有限元的matlab编程资料.pptxVIP

有限元编程示例;题目描述:;解题思路：;建立模型：;定义单元，即储存单元两端的节点号;集成总刚：;计算从局部坐标到整体坐标的坐标转换矩阵T;集成整体刚度矩阵K;集成总刚的对角线元素（这里的元素指2*2的小矩阵）;求解位移：;求解轴力：;输出求解结果：;例二：网架;思路分析;用户自定义输入;e=input(选择网架类型，0代表自由定义网架，1代表四角锥网架) %网架类型的选择; if Cont1==0 ;正放四角锥网架定义;for i=1:hu for j=1:lu Node((i-1)*lu+j,2)=(j-1)*dis_xu; Node((i-1)*lu+j,3)=(i-1)*dis_yu; Node((i-1)*lu+j,4)=dis_z; end end;Nsum=length(Node); %查询网架的节点数; for i=1:hd for j=1:ld-1 Men((i-1)*(ld-1)+(lu-1)*hu+(hu-1)*lu+j,2)=(i-1)*ld+j+hu*lu; Men((i-1)*(ld-1)+(lu-1)*hu+(hu-1)*lu+j,3)=(i-1)*ld+j+hu*lu+1; end end;网架腹杆单元的拓扑;单元编号录入单元储存矩阵;AE(1:m1,2)=A1; %将上层单元尺寸录入AE矩阵 AE((m1+1):m2,2)=A2; %将下层单元尺寸录入AE矩阵 AE((m2+1):Msum,2)=A3; %将腹杆单元尺寸录入AE矩阵;定义边界条件;elseif cont3==1 n1=2*(hd+ld-2); BC=zeros(n1,4); BC(1:ld-2,1)=2:ld-1; BC((ld-1):(2*ld-4),1)=ld*(hd-1)+2:ld*hd-1; BC((2*ld-3):(2*ld-4+hd),1)=1:ld:ld*(hd-1)+1; BC((2*ld-3+hd):n1,1)=ld:ld:hd*ld; for i=1:n1 BC(i,1)=BC(i,1)+hu*lu; end;Nsum=length(Node);Msum=length(Men); Psum=length(P);BCsum=length(BC); %提取各矩阵的行数;求解网架杆件的应力应变;Lt=L1-L0; %所有杆件长度的增量;网架节点位移求解函数; %定义转化矩阵 t=sqrt(l^2+n^2); if t==0 r=[0 m 0;-m 0 0;0 0 1]; else r=[l m n;-l*m/t t -m*n/t;-n/t 0 l/t]; end O=zeros(3,3); T=[r O;O r];;%单刚矩阵的扩充，使之行数、列数与总刚对应 G=zeros(6,3*Nsum); I=[1 0 0;0 1 0;0 0 1]; G(1:3,3*p-2:3*p)=I; G(4:6,3*q-2:3*q)=I; K=G*k*G;;具体算例;选择网架类型，0代表自由定义网架，1代表四角锥网架：1 输入网架上层节点行数：4 输入网架上层节点列数：6 输入网架上层节点列间距：2 输入网架上层节点行间距：2 输入网架上下层间距：2 请输入网架上层单元的截面面积：0.0003 请输入网架下层单元的截面面积：0.0003 请输入网架腹杆单元的截面面积：0.0004 定义节点荷载，若网架上层节点力与下层节点力均布，则输入0，否则输入1：0 请输入网架上层节点荷载：10000 请输入网架下层节点荷载：0 定义边界约束，若网架上层周边节点全约束，则输入0，若下层周边节点全约束，输入1，否则输入2：0 ;作业：;提交作业要求：