有限元杆单元讲解.pptVIP

第 2 章 杆单元与梁单元 第 2 章 杆单元与梁单元 第 2 章 杆单元与梁单元 § 2.1 等截面杆单元 2.1.1 一维等截面杆单元 杆单元 2.1.2 二维空间杆单元 ? 如何用直接法求杆单元特性？ ? 如何用公式法导出杆单元特性？ ? 什么是虚功原理？ ? 杆单元刚度矩阵的特点？ ? 什么叫坐标变换？ ? 如何对节点位移向量进行坐标变换？ ? 如何对刚度矩阵进行坐标变换？ ? 应用举例 第 2 章 杆单元与梁单元 § 2.1.1 一维等截面杆单元 ? ? ? ? ? ? ? j i u u d 单元节点位移： L — 杆长 A — 截面积 E — 弹性模量 考虑一个 2 节点一维等截面杆单元： ? ? ? ? ? ? ? j i f f f 单元节点力： 第 2 章 杆单元与梁单元 § 2.1.1 一维等截面杆单元 应力 — 应变关系： dx du ? ? ? ? E ? ) ( ) ( ) ( x x x u u ? ? ? ? ? ? ? —— 杆单元位移 —— 杆单元应变 —— 杆单元应力 应变 — 位移关系： 第 2 章 杆单元与梁单元 ? 应变： L ? ? ? ? 应力： L E E ? ? ? ? ? ? 杆内力： ? ? ? ? ? ? ? ? k L EA L EA A F ? ? 杆的轴向刚度： L EA k ? （一）直接法导出单元特性 杆单元伸长量 ： i j u u ? ? ? § 2.1.1 一维等截面杆单元 第 2 章 杆单元与梁单元 ? 轴向拉压变形模式下，该杆单元的行为与弹簧单元相同，因 此杆单元的刚度矩阵为： ? 比照弹簧元的刚度方程，写出杆单元的刚度方程为： ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? j i j i j i u u L EA u u k k k k f f 1 1 1 1 L EA k ? § 2.1.1 一维等截面杆单元 第 2 章 杆单元与梁单元 （二）公式法导出杆单元特性 ? 单元上假设近似位移函数 —— 位移模式 单元上位移假设为简单多项式函数： x a a u 1 0 ? ? 用插值法把多项式中的待定系数 转化为节点位移 , 从而得到插值形式的假设位移函数 —— 单元位移模式如下： 0 1 , a a ( ) ( ) ( ) i i j j u x N x u N x u ? ? ( ) 1 , ( ) i j x x N x N x L L ? ? ? 上式中： 形函数。 中称为形状函数，简称 是插值基函数，有限元 j i N N , § 2.1.1 一维等截面杆单元 第 2 章 杆单元与梁单元 单元位移模式写成矩阵形式： ? ? Nd ? ? ? ? ? ? ? ? j i j i u u N N u ( ) ( ) ( ) i i j j u x N x u N x u ? ? 注意：采用一次多项式是因为单元只有 2 个轴向位 移分量，对应 2 个多项式系数。 ? ? 称为单元形函数矩阵。 , j i N N ? N 式中 § 2.1.1 一维等截面杆单元 第 2 章 杆单元与梁单元 ? 单元应变： Bd d N ? ? ? ? ? ? ? ? ? dx d dx du ? —— 单元应变矩阵 B ? 单元应力： Bd E E ? ? ? ? ? 应用弹性体虚功原理导出单元刚度方程。 ? ? Nd ? ? ? ? ? ? ? ? j i j i u u N N u ? ? ( ) ( ) 1/ 1/ i j d N x N x L L dx ? ? ? ? ? ? ? B ? 虚功原理 § 2.1.1 一维等截面杆单元 第 2 章 杆单元与梁单元 弹性体受力平衡时，若发生 虚位移 ，则外力虚功等 于弹性体内的虚应变能。 —— 平衡条件 对于杆单元，定义虚位移如下： ? ? ? ? ? ? ? j i u u ? ? ? d 节点虚位移： 单元虚位移： d N ? ? ? u 节点力（外力）虚功： f d T ? 则单元虚应变： d B ? ? ?? ? ? ) ( u dx d § 2.1.1