【连续介质力学】张量的范数、各向同性和各向异性张量、同轴张量和极分解

news2024/11/23 11:14:06

张量的范数

张量的大小,使用Frobenius 范数:
∣ ∣ v ⃗ ∣ ∣ = v ⃗ ⋅ v ⃗ = v i v i (向量) ||\vec v|| = \sqrt{\vec v \cdot \vec v} = \sqrt{v_iv_i} (向量) ∣∣v ∣∣=v v =vivi (向量)
∣ ∣ T ∣ ∣ = T : T = T i j T i j (二阶张量) ||T|| = \sqrt{ T:T} = \sqrt{T_{ij}T_{ij}} (二阶张量) ∣∣T∣∣=T:T =TijTij (二阶张量)
∣ ∣ A ∣ ∣ = A : ⋅ A = A i j k A i j k (三阶张量) ||A|| = \sqrt{ A:\cdot A} = \sqrt{A_{ijk}A_{ijk}} (三阶张量) ∣∣A∣∣=A:A =AijkAijk (三阶张量)
∣ ∣ C ∣ ∣ = C : : C = C i j k l C i j k l (四阶张量) ||C|| = \sqrt{ C:: C} = \sqrt{C_{ijkl}C_{ijkl}} (四阶张量) ∣∣C∣∣=C::C =CijklCijkl (四阶张量)

在主空间,张量的特征值 T 1 , T 2 , T 3 T_1, T_2, T_3 T1,T2,T3,在主空间中:
∣ ∣ T ∣ ∣ = T : T = T i j T i j = T 1 2 + T 2 2 + T 3 3 = I T 2 − 2 I I T ||T|| = \sqrt{ T:T} = \sqrt{T_{ij}T_{ij}} =\sqrt{T_1^2+T_2^2+T_3^3}=\sqrt{I_T^2-2II_T} ∣∣T∣∣=T:T =TijTij =T12+T22+T33 =IT22IIT

所以, ∣ ∣ T ∣ ∣ ||T|| ∣∣T∣∣是一个不变量, ∣ ∣ T ∣ ∣ ||T|| ∣∣T∣∣表示的在主空间下主方向的长度
在这里插入图片描述

各向同性和各向异性张量

各向同性:在任意坐标系下,张量的分量都是一样的, T = T ′ T = T' T=T

各向同性一阶张量

在坐标系 ( x 1 , x 2 , x 3 ) (x_1, x_2, x_3) (x1,x2,x3) 的分量 ( v 1 , v 2 , v 3 ) (v_1, v_2, v_3) (v1,v2,v3)
在坐标系 ( x 1 ′ , x 2 ′ , x 3 ′ ) (x_1', x_2', x_3') (x1,x2,x3) 的分量 ( v 1 ′ , v 2 ′ , v 3 ′ ) (v_1', v_2', v_3') (v1,v2,v3)
根据变换定律:
v ⃗ = v i e ^ i = v j ′ e ^ j ′    ⟹    v i ′ = a i j v j \vec v = v_i \hat e_i = v_j' \hat e_j' \implies v_i' = a_{ij}v_j v =vie^i=vje^jvi=aijvj
由各向同性的定义, v i = v j ′ v_i = v_j' vi=vj,那么有: e ^ i = e ^ j ′ \hat e_i = \hat e_j' e^i=e^j

所以,要么是坐标系根本没有变换,要么满足各向同性的一阶张量只能是零向量 0 ⃗ \vec 0 0

各向同性二阶张量

例子:单位张量 1 1 1, δ k l \delta_{kl} δkl,根据二阶张量的变换定律:
δ i j ′ = a i k a j l δ k l = a i k a j k = δ i j \delta_{ij}' = a_{ik}a_{jl}\delta_{kl}=a_{ik}a_{jk}=\delta_{ij} δij=aikajlδkl=aikajk=δij

所以,如果一个二阶张量是各向同性的,那么这个张量是球形张量

各向同性三阶张量
Levi-Civita pseudo-tensor: ϵ i j k \epsilon_{ijk} ϵijk 不是真的张量
ϵ i j k ′ = a i l a j m a k n ϵ l m n = ∣ A ∣ ϵ i j k = ϵ i j k \epsilon _{ijk}'=a_{il}a_{jm}a_{kn}\epsilon_{lmn}=|A|\epsilon_{ijk}=\epsilon_{ijk} ϵijk=ailajmaknϵlmn=Aϵijk=ϵijk

各向同性四阶张量

