復(fù)合材料力學(xué)上機(jī)編程作業(yè)(計(jì)算層合板剛度)(共10頁(yè))

1、第 頁(yè)復(fù)合材料(f h ci lio)力學(xué)上機(jī)編程作業(yè)學(xué)院(xuyun)：School of Civil Engineering 專業(yè)(zhuny)：Engineering Mechanics小組成員信息：James Wilson（2012031890015）、Tau Young（2012031890011） 復(fù)合材料力學(xué)學(xué)了五個(gè)星期，這是這門(mén)課的第一次編程作業(yè)。我和楊濤結(jié)成一個(gè)小組，我用的是Fortran編制的程序，Tau Young用的是matlab編制。其中的算例以我的Fortran計(jì)算結(jié)果為準(zhǔn)。Matlab作為可視化界面有其獨(dú)到之處，在附錄2中將會(huì)有所展示。 作業(yè)的內(nèi)容是層合板的剛度的

2、計(jì)算和驗(yàn)算，包括拉伸剛度A、彎曲剛度D以及耦合剛度B。首先要給定層合板的各個(gè)參數(shù)，具體有：層合板的層數(shù)N；各單層的彈性常數(shù)E1、E2、G12；各單層對(duì)應(yīng)的厚度；各單層對(duì)應(yīng)的主方向夾角。然后就要計(jì)算每個(gè)單層板的二維剛度矩陣Q，具體公式如下：；得到Q矩陣后，根據(jù)課本上講到的得到。 然后根據(jù)z坐標(biāo)的定義求出到，接下來(lái)，最重要的一步，根據(jù)下式計(jì)算A、B、D。一、書(shū)上P110的幾個(gè)問(wèn)題可以歸納為以下幾個(gè)類型。（1）正交鋪設(shè)5層對(duì)稱層合板（T5-7）數(shù)據(jù)文檔層數(shù)5 層序數(shù)厚度mE1（Pa）E2（Pa）v12v21G12（Pa）角度（）1 1.00E-039.60E+102.40E+100.10 0.40

3、1.00E+100.00 2 1.00E-039.60E+102.40E+100.10 0.40 1.00E+1090.00 3 1.00E-039.60E+102.40E+100.10 0.40 1.00E+100.00 4 1.00E-039.60E+102.40E+100.10 0.40 1.00E+1090.00 5 1.00E-039.60E+102.40E+100.10 0.40 1.00E+100.00 結(jié)果文檔 拉伸剛度A= 3.5000E+08 5.0000E+07 -4.3711E-01 5.0000E+07 2.7500E+08 -6.1196E+00 -4.3711E-

4、01 -6.1196E+00 5.0000E+07 耦合剛度B= -9.0523E-12 -7.2831E-13 -3.8514E-21 -7.2831E-13 -2.2631E-12 -7.0304E-20 -3.8514E-21 -9.1480E-20 -9.0417E-13 彎曲剛度D= 8.7917E+02 1.0417E+02 -4.7354E-07 1.0417E+02 4.2292E+02 -6.6296E-06 -4.7354E-07 -6.6296E-06 1.0417E+02由此可以從課本上了解到的：A16=A26=0；D16=D26=0；相吻合。這里B顯然是等于零的。正交

5、鋪設(shè)(p sh)6層反對(duì)稱層合板（T5-8）數(shù)據(jù)文檔層數(shù)6 層序數(shù)厚度mE1（Pa）E2（Pa）v12v21G12（Pa）角度（）1 1.00E-039.60E+102.40E+100.10 0.40 1.00E+100.00 2 1.00E-039.60E+102.40E+100.10 0.40 1.00E+1090.00 3 1.00E-039.60E+102.40E+100.10 0.40 1.00E+100.00 4 1.00E-039.60E+102.40E+100.10 0.40 1.00E+1090.00 5 1.00E-039.60E+102.40E+100.10 0.40 1

6、.00E+100.00 6 1.00E-039.60E+102.40E+100.10 0.40 1.00E+1090.00 結(jié)果文檔 拉伸剛度A= 3.7500E+08 6.0000E+07 -6.5567E-01 6.0000E+07 3.7500E+08 -9.1794E+00 -6.5567E-01 -9.1794E+00 6.0000E+07 耦合剛度B= -1.1250E+05 1.8943E-10 -3.2784E-04 1.8943E-10 1.1250E+05 -4.5897E-03 -3.2784E-04 -4.5897E-03 1.8102E-10 彎曲剛度D= 1.125

