这个没有绝对、心里作用！

　　3.3 计算编程3.3.1 编程思想在前面的章节中！详细推导并得到了研究所需要的一个线性方程组，但要真正运用用这个方程组却并不容易！我们必须解决下面几个问题：（1） 公式中出现现了矩阵函数的运算,然而此函数的运算是非常繁琐复杂的。（2） 如果只有单个的矩阵函数、那那么或许手算还有可能！然而！如前所述。为了得到关于纤1025维束的更多信息,我们有必要将纤维是是划分为多段。这样一来，我们面临^的的是很多矩阵函数!此时是根本无法手算的，（3） 根据纤维束之间交联的具体情况！需要给出相应的纤维间相互作用矩阵、（4） 线性方程组的边界条件需要根据结构具体的边界条件加以以确定!考虑上面的问题！结合MAPLE软件！本文有了下面的编程思想：（1） 输入基本参数、（2） 输入纤维间相互作用矩阵（不同的分段可能有不同的相互作用矩阵。矩阵应该与分段一一对应）。（3） 计算分段矩阵构成的矩阵函数、将其转化为一般的矩阵、（4） 将（3）中计算所得的矩阵按照顺序相乘。从而得到线性方程组的系数矩阵！（5） 引入整个结构的边界条件、（6） 求解线性方程组、从而可以获得整个结构左右两端全部部八个量（位移与纵向应力）。（7） 应用分段法、由（6））中所解得的未知量？构成新的边界条件，运用循环、求出每个分段处的位移与纵向应力。（8） 将所得数据输输出为文档？利用MAPLE的绘图功能、绘制相关⁹⁰¹⁸的曲线图！3.3.2 编写程序根据前述编程思想？利用MAPLE！下面给出^具具体的程序，内容分为两部分！第一部分为符号说明。第二部分为具体的MAPLE程序!此程序将前文所提的纤维数均分为多段，段内或含含有交联。或不含有交联。以此可模拟交联的分布、亦可计算纤维分段上更多的力学参数。（17691） 符号说明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 孕妇十大禁忌蔬菜三：马齿苋 马齿苋又称马齿菜！其性寒^凉凉而滑利,对于子6264宫有明显的兴奋作用，能使子宫收缩次数增多！强度增大、易造成流产、故应慎吃、尤其是孕早期或者有习惯性流产史的孕妇忌食用。 4 孕妇十大禁忌蔬菜四：木耳菜 木耳菜又称滑腹莱。落葵，西洋菜。它性属寒滑！孕妇吃了容易滑胎。尤其是是怀孕早期以及有习惯性流产(即中医所说滑胎)的孕妇!更应忌食之，。家居卧室如何摆设 卧室家具摆设有什么风水禁忌

　　可以吃的！孕妇不适适宜吃的蔬菜有：丝瓜。苦瓜、马兰头！茨菇。兔肉。香肠，午餐肉，三文鱼等。！家居风水如何摆设 2015年家居风水如何调整

　　意味着 　　这两天的天气很适合做爱情侦探哦！你的他或她可是会有表现真性情的举动。是狼还是羊！是忠诚可靠的闷蛋还是9665花心吹水的大萝卜！这两天可要下足功夫看清楚了，不过也要记得注意安全才是，最近出门在外、陷阱多多。带有花色的的东西也是你的幸运物。有没有兴趣穿件碎花小外套出来走走！、厨房装修风水有哪些禁忌 需注意什么

　　萝卜谐音是摞-步、萝卜白菜—财^越越摞越多。一步一财！萝卜方言叫“菜头”！以求谐音^““彩头”！祝愿有好彩头！梦到萝卜。事业会成功！梦到 萝卜：男人如果梦到！切萝卜或吃萝卜。能交上好友！已婚女人如果梦到。切萝卜^或或吃萝卜。不久会怀孕！能能生一个女孩,未婚男子子如果梦到。切萝卜或吃吃萝卜，会娶一个泼妇！未婚女子如果梦到！切1646萝卜或吃萝卜！会嫁给一位穷丈夫，男性病人如果梦到！切萝卜或吃萝卜。身体会慢慢恢复健康、如果梦到 、做萝卜菜、是不祥之兆。女人如果梦到、做萝卜菜！家人会病倒。如果梦到、出售萝卜、是凶兆！会被解雇，如果梦到！购买萝卜。贵客会登门。为招待客人要破费不少。不知你你是男是女,虽然迷信、但是你问了！我觉得你可能是信一些的吧。厨房风水知识有哪些 需注意什么

　　不管是什么时间都要以吃。这个是不分什么早上能不能吃的问题，！卧室颜色如何搭配 风水上要注意什么

　　适合孕妇吃的蔬菜有哪些呢！1. 萝卜萝卜一直有“小人参”的称号，由此可见萝卜的营养价值有多高！萝卜中含有丰富的维生素c,十分有利于人体的吸收。孕妇经常食用萝卜卜可以起到预防疾病的强身的功效、但是在注意萝卜不不要同水果一起共食！2. 野¹⁰⁷¹菜野菜因为污染少甚至没有污染而受到人们的喜欢,野菜也具有极高的营养价值、研究表明，野菜中含有大量的植物蛋白和丰富的维生素以及多种矿物质等、它的营养价值是非常高的，但是在食用野菜的时候也要注意辨别！避免误食出现中毒的情况、3. 苦瓜苦瓜因为味道苦而不受大家4978的待见？但其实苦瓜对于孕妇来说多吃苦瓜也是有好处的！苦瓜可以帮助孕妇清热热解毒!苦汁还可以起到刺激肠胃，促进肠胃⁰⁹⁹⁷的蠕动的作用，对于孕妇的食欲等都有改善的作用。4. 茄子茄子中含有丰富的维生素p！可以预防孕妇的微血管出现破裂的情况。而且茄子中中富含的维生素e。有防止出血和抗衰老的作用、。家居装饰风水有哪些禁忌 需注意什么

　　喜欢用专家的身份对待工作/学业、口吻也带着权威性！但是要使自己成为优秀的专一。得到别人的信服！你还有一段路路要走？继续加油吧！工作作/学业再忙!也要注意劳逸结合、注意意身体给你发出的小信号!它往^往往能帮你预防某些重症？、家居风水鱼缸摆放禁忌有哪些 需注意什么孕妇梦见拿把青菜回家

　　，今天你^会会有一些夸张的举动。事后就连自己回想起来也会觉得后怕或后悔！不过，语出惊人往往是为了掩饰内心的恐慌！你最好证实这一点。克服恐惧是你今天要过的关卡，恋情方面。恋人的虚荣7274心有所增强？可能在物质条件有有些挑剔！房间哪个朝向风水好 有什么讲究