I I ‾ ‾ i j k l = δ i j δ k l \overline{\overline{II}}_{ijkl} = \delta_{ij}\delta_{kl} IIijkl=δijδkl
I I i j k l = δ i k δ j l II_{ijkl} = \delta_{ik}\delta_{jl} IIijkl=δikδjl
I I ‾ i j k l = δ i l δ j k \overline{II}_{ijkl} = \delta_{il}\delta_{jk} IIijkl=δilδjk

所以,任意一个四阶各向同性张量可以表示为以上张量的线性组合:
在这里插入图片描述

问题1.37 四阶张量: C i j k l = λ δ i j δ k l + μ ( δ i k δ j l + δ i l δ j k ) C_{ijkl}=\lambda\delta_{ij}\delta_{kl}+\mu (\delta_{ik}\delta_{jl} + \delta_{il} \delta_{jk}) Cijkl=λδijδkl+μ(δikδjl+δilδjk),证明C是各向同性

在这里插入图片描述

同轴张量 Coaxial tensor

如果张量 T T T S S S有相同的特征向量,那么它们是同轴张量,也就是它们之间点积是可以交换的:
T ⋅ S = S ⋅ T    ⟹    S , T 是同轴的 T \cdot S = S \cdot T \implies S,T 是同轴的 TS=STS,T是同轴的

如果张量 T T T S S S同轴且对称的,那么它们的谱表示:
T = ∑ a = 1 3 T a n ^ ( a ) ⨂ n ^ ( a ) ; S = ∑ a = 1 3 S a n ^ ( a ) ⨂ n ^ ( a ) ; T= \sum_{a=1}^3 T_a\hat n^{(a)}\bigotimes \hat n^{(a)} ; \quad S= \sum_{a=1}^3 S_a\hat n^{(a)}\bigotimes \hat n^{(a)} ; T=a=13Tan^(a)n^(a);S=a=13San^(a)n^(a);

S S S S − 1 S^{-1} S1是同轴的:
S − 1 ⋅ S = S ⋅ S − 1 = 1 S^{-1}\cdot S = S \cdot S^{-1} = 1 S1S=SS1=1
S = ∑ a = 1 3 S a n ^ ( a ) ⨂ n ^ ( a ) ; S − 1 = ∑ a = 1 3 1 S a n ^ ( a ) ⨂ n ^ ( a ) ; S= \sum_{a=1}^3 S_a\hat n^{(a)}\bigotimes \hat n^{(a)} ; \quad S^{-1}= \sum_{a=1}^3 \frac{1}{S_a}\hat n^{(a)}\bigotimes \hat n^{(a)} ; S=a=13San^(a)n^(a);S1=a=13Sa1n^(a)n^(a);

如果 S S S T T T是同轴且对称的张量,那么 S ⋅ T S \cdot T ST的结果是另一个对称张量:
T ⋅ S = S ⋅ T    ⟹    T ⋅ S − S ⋅ T = 0    ⟹    T ⋅ S − ( T ⋅ S ) T = 0    ⟹    2 ( T ⋅ S ) s k e e w = 0 T\cdot S = S \cdot T \\ \implies T \cdot S -S \cdot T = 0 \\ \implies T \cdot S -(T \cdot S)^T = 0 \\ \implies 2(T\cdot S)^{skeew} = 0 TS=STTSST=0TS(TS)T=02(TS)skeew=0

所以张量的反对称部分是零张量,那么张量只由对称张量组成:
( T ⋅ S ) ≡ ( T ⋅ S ) s y m (T \cdot S) \equiv (T \cdot S)^{sym} (TS)(TS)sym

极分解

F F F是非奇异的二阶张量, det ⁡ F ≠ 0    ⟹    ∃ F − 1 \det F \neq 0 \implies \exists F^{-1} detF=0F1

张量满足: F ⋅ N ^ = f ⃗ ( N ^ ) = ∣ ∣ f ⃗ ( N ^ ) ∣ ∣ n ^ = λ n ^ n ^ ≠ 0 ⃗ F \cdot \hat N = \vec f^{(\hat N)}=||\vec f^{(\hat N)}||\hat n = \lambda_{\hat n}\hat n \neq \vec 0 FN^=f (N^)=∣∣f (N^)∣∣n^=λn^n^=0

给定一个正交基 N ^ ( a ) \hat N^{(a)} N^(a), 可得:
在这里插入图片描述
以上 F F F 的表达式并不是 F F F 的谱表示,因为 λ a \lambda_a λa不是 F F F的特征值,并且 n ^ ( a ) \hat n^{(a)} n^(a) N ^ ( a ) \hat N^{(a)} N^(a)也都不是 F F F的特征向量

