土木工程中的几何非线性问题.pptVIP

非线性有限元 基本概念 1.物体运动的描述 拉格朗日描述 欧拉描述 2.大变形问题的应变描述 应变张量 3.大变形问题的应力描述 4.大变形问题的本构关系 5.大变形问题的有限元方程 6.大变形问题的载荷处理 UL法有限元方程的建立 特点：总以t 时刻（即现时构型）为参考构型，也就是说参考构型是变化的，因而，采用现时Kirchhoff应力（增量）和现时Green应变（增量）。 UL法：Updated Lagrangian Description （ULD） 仿照TL法的推导，可得虚功方程： 优点：可以处理加载方式更为复杂的问题，亦可处理边界非线性问题等。 UL法的增量型虚功方程： UL法有限元方程的建立（续） 将有限元位移插值、初始构型下的几何关系和本构关系引入后，得到 UL法的求解步骤及与TL法的比较: Step 1：利用有限元方程求出 间隔内的位移增量 ； Step 2：利用几何关系，计算现时Green应变增量 ； Step 3：利用本构关系，计算现时Kirchhoff应力增量 ； Step 4：更新当前时刻 ；更新当前应力, 根据 计算 ，并且使得 ；更新当前构型 ；计算当前刚度矩阵与载荷向量。 Step 5：转到Step 1，进入下一个时间间隔计算。 小结 大变形问题有限元方法与弹塑性问题有限元方法都是在增量意义上通过拟线性化，进而加以求解。 但弹塑性问题有限元方法在确定弹塑性状态时还应当进行迭代或按优化问题处理，这点与接触问题类似。 所以，从方法上说，弹塑性问题有限元方法包含了大变形问题有限元和接触问题有限元两类问题的所有特点。 载荷目前还没有考虑 重要区别 TL法的载荷项： UL法的载荷项： 体积力 表面力 体积力的处理 原则：物体的重力在变形过程中保持不变。 全量 增量 区别之处 * 几何非线性问题：位移与应变成非线性（微分意义上）关系。 物理现象：将位移（转动）和/或应变较大的问题统称为大变形问题，有时称为有限变形问题。这类问题又分为大位移（转动）小应变问题及大位移大应变问题两大类。 研究意义：和材料非线性问题一样重要。例如，平板的弯曲问题，大挠度理论分析结果更符合实际情况；薄壳的屈曲，非线性理论的预测值更好。又例如，对于橡皮型材料，大变形还必须考虑本构关系的变化，这与纯粹的材料非线性又有区别。 几何线性问题： 位移与应变成线性（微分）关系； 研究现状：大变形问题有限元分析的理论和方法存在不同学派间的争鸣，尚未得到一个权威性的结论。随之并发的其它问题，如解的稳定性、收敛性及收敛率等，都有待进一步深入研究。 t=0的坐标为Xi，t时刻位置为xi，质点运动可表为 对物体t时刻位置和变形的刻划称为构形(configuration)，如图示。该描述实质是给出初始位置坐标为Xi 的质点运动轨迹。 描述运动的参照基准称为参考位形，以初始位形作参考位形的描述称为物质描述或拉格朗日描述，Xi称为物质坐标, Xi和t称为拉格朗日变量。 现时构形中，经过空间位置xi的质点的速度为vi 研究不同时刻经过同一空间点xi的质点的运动状态 以现时位形作参考位形的描述称为空间描述或欧拉描述，xi称为空间坐标, xi和t称为欧拉变量。 注意：两种描述下对某个质点加速度的描述是不一样的 物体现时坐标xi对物质坐标Xi的偏导数 称为变形梯度，是非对称的二阶张量。(思考？什么时候其是对称的) 因此可以将变形梯度视作一种线性变换，它将参考位形中的线元dXi变换为现时位形中的线元dxi，这变换中既有伸缩，也有转动。变形梯度在大变形分析中很重要。 现时位形两邻点的距离为 点的变换 物体运动和变形是单值和连续的，也即在任一时刻， 和 是一一对应的，那么在参考位形的任意点Jacobi行列式J不为零。也即变形梯度可逆 设图示初始位形微元体体积为dV0，三线元为 运动变形后，现时位形三线元为 体积变换 体积变换公式 dV=JdV0 图示面元可表示为 如果记初始和现时位形的密度分别为 则由质量守恒，可得 因此对不可压缩物体 面积变换 面积变换公式 关注P、Q两点的距离 研究变形前后线段尺度的变化 可以获得变形的度量－应变 格林应变张量 阿尔曼西张量 格林应变张量用初始位形定义，也即用变形前的坐标定义，它是lagrange坐标的函数。 阿尔曼西应变张