第7讲非线性专题.pptVIP

北京航空航天大学 也称为运动的Jacobian矩阵 使用Jacobian矩阵行列式可以将当前构形和参考构形上的积分联系起来： 变形梯度定义 关于变形梯度 (反映变形特征) 北京航空航天大学 物体质量在变形过程中是不变的。 关于质量守恒定律 北京航空航天大学 应变的度量 有限变形情况下，小变形情况下的应变度量方法已经不再适用，必须重新定义。 刚体转动，应变度量必须为零。小变形情况下的应变度量方法不能满足这个要求 非线性连续介质力学中使用了多种不同的应变和应变率度量。在有限元方法中使用最普遍的有两种： Green应变张量Eij 变形率张量Dij（率形式的本构方程中） 北京航空航天大学 初始构形中P0Q0距离的平方： 现时构形中PQ距离的平方： Green应变张量的定义 北京航空航天大学 变形前后微段长度平方的变化： 格林(Green)应变张量——初始构形中定义 进一步引入位移： 北京航空航天大学 对于小变形问题 Green应变张量简化为无限小应变张量 ： 北京航空航天大学 Green应变不能直接反映变形的大小，它本身不是伸长度和角变度，而是与伸长度和角变度有一定的关系。Green应变和工程(名义)应变分量的关系为： 北京航空航天大学 变形率张量Dij的定义 为了定义变形率张量的表达式，首先定义速度梯度 变形率张量 转动张量 变形率张量Dij在现时构形中定义，其时间积分和真应变相对应。 北京航空航天大学 应变度量的客观性 Green应变张量Eij和变形率张量Dij都是不随刚体转动而变化的客观张量。 对于纯刚体运动，应变度量必须为零。否则预示着纯刚体运动中，将导致非零应力。这是线性应变位移方程在非线性理论中被放弃的主要原因。 举例具体说明如下： 北京航空航天大学 考虑物体的刚体运动： （绕原点的转动和平动） 北京航空航天大学 格林(Green)应变张量： 写成矩阵形式： 对于物体的刚体运动： 则变形梯度： 正交矩阵 所以： 北京航空航天大学 对比：线性应变位移方程 如果 较大，伸长应变不为零。因此，线应变张量不能用于大变形问题，即几何非线性问题。 北京航空航天大学 大变形情况下，是在变形后的物体中截取微元体建立平衡方程和与之等效的虚功(率)原理的。 如果应变是用变形前的坐标表示的Green应变张量，则需要定义与之对应的关于变形前构形的应力张量。 0 应力的度量 北京航空航天大学 Cauchy应力张量 (表征现时构形中质点的应力状态) Lagrange应力张量 Tij (第一类Piola-Kirchhoff应力张量) Kirchholf应力张量 Sij (第二类Piola-Kirchhoff应力张量) 三者之间关系 北京航空航天大学 Cauchy应力张量在现时构形中定义，具有明确的物理意义，代表真实应力。 Cauchy应力张量和变形率张量Dij在功率上共轭（或在功或能量上耦合）。（给出内部功率） Lagrange应力张量和Kirchholf应力张量在参考构形上定义，没有实际的物理意义，仅仅是为了构成能量上共轭的需要。 Lagrange应力张量是非对称张量，很难实际应用。因此进一步引入对称的Kirchholf应力张量。Kirchholf应力张量被广泛应用于与路径无关材料（如橡胶），并且和Green应变率在功率（能量）上耦合。 北京航空航天大学 弹塑性问题中本构方程的客观性 在大变形弹塑性有限元分析中，需要采用增量理论来分析依赖于材料变形历史的塑性变形过程，因而通常采用微分型或速率型的材料本构关系。采用速率型的本构关系也应该具有客观性。 可以证明，在固定坐标系下，Cauchy应力率张量不具有客观性（不能反映刚体运动）。因此定义一种不随刚体转动而变化的速率型的应力张量——Cauchy应力张量的Jaumann率 另外两种经常使用的率是Truesdell率和Green-Naghdi率。（构造方法参见相关文献） 北京航空航天大学 有限变形条件下的弹塑性本构方程与小变形弹塑性本构方程在形式上是相同的，但其中的应力速率和应变速率要具有客观性。 如果在有限元分析中采用随动共旋坐标系，本构方程可以直接使用Cauchy应力率张量和变形率张量。 习惯上 北京航空航天大学 7.4.3 几何非线性问题的表达格式 两种表达格式(根据参考构形的情况)： 完全的Lagrange格式——变量总是参考初始构形，一般采用全量型的本构关系。 更新的Lagrange格式——变量参考每一载荷步开始时的构形，即分析过程中参考构形是不断被更新的。一般采用速率型的本构关系。 可根据材料本构关系的特点选择最有效的格式，如下表所示： 北京航空航天大学 几何非线性 问题 (大位移、大转动) 材料本构关系 弹性问题 非弹性问题 小应变 大应变