ansys材料模型.doc
ANSYS邓肯张材料模型

ANSYS邓肯-张材料模型楼主给的在ANSYS上实现邓肯-张模型的方法很有用,但其中还有几点需要修正的,这也是楼上的兄弟们有疑问的原因。
我把楼主的代码运行了一下,然后对照作了修改,现在上传一下,有问题的兄弟可以仔细对照一下,在这里我对其中几个比较明显的问题说明一下:1.MP命令不能直接给单元加材料,这是对的。
在这里,楼主遗漏了一下命令:MPCHG,具体见下面的修改过的代码。
2.关于密度的问题。
这些要在宏中定义,每修改一种材料(即调用一次邓肯-张子程序)就要修改一次材料的密度,其他有关材料的问题可以类推。
3.关于施加重力的问题。
要在调用宏后,在同一个循环中重新定义一下重力。
以下是我修改过的楼主的代码,希望对兄弟们有所帮助。
!用APDL得到初步成果,贴于此供感兴趣的朋友参考,不当之处敬请指正,!欢迎加以完善。
!基本思路:!邓肯-张模型的关键点是材料的弹性模量随大小主应力差!及小主应力(围压)的变化而变化,用APDL实现之的基本思路是:!给每个单元定义一个材料号,分级施加荷载,在每个荷载步结束时提取出各!单元的大小主应力,据此计算出下个荷载步的弹性模量Et,修改各单元之MP,!用于下一步计算。
!以下是一个简单算例,copy出去可直接运行。
!!!常规三轴试验模拟!**********************************************************FINISH/CLEAR/TITLE,Numerical Simulation of three axes testing of soils/PREP7*dim,SUy,array,50!Settlement records*dim,MaxPs,array,120!Max history p1-p3*dim,MaxDs,array,120!Max history Ds!*dim,EEt,array,50!Et of elememt!!!Duncan-Chang Model!!!Symbols:c-粘滞力,Fai-内摩擦角,Sf-破坏强度(p1-p3)f,!Ds-应力水平,Pa-大气压,P3-围压!********************************************************************** *CREATE,Duncan-Chang!Creat Macro file*afun,deg!Unit of angle*set,Pa,1e5*set,P1,-ArrS3(i)!注意:岩土工程中应力为拉负压正*set,P3,-ArrS1(i)*if,P3,LT,0.1*Pa,thenP3=0.1*Pa!围压最小取值*endifSf0=2*(c0*cos(Fai)+P3*sin(Fai))/(1-sin(Fai))!Mohr-Coulomb破坏强度(p1-p3)fDs=(P1-P3)/Sf0!应力水平,*if,Ds,GT,0.95,thenDs=0.95!应力水平最大取值*endif!判断加卸荷,如果(P1-P3)小于历史最大值视为卸荷-再加荷过程*if,MaxPs(i),LT,P1-P3,thenEi=k0*Pa*(P3/Pa)**n0Et0=Ei*(1-Rf*Ds)**2!加荷情况的切线模量MaxPs(i)=P1-P3!保存历史最大应力*elseif,MaxPs(i),GE,P1-P3Et0=Kur*Pa*(P3/Pa)**n0!卸荷模量*endifmp,ex,i,Et0!修改材料i的Etmp,nuxy,i,Mu0mp,dens,I,1800!重新定义第I种材料的密度*END!***********************单元类型*********************************** et,1,42!平面四节点单元KEYOPT,1,3,2!平面应变!*********************以下定义材料初始模量************************* mp,ex,1,3.728e7!砂土的弹性模量mp,nuxy,1,0.33mp,dens,1,1800!****************************************************************** blc4,0,0,0.08,0.15!建立几何模型/PNUM,AREA,1/REPLOTaesize,all,0.01!网格划分mat,1amesh,allnsel,s,loc,y,0!边界条件d,all,Uy!底边界竖向约束nsel,s,loc,x,0d,all,Ux!左侧边界水平向约束nsel,all/replotfini/SOLUtime,0.01!施加围压ACEL,0,9.8,0sfl,all,pres,2e5!200kPasolve!********分级施加荷载,实现非线性计算,荷载增量10kPa,共50级******** *DO,ti,1,2!取出计算结果,修改弹性模量/POST1*get,SUy(ti),node,29,u,y!Settlement record of time ti ETABLE,EtabS1,S,1!取各单元第一主应力ETABLE,EtabS3,S,3!取各单元第三主应力*dim,ArrS1,array,120*dim,ArrS3,array,120*do,Num,1,120!Num为单元编号*get,ArrS1(Num),elem,Num,etab,EtabS1!将单元结果存入数组*get,ArrS3(Num),elem,Num,etab,EtabS3*enddo/PREP7c0=0Fai=35Rf=0.7k0=400n0=0.6Mu0=0.33Kur=326.7!修改砂土单元的Et,单元号1-120*do,i,1,120!各单元循环计算*use,Duncan-Chang,c0,Fai,Rf,k0,n0,Mu0,Kur!调用Duncan-Chang宏文件mpchg,i,i!!!!!!注意,这个命令是把第I种材料施加给第I个单元*enddo/SOLU!EEt(ti)=ET!保存第120单元之ETtime,tisfl,3,pres,2e5+1e4*ti!施加荷载,增量1e4ACEL,0,9.8,0!重新施加重力solve!对ti级荷载情况求解*ENDDO。
ANSYS超弹性、粘超弹性模拟

