有限元绝密讲义2.ppt

第2章 平面问题的有限元方法 2.1 弹性理论基础 Ⅰ、基本假设： 连续性－物质连续。相应的应力应变，位移等连续变量可以用坐标的连续函数表示； 均质各向同性——物体内部各点，各方向上物理性质相同，材料常数（弹性模量，泊松比）不随坐标方向而变； 完全弹性——材料服从Hooke定律； 小变形（几何假设）——略去二阶小量，所有微分方程为线性的； 无初应力——加载前物体内无初应力。 Ⅱ、基本方程 A、平衡方程： (2.1) 对于平面问题： (2.2) B、几何方程： (2.3) 对于平面问题： (2.4) (2.5) C：物理方程（应力－应变关系）： (2.6) 为弹性矩阵——由弹性模量和泊松比构成，为常数矩阵。 平面问题一般分为两类： ⑴：平面应力状态 仅在平面内存在应力分量，即 注意： (2.7) ⑵：平面应变状态 仅在平面内存在应变分量，即 注意： (2.8) 实际上，只需对（2.7）式中的 、 作如下变换，即能得到（2.8）式。 D、边界条件： 力的边界条件（ ）： (2.9) 为边界外法线的方向余弦 位移边界条件（ ）： Ⅲ、能量原理 应变能： 虚功原理：变形体中满足平衡的力系在任意满足协调条件的变形状态上作的虚功等于零。 体系外力的虚功与内力的虚功相等。 虚位移原理：如果弹性体处于平衡状态，当发生约束所允许的任意微小虚位移时，外力在虚位移上所作的功等于弹性体内的应力在相应虚应变上所作的功。 设弹性体在体积力 ，表面力 作用下处于平衡，体内引起相应的位移 ，应力 和应变 。 假定弹性体内有某约束条件所允许的虚位移 ，相应有虚应变 ，则虚位移原理可表述为： (2.10) 平面问题：设板厚t为常量，即 则， （2.11） 以后证明：上式为平衡方程和力边界条件的等效积分提法 虚应力原理：几何方程和位移边界条件的等效积分提法 应力、约束反力为虚拟量，而位移为确定量。 注意：虚位移或虚应力原理中，均未涉及物理方程。 也就是说，线弹性，非线弹性，弹塑性问题均适用； 但是，平衡方程、几何方程都是在小变形条件下建立的，因此，二原理不能直接应用于大变形问题。 更多的原理将在后续章节中阐述。 2.2 分析流程Ⅰ 对于任一平面问题，分析步骤的第一步是对结构离散化。 考虑图示一平板（平面应力或平面应变不限） 弹性理论求解时，取无限小微元体作为出发点，建立平板的力平衡、几何变形的微分方程形式，为一无限自由度问题。 而具体求解时可采用半逆解法，引入一多项式应力函数，定其系数，解出应力，积分得到位移，最后结果为一全域连续的分布函数。 当采用有限元方法求解时，第一步是将平板离散成有限个小单元。 A： 梁结构的离散：取一段梁为一单元 单元类型：简单直线段 离散原则：几何上真实模拟原结构及其变形 平板的离散：取一小面积板为一单元 单元类型：由最基本的平面图形构成 三角形、四边形（如正方形、长方形、梯形） 而五边形、圆、扇形不宜作为单元。 离散原则：几何上真实模拟原结构（无缺陷、重叠） 模拟变形状态 B：节点的选择 梁单元：取端点为节点，单元内任一点的几何位置可由节点坐标确定，单元内任一点的变形同样能由节点的变形确定。 对二节点梁单元，单元内任一点的坐标、变形可由节点的坐标、变形的线性函数插值确定。 平面单元：取角点为节点，单元内任一点的几何位置、变形可由节点的坐标几变形插值确定。 三角形单元节点：三个角点？ 矩形单元的节点：四个角点？ C：节点位移 梁单元中节点的位移分量对应直梁、平面梁、空间梁分别有1位移+1转动、2位移+1转动、3位移+3转动分量； 平面单元的节点只有2位移： 为什么？ 相应地，节点载荷： D：离散结果 图示平板可分离成18个三角形单元、16个节点（32个位移分量）。 无限自由度转换成32个自由度！ 计算结果：仅16个节点的位移数值； 可以将节点的位移拟合成一条曲线；当单元趋 小，节点增多，拟合曲线才逐渐逼近真解。 与弹性理论解的形式差别：不直观； 结果比较：不能直接反