【小呆的力学笔记】弹塑性力学的初步认知二:应力应变分析(2)-程序员宅基地

技术标签: 笔记  矩阵  线性代数  有限元  

1.4 主应力空间、八面体应力

一点的应力状态不论如何变化,其主应力和主方向一致的话,该点的应力状态就是唯一确定的。因此,我们用主应力方向建立一个三维坐标系来描述问题将不失一般性,该坐标系如下图4,我们称之为主应力空间。我们考察等倾面组成的八面体,图中O’P点为等倾面ABC上面的应力向量 ( p 1 , p 2 , p 3 ) (p_1,p_2,p_3) (p1,p2,p3),八面体为等倾面八面体,即面ABC的法线方向余弦为 ( 1 3 , 1 3 , 1 3 ) (\frac{1}{\sqrt 3},\frac{1}{\sqrt 3},\frac{1}{\sqrt 3}) (3 1,3 1,3 1)。将O’P分解
O ’ P ‾ = O ’ Q ‾ + O ’ N ‾ (25) \overline {O’P}=\overline {O’Q}+\overline{O’N}\tag{25} OP=OQ+ON(25)

在这里插入图片描述

图 4 八面体 图4八面体 4八面体
取等倾面和三个轴的坐标面组成的四面体为研究对象,如下图5所示。
在这里插入图片描述
图 5 等倾面四面体 图5等倾面四面体 5等倾面四面体
根据斜面应力公式 p j = σ i j n i p_j=\sigma_{ij}n_i pj=σijni,不难得到以下关系式(矩阵形式)
[ p 1 p 2 p 3 ] = [ σ 1 0 0 0 σ 2 0 0 0 σ 2 ] [ n 1 n 2 n 3 ] (26) \begin{bmatrix} p_1 \\ p_2\\p_3 \end{bmatrix}=\begin{bmatrix} \sigma_1 & 0 & 0\\ 0 & \sigma_2 & 0 \\0 & 0 & \sigma_2 \end{bmatrix}\begin{bmatrix} n_1 \\ n_2\\n_3 \end{bmatrix}\tag{26} p1p2p3 = σ1000σ2000σ2 n1n2n3 (26)

其中 ( n 1 , n 2 , n 3 ) = ( 1 3 , 1 3 , 1 3 ) (n_1 ,n_2,n_3)=(\frac{1}{\sqrt 3},\frac{1}{\sqrt 3},\frac{1}{\sqrt 3}) (n1,n2,n3)=(3 1,3 1,3 1)为等倾面的法线方向余弦。
那么,有
σ 8 = [ n 1 n 2 n 3 ] [ p 1 p 2 p 3 ] = σ 1 n 1 2 + σ 2 n 2 2 + σ 3 n 3 2 = 1 3 ( σ 1 + σ 2 + σ 3 ) = 1 3 I 1 (27) \sigma_8 = \begin{bmatrix} n_1 & n_2 & n_3 \end{bmatrix}\begin{bmatrix} p_1 \\ p_2\\p_3 \end{bmatrix}=\sigma_1n_1^2+\sigma_2n_2^2+\sigma_3n_3^2=\frac{1}{3}(\sigma_1+\sigma_2+\sigma_3)=\frac{1}{3}I_1 \tag{27} σ8=[n1n2n3] p1p2p3 =σ1n12+σ2n22+σ3n32=31(σ1+σ2+σ3)=31I1(27)
八面体相应的剪应力为
τ 8 = p 2 − σ 8 2 = p 1 2 + p 2 2 + p 3 2 − ( σ 1 n 1 2 + σ 2 n 2 2 + σ 3 n 3 2 ) 2 = σ 1 2 n 1 2 + σ 2 2 n 2 2 + σ 3 2 n 3 2 − ( σ 1 n 1 2 + σ 2 n 2 2 + σ 3 n 3 2 ) 2 = 1 3 ( σ 1 2 + σ 2 2 + σ 3 2 ) − 1 9 ( σ 1 + σ 2 + σ 3 ) 2 = 1 3 3 ( σ 1 2 + σ 2 2 + σ 3 2 ) − ( σ 1 2 + σ 2 2 + σ 3 2 + 2 σ 1 σ 2 + 2 σ 1 σ 3 + 2 σ 2 σ 3 ) = 1 3 ( σ 1 − σ 2 ) 2 + ( σ 1 − σ 3 ) 2 + ( σ 2 − σ 3 ) 2 = 2 3 J 2 = 1 3 s i j s i j (28) \tau_8 = \sqrt{p^2-\sigma_8^2}=\sqrt{p_1^2+p_2^2+p_3^2-(\sigma_1n_1^2+\sigma_2n_2^2+\sigma_3n_3^2)^2}\\ =\sqrt{\sigma_1^2n_1^2+\sigma_2^2n_2^2+\sigma_3^2n_3^2-(\sigma_1n_1^2+\sigma_2n_2^2+\sigma_3n_3^2)^2}\\ =\sqrt{\frac{1}{3}(\sigma_1^2+\sigma_2^2+\sigma_3^2)-\frac{1}{9}(\sigma_1+\sigma_2+\sigma_3)^2}\\ =\frac{1}{3}\sqrt{3(\sigma_1^2+\sigma_2^2+\sigma_3^2)-(\sigma_1^2+\sigma_2^2+\sigma_3^2+2\sigma_1\sigma_2+2\sigma_1\sigma_3+2\sigma_2\sigma_3)}\\ =\frac{1}{3}\sqrt{(\sigma_1-\sigma_2)^2+(\sigma_1-\sigma_3)^2+(\sigma_2-\sigma_3)^2}=\sqrt{\frac{2}{3}J_2}=\sqrt{\frac{1}{3}s_{ij}s_{ij}} \tag{28} τ8=p2σ82 =p12+p22+p32(σ1n12+σ2n22+σ3n32)2 =σ12n12+σ22n22+σ32n32(σ1n12+σ2n22+σ3n32)2 =31(σ12+σ22+σ32)91(σ1+σ2+σ3)2 =313(σ12+σ22+σ32)(σ12+σ22+σ32+2σ1σ2+2σ1σ3+2σ2σ3) =31(σ1σ2)2+(σ1σ3)2+(σ2σ3)2 =32J2 =31sijsij (28)

