最新大变形问题的有限元分析

合集下载

有限元实验报告

有限元实验报告

有限元实验报告一、实验目的本实验旨在通过有限元方法对一个复杂的工程问题进行数值模拟和分析,从而验证理论模型的正确性,优化设计方案,提高设计效率。

二、实验原理有限元方法是一种广泛应用于工程领域中的数值分析方法。

它通过将连续的求解域离散化为由有限个单元组成的集合,从而将复杂的偏微分方程转化为一系列线性方程组进行求解。

本实验将采用有限元方法对一个具体的工程问题进行数值模拟和分析。

三、实验步骤1、问题建模:首先对实际问题进行抽象和简化,建立合适的数学模型。

本实验将以一个简化的桥梁结构为例,分析其在承受载荷下的应力分布和变形情况。

2、划分网格:将连续的求解域离散化为由有限个单元组成的集合。

本实验将采用三维四面体单元对桥梁结构进行划分,以获得更精确的数值解。

3、施加载荷:根据实际工况,对模型施加相应的载荷,包括重力、风载、地震等。

本实验将模拟桥梁在车辆载荷作用下的应力分布和变形情况。

4、求解方程:利用有限元方法,将偏微分方程转化为线性方程组进行求解。

本实验将采用商业软件ANSYS进行有限元分析。

5、结果后处理:对求解结果进行可视化处理和分析。

本实验将采用ANSYS的图形界面展示应力分布和变形情况,并进行相应的数据处理和分析。

四、实验结果及分析1、应力分布:通过有限元分析,我们得到了桥梁在不同工况下的应力分布情况。

如图1所示,桥梁的最大应力出现在支撑部位,这与理论模型预测的结果相符。

同时,通过对比不同工况下的应力分布情况,我们可以发现,随着载荷的增加,最大应力值逐渐增大。

2、变形情况:有限元分析还给出了桥梁在不同工况下的变形情况。

如图2所示,桥梁的最大变形发生在桥面中央部位。

与理论模型相比,有限元分析的结果更为精确,因为在实际工程中,结构的应力分布和变形情况往往受到多种因素的影响,如材料属性、边界条件等。

通过对比不同工况下的变形情况,我们可以发现,随着载荷的增加,最大变形量逐渐增大。

3、结果分析:通过有限元分析,我们验证了理论模型的正确性,得到了更精确的应力分布和变形情况。

恒阻大变形锚杆冲击拉伸实验及其有限元分析

恒阻大变形锚杆冲击拉伸实验及其有限元分析

或15 m/s(撞击40次)的轴向初速度撞击托盘。

(7) 模型中所有接触面设定自动面接触。

(8) 其他相关参数设定。

上述前处理工作完成后投入LS-DYNA 软件进行计算并通过LSPREPOST 软件进行结果分析。

其中,网格模型缩略图和模型参数分别见图9和表3。

图9 1/4模型网格缩略示意图Fig.9 Schematic of partial meshing of 1/4 of model 表3 模型参数Table 3 Parameters of model介质 密度/ (kg ·m -3) 弹性模 量/GPa 泊松比屈服强度/MPa 单元 数量 单元尺寸节点数量冲击杆 7 850 206 0.300400 10 152 5 12 922托盘 7 850 206 0.300400 32 414235 380螺栓 7 850 206 0.300400 6 287 1,27 718恒阻套筒 7 800 210 0.269255 33 695 1 46 128杆体 7 850 206 0.300400 1 564 2 2 115子弹 7 850 206 0.300400 32 414 2 35 380 5 有限元结果分析与比较 5.1 膨胀量分析图10为单次冲击后选取的恒阻套筒外螺纹凹凸处参考单元示意图,图11为参考单元的径向位移(即膨胀量)曲线,其中,点A 为外螺纹凹处参考点,点B 为外螺纹凸处参考点。

图10 恒阻套筒参考单元示意图Fig.10 Sketch of selected elements of sleeve pipe0.00.10.20.30.40.50.60.70.80.9径向位移/m m时间/ms 图11 参考单元径向位移曲线Fig.11 Radial displacement of curves selected elements 从图11中可以看出:(1) 无论外螺纹凹处还是凸处,锚杆试样经受冲击后均表现出弹性回落,其中以凸处最为明显。

