有限元法理论及应用答案.doc

有限元法理论及应用大作业 试简要阐述有限元理论分析的基本步骤主要有哪些？ 答：有限元分析的主要步骤主要有： （1）结构的离散化，即单元的划分； （2）单元分析，包括选择位移模式、根据几何方程建立应变与位移的关系、根据虚功原理 建立节点力与节点位移的关系，最后得到单元刚度方程； （3）等效节点载荷计算； （4）整体分析，建立整体刚度方程； （5）引入约束，求解整体平衡方程。 2、有限元网格划分的基本原则是什么？指出图示网格划分中不合理的地方。 题2图 答：一般选用三角形或四边形单元，在满足一定精度情况，尽可能少一些单元。 有限元划分网格的基本原则: 拓扑正确性原则。即单元间是靠单元顶点、或单元边、或单元面连接 几何保持原则。即网络划分后，单元的集合为原结构近似 特性一致原则。即材料相同，厚度相同 单元形状优良原则。单元边、角相差尽可能小 密度可控原则。即在保证一定精度的前提下，网格尽可能的稀疏一些。 （a） (1). (2). 答：（1）不能，因为位移函数要满足几何各向同性，即单元的位移分布不应与人为选取的 坐标方位有关，即位移函数中的坐标x,y应该是能够互换的。所以位移多项式应按巴斯卡三角形来选择。 （2）不能，位移函数应该包括常数项和一次项。 6、设位移为线性变化，将图示各单元边上的载荷等效到相应的节点上去。 （1）集中力F平行于x轴，e点到i、j点的距离分别为ｌieｌjeｌijij边上有线性分布载荷，最大值为q。 题6图 答：（1） （2）i,j两节点受到的力分别为， 7、图示三角形ijm为等边三角形单元，边长为ｌF垂直作用于mj边的中点，集度为q的均布载荷垂直作用于im边。写出三角形单元的节点载荷向量。 题7图 题8图 答：将q移置到m,i节点： 将F移置到m,j两节点： 将重力移置到i,j,m点： 叠加后得： 8、如图所示为线性位移函数的三角形单元，若已知i、jij边上任意一点的位移都为零。 证：设ij边上任一点坐标为x,y，则其位移为： ∵i、j点位移为0 ∴所以ui,vi,uj,vj均为0 要证 {δ}=0，只需证 Nm=0 ∵Nm=(am+bm x +cmy)/2A ，am=xiyj-xjyi ，bm=yi-yj ，cm=xj-xi ∴Nm= [xiyj-xjyi+(yi-yj)x+(xj-xi)y]/2A=[xyi-yxi]/2A ∵该点为ij边上任一点 ∴yi/xi=y/x ∴Nm = 0 9、已知图示的三角形单元，其jm边和mi边边长均为a，单元厚度为t，弹性模量为E，泊松比为μ＝0，试求： （1）N； （2）应变矩阵B； （3）应力矩阵S； （4）K。 解：令m点为坐标原点，则m点坐标为（0,0），j点坐标为（0,a），i点坐标为（a,0） ，， ，，; ，，. ,, 题9图 题10图 10、如图所示，设桁架杆的长度为ｌA，材料弹性模量为E，单元的位移函数为u(x)=α1＋α2x 2点： x=l u=u2 [K] e=∫∫V[B]T[D][B]dv [D] -----为弹性矩阵(对于一维问题,为E) 11、如图为一悬臂梁，其厚度为1m，长度为2 m，高度为1 m，弹性模量为E，泊松比为μ＝1/3，在自由端面上作用有均匀载荷，合力为F，若用图示两个三角形单元进行有限元分析，试计算各个节点的位移；若将悬臂梁离散为四个平面三角形单元，令μ＝0，试求整体刚度矩阵。 解：离散为两个单元求各节点位移,假设t很小,则该问题为平面应力问题： 一、单元编号、节点坐标 单元号 节点号 ① ② i 1 2 j 2 3 m 4 4 各节点的坐标为：1(0,0),2(2,0),3(2,1),4(0,1) 面积A=1; 二、求单元刚度矩阵 （１）对单元① (i=1,j=2,m=4) 由ai=xjym-xmyj bi=yi-ym ci=xm-xj 得 b1= -1 c1=-2 b2=1 c2=0 b4=0 c4=2 由 r,s = i,j,m 令 得: （１）、对单元② (i=2,j=3,m=4) 同理求得：b2 = 0 c2 = -2 b3 = 1 c3 = 2 b4 = -1 c4=0 求得：