有限元法2011-应力张量
基于有限元法的复合材料结构的应力分析

基于有限元法的复合材料结构的应力分析复合材料是一种由两种或两种以上不同性质的材料组成的复合材料,具有优异的力学性能和轻质化的特点。
在现代工程设计中,复合材料广泛应用于航空航天、汽车、建筑等领域。
为了确保复合材料结构的可靠性和安全性,需要进行应力分析。
本文将介绍基于有限元法的复合材料结构的应力分析方法。
有限元法是一种数值分析方法,通过将结构离散为有限数量的单元,将结构的连续性问题转化为离散的代数方程组,从而求解结构的应力、位移等物理量。
在复合材料结构的应力分析中,有限元法是一种非常有效的工具。
首先,需要建立复合材料结构的有限元模型。
复合材料结构通常由多层纤维增强复合材料层和粘接剂层组成。
在有限元模型中,可以将复合材料层和粘接剂层分别建模为不同的单元类型。
对于复合材料层,可以采用壳单元或板单元进行建模,而对于粘接剂层,可以采用固体单元进行建模。
此外,还需要定义材料的力学性质,如弹性模量、剪切模量等。
其次,需要施加边界条件和加载条件。
边界条件是指结构的约束条件,用于限制结构的自由度。
加载条件是指施加在结构上的外部载荷。
在复合材料结构的应力分析中,边界条件和加载条件的选择对结果的准确性有重要影响。
需要根据实际情况选择合适的边界条件和加载条件。
然后,进行有限元分析。
有限元分析是通过求解离散的代数方程组来计算结构的应力和位移。
在复合材料结构的有限元分析中,通常采用迭代求解的方法,通过迭代计算得到结构的应力和位移。
在求解过程中,需要选择合适的求解方法和收敛准则,以确保结果的准确性和稳定性。
最后,对分析结果进行后处理。
分析结果包括结构的应力、位移等物理量。
可以通过可视化的方式将结果呈现出来,如应力云图、位移云图等。
此外,还可以通过对结果进行进一步的分析和评估,如应力分布的均匀性、结构的强度等。
综上所述,基于有限元法的复合材料结构的应力分析是一种有效的工具。
通过建立有限元模型、施加边界条件和加载条件、进行有限元分析以及对分析结果进行后处理,可以得到复合材料结构的应力分布情况,为工程设计提供有力支持。
构造应力场模拟——有限元理论、方法和研究进展

构造应力场模拟——有限元理论、方法和研究进展张胜利【摘要】采用有限元数值模拟方法对构造地质问题进行描述和定量化求解是当前地质学领域的研究的一个热点,在近10年以来取得了重要进展,形成了比较完整的理论和技术体系,并在一些典型的地质构造带获得了重要的研究成果.本文以有限元数值模拟方法理论作为出发点,总结分析了国内外有限元数值模拟方法在构造应力场领域的研究进展情况和技术方法,并讨论了其目前存在的问题和未来发展方向.【期刊名称】《地震工程学报》【年(卷),期】2010(032)004【总页数】6页(P405-410)【关键词】构造应力场;数值模拟;有限单元法【作者】张胜利【作者单位】五邑大学信息学院,广东,江门,529020;中科院广州地化所,广东,广州,510640【正文语种】中文【中图分类】P315.12Abstract:The Finite Element Method(FEM)has been used in the study of tectonic stress field for a long time,and the essence of numerical modeling has been adopted to the well-established numerical methods of multidisciplinary acknowledge including mathematics,physics andmechanics for studing characters of geological tectonics.In the last decade,great advances have been made on the numerical simulation method,not only an integrated theory has been built up,but also some significant results have been born from several typical tectonic belts.So the FEM becomes one of the most important numerical methods in the study of tectonic stress field.In this paper,taking theory of FEM as a springboard,the new progress and methods in this field at home and abroad is summarized and analyzed.Some problems and prospect of the researching on the field is also given.Key words:Tectonic stress field;Numerical model;Finite element method地壳中的各种地质构造都是岩石受力发生变形的产物,它们的产生和发展必然也受力学规律的支配。
应力张量的定义

