Kirchhoff应力张量
- 格式:pptx
- 大小:334.35 KB
- 文档页数:8
应力张量的认识(一)本文主要是对材料成形相关专业学习过程中对一些问题的思考,也许并不深刻,但却是自己从初学时的迷惑到后来逐渐认识的过程。
相关还有:Levy-Mises理论的思考从本科的材料成形原理教材上就认识了应力张量,然后一直出现在我们的视野里。
初始,以一个基本定义记住了它,进而有过疑惑,随着矩阵论的学习又有了新的认识。
曾经就有记录下对其理解的想法,但因思路尚未完善而一再搁置;直到今天重新想起,完成了方向余弦作为线性空间的证明,才终于开始详细记录。
我将这部分思考分为以下三部分:应力张量的认识(一)应力张量的认识(二)应力张量的认识(三)本文介绍第一部分应力的基本知识和常规认识。
应力初中物理就已知道,因外力作用而在物体内部产生的力成为内力。
单位面积上的内力即是应力,表征内力的强度。
为了研究某一点P处的应力,用某个截面在P点处切开物体,如下图所示。
根据定义可以得到P点的正应力σ、切应力τ,他们的合成即为全应力T。
需要注意的是,一个确定的截面对应了一组正应力和切应力。
但是过P点有无数的截面,那么如何才能真正描述P点的应力状态呢?应力状态点的应力状态是受力物体内某一点各个截面上应力的变化情况。
上面已经意识到过一点点有无数的截面,只有任意截面上的应力分量都可以确定,才可以说应力状态是确定的。
通常在无数的截面中,任意取三个互相垂直的截面,并以他们的法线方向建立笛卡尔坐标系。
也即在P点截取一个无限小的平行六面体,称为单元体。
单元体无限小,视为一点,因此单元体上相互平行的两个平面视为过该点的同一平面,也即他俩的应力是相同的。
这样就只用三个互相垂直的截面上的应力来分析问题。
由于单元体处于静力平衡状态,由绕各轴合力矩为零可以得到切应力互等定律。
问题:既然单元体上相互平行的两个平面视为过该点的同一平面,那为什么上图平行的平面上应力是相反的?单元体上相互平行的两个平面视为过该点的同一平面,但是分别是被截开的的两部分的平面,截开前他们是重合的,截开后成为了两部分各自的表面,而外表面是有方向的。
基于Odgen模型O型橡胶密封圈的大变形接触分析桑建兵;邢素芳;刘宝会;王静远;刘春阳【摘要】By utilizing nonlinear finite element analysis softwareMSC.Marc,the material parameters of Odgen model were determined through curve fitting by using the experimental data of uniaxial tensile experiments of rubber-like materi-als.The nonlinear finite element model of rubber O-ring seals was set up and the influence of different medium pressure on the mechanical properties of rubber O-ring seals was analyzed.The distribution of Von-mises stress and contact stress on primary contact surface and side contact surface was achieved.The results show that with the increase of medium pressure, the Von-mises stress is increased and transfered to the contact zone between rubber seal and groove,and its distribution on primary contact surface and side contact surface is approximate quadratic parabola.The maximum contact stress appears on the middle contact zone,and it is increased obviously with the increases of medium pressure.The maximum contact stress is greater than medium pressure,which can prevent the leakage of medium.%利用大型非线性有限元软件MSC.Marc,基于橡胶类材料单轴拉伸所得到的应力应变曲线,通过数据拟合确定Odgen模型的材料参数。
简析非线性有限元法游潇;苏小卒【摘要】In the analysis of reinforeced concrete structures subjected to general loading conditions,the realistic constitutive model and robust analytical procedure are two key preconditions to produce reasonably accurate simulations of nonlinear behaviors of such structures.Based on the FEM analysis,suggestions for further studies are given.%采用有限元法分析一般荷载作用下的钢筋混凝土结构时,要得到对结构性能的合理准确的模拟结果,除了需要合理的本构模型,还要有先进的数值分析方法。
文中将对非线性有限元的特点做出分析,为进一步研究提供参考。
【期刊名称】《江西科学》【年(卷),期】2012(030)001【总页数】4页(P75-78)【关键词】非线性有限元;数值方法;钢筋混凝土【作者】游潇;苏小卒【作者单位】同济大学土木工程学院,上海200092;同济大学土木工程学院,上海200092【正文语种】中文【中图分类】TU311目前,钢筋混凝土结构作为一种经济、实用的结构,是我国工业与民用建筑中最为广泛采用的一种结构形式。
这类结构在各类荷载作用下的反应特性,以及合理的设计方法和构造措施,一直以来是结构研究人员和工程师们经常研究的课题。
由于钢筋混凝土是由2种性质截然不同的材料—钢筋和混凝土组合而成,因此它的性能明显的依赖于这2种材料的性能。
尤其是在非线性阶段,钢筋和混凝土本身的各种非线性性能,都不同程度地在这种组合材料中充分反映出来。
目前钢筋混凝土非线性方面的分析和研究还存在着若干有待解决的问题。
橡胶材料硬化的本构模型与有限元分析朱艳峰;王红【摘要】针对橡胶类材料大应变时硬化现象,采用国家GB528标准,在室温下通过单轴拉伸本构实验,建立了基于主伸长的连续介质力学的新本构模型,并确定应变能密度函数中的本构参数,再利用简单剪切实验进行参数验证,表明新本构关系的可行性与有效性.最后将新本构方程加入通用有限元软件,利用非线性有限元对平面应力橡胶板进行了计算.【期刊名称】《武汉工程大学学报》【年(卷),期】2008(030)001【总页数】3页(P34-36)【关键词】橡胶类材料;材料硬化;本构模型;非线性有限元【作者】朱艳峰;王红【作者单位】广东工业大学建设学院工程力学研究所,广东,广州,510640;广东工业大学建设学院工程力学研究所,广东,广州,510640【正文语种】中文【中图分类】O3450 引言橡胶类材料产生大的应变时具有高度几何非线性、材料非线性,又呈现出硬化或软化现象,且体积不可压缩.橡胶材料的本构模型主要有[1]:a.基于分子链形式的统计学模型.b. 基于不变量形式的模型.c.基于主伸长的连续介质力学模型等.针对大应变硬化现象的本构模型目前有neo-Hookean型[2]、Mooney-Rivlin[3]型等.在对橡胶类材料的有限元分析过程中,由于其本构理论尚未成熟,导致分析结果的差别非常大.本文通过简单拉伸本构实验确定材料变形模式,用回归分析方法把实验得到的应力-应变数据拟合为一适当的应变能函数,建立了一种新的橡胶材料硬化模型,并推广到复杂的变形形式,通过验证并加入通用有限元软件,针对平面应力橡胶板进行了计算,为橡胶材料硬化时进一步的本构研究提供了基础.1 橡胶材料硬化的本构实验橡胶材料的本构实验试件采用硫化橡胶,在160℃下保持10 MPa压力,硫化6 min,配方及质量比如表1所示.表1 试件配方及质量比Table 1 The test-piece in parts by weight湛江农垦局天然橡胶(1#)SZnOSA防264N330CZTT100%1.0%4.0%2.0%1.5%3%1.0%0.2%1.1 单轴拉伸实验试件为哑铃型,如图1所示.规格尺寸符合GB528标准,室温26℃,采用岛津电子拉伸实验机,拉伸速度为(500±50) mm/min,共5组试件,试件拉力与伸长比实验曲线如图2.图1 单轴拉伸试件Fig.1 Rubber test piece for simple tension图2 拉应力t-伸长率λ曲线Fig.2 The curve of stress t versus extension1.2 简单剪切由于没有国家标准,参照Treloar1943年的实验[4],矩形试件,宽60 mm,厚2 mm,测试长度5 mm,用夹具分别夹住长边,室温160℃,拉伸速度为(500±50)mm/min,实验数据如表2.表2 简单剪切的最大应力与应变Table 2 The maximal stress and strain on simple shear试件最大应力/N·mm-2最大应变/%15.3333493.3325.6667543.3335.1333486.67平均5.3777507.782 橡胶类材料硬化时的新本构模型常温条件下,橡胶类材料为各项同性超弹性材料,本构模型以应变能函数的形式来表示,I1,I1,I3为右Cauchy-Green变形张量的第一、二、三基本不变量,在初始无应力构形且不考虑大应变硬化时应变能函数W可表示为[5]W=W(I1,I2,I3)I1=trC=C∶I=Ci iI3=detC.橡胶类材料在变形过程中近似认为体积不可压缩,变形后与变形前的体积比J=1.由本构实验,针对大应变时呈现出的明显硬化现象,本文提出一种新的应变能函数:W(λ1,λ2,λ3,φ)=μ(I1-3)(I2-3)+H(φ)φ=(λ*-λ1)(λ*-λ2)(λ*-λ3)λ1、λ2、λ3为主伸长,λ*为极限伸长.即φ为无穷时,又回到应变能不考虑硬化的传统形式,为所谓的根应变能.单轴拉伸时,第一类Piola-kirchhoff应力f=对单轴拉伸试验数据进行拟合,得:μ=0.004 8 MPa,λ*=9.47, k=9 233由此绘出拉伸曲线,如图3所示;将上述参数代入新本构方程,计算简单剪切时的f-λ曲线与实验值进行比较如图4所示.可知,本文提出的本构关系能够较好反映橡胶材料在单轴拉伸时的硬化现象,且与剪切硬化实验吻合良好,故可适用于对此类橡胶材料硬化进行力学分析.图3 单轴拉伸曲线Fig.3 The curve of simple extension图4 简单剪切曲线Fig.4 The curve of simple shear3 有限元分析第二类Piola-Kirchhoff 应力张量S与Green应变张量E存在下列关系与分别为应力张量与应变张量的率形式,因此⊗NαNα为原始构形中沿主方向的正交单位向量.本构方程中的Lagrangian弹性张量C=⊗Nα⊗Nβ⊗Nβ+⊗Nβ⊗Nα⊗Nβ因此,对于新本构函数C=Nα⊗Nα⊗Nβ-⊗Nβ⊗Nα⊗Nβ4 算例本文利用通用有限元软件,加入新本构方程,对矩形均匀伸长的橡胶薄板进行了计算,薄板长100 mm,宽50 mm,厚2 mm,中心开孔直径5 mm,近似认为平面应力.计算时利用二阶平面应力减缩积分单元,1/4橡胶薄板的网格划分与位移、应力计算结果如图5,6,7所示,此时,板伸长为15.52 mm,最大应力0.6 MPa.图5 中间开孔的1/4橡胶薄板网格划分Fig.5 Rubber elastic sheet with a circular hole - the geometry and the mesh for a quarter-sheet图6 位移分布图Fig.6 Final displaced configuration of the quarter-sheet 图7 应力分布图Fig.7 Final stress distributing of the quarter-sheet参考文献:[1]朱艳峰,刘锋,黄小清,等.橡胶材料的本构[J].橡胶工业.2006,53(1):119-125.[2]Horgan, Saccomandi. Constitutive modeling of rubber-like and biological materials with limiting chain extensibility[J]. Mathematics and mechanics of solids, 2002,(7):353-371.[3]朱艳峰,刘锋,黄小清,等.橡胶类材料大应变时明显硬化的本构分析[J].暨南大学学报.2005,26(1):98-99.[4]Treloar L R G. Stress-Strain data for vulcanized rubber under various of deformation[J]. Trans Faraday Soc, 1944,40(6):59-70.[5]Ogden R W. Non-Linear Elastic Deformations[M]. Chichester, UK:Ellis Horwood, 1984.。
第二皮奥拉-基尔霍夫应力张量第二皮奥拉-基尔霍夫应力张量是一个用于描述非均质连续介质内部应力分布的张量。
它通常用于解决固体力学、弹性力学和流变学等领域中的问题。
该张量以法向和切向压力的梯度为基础,用以表示介质内在应力场的行为。
而二阶张量由于其两个指标,能够同时考虑空间方向和应力的耦合,更加适用于对于应力场的描述。
一般而言,第二皮奥拉-基尔霍夫应力张量可以表示为:σ_ij = -pδ_ij + 2με_ij其中,σ_ij代表某一时刻物质内部的应力,p代表介质的压力,δ_ij代表克罗内克(Kronecker)δ符号,它用于将一个索引位置的张量分量与另一个相同的索引位置分量相乘并将它们加起来。
ε_ij代表应变张量,而μ则代表剪切模量。
实际上,这个方程中的前半部分代表了压力对于应力的贡献,而后半部分则是介质变形对于应力的贡献。
比如说,如果一个物质受到压缩,那么压力将导致介质扭曲,从而会改变物质内部的应力分布。
因此,这个张量的作用就在于将不同因素的贡献整合到一个方程中,从而准确地描述物质内部应力的行为。
需要注意的是,这个方程中只考虑了刚性体(没有流动性)中的应力情况。
如果涉及流动性质的话,就需要考虑黏度和非牛顿性等因素。
总而言之,第二皮奥拉-基尔霍夫应力张量是描述非均质连续介质内部应力分布的一个重要工具。
它能够将压力和介质变形对于应力的影响综合起来,从而准确地描述整个应力场的行为。
它在多个学科领域中都有广泛的应用,包括工程力学、地质学、流体力学等。
通过对于这个张量的深入理解,我们可以更好地研究介质内部应力的行为,为如何合理地利用自然资源,设想更高效的新材料等领域,提供更为深入的基础。
质量与检测聚氯乙烯膜材各向异性超弹性本构模型安晓通(新疆天域新实化工有限公司,新疆石河子832000)摘要:为有效且充分地反映和揭示聚氯乙烯膜材料在拉剪耦合应力作用条件下发生的大变形特性、非线性特性,以及各向异性等多样化的力学特性,本项研究基于纤维增强复合材料连续介质力学基本理论,借由将应变能解耦成经向纤维束与纬向纤维束的拉伸变形能,以及两向纤维束之间的角度变化产生的剪切变形能,建构形成了针对聚氯乙烯膜材的各项异性本构模型。
本文围绕聚氯乙烯膜材各向异性超弹性本构模型论题,择取两个具体方面展开了简要阐释。
关键词:聚氯乙烯膜材;各向异性;超弹性本构关系;拉剪耦合从材料技术特性角度展开阐释分析,聚氯乙烯膜材具备材料非线性特点,拉伸技术过程中编织结构卷曲转换几何非线性特点,以及能促使各组分材料之间在协同受力条件下摩搽接触非线性特点,在力学行为方面,该种材料整体上呈现大变形、非线性和各向异性的技术特点。
在现阶段常规组织开展的工程设计活动过程中,通常倾向于将聚氯乙烯膜材假设成一种线弹性材料类型,由于在主要技术参数项目的选取和确定过程中不能全面充分地反应和迎合上述的力学技术行为特征,在一定程度上给膜结构设计方案形成过程中的可靠性造成了不良影响。
有鉴于上述研究背景,本文将会围绕聚氯乙烯膜材各向异性超弹性本构模型论题,展开简要阐释。
1聚氯乙烯膜材本构模型已有研究回顾针对聚氯乙烯膜材本构模型展开的研究分析,主要基于宏观表征模型,以及细观机理模型两个具体方面加以展开。
细观机理模型在形成过程中,主要将聚氯乙烯膜材的细观结构特征和变形机理作为基点,在纱线模拟基础上建构形成编织结构的数值模型或者是理论模型,考量到经向纤维与纬向纤维之间存在的相互作用因素、纤维材料发生的卷曲变形因素、纤维材料与涂层之间发生的伸长因素,以及纤维材料与涂层之间的相互摩擦等,具体建构材料细观本构关系模型,之后遵照实验过程中获取的相关数据,反向推导获取模型中包含的各项相关参数。
kirchhoff应力转换cauchy应力
Kirchhoff应力和Cauchy应力都是表示材料中的应力分布的数学概念。
它们之间的转换关系可以通过应力张量的变换规律来描述。
在固体力学中,应力是描述物体内部受力情况的物理量。
通常用一个二阶的张量来表示应力。
在三维空间中,应力张量可以被表示为一个3x3的矩阵。
Kirchhoff应力是一种测量物体在不同形变状态下的应力的方法。
它是在物体的当前形变状态下测量的应力。
Kirchhoff应力张量可以用物体的当前形变矩阵直接计算得到。
Cauchy应力是一种测量物体在其未发生形变之前的初始应力的方法。
它是在物体未发生形变之前测量的应力。
Cauchy应力张量可以用物体的位移梯度矩阵及其转置矩阵计算得到。
Kirchhoff应力和Cauchy应力之间的转换关系可以用下面的公式表示:
Kirchhoff应力 = J * Cauchy应力 * J^T
其中,J是物体的形变梯度矩阵。
需要注意的是,Kirchhoff应力和Cauchy应力描述了同一个应力场,只是基于不同的测量方法和参考系来表示。
因此,它们的数值大小和分布可能会有所不同,但包含的物理信息是相同的。
柯西应力张量动量平衡方程
柯西应力张量与动量平衡方程是物理学中的两个概念。
柯西应力张量是一种描述物质点在运动过程中所受到的应力的量纲,通常用于描述流体力学中的应力状态。
其定义为Σ(a)(t)=12∑b:b≠aF(ab)(t)(x(b)(t)−x(a)(t)),其中F(ab)(t)表示物质点a和b之间的相互作用力,x(a)和x(b)表示物质点a和b的位置。
动量平衡方程则是描述一个封闭系统中的动量变化率等于该系统所受到的合外力。
在连续介质力学中,动量平衡方程可以表示为ρ*(udt)dt=∫(Tj)eρ*(u·Je)dVdt,其中ρ*表示质量密度,u表示速度矢量,t表示时间,∫(Tj)e表示对整个系统的积分,Je表示单位体积的体力密度。
柯西应力张量和动量平衡方程在流体力学、气体力学、固体力学等领域都有广泛的应用。