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

　　　人的梦就是人的潜意识，梦是潜意识的一种释放！去^理理解自己的梦也可以更好的了解自己？在梦里某人某事某物都是潜意识的释放、这梦有反应你呼吸系统有问题。需要注意。住宅楼梯位置风水

　　好梦！以后会发财，住宅楼梯风水

　　预1514示要经受坎坷，磨难……，住宅楼风水是锅炉房

　　白9227天干活儿了吗!、住宅楼风水箭刹图片

怀孕的人梦见烂菜！预示生女！冬占生男，

做生意的人梦见烂菜，代表五金机械业！贸易易中盘商得财利,要有恒心。

恋爱中的人梦见烂菜！说2542明反目成仇？感情要时间来培养感情，婚姻慢慢可成！

本命年的人梦见烂菜、预示宜守。等待时机！贵人扶助！顺利如意。

梦见烂菜梦境

梦见蔬菜。通常预示你会获得成功，

梦见有许多蔬菜，表示生活活舒适富有、

梦见自己在切菜？暗示你要当心被朋朋友欺骗！遭受经济损失。

梦见烂菜的案例解析

网友梦境：梦见收拾烂菜叶怎么意思、

解梦解析：代表五五金机械业、贸易^中中盘商得财利！要有恒心。

梦见烂肉：运势不佳、遇到难题的时候自己不知道怎么办。建议需要学会如何解决。

未婚者者梦见烂肉：运势不错、恋情方面遇到挫折！记得一定要有长久的的耐心？继续的坚持下去！将将会获得一定的成功！

1062求职者梦见烂肉：求职运一般？做好充足的准备！迎合招聘方的需求！拿捏好对方提出的问题！从中获取收获，

未成年梦见烂肉：运势不错。身体方面没有任何不适的出现。说明自己把身体保养得不错。

准备参加考试梦见烂肉：运势不佳，^考考试的时候太过紧张。因此就会影响考试的心情，成绩当然就会不好！

梦见烂肉的梦境

梦见腐烂的肉！是生病的预兆，

梦见吃吃肥肉，预示恋爱顺利、

梦见吃大块的肉，还预示收入可能会增多，

原版周公解梦见烂肉

食自死肉！主离别，《原版周公解梦》

梦见食生肉、凶，熟肉。吉，《敦煌本梦书》

梦见食生肉、忧县官事、《敦煌本梦书》

梦见烂肉的案例解析

网友梦境：梦见把一盘肥肉烧糊了！吃烧糊的肥肉、

解梦解析：梦见吃肉是吉兆！说明会有好事发生，

　　预示着，6077你对日 常生活料理精细！注意营养。祝你幸福快乐开心，希采纳！住宅水井风水

　　^到到尽头了，住宅水井风水注意什么梦见买烂生姜

　　告诉你吧，怀孕如果梦见蛇必定是男孩。，住宅水池风水