有限元方法与MATLAB程序设计 第3章 桁架和刚架.pptxVIP

第3章桁架和刚架§3.1 平面桁架§3.2 平面刚架 §3.1.4 求解结点位移未知，结点位移已知，结点力已知。结点力未知。①②③④11 大型线性方程组求解方法12 大型线性方程组求解方法13 例 3.1结点、单元编号，建立结构和局部坐标系；建立结点位移向量和结点力向量；①②③④14 3. 建立单元刚度矩阵单元①①②③④15 3. 建立单元刚度矩阵单元②①②③④16 3. 建立单元刚度矩阵单元③①②③④17 4. 形成总刚度矩阵和总刚度方程总刚度矩阵18 桁架的总刚度方程19 5 建立桁架的结构刚度方程6 计算结点位移20 7 计算杆端力①②③④21 7 计算杆端力①②③④22 8 计算支座反力23 9 校核整体平衡①②③④24 §3.1.5 程序设计function [ke,T]=TrusStif(xy,EA) dl=xy(2,:)-xy(1,:); l=sqrt(dl*dl); ca=dl(1)/l;sa=dl(2)/l;% 杆长% 角度T0=[ca,sa;-sa,ca]; T=[T0,zeros(2,2);zeros(2,2),T0];k0=[1,0;0,0;]; k1=[k0,-k0;-k0,k0];ke=EA/l*T*k1*T; end% 单元坐标刚度矩阵% 结构25 结点位移向量26结点号与自由度的对应关系自由度序数结点号结点i的位移在2i-1:2i位置如果单元的两个结点是i和j，它的位移在[2i-1:2i 2j-1:2j]位置 §3.1.5 程序设计①②③④function Truss gxy=[0,3;0,0;4,3;4,0];ndel=[1,3;2,3;4,3;2,4];EA=3e8*ones(1,4); nd=size(gxy,1); ne=size(ndel,1); F=zeros(2*nd,1); F([6,7])=[-2e4;3e4]; b=1:4;a=setdiff(1:2*nd,b);% 桁架% 结点坐标% 单元信息% 抗拉刚度% 结点总数% 单元总数% 置结点力% 约束自由度*% 无约束自由度27 刚度矩阵和求解K=sparse(2*nd,2*nd);阵for el=1:ne [ke,T]=TrusStif(gxy(ndel(el,:),:),EA(el));N(2:2:4)=2*ndel(el,:);N(1:2:4)=N(2:2:4)-1; K(N,N)=K(N,N)+ke;阵 endU=zeros(2*nd,1); U(b)=0;U(a)=K(a,a)\(F(a)-K(a,b)*U(b));% 结构总刚度矩% 单元刚度矩阵% 单元自由度% 安装总刚度矩%置约束位移28 29%for el=1:ne %输出单元号，单元信息，抗拉刚度和杆端力[ke,T]=TrusStif(gxy(ndel(el,:),:),EA(el));单元刚度矩阵 N(2:2:4)=2*ndel(el,:);N(1:2:4)=N(2:2:4)-1;单元自由度Fe=T*ke*U(N); % 杆端力fprintf(%4i%4i%4i%14.4g%14.4g\n ,el,ndel(el,:),EA(el),Fe(3));end % 输出单元号，单元信息，抗拉刚度和杆端力结果输出fprintf(%4s%7s%7s%12s%12s\n,结点,X坐标,Y坐标,u位移,v位移)for i=1:nd输出结点编号，坐标和位移 fprintf(%4i%10.4f%10.4f%14.4g%14.4g,i,gxy(i,:),U(2*i+(-1:0))) endfprintf(%4s%4s%12s%12s\n,单元,结点,面积,轴力)) 30运行结果结点X坐标Y坐标u位移v位移10.000030000000003.00000.0003556-0.001444.00000.00000.0004-0.0014单元结点面积轴力1133e+0082.667e+0042233e+008-3.333e+0043433e+0082.91e-0114243e+0083e+004 §3.2 平面刚架§3.2.1位移模式思考：这些形函数是近似地还是精确地描述变形？边界条件31 §3.2 平面刚架§3.2.1位移模式定义形函数矩阵32 形函数几何意义33 梁的弯曲应力ydxxdxxy34 §3.2.2 应变和应力几何矩阵几何方程35 应变矩阵（几何矩阵）36 位移应变应力单元结点位移列向量37 单元刚度矩阵虚应变虚位移应变38位移应力内力虚功刚架单元刚度矩阵 外力分布外力yxxy集中外力xyc39 杆端力杆端位移内力虚功虚功原理结点力与外力功位移虚位移 外力虚功等效结点力40 单元刚度矩阵41 §3.2.4 结构坐标系下单元刚度方