静力触探过程中的大变形弹塑性有限元分析

静力触探过程中的大变形弹塑性有限元分析

静力触探过程中的大变形弹塑性有限元分析摘要:本文利用有限元软件ABAQUS为工具,采用任意拉格朗日-欧拉(ALE)运动学方法描述,考虑土体的大变形特性,模拟静力触探过程中贯入阻力及沿探杆长度的应力场和位移场。

结果表明,贯入阻力随着贯入深度的增加线增大而后基本保持恒定,随土层深度的增加而增大,与土的内摩擦角和内聚力正相关。

通过对比大变形和小变形计算结果,表明对于静力触探,采用考虑大变形的ALE网格重划分技术可以较好的模拟锥尖处的应力集中现象。

关键词:静力触探;拉格朗日-欧拉(ALE)方法;大变形引言:本文以怀柔火工区铁路专用线一期工程、二期工程、张家口大唐电力铁路专用线工程、河北宣化钢铁集团铁路专用线工程、天津静海东方红路下穿京沪铁路立交桥工程的路基检测工程实践为课题基础,联合北京交大、首钢地质勘察院、北京金地通检测技术有限公司等多家单位进行研究,利用有限元软件ABAQUS为工具,采用任意拉格朗日-欧拉(ALE)方法描述的有限元方法,考虑土体的大变形特性,模拟静力触探过程中贯入阻力及沿探杆长度的应力场和位移场,并分析了各土性参数对应力场和位移场的影响。

1 有限元大变形问题的ALE描述简介在大变形分析中,位移与应变一般呈现非线性关系,使得以往的的小位移小应变分析方法不能适用于大变形分析计算。

目前,大变形分析的方法一般采用增量分析法,即在一般弹塑性分析方程中加入增量部分,随着加载变化而进行构型变化。

所以其求解方法以增量法为基础主要分为三种:Lagrange法、Euler法和Lagrange-Euler混合方法。

在拉格朗日描述(即物质描述或 L 描述)中,计算网格固定在物体上随物体一起运动,即网格点与物质点在物体的变形过程中始终保持重合,因此物质点与网格点之间不存在相对运动(即迁移运动,也称对流运动)。

对于物质的变形,一般采用应力和应变增量表示。

这大大地简化了控制方程的求解过程,而且能准确描述物体的移动界面,并可跟踪质点的运动轨迹。

有限元分析ppt

有限元分析ppt

分 片 近 似位
移 函 数
m(xm ym ) Fmy
vm um
vi i(xi yi )
Fmx ui
vj
y
Fix x
Fiy
uj
j(xj yj)
单 元 平 衡单
刚 方 程
整 体 平 衡总
刚 方 程


求 解
节 点 位



阶梯轴(梁)
A E (1)
(1)
A E (2) (2)
F
1
2
3
3
Φ1
Φ2
Φ3
l(1)
ui
vi
u
v
j j
um
vm
Fxi
Fyi
F
Fxj Fyj
Fxm Fym
y
vm
m
um vj
vi
j uj
i
ui
Fym
m
Fyi
i
Fxm Fyj
j Fxj Fxi
x
平面应变板单元
1.2.3 .1 单元刚度的概念 单元分析的主要工作是:通过研究单元力和单元位移
之间关系,建立单元刚度矩阵。 对任意单元而言,描述单元力和单元位移之间关系的
l(2)
F1
F2
F3
分为两个单元,共有三个节点。整体结构中,节点 载荷F及节点位移Φ都用大写。其脚标为节点在总体 结构中的编码,简称为总码。
1.1 有限元法概述
二.一个简单的应用实例
1. 离散化
① 局部码:各单元内,节点的编码; ② 各节点的位移分量及载荷分量分别用小写φ及f标记 ③ 所有节点位移的集合为该单元节点位移矢量{φ},节

压缩状态下橡胶件大变形有限元分析

压缩状态下橡胶件大变形有限元分析

