复合材料力学上机编程作业(计算层合板刚度).docxVIP

复合材料力学上机编程作业学院：School of Civil Engineering 专业：Engineering Mechanics小组成员信息：James Wilson（2012031890015）、Tau Young（2012031890011） 复合材料力学学了五个星期，这是这门课的第一次编程作业。我和杨涛结成一个小组，我用的是Fortran编制的程序，Tau Young用的是matlab编制。其中的算例以我的Fortran计算结果为准。Matlab作为可视化界面有其独到之处，在附录2中将会有所展示。 作业的内容是层合板的刚度的计算和验算，包括拉伸刚度A、弯曲刚度D以及耦合刚度B。首先要给定层合板的各个参数，具体有：层合板的层数N；各单层的弹性常数E1、E2、、G12；各单层对应的厚度；各单层对应的主方向夹角。然后就要计算每个单层板的二维刚度矩阵Q，具体公式如下：；；；；得到Q矩阵后，根据课本上讲到的得到。 然后根据z坐标的定义求出到，接下来，最重要的一步，根据下式计算A、B、D。一、书上P110的几个问题可以归纳为以下几个类型。（1）正交铺设5层对称层合板（T5-7）数据文档层数5 层序数厚度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.00E+100.00 结果文档 拉伸刚度A= 3.5000E+08 5.0000E+07 -4.3711E-01 5.0000E+07 2.7500E+08 -6.1196E+00 -4.3711E-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显然是等于零的。正交铺设6层反对称层合板（T5-8）数据文档层数6 层序数厚度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.00E+100.00 6 1.00E-039.60E+102.40E+100.10 0.40 1.00E+1090.00 结果文档 拉伸刚度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.1250E+03 1.8000E+02 -1.967