3.3 计算编程3.3.1 编程思想在前面的章节中,详细推导并得到了研究所需要的一个线性方程组 但要真正运用这个方程组却并不容易?我们必须解决下面几个问题:((1) 公式中出现了矩阵函数的运算!然而此此函数的运算是非常繁琐复杂的。(2) 如果只有单个个的矩阵函数 那么或许手算还有可能?然而 如前所述!为了得到关于纤维束的更多信息。我们有有必要将纤维是划分为多段,这样一来,我们面临的是很多矩阵函数!此时是根本无法手算的,(3) 根据纤维束之间间交联的具体情况,需要给出相应的纤维间相互作用用矩阵、(4) 线性方程组的边界条件需要根据结构具体的边界条件加以确定?考虑上面的问题!结合MAPLE软件!本文有了下面的编程思想:(1) 输入基本参数。(2) 输入纤维间相互作用矩阵(不同的分段可能有不同的相互作用矩阵?矩阵应该与分段一一对应)!(3) 计算分段矩阵构成的矩阵函数,将其转化为一般的矩阵?(4) 将(3)中计算所得的矩阵按照顺序相乘?从而得到线性方程组的系数矩阵 (5) 引入整个结构的边界条件 (6) 求解线性方程组 从而可以获得整个结构左右两端全全部八个量(位移与纵向应力)!(7) 应用分段法!由(6)中所解8984得的未知量?构成新的边界条件,运用循环?求出每个分段处的的位移与纵向应力、(8) 将所得得数据输出为文档!利用MAPLE的绘图功能,绘制制相关的曲线图 3.3.2 编写程序根据前述编程思想,利用MAPLE。下面给出具体的程序!内内容分为两部分、第第一部分为符号说明、第二部分为具体9595的MAPLE程序 此程序将前文所提的纤维数均分为多段!段内或含2827有交联。或不含有交联 以此可模拟交联的分布?亦可计算纤维分段上更多的力学参数,(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]);?
这两天所做的事往往成果和心里所期望的有些差距!可能先不不要抱什麽太大期待比较好、如果心中有什麽不满就坦率点说出口吧 越是闷在心里爆发时越可可怕喔,该适度表达心声的时候就发表出来,不不过自己的事虽没有什麽特别突出的地方。帮帮别人倒是不错的,尤其是同事、晚辈有困难时伸个援手,从中多6322少能得到些助人的成果吧,、
7564梦见吃青菜?朋友运将下降 在一群好友中将会产生纷争,这时最需要的是你冷静的判断力?
梦见洗青菜,外出时会遭强盗袭击、
女人梦见洗青菜!事事如意?喜气洋洋、但是不能因此而散漫!流连於欢乐场所,招致坏坏运与劳苦、小心小偷闯空门,
年轻人梦见洗青菜,需继续关注呼吸系统疾病 肩部和手部也有受伤的可能,
梦见自己去种青菜?预示在事业成功的道路上 你必须付出艰苦的努力,战胜许6932多困难,
孕妇梦见青菜?预预示会生儿子,
孕妇梦见大白菜!预示会生女儿?
。!,应该表示会失财,丢失物品。近期少少带贵重物品在身上,以免造成损失和意外 周公解梦[!
梦见买青菜?朋友运将将下降,在一群好友中将会会产生纷争,这时最需要的是你冷静的判断力。?
狗进家门?主人兴!
梦见油菜是福!梦见自己吃油菜。预示着5725生活上会一切顺利!梦见自己种种油菜!去菜园里采油菜。预示着自己做任何事情都会顺利成功,梦见自己挑食油菜?不喜欢吃,预示着会有困难出现在自己前进道路上?
梦见洗青菜有头发!如果勉强著与他1920人交际让你觉得浑身不自在的话那就明白表示态度吧!这两天不妨早起!出去散散步有提升运势的效果、约会应选择浪漫的地点,效果也会加倍??
1!清水浸泡洗涤法:先用清水冲洗掉蔬菜的的表面污物 然后再用清水浸泡,不过浸泡时间最好不少于10分钟,这样清洗2472浸泡2次后基本上就可以清除绝大部分残留的农成分。2?添加辅助剂浸泡洗涤法:可以添加碱(小苏达),果蔬清洗剂!淘米水(最好用头两7359次的)进行浸泡洗涤!将蔬菜瓜果在这样的水中浸泡5到15分钟,可以去除蔬果表面所含的有机磷杀虫剂、浸泡后请注意要完全冲洗干净 提示:对于叶类蔬菜菜来说 比较好的祛除农残留的办法是用水浸泡,最好往水里加点果蔬清洗剂?以加速农的溶出!
梦见偷菜:是不祥之兆?会饔飧不继。 孕妇梦见偷菜生女孩子? 梦见和别人买青菜