文章编号:100021506(2001)0120076204压缩状态下橡胶件大变形有限元分析郑明军,谢基龙(北方交通大学机械学院,北京100044)摘 要:分析了橡胶硬度与橡胶力学常数C 1和C 2的一般关系,通过单向压缩试样试验和有限元计算,确定了C 1和C 2.在此基础上,研究了压缩状态下不同硬度橡胶支座的大变形特点,进一步探讨了C 1和C 2与硬度的关系.关键词:橡胶;力学常数;非线性有限元中图分类号:O631.21;O343.5 文献标识码:AFinite E lement Analysis of Large Deform ationof Compressed Rubber ComponentZH EN G M i ng 2j un ,X I E Ji 2long(College of Mechanical and Manipulative Engineering ,Northern Jiaotong University ,Beijing 100044,China )Abstract :This paper analyses the general relation between rubber hardness and rubber mechanicalconstant symbolized by C 1and C 2,which are determined through uniaxial tension test and finiteelement computation.On the basis of it ,the large deformation of compressed rubber supportingabout different hardness is researched and the relation between the rubber mechanical constantsand the hardness is further discussed.K ey w ords :rubber ;mechanical constant ;non 2linear finite element橡胶具有良好的弹性且容易变形,被广泛地应用载重结构的座架、弹簧、密封件、减震衬垫、联轴器和轮胎,然而由于橡胶材料的非线性、不可压缩性和大变形特性,使得描述橡胶力学特性的常数C 1和C 2的确定比较烦琐,一般采用实验的方法来得到[1].本文根据文献[2,3]的橡胶硬度与弹性模量关系的试验数据,得到了硬度与C 1和C 2的一般关系式,这样将两个待定常数减少为一个.在此基础上,采用有限元法计算了压缩状态下橡胶支座的载荷—变形曲线,与已有的试验数据[4]相比,表明本文的方法是可靠的.文中利用有限元还进一步地分析了不同硬度下橡胶支座的变形特点,从而确定了橡胶在不同硬度下的力学常数C 1和C 2,这对橡胶件的力学特性分析和设计具有更广泛的指导意义.1 橡胶材料的本构关系1.1 橡胶弹性理论橡胶材料在较短时间内及恒定的环境温度下通常被处理为各向同性不可压缩材料,其应变能密度函数W 是变形张量不变量I 1、I 2、I 3的函数[5],即W =W (I 1,I 2,I 3),其中,I 1=λ21+λ22+λ23, I 2=λ21λ22+λ22λ23+λ21λ23, I 3=λ21λ22λ23(1)式中,λ1,λ2,λ3是3个主伸长比.根据橡胶的不可压缩性,有收稿日期:2000211212作者简介:郑明军(1971—),男,河南温县人,硕士生.em ail :zmj -l @ 第25卷第1期2001年2月 北 方 交 通 大 学 学 报JOURNAL OF NORTHERN J IAO TON G UN IV ERSIT Y Vol.25No.1Feb.2001I 3=λ21λ22λ23=1(2)从而W 可以用变形张量不变量的级数形式表示,该式由Rivlin 所推导[5]W =∑∞i ,j =0C ij (I i -3)i (I j -3)j (3)式中,C ij 是材料常数. 一般广泛采用的是Mooney 2Rivlin 模型,即W =C 1(I 1-3)+C 2(I 2-3)(4)该模型能很好地描述橡胶变形在150%内的特性[6].由K irchoff 应力张量t ij 和Green 应变量γij 间的关系得到t ij =5W 5I 15I 1γij +5W 5I 25I 2γij +5W 5I 35I 3γij (5) 利用式(1)和式(2)得出主应力t i 和主伸长比λi 之间关系为t i =2λ2i 5W 5I 1-1λ2i 5W 5I 2+P ,其中,P 为任意流体静压力.各式相减消去P ,得到3个主应力的差值,即t 1-t 2=2(λ21-λ22)5W 5I 1+λ235W 5I 2t 2-t 3=2(λ21-λ23)5W 5I 1+λ215W 5I 2t 3-t 1=2(λ23-λ21)5W 5I 1+λ225W 5I 211.2 C 1和C 2的实验确定方法[7]对于单向拉伸或压缩,有t 2=t 3=0,则λ22=λ23=λ-11.因此t 1=2λ1-1λ215W 5I 1+1λ15W 5I 2(6)考虑方程(4),可见5W I 1=C 1, 5W I 2=C 2(7)把式(7)代入式(6)得t 12(λ1-λ-21)=C 1+1λ1C 2(8)式(8)是单向拉伸或压缩试验确定橡胶材料常数C 1和C 2的基本公式.得到C 1和C 2的方法是根据试验测试出不同拉伸比λ1下的应力值t 1,然后以1λ1为横坐标,以t 12(λ1-λ-21)为纵坐标,把试验点描述在相应的坐标系中,并把这些试验点回归成一条直线,C 1为这条直线的截距,C 2为这条直线的斜率.1.3 橡胶材料的硬度与C 1和C 2的关系对于橡胶材料,其弹性模量E 0与剪切模量G 有下述关系G =E 02(1+μ),由橡胶的不可压缩性得泊松比μ=015,从而E 0=3G.G 或E 0与材料常数的关系为G =2(C 1+C 2), E 0=6C 11+C 2C 1(9)文献[2,3,8]给出了橡胶硬度H r (IRHD 硬度)与弹性模量E 0的试验数据,经拟合得77第1期 郑明军等:压缩状态下橡胶件大变形有限元分析log E 0=0.0198H r -0.5432(10)橡胶硬度很容易测得,根据式(9)和式(10),可见在已知橡胶硬度下,其力学常数C 1与C 2之和取决于H r .2 橡胶件大变形有限元分析2.1 橡胶柱的大变形分析一硬度为60(IRHD 硬度)的橡胶圆柱,受轴向压缩载荷,通过两块刚性的金属平板施加于橡胶上.橡胶圆柱及其所受载荷均为轴对称,故取一过轴线的剖面进行有限元建模(见图1),计算软件为Ansys5.6的轴对称4节点橡胶单元.有限元分析中所需常数C 1和C 2一般由试验确定,测试C 1和C 2需要专门加工试样,但这仅在橡胶组件可用的时候,或者橡胶老化导致材料性能发生变化等情况下,因此这一方法显得不切实际.在本研究中,在给定C 2/C 1不同比值的条件下,采用1.3节的方法,由有限元计算出不同C 2/C 1条件下的载荷—变形曲线,与橡胶柱压缩实际试验的载荷—变形曲线相比,确定合适的C 2/C 1值.分别取C 1为0.735、0.700、01490,相应的C 2值分别为0.035、0.245,即C 2/C 1值为0、0.05、0.5,受压橡胶柱载荷—变形计算结果与试验结果见图2.由图2可见,变形量小于5mm 时,C 2与C 1之比对计算结果影响很小;变形量大于5mm 时,对于C 2/C 1=0,计算结果与Rivlin [2]分析结果一致,对于C 2/C 1=0.5,曲线上移,对于C 2/C 1=0.05时,有限元计算结果与试验吻合最好.图1 受轴向载荷橡胶圆柱有限元模型图2 橡胶圆柱的载荷—变形曲线2.2 橡胶支座的大变形分析一受轴向压缩载荷作用下受剪的橡胶支座,其硬度与前述橡胶柱相同,在顶面钢板加载[3].采用轴对称条件,橡胶支座的有限元分析模型见图3,使用软件和单元类型与橡胶柱相同,使用2.1中的C 1和C 2值进行计算,所得载荷—变形结果见图4,将实测载荷—变形曲线绘于图4中.可见在C 2/C 1=0.05时,有限元计算值与实测值最为吻合,这表明由受压圆柱分析后得出的材料常数C 1和C 2同样适用于同硬度橡胶组件的力学特性分析.图3 橡胶支座有限元模型图4 硬度60的橡胶支座载荷—变形曲线87北 方 交 通 大 学 学 报 第25卷2.3 不同硬度下橡胶材料常数C 1和C 2的确定对于该橡胶支座,文献[4]给出了不同橡胶硬度下支座的载荷—变形曲线(见图5).利用前述分析方法和有限元建模,并与实测值进行比较确定不同硬度下材料常数C 1和C 2的最佳取值.由图5计算结果与实测结果的比较可见:当橡胶硬度分别为40、60、70时,C 2/C 1在分别取0.1、0.05、0.02下,计算值与实测值较吻合.根据分析结果,绘制了C 1、C 2和C 2/C 1随H r 的变化曲线(见图6),这表明对于不同硬度的橡胶,C 2/C 1的值也不相同,表现为硬度提高,比值下降.图5 不同硬度下橡胶支座的载荷—变形曲线图6 不同硬度下的橡胶力学常数曲线3 结论在橡胶以压缩状态为主的条件下,橡胶材料力学常数C 1和C 2之和由橡胶硬度决定,且随硬度的增大而增大;在已知橡胶硬度及其载荷—变形曲线时,采用有限元分析可得到可靠的橡胶力学常数C 1和C 2;不同硬度的橡胶材料,其C 2与C 1的比值不同,C 2/C 1随硬度的增加而下降.参考文献:[1]杨晓翔.非线性橡胶材料的有限单元法[M ].北京:石油出版社,1999.[2]Lee B S ,Rivin E I.Finite Element Analysis of Load 2Deflection and Characteristics of Com pressed Rubber Components for Vi 2bration Control Devices[J ].Journal of Mechanical Design ,1996,118:328-335.[3][英]弗雷克利K ,佩恩P K.橡胶在工程中应用的理论与实践[M ].杜承泽,唐宝华,罗东山,等译.北京:化学工业出版社,1985.[4]PAUL STRA 公司.橡胶支座产品介绍[Z].法:PAUL STRA 公司,1998.[5]于建华.魏泳涛.不可压缩超弹性材料的有限元应力分析[J ].西安交通大学学报,1998,33(1):41-45.[6][英]特雷劳尔L R G.橡胶弹性物理力学[M ].王梦蛟,王培国,薛广智译.北京:化学工业出版社,1982.[7]李洪升,张小朋,杨全生.橡胶大变形力学常数测试研究[J ].大连理工大学报,1989,29(6):629-634.[8]戚震华,方永明,张定贤.橡胶弹簧非线性刚度的有限元解[J ].上海力学,1994,15(4):33-41.97第1期 郑明军等:压缩状态下橡胶件大变形有限元分析。