2 3
单轴拉伸与压缩实验
11 2 12
1 W 1 W 2 I1 1 I 2
1 2 12 -
正交双轴拉伸实验
1 W 2 W 2 22 12 I1 I 2 1 W 2 W 1 22 12 I1 I 2
13
© 2011 ANSYS, Inc. September 2, 2013 Release 14.0
粘弹性模型 静态
其中剪切松弛模量的Prony级数表达式为
n t G t G0 i exp i 1 i
其中, G0——t = 0时的松弛模量 G∞——t =∞时的松弛模量
September 2, 2013
Release 14.0
粘弹性模型 动态 滞后
• 滞后:试样在交变应力作用下,应变变化落后于应力变化的 现象
(t ) 0 sin wt
σ(t) ε(t)
0 δ σ ε (粘弹性) π
(t ) 0 sin(wt )
3π ω t t
2π
材料的变形过程是可逆的,无其它不可逆伴随,变形过程中 的熵变为零,此种材料成为超弹性材料。
2
W W 1 B2 B pI I1 I 2
式中I──单位变形张量 p──球张量 Ii──为变形张量B的不变量 W──应变能函数 基于假设:各向同性、不可压缩
W W I1 , I 2 , I3
9
© 2011 ANSYS, Inc.
September 2, 2013
Release 14.0
粘弹性模型
ANSYS结构分析-材料模型

ANSYS 结构分析材料模型
1 材料模型的分类
a. ANSYS 结构分析材料属性:
线性(Linear)、非线性(Nolinear)、密度(Density)、热膨胀(Thermal Expansion)、阻尼(Damping)、摩擦系数( Friction Coefficient)、特殊材料(Specialized Materials) 等七种,可通过材料属性菜单分别定义。
b. 材料模型:
线性、非线性及特殊材料三类,每类材料中又可分为多种材料类型,而每种材料类型则有不同的属性。
2 材料模型的定义及特点
材料模型及其属性均可通过GUI 方式输入。
线弹性材料可通过MP 命
令输入,而非线性及特殊材料则通过TB 命令定义,其属性则通过TBDATA 表输入。
表中前几项是常用的塑性材料模型,其后部分的材料模型有专用材料模型和可与前几项组合使用的材料模型。
表中屈服准则列中的Mises/Hill,指针对不同的单元分别采用Mises 屈服准则或Hill屈服准则,凡是可以考虑塑性的所有单元均可采用二者。
常用的单元
杆单元:LINK8、LINK10、LINK180
梁单元:BEAM3、BEAM4、BEAM188、BEAM189
管单元:PIPE16、PIPE20
2D 实体单元:PLANE82、PLANE183
3D 实体单元:SOLID65、SOLID92/95、SOLID191
壳单元:SHELL63、SHELL93、SHELL181
弹簧单元:COMBIN14、COMBIN39
质量单元:MASS21
矩阵单元:MATRIX27
表面效应单元:SURF154。
ansys粘弹性maxwell模型全参数总结材料