7、0E+03 1.8000E+02 -1.9670E-06 1.8000E+02 1.1250E+03 -2.7538E-05 -1.9670E-06 -2.7538E-05 1.8000E+02由此可以和課本上了解到的：A11=A22；D11=D22；A16=A26=D16=D26=0；相吻合。（3）5層對(duì)稱(duchn)角鋪設(shè)層合板（T5-9）數(shù)據(jù)文檔層數(shù)5 層序數(shù)厚度mE1（Pa）E2（Pa）v12v21G12（Pa）角度（）1 1.00E-039.60E+102.40E+100.10 0.40 1.00E+1030.00 2 1.00E-039.60E+102.40E+100.10 0.

8、40 1.00E+10-30.00 3 1.00E-039.60E+102.40E+100.10 0.40 1.00E+1030.00 4 1.00E-039.60E+102.40E+100.10 0.40 1.00E+10-30.00 5 1.00E-039.60E+102.40E+100.10 0.40 1.00E+1030.00 結(jié)果文檔 拉伸剛度A= 3.4531E+08 1.1094E+08 2.3274E+07 1.1094E+08 1.5781E+08 9.2015E+06 2.3274E+07 9.2015E+06 1.1094E+08 耦合剛度B= -1.2975E-11 -

9、1.4779E-12 1.4282E-12 -8.4022E-12 1.3998E-12 -6.5192E-13 1.4282E-12 -8.2956E-13 -8.4022E-12 彎曲剛度D= 7.1940E+02 2.3112E+02 1.4159E+02 2.3112E+02 3.2878E+02 5.5976E+01 1.4159E+02 5.5976E+01 2.3112E+02由此可以和課本上了解到的：A16、A26相對(duì)要小；D16、D26相對(duì)要?。幌辔呛?。這里B顯然是等于零的。（4）6層反對(duì)(fndu)稱角鋪設(shè)層合板（T5-10）數(shù)據(jù)文檔層數(shù)6 層序數(shù)厚度mE1（Pa）E2（P

10、a）v12v21G12（Pa）角度（）1 1.00E-039.60E+102.40E+100.10 0.40 1.00E+1030.00 2 1.00E-039.60E+102.40E+100.10 0.40 1.00E+10-30.00 3 1.00E-039.60E+102.40E+100.10 0.40 1.00E+1030.00 4 1.00E-039.60E+102.40E+100.10 0.40 1.00E+10-30.00 5 1.00E-039.60E+102.40E+100.10 0.40 1.00E+1030.00 6 1.00E-039.60E+102.40E+100.1

11、0 0.40 1.00E+10-30.00 結(jié)果文檔 拉伸剛度A= 4.1437E+08 1.3313E+08 -1.3479E-09 1.3313E+08 1.8938E+08 3.0013E-10 -1.3479E-09 3.4470E-10 1.3313E+08 耦合剛度B= 1.2989E-11 -9.3543E-12 -6.9823E+04 7.4607E-12 -2.7711E-12 -2.7605E+04 -6.9823E+04 -2.7605E+04 7.4607E-12 彎曲剛度D= 1.2431E+03 3.9938E+02 2.9643E-14 3.9938E+02 5.

12、6813E+02 9.5583E-15 2.9643E-14 7.4281E-15 3.9938E+02由此可以和課本上了解到的：A16=A26=D16=D26=0;B11=B22=B12=B21=B66=0；相吻合。（5）我還想驗(yàn)證一個(gè)書(shū)上的例題，在課本P114。三層層合板，外層厚度t1，內(nèi)層10t1，正交鋪設(shè)比m=0.2,。玻璃(b l)/環(huán)氧單層板性能：E1=5.4E10Pa，E2=1.8E10Pa，v21=0.25，G12=8.8E9Pa。數(shù)據(jù)文檔層數(shù)3 層序數(shù)厚度mE1（Pa）E2（Pa）v12v21G12（Pa）角度（）1 1.00E+005.40E+101.80E+100.083

13、 0.250 8.80E+090.00 2 1.00E+015.40E+101.80E+100.083 0.250 8.80E+0990.00 3 1.00E+005.40E+101.80E+100.083 0.250 8.80E+090.00 結(jié)果文檔 The NORMALISED extension stiffness A* equals: 2.4509E+10 4.5954E+09 -1.3893E+02 4.5954E+09 4.9017E+10 -1.2002E+03 -1.3893E+02 -1.2002E+03 8.8000E+09 The NORMALISED coupling