大变形问题有限元分析

大变形问题有限元分析

解线性方程组
通过求解由刚度矩阵构成的线性方程 组,得到离散解。
后处理
对离散解进行后处理,如误差估计、 收敛性分析等。
04
大变形问题的有限
元分析
有限元模型的建立
确定问题类型
选择单元类型
根据实际问题,确定是弹性问题、塑性问 题还是流体动力学问题等。
根据问题类型和求解精度要求,选择合适 的单元类型,如四边形单元、六面体单元 等。
在大变形问题中,由于物体的位移和变形较大,传统的有限 元分析方法可能无法准确描述物体的变形行为,因此需要采 用更高级的有限元分析方法。
研究意义
大变形问题在工程实践中具有广泛的应用,如桥梁、建筑 、航空航天等领域的结构分析。因此,研究大变形问题的 有限元分析方法具有重要的实际意义。
通过研究大变形问题的有限元分析方法,可以更好地了解 物体的变形行为,提高工程结构的可靠性和安全性。此外 ,该研究还可以为其他复杂工程问题的有限元分析提供理 论支持和方法指导。
求解方程组
利用选定的求解方法,求解建立的方程组, 得到各节点的数值解。
有限元分析的步骤和流程
对计算结果进行可视化、分析和 解释等。
建立方程组、选择求解方法和求 解方程组等。
建立几何模型、划分网格和离散 化处理等。
前处理
求解过程
后处理
05
有限元分析的实例
实例一:简单大变形问题分析
模型描述
考虑一个简单的弹性体,在受到外力作用时 发生的变形。
建立几何模型
划分网格
根据实际问题,建立相应的几何模型,包 括形状、尺寸和边界条件等。
将几何模型划分为有限个小的单元,每个 单元由节点和边组成。
有限元模型的求解
离散化处理