在这里插入图片描述
可以发现,对于任意正交基 N ^ ( a ) \hat N^{(a)} N^(a),新基 n ^ ( a ) \hat n^{(a)} n^(a) 不一定是正交的

现在,我们希望找到一个基 N ^ ( a ) \hat N^{(a)} N^(a),使得新基 n ^ ( a ) \hat n^{(a)} n^(a)是正交的,也就是:
f ⃗ ( N ^ ( 1 ) ) ⋅ f ⃗ ( N ^ ( 2 ) ) = 0 f ⃗ ( N ^ ( 2 ) ) ⋅ f ⃗ ( N ^ ( 3 ) ) = 0 f ⃗ ( N ^ ( 3 ) ) ⋅ f ⃗ ( N ^ ( 1 ) ) = 0 \vec f^{(\hat N^{(1)})} \cdot \vec f^{(\hat N^{(2)})} = 0 \\ \vec f^{(\hat N^{(2)})} \cdot \vec f^{(\hat N^{(3)})} = 0 \\ \vec f^{(\hat N^{(3)})} \cdot \vec f^{(\hat N^{(1)})} = 0 f (N^(1))f (N^(2))=0f (N^(2))f (N^(3))=0f (N^(3))f (N^(1))=0

那么我们根据正交变换 n ^ ( a ) = R ⋅ N ^ ( a ) \hat n^{(a)} = R \cdot \hat N^{(a)} n^(a)=RN^(a) 来寻找一个空间,可以保证 n ^ ( a ) \hat n^{(a)} n^(a) 的正交性,因为正交变换既不改变向量之间的角度,也不改变其大小

所以,假设有一个从 N ^ ( a ) \hat N^{(a)} N^(a) n ^ ( a ) \hat n^{(a)} n^(a) 的变换,正交变换 n ^ ( a ) = R ⋅ N ^ ( a ) \hat n^{(a)} = R \cdot \hat N^{(a)} n^(a)=RN^(a),那么:
F = ∑ a = 1 3 λ a n ^ ( a ) ⨂ N ^ ( a ) = ∑ a = 1 3 λ a R ⋅ N ^ ( a ) ⨂ N ^ ( a ) = R ⋅ ∑ a = 1 3 λ a N ^ ( a ) ⨂ N ^ ( a ) = R ⋅ U F = \sum_{a =1}^3\lambda_a \hat n^{(a)}\bigotimes \hat N^{(a)} =\sum_{a =1}^3\lambda_a R \cdot \hat N^{(a)}\bigotimes \hat N^{(a)} \\ =R \cdot\sum_{a =1}^3\lambda_a \hat N^{(a)}\bigotimes \hat N^{(a)} = R \cdot U F=a=13λan^(a)N^(a)=a=13λaRN^(a)N^(a)=Ra=13λaN^(a)N^(a)=RU

F = R ⋅ U    ⟹    U = R T ⋅ F F = R \cdot U \implies U = R^T \cdot F F=RUU=RTF

其中, U = N ^ ( a ) ⨂ N ^ ( a ) U = \hat N^{(a)}\bigotimes \hat N^{(a)} U=N^(a)N^(a),是一个对称张量, U = U T U = U^T U=UT

反过来, N ^ ( a ) = R T ⋅ n ^ ( a ) = n ^ ( a ) ⋅ R \hat N^{(a)} = R^T \cdot \hat n^{(a)} = \hat n^{(a)} \cdot R N^(a)=RTn^(a)=n^(a)R, 那么:
F = ∑ a = 1 3 λ a n ^ ( a ) ⨂ N ^ ( a ) = ∑ a = 1 3 λ a n ^ ( a ) ⨂ n ^ ( a ) ⋅ R = V ⋅ R F = \sum_{a =1}^3\lambda_a \hat n^{(a)}\bigotimes \hat N^{(a)} =\sum_{a =1}^3\lambda_a \hat n^{(a)}\bigotimes \hat n^{(a)} \cdot R \\ =V \cdot R F=a=13λan^(a)N^(a)=a=13λan^(a)n^(a)R=VR

F = V ⋅ R    ⟹    V = F ⋅ R T F = V \cdot R \implies V = F \cdot R^T F=VRV=FRT

其中, V = n ^ ( a ) ⨂ n ^ ( a ) V = \hat n^{(a)}\bigotimes \hat n^{(a)} V=n^(a)n^(a)

U U U V V V有相同的特征值,但不同的特征向量

定义极分解
F = R ⋅ U = V ⋅ R T (极分解) \boxed{F = R \cdot U = V \cdot R^T} (极分解) F=RU=VRT(极分解)

F T F^T FT F F F进行点积:
在这里插入图片描述
在这里插入图片描述

