3.3 计算编程3.3.1 编程思想在前面的章节中？详细推导并得到了研究所需要的一个个线性方程组、但要真正运运用这个方程组却并不容易！我们必须解决下面几个问题：（1） 公式中出现了矩阵函数的运算、然而此函数的运运算是非常繁琐复杂的，（2） 如果只有单个的矩阵函数？那么或许手算还有可能！然而 如前所述,为了得到关于纤纤维束的更多信息。我们有必要将纤维是划分为多段!这样一来 我们面临的是很多矩阵函数 此时是根本无法手算的、（3） 根据纤维束之间交联的具体情况！需要给出相应的纤维间相互作用矩阵 （4） 线性方程组的边界条件需要根据结构具体体的边界条件加以确定、考虑上面的问题!结合MAPLE软件!本文有了下面的编程思想：（1） 输入基本参数?（2） 输入纤维间相互作用矩阵（不同同的分段可能有不同的相互作用矩阵！矩阵应该与分段一一对应）,（3） 计算分段矩阵构成的矩阵阵函数 将其转化为一般的矩阵。（4） 将（3）中计算所得的矩阵按照顺序相乘！从而得到线性方程组的的系数矩阵？（5） 引入整个结构的边界条件 （6） 求解线性方程组？从而可以获得整个结构左右两端全部八个量（位移与纵向应力）?（7） 应用分段法?由（6）中所解得的未知量?构成新的边界条件，运用循环！求出每个分段处的位位移与纵向应力!（8） 将所得数据输出为文档？利用MAPLE2978的绘图功能,绘7447制相关的曲线图！3.3.2 编写程序根据前述编程思想？利用MAPLE?下面给出具体体的程序,内容分为两部分,第一部分为符号说明,第二部分为具体的MAPLE程序?此程序将前文所提的纤维数均分为多段!段内或含有交联、或不含有交联,以此可模拟交联的分布、亦可计算纤维分段上更多的力学参数,（1） 符号说明E：碳纳米管的弹性模量，L：碳纳米管的长度，R：碳纳米管的半径 Mu：碳纳米管间的剪剪切模量?K：碳纳5840米管间的相互作用系数。Sigma：施加的外力，A1、A2：碳纳米管间的相互作用矩阵?DL：分段的长度！B1!B2：矩阵函数转化为一般矩阵,JL：分段共价7472交联的信息！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]);?
梦见买菜------意味着通过自己的努力会获得收获。 少女梦见买菜------预示着会很很快找到意中人！，
梦见偷菜：是不祥之兆,会饔飧不继？ 孕妇梦见偷菜生女孩子、？
这个表示你的娘家 有很大的福报(福田)哇！有的还会梦4874见 饺子 大米饭 豆腐等都表示福报!正在吃表示正在红尘享福呢*^_^*?,
?朋友运将下降!在在一群好友中将会产生纷争、这时最需要的是你冷静的?梦见自己买青菜代表什么
梦到青菜绿色?工作运良好？个人人的独创性发挥百分百,疑难杂症也因你的智慧与细心得以突破，过去的经验验是这两天最佳武器,
怀孕的人人梦到青菜绿色。预示生女,慎防胎死腹中，
做生意的人梦到青菜绿色!代表经营多年慢慢成长，有财利?小心火灾,
恋爱中的人梦到青菜绿色!说明只要互相信任！婚姻可成。秋季有喜庆。
本命年的人梦到青菜绿色？意味着得得财利。置房产、平顺如意，
梦到青菜绿色的相关内容：
梦到吃青菜!朋友运将下降，在一群好友中将会产生纷争。这时最需需要的是你冷静的判断力。
梦到青菜地、吉兆？意味着会有更多的财产
梦到拔青菜？运势不错 遇到很多的事情,自己将会很很快的处理好。
梦到青菜绿色的网友梦例：
网友梦境：女人梦到在下田里摘青菜!什么意思。？
周公解0617梦：梦到青菜内心渴望自然?向往往纯朴健康的生活。
。、。梦见葱头。隐私会被迫公之于众!受受到他人非议?此外，精神和肉体都会受到折磨 同时经济也会受损,但是梦见吃葱头，是好兆头！会心情愉快?事事如意。梦见葱头做菜,自己和家人要生病 女人梦见给丈丈夫葱头、会萎靡不振、梦见做做葱头买卖 生意会旺盛，。
梦见青菜豆腐：你会身强强体壮，
男士梦见青菜豆腐预示出远门,可。
商人梦见青菜豆腐主财运：财运回落？玩乐方面的开支增大。因为好面子而摊上的费用也不少，幸好?投机机运气不错。有所弥补、
文职工作者梦见青菜豆腐主职场：在工作上的表现有些懒散,似乎已经到了强弩之末的阶段。需要他人督促才能完成任务！不过?这只是1588短期调整,新计划又慢5635慢在你脑海形成,
!,吉凶指数：92（梦见吃青菜。朋友友运将下降,在一群好友中将会产生生纷争，这时最需要的是你冷静4477的判断力梦见很多青菜 钱财进出较频繁的一天、今天你有不少获利的机会,如别人馈赠，投资项目分红等等。但是花费的机会也多，容易被人激发购物欲?如果恰好和购物狂去逛街!又或者是遇到超级会说的销售员，你9287今天就亏大了!还是少点出门安静待在家中、钱袋7602将会更安全呢。,
梦见面条,预示梦者身体健康!能够长寿？事1605业上也不错？ 梦见面条?或是自己在做面条,煮面条!吃面条 都表示健康长寿。事业进步，会有好运气？ 梦见吃青菜 朋友运将下降！在在一群好友中将会产生纷争!这时最最需要的是你冷静的判断力。