abaqus中橡胶大变形问题

abaqus中橡胶大变形问题

Abaqus中橡胶大变形问题橡胶材料在工程中广泛应用,其特性之一就是其在受力时会产生大变形。

在工程实践中,需要对橡胶材料的大变形行为进行准确的预测和仿真,以便设计出更加可靠和安全的产品。

而Abaqus作为一款强大的有限元分析软件,可以帮助工程师们对橡胶材料的大变形问题进行深入研究和分析。

在Abaqus中,对橡胶材料的大变形问题进行仿真和分析通常需要考虑以下几个方面的内容:橡胶材料的本构模型、边界条件的设定、大变形时的网格变形和接触问题等。

在本文中,我将针对这些内容展开深入的讨论和分析,并结合个人的经验和理解,希望能为你带来有价值的信息和见解。

1. 橡胶材料的本构模型橡胶材料的大变形行为是非线性的,因此在Abaqus中对其进行仿真时,需要使用适当的本构模型来描述其力学行为。

常见的橡胶材料本构模型包括各向同性模型、各向异性模型、超弹性模型等。

在选择本构模型时,需要考虑橡胶材料的实际性能和实验数据,以及仿真的准确性和计算效率。

需要对本构模型的参数进行合理的设定和校准,以确保仿真结果的准确性和可靠性。