14、 tensor C equals: -9.3609E+09 1.4411E-05 -4.2450E+01 1.4411E-05 9.3609E+09 -3.6673E+02 -4.2450E+01 -3.6673E+02 1.6107E-05 The NORMALISED bending stiffness D* equals: 3.3869E+10 4.5954E+09 -9.6477E+01 4.5954E+09 3.9656E+10 -8.3347E+02 -9.6477E+01 -8.3347E+02 8.8000E+09由于課本上只是分析Nx的荷載，給出了A*的數(shù)值和課本上計(jì)算的結(jié)果

15、一致。二、驗(yàn)證Verchery的論文(lnwn)里給出的數(shù)值算例。這里一直(yzh)到Table5的數(shù)據(jù)都是從Verchery的論文中截獲。Verchery論文中的18層序列，第（21）式【laminates without bending-extension coupling】的排列有兩種材料(cilio)，一種是Boron-Epoxy，一種是Glass-Epoxy。而且都給出了最終的計(jì)算結(jié)果Q，A*，D*。下面是我的Fortran計(jì)算數(shù)據(jù)文檔和結(jié)果文檔。（1）Boron-Epoxy材料。（Boron-Epoxy）數(shù)據(jù)文檔層數(shù)18 層序數(shù) 厚度m E1（Pa）E2（Pa）v12v21G12（

16、Pa）角度（）1 1.00E-032.04E+111.85E+100.021 0.230 5.59E+090.00 2 1.00E-032.04E+111.85E+100.021 0.230 5.59E+090.00 3 1.00E-032.04E+111.85E+100.021 0.230 5.59E+0960.00 4 1.00E-032.04E+111.85E+100.021 0.230 5.59E+0960.00 5 1.00E-032.04E+111.85E+100.021 0.230 5.59E+09-60.00 6 1.00E-032.04E+111.85E+100.021 0.

17、230 5.59E+09-60.00 7 1.00E-032.04E+111.85E+100.021 0.230 5.59E+09-60.00 8 1.00E-032.04E+111.85E+100.021 0.230 5.59E+0960.00 9 1.00E-032.04E+111.85E+100.021 0.230 5.59E+090.00 10 1.00E-032.04E+111.85E+100.021 0.230 5.59E+09-60.00 11 1.00E-032.04E+111.85E+100.021 0.230 5.59E+09-60.00 12 1.00E-032.04E+

18、111.85E+100.021 0.230 5.59E+0960.00 13 1.00E-032.04E+111.85E+100.021 0.230 5.59E+0960.00 14 1.00E-032.04E+111.85E+100.021 0.230 5.59E+090.00 15 1.00E-032.04E+111.85E+100.021 0.230 5.59E+090.00 16 1.00E-032.04E+111.85E+100.021 0.230 5.59E+090.00 17 1.00E-032.04E+111.85E+100.021 0.230 5.59E+0960.00 18

19、 1.00E-032.04E+111.85E+100.021 0.230 5.59E+09-60.00 （Boron-Epoxy）結(jié)果文檔The stiffness of number 1 ply is: 2.0499E+11 4.2757E+09 0.0000E+00 4.2757E+09 1.8590E+10 0.0000E+00 0.0000E+00 0.0000E+00 5.5900E+09 The coupling stiffness B equals: -5.6581E-10 2.1876E-10 -1.4447E-10 2.1876E-10 6.5563E-10 -5.0525E

20、-10 -1.4447E-10 -4.5370E-10 2.5395E-10 The NORMALISED extension stiffness A* equals: 8.7706E+10 2.8359E+10 -6.5585E-07 2.8359E+10 8.7706E+10 -2.9973E-06 -6.5585E-07 -3.0072E-06 2.9674E+10 The NORMALISED coupling tensor C equals: -2.1719E+10 4.4599E+09 -1.4444E+09 4.4599E+09 1.2799E+10 -4.5344E+09 -1