应力张量的定义什么是应力张量应力是物体内部的力对单位面积的分布情况。
在固体力学中,应力张量是一个描述物质内部应力状态的重要工具。
应力张量是一个二阶的张量,用于描述各个方向的应力分量。
构成应力张量的分量应力张量由9个分量组成,它们分别是X方向上的正应力、Y方向上的正应力、Z 方向上的正应力以及XY平面上的剪应力、YZ平面上的剪应力和ZX平面上的剪应力。
这些分量可以用一个3x3的矩阵表示,矩阵的每个元素代表相应的应力分量。
应力张量的矩阵表示如下:[ σ11 σ12 σ13 ][ σ21 σ22 σ23 ][ σ31 σ32 σ33 ]其中,σ11、σ22和σ33分别代表X、Y和Z方向上的正应力,σ12、σ21代表XY平面上的剪应力,σ13、σ31代表XZ平面上的剪应力,σ23、σ32代表YZ平面上的剪应力。
应力张量的物理意义应力张量的物理意义是它能够给出物体内部应力状态的详细信息。
通过应力张量,我们可以了解到不同方向上的应力分量的大小和方向。
应力张量可以用来计算物体的应力分布、应力变化以及应力的大小。
在材料工程、土木工程和机械工程等领域中,应力张量是非常重要的。
应力张量的计算方法计算应力张量的方法可以通过实验测定和数值模拟两种方式进行。
实验测定方法实验测定方法需要使用一些仪器设备来测量物体内部的应力情况。
常用的实验测定方法有:1.应变计法:通过在物体表面安装应变计,测量应变分布从而计算得到应力分布。
2.简支梁法:通过在梁上进行弯曲试验,测量不同位置上的弯矩和曲率来计算得到应力分布。
3.角度测量法:通过在物体表面安装角度测量仪,测量不同位置上的角度变化来计算得到应力分布。
数值模拟方法数值模拟方法是利用计算机进行模拟计算,通过输入物体的几何形状、材料力学性质和加载条件等参数,通过数学计算得到应力分布。
常用的数值模拟方法有:1.有限元法:将物体分割为有限个小单元,利用力学原理和数学方法建立方程组,通过求解方程组得到应力分布。
弹性力学与有限元智慧树知到答案章节测试2023年武汉工程大学

第一章测试1.下列不属于弹性力学研究对象的是()。
A:板壳B:刚体C:杆件D:实体结构答案:B2.下列不属于弹性力学中基本未知量的是()。
A:位移分量B:应力分量C:面力分量D:应变分量答案:C3.在工程强度校核中起着重要作用的是()。
A:应力分量B:主应力C:正应力D:切应力答案:B4.已知物体内某点的应力张量(单位:Pa),则沿方向的正应力大小为()。
A:222.22 PaB:888.89 PaC:666.67 PaD:444.44 Pa答案:D5.下列关于应力分量的说法,正确的有()。
A:坐标面上的应力B:一点的9个应力分量可以完全确定该点的应力状态C:应力分量与面力分量的正负号规定相同D:正截面上的应力E:弹性力学中应力分量的正负号规定反映了作用力与反作用力原理以及“受拉为正、受压为负”的传统观念。
答案:ABDE6.理想弹性体满足的假设有()。
A:无初始应力假设B:均匀性假设C:连续性假设D:完全弹性假设E:各向同性假设答案:BCDE7.建立在基本假设上的弹性力学,也称为()。
A:弹性理论B:线性弹性力学C:应用弹性力学D:数学弹性力学答案:ABD8.弹性力学的主要任务是解决各类工程中所提出的问题,这些问题包括()。
A:稳定B:刚度C:强度D:动力答案:ABC9.弹性力学的研究方法是在弹性体的区域内严格考虑三方面条件,建立三套基本方程,这三方面条件包括()。
A:几何学B:物理学C:静力学D:动力学答案:ABC10.中国科学家胡海昌于1954年最早提出了三类变量的广义变分原理。
()A:错B:对答案:B11.物体内任意一点的应力分量、应变分量和位移分量,都不随该点的位置而变化,它们与位置坐标无关。
()A:对B:错答案:B12.在最大正应力的作用面上切应力为零,在最大切应力的作用面上正应力为零。
()A:对B:错答案:B13.应力张量的三个不变量是与坐标选择无关的标量。
()A:错B:对答案:B14.弹性力学与材料力学在研究方法上是完全相同的。
同济大学地下结构有限元11-应力磨平与处理

