复合材料的MATLAB计算M文件

这是关于用MATLAB编写M文件求复合材料层合板性能指标的一个作业，第一个M文件是对层合板性能的分析，第二个M文件是对层合板破坏的分析。

%MOCMexp01.m
%MECHANICS OF COMPOSITE MATERIALS EX1
%已知：层合板[θ/-θ]s其单层板性能指标为E1=140GPaE2=10GPaG12=5GPau12=0.3tp=0.125mm
%求弹性系数（弯曲和拉伸两组）ExGxyuxymx与theta的关系（θ为0~90°，Δθ=5°）
%对结果进行列表，与单层板结果进行对比论证
%清屏，清除原有变量，程序开始
echo off;
clc;
clear;
format short;
%…………………………………………………………………………………………………………
%1、预处理，定义已知材料性能指标及相关参数
E1=140;
E2=10;
G12=5;
u12=0.3;
tk=[0.1250.1250.1250.125];  %对应1-4层单层板的厚度，单位mm
zk=[-0.1875-0.06250.06250.187];  %对应1-4层单层板的中面坐标，单位mm
t=0.5;  %对应层合板的总厚度，单位mm
%由单层板性能指标求刚度系数Q
Q=zeros（33）;
u21=u12*E2/E1;
Q（11）=E1/（1-u12*u21）;
Q（22）=E2/（1-u12*u21）;
Q（33）=G12;
Q（12）=u12*Q（22）;
Q（21）=Q（12）;
Q（13）=0;Q（31）=0;Q（23）=0;Q（32）=0;
%此处单位为GPa
%…………………………………………………………………………………………………………
%2、计算，引入变量θ（0~90°）以Δθ=5°的步进循环求值
for theta=0:5:90
    c=theta/5+1;  %c为计数变量
    %转换成变换刚度矩阵Q（θ）
  

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       3845  2010-05-10 21:08  MOCMexp02.m

     文件       4679  2010-05-10 23:10  MOCMexp01.m

----------- ---------  ---------- -----  ----

                 8524                    2