最新大变形问题的有限元分析
有限元实验报告
有限元实验报告一、实验目的本实验旨在通过有限元方法对一个复杂的工程问题进行数值模拟和分析,从而验证理论模型的正确性,优化设计方案,提高设计效率。
二、实验原理有限元方法是一种广泛应用于工程领域中的数值分析方法。
它通过将连续的求解域离散化为由有限个单元组成的集合,从而将复杂的偏微分方程转化为一系列线性方程组进行求解。
本实验将采用有限元方法对一个具体的工程问题进行数值模拟和分析。
三、实验步骤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
分 片 近 似位
移 函 数
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中,对橡胶材料的大变形问题进行仿真和分析通常需要考虑以下几个方面的内容:橡胶材料的本构模型、边界条件的设定、大变形时的网格变形和接触问题等。
在本文中,我将针对这些内容展开深入的讨论和分析,并结合个人的经验和理解,希望能为你带来有价值的信息和见解。
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 为常数
[硕士]强夯加固地基大变形有限元分析初探_pdf
tamped pits
Based O/I the minor deformation
in classical mechanics
ቤተ መጻሕፍቲ ባይዱ
theory mat geometrical nonlinearity could not simulated the∞山d rigid body击印la∞m就呜
di曲l geometrical nonlinearity finite dement method is adopted in this paper to analyze the
在西安理工大学项士学位论文扰动荷载作用下随着塑性指数的增大扰动所诱致的孔隙水压力也相应增大由同一塑性指数重塑土样和原状土样的对比分析可知对于原状土而言由于其结构性的存在在冲击扰动荷载的初期扰动所诱致的孔隙水压力很小而后随着其结构破坏孔隙水压力有较大地增长当冲击次数达一定值后孔压增长才趋于平缓并与重塑土样有类似的变化趋势
actual由m锄血∞mp薯衄plastic血曲硼响蛀on Conclusions could be referenced in later similar projects. 2.ha the
would be∞∞e缸o being in soil
of蛐e血耐k top,on body.Rigid body di印la∞m咖on the point of compaction is the main托ason that forms the
1.Based On LS-DYNA3D∞nw鸸fla-ee-dimensional solid model is constructed usi|唱ela3to
ammih肌relation plastic
is弘DI鹏d Nonl/ncar theory
大变形问题的有限元分析概要课件
未来研究可以进一步优化有限元算法,提高计算效率 和精度,同时拓展有限元分析在大变形问题中的应用
范围。
有限元分析在大变形问题中具有重要应用价值 ,能够有效地模拟和分析复杂结构的变形和应 力散布。
有限元分析方法具有灵活性和通用性,可以广泛 应用于各种领域,如机械、航空、土木等。
研究展望
深入研究大变形问题的非 线性行为和复杂边界条件 ,提高有限元模型的精度 和可靠性。
随着科技的发展,大变形问题在许多工程领域中变得越来越重要,如航空航天、汽 车、建筑等。
有限元分析在大变形问题的研究中扮演着重要的角色,因为它能够模拟和分析复杂 的结构和材料行为。
目的和意义
目的
本课件旨在介绍有限元分析在大变形 问题中的应用,包括基本原理、方法 和实例。
意义
通过本课件的学习,学生可以了解有 限元分析在大变形问题中的重要性和 应用,为进一步研究和解决实际问题 打下基础。
实例三:流体大变形问题
总结词
材料属性
流固耦合问题,如流体对管道的冲刷和挤 压
流体密度、粘度、弹性模量等
变形特点
有限元分析
流体流动导致结构大变形,流固相互作用
采用流体动力学软件如ANSYS Fluent,考 虑流固耦合效应和流体对结构的冲刷
06
结论与展望
结论总结
通过对大变形问题的有限元分析,可以深入了 解材料的力学行为和变形机制,为工程设计和
有限元方法将连续的求解域离散化为一 系列的有限元,每个有限元都是一个简 单的、规则的几何形状,如四边形、六
面体等。
通过将问题离散化,有限元方法可以将 复杂的数学问题转化为求解一系列线性 方程组的问题,从而大大简化了理
建立模型并划分网格,将连 续的求解域离散化为有限个 简单元。
大变形问题的有限元分析课件
05
大变形问题的实例分析
总结词
弹性大变形问题是指物体在外力作用下产生的形变,形变程度超过弹性极限,导致材料内部应力散布产生变化的问题。
建立大变形问题的有限元模型需要选择合适的单元类型和材料模型。单元类型应能够准确地描述物体的变形行为,而材料模型则应能够反应物体的物理性质和行为。
在建立有限元模型时,还需要考虑边界条件、载荷和束缚条件等因素,以确保模型的准确性和有效性。
收敛性是评估有限元方法求解大变形问题的有效性的关键因素之一。收敛性分析涉及检查随着离散化程度的增加,解的近似值是否趋近于真实解。
方法是一种数值分析方法,通过将复杂的连续问题离散化为有限个简单的问题,从而求解复杂的数学模型。
它广泛应用于工程领域,如结构分析、流体动力学、电磁场等领域。
有限元方法的基本思想是将连续的求解区域离散化为有限个小的单元,并在每个单元上定义节点,通过节点之间的相互关系建立方程组,求解该方程组得到原问题的近似解。
总结词
大变形问题通常涉及到物体在外力或内力作用下的变形,这种变形可能是弹性变形、塑性变形、粘性变形等。根据变形的性质和程度,大变形问题可以分为不同的类型,每种类型都有其特定的数学模型和求解方法。
详细描述
大变形问题的基本方程包括运动方程、几何方程和本构方程。这些方程描述了物体的运动状态、变形几何和材料属性之间的关系。
详细描述
弹性大变形问题通常涉及到复杂的应力散布和应变状态,需要考虑材料的弹塑性行为和应变局部化现象。有限元分析在解决这类问题中发挥了重要作用,可以模拟材料的弹塑性行为和应力应变散布,为工程设计和优化提供根据。
大变形问题的有限元分析
目的:以大变形问题为例, 目的:以大变形问题为例,介绍几何非线性问题的有限元 方法。 方法。 特点:与线性有限元方法比较,几何关系不再是线性的。 特点:与线性有限元方法比较,几何关系不再是线性的。 内容: 内容: 引言 大变形问题的应变描述 大变形分析中的应力描述及本构关系 大变形问题有限元方程的建立 大变形分析中的载荷处理 小结
11
2011-10-19
现时构型时材料的密 随变形变化。 度-随变形变化。
大变形分析中的本构关系(4/5) 大变形分析中的本构关系(4/5)
次弹性材料 若应力率与变形率之间成线性变化规律,这类材料称为次弹性材料。但 应力率与变形率之间成线性变化规律,这类材料称为次弹性材料。 之间成线性变化规律 本构关系描述时要求“ 为与刚体转动无关的客观时间导数。 本构关系描述时要求“率”为与刚体转动无关的客观时间导数。
本构关系的客观性要求:需要选取合适的应力- 本构关系的客观性要求:需要选取合适的应力-应变共轭对 描述材料的本构关系。 描述材料的本构关系。 弹性材料:加载曲线与卸载曲线相同的材料。 弹性材料:加载曲线与卸载曲线相同的材料。
,
本构关系有三种形式
σ ij = Aijkl ε kl
∂W σ ij = ∂ε ij
初始构型时材料 的密度- 的密度-常数 Case-2
* W = W ( ε ij )
* ∂ ρ *W ( ε kl ) * Sij = * ∂ε ij
≜ AIJKLε KL
相 比 较
坐标变换
一阶近似
总之, 总之,对于一般的大变形 问题, 问题,在连续介质力学中 不能简化! 不能简化! 常用超弹性来表征材料的 本构关系。 本构关系。
有限元分析的力学基础
应用场景:流体 动力学分析广泛 应用于航空航天、 汽车、船舶、能 源等领域如飞机 机翼的气动性能 分析、汽车发动 机的流体动力学 分析等。
优势:有限元分 析能够处理复杂 的几何形状和边 界条件提供高精 度和可靠的分析 结果有助于优化 设计和改进产品 性能。
未来发展:随着 计算技术和数值 方法的不断进步 有限元分析在流 体动力学分析中 的应用将更加广 泛和深入有望在 解决复杂流体动 力学问题方面发 挥更大的作用。
特点:适用于大规模复杂问题的求解但需要设置合适的初值和解的精度要求。
有限元分析的精度与收敛性
精度:有限元分析的精度取决于网格划分的大小和形状以及所选择的近似函数。 收斂性:有限元分析的收敛性是指随着网格的细化解的近似值将逐渐接近真实解。 收敛速度:收敛速度取决于所选择的有限元类型和边界条件。 误差估计:通过误差估计可以确定所需的网格细化程度以确保解的精度。
弹性力学的 应用实例
塑性力学基础
定义:塑性力学是研究材料在达到屈服点后发生不可逆变形时行为规律的学科。 特点:塑性变形过程中外力的大小和方向可以发生变化而材料的内部结构保持不变。 塑性力学的基本方程:包括应力-应变关系、屈服准则、流动法则等。 应用:塑性力学在工程领域中广泛应用于金属成型、压力容器设计等领域。
局限性:塑性力 学模型忽略了材 料在塑性变形过 程中的微观结构 和相变行为因此 对于某些特定材 料或极端条件下 的应用可能存在 局限性。
流体动力学模型
简介:流体动力 学模型是有限元 分析中用于描述 流体运动的数学 模型包括流体压 力、速度、密度
等参数。
方程形式:流体 动力学模型通常 由一组偏微分方 程表示如NvierSkes方程描述了 流体的运动规律。
单元分析: 对每个单元 进行力学分 析包括内力、 外力、位移 等
L形管件弯曲成形有限元分析
L形管件弯曲成形有限元分析1. 引言- 描述L形管件在工业领域中的应用- 引入有限元分析作为一种计算方法- 概述本文将要讨论的内容和结论2. 有限元建模- 对L形管件进行三维建模- 定义有限元网格和边界条件- 描述管件材料和物理特性3. 弯曲成形模拟- 使用有限元分析模拟管件弯曲成形过程- 对比模拟结果与实验数据- 分析管件应力分布和形变情况4. 影响因素分析- 将影响管件弯曲成形的因素分为几类,如管件材料、管件壁厚等- 对每个影响因素进行有限元分析- 分析每个因素的影响程度和相互作用5. 结论和展望- 总结本文的研究成果和分析结论- 提出下一步的研究方向和拟解决的问题- 对应用前景和发展方向进行展望第1章:引言L形管件在工业领域中应用广泛,其具有良好的耐用性、耐腐蚀性和稳定性等特性,是各种管道系统中不可或缺的组成部分。
在实际应用中,L形管件的弯曲成形过程是至关重要的,因为这直接影响着管道的连接质量和稳定性。
传统的L形管件弯曲成形方法主要是通过人工操作机械或使用加热等手段进行成形,这种方法成本高且效率低,同时也容易影响到管件的质量。
为了提高管件成形的准确度和效率,有限元分析技术被广泛应用于管件成形过程的模拟和优化。
因此,本文旨在探讨如何通过有限元分析技术对L型管件弯曲成形过程进行模拟与分析,以更好地实现管件成形的精准度与质量。
本文将包括以下内容:首先介绍L形管件在工业中的应用背景和重要意义,然后介绍本文使用的有限元分析技术以及对L形管件三维建模的方法。
接着将介绍如何使用有限元分析模拟L形管件的弯曲成形过程,并对所得结果进行分析和讨论,最后探讨影响L形管件弯曲成形的因素,对因素进行分析和优化,得出结论并对未来发展进行展望。
第2章:有限元建模有限元分析技术是目前研究物理系统的一种最常用的数值分析方法。
在L形管件弯曲成形过程中,我们可以通过有限元分析技术对管件的形状变化,应力和形变等力学参数进行预测,从而优化管件的结构和形状设计。
边坡大变形弹塑性有限元分析
理状态, 滑坡体 分为四层, 第一层 为粘土 层( 滑 动带) 。 第二层为粉性土, 雨季常在地下水 位以下, 而第三层的 土性与第二层相同, 除了特大暴雨 外, 在通常情况下常 位于地下水位以上。第四 层土为粘性土, 主要分布在滑 坡的前部。这些土层的分布与物理力学性质见图 1 和 表 1, 从中可以看出土层的工程性质有明显的区别。与 通常情况 比较, 雨 季时第一 层土的凝聚 力和摩 擦角大 大下降, 所有土层的容重都有不同 程度的增加。滑坡的 这些地质条件十分有利于边坡的滑动。当下雨时, 地表 水可随着第二和三层土体中的孔隙和裂隙渗透到透水 性较差的 第一层和第 四层的边 界上, 结 果大量 的雨水 滞留在第二层和第三层中, 不易排 出, 造成滑坡体后部 荷重增加。而第一层粘土 由于被饱和其强度也降低, 这 些因素极易导致边坡的不稳定。
2. 3 变形 由 于本文 采用 了大 变形 弹塑 性模 型和 U pdat ed-
L agr ag ian 方法, 滑坡在荷载作用下的变形过程就 可以 进行模拟。图 3 展示了滑 坡初始和最终的变形形状。从 变形后的滑坡形态可以看出, 所有土层特别是第一层
图 2 滑坡体中屈服高斯点的分布 Fig. 2 Distribution of Gauss points
of yield in the landslide mass a—初始状态; b —最终状态
图 3 滑坡的初始和最终变形图形 Fig. 3 Diagram showing the initial and
f inal def ormation of landslide 土产生了较大的变形, 滑坡前部的 移动量比后部要大, 前部在 X 方向上的移动量达 2. 55m ( 见图 4) 。事实上, 在每一荷载作用下的变形也反映了滑坡在降雨过程中 的变形过程, 因为滑坡的净荷载主 要是来自于降雨。 2. 4 应力
悬臂梁大变形的向量式有限元分析
d s a e n s o t i d a c r i o i tr a o a o c n q i ae tma so a lm e t a d t e iplc me ti b ane c o d ngt n e ln d fr ea d e u v n s ff me e e n s n h n l l r
a c mp t t n p o e u e whih i r g a o u ai r c d r c sp o r mme ORTRAN.Th o o d by F e c mpu ain r s t e c n itn t tto e ul a o sse twi sr h t e t e r tc ou in.T e meho a s d t e f r t esmulto d a a y i fl g a t e e h h o eia s l to h t d c n beu e o p ro m h i l ain a n sso a e c n i v r n l r l
ee n meh d s s d o l me t to i u e t dic eie tu t r s nt p ril s se a d ik g ee n s mo g s r t sr cu e i o a t e y t ms n l a e l me t a n z c n
有限元分析—温度场与热变形问题专题
• 设单元只有三节点温度,jk为边界,将温度插值函数代 入前述的泛函,并求导得极值条件:
J e T T T T k[ ( ) ( )]dxdy Ti x Ti x y Ti y e T (T Ta ) T j jk 0
e e e e
• 两部分相加可得边界单元的温度刚阵:
1 ( H H ) T p 0
即
H T
J J e 0 m 1, 2,...n Tm e Tm
y
6 1 3 2 1 5
h11
h12 h1n T1 p1 T p h22 h2 n 2 2 hnn Tn pn
H T p
• 对于其他带热源的稳态温度场或三维温度场计算 其方法相似。
9-3 平面稳态温度场的有限元法
• 1、泛函与变分 • 函数 y=f(x) 求y 的极值,即求微分,由dy=0 可得。 • 泛函J=J [y(x)] 函数y(x)为自变量,J为函数y的函数,称J为y J 0 的泛函,求泛函的极值,即求变分, 由 可得。 • 例:平面上AB两点,连接AB的曲线很多,要求一条曲线使重 物靠自重由A沿此曲线滑到B所需的时间最短,即求最速下降 曲线。 • 显然,AB间直线路径最短,但重物运动的速度增长并不是最 大,即下滑的时间并非最短。
k T (Ta T ) n
• 1、三维瞬态热传导方程及边界条件
T T T T c (k x ) (k y ) (k z ) Q 0 在内 t x x y y z z T ( x, y, z, t ) T (1 , t ) T k (Ta T ) n 在1上 在 2上
大变形问题的有限元分析
1
引言
几何线性问题: 位移与应变成线性 (微分)关系;
几何非线性问题:位移与应变成非线性(微分意义上)关系。
物理现象:将位移(转动)和/或应变较大的问题统称为大变形 问题,有时称为有限变形问题。这类问题又分为大位移
(转动)小应变问题及大位移大应变问题两大类。
研究意义:和材料非线性问题一样重要。例如,平板的弯曲问题, 大挠度理论分析结果更符合实际情况;薄壳的屈曲,非线性理 论的预测值更好。又例如,对于橡皮型材料,大变形还必须考 虑本构关系的变化,这与纯粹的材料非线性又有区别。
Case-1 W ? W ??KL ?
例如
W
?
1
2? 0
? IJ AIJKL ?KL
(不限于这种形式)
增量形式 …
? ? SIJ
?
?0
?W ? KL ? ? IJ
初始构型时材料 的密度- 常数
一阶近似
? ? SIJ
?
?0
? 2W ? MN ?? IJ ?? KL
? KL
A ? IJKL KL
坐标变换
? ? ? ? ? ? ??K
S
t IJ
???? U ? ?
?F
SItJ
刚度矩阵
载荷向量 形式较复杂,因问题的类型而不同。
TL 法的求解步骤:
Step 1 :利用有限元方程求出 t ~ t ? ? t 间隔内的位移增量 ? UI ;
Step 2 :利用几何关系,计算 Green 应变增量 ? ? IJ ;
?
1 2
??
? KJ
? uK ,J
? uK,I ?
? KI ? uK , I
? uK ,J
?? ?
某门座起重机大车平衡梁局部起鼓变形的有限元分析
某门座起重机大车平衡梁局部起鼓变形的有限元分析发布时间:2022-09-20T05:24:44.172Z 来源:《工程建设标准化》2022年第10期作者:褚悦[导读] 介绍了某门座起重机行走机构大平衡梁结构出现的局部起鼓变形情况,应用有限元软件Ansys对平衡梁结构进行三维计算与分析褚悦中船第九设计研究院工程有限公司上海市杨浦区 200000摘要:介绍了某门座起重机行走机构大平衡梁结构出现的局部起鼓变形情况,应用有限元软件Ansys对平衡梁结构进行三维计算与分析,制定平衡梁结构的加强修护方案,并对优化后的方案进行有限元分析,论证加强结构的有效性。
关键词:门座式起重机;平衡梁;起鼓变形;有限元分析一、前言近年来,随着我国船舶业的蓬勃发展,门座起重机作为重要的装卸设备,在现代港口及船厂具有广泛的应用。
由于其使用频繁,且常年在露天或恶劣天气条件下工作,载荷工况复杂多变,极易引发变形及疲劳问题。
大车平衡梁作为起重机行走机构的重要部件,通过铰接的形式连接台车与主结构,起着支撑整机结构质量、传递轮压等重要作用,其结构与强度将直接影响门机的正常使用,因此对平衡梁结构变形的分析十分重要。
本文对某船厂一台门座起重机大车平衡梁结构起鼓变形维修进行过分析,通过结构受力分析,来判断起鼓变形原因,并根据计算分析结果制定相应的加强修护复方案,论证其对结构的有效性。
二、平衡梁受力情况及起鼓变形部位该门座式起重机额定起重量为75t,整机自重1335t,轨距12m,基距17m,海陆侧共布置4组大车行走机构,每组设大平衡梁一组,中平衡梁(I)1组,中平衡梁(II)3组,小平衡梁3组,车轮9只。
平衡梁钢结构采用Q345B焊接而成,如图一所示。
平衡梁所受垂直方向力主要由整机自重与工作载荷引起,水平方向力主要由工作时的惯性力及风载荷引起[1] 。
图1 行走机构示意图在日常检修中,发现在中平衡梁与小平衡梁连接铰点处发生起鼓变形,主要集中在中平衡梁腹板拐角处至重磅板焊缝上部,起鼓变形处圆弧的长度约为280mm,宽度约为200mm,起鼓最大处的高度约为46mm。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 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
大变形问题的分析方法:增量法。