【连续介质力学】张量的偏微分、球张量和偏张量

LAWKAWAI

于 2023-05-20 20:52:58 发布

阅读量1.2k

点赞数

分类专栏：连续介质力学文章标签：笔记

本文链接：https://blog.csdn.net/weixin_38362786/article/details/130703102

版权

连续介质力学专栏收录该内容

16 篇文章 8 订阅

订阅专栏

张量的偏微分

张量的一阶微分，定义：
$\frac{\partial A}{\partial A} = A_{,A}=\frac{\partial A_{ij}}{\partial A_{kl}}(\hat e_i\bigotimes\hat e_j\bigotimes\hat e_k\bigotimes\hat e_l )\\=\delta_{ik}\delta_{jl}(\hat e_i\bigotimes\hat e_j\bigotimes\hat e_k\bigotimes\hat e_l ) = I$

张量的迹的微分：
$\frac{\partial Tr(A)}{\partial A} = Tr(A)_{,A}=\frac{\partial A_{kk}}{\partial A_{ij}}(\hat e_i\bigotimes\hat e_j )\\=\delta_{ki}\delta_{kj}(\hat e_i\bigotimes\hat e_j ) = \delta_{ij}(\hat e_i\bigotimes\hat e_j ) = 1$

张量的迹平方的微分：
$\frac{\partial Tr(A)^2}{\partial A} = 2Tr(A)Tr(A)_{,A}=2Tr(A)1$

张量平方的迹的微分：
$\frac{\partial Tr(A^2)}{\partial A} =\frac{\partial A_{sr}A_{rs}}{\partial A_{ij}}(\hat e_i\bigotimes\hat e_j )\\=[A_{rs}\frac{\partial A_{sr}}{\partial A_{ij}} + A_{sr}\frac{\partial A_{rs}}{\partial A_{ij}}](\hat e_i\bigotimes\hat e_j )\\ =[A_{rs}\delta_{si}\delta_{rj}+A_{sr}\delta_{ri}\delta_{sj}](\hat e_i\bigotimes\hat e_j )\\= [A_{ji}+A_{ji}](\hat e_i\bigotimes\hat e_j )=2A_{ji}(\hat e_i\bigotimes\hat e_j )=2A^T$

张量立方的迹的微分：
$\frac{\partial Tr(A^3)}{\partial A} =\frac{\partial A_{pq}A_{qr}A_{rp}}{\partial A_{ij}}(\hat e_i\bigotimes\hat e_j )\\=[A_{qr}A_{rp}\frac{\partial A_{pq}}{\partial A_{ij}} + A_{pq}A_{rp}\frac{\partial A_{qr}}{\partial A_{ij}} + A_{pq}A_{qr}\frac{\partial A_{rp}}{\partial A_{ij}} ](\hat e_i\bigotimes\hat e_j )\\ =[A_{qr}A_{rp}\delta_{pi}\delta_{qj} + A_{pq}A_{rp} \delta_{qi}\delta_{rj} + A_{pq}A_{qr}\delta_{ri}\delta_{pj}](\hat e_i\bigotimes\hat e_j )\\= [A_{jr}A_{ri}+A_{pi}A_{jp}+A_{jq}A_{qi}](\hat e_i\bigotimes\hat e_j )\\ =3A_{jr}A_{ri}(\hat e_i\bigotimes\hat e_j ) =3(A^2)^T$

对于对称张量：
$\frac{\partial Tr(C)}{\partial C} = 1 \\ \frac{\partial [Tr(C)]^2}{\partial C}=2Tr(C) 1 \\ \frac{\partial [Tr(C^2)]}{\partial C}=2C^T=2C \\ \frac{\partial [Tr(C^3)]}{\partial C}=3(C^2)^T=3C^2$

张量的范数的微分：
$\frac{\partial ||C||}{\partial C}=\frac{\partial \sqrt{C:C}}{\partial C} = \frac{\partial (\sqrt{Tr(C\cdot C^T)})}{\partial C} \\= \frac{\partial \sqrt{Tr(C^2)}}{\partial C}=\frac{1}{2}[Tr(C^2)]^{-\frac{1}{2}}[Tr(C^2)]_{,C}\\=\frac{1}{2}[Tr(C^2)]^{-\frac{1}{2}}2C$