1.粘弹性:ANSYS中的粘弹性模型是Maxwell模型的通用积分形式,其松弛函数由Prony级数表示。
该模型功能全面,Maxwell、Kevin和标准线性实体都是其特殊形式,全面支持亚粘弹性和大应变超粘弹性。
大应变超粘弹性基于Simo建议的列式,粘弹性行为的定义分为超弹性和松弛两个部分,所有的ANSYS超弹性材料模型都可采用粘弹性选项(PRONY)。
2. 粘弹性是率相关行为, 材料特性可能与时间和温度都有关,粘弹性响应可看作由弹性和粘性部分组成。
–弹性部分是可恢复的, 且是瞬时的。
–粘性部分是不可恢复的, 且在整个时间范围内发生。
ANSYS 中能模拟线性粘弹性,这导致如下假设:
–应变率与瞬态应力成比例
–瞬态应变与瞬态应力也成比例
–限于小应变、小变形行为(NLGEOM,OFF)
–
C5=1
FICT TEMP可以从帮助文件里找到
注意密度。
ansys模型建立-材料设置-参数化过程

一.建立三维模型
8.标注如图尺寸
9.回到三维建模页面建立拉伸选项,operation选项该为去除材 料,生成特征
一.建立三维模型
10.建立阵列选项Greate→Pattern,选择阵列几何为圆孔内表面 点击Apply,在选择阵列轴为Z轴,生成如图特征。
一.建立三维模型
11.建立旋转特征,选择ZX平面并建立如图草图,尺寸标注如 图。
点击
三.参数化
建立变量之间的关系,此处使得小圆直径为大圆直径的1/6 如果含有更多变量可以建立更多变量关系,此处只建立大圆与
小圆直径的关系
三.参数化
2.双击Parameters,选择优化结果,如图在Table of Design poits
中添加不同变量的尺寸大小,点击
进行计算。
三.参数化
二.Mechanicals设置
6.施加载荷与约束(这里只添加简单的约束)如图,选择六个 圆孔面添加固定约束
选择凹槽面施加轴承载荷, 载荷设置如图
二.Mechanicals设置
7.设置结果选项,点击运算产生结果,如图。
三.参数化
1.返回DesignModeler,选择大圆的直径和孔的直径实现参数化。 如图点击尺寸前的小方框,并且命名。
ANSYS 培训举例
一.建立三维模型
1.点击XYplane
2.点击如图图标使得XY面正对屏幕
一.建立三维模型
3.选择草绘命令建立如图所示图形
一.建立三维模型
5.标注尺寸
一.建立三维模型
6.回到三维建模页面建立拉伸选项Extrude1,拉升长度为10mm
一.建立三维模型
7.建立孔特征,点击面选择工具,选择拉伸体的一个端面为建 立草图的平面,建立如图草图
压电材料ansys建模需要输入的属性

在定义完单元类型后,需要定义材料特性,压电模型需要输入的材料特性参数有:介电常数(又名电容率),压电矩阵弹性系数矩阵弹性基底材料及属性1、介电系数矩阵对于SOLID98 单元,可用MP 命令确定相对介电系数值作为PERX、PERY 和PERZ。
介电系数值分别代表介电系数矩阵的对角分量,即用MP 命令输入的介电系数总认为是常应变下的介电系数.2、压电矩阵可以用[ e ]形式-压电应力矩阵或用[ d ]形式—压电应变矩阵定义压电矩阵。
[ e ]矩阵通常与以刚度矩阵[ c ]形式表示各向异性弹性输入有关,而[ d ]矩阵和柔度矩阵相关.ANSYS 使用第一次界定温度下的弹性矩阵将压电应变矩阵[ d ]矩阵转换成压电应力矩阵[ e ]。
为确定这种变化所需要的弹性矩阵,可用TB,ANEL 命令.使用TB,PIEZ 和TBDATA 命令定义压电矩阵。
对于大多数已经公布的压电材料,压电矩阵数据都是按照x,y,z,yz,xz,xy 顺序的,基于IEEE 标准之上,而ANSYS 的输入数据是按照以上所示的x,y,z,xy,yz,xz 顺序。
也就是说,输入该参数是必须通过改变剪切项的行数据以转换到ANSYS 数据格式。
因此需将IEEE 常数61 62 63 [ e e e ]输入为ANSYS 的xy 行;将IEEE 常数41 42 43 [e e e]输入为ANSYS 的yz 行;将IEEE 常数51 52 53 [ e e e ]输入为ANSYS 的xz 行。
具体形式如矩阵(5.2)所示:3、弹性系数矩阵该矩阵为6 ×6对称矩阵(对于2-D 模型为4 ×4矩阵),它确定刚度系数([ c ]矩阵)或柔度系数([ s ]矩阵).弹性系数矩阵使用以下的数据表输入:使用TB,ANEL 和TB,DATA 命令定义系数矩阵[ c ]。
和压电矩阵的情况类似,大多数压电材料的[ c ]矩阵使用不同的参数顺序,需要通过如下方式改变剪切项的行和列数据以将IEEE矩阵转换到ANSYS 输入顺序。
(完整word版)ANSYS19.0帮助文件中材料参考MaterialReference

