07弹性平面问题程序设计研究.pptx

弹性平面问题程序设计;一、弹性平面问题的程序结构 [enter] 二、主控变量和数据描述 [enter] 三、弹性平面问题的功能函数 [enter] 四、弹性平面问题的计算程序 [enter] 五、弹性平面问题的数值算例 [enter]; 一、弹性平面问题的程序结构;二、主控变量和数据描述;以图示弹性平面应力问题的有限元模型为例，假设板厚均匀、上半部分和下半部分材料不同，其主控变量： 结点总数 N=9； 已知结点位移分量总数 LR=4； 三角形单元总数 MT=8； 单元的材料特征类型数 LTB=2； 外载荷总数 LPQ=4； 应力状态信息 XEV=0。;弹性平面问题的数据描述;2. 已知结点位移数组;3. 单元信息数组;4. 材料特征数组;5. 载荷数组;三、弹性平面问题的功能函数;根据上图所示单元刚度矩阵函数的程序结构图，并利用单元子方阵函数krst，用MATLAB语言编制平面3结点三角形单元刚度矩阵函数dan_gang。;%----element1-1,2,5----------- xye1=[0,150;0,75;100,75]; e1=1;t1=1;mu1=0.3; ke1=dan_gang(xye1,e1,t1,mu1) %----element2-2,3,6----------- xye2=[0,75;0,0;100,0]; e2=1;t2=1;mu2=0.3; ke2=dan_gang(xye2,e2,t2,mu2);%----element2-2,3,6----------- xye2=[0,75;0,0;100,0]; e2=1;t2=1;mu2=0.3; ke2=dan_gang(xye2,e2,t2,mu2) %----element3-5,4,1----------- xye3=[100,75;100,150;0,150]; e3=1;t3=1;mu3=0.3; ke3=dan_gang(xye3,e3,t3,mu3) %----element3-1,5,4,1----------- xye3=[0,150;100,75;100,150]; e3=1;t3=1;mu3=0.3; ke3=dan_gang(xye3,e3,t3,mu3);三、弹性平面问题的功能函数;位移约束函数;结点位移函数;单元应力函数;求得单元应力后，各结点应力采用绕结点单元的平均应力计算得到。计算过程可由如下结点应力函数node_stress来完成。;四、弹性平面问题的计算程序;前处理程序;后处理程序;五、弹性平面问题的数值算例;clear; clc; XY=[1, 0, 150e-2 2, 0, 75e-2 3, 0, 0 4, 100e-2, 150e-2 5, 100e-2, 75e-2 6, 100e-2, 0 7, 200e-2, 150e-2 8, 200e-2, 75e-2 9, 200e-2, 0]; EMUTW=[180e9, 0.35, 5e-2, 0, 0 150e9, 0.40, 10e-2, 0, 0]; ELT=[1, 1, 2, 5, 1 2, 2, 3, 6, 2 3, 1, 5, 4, 1 4, 2, 6, 5, 2 5, 4, 5, 8, 1 6, 5, 6, 9, 2 7, 4, 8, 7, 1 8, 5, 9, 8, 2]; ; 结点号 水平位移 竖直位;五、弹性平面问题的数值算??;运行程序FEM，得到如下结果：;再 见 !