有限元基本理论知识.ppt

2.7 材料非线性问题的有限元法 2.7.1 材料非线性本构关系 一、非线性弹性介质本构关系 2.5.5 坐标变换 母单元几何形状简单而且规则，便于进行运算，但难以适应实际工程中出现的各种结构的复杂形状。通过进行坐标变换，使(ξ，η，ζ)坐标系中形状简单的母单元，在(x，y，z)坐标系中变换为具有曲线(面)边界的形状复杂的单元，变换后的单元称为子单元。子单元在几何上可以适应各种实际结构的复杂外形。经过这样处理，单元具有双重特性：一方面，子单元的几何特征、荷载等等，都来自实际结构，充分反映了实际情况，另一方面，大量计算工作是在母单元内进行的，由于它的形状简单而且规则，计算比较方便，并便于循环，特别有利于在电子计算机上进行计算。因此兼有两方面的优点。 平面坐标变换 二维线性单元 坐标变换公式为 直线24的方程 形心坐标 子单元的4条边都是二次曲线，局部坐标系（ξ，η）是曲线坐标 空间坐标变换 经过空间坐标变换后，原来的直线将变成空间曲线，原来的平面将变成空间曲面。 母单元正六面体，将变为具有曲棱、曲面的六面体子单元。 例 相邻单元公共边连续性验证 2.5.6 位移函数 单元位移用位移函数表示如下： 如果单元坐标变换和位移函数所采用形函数的阶次相等，那么用以规定单元形状的结点数应等于用以规定单元位移的结点数，这种单元称为等参单元。如果坐标变换所用形函数的阶次高于位移函数中的形函数，坐标变换的结点数应超过用以规定单元位移的结点数，这种单元称为超参数单元。反之，如果坐标变换所用形函数的阶次低于位移函数中形函数的阶次，则称为逊参数单元。 写成矩阵形式 等参单元的位移函数满足刚体位移好常应变条件，满足连续性条件，也满足收敛条件 2.5.7 单元应变 空间问题的应变可表示为： 代入位移函数 形函数 是用局部坐标表示的，根据偏微分法则，有 雅可比矩阵计算示例 单元a 单元b 单元c 2.5.8 刚度矩阵 单元刚度矩阵的一般公式 令 将整体坐标系中的积分转换成局部坐标系中积分（注意是矢量乘积） 2.5.9 结点荷载 体积力 产生的结点荷载为 作用于单元边界 上的分布力 所产生的结点荷载按下式计算 式中：[N]为表面 Ω 的形函数矩阵，是3×3s阶矩阵，s为单元表面Ω 的结点数 设 所作用的表面Ω是ζ＝1的面，在坐标变换公式中令ζ＝1可得到Ω的方程 表面Ω上任一点的压强为 设表面Ω任一点的法线的方向余弦为l，m，n，则表面力在各方向的分力为 2.5.10 等参数单元的退化 4结点平面等参数单元的退化 令 8结点空间等参单元的退化 高次单元的退化 高次单元退化后必须修改形函数 2.5.11 数值积分 在求解刚度矩阵和结点荷载时，需计算如 的积分。但 一般是很复杂的，通常难以用显式表示其积分，一般都用数值积分方法计算积分 值，即在单元内选出某些点，称为积分点，求出被积函数 在这些点的 值，然后根据这些数值求出积分值。 数值积分有两类方法，一类方法积分点是等间距的，如辛普生方法；另一类方法 积分点是不等间距的，如高斯方法。 一维高斯积分公式 0.3478548451 0.6521451549 n＝4 0.8611363116 0.3399810436 0.5555555556 0.8888888889 n＝3 0.7745966692 0.0000000000 1.0000000000 n＝2 0.5773502692 Hi ±ξi 和 是根据计算精度最高而选定的，积分点 应是勒让德多项式 的根。 加权系数 按下式计算。 二维及三维高斯积分公式 先令 保持常数，计算沿 方向的积分 再沿 方向积分 对三重积分有 一般采用2×2×2高斯积分 2.6 非线性有限元分析方法 线性弹性力学采用两个基本假定：1。材料的应力应变关系是线性的，即假定材料符合胡克定律 2。应变位移关系是线性的，即小位移假定。 例如当钢材的应力超过其比例极限后，应力应变关系便是非线性的。又如土壤和岩石 的应力应变关系也是非线性的。这些称为材料非线性。 又如梁、板及薄壳等结构失稳后，由于产生了大位移，其应变位移关系是非线性的， 这些称为几何非线性。 在热传导问题中，某些情况下材料的导温系数及内部热源与温度有关。在流体力学中， 粘滞系数与流速有关，或者由于出现紊动，达西定律不再适应。这些问题都是非线性的。 当材料的应力应变关系是非线性的时，刚度矩阵不是常数，而与应变和变位值有关。 可记为 。这时结构的