有限元法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点,并以其法向建立笛卡尔坐标系。
这样的截面有无数多对,这样的坐标系也有无数多个。
有限元_塑性力学
rg
连续体
三结点三角形单元
六结点三角形单元
矩形单元 任意四边形单元 8结点曲边 四边形单元
2. 单元分析
建立: {F}e=[k]{d}e
[k]:单元刚度矩阵
y vj j
vi , (Vi) i ui , (Ui)
uj
vm
m um
x
单元结点力向量{F}e=[Ui Vi Uj Vj Um Vm]T
利用单元的结点位移将整个单元的位移分量表
示为坐标的函数。
y
i
2. 三结点三角形单元的位移模式 j P
设:u=a1+a2 x+a3 y v=a4+a5 x+a6 y
m x
系数a1~a6由结点位移 ui , vi , uj , vj , um , vm确定。
• 将位移模式写成结点位移的显式:
i
u= Niui+ Njuj +Nmum
单元结点位移向量{d}e=[ui vi uj vj um vm]T
•
体力、面力
静力等效 ——
等效结点荷载
3. 整体分析
建立: {F}=[K]{d}, [K]:整体刚度矩阵
由各结点平衡{F}={R},得有限元方程:
[K]{d}= {R}
0.3 位移模式与解答的收敛性
1. 什么是位移模式(位移函数)
下标1,2,3表示坐标 x1, x2, x3(即x, y, z)方向.
2. 一点斜面上的应力
3
SNi i1l1 i2l2 i3l3 ijl j j1 (i 1, 2, 3)
x3 N
SN3
SN
SN1 O
用张量分析方法推导含偶应力弹性力学有限元理论
河 北 水 利 电 力 学 院 学 报JournalofHebeiUniversityof WaterResourcesandElectricEngineering2021 年3 月第31卷第1期Mar2021Vol31 No1文章编号:2096 — 5680(2021)01 — 0075 — 06用张量分析方法推导含偶应力弹性力学有限元理论孙晓勇1 2 ,宋兴海2,侯娜娜12,付建航2,刘立悦1,2(1.河北省数据中心相变热管理技术创新中心,河北省沧州市重庆路1号061001;2.河北水利电力学院土木工程学院,河北省沧州市重庆路1号061001)摘要:经典弹性力学理论用位移梯度表示无限小变形,不考虑旋转变形,把微元体的旋转视为刚体旋转。
含偶应力弹性力学理论将旋转变形以旋转张量表示,微元体旋转和微元体平动位移同量级,而旋转张量和应变张量同量级,旋转张量与旋转矢量一一对应,用旋转矢量的梯度表示旋转变形。
含偶应力弹性力学理论本构关系包括应力-应变关系和偶应力-曲率张量关 系,用等参变换方法离散单元位移到节点上,从虚功原理出发,增加罚函数项以降低有限元方程对高阶单元的需求,推导了拟 解决三维及二维问题的含偶应力弹性线力学有限元理论,可得三维及二维问题中位移、应力、应变等分布情况,对结构进行力 学评价。
关键词:偶应力;旋转变形;旋转张量;张量分析中图分类号:O343文献标识码:A DOI : 10. 16046/j. cnki. issn2096-5680. 2021. 01. 0151经典线弹性理论与考虑偶应力线弹 性理论在经典弹塑性力学理论中,物体内任意一点的 应力状态只和应变或应变的历史有关,其基本变量为位移,对位移求梯度得到应变张量,用位移梯度描述无限小的变形,然后再由一点的应变张量分析得 到应力张量[1]。
含偶应力的线弹性力学理论认为, 物体内任意一点的微元体,除有各个方向的位移外,还有本身的旋转变形,而这种旋转变形并非单纯的 以旋转角表达,而是用和应变张量一个量级的旋转张量来表示[]。
有限元三大方程
有限元三大方程有限元方法(Finite Element Method, FEM)是一种将连续介质离散为有限数量的小单元来求解力学问题的数值方法。
有限元方法通过三大方程来描述系统的力学行为:平衡方程、运动方程和本构关系。
1. 平衡方程:平衡方程是描述系统在受力平衡状态下的行为。
对于一个连续体,平衡方程可以用微分形式表示为:∇·σ + f = 0其中,∇·表示散度算子,σ是应力张量,f是体力。
在有限元方法中,将连续体离散为小单元后,平衡方程可以用积分形式表示为:∫(∇·σ)dV + ∫fdV = 0这里积分是对整个区域求和,dV表示体积元。
有限元方法的目标是通过对小单元的离散近似求解平衡方程来得到连续体的应力分布。
2. 运动方程:运动方程描述了系统的运动行为。
对于一个静力学问题,运动方程为:ρ∂²u/∂t² + ∇·(ρv) = 0其中,ρ是介质的密度,u是位移场,v是速度场。
在有限元方法中,将连续体离散为小单元后,运动方程可以写为:∑(∫(ρN∂²u/∂t²)dV) + ∑(∫(N∇u·ρvdV)) = 0这里,∑表示对所有小单元求和,N是形状函数。
有限元方法的目标是通过对小单元的离散近似求解运动方程来得到连续体的位移场。
3. 本构关系:本构关系描述了物质的力学性质。
对于一个线弹性材料,本构关系为:σ = Eε其中,σ是应力张量,E是弹性模量,ε是应变张量。
在有限元方法中,将连续体离散为小单元后,本构关系可以写为:∑(∫(NσdV)) = ∑(∫(B∇u)dV)这里,B是形状函数的导数矩阵。
有限元方法的目标是通过对小单元的离散近似求解本构关系来得到连续体的应力分布。
有限元方法通过将连续体离散为小单元,并在每个小单元内近似求解平衡方程、运动方程和本构关系,来得到连续体的力学行为。
通过不断迭代小单元的解,最终可以得到整个体系的力学行为。
第二章第2节作用于流体的力应力张量
24
4、牛顿流体 流体中流点的应力和变形速度间的关系满足广义
牛顿公式(2.36)的流体称为牛顿(粘性)流体。 如水和空气。
还有一些流体不满足(2.36)式,称为非牛顿流体, 如颜料、低温下的润滑油等,这些属于胶体化学。
25
总结
26
精品课件!
27
精品课件!
28
End
29
携手共进,齐创精品工程
3
问题
那么,要描写某一点的应力就需要知道所有通过该点 的面上所受的应力。-------是否一定要这样做? ----------不必,后面就会看到,过同一点不同面上所受到的 应力并不是处处相互独立,事实上,只要知道三个与坐标 面平行面上的应力,则任一以 为法向的面上的应力都可 以通过它们及 表示出来。
随着受力面元取向的不同而变化,即:
是空间某一点的位置, 是该点某一个受力面元的法向单位矢。
这段话可以这样理解:流体中有各个位置的点,不同点用 确定, 对于某一点 ,过这一点可以做无数个不同方向的面元,这些 面元就用 区别开来了,作用在这些面元上的面力 一般来说是不 同的,因此, 是 位置 和表面法向 的函数了,另外还 随着时间变化。
pzx
pxy pyy pzy
p p px y zzzzp1 0 01 0 01 0 02exx e ez1 3 yd xxxV iveyy e e1 3 x zd yy V ivezz e e1 3 x ydzz V iv
22
p px yx x
pzx
pxy pyy pzy
p p px y zzzzp1 0 01 0 01 0 02exx e ez1 3 yd xxxV iveyy e e1 3 x zd yy V ivezz e e1 3 x ydzz V iv
有限元法理论基础弹力变分原理教学内容
ij
Dijkl kl
W
ij
ij
Cijkl kl
➢ 对于线弹性材料,有:
其中:
Cijkl
D1 ijkl
W W ijij2
7
有限元法理论基础
应变能和应变余能
➢ 互余关系:
W W ijij
全功
➢ 对线弹性体应变能密度和应变余能密度可分 别表示为应变和应力的二次齐函数: W 12Dijkl ijkl
W 12Cijkl ijkl
8
有限元法理论基础
虚功原理
➢ 变形体虚功原理:变形体中,任意平衡力系 在任意满足协调条件的变形状态上作的虚功 等于零。即体系外力的虚功等于内力的虚功。
➢ 虚功原理:
➢ 虚位移原理、虚应力原理
9
有限元法理论基础
虚位移原理
➢ 虚位移:变形体几何约束所允许的位移称为可 能位移,取其任意微小的变化量即是虚位移。
15
有限元法理论基础
虚应力原理
➢ 虚应力:满足平衡方程和力边界条件的应力的变分 (微小的变化)。
➢ 虚应力原理:位移边界处给定位移在虚反力上所做的 余虚功等于应变在虚应力上的余虚功:
uipidS ijijdV
➢ 按力法建立有限元方Su 程的基V本方程。
➢ 几何方程和位ij移12边(u界i,j 条u件j,i)
有限元法理论基础
虚位移原理
➢ 几何方程:
ij 12(ui,j uj,i)
➢ 将应变的变分及上述分部积分结果代入等效积 分形式,得到等效积分形式:
V ijid j VS u iT id SVu ifidV
➢ 上式右边是变形体内的应力在虚应变上所作的虚功,即 内力虚功;上式左边是外力(体力+面力)在虚位移上 所做功,即外力虚功。即内力虚功等于外力虚功。
张量和应力张量PPT课件
1.1 角标符号
• 带有下角标的符号称为角标符号,可用来表 示成组的符号或数组。
• 例:
– 直角坐标系的三根轴
• x、y、z→x1、x2、x3 → xi(i=1,2,3);
– 空间直线的方向余弦
• l、m、n → lx、ly、lz → li(i=x,y,z);
– 表示一点应力状态的九个应力分量
• σxx、σxy… → σij(i,j=x,y,z);
a11b32
x13
a12b32 x23
a13b32
x33
y13
y33
1.3 张量的基本概念
• 只需一个实数就可以表示出来简单的物理量称为 标量。例如距离、时间、温度等。
• 需用空间坐标系中的三个分量来表示的物理量称 为矢量。 例如位移、速度、力等。
• 对于复杂的物理量,例如应力状态、应变状态等, 需要用空间坐标系中的三个矢量(也即九个分量) 才能完整地表示出来,这就是张量。
• 张量是矢量的推广,与矢量相类似,可以定义为: 由若干个当坐标系改变时满足转换关系的分量所 组成的集合称为张量。
• 物理量P
– 在j=1空,间2,坐3标);系xi (i=1,2,3)中存在九个分量Pij (i,
– 在新空间坐标系 xk(k=1’,2’,3’)中存在九个新 分量Pkr(k,r=1’,2’,3’)。
1.5 应力张量
• 外力确定后,受力物体内任意点的应力状态即已 确定。但表示该点应力状态的各个分量在不同坐 标系中将有不同的数值,因此在不同坐标系中该 点的应力分量之间应该存在一定的关系。
• 设受力物体内一点的应力状态为: – 在xi(i=x,y,z)坐标系中为σij(i,j=x,y,z);
– 在xk (k=x’,y’,z’)坐标系中为σkr(k,r=x’,y’,
弹性理论第二讲—张量理论_211207049
如何表示???
我们在物体内兴趣点周围截出一微单元,该点的 应力状态由各微分面上的应力情况表示。 应力状态由各微分面上的应力情况表示
哑指标:该重复指标称为哑指标。 自由指标:除哑标外,在表达式或者方程的某项中 非成对出现的其它指标。
a11 x1 a12 x2 a13 x3 x1
10
附录:指标符号
求和约定
a11 x1 a12 x2 a13 x3 x1
利用哑标简化方 a21 x1 a22 x2 a23 x3 程中的表达式 x2
u e
i 1
3
i i
u1e1 u2 e2 u3e3
9
附录:指标符号
求和约定 求 约定
定义:如果在表达式的某项中,某指标重复 重复地出现 两次,则表示要把该项 两次 该项在该指标 该指标的取值范围 取值范围内遍历 遍历 求和。 求和
u e
i 1
3
i i
u1e1 u2 e2 u3e3 ui ei
19
eijk ai b j e k
A- 2
二者在矢量运算中的应用
[ a , b, c ] a b c a b c
三个矢量的混合积: a am e m , b b j e j , c ck e k ,
a b c (am e m ) (b j e j ck e k )
弹性力学基础及有限元
第二讲 张量理论
有限元法2011应力张量
V
V
[r ( σ f u&&) kjekjiei ]dV 0
V
0
kjekjieidV 0
V
V的任意性
kjekjiei ( 23 32 )e1 ( 31 13 )e2 (12 21)e3 0
即
ij ji
(13)
应力张量是对称的,九个应力分量中只有六个是独立的。这一结
2 2j n j
4 n
jnj
2n
j
(
2 j
2
j n )
式(c)为 即
n j
(nini 1) n j
2nj
nj
(
2 j
2
j
n
)
0
(对j不求和)
n1
(
2 1
21
n
)
0
n2
(
2 2
2
2
n
)
0
(d)
n3
(
2 3
2 3
n
)
0
n是单位矢量, 若ni全不为零
12 21n 0,
2 2
2
2
n
0,
2 3
2 n
12 (n12
n22
)
2 3
n32
[1(n12
n22 ) 3n32 ]2
12 (1
n32 )
n2 2
33
[1 (1
n32 )
3n32 ]2
(a)x3令ຫໍສະໝຸດ 2 n0得
n3 1 2
45 45
n3
由 nini 1
n12 n22 1 2
x2
n
1
3
2
(完整版)应力张量的认识(一)
应力张量的认识(一)本文主要是对材料成形相关专业学习过程中对一些问题的思考,也许并不深刻,但却是自己从初学时的迷惑到后来逐渐认识的过程。
相关还有:Levy-Mises理论的思考从本科的材料成形原理教材上就认识了应力张量,然后一直出现在我们的视野里。
初始,以一个基本定义记住了它,进而有过疑惑,随着矩阵论的学习又有了新的认识。
曾经就有记录下对其理解的想法,但因思路尚未完善而一再搁置;直到今天重新想起,完成了方向余弦作为线性空间的证明,才终于开始详细记录。
我将这部分思考分为以下三部分:应力张量的认识(一)应力张量的认识(二)应力张量的认识(三)本文介绍第一部分应力的基本知识和常规认识。
应力初中物理就已知道,因外力作用而在物体内部产生的力成为内力。
单位面积上的内力即是应力,表征内力的强度。
为了研究某一点P处的应力,用某个截面在P点处切开物体,如下图所示。
根据定义可以得到P点的正应力σ、切应力τ,他们的合成即为全应力T。
需要注意的是,一个确定的截面对应了一组正应力和切应力。
但是过P点有无数的截面,那么如何才能真正描述P点的应力状态呢?应力状态点的应力状态是受力物体内某一点各个截面上应力的变化情况。
上面已经意识到过一点点有无数的截面,只有任意截面上的应力分量都可以确定,才可以说应力状态是确定的。
通常在无数的截面中,任意取三个互相垂直的截面,并以他们的法线方向建立笛卡尔坐标系。
也即在P点截取一个无限小的平行六面体,称为单元体。
单元体无限小,视为一点,因此单元体上相互平行的两个平面视为过该点的同一平面,也即他俩的应力是相同的。
这样就只用三个互相垂直的截面上的应力来分析问题。
由于单元体处于静力平衡状态,由绕各轴合力矩为零可以得到切应力互等定律。
问题:既然单元体上相互平行的两个平面视为过该点的同一平面,那为什么上图平行的平面上应力是相反的?单元体上相互平行的两个平面视为过该点的同一平面,但是分别是被截开的的两部分的平面,截开前他们是重合的,截开后成为了两部分各自的表面,而外表面是有方向的。
应力张量
从一点应力状态分析,Cauchy应力张量描述的 是变形构形中点的应力状态,而Lagrange、 Kirchhoff则是描述初始构形中点的应力状态,不过 这个点是由变形构形中的点对应过来的,故是特定 的。 只有Cauchy应力张量真实反映了客观物理现象。 而Lagrange、Kirchhoff应力张量只是为了便于应用 的“近似计算”。
dT F dT
K
1
三种应力张量的区别和联系
Cauchy应力张量定义于Euler描述,应力是空间 坐标的函数,力矢作用在变形后的面元上,由此算 得的应力就是真实应力,但由于应用时边界条件常 用物质坐标给出,用空间坐标给出边界条件十分困 难,加之本构方程也常采用Lagrange描述,故用 Cauchy应力张量求解边值问题还须作进一步变换。谢谢!源自Kirchhoff应力张量
在大变形问题中,是从变形后的物体中截取出 微元体建立平衡方程和与之相等效的虚功原理,如 果应变是用变形前的坐标表示的Green应变张量,则 需要定义与之对应的关于变形前位形的应力张量。
Kirchhoff规定:变形前面积元上的内力分量和 变形后面积元上的内力分量的变换与标变换一致。
应力张量的概念及其应用PPT课件
广义胡克定律,应变能密度 应变能密度
vd
1 6 E (1 2 )2 (2 3 )2 (3 1 )2
v v
16 E 2(123)2
vd vv v
重要应用实例
承受内压薄壁容器任意点的应力状态
重要应用实例
m t l
m
m(p D)
D
m
p
pp D 2
4
D
pDl
p
t t (2 l ) t
应力张量的概念及其应用
1。应力张量及其不变量 2。应变张量及其不变量 3。广义胡克定理
1。应力张量及其不变量
一、张量的概念
自然界的物质的性质和规律是一种客观存在, 不受描述它的方法的影响。
数学方法描述时,会引入坐标系。不同的坐 标系的选择,会使问题简单化或复杂化。
希望有某些数学量在描述物理现象时,尽量 摆脱具体坐标系的影响。
应力和应变是二阶张量
1。应力张量及其不变量
二、一点的应力状态表示
用二阶张量在x, y, z 坐标系表示
xx yx
xy yy
xz yz
zx zy zz
或写成:
xx yx
xy yy
xz yz
zx zy zz
1。应力张量及其不变量
仍选用直角坐标系,坐标轴写成 x1, x2, x3 采用张量下标记号法:
)
广义胡克定律,应变能密度
应变能密度
广义胡克定律,应变能密度 应变能密度
1、微元应变能(Strain Energy)
2
1dydz~1dx
dy
1 2dxdz~2dy
dz 3 dx
3dydx~3dz
广义胡克定律,应变能密度 应变能密度
- 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,作用在 动量定理和动量矩定理。 对任一块体积 表面为 作用在 上的 体积力、惯性力和面力的合力必须为零。 体积力、惯性力和面力的合力必须为零。即