NOTE: 由于 det ⁡ F ≠ 0 \det F \neq 0 detF=0,所以 C C C b b b是正定对称张量,表示 C C C b b b的特征值是正实数,但是 det ⁡ F ≠ 0 \det F \neq 0 detF=0 有以下两种情况:

  • 如果 det ⁡ F > 0 \det F > 0 detF>0
    det ⁡ F = det ⁡ R det ⁡ U = det ⁡ V det ⁡ R T > 0 \det F = \det R \det U = \det V \det R^T > 0 detF=detRdetU=detVdetRT>0,那么:
    在这里插入图片描述
  • 如果 det ⁡ F < 0 \det F < 0 detF<0
    det ⁡ F = det ⁡ R det ⁡ U = det ⁡ V det ⁡ R T < 0 \det F = \det R \det U = \det V \det R^T < 0 detF=detRdetU=detVdetRT<0,那么:
    在这里插入图片描述

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

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

相关文章

okhttp篇2:Dispatcher

Dispatchers维护着一个线程池&#xff0c;3个双端队列&#xff0c;准备执行的AsynCall&#xff0c;正在执行的AsynCall&#xff0c;正在执行的同步Call&#xff08;RealCall&#xff09;。 同时规定每个Host最多同时请求5个Request&#xff0c;同时可最多执行64个Request。 p…

玩转Google开源C++单元测试框架Google Test系列(gtest)之一 - 初识gtest

一、前言 本篇将介绍一些gtest的基本使用&#xff0c;包括下载&#xff0c;安装&#xff0c;编译&#xff0c;建立我们第一个测试Demo工程&#xff0c;以及编写一个最简单的测试案例。 二、下载 如果不记得网址&#xff0c; 直接在google里搜gtest&#xff0c;第一个就是。目…

Docker下Gitlab配置Let’s Encrypt证书

Docker下Gitlab配置Let’s Encrypt证书 1 参考文档2 常见问题2.1 前置条件2.2 不支持ip2.3 重复签发2.4 外网无法访问 ※3 内网穿透配置&#xff08;可选&#xff09;4 Gitlab 创建并配置Let’s Encrypt证书4.1 开放Let’s Encrypt签发所需端口4.2 新增存储HTTPS证书文件夹4.3 …

多态与虚函数(补)

多态与虚函数&#xff08;补&#xff09; 静态联编与动态联编的深层次理解多态底层原理 示例示例一示例二示例三示例四 对象与内存虚析构函数构造函数为什么不能是虚函数&#xff1f; 静态联编与动态联编的深层次理解 我们首先看下面一段代码 class object { private: int va…

C learning_12 操作符前篇(算术操作符、移位操作符、位操作符、赋值操作符、单目操作符、关系操作符、逻辑操作符)

目录 算术操作符 移位操作符 移位规则 位操作符 交换两个整形变量的写法 赋值操作符 单目操作符 sizeof和数组的纠缠 和--运算符 多组输入的方案 关系操作符 逻辑操作符 算术操作符 -- 加法操作符&#xff08;&#xff09;&#xff1a;用于将两个值相加。 -- 减法操…

Python爬虫(二):Requests库

所谓爬虫就是模拟客户端发送网络请求&#xff0c;获取网络响应&#xff0c;并按照一定的规则解析获取的数据并保存的程序。要说 Python 的爬虫必然绕不过 Requests 库。 1 简介 对于 Requests 库&#xff0c;官方文档是这么说的&#xff1a; Requests 唯一的一个非转基因的 P…

存储知识点:RAID0、RAID1、RAID5、RAID10特点是什么?所需的硬盘数量分别为多少?

RAID&#xff08;Redundant Array of Independent Disks&#xff09;是一种将多个独立的硬盘组合成一个逻辑磁盘的技术&#xff0c;目的是提高性能或容错能力。RAID有不同的级别&#xff0c;常见的有RAID0、RAID1、RAID5、RAID10等。下面我们来介绍这些级别的特点和所需的硬盘数…

套接字编程简介

作者&#xff1a;V7 博客&#xff1a;https://www.jvmstack.cn 一碗鸡汤 少年辛苦终身事&#xff0c;莫向光阴惰寸功。 —— 杜荀鹤 Socket概述 在计算机中产生和接受IO流的数据源是多种多样的&#xff0c;在网络编程中&#xff0c;有一个特殊的数据源就是socket。通俗点soc…

linux的系统日志

目录 一、日志文件的产生 二、日志文件存放在哪儿 &#xff08;1&#xff09;文本日志 &#xff08;2&#xff09;二进制日志 三、日志存放规则的配置文件 四、日志轮转 五、分析和监控日志 一、日志文件的产生 日志内容&#xff1a;内核、开机引导、守护进程启动运行的…

