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

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; %单元编号 单元 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 的力作用在 matlab 环境下，输入则程序运行结果如下：该程序求解的结点位移结果和结点应力结果与课本给出的结果一致。附录：%%-------------平面三角形单元有限元法---------------------%%function [x strain stress]=tri_fem()[E,v,t,EN,ecode,NN,node,RN,RC,PN,PC]=input; %调入已定材料、几何尺寸以及单元和结点编号及约束和载荷分布[n m]=size(ecode);if EN~=n error('Wrong elementnumber EN or wrong elementcode ecode!'); return;end[n m]=size(node);if NN~=n error('Wrong nodenumber NN or wrong node-coordinate node!'); return;ende=zeros(EN,6); A=zeros(EN,1); %面积for i=1:EN e(i,:)=[node(ecode(i,1),:),node(ecode(i,2),:),node(ecode(i,3),:)]; %各三角形单元的节点坐标 D=[1,node(ecode(i,1),:) 1,node(ecode(i,2),:) 1,node(ecode(i,3),:)]; A(i,1)=1/2*det(D);end%% 形成单元参数 b=zeros(EN,3); c=zeros(EN,3); %各单元参数初始化for i=1:EN b(i,1)=e(i,4)-e(i,6); b(i,2)=e(i,6)-e(i,2); b(i,3)=e(i,2)-e(i,4); c(i,1)=-e(i,3)+e(i,5); c(i,2)=-e(i,5)+e(i,1); c(i,3)=-e(i,1)+e(i,3);end%% 求得总...