庄晓莹
弹性力学中的有限元
9
一些符号约定
精确解
近似解 近似解 (改善后) (改善前)
庄晓莹
弹性力学中的有限元
10
FEM应力计算方法
庄晓莹
弹性力学中的有限元
11
FEM应力计算方法
整体刚度矩阵求解位移 复习!得到的是什么变量?
对于每个单元则有
这里的上标 i 表示单元编号
我们首先以三结点常应变单元为例,单元结点应力计算方法为
原问题
等价问题
庄晓莹
弹性力学中的有限元
27
FEM应力结果改善的主要方法
庄晓莹
弹性力学中的有限元
28
FEM应力改善的主要方法
应力磨平的方法多种多样:
整体磨平 整体意义上对应力值的“近似平均”处理
局部磨平 单元内部对应力值的“近似平均”处理 分片磨平 绕结点周围的单元应力值的“近似平均”处理
地震作用下钢结构损伤过程数值模拟_段红霞

性硬化参数,还是随动硬化参数,难以真正通过试
验来测定,这时可以利用这些材料已有的常规试验
数据,比如对称应变循环试验数据、半循环应力-
应变数据、稳定循环应力-应变数据等,对这些数据
进行处理或通过计算模拟来获得所需用值。
在弹塑性损伤模型实际应用时,首先要判断损
伤的出现。塑性变形发展到一定阶段,微孔洞、微
裂纹不断累积,当满足下列准则时,认为结构出现
延性损伤[16]:
∫ ωD =
dε pl
ε
pl D
(η
,
ε
pl
)
=
1
(9)
式中: ωD 为损伤状态变量,随着塑性变形单调递
增;
ε
pl D
为延性损伤出现时的相当塑性应变,是三
轴应力度 η 和相当塑性应变率的函数,η = − p / q ,
摘 要:采用弹塑性损伤本构模型,该模型结合了非线性各向同性和随动强化准则以及基于塑性位移的损伤演化 规律,利用 ABAQUS 对一个 9 层钢结构在 EL-Centro 地震波作用下塑性变形和损伤区域的发展过程进行了数值 模拟。结果表明,上部楼层的梁端产生较明显的塑性变形并形成损伤部位。这与在 Northridge 地震中观测到的高 层钢结构的地震损伤情况是一致的。 关键词:地震荷载;钢结构;损伤演化;数值模拟;ABAQUS 中图分类号:TU393.2 文献标识码:A
李守巨(1960―),男,辽宁人,教授,博士,从事工程力学和参数反演研究(E-mail: lishouju@); 刘迎曦(1944―),男,四川人,教授,博士,从事有限元研究(E-mail: yxliu@).
工程力学
199
材的本构关系),对结构的动力方程采用逐步积分的 方法计算地震过程中每一瞬时结构的位移、速度、 加速度、内力反应,从而可以分析结构在弹性和非 弹性阶段的内力变化及构件破坏、直至倒塌的全过 程。进行结构弹塑性地震反应分析的关键问题是建 立循环荷载下精确的材料本构模型和计算模型。对 于建筑用钢材,常用的本构模型有理想弹塑性模 型、各向同性强化模型及随动强化模型。理想弹塑 性模型无法描述材料的硬化特性,与材料的实验结 果有一定偏差。同性强化理论允许屈服面膨胀、收 缩,适用于单调加载,对于循环塑性,同性强化不 能反映包辛格(Bauschinger)效应。随动强化理论假 定屈服面在应力空间中可以平移,但不能转动、膨 胀和收缩。随动强化理论比同性强化理论前进了一 步,但随动强化只适应于小应变的情况。对复杂荷 载历史工况,同性强化、随动强化都不能真实描述 钢材的循环特性。由 Hodge 提出,并由 Axelsson 和 Samuelsson[1]进一步发展的混合强化模型将各向 同性强化和随动强化结合起来,屈服面既能膨胀(收 缩)又能平移,可考虑钢材的包辛格(Bauschinger)效 应和屈服平台,模型简单而实用。许红胜介绍了一 种新的双曲面混合模型,对非比例加载有良好的适 用性,可用于钢结构在复杂动力荷载作用下的塑性 分析[2]。但是上述这些模型没有涉及到钢材损伤的 影响,实际上原生材料就存在损伤[3],随着荷载的 循环作用,在材料中的微裂纹、微孔洞、剪切带等 细观损伤萌生、串接、汇合、扩展,从而形成损伤 的动态演化过程,直接影响着钢结构的性能。因此, 要客观地描述建筑用钢材的循环本构关系,不可忽 略损伤对材性的不可逆劣化影响[4―7]。郑宏提出了 结构钢弹塑性各向异性损伤本构模型,采用非线性 有限元方法,分析了箱形柱在单轴循环荷载作用下 的滞回性能[8]。丁阳考虑了钢材的损伤累积效应和 应变强化效应,应用塑性应变和能量耗散理论建立 了钢材的损伤力学模型[9]。王连坤基于钢材各向同 性塑性累积损伤本构关系,推导了考虑材料损伤和 混合强化本构关系的弹塑性刚度矩阵,算例证明文 中方法可以达到较高精度[10]。Mashayekhi 采用连续 损伤本构模型,计算了带槽口矩形截面试件的损伤 参数,并通过试验证明了该模型的有效性[11]。为了 能够真实描述钢结构经受循环荷载时的损伤破坏 特性,本文同时考虑了钢材的混合强化模型和损伤 演化规律,利用高效有限元 ABAQUS 模拟了一个 9
超弹性仿真计算公式

