三角形有限元matlab,3结点三角形单元有限元程序MATLAB语言

3结点三角形单元有限元程序(MATLAB语言)

学号：2011100290 吴晴晴

该程序包括以下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; %单元编号单元1 3-1-2;单元2 1-3-4

1 3 4];

NN=4; %结点数

node=[0 0; %各结点坐标

2 0;

2 1;

0 1];

RN=2; %被约束的位移数

RC=[1 4]; %被约束的结点

PN=2; %有载荷的结点数

%PC(1)表示有载荷的结点，PC(2)表示各结点的力，PC(3)表示载荷方向，0为x方向，1为y方向

PC=[2 3;

-1/2 -1/2;

1 1]; %结点2、3分别有y负方向上-1/2N的力作用