弹性力学-06(gai).ppt

单元号 1 2 局部编码 整体编码 ui vi 2,0 0,0 uj vj 0,0 2,0 um vm 0,0 0,1 三、建立结构的整体结点荷载列阵 单元号 1 2 局部编码 整体编码 ui vi 2,0 0,0 uj vj 0,0 2,0 um vm 0,0 0,1 四、解平衡方程，求结点位移 得 ， 五、代回各单元，求得应力分布，应变分布。 §6-8 解题的具体步骤 单元的划分 求解结点位移的解题思路： 1．将计算对象划分成许多三角形单元，并按一定的规律将所有的结点 和单元分别编上号码； 2．选定一个直角坐标系，写出整体结点位移矩阵 3．写出整体结点荷载矩阵 4．算出整体劲度矩阵K； 5．列出整体平衡方程组 ，解出未知的结点位移分量 例2 在上题基础上，若 处为固定铰链支座， 处作用有水平力P和竖直力P。 解：因为： 且 其中： 又由于只有节点 上有力和位移，其它节点均为零，故只需求 ， 而由上题结果知： ? 所以 于是 * 第六章 用有限单元法解平面问题 第一节 基本量及基本方程的矩阵表示 第二节 有限单元法的概念 第三节 单元的位移模式与解答的收敛性 第四节 单元的应变列阵和应力列阵 第五节 单元的结点力列阵和劲度列阵 第六节 荷载向结点移置 单元的结点荷载列阵 第七节 结构的整体分析 结点的平衡方程 第八节 解题的具体步骤 单元的划分 第九节 计算成果的整理 第十节 计算实例 §6-1 基本量及基本方程的矩阵表示 1．体力列阵： 2．面力列阵： 3．应力列阵： 4．应变列阵： 5．位移列阵： 6．平面应力问题的物理方程： 称为弹性矩阵。 对于平面应变问题，只需 将弹性矩阵D 中的E、? 分别换成 即 可。 写成矩阵表示： 式中： 用 表示虚位移，用 表示与该虚位移相应的 虚应变，由虚功方程： 在有限单元法中，作用于弹性体的各种外力常以作用于结点的等效集中力来代替 i j m ui vi Fxi Fyi um vm Fym Fxm uj vj Fyj Fxj 集中力作用下的虚功方程可写成： 1、将连续体变换成为离散化结构。 §6-2 有限单元法的概念 一．有限单元法的概念 简单的说，就是用结构力学方法求解弹性力学问题。即首先将连续体变换成为 离散化结构，然后再用结构力学方法进行求解。 二．有限单元法的解题思路： 2、应用结构力学的位移法，求解离散化结构。 （1）取三角形单元的结点位移为基本未知量，它们是 （2）应用插值公式，求出单元的位移函数，即： 其中N 称为形函数矩阵。 （3）应用几何方程，求出单元的应变，即： （4）应用物理方程，求出单元的应力，即： 其中S 称为应力转换矩阵。 i j m ui vi Fxi Fyi um vm Fym Fxm uj vj Fyj Fxj （5）应用虚功方程，求出单元的结点力，即： 其中k 称为单元劲度矩阵。 （6）应用虚功方程，求出单元的结点荷载，即： （7）列出各结点的平衡方程，即： （i=1,2,3,…） 经整理，得整体平衡方程组 其中K 称为整体劲 度矩 阵。 (2)至(6) 步是针对每个单元进行的, 称为单元分析，第( 7) 步是对整个结构进行 的，称为整体分析。 i j m ui vi Fxi Fyi um vm Fym Fxm uj vj Fyj Fxj §6-3 单元的位移模式与解答的收敛性 在三结点三角形单元中，可以假设位移分量只是坐标的线性函数，即： 由左边三个方程求解 ，由右边三个方程求 解 。再代回u、v 中，得 在i、j、m三个结点，位移函数应当等于该结点位移值，即： 其中： 而 （ 轮换） 是单元的结点位移列阵。 而N 称为形函数矩阵。 其中： 位移模式的表达式可用矩