[工学]ANSYS平面问题单元分析.ppt

Plane problems：Element analysis Step2几何关系——应变与位移 应变：单位长度的改变 A dy dx P O A’ B’ B P’ 即 将位移模式 代入几何关系 得： 或写为： 几何矩阵(应变矩阵） 常数 因此，三节点三角形单元也叫常应变单元。 Step3 物理关系 For linear elastic isotropic materials (Hookes law) ： Plane stress 得： ： 平面应力 弹性矩阵 Plane strain 得： 平面应变弹性矩阵： 应力矩阵 ： Hooke‘s law 胡克定律 平面应力问题的应力矩阵 平面应变问题的应力矩阵将上式中的弹性常数E、μ换成 即可。 Step4 虚功原理 the Principle of Virtual Work ——单元节点力向量 ——单元节点虚位移向量 ——单元内一点的应力向量 ——单元内一点的虚应变向量 虚功原理： 厚度= t 外力虚功 内力虚功 因为有 所以有 代入 因虚位移是任意的，等号两边可消去 至此,完成平面问题单元分析！ 即 三节点三角形单元的单元刚度矩阵： ???? Summary---单元刚度矩阵特点 以三节点三角形单元的单刚为例 物理意义：单位节点位移分量所引起的节点力分量。如 k16 是m节点有单位垂直位移而其它节点无位移时，在i 节点上所引起的水平节点力。 * 什么是平面问题 ？What is a plane problem? 什么是单元？What is an element ? 什么是单元分析？What is an element analysis? 总结 Summary 弹性力学研究对象都为空间问题，但在一定条件下，有些可简化为平面问题来研究，从而大大减少计算工作量。 ???? 什么是平面问题? 研究对象 例：齿轮、发动机连杆等 沿厚度均布 且平行板面 t a b F1 F2 几何形状：薄板 tb 受力状况：平行于板面且沿厚度均布，板面上不受力。 应力特点：只平行板面的三个应力分量不为零 平面应力问题 FR FR 齿轮 平行xoy平面 例：花键、拦水大坝等 平面应力问题中 σz=0 但ε z≠0 平面应变问题中ε z=0 但σ z≠0 平面应力问题和平面应变问题统称为平面问题。 平面应变问题 注意 几何形状：长柱体l a 受力状况：平行于柱体横截面且沿长度均布. 应变特点：只平行板面的三个应变分量不为零 练习 判断下述问题是否平面应力或平面应变问题? P P P P q q q q q F q q F p q W 6 Stress Components 正应力拉为正压为负 正面上的剪应力沿着坐标轴正向时为正，反之为负。负面上的剪应力相反 什么是单元？ ???? 一个单元是组成弹性体的部分连续质点的集合。 discretization Geometry Finite Element Mesh Elements Nodes 说明：1、单元之间仅在节点处连接； 2、外力被等效到有关节点上； 3、各单元内的几何、物理特性是均匀的（uniform） 无 限 个 自 由 度 有 限 个 自 由 度 不是越细越好。 应力集中处单元应划细些。 单元几何形状应避免畸形。 不同厚度、不同材料划在不同单元。 有集中载荷处或分布载荷突变点一般应安排节点。 划分单元应注意： 平面单元的形状 不是越细越好 O x y A O x A 1 2 3 4 5 6 7 8 9 y O x y A 1 5 25 21 O x A y 1 9 81 73 节点数 100 % 90 95 9 25 81 精确解 应力集中处划细些 几何形状有变化（有应力集中）的地方单元要密由密到疏 平缓过渡 避免畸形单元 O x y A 1 12 x O y A 1 12 A 1 12 O x y O x y 1 12 A I II III Ⅳ Ⅴ x O y A 1 12 长宽比 100 % 90 95 1 4 8 精确解 2 3 5 6 7 85 Ⅴ Ⅳ II III I 不同厚度、不同材料划在不同单元 有集中载荷处或分布载荷突变点一般应安排节点 能剖分的单元形状由程序决定，非边缘节点一般不能只属于一个单元 单元剖分完后，将全部单元和节点编号，从1开始不得遗漏或重复。 Fig.1 Fig.2 Fig.3 Fig.5 Fig.4 i j m p i j m p i j m p i j m p i j m p 剖分练习 q 孔 t1 t2 t2 E2 μ2 E1 μ1 p F ???? 什么是单元分析？ 单元分析：求单元节点力与节点位