1.5 应变分析

应变分析的内容同应力分析内容,只是注意一点,应变张量和工程应变在剪应变分量是不同的,定义如下。
[ ε x x ε y x ε z x ε x y ε y y ε z y ε x z ε y z ε z z ] = [ ε x x 1 2 γ y x 1 2 γ z x 1 2 γ x y ε y y 1 2 γ z y 1 2 γ x z 1 2 γ y z ε z z ] (29) \begin{bmatrix} \varepsilon_{xx} & \varepsilon_{yx} & \varepsilon_{zx}\\ \varepsilon_{xy} & \varepsilon_{yy} & \varepsilon_{zy}\\ \varepsilon_{xz} & \varepsilon_{yz} & \varepsilon_{zz} \end{bmatrix}= \begin{bmatrix} \varepsilon_{xx} & \frac{1}{2}\gamma_{yx} & \frac{1}{2}\gamma_{zx}\\ \frac{1}{2}\gamma_{xy} & \varepsilon_{yy} & \frac{1}{2}\gamma_{zy}\\ \frac{1}{2}\gamma_{xz} & \frac{1}{2}\gamma_{yz} & \varepsilon_{zz} \end{bmatrix}\tag{29} εxxεxyεxzεyxεyyεyzεzxεzyεzz = εxx21γxy21γxz21γyxεyy21γyz21γzx21γzyεzz (29)
同样定义应变偏张量,有如下形式
[ e x x e y x e z x e x y e y y e z y e x z e y z e z z ] = [ ε x x ε y x ε z x ε x y ε y y ε z y ε x z ε y z ε z z ] − [ ε m 0 0 0 ε m 0 0 0 ε m ] (30) \begin{bmatrix} e_{xx} & e_{yx} & e_{zx}\\ e_{xy} & e_{yy} & e_{zy}\\ e_{xz} & e_{yz} & e_{zz} \end{bmatrix}= \begin{bmatrix} \varepsilon_{xx} & \varepsilon_{yx} & \varepsilon_{zx}\\ \varepsilon_{xy} & \varepsilon_{yy} & \varepsilon_{zy}\\ \varepsilon_{xz} & \varepsilon_{yz} & \varepsilon_{zz} \end{bmatrix}-\begin{bmatrix} \varepsilon_{m} & 0 & 0\\ 0 & \varepsilon_{m} & 0\\ 0 & 0 & \varepsilon_{m} \end{bmatrix}\tag{30} exxexyexzeyxeyyeyzezxezyezz = εxxεxyεxzεyxεyyεyzεzxεzyεzz εm000εm000εm (30)
其中 ε m = 1 3 ( ε x x + ε y y + ε z z ) \varepsilon_{m}=\frac{1}{3}(\varepsilon_{xx}+\varepsilon_{yy}+\varepsilon_{zz}) εm=31(εxx+εyy+εzz)

1.6 特殊应力、应变定义

