有限元与数值方法.docxVIP

1 加权残数法的基本思想及加权参数法与有限元方法的联系 基本思想：加权残数法是直接从应用科学中的控制微分方程式求得近似解的数学方法。联系：有限元法和加权残数法都是将微分方程转化为代数方程，但转化的方法不同。 2 编写一个用三角形单元分析弹性力学平面问题的程序 编写程序流程图： 阵 采用matlab编程 子函数 function [ v1 ] = fun( x1,x2,x3,y1,y2,y3,E) syms x; syms y; v = 0.3; B = zeros(3,6); zuobiao = [1 x1 y1;1 x2 y2;1 x3 y3]; area = det(zuobiao)/2; a1 = x2*y3-x3*y2; a2 = x3*y1-x1*y3; a3 = x1*y2-x2*y1; b1 = -(y3-y2); b2 = -(y1-y3); b3 = -(y2-y1); c1 = x3-x2; c2 = x1-x3; c3 = x2-x1; N(1) = 1/2/area*(a1+b1*x+c1*y); N(2) = 1/2/area*(a2+b2*x+c2*y); N(3) = 1/2/area*(a3+b3*x+c3*y); dif = [x,0;0,y;y,x]; N = [N(1),0,N(2),0,N(3),0;0,N(1),0,N(2),0,N(3)]; for i=1:3 for j=1:6 for k = 1:2 if dif(i,k)~=0 B(i,j) = B(i,j) + diff(N(k,j),dif(i,k)); end end end end D = [E/(1-v^2) v*E/(1-v^2) 0;v*E/(1-v^2) E/(1-v^2) 0;0 0 E/(1-v^2)*(1-v)/2]; K = B*D*B*area；%总刚度矩阵 end 主函数 clc; clear; E = 210e3;%弹性模量 t = 1;%壳单元厚度 x1=2;x2=0;x3=0;y1=0;y2=1;y3=0;%节点坐标 fun( x1,x2,x3,y1,y2,y3,E) 此函数返回矩阵总刚为×10e5 采用改0置1法，将约束条件引入方程 [K]×{u}={P} 其中{u}=[0 v1 0 0 0 0]T ，{P}=[0 -1000 0 0 0 0]T 得v1=0.04953mm 采用ansys进行计算，最终得到如图的云图，得节点1的位移为0.049376mm 与用matlab编写的程序求解一致，则证明程序编写的正确。 ansys命令流/CLEAR /PREP7 ******单元类型和材料参数****** ET,1,SHELL181 MP,EX,1,,210E3 MP,PRXY,1,,0.3 sect,1,shell,, secdata, 1,1,0.0,3 secoffset,MID seccontrol,,,, , , , *******建立有限元模型******** N,1,2,,,,,, N,2,0,1,,,,, N,3,0,0,,,,, FLST,2,3,1 FITEM,2,1 FITEM,2,3 FITEM,2,2 E,P51X FINISH ******施加约束和载荷****** /SOL FLST,2,2,1,ORDE,2 FITEM,2,2 FITEM,2,-3 D,P51X, , , , , ,ALL, , , , , FLST,2,1,1,ORDE,1 FITEM,2,1 F,P51X,FY,-1000 /STATUS,SOLU *****求解********* SOLVE FINISH /POST1 PLNSOL, U,Y, 0,1.0