2. 边界条件的设定橡胶材料在实际工程中往往处于复杂的受力和约束条件下。

在Abaqus中进行橡胶材料的大变形仿真时,需要对边界条件进行合理的设定。

这包括加载条件的设定、约束条件的设定以及边界条件的处理等。

合理的边界条件设置能够更好地模拟橡胶材料的受力和变形行为,从而得到准确的仿真结果。

3. 大变形时的网格变形和接触问题橡胶材料在受力过程中会产生较大的变形,这需要在Abaqus中进行合适的网格变形和接触处理。

在进行橡胶材料大变形仿真时,需要对网格进行合理的划分和调整,以适应材料的大变形,同时需要对接触问题进行有效的处理,保证仿真的准确性和稳定性。

总结回顾通过以上对Abaqus中橡胶材料大变形问题的讨论和分析,我们可以得出以下几点结论:在进行橡胶材料大变形仿真时,需要选择合适的本构模型,并对模型参数进行准确的设定和校准;在边界条件的设定上,需要考虑橡胶材料的受力和约束情况,以得到真实可靠的仿真结果;在进行大变形仿真时,需要合理处理网格变形和接触问题,以确保仿真的准确性和稳定性。

大变形问题有限元分析

大变形问题有限元分析

uK
,J
1 2
uK
,I
uK ,J
eIJ IJ
二者之间满足张
线性部分 非线性部分
量变换关系!
现时(Updated)Green应变增量:
* ij
1
ui
2 x j
u
j
xi
1 2
uk xi
uk x j
*eij *ij
IJ
xm X I
xn X J
* mn
非线性部分
2020/4/5
线性部分
4
kl
1 2
uk,l ul,k um,kum,l
注意:我们用下标的大小写表示坐标的大小写,对应于不同的构型。
大变形分析由于采用增量方法,需经常用到它们的增量形式。
2020/4/5
3
第3页/共26页
大变形问题的应变描述(3/4)
应变增量:
Green应变增量:
IJ
1 2
KJ
uK,J
uK ,I KI uK ,I
2020/4/5
9
第9页/共26页
大变形分析中的本构关系(3/5)
超弹性材料
假定材料具有单位质量的应变能函数,再根据能量原理来定义本构
关系,这类材料称为超弹性材料。
Case-1 W W KL
例如
W
1 20
IJ AIJKL KL
(不限于这种形式)
增量形式 …
SIJ
0
W KL IJ
初始构型时材料 的密度-常数
弹性材料:加载曲线与卸载曲线相同的材料。

本构关系有三种形式
(大变形分析中)
ij Aijkl kl
Aijkl 为常数
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