Material ReferenceTable of Contents1。
Introduction to Material Models1。
1。
Material Models for Displacement Applications1。
2. Material Models for Temperature Applications1.3. Material Models for Electromagnetic Applications1。
4。
Material Models for Coupled Applications1。
5. Material Parameters1。
6. How Material Properties Are Evaluated2. Material Model Element Support3。
Linear Material Properties3。
1. Defining Linear Material Properties3。
2. Stress-Strain Relationships3。
3。
Anisotropic Elasticity3.4. Damping3.5。
Thermal Expansion3.6。
Emissivity3。
7. Specific Heat3.8. Film Coefficients3。
9。
Temperature Dependency4。
Nonlinear Material Properties4.1. Understanding Material Data Tables4.2. Experimental Data4.3。
Porous Elasticity4。
3。
1。
Defining the Porous Elasticity Model4.4. Rate-Independent Plasticity4.4。
1。
Understanding the Plasticity Models4.4。
ANSYS材料模型

第七章材料模型ANSYS/LS-DYNA包括40多种材料模型,它们可以表示广泛的材料特性,可用材料如下所示。
本章后面将详细叙述材料模型和使用步骤。
对于每种材料模型的详细信息,请参看App endix B,Materi al ModelExampl es或《LS/DYNA Theore tical Manual》的第十六章(括号内将列出与每种模型相对应的L S-DYNA材料号)。
线弹性模型·各向同性(#1)·正交各向异性(#2)·各向异性(#2)·弹性流体(#1)非线弹性模型·Blatz-ko Rubber(#7)·Mooney-Rivlin Rubber(#27)·粘弹性(#6)非线性无弹性模型·双线性各向同性(#3)·与温度有关的双线性各向同性(#4)·横向各向异性弹塑性(#37)·横向各向异性FLD(#39)·随动双线性(#3)·随动塑性(#3)·3参数Bar lat(#36)·Barlat各向异性塑性(#33)·与应变率相关的幂函数塑性(#64)·应变率相关塑性(#19)·复合材料破坏(#22)·混凝土破坏(#72)·分段线性塑性(#24)·幂函数塑性(#18)压力相关塑性模型·弹-塑性流体动力学(#10)·地质帽盖材料模型(#25)泡沫模型·闭合多孔泡沫(#53)·粘性泡沫(#62)·低密度泡沫(#57)·可压缩泡沫(#63)·Honeyc omb(#26)需要状态方程的模型·Bamman塑性(#51)·Johnso n-Cook塑性(#15)·空材料(#9)·Zerill i-Armstr ong(#65)·Steinb erg(#11)离散单元模型·线弹性弹簧·普通非线性弹簧·非线性弹性弹簧·弹塑性弹簧·非弹性拉伸或仅压缩弹簧·麦克斯韦粘性弹簧·线粘性阻尼器·非线粘性阻尼器·索(缆)(#71)刚性体模型·刚体(#20)7.1定义显示动态材料模型用户可以采用ANSYS命令 MP,MPTEMP,MPDATA,TB, TBTEMP和 TBDATA以及ANSY S/LS-DYNA命令 EDMP来定义材料模型。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
! Failure strain
TBDATA,4,40.0
! C (strain rate parameter)
TBDATA,5,5.0
! P (strain rate parameter)
TBDATA,6,1
! LCID for true stress vs. true strain (see EDCURVE below)
TB,BISO
TBDATA,1, (屈服应力)
TBDATA,2, (切线模量)
例题参看B.2.7,Bilinear Isotropic Plasticity Example:Nickel Alloy。
B.2.7. Bilinear Isotropic Plasticity Example: Nickel Alloy
注--如果采用载荷曲线LCID1,则用TBDATA命令输入的屈服应力和切线模量将被忽略。另外,如果C和P设为0,则略去应变率影响。如果使用LCID2,用TBDATA命令输入的应变率参数C和P将被覆盖。只考虑真实应力和真实应变数据。在数据曲线一节中讲述了此种类型的例题。
注--例题参看B.2.16,Piecewise Linear Plasticity Example:High Carbon Steel。
应力应变特性只能在一个温度条件下给定。用MP命令输入弹性模量(Exx),密度(DENS)和泊松比(NUXY)。用TB,PLAW,,,,1和TBDATA命令中的1-6项输入屈服应力,切线斜率,硬化参数,应变率参数C和P以及失效应变:
如下所示,可以用TB,PLAW,,,,10和TBDATA命令中的1-5项定义其它参数。
MP,ex,1,100e9
! Pa
MP,nuxy,1,.36
! No units
MP,dens,1,4650
! kg/m3
TB,BKIN,1
TBDATA,1,70e6
! Yield stress (Pa)
TBDATA,2,112e6
! Tangent modulus (Pa)
7.2.3.6塑性随动模型
使用EDMP命令的同时,必须用MP命令定义刚体材料类型的杨氏模量(Ex),泊松比(NUXY)和密度(DENS)。必须指定实际的材料特性值,从而使程序能计算接触表面的刚度。基于此原因,在显动态分析中,刚性体不要用不切实际的杨氏模量或密度,刚体不能再变硬因为它已是完全刚硬的。
因为刚性体的质量中心的运动传递到节点上,所以不能用D命令在刚体上施加约束。刚体的一个节点上的约束和初始速度将转换到物体的质心。但是,如果约束了多个节点,就很难确定使用哪种约束。要正确在刚体上施加约束,使用EDMP命令的平移(VAL1)和转动(VAL2)约束参数域,表示如下:
3约束Y方向的旋转
4约束Z方向的旋转
5约束X,Y方向的旋转
6约束Y和Z方向的旋转
7约束Z和X方向的旋转
8约束X,Y和Z方向的旋转
例如,命令EDMP,IGID,2,7,7将约束材料的刚体单元的所有自由度。
在定义刚体之后,可以用EDIPART命令指定惯性特性、质量和初始速度矢量。如果没有定义刚性体的惯性特性,程序将会依据有限元模型计算它们。
TBDATA,2,763e6
! Tangent modulus (Pa)
TBDATA,4,40.0
! C (s-1)
TBDATA,5,5.0
! P
TBDATA,6,.75
! Failure strain
7.2.3.13分段线性塑性模型
多线性弹塑性材料模型,可输入与应变率相关的应力应变曲线。它是一个很常用的塑性准则,特别用于钢。采用这个材料模型,也可根据塑性应变定义失效。采用Cowper-Symbols模型考虑应变率的影响,它与屈服应力的关系为:
B.2.11. Plastic Kinematic Example: 1018 Steel
MP,ex,1,200e9
! Pa
MP,nuxy,1,.27
! No units
MP,dens,1,7865
! kg/m3
TB,PLAW,,,,1
TBDATA,1,310e6
! Yield stress (Pa)
这里 ——有效应变率,C和P——应变率参数, ——常应变率处的屈服应力,而 是基于有效塑性应变的硬化函数。用MP命令输入弹性模量(Exx),密度(DENS)和泊松比(NUXY)。用TB,PLAW,,,,8和TBDATA命令的1-7项输入屈服应力、切线模量、失效的有效真实塑性应变、应变率参数C、应变率参数P、定义有效全应力相对于有效塑性真应变的载荷曲线ID 以及定义应变率缩放的载荷曲线ID。
例题参看B.2.25,Rigid Material Example:Steel。
B.2.25. Rigid Material Example: Steel
MP,ex,1,207e9
! Pa
MP,nuxy,1,.3
! No units
MP,dens,1,7580
! kg/m3
EDMP,rigid,1,7,7
VAL1-平移约束参数(相对于整体笛卡尔坐标系)
0 没有约束(缺省)
1 约束X方向的位移
2 约束Y方向的位移
3 约束Z方向的位移
4 约束X和Y方向的位移
5 约束Y和Z方向的位移
6 约束Z和X方向的位移
7 约束X,Y,Z方向的位移
VAL2-转动约束参数(相对于整体笛卡尔坐标系)
1没有约束(缺省)
2约束X方向的旋转
各向同性弹性模型
各向同性弹性模型。使用MP命令输入所需参数:
MP,DENS—密度
MP,EX—弹性模量
MP,NUXY—泊松比
此部分例题参看B.2.1,Isotropic Elastic Example:High Carbon Steel。
B.2.1. Isotropic Elastic Example: High Carbon Steel
各向同性、随动硬化或各向同性和随动硬化的混合模型,与应变率相关,可考虑失效。通过在0(仅随动硬化)和1(仅各向同性硬化)间调整硬化参数β来选择各向同性或随动硬化。应变率用Cowper-Symonds模型来考虑,用与应变率有关的因数表示屈服应力,如下所示:
这里 —初始屈服应力, —应变率,C和P-Cowper Symonds为应变率参数。 —有效塑性应变, —塑性硬化模量,由下式给出:
TB,PLAW,,,, 8
TBDATA,1, (屈服应力)
TBDATA,2, (切线模量)
TBDATA,3, (失效时的有效塑性真应变)
TBDATA,4,C(应变率参数)
TBDATA,5,P(应变率参数)
TBDATA,6,LCID1(定义全真应力相对于塑性真实应变的载荷曲线)
TBDATA,7,LCID2(关于应变率缩放的载荷曲线)
TB,PLAW,,,,1
TBDATA,1, (屈服应力)
TBDATA,2, (切线模量)
TBDATA,3,β(硬化参数)
TBDATA,4,C(应变率参数)
TBDATA,5,P(应变率参数)
TBDATA,6, (失效应变)
例题参看B.2.11,Plastic Kinematic Example:1018 Steel。
用EDMP命令定义刚性体,例如,定义材料2为刚性体,执行:EDMP,RIGIS,2。用指定材料号定义的所有单元都认为是刚性体的一部分。材料号以及单元的单元类型和实常数类型号用来定义刚体的PART ID。这些 PART ID用于定义刚性体的载荷和约束(如第4章所述,Loading)。刚体内的单元不必用连接性网格连接。因此,为了在模型中表示多个独立的刚性体。必须定义多个刚体类型。但是,两个独立刚体不能共同使用一个节点。
MP,ex,1,210e9
! Pa
MP,nuxy,1,.29
! No units
MP,dens,1,7850
! kg/m3
双线性各向同性模型
使用两种斜率(弹性和塑性)来表示材料应力应变行为的经典双线性各向同性硬化模型(与应变率无关)。仅可在一个温度条件下定义应力应变特性。(也有温度相关的本构模型;参看Temperature Dependent Bilinear Isotropic Model)。用MP命令输入弹性模量(Exx),泊松比(NUXY)和密度(DENS),程序用EX和NUXY值计算体积模量(K)。用TB和TBDATA命令的1和2项输入屈服强度和切线模量:
*DIM,TruStran,,5
*DIM,TruStres,,5
TruStran(1)=0,.08,.16,.4,.75
TruStres(1)=207e6,250e6,275e6,290e6,3000e6
EDCURVE,ADD,1,TruStran (1),TruStres(1)
7.2.8.1刚性体模型
TB,BKIN
TBDATA,1, (屈服应力)
TBDATA,2, (切线模量)
例题参看B.2.10,Bilinear Kinematic Plasticity Example :Titanium Alloy。
B.2.10. Bilinear Kinematic Plasticity Example: Titanium Alloy
B.2.16. Piecewise Linear Plasticity Example: High Carbon Steel
MP,ex,1,207e9
! Pa
MP,nuxy,1,.30
! No units
MP,dens,1,7830
! kg/m3
TB,PLAW,,,,8
TBDATA,1,207e6
! Yield stress (Pa)