21、.4444E+09 -4.5344E+09 4.4599E+09 The NORMALISED bending stiffness D* equals: 1.0943E+11 2.3899E+10 1.4444E+09 2.3899E+10 7.4907E+10 4.5344E+09 1.4444E+09 4.5344E+09 2.5214E+10這里的結(jié)果顯然是B=0，而且得到的Q，A*，D*與論文上的數(shù)據(jù)一致，Glass-Epoxy材料(cilio)。（Glass-Epoxy）數(shù)據(jù)文檔層數(shù)18 層序數(shù) 厚度m E1（Pa）E2（Pa）v12v21G12（Pa）角度（）1 1.00E-033

22、.86E+108.27E+090.056 0.260 4.14E+090.00 2 1.00E-033.86E+108.27E+090.056 0.260 4.14E+090.00 3 1.00E-033.86E+108.27E+090.056 0.260 4.14E+0960.00 4 1.00E-033.86E+108.27E+090.056 0.260 4.14E+0960.00 5 1.00E-033.86E+108.27E+090.056 0.260 4.14E+09-60.00 6 1.00E-033.86E+108.27E+090.056 0.260 4.14E+09-60.00

23、 7 1.00E-033.86E+108.27E+090.056 0.260 4.14E+09-60.00 8 1.00E-033.86E+108.27E+090.056 0.260 4.14E+0960.00 9 1.00E-033.86E+108.27E+090.056 0.260 4.14E+090.00 10 1.00E-033.86E+108.27E+090.056 0.260 4.14E+09-60.00 11 1.00E-033.86E+108.27E+090.056 0.260 4.14E+09-60.00 12 1.00E-033.86E+108.27E+090.056 0.

24、260 4.14E+0960.00 13 1.00E-033.86E+108.27E+090.056 0.260 4.14E+0960.00 14 1.00E-033.86E+108.27E+090.056 0.260 4.14E+090.00 15 1.00E-033.86E+108.27E+090.056 0.260 4.14E+090.00 16 1.00E-033.86E+108.27E+090.056 0.260 4.14E+090.00 17 1.00E-033.86E+108.27E+090.056 0.260 4.14E+0960.00 18 1.00E-033.86E+108

25、.27E+090.056 0.260 4.14E+09-60.00 （Glass-Epoxy）結(jié)果文檔The stiffness of number 1 ply is: 3.9170E+10 2.1820E+09 0.0000E+00 2.1820E+09 8.3922E+09 0.0000E+00 0.0000E+00 0.0000E+00 4.1400E+09 The coupling stiffness B equals: -3.2222E-10 8.5606E-11 -3.2021E-11 8.5606E-11 2.1697E-10 -8.1172E-11 -3.4829E-11 -8

26、.1172E-11 -2.0606E-13 The The NORMALISED extension stiffness A* equals: 2.0451E+10 5.5118E+09 -1.4616E-07 5.5118E+09 2.0451E+10 -4.3469E-07 -1.4491E-07 -4.3469E-07 7.4698E+09 The NORMALISED coupling tensor C equals: -3.4665E+09 6.1663E+08 -2.8000E+08 6.1663E+08 2.2332E+09 -7.0721E+08 -2.8000E+08 -7.

27、0721E+08 6.1663E+08 The NORMALISED bending stiffness D* equals: 2.3918E+10 4.8952E+09 2.8000E+08 4.8952E+09 1.8218E+10 7.0721E+08 2.8000E+08 7.0721E+08 6.8532E+09這里的結(jié)果顯然是B=0，而且得到的Q，A*，D*仍然與論文上的數(shù)據(jù)一致。（3）當(dāng)然我也驗(yàn)證了第（22）【laminates with equal elastic properties in bending and extension】、（23）【quasi-homogeneo

28、us laminates】的排序(pi x)，材料是Boron-Epoxy，下面給出計(jì)算的結(jié)果。從下面的兩個(gè)(lin )結(jié)果表中可以知道，（22）排序的確是C=0，（23）的排序的確是B=0且C=0。驗(yàn)證成功。（A）第（22）排序(pi x)。（Boron-Epoxy）結(jié)果文檔 The ACTUAL stiffness tensor of the laminate: The extension stiffness A equals: 1.5787E+09 5.1047E+08 -1.9256E-08 5.1047E+08 1.5787E+09 -3.9050E-08 -1.9256E-08 -

29、5.4129E-08 5.3412E+08 The coupling stiffness B equals: -4.0463E+06 8.3088E+05 4.8750E+05 8.3088E+05 2.3845E+06 1.5303E+06 4.8750E+05 1.5303E+06 8.3088E+05 The bending stiffness D equals: 4.2625E+04 1.3783E+04 -7.3708E-13 1.3783E+04 4.2625E+04 -4.1602E-12 -7.3708E-13 -3.8027E-12 1.4421E+04 The NORMAL

