3.3 计算编程3.3.1 7320编程思想在前面的章节中,详细推导8372并得到了研究所需要的一个线性方程组！但要真正运用这个方程组却并不容易、我们必须解决下面几个问题：（1） 公式中出现了矩阵函数的运算！然而此函数的运算是非常繁琐复杂的！（2） 如果只有^单单个的矩阵函数。那么或许手算还有可能，然而、如前所述、1018为了得到关于纤维束的更多信息、我们有必要将纤维是划分为多段。这样一来、我们面临的^是是很多矩阵函数、此时是根本2772无法手算的？（3） 根据纤维束之间交联的具体情况。⁰²³³需要给出相应的纤维间相互作用矩阵!（4） 线性方程组的边界条件需要根据结构具体的边界条件加以确定，考虑上面的问题、结合MAPLE软件！本文有了下面的编程思想：（1） 输入基本参数！（2） 输入纤维间相互作用矩阵（不同的分段可能有不同的相互作用矩阵！矩矩阵应该与分段一一对应）？（3） 计算分段矩阵构成的矩阵函数，将其转化为一般的矩阵、（4） 将（3）中计算所得的矩阵按照顺序相乘，⁸⁷¹²从而得到线性方程组的系数矩阵,（5） 引入整个结构的边界条件！（6） 求解线性方程组、从而可以获得整个结构左右两端全部八个量（位移与纵向应力），（7） 应用分段法。由（6）中所解得的未知量，构成新的边界条件，运用循环、求出每个分段处的6386位移与纵向应力，（8） 将所得数据输出为文档、利用MAPLE的绘图功能。绘制相关的曲线图，3.3.2 编写程序根据前述编程思想！利用MAPLE。下面给出具体的程序！内容2158分为两部分。第一部分为符号说明。第二部分为具体的MAPLE程序、此程序将前文所提的8704纤维数均分为多段！段内或含有有交联、或不含有交联，以此可模拟交联的分布！亦可计算纤维分段上更多的力学参数！（1） 符号说明E：碳纳米管的弹性性模量,L：碳纳米管的长度。R：碳纳米管的半径！Mu：碳纳米管间的剪切模量！K：碳纳米管间的相互作用系数！Sigma：：施加的外力,A1！A2：碳纳米管间的相互作用矩阵，DL：分段的长度、B1、B2：⁶⁹²⁶矩阵函数转化为一般矩阵,JL：分段共价交联的信息,C：线性方程组系数矩阵、（2） 详细程序E := .46*10^12;L := 19.84*10^(-6);R := 1.5*10^(-6);Mu := .24*10^12;d := 3*R;k := mu/(R^2*ln(d/(2*R)+sqrt(d^2/(4*R^2)-1)));sigma := 10*10^9;A1 := Matrix(4, 4, [[0, k, 0, -k], [1/E, 0, 0, 0], [0, -k, 0, k], [0, 0, 1/E, 0]]);A2 := Matrix(4, 4, [[0, 0, 0, 0], [1/E, 0, 0, 0], [0, 0, 0, 0], [0, 0, 1/E, 0]]);with(LinearAlgebra);DL:= (1/100)*L;B1 := MatrixFunction(A1*DL, exp(x), x);B2 := MatrixFunction(A2*DL, exp(x), x);JL := readdata("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\JL.txt”,1 )C := Matrix(4, 4, [[1, 0, 0, 0], [0, 1, 0, 0], [0, 0, 1, 0], [0, 0, 0, 1]]):for i from 1 by 1 to 100 doif JL[i] = 1 thenC := B1 . CelseC := B2 . Cend ifend do:XS := evalf(C):Y := Vector[column](4, [t10, 0, 0, u20]):M := evalf(XS . Y):eqns := {0 = M[1], sigma = M[3], u110 = M[2], u210 = M[4]}:sols := evalf(solve(eqns, {t10, u110, u20, u210})):Y[1] := op(2, op(1, sols))：Y[4] := op(2, op(3, sols))：XSBL := Matrix(4, 4, [[1, 0, 0, 0], [0, 1, 0, 0], [0, 0, 1, 0], [0, 0, 0, 1]]):Z := Vector[column](4, [0, 0, 0, 0]):for i from 1 by 1 to 100 doIf JL[i]=1 thenXSBL:=B1.XSBL;elseXSBL:=B2.XSBL;end if;M := evalf(XSBL . Y);eqns := {z1 = M[1], z2 = M[2], z3 = M[3], z4 = M[4]};sols := evalf(solve(eqns, {z1, z2, z3, z4}));Z[1] := op(2, op(1, sols));Z[2] := op(2, op(2, sols));Z[3] := op(2, op(3, sols));Z[4] := op(2, op(4, sols));YL1 := array([[i*DL, Z[1]]]);writedata[APPEND]("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\YL1.txt",YL1 );YL2 := array([[i*DL, Z[3]]]):writedata[APPEND]("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\YL2.txt”,YL2 );WY1 := array([[i*DL, Z[2]]]);writedata[APPEND]("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\WY1.txt",WY1);WY2:=array([[i*DL,Z[4]]]);writedata[APPEND]("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\WY2.txt",WY2);If i=1 thenQYL1 := array([[i*DL, Z[1]-Y[1]]]);writedata[APPEND]("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\QYL1.txt",QYL1 );QYL2 := array([[i*DL, Z[3]-Y[3]]]);writedata[APPEND]("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\QYL2.txt",QYL2G);next end if;XSBL2 := Matrix(4, 4, [[1, 0, 0, 0], [0, 1, 0, 0], [0, 0, 1, 0], [0, 0, 0, 1]]);for j from 1 by 1 to i-1 doif JL[j]=1 thenXSBL2 := B1 . XSBL2;ElseXSBL2 := B2 . XSBL2;end if;end do;MM := evalf(XSBL2 . Y);eqns := {qz1 = MM[1], qz2 = MM[2], qz3 = MM[3], qz4 = MM[4]};sols := evalf(solve(eqns, {qz1, qz2, qz3, qz4}));QZ[1] := op(2, op(1, sols));QZ[2] := op(2, op(2, sols));QZ[3] := op(2, op(3, sols));QZ[4] := op(2, op(4, sols));QYL1 := array([[i*DL, Z[1]-QZ[1]]]);writedata[APPEND]("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\QYL1.txt",QYL1);QYL2 := array([[i*DL, Z[3]-QZ[3]]]);writedata[APPEND]("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\QYL2.txt",QYL2 );end do:YL1 := readdata("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\YL1.txt6",YL1 );YL2 := readdata("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\YL2.txt",YL2) ;QYL1 := readdata("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\QYL1.txt",QYL1) ;QYL2 := readdata("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\QYL2.txt",QYL2);WY1 := readdata("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\WY1.txt" ,WY1) ;WY2 := readdata("D:\\编程\\应力分布与位移分布的研究\\均匀分布\\WY2.txt",WY2);plot(YL1);plot(YL2);plot(QYL1);plot(QYL2);plot(WY1);plot(WY2);plot([YL1, QYL1]);plot([YL2, QYL2]);!

