ANSYS 模拟超弹性球压缩
- 格式:pdf
- 大小:2.85 MB
- 文档页数:13
图1为一个承受内压的薄板,在其中心位置有一个小圆孔,相关的结构尺寸参考图1所示。
材料属性:弹性模量E=2e11Pa,泊松比为0.3。
拉伸载荷为:q=3000Pa。
平板的厚度为:t=0.01mm。
通过简单力学分析,该问题属于平面应力问题,又因为平板结构的对称性,所以只要分析其中的1/4即可,如图2所示。
图1 板的结构示意图图2 有限元分析见图一、前处理(1)定义工作文件名:Utility Menu>File>Change Jobname,弹出如图3所示的Change Jobname对话框,在Enter new Jobname后面的输入栏中输入Plate,并将New Log and error files复选框选为yes,单击OK。
图3 定义工作文件名对话框(2)定义工作标题:Utility Menu>File>Change Title,在出现的对话框中输入The Analysis of Plate Stress with small Circle,单击OK。
图4 定义工作标题对话框(3)重新显示:Utility Menu>Plot>Replot。
(4)关闭三角坐标符号:Utility Menu>PlotCtrls>Window Controls>Window options,弹出一个对话框,在Location of triad 后面的下拉式选择框中,选择Not Shown,单击OK。
(5)选择单元类型:Main Menu>Preprocessor>Element Type>Add/Edit/Delete,弹出Element Type对话框,单击Add按钮,又弹出如图5所示的Library of Element Types对话框,在选择框中分别选择Structural Solid和Quad 8node 82,单击OK,然后单击Close。
ANSYS中超弹性模型及其在橡胶工程中的应用张振秀,聂军,沈梅,辛振祥(青岛科技大学橡塑材料与工程教育部重点实验室,山东 青岛 266042)摘要:本文主要介绍了超弹性理论和有限元分析软件A N S YS 非线性超弹性模型,其中主要叙述了各种超弹性模型的本构方程(Polyn omid F orm 模型、Mooney -Rivlin 模型、Neo -Hookean 模型、Yeoh 模型、Arruda -Boyce 模型、Gen t 模型、Ogden 模型、Hyperf oam 模型、Bla tz -K o 模型),在不同情况下模型的选取、材料数据的获得、参数拟合及不同橡胶材料选用不同模型的应用实例。
关键词:有限元分析软件A N S YS ;橡胶;超弹性;模型选用中图分类号:TQ330.1 文献标识码:B 文章编号:1009-797X(2005)09-0001-05作者简介:张振秀(1975-),男,现为青岛科技大学高分子科学与工程学院材料加工工程专业在读硕士研究生,主要从事高分子制品及模具的计算机辅助设计及分析方面的工作。
收稿日期:2004-06-11为了适应市场经济激烈竞争的挑战,企业必须快速推出高质量、低成本的产品,计算机辅助工程分析(C A E )软件将成为关键的工具和手段,而有限元分析是其中应用最为广泛的技术。
橡胶制品被广泛应用于国民经济各行各业,如轮胎、密封圈、软管、垫圈、震动隔离器、护舷、轴承套等,其产品种类多,使用的材料性能和产品生产工艺非常复杂,产品性能质量难以确定。
应用传统技术及原有经验,产品的开发和结构设计时间长、成本高、费用大、与国际先进水平有所差距。
在现代橡胶制品生产企业中,使用有限元技术进行产品开发是一重要手段,对于缩短新品开发周期、提高产品质量及降低产品开发与制造成本具有重要意义。
1 超弹性理论[1,4,7]橡胶又称为弹性体,包括天然橡胶及合成橡胶,是无定形的高聚物。
橡胶是一种超弹性材料,具有良好的伸缩性和复原性。
目录什么是塑性 (1)路径相关性 (1)率相关性 (1)工程应力、应变与真实应力、应变 (1)什么是激活塑性 (2)塑性理论介绍 (2)屈服准则 (2)流动准则 (3)强化准则 (3)塑性选项 (5)怎样使用塑性 (6)ANSYS输入 (7)输出量 (7)程序使用中的一些基本原则 (8)加强收敛性的方法 (8)查看结果 (9)塑性分析实例(GUI方法) (9)塑性分析实例(命令流方法) (14)弹塑性分析在这一册中,我们将详细地介绍由于塑性变性引起的非线性问题--弹塑性分析,我们的介绍人为以下几个方面:∙什么是塑性∙塑性理论简介∙ANSYS程序中所用的性选项∙怎样使用塑性∙塑性分析练习题什么是塑性塑性是一种在某种给定载荷下,材料产生永久变形的材料特性,对大多的工程材料来说,当其应力低于比例极限时,应力一应变关系是线性的。
另外,大多数材料在其应力低于屈服点时,表现为弹性行为,也就是说,当移走载荷时,其应变也完全消失。
由于屈服点和比例极限相差很小,因此在ANSYS程序中,假定它们相同。
在应力一应变的曲线中,低于屈服点的叫作弹性部分,超过屈服点的叫作塑性部分,也叫作应变强化部分。
塑性分析中考虑了塑性区域的材料特性。
路径相关性:即然塑性是不可恢复的,那么这种问题的就与加载历史有关,这类非线性问题叫作与路径相关的或非保守的非线性。
路径相关性是指对一种给定的边界条件,可能有多个正确的解—内部的应力,应变分布—存在,为了得到真正正确的结果,我们必须按照系统真正经历的加载过程加载。
率相关性:塑性应变的大小可能是加载速度快慢的函数,如果塑性应变的大小与时间有关,这种塑性叫作率无关性塑性,相反,与应变率有关的性叫作率相关的塑性。
大多的材料都有某种程度上的率相关性,但在大多数静力分析所经历的应变率范围,两者的应力-应变曲线差别不大,所以在一般的分析中,我们变为是与率无关的。
工程应力,应变与真实的应力、应变:塑性材料的数据一般以拉伸的应力—应变曲线形式给出。
ANSYS命令流学习笔记17-超弹性材料分析及WB-ABAQUS 分析对⽐! ANSYS 命令流学习笔记17-超弹性材料分析及WB-ABAQUS 分析对⽐ !学习重点:⾮线性材料建⽴在线性材料的基础上,理解好线性才⾏,在概念上就能理解好⾮线性材料。
但是⾮线性的计算⼜是另外⼀个概念,先学习材料部分知识吧。
理解应⼒应变的张量形式、应变能函数、⾼度⾮线性下应变能函数形式。
!1、应变张量张量最初是⽤来表⽰弹性介质中各点应⼒状态的,在三维坐标下,应⼒和应变的状态可以⽤9个分量来表⽰,超弹性材料主要使⽤应变张量及应变张量不变量这两个概念。
任意⼀点的应变状态可由矩阵表⽰:??z zy zxyz y yx xz xy x εγγγεγγγε存在三个相互垂直的⽅向。
在这三个⽅向上没有⾓度偏转,只有轴向的应变,该正应变称为主应变,此三⽅向成为主⽅向。
此时,该点应⼒状态由矩阵表⽰:但是应变张量表达中,某⼀点的应变状态矩阵,和坐标⽅向的选取有着很⼤关系。
为了表达坐标⽆关的某点应变状态,定义应变张量不变量I 1、I 2、I 3 ,分别为应变张量的第⼀,第⼆和第三不变量。
由下式表⽰:取= 1/3*I 1,将应变张量可以分解为应变球张量和应变偏张量,分别对应应变的形状改变部分和体积改变部分。
+---=m mm m z zyzx yzmy yxxzxy m x ijεεεεεγγγεεγγγεεε000000?=321000000εεεεij m ε!2、应变能函数⼀维应变能函数:⼀维应变能密度函数:W 或U 函数形式能够确定的话,应⼒与应变之间的关系也就完全确定了,反之应变应⼒关系确定可以反推应变能密度函数。
可以认为应变能密度函数是材料本构关系的⼀种表达形式。
!3、应变能函数形式(1)延伸率、不变量、体积⽐在确定应变能函数形式之前,⾸先要确定应变能函数的变量。
⾸先定义延伸率λ:其中,E ε⼀般称为⼯程应变或名义应变。
(此外,⼀般说的⼯程应⼒,真实应⼒)。
压缩过程分析实例说明:有一输气管道,管道内径为1200mm,外径为1600mm,其横截面结构为一圆环,输送气体压力为0.5MPa。
求管壁的应力场分布。
管道材料弹性模量为260GPa,泊松比为0.26。
分析:管道沿长度方向尺寸较大,一般应远大于管道直径,该问题属于平面应变问题,可取其横截面的1/4建立有限元分析模型进行求解。
操作步骤第1步:建立工作文件名和工作标题(1)选择Utility Menu | File | Change Jobname命令,出现Change Jobname对话框。
(2)在Enter new Jobname文本框中输入工作文件名,并将New log and error files设置为Yes。
(3)选择Utility Menu | File | Change Title命令,出现Change Title对话框,在文本文档中输入文件名。
第2步:创建几何模型(1)选择Main Menu | Preprocessor | Modeling | Create |Areas | Circle| Partial Annulus命令。
(2)在Rad-1中输入0.6,在Theta-1中输入0,在Rad-2中输入0.8,在Theta-2中输入90。
(3)单击OK关闭对话框。
第3步:定义单元类型(1)选择Main Menu | Preprocessor | Element Type | Add/Edit/Delete命令,出现Element Types对话框。
(2)单击Add按钮,出现Library of Element Types对话框。
(3)在其中第一个列表框中选择Structural Solid,在第二个列表框中选择Quad 8node82,在下面的文本框中输入1。
(4)单击OK关闭Library of Element Types对话框。
(5)单击Elemenl Types对话框上的Options按钮,出现PLANE82 Element type options 对话框。
ANSYS 模拟超弹性球压缩
采用Ogden 三对材料常数模型
分析橡胶球的压缩,球的直径40mm。
将橡胶球压缩其直径的1/2。
几乎不可压缩的Ogden。
由二维轴对称PLANE182单元和刚-柔接触面组成的二维轴对称模型,该接触对考虑了板的厚度变化效应,接触对摩擦系数指定
为0.35。
1,选择结构分析类型,选择单元类型
Main Menu: Preprocessor →Element Type→Add/Edit/Delete→Add
→Solid - 4 Node 182 →OK→Option→K3:Axisymmetric →Ad d →Contact →刚性接触单元2D Targe169(显
示TARGE169) →Apply →
Contact →柔性接触单元 2 nd surf 171 显示
(CONTA171) →OK →[Close]
提示: 单元类型 1 采用PLANE182 单元,因为这是一个体积变形问
题, 所以采用缺省公式, 即B-Bar 方法。
采用几乎不可压缩超弹性材料属,所以不需要混合U-P 公式。
2. 定义材料参数--输入Ogden 模型参数
Main Menu →Preprocessor →Material Props →Material Mode
ls →Structural →Nonlinear →Elastic→Hyperelastic →O gden → 3 terms”→mu_1:6.3, a_1:1.3, mu_2:0.012, a_2:5.0, m u_3:-0.1, a_3:-2.0, d_1:2e-
4 →[OK] →New Material Model →Structural →Friction Coefficient →0.3
5 Exit
3.生成几何模型
生成特征点
Main Menu: Preprocessor →Modeling→Create→Key points →In
Active CS →依次输入 2 个点的坐标:input:1(0,0.02,0),
2(0.03,0.02,0) →OK生成一条横线
Main Menu: Preprocessor→Modeling→Create→Lines→Straight Line →连接 2 个特征点,1→2→OK
生成一个1/4 圆形面
Main Menu: Preprocessor →Modeling→Create→Areas
→Circle→Part ial Annulus(WPX=0,WPY=0,Rad1=0.02,Theta=0,
Rad2=0,Theta=90)→OK
4.网格划分
Main Menu: Preprocessor →Meshing→Mesh Tool→Element Attribut e →Area→Material num:1 (PLANE182)→Smart size (选6) →Mesh: Area, Shape: Free, Quad →Mesh→拾取圆弧形
面→Close
注意:由于要生成线刚-柔接触对,表示模具的横线要指定属性。
Main Menu: Preprocessor →Meshing→Mesh Tool →Element Attrib ute →Line→拾取横直线→Material num:2 →单元类型:
TARGE169 →Apply →Element Attribute →Line→拾取拾取圆表面线
→单元类型:CONTAC171 →OK
5.生成线刚-柔接触对
由于橡胶球和模具在接触时是无过盈配合,橡胶球外表面和模具的表面之间将构成线-线接触对。
在生成接触对的同时,ANSYS 程序将自动给接触对分配实常数号。
1).打开接触管理器
Main Menu→Preprocessor→Modeling→Create→Contact Pair
2).单击接触管理器中的工具条上的最上左边按钮,将弹出
Add Contact Pair 对话框。
3).指定接触目标表面。
模具看成是刚体,Target Surface:用默认的“Line”,在Target Type:选:Rigid w/Pilot,单
击Pick Target:弹出Select Line for Target 对话框,在图形输出窗口中单击那条横线,模具选定→OK→单击NEXT 进入下一
步→选Pick existing free keypoint →Pick Entity →选特征点 1 用作接触对的Pilot 点→OK →NEXT →将下一个画
面
4).弹出选中接触面的对话框,单击Pick Contact:弹
出Select Line for Target 对话框,在图形输出窗口中单击球的外环面将其选定→OK →NEXT →对接触对属性进行设置。
5).单击Material ID 下拉框中的“1”,指定接触材料属性为定义的一号材料。
并在Coefficient Of Friction文本框中输入“0.35”,指定摩擦系数为0.35。
单击按钮,来对接触问题的其它选项进行设置。
6).打开对摩擦选项设置的选项卡Optional setting,在对话框中
的Normal Penalty Stiffness 的框中输入“0.1”,指定接触刚度的罚系数为0.1。
在Friction 选项中的Stiffness matrix 下拉框中选:“Unsymmetric”选项,指定本实例的接触刚度为非对称矩阵。
其余的设置保持缺省,单击Create 按钮关闭对话框,完成对接触选项的设置。
6.模型施加约束
给球的左、下端施加对称约束
Main Menu: Solution →Define Loads →Apply→Structural→Displa
cement →Symmetry B.C. →选球的左、下端→OK
给模具 1 点施加位移-0.01m
Main Menu: Solution →Define Loads →Apply→Structural→Displa
cement →On Keypoints →选横线最左端点,
Keypoint 1 →OK →Uy: →VALUE 选-0.01 →OK
7.分析参数设定与计算
Main Menu: Solution → Analysis Type → New analysis → Sta tic → Sol’n Controls → Large deformation static → Number Of substeps =50; Max No substeps =100; Min No substeps =10 → OK
8.求解
Main Menu →Solution →-Solve- Current LS
9查看结果
面沿旋转轴旋转270 度,形成实体剖面图
Utility Menu →PlotCtrls →Style →Symmetry Expansion →2D Axi-Symmetic →3/4Expansion →OK
Main Menu →General Postproc →Plot Results →Contour Pl ot- Nodal Solu →Stress →Von Mises Stress →[OK]
Main Menu →General Postproc →Plot Results →Contour Pl ot- Element Solu →Stress →Hydrostatic Pressure
(HPRE)→[OK]
更多精彩关注公众号:工科小学生。