清华大学计算固体力学第三次课件连续介质力学.ppt

5 守恒方程 4 能量守恒 在域内由体积力,和在表面上由面力做的功率为 在物体中总能量的变化率为 由热源s和热流q提供的功率为 其中热流一项的符号是负的，因为正的热流是向物体外面流出的 能量守恒 5 守恒方程 4 能量守恒 即物体内总能量的变化率（包括内能和动能）等于外力的功率和由热流及热能源提供的功率。这是已知的热力学第一定律。 内能的支配依赖于材料。在弹性材料中，它以内部弹性能的形式存储起来，并在卸载后完全恢复；在弹塑性材料中，部分内能转化为热，部分由于材料内部结构的变化而耗散了。 应用Reynold定理将求导数移入积分内，然后将面积分转换为域积分 3 应变度量 Green应变张量E 以位移的形式使用指标写法 代入上式，表示为位移梯度的形式 3 应变度量 在任何刚体运动中，Green应变张量为零，满足了应变度量的一个重要要求。 考虑刚体运动 由变形梯度F 定义，绕原点纯转动时，给出为F＝R （证明见例3.2) 式中转动张量满足正交性，R是正交矩阵 Green应变张量E 第二个运动度量D，称为速度应变, 是变形的率度量, 定义速度梯度 3 应变度量 变形率张量D 速度梯度张量可以分解为对称部分和反对称部分为 令 变形率(对称) 转动(反对称) 二阶张量或方阵的标准分解：以上面的方式，任何一个二阶张量都可以表示为它的对称部分和反对称部分的和 所以 没有变形，转动张量和角速度张量相等：W＝Ω。 由速度梯度定义，在刚体运动中变形率D＝0，所以L＝W ，积分 其中xT和vT是积分常数，对比刚体动力学公式： 得到 在刚体转动中，转动和角速度张量是相同的。当刚体除了转动之外还有变形时，转动张量一般区别于角速度张量。 3 应变度量 变形率张量D 变形率是微小材料线段长度的平方的变化率度量 证明在刚体运动中变形率D＝0 3 应变度量 变形率张量D 3 应变度量 变形率的Green应变率形式 将变形率与Green应变张量的率联系起来，首先得到速度场的材料梯度，并通过链规则表示为空间梯度的形式 取变形梯度 的材料时间导数 应用链规则展开恒等式 得到 代入上面公式，有 3 应变度量 变形率的Green应变率形式 将变形率与Green应变张量的率联系起来 将变形率D前面点积FT，后面点积F，得到 这两种度量是看待相同过程的两种方式：Green应变率是在参考构形中表达的，变形率是在当前构形中表达的。 两种形式的性质的区别是，在例3.7中将会看到Green应变率对时间积分是与路径无关的，而变形率对时间积分是与路径有关的。 逆变换得到 前推运算 后拉运算 例3.5 拉伸和转动联合作用下的应变度量 考虑运动 其中a和b是正常数。计算作为时间函数的变形梯度F，Green应变 和变形率张量，并验证在t＝0与t＝1时的值。定义 计算变形梯度F 以上变形包括同时沿着X和Y轴材料线的拉伸和单元转动。在任何时刻在单元中的变形梯度是常数，应变度量也是常数。得到Green应变张量，由公式给出F，这样得到： 得到Green应变张量 当t＝0时，有x＝X和E＝0， 计算变形率，先获得速度，取运动的材料时间导数 在t＝0时，x＝X，y＝Y，c＝1，s＝0，A＝B＝1，速度梯度在t＝0时为 例3.5 拉伸和转动联合作用下的应变度量 为了确定变形率的时间历史，计算变形梯度的时间导数和逆 等式右边的第一项是变形率，因为它是速度梯度的对称部分， 而第二项是转动，它是反对称部分。变形率在t＝1时给出为 因此，当在中间步骤中，剪切速度－应变是非零的， 在t＝1时刻的构形中只有伸长的速度－应变是非零的。 当t＝1时刻的Green应变率通过对变形率后拉运算给出 例3.5 拉伸和转动联合作用下的应变度量 一个单元经历了图示的变形阶段。在这些阶段之间的运动是时间的线性函数。计算每一阶段的变形率张量D，对于回到未变形构形的整个变形循环，获得变形率的时间积分。 例3.7 计算变形率的时间积分 假定变形的每个阶段都发生在一个单位时间间隔内。时间标定与结果无关，从构形1到构形2的运动为： 确定变形梯度 得到速度梯度和变形率为 例3.7 计算变形率的时间积分 这样，变形率就是一个纯剪切，即两个拉伸分量都为零。由公式(3.3.5)得到Green应变为： 比较上面两式，E22非零，而D22＝0，当a为小量时，E22也小。 从构形2到构形3－剪切与y向拉伸的联合运动： 例3.7 计算变形率的时间积分 从构形3到构形4－纯剪切运动： 从构形4到构形5－y向拉伸（压缩）运动： 在构形