loading.
%生成弹性矩阵 D
D= [1 POISS 0; POISS 1 0;
0 0 (1-POISS)/2]*YOUNG/(1-POISS^2) %********************************************************** %计算当前单元的面积
A=-det([1 COORD(LNODS(i,1),1) COORD(LNODS(i,1),2); 1 COORD(LNODS(i,2),1) COORD(LNODS(i,2),2); 1 COORD(LNODS(i,3),1) COORD(LNODS(i,3),2)])/2
%********************************************************** %生成应变矩阵 B for j=0:2
b(j+1)=COORD(LNODS(i,(rem((j+1),3))+1),2)-COORD(LNODS(i,(rem((j+2),3))+1),2);
c(j+1)=-COORD(LNODS(i,(rem((j+1),3))+1),1)+COORD(LNODS(i,(rem((j+2),3))+1),1); end
B=[b(1) 0 b(2) 0 b(3) 0; 0 c(1) 0 c(2) 0 c(3);
c(1) b(1) c(2) b(2) c(3) b(3)]/(2*A); B1( :,:,i)=B;
%********************************************************** %求应力矩阵 S=D*B S=D*B;
ESTIF=B'*S*THICK*A; %求解单元刚度矩阵
a=LNODS(i,:); %临时向量,用来记录当前单元的节点编号 for j=1:3
for k=1:3
ASTIF((a(j)*2-1):a(j)*2,(a(k)*2-1):a(k)*2)=ASTIF((a(j)*2-1):a(j)*2,(a(k)*2-1):a(k)*2)+ESTIF(j*2-1:j*2,k*2-1:k*2); %根据节点编号对应关系将单元刚度分块叠加到总刚 %度矩阵中 end end end
%********************************************************** %将约束信息加入总体刚度矩阵（对角元素改一法）
. .
.
for i=1:NVFIX
if FIXED(i,2)==1
ASTIF(:,(FIXED(i,1)*2-1))=0; %一列为零 ASTIF((FIXED(i,1)*2-1),:)=0; %一行为零 ASTIF((FIXED(i,1)*2-1),(FIXED(i,1)*2-1))=1; %对角元素为 1 end
%********************************************************** %生成单元刚度矩阵并组成总体刚度矩阵
%**********************************************************
if FIXED(i,3)==1
ASTIF( :,FIXED(i,1)*2)=0; %一列为零 ASTIF(FIXED(i,1)*2,:)=0; %一行为零 ASTIF(FIXED(i,1)*2 ,FIXED(i,1)*2)=1; %对角元素为 1
end end
%********************************************************** %生成荷载向量
ASLOD(1:2*NPION)=0; %总体荷载向量置零 for i=1:NFORCE
ASLOD((FORCE(i,1)*2-1):FORCE(i,1)*2)=FORCE(i,2:3); end
%********************************************************** %求解内力
ASDISP=ASTIF\\ASLOD' %计算节点位移向量 ELEDISP(1:6)=0; %当前单元节点位移向量 for i=1:NELEM for j=1:3
ELEDISP(j*2-1:j*2)=ASDISP(LNODS(i,j)*2-1:LNODS(i,j)*2); %取出当前单元的节点位移向量 end i
STRESS=D*B1(:, :, i)*ELEDISP' %求内力 end
（程序计算结果和有限元软件得出的结果稍有偏差，可能是程序某些地方数据输入时出了问题，还在寻找具体原因）
. .
.
二、有限元软件分析
设置材料参数 建模
. .
.
. 边界约束 .
网格划分
添加载荷 .
欢迎您的光临，Word文档下载后可修改编辑双击可删除页眉页脚谢谢！希望您提出您宝贵的意见，你的意见是我进步的动力。赠语； 1、如果我们做与不做都会有人笑，如果做不好与做得好还会有人笑，那么我们索性就做得更好，来给人笑吧！ 2、现在你不玩命的学，以后命玩你。、我不知道年少轻狂，我只知道胜者为王。、不要做金钱、权利的奴隶；应学会做“金钱、权利”的主人。、什么时候离光明最近？那就是你觉得黑暗太黑的时候。、最值得欣赏的风景，是自己奋斗的足迹。、压力不是有人比你努力，而是那些比你牛×几倍的人依然比你努力。
. .