超弹性仿真计算公式超弹性材料是一类具有非线性、大变形和大应变能力的材料,常见于橡胶、聚合物等材料中。
超弹性材料的力学行为与普通材料存在很大的差异,因此需要特殊的计算方法来描述其力学性能。
超弹性仿真计算公式是描述超弹性材料力学行为的重要工具,本文将介绍超弹性仿真计算公式的基本原理和应用。
超弹性材料的力学行为可以用应力-应变关系来描述,而超弹性材料的应力-应变关系通常不遵循胡克定律,因此需要使用特殊的公式来描述。
在超弹性仿真计算中,常用的描述超弹性材料力学行为的公式包括Mooney-Rivlin模型、Ogden模型、Yeoh模型等。
这些模型都是基于应变能密度函数来描述超弹性材料的应力-应变关系,其基本形式如下:Mooney-Rivlin模型:$$W = C_1(I_1-3) + C_2(I_2-3)$$。
Ogden模型:$$W = \sum_{i=1}^{N} \frac{\mu_i}{\alpha_i}(\lambda_1^{\alpha_i} +\lambda_2^{\alpha_i} + \lambda_3^{\alpha_i} 3)$$。
Yeoh模型:$$W = \sum_{i=0}^{N} C_i(I_1-3)^i$$。
其中,W表示应变能密度函数,C1、C2、μi、αi、Ci等参数是需要通过实验或拟合得到的材料参数,I1、I2、I3分别是应变张量的主应变,λ1、λ2、λ3分别是应变张量的主应变比。
这些公式通过描述应变能密度函数来描述超弹性材料的应力-应变关系,可以很好地描述超弹性材料的非线性和大变形行为。
超弹性仿真计算公式的应用可以帮助工程师和科研人员更好地理解和预测超弹性材料的力学行为。
在工程设计中,超弹性仿真计算可以用来预测超弹性材料在复杂加载条件下的应力分布和变形情况,从而指导材料选择和结构设计。
在科学研究中,超弹性仿真计算可以用来研究超弹性材料的力学行为和性能,为新材料的设计和开发提供重要参考。
应力张量的认识(二)

