有限元理论题答案.docVIP

课 程 名 称 有限元单元法 学 院 机械工程学院 专 业 机械制造及其自动化 学 生 姓 名 周祥态 学 号 110010021 教 师 陈爱军 理论部分习题及解答： 1．已知平面应变下三节点三角形单元的节点坐标、和，单元的节点位移分量、，材料的弹性模量，泊松比为。试求单元的应变和应力分量？ 解：单元面积： 单元应变 平面应变下： 2．图示为一个三个节点的杆单元，为坐标原点，其位移模式取为。设其为常数，试求其单元的刚度矩阵。 解： 形函数： 3．已知图示正方形薄板的边长为，厚度为，弹性模量为，泊松比。现将其分成两个三角形单元，设节点2、3和4为不动点，在节点1处受到向上的集中载荷P。试求节点位移，支座反力以及单元A和单元B内的应力？ 解： A单元： B单元： 整体刚度矩阵： 单元A内的应力： 单元B内的应力： 4．已知集中载荷，试求图示六节点三角形单元的等效节点载荷列阵。 解： 在点（4，4）处，有，， ，， ，， 5．已知各向同性材料的常数和。计算图示四面体单元的形函数；单元的刚度矩阵中的元素和。 解：i(0,0,0) j(a,0,0) m(0,b,0) p(0,0,c) 6、图示矩形单元受到均匀重力载荷（每单位体积载荷）作用，试用阶 高斯积分计算节点2处的等效节点载荷。已知材料的弹性模量、泊松比和厚度。 解：令， 得到形函数： 7．写出线性动力学有限元方程，并说明方程中各个符号的力学含义。简介线性动力学有限元方程的求解方法，并且说明方法的思想和分析步骤。 答：其中，M、C、K分别为整体质量矩阵、阻尼矩阵和刚度矩阵；、、、分别为整体节点加速度、速度、位移和动载荷向量。 求解方法：为简便推导动力有限元方程式，可以先把动力问题化成静力问题，即将惯性力和阻尼力看作体积力施加在结构物上，这样就可以按照静力问题的有限元分析推到有关方程。问题的关键在于将惯性力和阻尼力转化成单元等效节点载荷。 （1）位移、速度、加速度 将结构离散以后，单元位移函数为： 其中单元节点位移向量为： ， 形函数矩阵为 I为三阶单位矩阵。 单元内部各点速度和加速度分别可用单元节点的速度和加速度表示为： （2）惯性力的等效节点载荷 表材料的质量密度为ρ，则根据D’Alembert原理，单元内单位体积上作用的惯性力为： 下标i表示惯性。 惯性力的等效节点载荷： （3）阻尼力的等效节点载荷 假设阻尼力正比于运动速度，比例常数为α，则单元内单位体积上作用的阻尼力可表示为： 下标d表示阻尼。 阻尼力的等效节点载荷为： （4）动力有限元方程的建立 若用P表示外部动载荷产生的整体节点载荷向量，则总的整体节点载荷向量为：。将虚位移原理应用于整个机构，可得：，即，。 系数矩阵： （1）质量矩阵M 常用方法：一致质量法和集中质量法。 （2）阻尼矩阵C 采用Royleigh阻尼：，其中α、β为不依赖于频率的常数。 Newmark方法的思想和分析步骤 假设： （a） （b） 其中，r和β是按积分精度和稳定性要求而决定的参数。 当r=1/2和β=1/2时，式（a）和（b）相应于线性加速度。 在Newmark法中，时刻的位移解答也是通过满足时刻的运动方法： (*)而得到的。 首先从式（b）中解得 代入式（a），然后再一并代入式（*），从而得到由、、计算的公式： 其中、分别为： 8．图示两个三角形单元的集合体，边长为，厚度为，材料的密度为，弹性模量为，泊松比，并且已知阻尼的常数和。写出单元集合体的整体一致质量矩阵和阻尼矩阵。 解：首先对局部节点进行编号 对单元①： 单元①的一致性质量矩阵： 单元①的刚度矩阵： 对单元②： 单元②的一致性质量矩阵： 单元②的刚度矩阵： 整体一致性质量矩阵为： 整体刚度矩阵为： Rayleigh阻尼矩阵为： 上机部分习题答案： a ② ① a x y 4 3 1 2 B A a a L/2 L/2 3 2 O 1 x y P 4 3 1 2 3 45o m(5,6) j(6,3) i(2,2) (4,4) x y P 1 2 z y x m p j i c a b b b