“有限元法基础及应用”讲义-2014.pdf

有限元法基础及应用 补充讲义 2014 年 2 月 “有限元法基础及应用”补充讲义 一、弹簧单元与弹簧系统 1、弹簧单元分析 1）单元描述 弹簧系统受力平衡时，从中隔离出一个典型弹簧单元进行分析。 图 1-1 弹簧单元的端点i, j 设置为单元节点。 基本未知量为节点位移：u ,u i j 单元节点力（单元在节点处受到的作用力）：f i , f j 已知弹簧的物理特性：F k ⋅∆ 其中：k为弹簧刚度 u u 为弹簧伸长量，F 为弹簧力（拉伸为正） , ∆ j − i 2）建立弹簧单元的单元刚度方程 考虑弹簧单元在弹簧系统中受力平衡和弹簧物理特性，得到下列方程： f −F −k (u −u ) ku −ku i j i i j （1-1） f F k (u −u ) −ku +ku j j i i j 写成矩阵形式： f u    k −k   i i （1-2 ）      f j −k k  uj     上式的矩阵符号形式为： f kd （1-3 ） 方程（1-2 ）或（1-3 ）称为弹簧单元的刚度方程，反映了单元的力学特性， 即节点力~节点位移之间的关系。 式 （1-3 ）中： 1  k −k  k  , 称为单元刚度矩阵  −k k  u   i d , 称为单元节点位移列阵   u  j  f   i f , 称为单元节点力列阵   f j   3）弹簧单元刚度方程的讨论 a. k有何特点？ 对称、奇异、主对角元素恒正。 b. k中元素的物理意义是什么？ 刚度矩阵元素的大小等于弹簧刚度。 从对方程（1-2 ）的分析可以看出，矩阵中某列的各元素代表列序号对应节 点有单位位移，其它节点位移为零时，单元各节点上的节点力；某行的各元素分 别是单元各节点位移对行序号对应节点的节点力的贡献系数。 因此，矩阵中任意一个元素k 的物理意义是： 节点的位移对 节点的节点 ij j i 力贡献系数，或者 节点有单位位移，其他节点位移为零时， 节点上的节点力。 j i c. 单元刚度方程可以求解吗？为什么？ 不可以。单元刚度方程仅仅表征一个单元的力学特性，单元水平上无法确定 单元节点位移。只有把系统中所有单元特性集成后，在系统水平上才可