ANSYS 模拟超弹性球压缩
弹性力学ansys分析

图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中超弹性模型及其在橡胶工程中的应用

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]橡胶又称为弹性体,包括天然橡胶及合成橡胶,是无定形的高聚物。
橡胶是一种超弹性材料,具有良好的伸缩性和复原性。
ANSYS13.0 Workbench 结构非线性培训 超弹性

3. O.H. Yeoh, “Phenomenological Theory of Rubber Elasticity,” Comprehensive Polymer Science, ed. G. Allen, Elsevier, Oxford, 1996, Chapter 12.
– di 反比于体积模量. 默认地, 如曲线拟合(下一部分)中没引入体积试验数据, 则材料
假定为完全不可压缩的 (di=0).
N
iai
o
i 1
2
o
2 d1
... 体积容差
• 体积协调约束中的容差(vtol)可通过 Command Objects放松.
为接受后续的solc,,,vtol手动激活 Mixed u-P 是必要的
参考文献
一些关于橡胶机理的参考文献:
1. R.S. Rivlin, “Large Elastic Deformations,” Rheology: Theory & Applications - Vol. 1, ed. F.R. Eirich, Academic Press, Inc., New York, 1956, Chapter 10.
• 高弹体是一种聚合物, 具有如下性能
– 高弹体包括天然和合成橡胶, 它是非晶态的, 由 长的分子链组成
• 分子链高度扭转、卷曲, 且在未变形状态下取向任 意
• 在拉伸载荷作用下, 这些分子链部分变得平直、不 扭曲
• 去除载荷后, 这些分子链恢复最初的形态
ansys超弹性分析练习十五

Workshop Supplement
Advanced Structural Nonlinearities 6.0
–
Utility Menu > File > Read Input From… • 选择文件 “ANL_W11_Oring.inp” • 点击 [OK]
或命令输入方法:
–
/INPUT,
ANL_W11_Oring,
练习 8C: Ogden 模型
…超弹性O型环
3. 输入 Ogden 模型参数
使用 GUI 菜单方法: – Main Menu > Preprocessor > Material Props > Material Models …
Workshop Supplement
Advanced Structural Nonlinearities 6.0
inp
提示: 这将读入一个输入文件, 产生当前练习的几何体 、载荷和边界条件,得到一个由 PLANE182 单 元和刚-柔接触面组成的二维轴对称模型。
当前数据库将被删除, 工作名改为
“Exercise_11”。 图形和其它设置将与这些幻灯片一致。
September 30, 2001 Inventory #001492 W11-3
练习 8C: Ogden 模型
…超弹性O型环
继续进行该练习前, 也可检查模型的网格和边 界条件, 以便更熟悉它。 O 型环是一个轴对称 PLANE182 模型, 具有 几乎不可压缩的 Ogden 材料属性。 定义三个刚-柔接触对,在第一载荷步中, 最外 面的刚性目标面的拖拉距离为直径的1/3。 在第二载荷步中, 在外露面施加一个 20 单位 的压力。
刚性目标面
08 ANSYS13.0 Workbench 结构非线性培训 超弹性汇总

... 模型分类
其他模型 (Mooney-Rivlin, Yeoh 和 Neo-Hookean) 都是多项式模型的缩减形式。
• 它们的应力-应变关系是高度非线性的
• 通常, 拉伸状态下, 材料先软化再硬化,而压缩时材料急剧硬化. F
Tension
u
Compression
B.超弹性背景
ANSYS中关于超弹性本构模型有一些关键假设
• 材料响应是各向同性、等温和弹性的 • 热膨胀是各向同性 • 变形完全可恢复(保守的)
• 材料是完全或几乎不可压缩的 • 更复杂的真实橡胶行为理想化
Example of Rubber ot, o-rings/seals
...高弹体背景
宏观上,橡胶行为呈现下列特征
• 能承受大弹性(可恢复)变形,任何地方都可达100-700%。正如前面提及的,这是由于交联分子链 拆开的原因.
• 由于分子链的拉直引起变形, 所以在外加应力作用下, 体积变化很小。 因此, 高弹体几乎不 可压缩.
lp 和 I p
因而, 使用
可以将W分为偏差项和体积项。
... 应力和应变的计算
通过应变能密度函数,可计算应力应变.
• •
需基要于采W,用确应定力第-二应P变io共la轭-测Ki量rchoffS应ij力(和ddGEWreijen-Lagrange 应变) :
• 注意ANSYS结果以真实应力和应变输出。超弹性曲线拟合(稍后描述) 要求工程应力和工程应变.
... 多项式
多项式形式 基W于第一i和Nj第1二c应ij 变I不1 变量3,i 它I是2 如下3形式j 的现kN象1学d模1k型 J 1 2k
ansys超弹性分析练习十三资料