定义应力强度或等效应力 σ ‾ \overline\sigma σ
σ ‾ = 3 J 2 = 3 2 s i j s i j = 1 2 [ ( σ 1 − σ 2 ) 2 + ( σ 1 − σ 3 ) 2 + ( σ 2 − σ 3 ) 2 ] = 1 2 [ ( σ x x − σ y y ) 2 + ( σ x x − σ z z ) 2 + ( σ y y − σ z z ) 2 + 6 ( τ x z 2 + τ x y 2 + τ y z 2 ) ] (31) \overline\sigma=\sqrt{3J_2}=\sqrt{\frac{3}{2}s_{ij}s_{ij}}\\ =\sqrt{\frac{1}{2}[(\sigma_{1}-\sigma_{2})^2+(\sigma_{1}-\sigma_{3})^2+(\sigma_{2}-\sigma_{3})^2]}\\ =\sqrt{\frac{1}{2}[(\sigma_{xx}-\sigma_{yy})^2+(\sigma_{xx}-\sigma_{zz})^2+(\sigma_{yy}-\sigma_{zz})^2+6(\tau_{xz}^2+\tau_{xy}^2+\tau_{yz}^2)]} \tag{31} σ=3J2 =23sijsij =21[(σ1σ2)2+(σ1σ3)2+(σ2σ3)2] =21[(σxxσyy)2+(σxxσzz)2+(σyyσzz)2+6(τxz2+τxy2+τyz2)] (31)
定义应变强度或等效应变 ε ‾ \overline \varepsilon ε
ε ‾ = 2 3 e i j e i j (32) \overline \varepsilon=\sqrt{\frac{2}{3}e_{ij}e_{ij}} \tag{32} ε=32eijeij (32)

定义剪切等效应力 T ‾ \overline T T
T ‾ = 1 2 s i j s i j (33) \overline T=\sqrt{\frac{1}{2}s_{ij}s_{ij}} \tag{33} T=21sijsij (33)
定义剪切等效应变 Γ ‾ \overline\Gamma Γ
Γ ‾ = 2 e i j e i j (34) \overline\Gamma=\sqrt{2e_{ij}e_{ij}} \tag{34} Γ=2eijeij (34)
加上上面定义的八面体剪应力、八面体剪应变
τ 8 = 1 3 s i j s i j γ 8 = 4 3 e i j e i j (35) \tau_8=\sqrt{\frac{1}{3}s_{ij}s_{ij}}\\ \gamma_8=\sqrt{\frac{4}{3}e_{ij}e_{ij}}\tag{35} τ8=31sijsij γ8=34eijeij (35)

至于为什么定义这些应力应变,我们在后面再介绍。

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/u012915522/article/details/135119056

智能推荐

874计算机科学基础综合,2018年四川大学874计算机科学专业基础综合之计算机操作系统考研仿真模拟五套题...-程序员宅基地

文章浏览阅读1.1k次。一、选择题1. 串行接口是指( )。A. 接口与系统总线之间串行传送,接口与I/0设备之间串行传送B. 接口与系统总线之间串行传送,接口与1/0设备之间并行传送C. 接口与系统总线之间并行传送,接口与I/0设备之间串行传送D. 接口与系统总线之间并行传送,接口与I/0设备之间并行传送【答案】C2. 最容易造成很多小碎片的可变分区分配算法是( )。A. 首次适应算法B. 最佳适应算法..._874 计算机科学专业基础综合题型

XShell连接失败:Could not connect to '192.168.191.128' (port 22): Connection failed._could not connect to '192.168.17.128' (port 22): c-程序员宅基地

文章浏览阅读9.7k次,点赞5次,收藏15次。连接xshell失败,报错如下图,怎么解决呢。1、通过ps -e|grep ssh命令判断是否安装ssh服务2、如果只有客户端安装了,服务器没有安装,则需要安装ssh服务器,命令:apt-get install openssh-server3、安装成功之后,启动ssh服务,命令:/etc/init.d/ssh start4、通过ps -e|grep ssh命令再次判断是否正确启动..._could not connect to '192.168.17.128' (port 22): connection failed.

杰理之KeyPage【篇】_杰理 空白芯片 烧入key文件-程序员宅基地

文章浏览阅读209次。00000000_杰理 空白芯片 烧入key文件

一文读懂ChatGPT,满足你对chatGPT的好奇心_引发对chatgpt兴趣的表述-程序员宅基地