30、ISED stiffness tensor of the laminate: The NORMALISED extension stiffness A* equals: 8.7706E+10 2.8359E+10 -1.0698E-06 2.8359E+10 8.7706E+10 -2.1694E-06 -1.0698E-06 -3.0072E-06 2.9674E+10 The NORMALISED coupling tensor C equals: 1.4275E-05 5.7798E-06 4.4685E-07 5.7798E-06 8.9034E-06 6.3907E-06 4.468

31、5E-07 4.8174E-06 2.6114E-06 The NORMALISED bending stiffness D* equals: 8.7706E+10 2.8359E+10 -1.5166E-06 2.8359E+10 8.7706E+10 -8.5601E-06 -1.5166E-06 -7.8246E-06 2.9674E+10（B）第（23）排序。（Boron-Epoxy）結(jié)果文檔 The ACTUAL stiffness tensor of the laminate: The extension stiffness A equals: 1.5787E+09 5.1047E

32、+08 -2.2981E-08 5.1047E+08 1.5787E+09 -5.3951E-08 -2.2981E-08 -5.4129E-08 5.3412E+08 The coupling stiffness B equals: -1.0897E-09 2.1876E-10 -1.7357E-10 2.1876E-10 4.2280E-10 -4.4705E-10 -1.7357E-10 -4.5370E-10 3.7036E-10 The bending stiffness D equals: 4.2625E+04 1.3783E+04 -7.3708E-13 1.3783E+04 4

33、.2625E+04 -2.3412E-12 -7.3708E-13 -2.8932E-12 1.4421E+04 The NORMALISED stiffness tensor of the laminate: The NORMALISED extension stiffness A* equals: 8.7706E+10 2.8359E+10 -1.2767E-06 2.8359E+10 8.7706E+10 -2.9973E-06 -1.2767E-06 -3.0072E-06 2.9674E+10 The NORMALISED coupling tensor C equals: -6.9

34、290E-07 5.3458E-06 2.3989E-07 5.3458E-06 1.0632E-05 1.8201E-06 2.3989E-07 2.9460E-06 1.0096E-05 The NORMALISED bending stiffness D* equals: 8.7706E+10 2.8359E+10 -1.5166E-06 2.8359E+10 8.7706E+10 -4.8174E-06 -1.5166E-06 -5.9532E-06 2.9674E+10附件1：計(jì)算所用的程序代碼。PROGRAM COMPOSITE!Coded by James WilsonIMPLI

35、CIT NONEREAL(8):A(3,3),B(3,3),D(3,3),MC(5),TEMP,ROT(3,3)!A拉伸剛度；B耦合剛度；D彎曲剛度；!MC讀入材料常數(shù)；ROT旋轉(zhuǎn)矩陣REAL(8):TOTAL_TH,HALF_TH !總厚度；半厚度REAL(8),ALLOCATABLE:Q(:,:,:),AL(:),T(:),Z(:),Z1(:),Z2(:),Z3(:)!Q每層板相應(yīng)剛度；AL轉(zhuǎn)角；T每層板的厚度；Z坐標(biāo)量INTEGER(4):N,I,J,K,SEQ,L!_IJK循環(huán)變量；N板的層數(shù)；SEQ序數(shù)CHARACTER(50):CHR(8),TEMPC,filename1,file