线性部分
非线性部分
二者之间满足张 量变换关系!
xm xn * mn X I X J
现时(Updated)Green应变增量:
* ij 1 ui u j 1 uk uk 2 x x i j 2 xi x j
IJ
1
引言
几何线性问题: 位移与应变成线性 (微分)关系;
几何非线性问题:位移与应变成非线性(微分意义上)关系。 物理现象:将位移(转动)和/或应变较大的问题统称为大变形 问题,有时称为有限变形问题。这类问题又分为大位移
(转动)小应变问题及大位移大应变问题两大类。 研究意义:和材料非线性问题一样重要。例如,平板的弯曲问题, 大挠度理论分析结果更符合实际情况;薄壳的屈曲,非线性理 论的预测值更好。又例如,对于橡皮型材料,大变形还必须考 虑本构关系的变化,这与纯粹的材料非线性又有区别。
*eij *ij
2019/1/15
线性部分
非线性部分是高阶小量
6
大变形问题的应力描述(1/2)
应力是借助于微元体来定义的,但在大变形分析中,必须注意 微元体所在的构型。 与应变类似,连续介质力学理论具有严格的应力定义和多 种不同的应力概念。这里也只介绍后面将要用到的几种。
Euler应力:
从当前构型中取出微元体,在其上定义的应力称为Euler应力,用 表示。Euler应力代表物体的真实应力。然而,当前构型是待求的未知构型, 因而,有必要通过已知构型上的微元体再对应力进行描述。
Kirchhoff应力:
通过初时构型上的微元体定义的应力称为Kirchhoff应力,用 S 表示; 通过现时构型的微元体定义的应力称为现时(Updated)Kirchhoff 应力, * 用 表示。 S
kl
1 uk ,l ul ,k um,k um,l 2
注意:我们用下标的大小写表示坐标的大小写,对应于不同的构型。 大变形分析由于采用增量方法,需经常用到它们的增量形式。
2019/1/15 4
大变形问题的应变描述(3/4)
应变增量: Green应变增量:
IJ 1 1 KJ uK , J uK , I KI uK , I uK , J uK , I uK , J 2 2 eIJ IJ
3
大变形问题的应变描述(2/4)
描述的出发点:物体的变形描述建立在确定的参考构型上。 Green应变张量:以初始构型为参考构型所定义的应变,数学 表示为
KL
1 u K , L u L , K uM , K uM , L 2
现时(Updated)Green应变张量:以现时构型为参考构 型所定义的应变,数学表示为
研究现状:大变形问题有限元分析的理论和方法存在不同学派间的 争鸣,尚未得到一个权威性的结论。随之并发的其它问题,如 解的稳定性、收敛性及收敛率等,都有待进一步深入研究。
2019/1/15 2
大变形问题的应变描述(1/4)
问题的特点:由于变形较大,使得不同时刻物体具有差别不能 忽略的不同构型,这是大变形问题分析的基本出发点。
2019/1/15 7
大变形问题的应力描述(2/2)
Kirchhoff、现时Kirchhoff及Euler应力(增量)间的关系:
*Sij ij * Sij
现时Kirchhoff应力增量
现时Kirchhoff应力 t t 时刻
Euler应力 t 时刻
特点:以现时构型为参考。
根据张量的坐标变换规则,它们之间还有以下关系
D
* N 1

y1 , y2 , y3 yi x1 , x2 , x3 x j
2019/1/15
8
大变形分析中的本构关系(1/5)
本构关系的客观性要求:需要选取合适的应力-应变共轭对描 述材料的本构关系。 弹性材料:加载曲线与卸载曲线相同的材料。

本构关系有三种形式
ij Aijkl kl
W ij ij
第三章 大变形问题的有限元分析
目的:以大变形问题为例,介绍几何非线性问题的有限元 方法。
特点:与线性有限元方法比较,几何关系不再是线性的。 内容:

2019/1/15
引言 大变形问题的应变描述 大变形分析中的应力描述及本构关系 大变形问题有限元方程的建立 大变形分析中的载荷处理 小结
线性部分
非线性部分是高阶小量 对于小变形情形
IJ * ij 1 ui u j 2 X i X j ij
现时(Updated)Green应变增量退化成:
* ij 1 ui u j 1 uk uk 2 x x j i 2 xi x j
1 xi x j Sij N Skl D X K X L
*
ij ij
1 D
* N 1
yi y j kl *Skl xk xl
D
N

x1 , x2 , x3 xi X 1 , X 2 , X 3 X J
*eij *ij
2019/1/15
线性部分
非线性部分
5
大变形问题的应变描述(4/4)
应变增量:(续)-对于大变形小应变情形 Green应变增量退化成:
IJ 1 1 KJ uK , J uK , I KI uK , I uK , J uK , I uK , J 2 2 eIJ IJ
xi
XI
yi
(a)
(b)
(c)
初始构型(0时刻)
现时构型(t 时刻)
当前构型( t t 时刻)
连续介质力学理论对物体经历大变形后的变形有严格的定义 和推导。这里不准备过多引入复杂的概念和符号,而是与小变形 理论对照,介绍进行大变形分析时必需的几个概念和术语。
2019/1/15
大变形问题的分析方法:增量法。
相关文档
最新文档