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

TL法有限元方程的建立 特点：始终以初始（0时刻）构型做为应力与应变描述的参考构型，因而，采用Kirchhoff应力（增量）和Green应变（增量）。 t 时刻: TL法：Total Lagrangian Description （TLD） 虚功方程： 优点：参考构型不发生变化，本构关系与虚功方程描述形式简单。 时刻: 两式相减，得增量型虚功方程： 第三十一页，共三十九页，2022年，8月28日 TL法有限元方程的建立（续） 将有限元位移插值、初始构型下的几何关系和本构关系引入后，得到 刚度矩阵 形式较复杂，因问题的类型而不同。 载荷向量 TL法的求解步骤: Step 1：利用有限元方程求出 间隔内的位移增量 ； Step 2：利用几何关系，计算Green应变增量 ； Step 3：利用本构关系，计算Kirchhoff应力增量 ； Step 4：更新当前时刻 ；更新当前应力 ；计算当前刚度矩阵和载荷向量。 Step 5：转到Step 1，进入下一个时间间隔计算。 第三十二页，共三十九页，2022年，8月28日 UL法有限元方程的建立 特点：总以t 时刻（即现时构型）为参考构型，也就是说参考构型是变化的，因而，采用现时Kirchhoff应力（增量）和现时Green应变（增量）。 UL法：Updated Lagrangian Description （ULD） 仿照TL法的推导，可得虚功方程： 优点：可以处理加载方式更为复杂的问题，亦可处理边界非线性问题等。 UL法的增量型虚功方程： 第三十三页，共三十九页，2022年，8月28日 UL法有限元方程的建立（续） 将有限元位移插值、初始构型下的几何关系和本构关系引入后，得到 UL法的求解步骤及与TL法的比较: Step 1：利用有限元方程求出 间隔内的位移增量 ； Step 2：利用几何关系，计算现时Green应变增量 ； Step 3：利用本构关系，计算现时Kirchhoff应力增量 ； Step 4：更新当前时刻 ；更新当前应力, 根据 计算 ，并且使得 ；更新当前构型 ；计算当前刚度矩阵与载荷向量。 Step 5：转到Step 1，进入下一个时间间隔计算。 第三十四页，共三十九页，2022年，8月28日 小结 大变形问题有限元方法与弹塑性问题有限元方法都是在增量意义上通过拟线性化，进而加以求解。 但弹塑性问题有限元方法在确定弹塑性状态时还应当进行迭代或按优化问题处理，这点与接触问题类似。 所以，从方法上说，弹塑性问题有限元方法包含了大变形问题有限元和接触问题有限元两类问题的所有特点。 第三十五页，共三十九页，2022年，8月28日 6.大变形问题的载荷处理 第三十六页，共三十九页，2022年，8月28日 载荷目前还没有考虑 重要区别 TL法的载荷项： UL法的载荷项： 体积力 表面力 第三十七页，共三十九页，2022年，8月28日 土木工程中的几何非线性问题 第一页，共三十九页，2022年，8月28日 基本概念 几何非线性问题：位移与应变成非线性（微分意义上）关系。 物理现象：将位移（转动）和/或应变较大的问题统称为大变形问题，有时称为有限变形问题。这类问题又分为大位移（转动）小应变问题及大位移大应变问题两大类。 研究意义：和材料非线性问题一样重要。例如，平板的弯曲问题，大挠度理论分析结果更符合实际情况；薄壳的屈曲，非线性理论的预测值更好。又例如，对于橡皮型材料，大变形还必须考虑本构关系的变化，这与纯粹的材料非线性又有区别。 几何线性问题： 位移与应变成线性（微分）关系； 研究现状：大变形问题有限元分析的理论和方法存在不同学派间的争鸣，尚未得到一个权威性的结论。随之并发的其它问题，如解的稳定性、收敛性及收敛率等，都有待进一步深入研究。 第二页，共三十九页，2022年，8月28日 1.物体运动的描述 第三页，共三十九页，2022年，8月28日 拉格朗日描述 t=0的坐标为Xi，t时刻位置为xi，质点运动可表为 对物体t时刻位置和变形的刻划称为构形(configuration)，如