应力张量的认识(二)本文主要是对材料成形相关专业学习过程中对一些问题的思考,也许并不深刻,但却是自己从初学时的迷惑到后来逐渐认识的过程。
相关还有:Levy-Mises理论的思考应力张量的基础知识参见应力张量的认识(一)这一部分主要是对应力张量本质的理解。
相似矩阵通过基础知识我们已经认识到,应力张量代表点的应力状态,它不依赖于坐标系的选取,并且对应着同一主应力状态。
用矩阵的观点理解为:同一点的应力张量矩阵是相似矩阵,并且可以对角化。
然而问题是为什么会这样?我们可以这么理解,同一点不同截面下的应力张量描述的都是相同的应力状态,因此他们有着内在的联系(例如满足静力平衡方程)。
由切应力互等定律可知应力张量矩阵是实对称矩阵,由矩阵论可知实对称矩阵必定可以对角化,即不同截面应力张量对相似于一个主应力张量;而同一点的主应力状态是确定的。
于是由相似矩阵的传递性可知,不同截面下的应力张量矩阵是相似的。
线性变换从应力张量是相似矩阵再进一步——相似矩阵的本质是同一线性变换在不同基下表示的矩阵,我们就可以从更根本的角度看待应力张量了:应力张量代表一个线性变换!这是一个抽象的认识,或者说是从相似矩阵推论出来的。
那么接下来让他具体化:从前面已经知道,利用三个相互垂直的截面截取一点P,将各个截面上的所有应力分量组成一个整体的物理量——应力张量。
以这三个截面的法线方向为正方向建立笛卡尔坐标系,如下图所示。
那么对于任意法向为n=(n1,n2,n3)截面,可以得到面上作用的应力分量将上式改写成矩阵表达式其中,根据线性变换与矩阵的一一对应的关系,可知应力张量代表一个线性变换(确切的说是线性映射)——应力张量将截面位置映射到截面应力。
数学表述为:U是截面方向余弦组成的线性空间,V是截面应力作成的线性空间。
举例说明下面说明这个线性映射是如何与力学描述相一致的。
前面多次提到,用三个互相垂直的截面截取P点,并以其法向建立笛卡尔坐标系。
这样的截面有无数多对,这样的坐标系也有无数多个。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
σ 33
σ 31
3.4
主应力
如果作用在某一微分面上的应力矢量和这一微分面垂直, 如果作用在某一微分面上的应力矢量和这一微分面垂直,即这 一微分面上只有正应力而无剪应力,则这一微分面称为主平面 一微分面上只有正应力而无剪应力,则这一微分面称为主平面 其法线方向称为应力主方向 其上的应力称为主应力 应力主方向, 主应力。 ,其法线方向称为应力主方向,其上的应力称为主应力。如果 三个坐标轴方向都是主方向,则称这一坐标系为主坐标系 主坐标系。 三个坐标轴方向都是主方向,则称这一坐标系为主坐标系。 在主平面上,正应力取极值、剪应力为零。 在主平面上,正应力取极值、剪应力为零 n主平面上的单位法向量,σ主应力,主平面上的应力矢量 主平面上的单位法向量, 主应力 主应力, 主平面上的单位法向量
T1 , T2 , T3
T = T1i + T2 j + T3k
分别为应力矢量T(n)沿三个坐标方向的分量 沿三个坐标方向的分量 分别为应力矢量
3.2 应力张量
通过一点P 通过一点 的各个面上应力状况的集合 —— 称为一点的应力状态 x面的应力: σx ,τ xy ,τ xz 面的应力: 面的应力 y面的应力: 面的应力: 面的应力 z面的应力: 面的应力: 面的应力
三个微分面的外法向单位矢量e 三个微分面的外法向单位矢量 1, e2,e3.对应的应力矢量为 对应的应力矢量为
T, T , T 1 2 3
e3
σ 11
σ 13 σ 12
σ 23
σ 33
e3
σ 21
σ 22
e3
σ 32 σ 31
e2 e1
e2
e1
e2
e1
e1面的应力: 11,σ12,σ13 面的应力: σ
z
∆F
T 3
T 1
2 ∆S T
(1) T 是坐标的连续分布函数; 说明: 说明: (2) T 的加载方式是任意的;
k
i
O j
y
x (3) T T T的正负号由坐标方向确定。 的正负号由坐标方向确定。 3 1 2
3 、应力 (1) 一点应力的概念
(1) 物体内部分子或原子间的相互 内力 不考虑) 不考虑 作用力; (不考虑 作用力 (2) 由于外力作用引起的相互作用力 由于外力作用引起的相互作用力.
力分量的方向。 力分量的方向。 •以上 个分量,构成应力张量在笛卡儿坐标系下的分量 以上9个分量 构成应力张量 应力张量在笛卡儿坐标系下的分量 以上 个分量,
σ x τ xy τ xz σ ij = τ yx σ y τ yz τ τ σ z zx zy
取代下标x、 、 , 张量表示 用1、2、3取代下标 、y、z, 、 、 取代下标
即
σ ij = σ ji
(13)
应力张量是对称的,九个应力分量中只有六个是独立的。 应力张量是对称的,九个应力分量中只有六个是独立的。这一结 论称为剪应力互等定理 剪应力互等定理。 论称为剪应力互等定理。 剪应力互等定理: 剪应力互等定理:过物体内任意一点的两个相互垂直的微分面 和这两个微分面的交线垂直的两个剪应力相等。 上,和这两个微分面的交线垂直的两个剪应力相等。
σ e2面的应力: 22,σ21,σ23 面的应力:
σ e3面的应力: 33,σ31,σ32 面的应力:
T =σ11e1 +σ12e2 +σ13e3 1 T =σ21e1 +σ22e2 +σ23e3 2 T =σ31e1 +σ32e2 +σ33e3 3
Ti = σ ij e j
(8) )
σ ij 的第一个下标表示应力分量的作用面,第二个下标表示应 的第一个下标表示应力分量的作用面,
第3讲 应力张量 讲
3.1
外力与应力矢量
面力—— 作用于物体表面单位面积上的外力 面力
∆F —— 面力分布集度(矢量) 面力分布集度(矢量) T = lim S ∆S→ ∆ 0 (2) ) T =T i +T2j +T3k 1 T T T —— 面力矢量在坐标轴上投影 1 2 3
单位: 单位: 1N/m2 =1Pa (帕) 1MN/m2 = 106Pa = 1MPa (兆帕)
主应力的确定: 主应力的确定:
设主平面存在,其外法线为 , 设主平面存在,其外法线为n, 方向余弦: 方向余弦:n1,n2,n3
C
n
σn =σ 其上应力: 则:其上应力: τn = 0
T1 = n 1σ T 2 = n 2σ T 3 = n 3σ
代入
τ yx τ xy
z
σ
xz ττzx
σx
B
轴上的投影为: 在x,y,z轴上的投影为: , , 轴上的投影为
(b)
(a )
T = Ti ei
T = σ nn + τ ns
(5)
T
τn
P
σn
n
式中n 式中n和s分别为微分面的单位法向量和 单位切向量 σ n = T ⋅ n = Ti ni (6)
2 τn = T ⋅s = T 2 −σ n
(7) T为应力矢量 为应力矢量T(n)的大小,称为总应力 的大小, 为应力矢量 的大小 称为总应力
应力张量也满足张量的变换规律
弹性体的应力边界条件: 弹性体的应力边界条件: 当面abc为物体的边界面时,则其应力分量 为物体的边界面时, 当面
Tx ,Ty ,Tz
成为面力分量
Tx ,Ty ,Tz
e3
c
P
T x = n1σ x + n 2τ yx + n 3τ zx
由
n
b e2
T y = n1τ xy + n 2σ y + n 3τ zy T z = n1τ xz + n 2τ yz + n 3σ z
i i k kj
∫ r × TdS = ∫ r × (n ⋅σ )dS = ∫ x e × n σ
S S S
V
e j dS = ∫ ( xi ei × σ kj e j ),k dV
V
V
= ∫ (e k × σ kj e j= nkei k×⋅σ kjmje jm dV e j ∫=σ kjkekjimjiδ km e j∇ ⋅ σ )dV n ⋅ σ + xi e σ ,k e ) ⊗ = ( n σ e + r ×
n
σn
(法线 法线) 法线
应力分量
应力的切向分量
σn —— 正应力
τn —— 剪应力
单位: 单位
与面力相同
MPa (兆帕)
应力关于坐标连续分布的
σ =σ (x, y, z) τ =τ (x, y, z)
∆F
T ( r , − n) = −T ( r , n)
P
∆S
(4) )
P
∆S
n
−n
−∆F
应力矢量T(n)的下标 表示微分面 的下标n表示微分面 应力矢量 的下标 的外法线方向, 的外法线方向,它用于反映应力作 用面的方向。 用面的方向。
σ ij σ 11 = σ 21 σ 31 σ 12 σ 22 σ 32 σ 13 σ 23 σ 33
(9) )
过一点任意微分面上的应力矢量可由三个相互垂直的 微分面上的应力矢量表示出来, 微分面上的应力矢量表示出来,也即可由应力张量表 示出来。故应力张量完全确定了一点的应力状态。 示出来。故应力张量完全确定了一点的应力状态。
∆F —— 体力分布集度 f = lim 矢量) (矢量) ∆V→ ∆ 0 V f = f1i + f2j + f3k (1) )
f1、f2、f3为体力矢量在坐标轴上的投影 单位: 单位: N/m3 kN/m3
z
f3
∆F
k i
O j
f1 ∆V f2
y
(1) f 是坐标的连续分布函数; x 说明: 说明:(2) f 的加载方式是任意的 (如:重力,磁场力、惯性力等 如 重力,磁场力、惯性力等) (3) f1、f2、f3的正负号由坐标方向确定。 的正负号由坐标方向确定。
&& ∫ ( f − ρ u)dV + ∫ TdS = 0
V S
T = Ti ni = n ⋅ σ
(10) )
∫ TdS = ∫ n ⋅ σdS = ∫ ∇ ⋅ σdV
S S V
&& ∫ (∇ ⋅ σ + f − ρ u)dV = 0
V
V的任意性 的任意性
&& ∇ ⋅σ + f = ρu
运动方程
(11a) )
n⋅σ = σ n 即 σ ⋅n = σ n
T =σn
(14)
主应力σ是应力张量σ的特征值, 是 主应力 是应力张量σ的特征值, n是σ的特征矢量 是应力张量 谱定理可知,必有三个相互垂直的应力主方向, 谱定理可知,必有三个相互垂直的应力主方向,对应的 有三个主应力 σ 1 ≥ σ 2 ≥ σ 3 因此必定存在一个主坐标系。 因此必定存在一个主坐标系。
1 1 x
)s
s
τ + n2 (
2
yx y
s
3
zx
s
x
a e1
xy
s
3
zy
s
y
xz
s
2
yz
s
3
n
s
z
x s
y s
z s
xy s
xz s
yz s
3.3
平衡方程和运动方程
应力的变化并不是任意的, 应力的变化并不是任意的,应力张量的变化必须满足平衡条件或 动量定理和动量矩定理。 对任一块体积V,表面为 作用在V上的 表面为S,作用在 动量定理和动量矩定理。 对任一块体积 表面为 作用在 上的 体积力、惯性力和面力的合力必须为零。 体积力、惯性力和面力的合力必须为零。即