　　可以、黄花菜⁴³⁵⁰有消肿，利尿、解热。补血、健脑。通乳等作用,孕妇食用有好处。不过孕妇吃鲜黄花菜,需要经过水泡与充分的加热,。卫生间用什么灯具好 卫生间灯光风水学介绍

　　孕妇梦见兔子。小兔子。抓兔子！抓到一只兔子、追赶兔子都是好的征兆;但孕妇梦见多只兔子在睡觉，杀兔子则是不好的征兆、孕妇梦见兔子、更是预示着自己和孩子都会很健康、孕妇梦见小兔子。怀的如果是女孩多貌美！怀的是男孩多清秀、异性缘好、孕妇梦见抓兔子。预示着会生个活泼健康的的宝宝！孕妇梦见多只兔子在睡觉、提醒梦梦者在怀孕或生孩子方面，应该多加注意、家人也一定要照顾陪伴好梦者！尽量减少梦者独自一人行动的几率。孕妇梦3654见杀兔子，提提醒梦者要更多注意自己和孩子的身体健康！孕妇则要要当心流产，孕妇梦见抓到一只兔子！不久之后便会给家中添一个聪明的小宝宝，孕妇⁷⁰¹⁸梦见追赶兔子,预1969示着母子平安、是个好兆头，，卫生间的风水禁忌有哪些 6大禁忌介绍

