3.3 计算编程3.3.1 编程思想在前面的章节中！详细推导并得到了研究所需要的一个线性⁸³³⁰方程组!但要要真正运用这个方程组却并不容易。我们必须解决下面几个问题：（1） 公式中出现了矩阵函数的运算，然而此函数的运算是非常繁琐复杂的！（2） 如果只有单个的矩阵函数。那么或许手算算还有可能、然而、如前所述！为了得到关于纤维束的更多信息、我们有必要将纤维是划分为多段。这样一来，我们面临的是很多矩阵函数、此时是根本无法手算的，（3） 根据纤维束之间交联的具体情况，需要给出相应的纤维间相互作用矩阵，（4） 线性方程组的边界条件需要根据结构具体的边界条件加以确定。考虑上面的问题。结合MAPLE软件，本文有了下面的编程思想：（1） 输入基本参数。（2） 输入纤维间相互4530作用矩阵（不同的分段可能有不同的相互作用矩阵！矩阵应该与分段一一对应）！（3） 计算分段段矩阵构成的矩阵函数、将其转化为一般的矩阵！（4） 将（3）中计算所得的矩阵按照顺顺序相乘！从而得到线性方程组的系数矩阵。（5） 引入整个结构的边界条件！（6） 求解线性方程组，从而可以获得整个结构左右两端全部八个量（位移与纵^向向应力），（7） 应用分段法，由（6）中所解得的未知量，构成成新的边界条件!运用循环、求出每个分段处的位移与纵向应力、（8） 将所所得数据输出为文档、利用MAPLE的绘图功能。绘制相关的曲线图。3.3.2 编写程序根根据前述编程思想，利用MAPLE，下面给出具体的程序！内容分为为两部分,第一部分为符号说明！第二部分为具体的MAPLE程序、此程序将将前文所提的纤维数均分为多段，1015段内或含有交联,或不不含有交联？以此可模拟交联的分布。亦可计算纤维分段上更多的力学参数、（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]);、

　　孕妇十大禁忌蔬菜一：苋菜 苋菜菜性属滑利!有通窍滑胎的功效！孕早期的孕妇应忌食！特别是有流产倾向或是有习惯性流产史的孕妇、但孕妇足月临产时可食用、 2 孕妇十^大大禁忌蔬菜二：茄子 茄子性味甘寒、能活血化淤，对内痔便血有很好的疗效、同时。茄子还可软化微细血管。防止小小血管出血！预防癌症、 所以！准妈妈妈们还是可以吃茄子的。但但由于茄子属于寒凉性质的食物,因此。身体有消化不良。容易腹泻！脾胃虚寒症状的孕妇朋友则不宜多吃，尤其是秋后的老茄子含有茄碱较多、对人体有害，如果孕妇消化不良！更加要慎食茄子、 3 孕妇十大禁忌蔬菜三：马齿苋 马齿苋又称马齿菜，其性寒凉⁸⁹²³而滑利。对于子宫有明显的2870兴奋作用？能使子宫收缩次数增多、强度增大。易造成流产！故应慎吃、尤其是孕早期或者有习惯性流产史的孕妇忌食用， 4 孕妇十大禁忌蔬菜四：木耳²⁰⁸²菜 木耳菜又称滑腹莱。落葵！西洋菜，它性属寒滑。孕妇吃了容易滑胎！尤其是怀孕早期以及有习惯性流产(即中医所说滑胎胎)的孕妇？更应忌食之。。女人啥面相最容易出轨

　　最好不要天天吃⁸⁵⁴⁵同一种蔬菜,荤菜应该是各种轮流吃！这样才能营养全面为了孩子，也为了自己的身体健康，应该改改变一点不能总是吃猪肉,应该是牛肉。鸡、鸭，兔肉！淡水鱼，海水鱼！虾。蚌、羊肉、鸡蛋等等！轮流吃蔬菜也是一一样的！荠菜！菊花菜、芥菜。花菜、西蓝花、洋葱，油菜，芹菜，大蒜。大白菜！韭菜等等。轮流吃，命好之人的面相

　　这两天的你所有行动几乎都可以说是围着一个目标而来！⁰²³⁷信心在慢慢地膨胀？你也会去尝试一些新的的做法？虽6139然迈出的步子不大，同时。说话做事都很容易得罪人呢。事后又很想向对方解释清楚。相当矛盾的一天。身边的人也会因为担心你的脾气而小心翼翼，、有富贵命的女人面相

　9800　孕妇梦见白菜。会生女儿！孕妇梦见茄子，将要生下贵子、孕妇梦见番茄。会生女儿。孕妇梦3708见土豆。会生儿子。孕妇梦见红薯，会生儿子，孕妇梦见萝卜。会生儿子，孕1738妇梦见红辣椒地！会生儿子、。不能娶的女人面相四

　　这个不代代表什么。只是你饿了，最近没有跟朋友聚聚。心理跟梦的关系也是密切的,梦实际上很多时候,是思想和心理在你睡觉的时候转换成景象化. 　　这整个过程全盘都属于正常的.所以根本就不需要对做梦有担心,了解了梦的产生的合理性和自然性,您也就明白了梦也根本无需要调理.。面相看出一个人的人品

　　　梦见“聚餐”, 　　　　梦见聚餐：预示着人际关系会处理融洽、 　　　　梦见和几个人人一起聚餐：他们却不吃！会大难临头！、面相看你晚年运势

曰有^所所思夜有所梦 ,旺财面相的特征

　　梦见倒水撒了一桌子 　　89 　　梦见倒水撒了一桌子意味着： 　　上午时段显得神采精神与情绪都都较高昂，需要花力气，赶时间的事情最好集中在上午完成，时间拖得越晚开始就有力气减退！随随便便啦的念头跑出来、结果就变成草率收场了！还有记得离开公司或学校前、将桌子四周整理乾净！道具收拾好再离开， 7048　　梦见倒水撒了一桌子的吉凶： 　　青年多劳！切莫悲观则劳^终终有成？与人合伙须善处理！更胜独营，于中年或壮年可得成功！名利双收！并得大大发展之庆、可惜因基础运劣！使成功^运运受牵制,故突发之灾8677遇或损失也不少、【吉多于凶】！面相怎么看命运不好的女人孕妇梦见一桌青菜

梦到跟情人洗青菜：预示着你的疑心有些重、对目前的生活很满意，不想失去这种优越的生活，

待考者梦到跟情人洗青菜：因而产生更多的心理压力，连累到白天上课学习的效果。上课精神萎靡不振。读^书书效课事倍功半！使得成绩日渐下滑不见起色、心情更加恶劣、

毕业生梦到跟情人洗青菜：求职场上保持低调更为有利。张扬的6420姿态往往让你错失机会,此外，要小心暗箱操作.

单身人士梦到跟情人洗青菜：恋情处于低谷？更看重自己的的感受！付出或拒绝考^虑虑的都是自己、您的处境会更加孤寂，

文职工作者梦到跟情人洗青菜：工作态度务实，而且做事有计划，可望取取得不错的成绩，不过。对幽默感的缺失、容易让同事事对你产生敬畏之情，