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 的力作用
在
评论1
最新资源