　7810　今年一定是丰收年,如何看住宅风水 家居风水大师教你看住宅风水

　　〖梦见自己抱稻谷回家意味着〗：周围充满满了玩乐的气息，朋友会打电话来约你出去。家人也会有游玩的计划、不过随着夜幕的降临！后悔的感觉也在悄悄地迫近。让名人来告诫一下你吧。“放肆的生活！不是自自由的生活。””【佛滔算命】！卫生间门对客厅风水好吗 有什么化解方法

梦见母亲抱小孩：整体运势走低，现在正是休整的好时刻、不久后整装出发！会获得得不错的成绩！

梦见妈妈抱孩子！惊奇有趣，令你期待不已！

毕业生梦见母亲抱小孩预示找工作：求职运势良好？的具体执行能力增强。有计划也有行动！只要要不是要求太高、很有把握赢得考官赏识、争取到心水职位、

梦见女人抱小男孩，预示着你正在计划的事情或愿望会达成，是吉兆！

　　你好、很高兴能回答这位朋友的的题目!此梦说明6260：你要小心体能和精力的枯竭？瑜伽是可以帮你静心的运动！今天不妨一试！、风水知识大揭秘 客餐厅风水问题有哪些

梦见被抱：预示着你必须要善用智慧去经营与争取，只要要累积足够实力？是可以因自己的积极努力而赢得机会。遇到困难时别轻言放弃！

创业者者梦见被情人拥抱：预示着你最近的财运继续走高。家人对你的财务状况起支持作用。如果有较大开支。也往往因为家庭事务，久追未回的债款、行动会有所收获！

已婚者梦见被女人拥抱：预示着不久会出远门，但要防止节外生枝。可以避免一些挫折。是祥兆！

单身者梦见被情人拥抱：预示着近期在感情上会遇到一些小困难，只要能耐心心的去解决,不久会会得到好转!

职员⁵⁶²⁸梦见被女人拥抱：预示着你最近的工作状态很好！财运也会高涨！有聚财的的好现象,

求学者梦见被抱：预示着你的考^试试成绩很好？会去掉6366不错的进步，但不能太过骄傲，需要耐心！细心的学习才不会太过失落、

梦见被陌生人拥抱：预示着近期会有一些不受欢迎的可人上门！自己还需要虚情假意的应酬着，建议你你要做好心理准备和应对措施、

梦见被亲人拥抱：预示着家中会有人生病或生活不够快乐。自己很希望能够得到到别人更多的关心和爱护,

梦见被陌生人拥抱：预示着你多利用时间来审视下最近的得失。然后找到改进的方向再付诸实行。可望受到肯定。

梦见被抱的梦境

梦见自己拥抱陌生人！预示将有有不受欢迎的客人、

梦见拥抱亲戚、预示他们的疾病和不快。

男人梦见和自己女人抱在一起的情景、表示内心有强烈的欲望，

女人梦见和男人拥抱。预示会得到到异性朋友？但如果梦里是躺着拥抱男人，则暗示私生活不严谨！并惹来闲言碎语，

原版周公解梦见被抱

卧抱男人之女,凶！《原版周公解梦》

与人拥抱。主不祥。《原版周公解梦》

梦见被抱的案例解析

网友梦境：梦见过去情人带小孩我看小孩脸脏去给小孩洗脸我和情人拥抱亲吻，

解梦解析：梦^境境中的情人,是情感折磨的象征。

　　　说明最近会有艳遇，家居风水知识知多少 室内家居风水的三大讲究孕妇梦见抱花菜

　　孕妇梦见豹子、有时也提醒你周0460围有不忠实的朋友，小心别人散布谣言！对你恶意意中伤。！餐厅灯具风水禁忌不可忽略 选一款好灯很重要