华为和思科两种常见的网络设备如何进行ospf配置?

概述 ospf&#xff08;开放最短路径优先&#xff09;是一种基于链路状态的动态路由协议&#xff0c;它可以在网络中自动发现和维护最优的路由路径。ospf广泛应用于大型和复杂的网络环境&#xff0c;因为它具有以下优点&#xff1a; 支持分层路由&#xff0c;可以将网络划分为…

WebAssembly黑暗的一面

案例1&#xff1a;技术支持诈骗 什么是技术支持诈骗&#xff1f; 技术支持诈骗是一种电话欺诈&#xff0c;其中诈骗者声称可以提供合法的技术支持服务。该骗局可能以陌生电话开始&#xff0c;骗子通常会声称来自合法的第三方的员工&#xff0c;如“微软”或“Windows部门”。他…

YOLOv5实现目标分类计数并显示在图像上

有同学后台私信我&#xff0c;想用YOLOv5实现目标的分类计数&#xff0c;因此本文将在之前目标计数博客的基础上添加一些代码&#xff0c;实现分类计数。阅读本文前请先看那篇博客&#xff0c;链接如下&#xff1a; YOLOv5实现目标计数_Albert_yeager的博客 1. 分类实现 以co…

web 实验一 HTML基本标签实验

实验原理 通过创建HTML5网页&#xff0c;验证form内多种元素标签及其属性的作用及意义。 实验目的 理解并掌握Form表单提交必须声明的内容 理解并掌握Input元素中多种类型属性的使用方法及使用场景 理解并掌握Label元素的使用方法 理解并掌握Datalist元素的使用方法 理解并掌握…

软件测试学习——笔记一

一、软件和软件测试 1、软件和软件分类 &#xff08;1&#xff09;软件&#xff1a;程序、数据、文档——用户手册 &#xff08;2&#xff09;软件的分类 按层次划分&#xff1a;系统软件、应用软件按组织划分&#xff1a;开源软件&#xff08;代码公开&#xff09;、商业软…

RSA 加密算法在C++中的实现 面向初学者(附代码)

概述 博文的一&#xff0c;二部分为基础知识的铺垫。分别从密码学&#xff0c;数论两个方面为理解RSA算法做好了准备。第三部分是对RSA加密过程的具体介绍&#xff0c;主要涉及其密钥对&#xff08;key-pair&#xff09;的获取。前三个部分与编程实践无关&#xff0c;可以当作…

C# | 内存池

内存池 文章目录 内存池前言什么是内存池内存池的优点内存池的缺点 实现思路示例代码结束语 前言 在上一篇文章中&#xff0c;我们介绍了对象池的概念和实现方式。对象池通过重复利用对象&#xff0c;避免了频繁地创建和销毁对象&#xff0c;提高了系统的性能和稳定性。 今天我…

你真的了解索引吗

当我们学习存储算法和索引算法时&#xff0c;他们可以深入了解如何在系统中存储和查询数据。因为存储和查询数据是许多系统的核心功能之一&#xff0c;例如数据库、搜索引擎等。理解这些算法可以帮助程序员更好地设计和优化系统架构&#xff0c;提高系统的可扩展性、可用性和性…

玩转Google开源C++单元测试框架Google Test系列(gtest)之二 - 断言

一、前言 这篇文章主要总结gtest中的所有断言相关的宏。 gtest中&#xff0c;断言的宏可以理解为分为两类&#xff0c;一类是ASSERT系列&#xff0c;一类是EXPECT系列。一个直观的解释就是&#xff1a; 1. ASSERT_* 系列的断言&#xff0c;当检查点失败时&#xff0c;退出当前…

大数据之光:Apache Spark 实用指南 大数据实战详解【上进小菜猪大数据】

上进小菜猪&#xff0c;沈工大软件工程专业&#xff0c;爱好敲代码&#xff0c;持续输出干货。 本文将深入探讨Apache Spark作为一种强大的大数据处理框架的基本概念、特点和应用。我们将详细介绍Spark的核心组件&#xff0c;包括Spark Core、Spark SQL、Spark Streaming和Spa…

百子作业 —— 中国邮递员问题

题目 严老师和宋老板去勘测武威市区的道路网&#xff0c;每一条路都需要勘测&#xff0c;且需要两人合作.武威市区可以近似地看成六横六纵组成的道路网&#xff0c;自西向东依次为学府路、民勤路、西关路、中关路、富民路、滨河路&#xff1b;自北向南依次为雷海路、宣武路、祁…