36、name2!CHR、TEMPC:character variablesWRITE(*,*)Please insert the INP file name(a.txt for example):READ(*,*)filename1OPEN(8,file=filename1)!Open data file!Read in dataREAD(8,*)TEMPC,NALLOCATE(Q(3,3,N),AL(N),T(N),Z(N+1),Z1(N),Z2(N),Z3(N)READ(8,*)CHR(1:8)DO I=1,NREAD(8,*)SEQ,T(I),MC(1:5),AL(I)Q(:,:,I)=0!

37、Calculate stiffness of each layer for the principal axisTEMP=1./(1-MC(3)*MC(4)Q(1,1,I)=MC(1)*TEMPQ(2,2,I)=MC(2)*TEMPQ(3,3,I)=MC(5)Q(1,2,I)=MC(4)*MC(2)*TEMPQ(2,1,I)=Q(1,2,I)AL(I)=AL(I)*3.1415926535898/180ROT(1,1)=(cos(AL(I)*2!Work out Rot MatrixROT(2,2)=ROT(1,1)ROT(3,3)=cos(2*AL(I)ROT(2,1)=1-ROT(1,1)

38、ROT(1,2)=ROT(2,1)ROT(3,1)=0.5*sin(2*AL(I)ROT(3,2)=-ROT(3,1)ROT(1,3)=-2*ROT(3,1)ROT(2,3)=-2*ROT(3,2)Q(:,:,I)=MATMUL(MATMUL(ROT,Q(:,:,I),TRANSPOSE(ROT)ENDDOTOTAL_TH=sum(T)HALF_TH=TOTAL_TH/2Z(1)=-HALF_TH!Work out ZDO I=1,NZ(I+1)=Z(I)+T(I)ENDDO!calculate tensor A、B and DDO K=1,NZ1(K)=(Z(K+1)-Z(K)Z2(K)=(

39、Z(K+1)-Z(K)*(Z(K+1)+Z(K)/2Z3(K)=(Z(K+1)*3-Z(K)*3)/3ENDDOA=0;B=0;D=0WRITE(*,*)Please insert the OUP file name(b.txt for example):READ(*,*)filename2OPEN(9,file=filename2)!Write in stiffness tensor for each single plyDO K=1,NWRITE(9,100)K100 FORMAT(The stiffness of number,1X,I2,2X,ply is:)DO I=1,3WRITE

40、(9,200)Q(I,:,K)200 FORMAT(ES12.4,6X,ES12.4,6X,ES12.4)ENDDOWRITE(9,(/)A=A+Q(:,:,K)*Z1(K)B=B+Q(:,:,K)*Z2(K)D=D+Q(:,:,K)*Z3(K)ENDDO!Output the actual stiffness of the laminateWRITE(9,(/);WRITE(9,(/)WRITE(9,*)The ACTUAL stiffness tensor of the laminate:WRITE(9,(/)WRITE(9,*)The extension stiffness A equa

41、ls:DO I=1,3WRITE(9,200)A(I,1:3)ENDDOWRITE(9,(/)WRITE(9,*)The coupling stiffness B equals:DO I=1,3WRITE(9,200)B(I,1:3)ENDDOWRITE(9,(/)WRITE(9,*)The bending stiffness D equals:DO I=1,3WRITE(9,200)D(I,1:3)ENDDO!Normalised tensor outputWRITE(9,(/);WRITE(9,(/)WRITE(9,*)The NORMALISED stiffness tensor of

42、the laminate:WRITE(9,(/)WRITE(9,*)The NORMALISED extension stiffness A* equals:DO I=1,3WRITE(9,200)A(I,1:3)/TOTAL_THENDDOWRITE(9,(/)WRITE(9,*)The NORMALISED coupling tensor C equals:DO I=1,3WRITE(9,200)A(I,1:3)/TOTAL_TH-12*D(I,1:3)/TOTAL_TH*3ENDDOWRITE(9,(/)WRITE(9,*)The NORMALISED bending stiffness

43、 D* equals:DO I=1,3WRITE(9,200)12*D(I,1:3)/TOTAL_TH*3ENDDOWRITE(*,*)OUTPUT successfully,please press any key to end program!READ(*,*)END PROGRAM COMPOSITE附2楊濤同學(xué)(tng xu)的MATLAB（GUI）計(jì)算程序。主要程序：（編了個(gè)小界面，程序略長(zhǎng)，刪掉一些程序自帶解釋(jish)語(yǔ)句，添加了一些對(duì)關(guān)鍵語(yǔ)句的解釋。）界面是：作的一個(gè)(y )算例如下：該算例結(jié)果(ji gu)與組內(nèi)同伴James Wilson同學(xué)基本一致，其余算例結(jié)果(ji

44、gu)也基本一致，僅僅在趨近于零時(shí)有略微差異，在此不贅于。 后邊(hu bian)附上源代碼：function varargout = composit_plate(varargin)gui_Singleton = 1;gui_State = struct(gui_Name, mfilename, . gui_Singleton, gui_Singleton, . gui_OpeningFcn, composit_plate_OpeningFcn, . gui_OutputFcn, composit_plate_OutputFcn, . gui_LayoutFcn, , . gui_Callb

45、ack, );if nargin & ischar(varargin1) gui_State.gui_Callback = str2func(varargin1);end if nargout varargout1:nargout = gui_mainfcn(gui_State, varargin:);else gui_mainfcn(gui_State, varargin:);endfunction composit_plate_OpeningFcn(hObject, eventdata, handles, varargin)handles.output = hObject;guidata(

46、hObject, handles);ha=axes(units,normalized,position,0 0 1 1);%嵌入坐標(biāo)，為嵌入背景圖片準(zhǔn)備uistack(ha,down)%作為背景II=imread(武漢大學(xué).jpg);%讀入圖片信息image(II)colormap hsvset(ha,handlevisibility,off,visible,off)function varargout = composit_plate_OutputFcn(hObject, eventdata, handles) varargout1 = handles.output;function edi

47、t1_Callback(hObject, eventdata, handles)function edit1_CreateFcn(hObject, eventdata, handles)if ispc & isequal(get(hObject,BackgroundColor), get(0,defaultUicontrolBackgroundColor) set(hObject,BackgroundColor,white);endfunction edit2_Callback(hObject, eventdata, handles)function edit2_CreateFcn(hObje

48、ct, eventdata, handles)if ispc & isequal(get(hObject,BackgroundColor), get(0,defaultUicontrolBackgroundColor) set(hObject,BackgroundColor,white);endfunction edit3_Callback(hObject, eventdata, handles)function edit3_CreateFcn(hObject, eventdata, handles)if ispc & isequal(get(hObject,BackgroundColor),

49、 get(0,defaultUicontrolBackgroundColor) set(hObject,BackgroundColor,white);endfunction edit4_Callback(hObject, eventdata, )function edit4_CreateFcn(hObject, eventdata, handles)if ispc & isequal(get(hObject,BackgroundColor), get(0,defaultUicontrolBackgroundColor) set(hObject,BackgroundColor,white);en

50、d function pushbutton1_Callback(hObject, eventdata, handles)syms e1 e2 v21 g12 ang %本程序采用符號(hào)(fho)運(yùn)算v12=v21*e2/e1;q=e1/(1-v12*v21),v21*e2/(1- v12*v21),0v21*e2/(1-v12*v21),e2/(1-v12*v21),00,0,g12;tran= cos(ang)2, sin(ang)2, -sin(2*ang)sin(ang)2, cos(ang)2, sin(2*ang)sin(2*ang)/2, -sin(2*ang)/2, cos(2*a

51、ng);q1=tran*q*tran;%得到(d do)Qn=str2num(get(handles.edit2,string);%讀入層數(shù)n=floor(n);nn=0;A=0;B=0;D=0;t=str2num(get(handles.edit3,string);%讀入每層厚度(hud) t1=zeros(1,n+1);t1(1)=0; for nn=1:n t1(nn+1)=t1(nn)+t(nn); endang1=str2num(get(handles.edit4,string);%讀入每層角度const=str2num(get(handles.edit1,string);%讀入材料

52、系數(shù)t0=sum(t)/2;t1=t1-t0;e1=const(1);e2=const(2);v21=const(3);g12=const(4);q11=subs(q1);for nn=1:nang=ang1(nn);nn=nn+1;A=A+subs(q11*(t1(nn)-t1(nn-1);B=B+subs(0.5*q11*(t1(nn)2-t1(nn-1)2);D=D+subs(1/3*q11*(t1(nn)3-t1(nn-1)3);end %累加計(jì)算set(handles.edit5,string,num2str(A(1,:);%以下為輸出結(jié)果set(handles.edit6,stri

53、ng,num2str(A(2,:);set(handles.edit7,string,num2str(A(3,:);set(handles.edit8,string,num2str(B(1,:);set(handles.edit9,string,num2str(B(2,:);set(handles.edit10,string,num2str(B(3,:);set(handles.edit11,string,num2str(D(1,:);set(handles.edit12,string,num2str(D(2,:);set(handles.edit13,string,num2str(D(3,:

54、);function edit5_Callback(hObject, eventdata, handles)function edit5_CreateFcn(hObject, eventdata, handles)if ispc & isequal(get(hObject,BackgroundColor), get(0,defaultUicontrolBackgroundColor) set(hObject,BackgroundColor,white);endfunction edit6_Callback(hObject, eventdata, handles)function edit6_CreateFcn(hObject, eventdata, handles)if ispc & isequal(get(hObject,BackgroundColor), get(0,defaultUicontrolBackgroundColor) set