$\frac{\partial ||C||}{\partial C} = \frac{C}{\sqrt{Tr(C^2)}}=\frac{C}{||C||}$

有趣的微分：
在这里插入图片描述

张量的逆的微分：
由于： $\frac{\partial 1}{\partial C}=\frac{\partial (C^{-1}\cdot C)}{\partial C}= 0$
在这里插入图片描述
由于： $C_{qj}= \frac{1}{2}(C_{qj}+C_{jq})$

张量形式：
$\frac{\partial C^{-1}}{\partial C}=\frac{1}{2}[C^{-1}\overline \bigotimes C^{-1}+C^{-1} \underline \bigotimes C^{-1}]$

NOTE: 如果C不是对称的，那么
$\frac{\partial C_{iq}^{-1}}{\partial C_{kl}}\delta_{qr}=-C_{iq}^{-1}\frac{\partial C_{qj}}{\partial C_{kl}}C_{jr}^{-1}=-C_{iq}^{-1}\delta_{qk}\delta_{jl}C_{jr}^{-1} =-C_{ik}^{-1}C_{lr}^{-1}$
不是对称的

不变量的偏微分

$I_T$ 的微分：
$\frac{\partial I_T}{\partial T}=\frac{\partial Tr(T)}{\partial T} =Tr(T)_{,T}=1$

$II_T$ 的微分：
在这里插入图片描述

应用Cayley-Hamilton定理：
在这里插入图片描述
将上式T表达式，代入 $II_T$ 的表达式：

第三不变量 $III_T$ 的微分：
在这里插入图片描述
再次应用Cayley-Hamilton定理：

转置：

通过比较，可以求出另一种表示 $III_T$ 的表达式：
$\frac{\partial III_T}{\partial T}=(III_T T^{-1})^T=III_T T^{-T}$

张量的时间偏导

定义：
$\frac{D}{Dt}T = \dot T \quad \frac{D^2}{Dt^2}=\ddot T$

张量的行列式的时间偏导：
$\frac{D}{Dt}[\det T] = \frac{DT_{ij}}{Dt}cof (T)$

其中， $co f (T)$ 是T的余子式， $cof[T_{ij}]]^T = \det (T )(T^{-1})_{ij}$

问题1.38 考虑 $[\det b]^{\frac{1}{2}} = (III_b)^{\frac{1}{2}}$ ， $b$ 是二阶对称张量，求出 $J$ 和 $\ln J$ 的关于 $b$ 的偏导

在这里插入图片描述

球张量和偏张量

任意一个张量都可以分解成球张量和偏张量：
$T^{sph}+T^{dev}=\frac{Tr(T)}{3}1+T^{dev}=\frac{I_T}{3}1+T^{dev}=T_m1+T^{dev}$

所以，偏张量的定义：
$T^{dev} = T - \frac{Tr(T)}{3}1=T - T_m1$

由于张量T是对称的， $T = T^T$ ，所以：
在这里插入图片描述
在笛卡尔坐标系表示球张量和偏张量：

接下来介绍根据张量T主不变量的偏张量不变量

偏张量的第一不变量

在这里插入图片描述
任意偏张量的迹都为0

偏张量的第二不变量

在主空间中，张量T的分量是：
在这里插入图片描述
主不变量为： $I_T = T_1+T_2+T_3; \quad II_T = T_1T_2+T_2T_3+T_3T_1; \quad III_T=T_1T_2T_3$

那么偏张量 $T^{dev} = T - T_m1$ 在主空间的分量为：
在这里插入图片描述
所以，偏张量的第二不变量是：

同样地，也可以用从第二不变量的定义出发证明：
定义：

那么偏张量的第二不变量为：

因为不变量不随坐标系的改变而改变，所以在主空间和一般的笛卡尔坐标系中的表达式一样

