《大变形问题的有限元分析.ppt》由会员分享，可在线阅读，更多相关《大变形问题的有限元分析.ppt（27页珍藏版）》请在三一文库上搜索。
1、2019/2/12,1,第三章 大变形问题的有限元分析,目的：以大变形问题为例，介绍几何非线性问题的有限元方法。 特点：与线性有限元方法比较，几何关系不再是线性的。,内容：,引言 大变形问题的应变描述 大变形分析中的应力描述及本构关系 大变形问题有限元方程的建立 大变形分析中的载荷处理 小结,2019/2/12,2,引言,几何非线性问题：位移与应变成非线性（微分意义上）关系。,物理现象：将位移（转动）和/或应变较大的问题统称为大变形问题，有时称为有限变形问题。这类问题又分为大位移（转动）小应变问题及大位移大应变问题两大类。,研究意义：和材料非线性问题一样重要。例如，平板的弯曲问题，大挠度理论分
2、析结果更符合实际情况；薄壳的屈曲，非线性理论的预测值更好。又例如，对于橡皮型材料，大变形还必须考虑本构关系的变化，这与纯粹的材料非线性又有区别。,几何线性问题： 位移与应变成线性（微分）关系；,研究现状：大变形问题有限元分析的理论和方法存在不同学派间的争鸣，尚未得到一个权威性的结论。随之并发的其它问题，如解的稳定性、收敛性及收敛率等，都有待进一步深入研究。,2019/2/12,3,大变形问题的应变描述(1/4),问题的特点：由于变形较大，使得不同时刻物体具有差别不能忽略的不同构型，这是大变形问题分析的基本出发点。,初始构型（0时刻）,现时构型（t 时刻）,当前构型（ 时刻）,连续介质力学理论对
3、物体经历大变形后的变形有严格的定义和推导。这里不准备过多引入复杂的概念和符号，而是与小变形理论对照，介绍进行大变形分析时必需的几个概念和术语。,大变形问题的分析方法：增量法。,2019/2/12,4,大变形问题的应变描述(2/4),描述的出发点：物体的变形描述建立在确定的参考构型上。,大变形分析由于采用增量方法，需经常用到它们的增量形式。,Green应变张量：以初始构型为参考构型所定义的应变，数学表示为,现时（Updated）Green应变张量：以现时构型为参考构型所定义的应变，数学表示为,注意：我们用下标的大小写表示坐标的大小写，对应于不同的构型。,2019/2/12,5,大变形问题的应变描
4、述(3/4),应变增量：,Green应变增量：,现时（Updated）Green应变增量：,线性部分,非线性部分,线性部分,非线性部分,二者之间满足张量变换关系！,2019/2/12,6,大变形问题的应变描述(4/4),应变增量：（续）对于大变形小应变情形,Green应变增量退化成：,现时（Updated）Green应变增量退化成：,线性部分,非线性部分是高阶小量,线性部分,非线性部分是高阶小量,对于小变形情形,2019/2/12,7,大变形问题的应力描述(1/2),应力是借助于微元体来定义的，但在大变形分析中，必须注意微元体所在的构型。,Euler应力：,与应变类似，连续介质力学理论具有严格
5、的应力定义和多种不同的应力概念。这里也只介绍后面将要用到的几种。,从当前构型中取出微元体，在其上定义的应力称为Euler应力，用 表示。Euler应力代表物体的真实应力。然而，当前构型是待求的未知构型，因而，有必要通过已知构型上的微元体再对应力进行描述。,Kirchhoff应力：,通过初时构型上的微元体定义的应力称为Kirchhoff应力，用 表示；通过现时构型的微元体定义的应力称为现时（Updated）Kirchhoff 应力，用 表示。,2019/2/12,8,大变形问题的应力描述(2/2),Kirchhoff、现时Kirchhoff及Euler应力（增量）间的关系：,根据张量的坐标变换规
6、则，它们之间还有以下关系,现时Kirchhoff应力,Euler应力,现时Kirchhoff应力增量,时刻,t 时刻,特点：以现时构型为参考。,2019/2/12,9,大变形分析中的本构关系(1/5),本构关系的客观性要求：需要选取合适的应力应变共轭对描述材料的本构关系。,弹性材料：加载曲线与卸载曲线相同的材料。,本构关系有三种形式,，,为常数,线弹性材料 (elasticity),超弹性材料 (hyperelasticity),次弹性材料 (hypoelasticity),（大变形分析中）,2019/2/12,10,大变形分析中的本构关系(2/5),弹性材料,若Kirchhoff应力与Gre
7、en应变之间存在一一对应关系，则称这类材料为弹性材料,不依赖于构型变化,弹性本构关系多用于大位移（转动）小应变的情形。,特殊情形,2019/2/12,11,大变形分析中的本构关系(3/5),超弹性材料,假定材料具有单位质量的应变能函数，再根据能量原理来定义本构关系，这类材料称为超弹性材料。,(不限于这种形式）,总之，对于一般的大变形问题，在连续介质力学中常用超弹性来表征材料的本构关系。,例如,一阶近似,初始构型时材料的密度常数,增量形式 ,坐标变换,现时Kirchhoff应力或增量形式 ,Case-1,Case-2,不能简化！,一阶近似,现时构型时材料的密度随变形变化。,相比较,2019/2/
8、12,12,大变形分析中的本构关系(4/5),次弹性材料,若应力率与变形率之间成线性变化规律，这类材料称为次弹性材料。但本构关系描述时要求“率”为与刚体转动无关的客观时间导数。,同乘以时间增量,增量形式 ,Case-2,Case-1,可以证明，这两个率都与转动无关,Jaumann应力率,现时Green应变的线性部分,可以证明，这两个率都与转动无关,旋转率,2019/2/12,13,大变形分析中的本构关系(5/5),三种本构关系间的关系,对于实际的大变形问题，上述三种本构关系并不等价。可以证明，弹性材料是一种特殊的次弹性材料，超弹性材料是一种特殊的弹性材料。,实际材料所遵守的本构关系，只有通过实
9、验测试才能得以确定。,次弹性材料,弹性材料,超弹性材料,2019/2/12,14,大变形问题有限元方程的建立 (1/6),与塑性力学有限元方法的异同,区别：塑性力学的本构关系随加载变化，而大变形问题的构型随加载变化。,TL？UL? 本节讨论,相似：都采用增量方法，都不显含时间。,导致分析方法、应力应变描述、本构关系、控制方程的变化。,构型对应,构型相关，本节讨论 。,客观性描述,2019/2/12,15,大变形问题有限元方程的建立 (2/6),TL法有限元方程的建立,特点：始终以初始（0时刻）构型做为应力与应变描述的参考构型，因而，采用Kirchhoff应力（增量）和Green应变（增量）。,
10、t 时刻:,TL法：Total Lagrangian Description （TLD）,虚功方程：,优点：参考构型不发生变化，本构关系与虚功方程描述形式简单。,时刻:,两式相减，得增量型虚功方程：,2019/2/12,16,大变形问题有限元方程的建立 (3/6),TL法有限元方程的建立（续）,将有限元位移插值、初始构型下的几何关系和本构关系引入后，得到,刚度矩阵,形式较复杂，因问题的类型而不同。,载荷向量,TL法的求解步骤:,Step 1：利用有限元方程求出 间隔内的位移增量 ；,Step 2：利用几何关系，计算Green应变增量 ；,Step 3：利用本构关系，计算Kirchhoff应力增
11、量 ；,Step 4：更新当前时刻 ；更新当前应力 ；计算当前刚度矩阵和载荷向量。,Step 5：转到Step 1，进入下一个时间间隔计算。,2019/2/12,17,大变形问题有限元方程的建立 (4/6),UL法有限元方程的建立,特点：总以t 时刻（即现时构型）为参考构型，也就是说参考构型是变化的，因而，采用现时Kirchhoff应力（增量）和现时Green应变（增量）。,UL法：Updated Lagrangian Description （ULD）,仿照TL法的推导，可得虚功方程：,优点：可以处理加载方式更为复杂的问题，亦可处理边界非线性问题等。,TL法的增量型虚功方程：,2019/2/
12、12,18,大变形问题有限元方程的建立 (5/6),UL法有限元方程的建立（续）,将有限元位移插值、初始构型下的几何关系和本构关系引入后，得到,UL法的求解步骤及与TL法的比较:,Step 1：利用有限元方程求出 间隔内的位移增量 ；,Step 2：利用几何关系，计算现时Green应变增量 ；,Step 3：利用本构关系，计算现时Kirchhoff应力增量 ；,Step 4：更新当前时刻 ；更新当前应力, 根据 计算 ，并且使得 ；更新当前构型 ；计算当前刚度矩阵与载荷向量。,Step 5：转到Step 1，进入下一个时间间隔计算。,2019/2/12,19,大变形问题有限元方程的建立 (6/
13、6),小结,大变形问题有限元方法与弹塑性问题有限元方法都是在增量意义上通过拟线性化，进而加以求解。 但弹塑性问题有限元方法在确定弹塑性状态时还应当进行迭代或按优化问题处理，这点与接触问题类似。 所以，从方法上说，弹塑性问题有限元方法包含了大变形问题有限元和接触问题有限元两类问题的所有特点。,2019/2/12,20,大变形分析中的载荷处理 (1/4),载荷目前还没有考虑,重要区别,TL法的载荷项：,UL法的载荷项：,体积力,表面力,2019/2/12,21,大变形分析中的载荷处理 (2/4),体积力的处理,原则：物体的重力在变形过程中保持不变。,全量,增量,在TL法中，原来的计算方法是正确的；
14、在UL法中，需要按本节的方法计算，可以看出，其差别不能被忽略。,区别之处,2019/2/12,22,大变形分析中的载荷处理 (3/4),表面力的处理,表面力的处理较为复杂，不但与构型变化有关，还与表面力的施加方式有关。以常见的集中力和均布力为例。,1）集中力,第一种情形：集中力的方向在整个变形过程中保持不变；,第二种情形：集中力的方向与所作用表面的夹角不变。,全量,增量,可按表面分布力的特殊情形加以处理，详见下面的分析。,2019/2/12,23,大变形分析中的载荷处理 (4/4),表面力的处理（续）,大变形分析中一般分布载荷随变形的变化，是一个复杂的问题，很难进行定量研究。这里对均匀表面分布
15、力随大变形的变化进行分析。,2）表面分布力,全量,增量,原则：不同时刻，这类载荷的合力一般保持不变。（其它情形仿此进行）,在 t 时刻：,在 时刻：,近似处理后,2019/2/12,24,本章总结 (1/2),由于发生了不可忽略的较大变形，大变形问题的分析更加困难，所涉及内容更加丰富。一般来说，相对于小变形问题，大变形分析具有以下特点：,1）应变定义发生了变化。大变形问题的最鲜明特征就是描述应变位移的几何关系发生了变化。在一些特殊问题中，大变形并没有产生较大的应变，这类问题可略去几何关系中的高阶部分，但必须采用大变形分析方法。实际上，这种对几何关系简化处理，是经过对更一般的几何非线性关系检验后
16、而确立的。,2）构型发生了变化。由于具有较大的变形（位移和/或应变），大变形问题在不同时刻其构型差异较大，必须区别对待 。,3）应力、应变描述发生了变化。由于不同时刻具有不同的构型，并且可以选取不同的参考构型，大变形问题中的应力描述、应变描述以及变分方程都发生了相应的变化。,2019/2/12,25,小结 (2/2),4）求解方法发生了变化。以增量法为基础，大变形问题的求解方法分为Lagrange方法、Euler方法以及Lagrange-Euler混合方法。,5）其它相应的变化。由于具有较大的变形，即使对于弹性材料，本构关系也需要做相应变化，以正确描述大变形情形下材料的本构规律；载荷的作用方式
17、也必须考虑构型的差异，以获得高精度的大变形问题的分析结果。,实际上，当物体发生大变形时，一般都会伴随材料的非线性弹性或非弹性行为，因而必须同时考虑材料和几何两种非线性行为，综合运用两种非线性问题的分析方法。另外，在很多情况下，大变形分析可用来估算结构在失稳前所能承受的最大载荷，如结构屈曲问题等，这些仍是许多学者目前研究的课题。,2019/2/12,26,思考题,举出大位移（转动）小应变和大应变问题这两类大变形问题的典型实例。 根据定义，推导Green应变增量和现时Green应变增量的具体表达式。 解释下列应力变换的物理意义,仿照TL法增量型变分方程，推导UL法的增量型变分方程。 为什么说大变形一般都引起材料性能的非线性行为？能否利用公式加以说明？,2019/2/12,27,谢谢！ 请提问、进行讨论！,