【连续介质力学】张量分量的变换定律

news2024/11/26 7:51:09

张量分量的变换定律

张量的分量是依赖于坐标系的,所以当坐标系发生旋转,张量分量也会发生改变,张量分量与坐标系通过分量变换规律互相关联起来的。
在这里插入图片描述

考虑在正交基 ( e ^ 1 , e ^ 2 , e ^ 3 ) (\hat e_1, \hat e_2, \hat e_3) (e^1,e^2,e^3)的坐标系 ( x 1 , x 2 , x 3 ) (x_1, x_2, x_3) (x1,x2,x3),任意向量可以表示成:
v ⃗ = v i e ^ i = v 1 e ^ 1 + v 2 e ^ 2 + v 3 e ^ 3 \vec v = v_i\hat e_i = v_1\hat e_1 + v_2\hat e_2 + v_3\hat e_3 v =vie^i=v1e^1+v2e^2+v3e^3
在这里插入图片描述
矩阵形式:
在这里插入图片描述
考虑一个新的坐标系 ( x 1 ′ , x 2 ′ , x 3 ′ ) (x_1', x_2', x_3') (x1,x2,x3), 基于正交基 ( e ^ 1 ′ , e ^ 2 ′ , e ^ 3 ′ ) (\hat e_1', \hat e_2', \hat e_3') (e^1,e^2,e^3),在这个新的坐标系种,向量可以表示成 v ⃗ = v j ′ e ^ j ′ \vec v = v'_j \hat e_j' v =vje^j, 所以:
v ⃗ = v k ′ e ^ k ′ = v j e ^ j \vec v = v_k'\hat e_k' = v_j \hat e_j v =vke^k=vje^j

为了计算得到在某个坐标系的张量的分量,只需要将张量和坐标基进行点积即可
在这里插入图片描述
或者矩阵形式:
在这里插入图片描述
在这里插入图片描述
所以,坐标变换张量的分量的计算公式
a i j = e ^ j ⋅ e ^ i ′ = e ^ i ′ ⋅ e ^ j a_{ij} = \hat e_j \cdot \hat e_i' = \hat e_i' \cdot \hat e_j aij=e^je^i=e^ie^j