可以观察到： $Tr(T^2) = T_1^2 + T_2^2 + T_3^2=I_T^2-2II_T$ ，（问题1.31）,偏张量的第二不变量公式变为用张量T的第一和第二不变量表示的形式：
$II_{T^{dev}}=\frac{1}{2}[-I_T^2+2II_T+\frac{I_T^2}{3}]=\frac{1}{2}[2II_T-\frac{2I_T^2}{3}]=\frac{1}{3}(3II_T-I_T^2)$

另一种形式是用偏张量分量表示的形式：
$II_{T^{dev}}=-\frac{1}{2}Tr[(T^{dev})^2]=-\frac{1}{2}Tr[(T^{dev}\cdot T^{dev})]=-\frac{1}{2}T^{dev}\cdot \cdot T^{dev}=-\frac{1}{2}T_{ij}^{dev}T_{ji}^{dev}$

展开，得：
$II_{T^{dev}}=-\frac{1}{2}[(T_{11}^{dev})^2+(T_{22}^{dev})^2+(T_{33}^{dev})^2+2(T_{12}^{dev})^2+2(T_{13}^{dev})^2+2(T_{23}^{dev})^2]$

变换一下：
$T_{11}^{dev})^2+(T_{22}^{dev})^2+(T_{33}^{dev})^2=-2II_{T^{dev}}-2(T_{12}^{dev})^2-2(T_{13}^{dev})^2-2(T_{23}^{dev})^2$

另外，在主空间的分量：
$II_{T^{dev}}=-\frac{1}{2}T_{ij}^{dev}T_{ji}^{dev}=-\frac{1}{2}[(T_{1}^{dev})^2+(T_{2}^{dev})^2+(T_{3}^{dev})^2]$

或者：
在这里插入图片描述
或者：

将 $II_{T^{dev}}=-\frac{1}{2}T_{ij}^{dev}T_{ji}^{dev}=-\frac{1}{2}[(T_{1}^{dev})^2+(T_{2}^{dev})^2+(T_{3}^{dev})^2]$ 代入到以上式子：
$II_{T^{dev}}=-\frac{1}{6}[(T_{22}^{dev}-T_{33}^{dev})^2+(T_{11}^{dev}-T_{33}^{dev})^2+(T_{11}^{dev}-T_{22}^{dev})^2]-(T_{12}^{dev})^2-(T_{23}^{dev})^2-(T_{13}^{dev})^2$

如果在主空间，则：
$II_{T^{dev}}=-\frac{1}{6}[(T_{2}^{dev}-T_{3}^{dev})^2+(T_{1}^{dev}-T_{3}^{dev})^2+(T_{1}^{dev}-T_{2}^{dev})^2]$

偏张量的第三不变量

偏张量的第三部变量：
在这里插入图片描述
另一个形式:
$III_{T^{dev}}=T_{1}^{dev}T_{2}^{dev}T_{3}^{dev}=\frac{1}{3}T_{ij}^{dev}T_{jk}^{dev}T_{kl}^{dev}$

问题1.39 $\sigma$ 是对称二阶张量， $\sigma^{dev}$ 是一个偏张量，证明： $s:\frac{\partial s}{\partial \sigma}=s$ ，并证明 $\sigma$ 和 $\sigma^{dev}$ 是同轴张量

在这里插入图片描述
参考教材：
Eduardo W.V. Chaves， Notes On Continuum Mechanics

LAWKAWAI

关注

0
点赞
踩
2

收藏

觉得还不错? 一键收藏
0
评论
【连续介质力学】张量的偏微分、球张量和偏张量

张量的偏微分张量的一阶微分，定义：∂A∂A=A,A=∂Aij∂Akl(e^i⨂e^j⨂e^k⨂e^l)=δikδjl(e^i⨂e^j⨂e^k⨂e^l)=I\frac{\partial A}{\partial A} = A_{,A}=\frac{\partial A_{ij}}{\partial A_{kl}}(\hat e_i\bigotimes\hat e_j\bigotimes\hat e_k\bigotimes\hat e_l )\\=\delta_{ik}\delta_{jl}(\hat e_i
复制链接

扫一扫