结点三角形单元有限元程序MATLAB语言(共7页).doc

精选优质文档-倾情为你奉上3结点三角形单元有限元程序（MATLAB语言）学号： 吴晴晴该程序包括以下6个部分：1 主程序tri_fem：用于数据的录入和其他程序的调用；2 总刚程序Kf：计算结构的总体刚度；3 各结点位移求解程序xf：求解各结点的位移；4 线性方程组求解程序Jordan：Gauss-Jordan法求解非约束结点的位移；5 应力应变程序ss：由各结点位移求解各单元内的三个结点的应力stress和应变strain；6 数据录入程序input：录入材料、几何尺寸、单元编号和结点编号、位移约束和已知载荷等。以课本P25页例2.2为例，其input程序为function E,v,t,EN,ecode,NN,node,RN,RC,PN,PC=input()E=2.1e11; v=1/3; t=1; %杨氏模量Pa，泊松比，厚度EN=2; %单元数ecode=3 1 2; %单元