指标形式:
v i ′ = a i j v j \boxed{v_i'= a_{ij}v_j} vi=aijvj

坐标变换矩阵:张量在原坐标系变换到新坐标系的变换矩阵,新坐标基控制矩阵的行
在这里插入图片描述
矩阵不是对称的

由于 e ^ i ′ ⋅ e ^ j = ∣ ∣ e ^ i ′ ∣ ∣ ∣ ∣ e ^ j ∣ ∣ cos ⁡ ( x i ′ , x j ) \hat e_i' \cdot \hat e_j = ||\hat e_i'|| ||\hat e_j|| \cos(x_i', x_j) e^ie^j=∣∣e^i∣∣∣∣e^j∣∣cos(xi,xj),所以可以把变换矩阵表示成方向余弦的形式:
在这里插入图片描述
方向余弦的角度是向量与原坐标系的角度
cos ⁡ α 1 = cos ⁡ ( x 1 ′ , x 1 ) \cos \alpha_1 = \cos (x_1', x_1) cosα1=cos(x1,x1)
cos ⁡ β 1 = cos ⁡ ( x 1 ′ , x 2 ) \cos \beta _1= \cos (x_1', x_2) cosβ1=cos(x1,x2)
cos ⁡ γ 1 = cos ⁡ ( x 1 ′ , x 3 ) \cos \gamma_1 = \cos (x_1', x_3) cosγ1=cos(x1,x3)
在这里插入图片描述

上面已经讨论了向量在新坐标系的投影,下面是向量在原坐标系的投影
在这里插入图片描述
因此: e ^ i = a j i e ^ j ′ \boxed{\hat e_i = a_{ji} \hat e_{j}'} e^i=ajie^j

逆变换矩阵:
A − 1 v ⃗ ′ = A − 1 A v ⃗    ⟹    v = A − 1 v ⃗ ′ A^{-1}\vec v' = A^{-1}A \vec v \implies v = A^{-1}\vec v' A1v =A1Av v=A1v

由于 v = A − 1 v ⃗ ′ v = A^{-1}\vec v' v=A1v v ⃗ = A T v ⃗ ′ \vec v = A^T \vec v' v =ATv ,可以得出 A A A 是正交矩阵

A − 1 = A T    ⟹    A T A = 1 → a k i a k j = δ i j A^{-1} = A^T \implies A^T A = 1 \rightarrow a_{ki}a_{kj} = \delta_{ij} A1=ATATA=1akiakj=δij

二阶张量
从原坐标系的正交基 e ^ i \hat e_i e^i 变换到新坐标系得正交基 e ^ i ′ \hat e_i' e^i,遵循变换定律:
e ^ k = a i k e ^ i ′ \hat e_k = a_{ik} \hat e_i' e^k=aike^i

这可以将二阶张量表示为:
在这里插入图片描述
其中:
T i j ′ = T k l a i k a j l = a i k T k l a j l → T ′ = A T A ′ T_{ij}' = T_{kl}a_{ik}a_{jl}=a_{ik}T_{kl}a_{jl} \rightarrow T' = A T A' Tij=Tklaikajl=aikTklajlT=ATA

三阶张量
可以将三阶张量表示为正交基 e ^ i \hat e_i e^i 和正交基 e ^ i ′ \hat e_i' e^i:
在这里插入图片描述
所以,三阶张量在新正交基的分量:
S i j k ′ = S l m n a i l a j m a k n S'_{ijk} = S_{lmn}a_{il}a_{jm}a_{kn} Sijk=Slmnailajmakn

根据张量阶数总结变换定律
在这里插入图片描述

问题1.27 给定 T ′ = A ⋅ T ⋅ A ′ T' = A \cdot T \cdot A' T=ATA, 计算 T ′ T' T的分量

思路: 先将其转换成指标形式,然后为了求二阶张量的分量表示,需要同时作用双缩并
: ( e ^ i ⨂ e ^ j ) :(\hat e_i \bigotimes \hat e_j) :(e^ie^j)
在这里插入图片描述
在这里插入图片描述

问题1.128 T T T是一个对称二阶张量, I T , I I T , I I I T I_T, II_T, III_T IT,IIT,IIIT是标量,其中 I T = T r ( T ) = T i i , I I T = 1 2 [ I T 2 − T r ( T 2 ) ] , I I I T = det ⁡ T I_T = Tr(T) = T_{ii}, \quad II_T = \frac{1}{2}[I_T^2-Tr(T^2)], \quad III_T = \det T IT=Tr(T)=Tii,IIT=21[IT2Tr(T2)],IIIT=detT, 证明 I T , I I T , I I I T I_T, II_T, III_T IT,IIT,IIIT是做坐标基变换的不变量

在这里插入图片描述
其中,用到正交张量的定义:
由于A是正交张量,所以有: a i k a i l = δ k l → A ⋅ A T = 1 a_{ik}a_{il} = \delta_{kl} \rightarrow A\cdot A^T = 1 aikail=δklAAT=1
其中,还用到张量的迹:
A : B = T r ( A ⋅ B T ) A:B = Tr(A \cdot B^T) A:B=Tr(ABT)
由于 T ′ T' T 是正交张量,所以 T ′ = ( T ′ ) T T' = (T')^T T=(T)T,所以
T r ( T ′ ⋅ T ′ ) = T r ( T ′ ⋅ ( T ′ ) T ) = T ′ : T ′ Tr(T' \cdot T') = Tr(T' \cdot (T')^T) = T' : T' Tr(TT)=Tr(T(T)T)=T:T

所以,目前提到了三个不变量
张量的迹;
张量的迹的平方 - 张量平方的迹;
张量的行列式

有四个坐标系
( x 1 , x 2 , x 3 ) (x_1, x_2, x_3) (x1,x2,x3)
( x 1 ′ , x 2 ′ , x 3 ′ ) (x_1', x_2', x_3') (x1,x2,x3)
( x 1 ′ ′ , x 2 ′ ′ , x 3 ′ ′ ) (x_1'', x_2'', x_3'') (x1′′,x2′′,x3′′)
( x 1 ′ ′ ′ , x 2 ′ ′ ′ , x 3 ′ ′ ′ ) (x_1''', x_2''', x_3''') (x1′′′,x2′′′,x3′′′)

考虑以下变换矩阵
A A A: 从 ( x 1 , x 2 , x 3 ) (x_1, x_2, x_3) (x1,x2,x3) 变换到 ( x 1 ′ , x 2 ′ , x 3 ′ ) (x_1', x_2', x_3') (x1,x2,x3)
B B B: 从 ( x 1 ′ , x 2 ′ , x 3 ′ ) (x_1', x_2', x_3') (x1,x2,x3) 变换到 ( x 1 ′ ′ , x 2 ′ ′ , x 3 ′ ′ ) (x_1'', x_2'', x_3'') (x1′′,x2′′,x3′′)
C: 从 ( x 1 ′ ′ , x 2 ′ ′ , x 3 ′ ′ ) (x_1'', x_2'', x_3'') (x1′′,x2′′,x3′′) 变换到 ( x 1 ′ ′ ′ , x 2 ′ ′ ′ , x 3 ′ ′ ′ ) (x_1''', x_2''', x_3''') (x1′′′,x2′′′,x3′′′)

在这里插入图片描述
如果有一个向量 v ⃗ \vec v v , 在每个坐标系之间的转换如下所示:

( x 1 , x 2 , x 3 ) (x_1, x_2, x_3) (x1,x2,x3) 变换到 ( x 1 ′ , x 2 ′ , x 3 ′ ) (x_1', x_2', x_3') (x1,x2,x3)
v ⃗ ′ = A v ⃗ \vec v' = A\vec v v =Av
( x 1 ′ , x 2 ′ , x 3 ′ ) (x_1', x_2', x_3') (x1,x2,x3) 变换到 ( x 1 , x 2 , x 3 ) (x_1, x_2, x_3) (x1,x2,x3)
v ⃗ = A T v ⃗ ′ \vec v = A^T\vec v' v =ATv

( x 1 ′ , x 2 ′ , x 3 ′ ) (x_1', x_2', x_3') (x1,x2,x3) 变换到 ( x 1 ′ ′ , x 2 ′ ′ , x 3 ′ ′ ) (x_1'', x_2'', x_3'') (x1′′,x2′′,x3′′):
v ⃗ ′ ′ = B v ⃗ ′ \vec v'' = B \vec v' v ′′=Bv
( x 1 ′ ′ , x 2 ′ ′ , x 3 ′ ′ ) (x_1'', x_2'', x_3'') (x1′′,x2′′,x3′′)变换到 ( x 1 ′ , x 2 ′ , x 3 ′ ) (x_1', x_2', x_3') (x1,x2,x3) :
v ⃗ ′ = B T v ⃗ ′ ′ \vec v' = B^T \vec v'' v =BTv ′′
代入 v ⃗ ′ = A v ⃗ \vec v' = A\vec v v =Av , 从 ( x 1 , x 2 , x 3 ) (x_1, x_2, x_3) (x1,x2,x3) 变换到 ( x 1 ′ ′ , x 2 ′ ′ , x 3 ′ ′ ) (x_1'', x_2'', x_3'') (x1′′,x2′′,x3′′) :
v ⃗ ′ ′ = B v ⃗ ′ = B A v ⃗ \vec v'' = B \vec v' = BA\vec v v ′′=Bv =BAv

( x 1 ′ ′ , x 2 ′ ′ , x 3 ′ ′ ) (x_1'', x_2'', x_3'') (x1′′,x2′′,x3′′) 变换到 ( x 1 , x 2 , x 3 ) (x_1, x_2, x_3) (x1,x2,x3) :
v ⃗ = A T B T v ⃗ ′ ′ \vec v = A^T B^T\vec v'' v =ATBTv ′′

( x 1 , x 2 , x 3 ) (x_1, x_2, x_3) (x1,x2,x3) 变换到 ( x 1 ′ ′ ′ , x 2 ′ ′ ′ , x 3 ′ ′ ′ ) (x_1''', x_2''', x_3''') (x1′′′,x2′′′,x3′′′) :
v ⃗ ′ ′ ′ = C B A v ⃗ \vec v''' = CBA\vec v v ′′′=CBAv
( x 1 ′ ′ ′ , x 2 ′ ′ ′ , x 3 ′ ′ ′ ) (x_1''', x_2''', x_3''') (x1′′′,x2′′′,x3′′′)变换到 ( x 1 , x 2 , x 3 ) (x_1, x_2, x_3) (x1,x2,x3) :
v ⃗ = A T B T C T v ⃗ ′ ′ ′ \vec v = A^TB^TC^T\vec v''' v =ATBTCTv ′′′

二维分量变换定律

考虑两个坐标系:
在这里插入图片描述
从(x, y) 到(x’, y’):
v ⃗ ′ = A v ⃗ \vec v' = A\vec v v =Av
在这里插入图片描述
利用三角恒等式:
α x ′ x = α y ′ y    ⟹    cos ⁡ ( α x ′ x ) = cos ⁡ ( α y ′ y ) = cos ⁡ ( α ) cos ⁡ ( α x ′ y ) = cos ⁡ ( π 2 − α ) = sin ⁡ ( α ) cos ⁡ ( α y ′ x ) = cos ⁡ ( π 2 + α ) = − sin ⁡ ( α ) \alpha_{x'x}=\alpha_{y'y} \implies \cos(\alpha_{x'x}) = \cos(\alpha_{y'y})=\cos(\alpha) \\ \cos(\alpha_{x'y}) = \cos(\frac{\pi}{2}-\alpha) = \sin(\alpha) \\ \cos(\alpha_{y'x}) = \cos(\frac{\pi}{2}+\alpha) = -\sin(\alpha) αxx=αyycos(αxx)=cos(αyy)=cos(α)cos(αxy)=cos(2πα)=sin(α)cos(αyx)=cos(2π+α)=sin(α)

所以,二维的变换矩阵:
在这里插入图片描述
另一个证明方法

考虑在两个坐标系的一个点P的向量方向:
在这里插入图片描述
在这里插入图片描述

矩阵形式:
在这里插入图片描述
又因为 A − 1 = A A^{-1} = A A1=A, 所以:
在这里插入图片描述

问题1.29 求出坐标系 ( x , y , z ) (x, y, z) (x,y,z) ( x ′ ′ ′ , y ′ ′ ′ , z ′ ′ ′ ) (x''', y''', z''') (x′′′,y′′′,z′′′)之间的变换矩阵

在这里插入图片描述
![
在这里插入图片描述
在这里插入图片描述
α , β , γ \alpha, \beta, \gamma α,β,γ 是欧拉角,用来计算刚体运动的方向

问题1.30 T T T是二阶张量

在这里插入图片描述

NOTE: 以上例子可以看出,张量 T T T 在新坐标基的分量有一个特性, 非对角元素是零,现在问题是:给定任意一个张量 T T T ,是否存在使得非对角元素为零的变换?
这种问题的类型是特征值问题。

参考教材:
Eduardo W.V. Chaves, Notes On Continuum Mechanics

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/525426.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

【51单片机】串口通信使用串口通信控制LED灯

🎊专栏【51单片机】 🍔喜欢的诗句:更喜岷山千里雪 三军过后尽开颜。 🎆音乐分享【51单片机】 🥰大一同学小吉,欢迎并且感谢大家指出我的问题🥰 目录 🍔串口通信 😎代码…

HTTP1.1(四)URI

一 URI的基本格式以及与URL的区别 ① URI直观感受 1)URI是HTTP协议的重要组成部分2)ABNF如何严格定义URI ② URL、URN、URI ③ URI细讲 ④ URI的组成 ⑤ ABNF关于URI定义 blob:http 地址的原理及生成方法 video 标签 src 带有blob:http data…

NetCore发布到Docker-单机

一.git安装 1.下载 https://github.com/git/git/releases/tag/v2.40.1 2.将本地的安装包上传到 linux 服务器上,我这里放在 /opt/git/ 目录下 3.解压压缩包,得到目录 git-2.40.1,位置在 /opt/git/git-2.40.1 tar -zxvf git-2.40.1.tar.gz 4.…

【刷题之路】LeetCode 225. 用队列实现栈

【刷题之路】LeetCode 225. 用队列实现栈 一、题目描述二、解题1、主要思路解析2、先实现栈3、实现各个接口3.1、初始化接口3.2、push接口3.3、pop接口3.4、myStackTop接口3.5、myStackEmpty接口3.6、myStackFree接口 一、题目描述 原题连接: 225. 用队列实现栈 题…

如何选择合适的自动化测试工具? 我教你

目录 前言 选择最佳自动化测试工具的重要性 自动化测试工具的种类 选择自动化测试工具时要考虑的因素 自动化测试工具的选择标准 前言 自动化测试是高质量软件交付领域中最重要的实践之一。在今天的敏捷开发方法中,几乎任一软件开发过程都需要在开发阶段的某个…

测试工具-JIRA

软件测试工具—JIRA使用 1.创建 项目:一般为业务线名称 问题类型:任务、故事、缺陷等 概要:bug的简要描述,一般写出定位地方和目前现状 Sprint:一般为业务线迭代简称 测试阶段:按照各公司的测试流程对应标…

vlc搭建rtsp直播Demo

文章目录 学习链接本地视频文件作为数据源推流步骤拉流步骤 本地摄像头作为数据源拉流步骤 学习链接 ffmpeg 使用VLC亲手搭建RTSP直播点播 本地视频文件作为数据源 推流步骤 点击媒体->流 点击添加,选择视频文件,后点击串流 点击串流后&#xf…

linux异步通知实验

一、异步通知简介 中断是处理器提供的一种异步机制,配置好中断后就可以让处理器去处理其他的事情,当中断发生以后会执行中断服务函数,在中断服务函数中做具体的处理。 Linux 应用程序可以通过阻塞或者非阻塞两种方式来访问驱动设备&#xff0…

04:mysql--DCL

目录 1:介绍 2:语法 3:DCL的使用 4:权限 1:介绍 DCL英文全称是Data ControlLanguage(数据控制语言),用来管理数据库 用户、控制数据库的访问 权限。 2:语法 DCL-管理用户 查询用户 use mysql; select * from user; 创建用户 create user "用户名&quo…

缓存服务器概述

一、概述 主要功能:应用解耦,异步消息,流量削锋等问题 架构设计:实现高性能,高可用,可伸缩和最终一致性架构 常用消息队列:ActiveMQ,RabbitMQ,ZeroMQ,Kafka…

谷歌chrome浏览器升级新版后字体显示不清楚解决方案

谷歌chrome浏览器升级新版后字体显示不清楚解决方案 参考图片: Chrome更新至版本Chrome 109.0.5414.120 字体看不清 浏览器症状与表现 Chrome更新至版本Chrome 109.0.5414.120 字体看不清;会很细,在设置中选择自定义的字体,仍无法…

MySQL数据库笔记-基础篇

文章目录 MySQL概述数据库相关概念小总结 MySQL数据库版本:启动与停止客户端连接关系型数据库数据模型小总结 SQLSQL通用语法SQL分类DDL数据库操作查询创建数据类型数据类型的使用练习 修改删除小总结 DML添加数据修改数据删除数据小总结 DQLDCL MySQL概述 数据库相…

BI技巧丨矩阵甘特图

BOSS:白茶,我们最近有很多项目,能不能做个甘特图? 白茶:没问题!(应用市场上随便找个甘特图就OK了嘛!) BOSS:对了,不要那种点进来会有点卡&#xf…

二、数据结构1:单链表 模板题+算法模板(单链表)

文章目录 算法模板单链表题目模板 模板题单链表原题链接题目思路题解 算法模板 单链表题目模板 // head存储链表头,e[]存储节点的值,ne[]存储节点的next指针,idx表示当前用到了哪个节点 int head, e[N], ne[N], idx;// 初始化 void init() …

【001】C++的第一个程序Hello World

C的第一个应用程序(Hello World程序) 引言一、代码二、代码解释三、注意事项四、注释语句总结 引言 💡 作者简介:专注于C/C高性能程序设计和开发,理论与代码实践结合,让世界没有难学的技术。 👉…

《JavaEE》网络编程套接字

文章目录 Socket我们来解释一下上面叫做有无连接我们再来理解一下上面是可靠传输和不可靠传输面向字节流与面向数据报的区别(后期会具体分析这个)全双工 单双工 UDPDatagramSocketDatagramPacket我们来尝试写一下UDP版本的代码 TCPServerSocketSocket我们…

什么年代了?不会还有人不会插件化吧?

一.到底什么是插件化开发 插件化开发是将整个app拆分成很多模块,这些模块包括一个宿主和多个插件,每个模块都是一个apk,最终打包的时候将宿主apk和插件apk分开或者联合打包。 为什么要将一定要使用插件化开发呢? 宿…

【Linux】一文带你掌握Linux权限!

Linux权限 问题一:Linux下的用户分类有哪些?问题二:什么叫做权限?问题三:见一见没有权限是什么现象(看看就好了)问题四:权限的修改问题。chmodchownchgrp 问题五:其他问题(三个小问题)文件的初始…

51单片机(十五)直流电机驱动(PWM)

❤️ 专栏简介:本专栏记录了从零学习单片机的过程,其中包括51单片机和STM32单片机两部分;建议先学习51单片机,其是STM32等高级单片机的基础;这样再学习STM32时才能融会贯通。 ☀️ 专栏适用人群 :适用于想要…

【并发编程】深入理解Java并发之synchronized实现原理

文章目录 一、synchronized 实现同步示例二、synchronized 典型错误示例三、Java 对象头与 Monitor四、synchronized代码块底层原理五、synchronized方法底层原理六、Java虚拟机对synchronized的优化 一、synchronized 实现同步示例 public class MyThread extends Thread{pri…