1、4 1.00E-03 9.60E+10 2.40E+10 0.10 0.40 1.00E+10 90.005 1.00E-03 9.60E+10 2.40E+10 0.10 0.40 1.00E+10 0.00（2）正交铺设6层反对称层合板（T5-8）层数 66 1.00E-03 9.60E+10 2.40E+10 0.10 0.40 1.00E+10 90.00（3）5层对称角铺设层合板（T5-9）1 1.00E-03 9.60E+10 2.40E+10 0.10 0.40 1.00E+10 30.002 1.00E-03 9.60E+10 2.40E+10 0.10 0.40 1.00E+
2、10 -30.003 1.00E-03 9.60E+10 2.40E+10 0.10 0.40 1.00E+10 30.004 1.00E-03 9.60E+10 2.40E+10 0.10 0.40 1.00E+10 -30.005 1.00E-03 9.60E+10 2.40E+10 0.10 0.40 1.00E+10 30.00（4）6层反对称角铺设层合板（T5-10）6 1.00E-03 9.60E+10 2.40E+10 0.10 0.40 1.00E+10 -30.00（5）我还想验证一个书上的例题，在课本P114。三层层合板，外层厚度t1，内层10t1，正交铺设比m=0.2,。
3、玻璃/环氧单层板性能：E1=5.4E10Pa，E2=1.8E10Pa，v21=0.25，G12=8.8E9Pa。层数 31 1.00E+00 5.40E+10 1.80E+10 0.083 0.250 8.80E+09 0.002 1.00E+01 5.40E+10 1.80E+10 0.083 0.250 8.80E+09 90.003 1.00E+00 5.40E+10 1.80E+10 0.083 0.250 8.80E+09 0.00二、验证Verchery的论文里给出的数值算例。这里一直到Table5的数据都是从Verchery的论文中截获。Verchery论文中的18层序列，第（2
4、1）式【laminates without bending-extension coupling】的排列有两种材料，一种是Boron-Epoxy，一种是Glass-Epoxy。而且都给出了最终的计算结果Q，A*，D*。下面是我的Fortran计算数据文档和结果文档。（1）Boron-Epoxy材料。（Boron-Epoxy）数据文档层数 18层序数 厚度m E1（Pa） E2（Pa） v12 v21 G12（Pa） 角度（1 1.00E-03 2.04E+11 1.85E+10 0.021 0.230 5.59E+09 0.002 1.00E-03 2.04E+11 1.85E+10 0.02
5、1 0.230 5.59E+09 0.003 1.00E-03 2.04E+11 1.85E+10 0.021 0.230 5.59E+09 60.004 1.00E-03 2.04E+11 1.85E+10 0.021 0.230 5.59E+09 60.005 1.00E-03 2.04E+11 1.85E+10 0.021 0.230 5.59E+09 -60.006 1.00E-03 2.04E+11 1.85E+10 0.021 0.230 5.59E+09 -60.007 1.00E-03 2.04E+11 1.85E+10 0.021 0.230 5.59E+09 -60.008
6、 1.00E-03 2.04E+11 1.85E+10 0.021 0.230 5.59E+09 60.009 1.00E-03 2.04E+11 1.85E+10 0.021 0.230 5.59E+09 0.0010 1.00E-03 2.04E+11 1.85E+10 0.021 0.230 5.59E+09 -60.0011 1.00E-03 2.04E+11 1.85E+10 0.021 0.230 5.59E+09 -60.0012 1.00E-03 2.04E+11 1.85E+10 0.021 0.230 5.59E+09 60.0013 1.00E-03 2.04E+11 1
7、.85E+10 0.021 0.230 5.59E+09 60.0014 1.00E-03 2.04E+11 1.85E+10 0.021 0.230 5.59E+09 0.0015 1.00E-03 2.04E+11 1.85E+10 0.021 0.230 5.59E+09 0.0016 1.00E-03 2.04E+11 1.85E+10 0.021 0.230 5.59E+09 0.0017 1.00E-03 2.04E+11 1.85E+10 0.021 0.230 5.59E+09 60.0018 1.00E-03 2.04E+11 1.85E+10 0.021 0.230 5.5
8、9E+09 -60.00（2）Glass-Epoxy材料。（Glass-Epoxy）数据文档层数18层序数厚度m E1（Pa）E2（Pa）v12v21G12（Pa）角度（11.00E-033.86E+108.27E+090.0560.2604.14E+090.0021.00E-033.86E+108.27E+090.0560.2604.14E+090.0031.00E-033.86E+108.27E+090.0560.2604.14E+0960.0041.00E-033.86E+108.27E+090.0560.2604.14E+0960.0051.00E-033.86E+108.27E+09
9、0.0560.2604.14E+09-60.0061.00E-033.86E+108.27E+090.0560.2604.14E+09-60.0071.00E-033.86E+108.27E+090.0560.2604.14E+09-60.0081.00E-033.86E+108.27E+090.0560.2604.14E+0960.0091.00E-033.86E+108.27E+090.0560.2604.14E+090.00101.00E-033.86E+108.27E+090.0560.2604.14E+09-60.00111.00E-033.86E+108.27E+090.0560.
10、2604.14E+09-60.00121.00E-033.86E+108.27E+090.0560.2604.14E+0960.00131.00E-033.86E+108.27E+090.0560.2604.14E+0960.00141.00E-033.86E+108.27E+090.0560.2604.14E+090.00151.00E-033.86E+108.27E+090.0560.2604.14E+090.00161.00E-033.86E+108.27E+090.0560.2604.14E+090.00171.00E-033.86E+108.27E+090.0560.2604.14E
11、+0960.00181.00E-033.86E+108.27E+090.0560.2604.14E+09-60.00（3）当然我也验证了第（22）【laminates with equal elastic properties in bending and extension】、（23）【quasi-homogeneous laminates】的排序，材料是Boron-Epoxy，下面给出计算的结果。从下面的两个结果表中可以知道，（22）排序的确是C=0，（23）的排序的确是B=0且C=0。验证成功。（A）第（22）排序。（B）第（23）排序。附件1：计算所用的程序代码。PROGRAM COM
12、POSITE!Coded by James WilsonIMPLICIT NONEREAL（8）:A（3,3）,B（3,3）,D（3,3）,MC（5）,TEMP,ROT（3,3）A拉伸刚度；B耦合刚度；D弯曲刚度；MC读入材料常数；ROT旋转矩阵TOTAL_TH,HALF_TH !总厚度；半厚度REAL（8）,ALLOCATABLE:Q（:,:）,AL（:）,T（:）,Z（:）,Z1（:）,Z2（:）,Z3（:Q每层板相应刚度；AL转角；T每层板的厚度；Z坐标量INTEGER（4）:N,I,J,K,SEQ,L_IJK循环变量；N板的层数；SEQ序数CHARACTER（50）:CHR（8）,TE
13、MPC,filename1,filename2CHR、TEMPC:character variablesWRITE（*,*）Please insert the INP file name（a.txt for example）:READ（*,*）filename1OPEN（8,file=filename1）!Open data fileRead 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）,
14、AL（I）,I）=0!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）
15、=1-ROT（1,1）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）,I）=MATMUL（MATMUL（ROT,Q（:,I）,TRANSPOSE（ROT）ENDDOTOTAL_TH=sum（T）HALF_TH=TOTAL_TH/2Z（1）=-HALF_THWork out ZZ（I+1）=Z（I）+T（I）calculate tensor A、B and DDO K=1,NZ1（K）=（Z（K+1）-Z（K）Z2（K）=（Z（K+1）-Z（K）
16、*（Z（K+1）+Z（K）/2Z3（K）=（Z（K+1）*3-Z（K）*3）/3A=0;B=0;D=0Please insert the OUP file name（b.txt for example）:READ（*,*）filename2OPEN（9,file=filename2）Write in stiffness tensor for each single plyWRITE（9,100）K100 FORMAT（The stiffness of number,1X,I2,2X,ply is:DO I=1,3WRITE（9,200）Q（I,:,K）200 FORMAT（ES12.4,6X,
17、ES12.4,6X,ES12.4）WRITE（9,（/）A=A+Q（:,K）*Z1（K）B=B+Q（:,K）*Z2（K）D=D+Q（:,K）*Z3（K）Output the actual stiffness of the laminate）;WRITE（9,*）The ACTUAL stiffness tensor of the laminate:The extension stiffness A equals:WRITE（9,200）A（I,1:3）The coupling stiffness B equals:WRITE（9,200）B（I,1:The bending stiffness
18、D equals:WRITE（9,200）D（I,1:Normalised tensor outputThe NORMALISED stiffness tensor of the laminate:The NORMALISED extension stiffness A* equals:3）/TOTAL_THThe NORMALISED coupling tensor C equals:3）/TOTAL_TH-12*D（I,1:3）/TOTAL_TH*3The NORMALISED bending stiffness D* equals:WRITE（9,200）12*D（I,1:OUTPUT
19、successfully,please press any key to end program!READ（*,*）END PROGRAM COMPOSITE附2杨涛同学的MATLAB（GUI）计算程序。主要程序：（编了个小界面，程序略长，删掉一些程序自带解释语句，添加了一些对关键语句的解释。）界面是：作的一个算例如下：该算例结果与组内同伴James Wilson同学基本一致，其余算例结果也基本一致，仅仅在趋近于零时有略微差异，在此不赘于。后边附上源代码：functionvarargout = composit_plate（varargin）gui_Singleton = 1;gui_Stat
20、e = struct（gui_Name, mfilename,.gui_Singleton, gui_Singleton,.gui_OpeningFcn, composit_plate_OpeningFcn,.gui_OutputFcn, composit_plate_OutputFcn,.gui_LayoutFcn, ,.gui_Callback, ）;ifnargin & ischar（varargin1）gui_State.gui_Callback = str2func（varargin1）;endifnargoutvarargout1:nargout = gui_mainfcn（gui
21、_State, varargin:）;elsegui_mainfcn（gui_State, varargin:functioncomposit_plate_OpeningFcn（hObject, eventdata, handles, varargin）handles.output = hObject;guidata（hObject, handles）;ha=axes（units,normalizedposition,0 0 1 1）;%嵌入坐标，为嵌入背景图片准备uistack（ha,down）%作为背景II=imread（武汉大学.jpg%读入图片信息image（II）colormaphsvset（ha,handlevisib