模型描述 • • 带厚度的二维平面应力 PLANE182 单元 刚-柔接触对
September 30, 2001 Inventory #001492 W10-2
练习 13: 超弹性键盘
…超弹性模型
1. 读入输入文件“ASNL_W08A_Keyboard.inp”
使用 GUI 菜单方法:
Workshop Supplement
或命令输入方法:
– –
/PREP7 ETLIST
提示:
因为预料该问题是弯曲占主要优势,将使用增强 应变公式。(如果时间允许,用B-Bar再运行该问 题并对比反作用力). 而且,前面所述,假定了具 有厚度的二维平面应力。
September 30, 2001 Inventory #001492 W10-5
Advanced Structural Nonlinearities 6.0
–
Utility Menu > File > Read Input From… • 选择文件 “ASNL_W08A_Keyboard.inp” • 点击 [OK]
或命令输入方法:
–
/INPUT,ASNL_W08A_Keyboard,inp
Advanced Structural Nonlinearities 6.0
–
Main Menu > Preprocessor > Element Type > Add/Edit/Delete … • 选择 “Type 1 PLANE182” • 点击 [Options] • 验证单元选项, 然后点击[OK] • 选择 [Close]
1 1 2 3 3 1
1 1 1
或命令输入菜单:
ansys超弹性分析练习十三

Advanced Structural Nonlinearities 6.0
–
Main Menu > Preprocessor > Material Props > Material Models … • 选择“Structural > Nonlinear > Elastic > Hyperelastic > Ogden > 2 terms” • “mu_1” 输入“80.194”, “a_1”输入“2”, “a_2” 输入“-2”, 其它空白。 • 点击 [OK] • 选择“Material > Exit” TB, HYPER, TBDATA, 1, 1, 1, 2, OGDEN 80.194, 2, 0, -2
提示: 记录最大 von Mises 应力(约149),画出其它 感兴趣的量, 如静水压力。
September 30, 2001 Inventory #001492 W10-9
练习 13: 超弹性键盘
…超弹性模型
7. 图示等效应变模型
使用GUI菜单方法:
Workshop Supplement
Advanced Structural Nonlinearities 6.0
–
Main Menu > General Postproc > Plot Results > Contour Plot > Element Solu • • 选择 “Elastic Strain > von Mises elastic strain” 点击 [OK] 图示等效应变
和命令输入方法:
–
PLESOL,EPEL,EQV
或命令输入方法:
– – – –
10 ANSYS13.0 Workbench 结构非线性培训 作业 超弹曲线拟合

Workbench Mechanical – Structural Nonlinearities
…作业 6A – 超弹曲线拟合
Workshop Supplement
4. 对分别为excel 文件 ‘tension-eb.xls‟ 和 ‘tension-pt.xls‟的二轴试验数 据和剪力试验数据重复步骤 3的过程
…作业 6A – 超弹曲线拟合
Workshop Supplement
15. 选择如图的顶点,并插入用户定义结果对x向应力应变进行后处理
Vertex for user defined post processing
– x 向应力表示为 “SX” – x 向应变表示为 “EPELX”
– 用户定义结果语法参见 PRNSOL 帮助命令
Workshop Supplement
WS6A-5
Workbench Mechanical – Structural Nonlinearities
…作业 6A – 超弹曲线拟合
• 点击 ‘Return to Project‟ • 双击 Model 打开 FE模型(Mechanical Session) (或 RMB=>Edit…)
– Biaxial Test Data – Shear Test Data
Workshop Supplement
WS6A-9
Workbench Mechanical – Structural Nonlinearities
…作业 6A – 超弹曲线拟合
2. 熟悉属性表和图
2a Properties 对话框对超弹试验数据应包含三个额外组件 2b 属性表描绘的是2a中高亮显示的组件数据. 2c Engineering Data 制图是 2b中列表数据的图形显示.
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
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]
更多精彩关注公众号:工科小学生。