[高二数学]第4章 平面问题有限单元法3.pptVIP

如以无量纲坐标系来表达，则上式可以写成 其中： 单元应变场 　　根据单元的位移场函数式，由几何方程可以得到单元的应变场表达式， 记为： 　　这里，B矩阵称为几何矩阵。B矩阵可以表示为分块矩阵的形式 其中 单元应力场 　　由物理方程及应变矩阵，可以得到单元的应力场表达式， 其中　　　为应力矩阵，D称为弹性矩阵，对于平面应力问题， 将应力矩阵表示为分块矩阵的形式 其中： 对于平面应变问题，只需将E换为　　　，　换为　　　。 单元刚度矩阵 　　和三角形单元一样，可以根据虚功理导出节点位移向量和节点力向量之间关系，即单元的刚度矩阵　　，可以将其写成分块的形式。 其中 　　对于平面应力问题，如果单元厚度t为常数，则可得刚度矩阵的显式形式： 课后作业 如图（a）所示一高深悬臂梁，在右端部受集中力F作用，材料弹性模量E、泊松比v＝1/3，悬臂梁的厚度（板厚）为t，如图（b）所示有限元模型，试按平面应力问题，求结构总体刚度矩阵和总体载荷向量，引入边界条件，求解节点1，2的位移。 再 见 * * 李建宇 天津科技大学 工程中的有限元方法 Finite Element Method in Engineering 内容 4 平面问题的有限元法 3 1 计算实例 2 位移边界条件的处理 3 平面矩形单元分析 要求 理解：有限元分析的基本步骤 引入位移边界条件的原理 平面矩形单元分析的原理 掌握： 引入位移边界条件的置“1”法和乘大数法 平面问题三角单元有限元分析整体流程 课后作业 完成边界条件的引入和刚度方程的求解 回顾 连续体有限元分析的基本流程 整体 离散 单元分析 单元 组装 引入边界条件整体解算 连续体结构 人工节点 逼近离散 单元刚度方程 总体刚度方程 总体刚度矩阵的组装 组装原理： 节点位移协调条件 节点平衡条件 i j m ③ ④ ① ② ③ ③ ① ① ④ ④ i i i i ) 1 ( i V yi P xi P ) e ( i U ? ) 3 ( i V ) 3 ( i U ) 1 ( i U ) 4 ( i V ) 4 ( i U ) e ( i V ? (a) (b) (c) 位移协调条件：各单元共享节点的位移相等 节点平衡条件：各单元内力与节点外力构成平衡力系 回顾 回顾 回顾平面连续体的虚功原理： 总体刚度矩阵组装的另一种解释 内力总虚功＝外力总虚功 将连续体区域离散化为单元的集合，即 对各单元，引入位移形函数[N]，则有 位移： 应变： 应力： 代入虚功方程，得 内力虚功＝ 外力虚功＝ 回顾 总体刚度矩阵组装的另一种解释 内力虚功＝ 外力虚功＝ 对各单元，引入单元节点位移提取矩阵[T]e，则有 引入总体节点位移向量： 单刚： 总刚： 总体外载荷向量： 回顾 关于总体外载荷向量的注解 总体外载荷向量： 单元外载荷向量： 1×6 (1×2) (2×6) (1×2) (2×6) 上式也是计算单元等效节点载荷的普遍公式。 1、集中力的等效 关于总体外载荷向量的注解 2、体力的等效 关于总体外载荷向量的注解 如重力、惯性力、磁力。将体积力的微元即可以看作集中力。设体力集度 微元体上的合力 积分 3、面力等效 如均布、线性分布，作用平面内的力 单位面力（面力度） 微元体上 面积分求和 注意： 1 ．同时作用面力、体力、集中力移置可单独移置，最后相加（小变形、弹性体）； 2 ．特例：线性分布的力可按静力学平行力系的分解原理直接求出。 关于总体外载荷向量的注解 ①均质等厚的等边三角形重力（体力）w ：1/3 w 移置到各节点上； ②如图，i ，j 边L 上有分布力q ，假设单元厚度为t，则合力为qtL，可将qtL/2 移置到各节点上； ③ i ，j 边L上有线性分布力，i 点大小为0，j 点大小为q ，假设单元厚度为t，则合力为qtL/2，将1/3*qtL/2移到i 点，2/3* qtL/2 移到j 点。 q i j i j q qtL/2 qtL/2 1/3*qtL/2 2/3* qtL/2 回顾 常用节点等效外载荷向量 计算实例 如图为一边长为a、厚度为t的正方形薄板。其中AB边固定，BC、CD边自由，AD边作用均