文章浏览阅读475次。2023年初,“ChatGPT”一词在社交媒体上引起了热议,人们纷纷探讨它的本质和对社会的影响。就连央视新闻也对此进行了报道。作为新传专业的前沿人士,我们当然不能忽视这一热点。本文将全面解析ChatGPT,打开“技术黑箱”,探讨它对新闻与传播领域的影响。_引发对chatgpt兴趣的表述

中文字符频率统计python_用Python数据分析方法进行汉字声调频率统计分析-程序员宅基地

文章浏览阅读259次。用Python数据分析方法进行汉字声调频率统计分析木合塔尔·沙地克;布合力齐姑丽·瓦斯力【期刊名称】《电脑知识与技术》【年(卷),期】2017(013)035【摘要】该文首先用Python程序,自动获取基本汉字字符集中的所有汉字,然后用汉字拼音转换工具pypinyin把所有汉字转换成拼音,最后根据所有汉字的拼音声调,统计并可视化拼音声调的占比.【总页数】2页(13-14)【关键词】数据分析;数据可..._汉字声调频率统计

linux输出信息调试信息重定向-程序员宅基地

文章浏览阅读64次。最近在做一个android系统移植的项目,所使用的开发板com1是调试串口,就是说会有uboot和kernel的调试信息打印在com1上(ttySAC0)。因为后期要使用ttySAC0作为上层应用通信串口,所以要把所有的调试信息都给去掉。参考网上的几篇文章,自己做了如下修改,终于把调试信息重定向到ttySAC1上了,在这做下记录。参考文章有:http://blog.csdn.net/longt..._嵌入式rootfs 输出重定向到/dev/console

随便推点

uniapp 引入iconfont图标库彩色symbol教程_uniapp symbol图标-程序员宅基地

文章浏览阅读1.2k次,点赞4次,收藏12次。1,先去iconfont登录,然后选择图标加入购物车 2,点击又上角车车添加进入项目我的项目中就会出现选择的图标 3,点击下载至本地,然后解压文件夹,然后切换到uniapp打开终端运行注:要保证自己电脑有安装node(没有安装node可以去官网下载Node.js 中文网)npm i -g iconfont-tools(mac用户失败的话在前面加个sudo,password就是自己的开机密码吧)4,终端切换到上面解压的文件夹里面,运行iconfont-tools 这些可以默认也可以自己命名(我是自己命名的_uniapp symbol图标

C、C++ 对于char*和char[]的理解_c++ char*-程序员宅基地

文章浏览阅读1.2w次,点赞25次,收藏192次。char*和char[]都是指针,指向第一个字符所在的地址,但char*是常量的指针,char[]是指针的常量_c++ char*

Sublime Text2 使用教程-程序员宅基地

文章浏览阅读930次。代码编辑器或者文本编辑器,对于程序员来说,就像剑与战士一样,谁都想拥有一把可以随心驾驭且锋利无比的宝剑,而每一位程序员,同样会去追求最适合自己的强大、灵活的编辑器,相信你和我一样,都不会例外。我用过的编辑器不少,真不少~ 但却没有哪款让我特别心仪的,直到我遇到了 Sublime Text 2 !如果说“神器”是我能给予一款软件最高的评价,那么我很乐意为它封上这么一个称号。它小巧绿色且速度非

对10个整数进行按照从小到大的顺序排序用选择法和冒泡排序_对十个数进行大小排序java-程序员宅基地

文章浏览阅读4.1k次。一、选择法这是每一个数出来跟后面所有的进行比较。2.冒泡排序法,是两个相邻的进行对比。_对十个数进行大小排序java

物联网开发笔记——使用网络调试助手连接阿里云物联网平台(基于MQTT协议)_网络调试助手连接阿里云连不上-程序员宅基地

文章浏览阅读2.9k次。物联网开发笔记——使用网络调试助手连接阿里云物联网平台(基于MQTT协议)其实作者本意是使用4G模块来实现与阿里云物联网平台的连接过程,但是由于自己用的4G模块自身的限制,使得阿里云连接总是无法建立,已经联系客服返厂检修了,于是我在此使用网络调试助手来演示如何与阿里云物联网平台建立连接。一.准备工作1.MQTT协议说明文档(3.1.1版本)2.网络调试助手(可使用域名与服务器建立连接)PS:与阿里云建立连解释,最好使用域名来完成连接过程,而不是使用IP号。这里我跟阿里云的售后工程师咨询过,表示对应_网络调试助手连接阿里云连不上

<<<零基础C++速成>>>_无c语言基础c++期末速成-程序员宅基地

文章浏览阅读544次,点赞5次,收藏6次。运算符与表达式任何高级程序设计语言中,表达式都是最基本的组成部分,可以说C++中的大部分语句都是由表达式构成的。_无c语言基础c++期末速成