【小呆的力学笔记】弹塑性力学的初步认知二:应力应变分析(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

智能推荐

文件夹权限问题和linux下搭建FTP服务器_linux 文件夹权限 影响ftp listfles-程序员宅基地

文章浏览阅读2.6w次。文件夹权限问题Linux、Fedora、Ubuntu修改文件、文件夹权限的方法差不多。很多人开始接触Linux时都很头痛Linux的文件权限问题。这里告诉大家如何修改Linux文件-文件夹权限。以主文件夹下的一个名为cc的文件夹为例。 下面一步一步介绍如何修改权限: 1.打开终端。输入su(没 Linux、Fedora、Ubuntu修改文件、文件夹权限的方法差不多。很多人开始接触Lin_linux 文件夹权限 影响ftp listfles

JavaScript调用Python程序_javascript 调用python-程序员宅基地

文章浏览阅读186次。JavaScript调用Python程序_javascript 调用python

这个为生信学习打造的开源Bash教程真香!!(目录更新)!-程序员宅基地

文章浏览阅读300次。生物信息学习的正确姿势NGS系列文章包括NGS基础、在线绘图、转录组分析(Nature重磅综述|关于RNA-seq你想知道的全在这)、ChIP-seq分析(ChIP-seq基本分析流..._bash教程 国外

oracle详解-程序员宅基地

文章浏览阅读4.3k次,点赞11次,收藏65次。首先看张图:对于一个数据库系统来说,假设这个系统没有运行,我们所能看到的和这个数据库相关的无非就是几个基于操作系统的物理文件,这是从静态的角度来看,如果从动态的角度来看呢,也就是说这个数据库系统运行起来了,能够对外提供服务了,那就意外着数据库系统启动了自己的一个实例,综合以上2个角度,Oracle如何定义上述描述呢?我们来引入第一个概念,Oracle服务器,所谓Oracle服务器是一个数据库管理系统,它包括一个Oracle实例(动态)和一个Oracle数据库(静态)。Oracle实例是一个运行的概念(_oracle

Cadence学习记录_layout如何提取参数-程序员宅基地

文章浏览阅读1.3k次,点赞2次,收藏14次。基于模拟集成电路仿真软件cadence零基础入门的一点学习记录_layout如何提取参数

Obsidian多端同步化方案_obsidian同步方案-程序员宅基地

文章浏览阅读1.7k次。Notion与Obsidian是目前市场上功能最强大,使用用户最多,门槛也最高的两款笔记软件。本文不会对Notion和Obsidian两款笔记软件的优劣进行对比,但从两者的存储方式,我们可以看到Notion和Obsidian是两个极端——前者只支持云存储,后者只支持本地存储。_obsidian同步方案

随便推点

java爬虫之登录到教务系统抓取成绩_java实现从学校教务网上爬取数据(-程序员宅基地

文章浏览阅读2.2k次。最近使用java写了个爬虫,可能我对java比较熟悉,所以相对于python来说,我觉得用java写更得心应手些。我采用的是java的jsoup,以及解析用到的json先放上学校教务系统的url http://222.200.98.147首先可以看到,这里是需要验证码输入的,所以我使用了以下的思路:第一步,先访问验证码所在的url,把图片下载到本地,然后保存cookie_java实现从学校教务网上爬取数据(

[ERROR] Can't find error-message file '/data/mysql/share/errmsg.sys'. Check error-message file locat...-程序员宅基地

文章浏览阅读5.6k次。1. MySQL5.7.21启动时报错:[ERROR] Can't find error-message file '/data/mysql/3307/share/errmsg.sys'. Check error-message file location and 'lc-messages-dir' configuration directive.2. 登录MySQL查看系统全局参数:..._can't find error-message file '/data/mysql/share/errmsg.sys'. check error-me

51单片机入门——数字时钟_51单片机 时钟-程序员宅基地

文章浏览阅读2.5w次,点赞58次,收藏419次。文章目录1. 前言1.1. 设计要求2. 硬件原理2.1. 时钟信号(晶振)2.2. 按键开关2.3. 数码管显示3. 原理图3.1. 仿真原理图3.2. AD原理图3.3. PCB图4. 软件设计4.1. 初版代码(无年月日)4.2. 终版代码5. 元器件清单5.1. 仿真软件5.2. 实物1. 前言在此之前我们已经学习了单片机的定时器、中断、数码管。这篇文章主要讲述如何用上述的知识自己制作一个基于51单片机的数字时钟。1.1. 设计要求(1)主电路由秒信号发生器、“时、分、秒”计数器_51单片机 时钟

【华为OD机考 统一考试机试】区间交集(C++ Java JavaScript Python)_给定一组闭区间,其中部分区间存在交集,任意两个给定区间的交集,称为公共区间-程序员宅基地

文章浏览阅读3.4k次。任意两个给定区间的交集,称为公共区间(如:[1,2],[2,3]的公共区间为[2,2],[3,5],[3,6]的公共区间为[3,5])。公共区间之间若存在交集,则需要合并(如:[1,3],[3,5]区间存在交集[3,3],需合并为[1,5])。区间元素为 X: -10000_给定一组闭区间,其中部分区间存在交集,任意两个给定区间的交集,称为公共区间

yum下载Zabbix4.0失败的解决方法_warning: failed loading '/etc/yum.repos.d/zabbix.r-程序员宅基地

文章浏览阅读1.7k次。根据官网说明配置的yum源,今天用yum下载Zabbix时莫名的报错,经过几番折腾,找到了解决方法。一、报错如下:二、 解决方法:[root@VM_0_6_centos ~]# cat /etc/yum.repos.d/zabbix.repo # 这个yum源是官方地址,服务器应该是在国外的,下载的时候报错[zabbix] ..._warning: failed loading '/etc/yum.repos.d/zabbix.repo', skipping.

人机交互实验五_人机交互界面设计演练-程序员宅基地

文章浏览阅读6.4k次,点赞2次,收藏10次。5.9.1实验目的本章实验的目的是:1)了解和熟悉人机界面设计过程管理的相关知识;2)了解和评价游戏软件的人机交互设计,提高自己的评价能力,提高自己对设计水平的鉴赏能力。5.9.2工具/准备工作在开始本实验之前,请回顾课文的相关内容。需要准备一台能够访问因特网的计算机。5.9.3实验内容与步骤1.概念理解1)成功的用户界面开发有4个支柱,它们能够帮助用户界面架构师将好的思想转化为成功的系统。经验表明,每个支柱都能在此过程中产生..._人机交互界面设计演练