ansys接触分析命令流
ANSYS命令流总结(全)
ANSYS结构分析单元功能与特性/可以组成一一些命令,一般是一种总体命令(session),三十也有特殊,比如是处理/POST1! 是注释说明符号,,与其他软件的说明是一样的,ansys不作为命令读取,* 此符号一般是APDL的标识符,也就是ansys的参数化语言,如*do ,,,*enddo等等NSEL的意思是node select,即选择节点。
s就是select,选择。
DIM是定义数组的意思。
array 数组。
MP命令用来定义材料参数。
K是建立关键点命令。
K,关键点编号,x坐标,y坐标,z坐标。
K, NPT, X, Y, Z是定义关键点,K是命令,NPT是关键点编号,XYZ是坐标。
NUMMRG, keypoint 用这个命令,要保证关键点的位置完全一样,只是关键点号不一样的才行。
这个命令对于重复的线面都可以用。
这个很简单,压缩关键。
Ngen 复制节点e,节点号码:这个命令式通过节点来形成单元NUMCMP,ALL:压缩所有编号,这样你所有的线都会按次序重新编号~你要是需要固定的线固定的标号NSUBST,100,500,50:通过指定子步数来设置载荷步的子步LNSRCH线性搜索是求解非线性代数方程组的一种技巧,此法会在一段区间内,以一定的步长逐步搜索根,相比常用的牛顿迭代法所要耗费的计算量大得多,但它可以避免在一些情况下牛顿迭代法出现的跳跃现象。
LNSRCH激活线性搜索PRED 激活自由度求解预测NEQIT指定一个荷载步中的最大子步数AUTOTS 自动求解控制打开自动时间步长.KBC -指定阶段状或者用跳板装载里面一个负荷步骤。
SPLINE:P1,P2,P3,P4,P5,P6,XV1,YV1,ZV1,XV6,YV6,ZV6(生成分段样条曲线)*DIM,Par,Type,IMAX,JMAX,KMAX,Var1,Var2,Var3(定义载荷数组的名称)【注】Par: 数组名Type:array 数组,如同fortran,下标最小号为1,可以多达三维(缺省)char 字符串组(每个元素最多8个字符)tableIMAX,JMAX,KMAX各维的最大下标号Var1,Var2,Var3 各维变量名,缺省为row,column,plane(当type为table时)/config是设置ansys配置参数的命令格式为/CONFIG, Lab, V ALUELab为参数名称value为参数值例如:/config,MXEL,10000的意思是最大单元数为10000杆单元:LINK1、8、10、11、180梁单元:BEAM3、4、23、24,44,54,188,189管单元:PIPE16,17,18,20,59,602D实体元:PLANE2,25,42,82,83,145,146,182,1833D实体元:SOLID45,46,64,65,72,73,92,95,147,148,185,186,187,191壳单元:SHELL28,41,43,51,61,63,91,93,99,143,150,181,208,209弹簧单元:COMBIN7,14,37,39,40质量单元:MASS21接触单元:CONTAC12,52,TARGE169,170,CONTA171,172,173,174,175,178矩阵单元:MATRIX27,50表面效应元:SURF153,154粘弹实体元:VISCO88,89,106,107,108,超弹实体元:HYPER56,58,74,84,86,158耦合场单元:SOLID5,PLANE13,FLUID29,30,38,SOLID62,FLUID79,FLUID80,81,SOLID98,FLUID129,INFIN110,111,FLUID116,130界面单元:INTER192,193,194,195显式动力分析单元:LINK160,BEAM161,PLANE162,SHELL163,SOLID164,COMBI16杆单元单元名称简称节点数节点自由度特性备注LINK1 2D杆 2 Ux,Uy EPCSDGB常用杆元LINK8 3D杆Ux,Uy,Uz EPCSDGBLINK103D仅受拉或仅受压杆EDGB模拟缆索的松弛及间隙LINK11 3D线性调节器EGB模拟液压缸和大转动LINK180 3D有限应变杆EPCDFGB 另可考虑粘弹塑性E-弹性(Elasticity),P-塑性(Plasticity),C-蠕变(Creep),S-膨胀(Swelling),D-大变形或大挠度(Large deflection),F-大应变(Large strain)或有限应变(Finite strain),B-单元生死(Birth and dead),G-应力刚化(Stress stiffness)或几何刚度(Geometric stiffening),A-自适应下降(Adaptive descent)等。
本文介绍了ANSYS APDL命令和轴承接触分析相关实例
本文介绍了ANSYS APDL命令和轴承接触分析相关实例。
finish/clear/filename,ex5-3 !设定工作名/title,the contact analysis of Bearing !设定工作标题!************************************!进入前处理模块!************************************/prep7*set,roll_r,5.7 !定义变量参数*set,inside_r1,15*set,inside_r2,18.5*set,outside_r1,29.8*set,outside_r2,33.3*set,deep,8*set,axes_r3,15.5*set,fillet_r4,1*set,axes_deep,20*set,fillet_r5,1!定义单元类型和属性et,1,solid45 !定义单元类型solid45MP,EX,1,2.06e5 !定义弹性模量MP,PRXY,1,0.3 !定义泊松比MP,MU,1,0.3 !定义摩擦系数!建立模型sph4,,-inside_r2-roll_r,roll_r+0.9 !生成轴承滚珠cyl4,,,inside_r1,-120,inside_r2,-60,deep !生成+z向的轴承内圈cyl4,,,outside_r1,-120,outside_r2,-60,deep !生成+z向的轴承外圈cyl4,,,inside_r1,-120,inside_r2,-60,-deep !生成-z向的轴承内圈cyl4,,,outside_r1,-120,outside_r2,-60,-deep !生成-z向的轴承外圈wpoff,0,0,-deep !偏移工作平面cyl4,,,0,-120,axes_r3,-60,-axes_deep !生成装配轴lfillt,57,59,fillet_r4,, !生成倒角线al,57,62,59,54,61 !由线生成面vrotat,32,,,,,,42,37,60,1, !旋转面生成体vdele,6,,,1 !删除编号为6的体!lfillt,31,39,fillet_r5,, !生成倒角线al,34,39,53,31,37 !由线生成面vdrag,27,,,,,,35 !拖拉面生成体vdele,4,,,1 !删除编号为4的体!vgen,,7,,,,,0.27,,,1 !平移编号为7的体!wpoff,0,0,deep !偏移工作平面csys,1 !激活柱坐标系asel,s,loc,x,inside_r2 !选择x=inside_r2的面asel,a,loc,x,outside_r1 !选择x=ouside_r1的面vsba,1,all !体被面分割vdele,4,,,1 !删除编号为4的体vdele,8,,,1 !删除编号为8的体allsel,all !选择全部图元vsel,u,volu,,7 !不选编号为7的体vglue,all !粘接全部的体!以下通过一些布尔操作以方便网格划分wpoff,0,-inside_r2-roll_r,0 !偏移工作平面vsbw,1 !用工作平面分割体1wpro,,-90, !旋转工作平面vsbw,2 !用工作平面分割体2vsbw,3 !用工作平面分割体3wpro,,,-90 !旋转工作平面vsbw,1 !用工作平面分割体1vsbw,2 !用工作平面分割体2vsbw,5 !用工作平面分割体5vsbw,6 !用工作平面分割体6!voffst,2,-4 !沿面的法向平移面2生成体voffst,9,-4 !沿面的法向平移面9生成体voffst,23,-4 !沿面的法向平移面23生成体voffst,53,-4 !沿面的法向平移面53生成体!voffst,3,4 !沿面的法向平移面3生成体voffst,25,4 !沿面的法向平移面25生成体voffst,38,4 !沿面的法向平移面38生成体voffst,58,4 !沿面的法向平移面58生成体!vovlap,all !对体进行搭接操作vdele,25,,,1 !删除编号为25的体及其所属的低阶图元vdele,32,,,1 !删除编号为32的体及其所属的低阶图元vdele,33,,,1 !删除编号为33的体及其所属的低阶图元vdele,34,,,1 !删除编号为34的体及其所属的低阶图元!vdele,31,,,1 !删除编号为31的体及其所属的低阶图元vdele,35,,,1 !删除编号为35的体及其所属的低阶图元vdele,36,,,1 !删除编号为36的体及其所属的低阶图元vdele,37,,,1 !删除编号为37的体及其所属的低阶图元vglue,all !对体进行粘接操作!划分网格esize,2 !设定网格单元尺寸mshape,0,3d !设定网格形状为六面体单元mshkey,1 !设定为映射网格划分方式vsel,s,volu,,1,3,2 !选择编号为1、3 的体vsel,a,volu,,4,5 !同时选择编号为4,5的体vsel,a,volu,,9 !同时选择编号为9的体vsel,a,volu,,12,14 !同时选择编号为12、13、14的体cm,sphere,volu !生成体的组件spherevmesh,all !对体进行网格划分!esize,1 !设定网格单元尺寸!vsel,inve,volu !对当前体选择集进行反选vsel,s,volu,,6vsel,a,volu,,22,23vsel,a,volu,,26,30vsel,a,volu,,38,40vsweep,all !对体sweep网格划分esize,1.5 !设定网格单元尺寸allsel,allvsweep,8,50,49 !设定源面和目标面并进行sweep网格划分vsweep,7,32,37 !设定源面和目标面并进行sweep网格划分!!生成耦合设置cmsel,s,sphere,volu !选择名称为sphere的组件vgen,2,all,,,,,,,0 !复制该组件cmsel,s,sphere,volu !选择名称为sphere的组件vclear,all !清除该组件包含图元的网格vdele,all,,,1 !删除该组件包含的图元!csys,1 !激活柱坐标系asel,s,loc,x,inside_r2 !选择x=inside_r2的面asel,a,loc,x,outside_r1 !同时选中x=outside_r1的面asel,u,loc,y,-90 !从当前选择集中不选y=-90的面nsla,s,1 !选择面所属的节点nrotat,all !旋转节点坐标系与当前激活坐标系平齐cpintf,ux !在重合节点生成自由度ux的耦合设置cpintf,uy !在重合节点生成自由度uy的耦合设置cpintf,uz !在重合节点生成自由度uz的耦合设置!!设定接触参数/PREP7ALLSEL,ALL !选择全部图元/COM, CONTACT PAIR CREATION - START !接触对设置开始/GSA V,cwz,gsav,,temp !将当前的图形设置保存在cwz.gsav文件中!MP,MU,1,0.3 !定义摩擦系数MAT,1 !激活材料属性1R,3 !定义实常数3REAL,3 !激活实常数3ET,2,170 !定义单元类型2ET,3,174 !定义单元类型3KEYOPT,3,9,0 !设定单元类型3的关键项9KEYOPT,3,10,1 !设定单元类型3的关键向10R,3,,,0.1, !设定法向接触刚度为0.1!生成目标面ASEL,S,,,30 !选择编号为30的面ASEL,A,,,90 !同时选中编号为90的面ASEL,A,,,98 !同时选中编号为98的面ASEL,A,,,104 !同时选中编号为104的面ASEL,A,,,113 !同时选中编号为113的面ASEL,A,,,138 !同时选中编号为138的面ASEL,A,,,143 !同时选中编号为143的面CM,AREA_TARGET,AREA !生成目标面组件target TYPE,2 !激活单元类型2NSLA,S,1 !选择面所属的节点ESLN,S,0 !选择节点依附的单元ESURF !在当前选择的单元上覆盖生成单元ESEL,ALL !选择所有的单元!生成接触面ASEL,S,,,35 !选择编号为35的面ASEL,A,,,36 !同时选中编号为36的面CM,AREA_CONTACT,AREA !生成接触面组件contact TYPE,3 !激活单元类型3NSLA,S,1 !选择面所属的节点ESLN,S,0 !选择节点依附的单元ESURF !在当前选择的单元上覆盖生成单元ALLSEL !选择全部图元ESEL,ALL !选择全部单元ESEL,S,TYPE,,2 !选择单元类型为2的单元ESEL,A,TYPE,,3 !同时选中单元类型为3的单元ESEL,R,REAL,,3 !在当前选择集中选出实常数为3的单元/PSYMB,ESYS,1 !打开单元坐标系显示/PNUM,TYPE,1 !打开单元类型编号/NUM,1 !打开颜色显示EPLOT !图形显示单元ALLSEL,ALL !选择全部图元/GRES,cwz,gsav !从cwz.gsav文件中恢复图形设置/COM, CONTACT PAIR CREATION - END !接触对结束!**********************************!进入求解模块!**********************************/solu !进入求解模块csys,1 !激活柱坐标系nsel,s,loc,x,outside_r2 !选择x=outside_r2的节点d,all,all !在节点上施加全部自由度约束asel,s,loc,y,-60 !选择y=-60的面asel,a,loc,y,-120 !同时选中y=-120的面da,all,symm !施加对称边界条件!施加装配轴的移动位移da,33,uz,2*deep !在编号为33的面上施加位移约束!非线性求解设置lnsrch,on !打开线性搜索pred,on !打开预测矫正autot,on !打开自动时间步nsubst,40,100,10 !设定子步数outres,all,all !输出所有子步上的全部数据allsel,all !选择所有图元solve !开始求解!**********************************!进入后处理模块!**********************************/POST1 !进入通用后处理器PLDISP,2 !图形显示结构变形图/DSCALE,1,1.0 !设定显示比例为1.0/EXPAND,6,POLAR,FULL,0,60,0,, !将结果扩展到360度/REPLOT !重绘当前图形PLNSOL,S,Z !图形显示z方向的应力PLNSOL,S,EQV !图形显示平均等效应力PLNSOL, CONT,PRES, 0,1.0 !接触应力等值线图PLNSOL, CONT,STA T, 0,1.0 !接触状态等值线图PLDI,2,ANMODE,10,0.5, ,0 !变形前后动画效果FINISH !后处理模块结束/EXIT,ALL !退出并保存全部数据。
ANSYS命令流解释大全
一、定义材料号及特性mp,lab, mat, co, c1,…….c4lab: 待定义的特性项目(ex,alpx,reft,prxy,nuxy,gxy,mu,dens)ex: 弹性模量nuxy: 小泊松比alpx: 热膨胀系数reft: 参考温度reft: 参考温度prxy: 主泊松比gxy: 剪切模量mu: 摩擦系数dens: 质量密度mat: 材料编号(缺省为当前材料号)c 材料特性值,或材料之特性,温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项,2次项,3次项,4次项的系数二、定义DP材料:首先要定义EX和泊松比:MP,EX,MAT,……MP,NUXY,MAT,……定义DP材料单元表(这里不考虑温度):TB,DP,MAT进入单元表并编辑添加单元表:TBDATA,1,CTBDATA,2,ψTBDATA,3,……如定义:EX=1E8,NUXY=,C=27,ψ=45的命令如下:MP,EX,1,1E8MP,NUXY,1,TB,DP,1TBDATA,1,27TBDATA,2,45这里要注意的是,在前处理的最初,要将角度单位转化到“度”,即命令:*afun,deg三、单元生死载荷步!第一个载荷步TIME,... !设定时间值(静力分析选项)NLGEOM,ON !打开大位移效果NROPT,FULL !设定牛顿-拉夫森选项ESTIF,... !设定非缺省缩减因子(可选)ESEL,... !选择在本载荷步中将不激活的单元EKILL,... !不激活选择的单元ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVE !选择所有非活动结点(不与活动单元相连的结点)D,ALL,ALL,0 !约束所有不活动的结点自由度(可选)NSEL,ALL !选择所有结点ESEL,ALL !选择所有单元D,... !施加合适的约束F,... !施加合适的活动结点自由度载荷SF,... !施加合适的单元载荷BF,... !施加合适的体载荷SAVESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D,F,SF和BF命令得到更详细的解释。
Ansys命令流大全(整理)
Ansys命令流大全(整理)1、A,P1,P2,P3,P4,P5,P6,P7,P8,P9此命令用已知的一组关键点点(P1~P9)来定义面(Area),最少使用三个点才能围成面,同时产生转围绕些面的线。
点要依次序输入,输入的顺序会决定面的法线方向。
如果超过四个点,则这些点必须在同一个平面上。
Menu Paths:Main Menu>Preprocessor>Create>Arbitrary>Through KPs2、*ABBR,Abbr,String--定义一个缩略语.Abbr:用来表示字符串"String"的缩略语,长度不超过8个字符.String:将由"Abbr"表示的字符串,长度不超过60个字符.3、ABBRES,Lab,Fname,Ext-从一个编码文件中读出缩略语.Lab:指定读操作的标题,NEW:用这些读出的缩略语重新取代当前的缩略语(默认)CHANGE:将读出的缩略语添加到当前缩略语阵列,并替代现存同名的缩略语.Ext:如果"Fname"是空的,则缺省的扩展命是"ABBR".4、ABBSA V,Lab,Fname,Ext-将当前的缩略语写入一个文本文件里Lab:指定写操作的标题,若为ALL,表示将所有的缩略语都写入文件(默认)5、add, ir, ia,ib,ic,name,--,--,facta, factb, factc将ia,ib,ic变量相加赋给ir变量ir, ia,ib,ic:变量号name: 变量的名称6、Adele,na1,na2,ninc,kswp !kswp=0时只删除掉面积本身,=1时低单元点一并删除。
7、Adrag, nl1,nl2,nl3,nl4,nl5,nl6, nlp1,nlp2,nlp3,nlp4,nlp5,nlp6 !面积的建立,沿某组线段路径,拉伸而成。
8、Afillt,na1,na2,rad !建立圆角面积,在两相交平面间产生曲面,rad为半径。
Ansys接触分析帮助文档-翻译
Ansys帮助文档-接触分析侯峰整理1.接触分析overview接触问题是高度非线性的,需要大量的计算机资源来解决这类问题。
解决这类问题时,需要你对物理问题有足够的了解,花足够的时间建立模型,再用尽量好的计算资源进行求解。
接触问题有两个很明显的难点。
第一,在对问题求解之前,我们是不知道接触区域的位置的。
在不同的载荷、材料、边界条件以及其他一些因素的情况下,表面之间可能以不可预知的奇怪的方式在较大尺度上彼此侵入与分离接触面。
第二,大多数接触问题需要考虑到摩擦。
有几种摩擦准则与模型可供选择,他们都是非线性的。
由摩擦产生的反应可能很复杂,导致求解的收敛困难。
除上述两种困难之外,在许多接触问题中,我们不得不强调多域情况下的影响,例如材料的导热率,电流强度以及在接触区域内的磁通量等。
如果在你的模型中不需要考虑摩擦的影响,且体之间的交互影响是确定的,那么,你就可以采用内部的多点约束来对模型进行约束。
另外一个选择是使用约束等式或者成对的自由度约束来进行约束。
这些外部约束方程或者耦合方程仅仅适用于小应变的情况。
除在这个guide中讨论的间接的接触问题外,ansys也能够提供采用ansys ls-dyna动态分析的系列产品进行分析。
直接分析套件对于分析暂态问题非常有用。
1.1一般的接触分类方法接触问题分为两类:刚体-刚体与柔体-柔体问题。
在刚体-柔体接触问题中,一个或多个的接触面被认为是刚性的。
一般来说,任何时候,只要是分析一个硬质材料与一个较软材料的接触问题,都被假设为刚体-柔体问题。
另外一类的问题,即柔体-柔体分析,是更加常见的一类问题。
在这类问题中,两个接触面都被认为是可变形的。
1.2接触分析的能力1.3面-面接触分析单元Ansys提供刚体-柔体接触、柔体-柔体接触的面-面分析单元。
这类分析单元采用一个“目标面”与一个“接触面”来组成一个接触对。
●目标面用TARGE 169(2-D) 与TARGE 170(3-D)两类单元●接触面用CONTA 171、CONTA 172、CONTA 173 、CONTA 174四类单元。
全面详细讲解ansys接触分析
搜集整理的资料,详细讲解ansys接触分析,有实例命令流,亲测过。
包括以下部分。
一般的接触分类ANSYS接触能力点─点接触单元点─面接触单元面─面的接触单元执行接触分析面─面的接触分析.接触分析的步骤:步骤1:建立模型,并划分网格步骤二:识别接触对步骤三:定义刚性目标面步骤4:定义柔性体的接触面步骤5:设置实常数和单元关键字步骤六:步骤7:给变形体单元加必要的边界条件步骤8:定义求解和载步选项第十步:检查结果点─面接触分析点─面接触分析的步骤点-点的接触接触分析实例(GUI方法)非线性静态实例分析(命令流方式)接触分析接触问题是一种高度非线性行为,需要较大的计算资源,为了进行实为有效的计算,理解问题的特性和建立合理的模型是很重要的。
接触问题存在两个较大的难点:其一,在你求解问题之前,你不知道接触区域,表面之间是接触或分开是未知的,突然变化的,这随载荷、材料、边界条件和其它因素而定;其二,大多的接触问题需要计算摩擦,有几种摩擦和模型供你挑选,它们都是非线性的,摩擦使问题的收敛性变得困难。
一般的接触分类接触问题分为两种基本类型:刚体─柔体的接触,半柔体─柔体的接触,在刚体─柔体的接触问题中,接触面的一个或多个被当作刚体,(与它接触的变形体相比,有大得多的刚度),一般情况下,一种软材料和一种硬材料接触时,问题可以被假定为刚体─柔体的接触,许多金属成形问题归为此类接触,另一类,柔体─柔体的接触,是一种更普遍的类型,在这种情况下,两个接触体都是变形体(有近似的刚度)。
ANSYS接触能力ANSYS支持三种接触方式:点─点,点─面,平面─面,每种接触方式使用的接触单元适用于某类问题。
为了给接触问题建模,首先必须认识到模型中的哪些部分可能会相互接触,如果相互作用的其中之一是一点,模型的对立应组元是一个结点。
如果相互作用的其中之一是一个面,模型的对应组元是单元,例如梁单元,壳单元或实体单元,有限元模型通过指定的接触单元来识别可能的接触匹对,接触单元是覆盖在分析模型接触面之上的一层单元,至于ANSTS 使用的接触单元和使用它们的过程,下面分类详述。
ANSYS接触和出图技巧
ANSYS接触和出图技巧1.ANSYS后处理时如何按灰度输出云图?1)你可以到utilitymenu-plotctrls-style-colors-window colors 试试2)直接utilitymenu-plotctrls-redirect plots2 将云图输出为JPG菜单->PlotCtrls->Redirect Plots->To JPEG Files3.怎么在计算结果实体云图中切面?命令流/cplane/type图形界面操作<1.设置工作面为切面<2.PlotCtrls-->Style-->Hidden line Options将[/TYPE]选项选为section将[/CPLANE]选项选为working plane4.非线性计算过程中收敛曲线实时显示solution>load step opts>output ctrls>grph solu track>on5.运用命令流进行计算时,一个良好的习惯是:使用SELECT COMMEND后.........其后再加上ALLSEL.........6.应力图中左侧的文字中,SMX与SMN分别代表最大值和最小值如你plnsolv,s,eqv则SMX与SMN分别代表最大值等效应力和最小值等效应力如你要看的是plnsolv,u则SMX与SMN分别代表位移最大值和位移最小值不要被S迷惑mx(max)mn(min)7.在非线性分析中,如何根据ansys的跟踪显示来判断收敛?在ansys output windows 有force convergenge valu 值和criterion 值当前者小于后者时,就完成一次收敛你自己可以查看两条线的意思分别是:F L2:不平衡力的2范数F CRIT:不平衡力的收敛容差,如果前者大于后者说明没有收敛,要继续计算当然如果你以弯矩M为收敛准则那么就对应M L2 和M CRIT希望你现在能明白8.两个单元建成公共节点,就成了刚性连接,不是接触问题了。
ansys命令流中文说明
ansys命令流中文说明展开全文KB、KE: 待划分线的定向关键点起始、终止号SECNUM: 截面类型号u SECPLOT,SECID,MESHKEY 画梁截面的几何形状及网格划分SECID:由SECTYPE命令分配的截面编号MESHKEY:0:不显示网格划分1:显示网格划分u /ESHAPE, SCALE 按看似固体化分的形式显示线、面单元SCALE: 0:简单显示线、面单元1:使用实常数显示单元形状u esurf, xnode, tlab, shape 在已存在的选中单元的自由表面覆盖产生单元xnode: 仅为产生surf151 或surf152单元时使用tlab: 仅用来生成接触元或目标元top 产生单元且法线方向与所覆盖的单元相同,仅对梁或壳有效,对实体单元无效Bottom产生单元且法线方向与所覆盖的单元相反,仅对梁或壳有效,对实体单元无效Reverse 将已产生单元反向Shape: 空与所覆盖单元形状相同Tri 产生三角形表面的目标元注意:选中的单元是由所选节点决定的,而不是选单元,如同将压力加在节点上而不是单元上u Nummrg,label,toler, Gtoler,action,switch 合并相同位置的itemlabel: 要合并的项目node: 节点, Elem,单元,kp: 关键点(也合并线,面及点)mat: 材料,type: 单元类型,Real: 实常数cp:耦合项,CE:约束项,CE: 约束方程,All:所有项toler: 公差Gtoler:实体公差Action: sele 仅选择不合并空合并switch: 较低号还是较高号被保留(low, high)注意:可以先选择一部分项目,再执行合并。
如果多次发生合并命令,一定要先合并节点,再合并关键点。
合并节点后,实体荷载不能转化到单元,此时可合并关键点解决问题。
u Lsel, type, item, comp, vmin, vmax, vinc, kswp 选择线type: s 从全部线中选一组线r 从当前选中线中选一组线a 再选一部线附加给当前选中组aunoneu(unselect)inve: 反向选择item: line 线号loc 坐标length 线长comp: x,y,zkswp: 0 只选线1 选择线及相关关键点、节点和单元u Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备Type: S: 选择一组新节点(缺省)R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0” 使用正负号“1”仅用绝对值u NSLL,type, nkey 选择与所选线相联系的节点u nsla, type, nkey: 选择与选中面相关的节点type:s 选一套新节点r 从已选节点中再选a 附加一部分节点到已选节点u 从已选节点中去除一部分nkey: 0 仅选面内的节点1 选所有和面相联系的节点(如面内线,关键点处的节点)u esel, type, item, comp, vmin, vmax, vinc, kabs 选择一组单元Type: S: 选择一组单元(缺省)R: 在当前组中再选一部分作为一组A: 为当前组附加单元U: 在当前组中不选一部分单元All: 选所有单元None: 全不选Inve: 反向选择当前组(?)Stat: 显示当前选择状态Item: Elem: 单元号Type: 单元类型号Mat: 材料号Real: 实常数号Esys: 单元坐标系号u ALLSEL, LABT, ENTITY 选中所有项目LABT: ALL: 选所有项目及其低级项目BELOW: 选指定项目的直接下属及更低级项目ENTITY: ALL: 所有项目(缺省)VOLU:体高级AREA:面LINE :线KP:关键点ELEM:单元NODE:节点低级u Tshap,shape 定义接触目标面为2D、3D的简单图形Shape: line:直线Arc:顺时针弧Tria:3点三角形Quad:4点四边形………….2.6 根据需要耦合某些节点自由度u cp, nset, lab,,node1,node2,……node17nset: 耦合组编号lab: ux,uy,uz,rotx,roty,rotznode1-node17: 待耦合的节点号。
简单易懂的ANSYS接触分析教程
一般的接触分类 (1)ANSYS接触能力 (2)点─点接触单元 2点─面接触单元 2面─面的接触单元 2执行接触分析 (3)面─面的接触分析 3接触分析的步骤: 3步骤1:建立模型,并划分网格 (3)步骤2:识别接触对 (4)步骤2:指定接触面和目标面 3步骤4:定义刚性目标面 (3)步骤5:定义柔性体的接触面 (5)步骤6:设置实常数和单元关键字 (7)步骤7:控制刚体目标的运动 (13)步骤8:给变形体单元加必要的边界条件 (14)步骤9:定义求解和载荷步选项14第十步:检查结果 (15)点─面接触分析 (16)点─面接触分析的步骤 (16)点-点的接触 (22)接触分析实例(GUI方法) (23)非线性静态实例分析(命令流方式) (25)接触分析接触问题是一种高度非线性行为,需要较大的计算资源,为了进行实为有效的计算,理解问题的特性和建立合理的模型是很重要的。
接触问题存在两个较大的难点:其一,在你求解问题之前,你不知道接触区域,表面之间是接触或分开是未知的,突然变化的,这随载荷、材料、边界条件和其它因素而定;其二,大多的接触问题需要计算摩擦,有几种摩擦和模型供你挑选,它们都是非线性的,摩擦使问题的收敛性变得困难。
一般的接触分类接触问题分为两种基本类型:刚体─柔体的接触,半柔体─柔体的接触,在刚体─柔体的接触问题中,接触面的一个或多个被当作刚体,(与它接触的变形体相比,有大得多的刚度),一般情况下,一种软材料和一种硬材料接触时,问题可以被假定为刚体─柔体的接触,许多金属成形问题归为此类接触,另一类,柔体─柔体的接触,是一种更普遍的类型,在这种情况下,两个接触体都是变形体(有近似的刚度)。
ANSYS接触能力ANSYS支持三种接触方式:点─点,点─面,平面─面,每种接触方式使用的接触单元适用于某类问题。
为了给接触问题建模,首先必须认识到模型中的哪些部分可能会相互接触,如果相互作用的其中之一是一点,模型的对立应组元是一个结点。
(完整版)史上最全的ANSYS命令流大全
《史上最全的ANSYS命令流查询与解释》【1】*************************************************************************************对ansys主要命令的解释1,/PREP7 ! 加载前处理模块2,/CLEAR,NOSTART ! 清除已有的数据, 不读入启动文件的设置(不加载初始化文件)初始化文件是用于记录用户和系统选项设置的文本文件/CLEAR, START !清除系统中的所有数据,读入启动文件的设置/FILENAME, EX10.5 ! 定义工程文件名称/TITLE, EX10.5 SOLID MODEL OF AN AXIAL BEARING ! 指定标题4,F,2,FY,-1000 ! 在2号节点上施加沿着-Y方向大小为1000N的集中力6,FINISH ! 退出模块命令7,/POST1 ! 加载后处理模块8,PLDISP,2 ! 显示结构变形图,参数“2”表示用虚线绘制出原来结构的轮廓9,ETABLE,STRS,LS,1 ! 用轴向应力SAXL的编号”LS,1”定义单元表STRSETABLE, MFORX,SMISC,1 ! 以杆单元的轴力为内容, 建立单元表MFORXETABLE, SAXL, LS, 1 ! 以杆单元的轴向应力为内容, 建立单元表SAXLETABLE, EPELAXL, LEPEL, 1 ! 以杆单元的轴向应变为内容, 建立单元表EPELAXLETABLE,STRS_ST,LS,1 !以杆件的轴向应力“LS,1”为内容定义单元表STRS_STETABLE, STRS_CO, LS,1 !以杆件的轴向应力“LS,1”定义单元表STRS_COETABLE,STRSX,S,X ! 定义X方向的应力为单元表STRSXETABLE,STRSY,S,Y ! 定义Y方向的应力为单元表STRSY*GET,STRSS_ST,ELEM,STEEL_E, ETAB, STRS_ST !从单元表STRS_ST中提取STEEL_E单元的应力结果,存入变量STRSS_ST;*GET, STRSS_CO,ELEM,COPPER_E,ETAB,STRS_CO”从单元表STRS_CO中提取COPPER_E单元的应力结果,存入变量STRSS_CO10 FINISH !退出以前的模块11, /CLEAR, START ! 清除系统中的所有数据,读入启动文件的设置12 /UNITS, SI !申明采用国际单位制14 /NUMBER, 2 !只显示编号, 不使用彩色/NUMBER, 0 ! 显示编号, 并使用彩色15 /SOLU ! 进入求解模块:定义力和位移边界条件,并求解ANTYPE, STATIC ! 申明分析类型是静力分析(STA TIC或者0)OUTPR, BASIC, ALL ! 在输出结果中, 列出所有荷载步的基本计算结果OUTPR,BASIC,ALL !指定输出所有节点的基本数据OUTPR,BASIC,LAST ! 选择基本输出选项,直到最后一个荷载步OUTPR,,1 ! 输出第1个荷载步的基本计算结果OUTPR,BASIC,1 ! 选择第1荷载步的基本输出项目OUTPR,NLOAD,1 ! 指定输出第1荷载步的内容OUTRES,ALL,0 !设置将所有数据不记录到数据库。
Ansys命令流大全
1、A,P1,P2,P3,P4,P 5,P6,P7,P8,P9欧阳学文此命令用已知的一组关键点点(P1~P9)来定义面(Area),最少使用三个点才能围成面,同时产生转围绕些面的线。
点要依次序输入,输入的顺序会决定面的法线方向。
如果超过四个点,则这些点必须在同一个平面上。
Menu Paths:MainMenu>Preprocessor>Create>A rbitrary>Through KPs 2、*ABBR,Abbr,String--定义一个缩略语.Abbr:用来表示字符串"String"的缩略语,长度不超过8个字符.String:将由"Abbr"表示的字符串,长度不超过60个字符.3、ABBRES,Lab,Fname,Ext-从一个编码文件中读出缩略语.Lab:指定读操作的标题,NEW:用这些读出的缩略语重新取代当前的缩略语(默认)CHANGE:将读出的缩略语添加到当前缩略语阵列,并替代现存同名的缩略语.Ext:如果"Fname"是空的,则缺省的扩展命是"ABBR".4、ABBSAV,Lab,Fname,Ext-将当前的缩略语写入一个文本文件里Lab:指定写操作的标题,若为ALL,表示将所有的缩略语都写入文件(默认)5、add, ir, ia,ib,ic,name,,,facta, factb, factc将ia,ib,ic变量相加赋给ir变量ir, ia,ib,ic:变量号name: 变量的名称6、Adele,na1,na2,ninc,kswp !kswp=0时只删除掉面积本身,=1时低单元点一并删除。
7、Adrag, nl1,nl2,nl3,nl4,nl5,nl6, nlp1,nlp2,nlp3,nlp4,nlp5,nlp6 !面积的建立,沿某组线段路径,拉伸而成。
8、Afillt,na1,na2,rad !建立圆角面积,在两相交平面间产生曲面,rad为半径。
9、*AFUN,Lab在参数表达式中,为角度函数指定单位.Lab:指定将要使用的角度单位.有3个选项.RAD:在角度函数的输入与输出中使用弧度单位(默认)DEG:在角度函数的输入与输出中使用度单位.STAT:显示该命令当前的设置(即是度还是弧度).10、Agen,itime,na1,na2,ninc,dx,dy,dz,kinc, noelem,imove !面积复制命令。
最新ANSYS命令流学习笔记3-Solid单元的接触分析资料
!ANSYS命令流学习笔记3-Solid单元的接触分析背景说明:两个半径均为100mm,的正交圆柱体发生正接触,作用在两圆柱接触体法线方向上的压力总和1000N,两圆柱体均为钢制分析两圆柱的接触情况。
分析思路:对模型进行对称约束,施加载荷后要进行节点自由度耦合,使加载面的所有节点能在承载后具有相同的Y方向位移,并且在变形后仍然保持为水平平面。
总结:1. 如果不进行节点耦合,结构变形会产生刚性位移,造成不收敛。
但是节点耦合造成了变形不符合实际情况,即忽略了施力面的泊松比造成的影响。
2. 由面网格生成体网格,几何只有面,有限元模型是体,需要定义网格单元,定义拉伸的目标单元,对拉伸选项做出设置。
3. 接触部分的细分十分重要,也是收敛的影响因素。
可以把下列文字,直接粘贴到TXT文档修改,更为方便。
! 本次学习重点:!1、命令流的复杂建模!重点学习下建模。
因为前处理尽可能用WB,所以这里也是规则模型。
还是太不好用了。
我选择狗带。
!2、接触的定义。
!Ansys依据实常数识别接触对,接触对要具有同样的实常数定义。
而且命令流是根据选择的单元来定义接触单元。
GUI操作,也不麻烦。
!3、分析步数的设定!非线性分析时步数设置,线性搜索,自动步长,定义子步。
!求解解析解finish/clearf=1000e=200r=100 !命令不区分大小写,参数也不区分大小写。
p=2.45*(f*e**2*((r+r)/(r*r))**2)**(1/3) !接触应力的解析解61.73,有限元法结果59.9MPa !正式分析命令流finish/clear/prep7 !进入前处理et,1,solid185et,2,mesh200keyopt,2,1,6 !查看help中关于mesh200和keyopt的解释,KEYOPT, ITYPE, KNUM, VALUE。
此处意思为定义了mesh(2)为4 nodes的四边形单元(1,6)。
mp,ex,1,2e11mp,prxy,1,0.3 !定义et 1的材料属性,Elastic moduli 和Poisson's ratios,(lab,mat,co)!lab这里是系统规定的,输出结果表格时lab是可以自定义的。
ANSYS Mechanical 接触分析
基本概念
接触是状态改变 非线性。也就是说, 系统的刚度依赖于接触状态:
状态 = 开
状态 = 闭合和粘结
状态= 闭合和滑动
基本概念
接触是强非线性,因为随着接触状态的改变,接触表面的法向和切 向刚度都有显著的变化。 刚度方面大的突变通常会导致严重的收敛困难。
开接触
闭接触
F
u
基本概念
使接触分析复杂化的因素包括:
基本概念
接触分类
接触问题通常分为两类:
刚性体 -柔性体 和柔性体 -柔性体
刚性体 - 柔性体
一个或更多的接触表面看作刚性体(一个表面与其它表面相比, 刚度显 然要大得多),许多金属成形问题可归为此类。 不计算刚性体内的应力
柔性体 - 柔性体
两个或所有的接触体都可变形(所有表面的刚度相近),螺栓结合凸缘连 接就是柔性体 - 柔性体接触的例子。
接触刚度例子
对Normal Stiffness factors 为0.10, 1.0, 和10.0重复以上分析.
将接触算法 改为 Normal Lagrange进行进一步对比.
接触刚度例子
结论
刚度增加, 渗透减少,而最大压力增加,并且通常会有更多的迭代和 更长运行时间。 指定正确的法向刚度不是一个独立问题,它总是需在结果质量 (精度) 和耗费(运行时间)间寻找平衡。基于这个研究, 法向接触刚度因子为 1.0 是满意的。
SEQV
寻找 FKN 值, 超过此值, 结果 不发生显著的改变
7. 重复步骤 5和6, 直到达到所预期的收敛
FKN
接触刚度
作为一个例子, 对于一个轴上套环的过盈配合分析, FKN 是变 化的,等效von Mises应力的最大值被监控。
ANSYS高级接触问题处理
• -开始估计时,选用 • FKN = 1.0 大面积实体接触 • FKN = 0.01-0.1 较柔软(弯曲占主导的)部分 • -另外,也可以指定一个绝对刚度值,单位:(力/长度)
/ 面积。 • ·点一点(除CONTA178)和点-面接触单元需要为罚刚
度KN输入绝对值: • -初始估计时: • 对于大变形: 0.1*E < KN < 1.0*E • 对于弯曲: 0.01*E < KN < 0.1*E • E 为弹性模量
ANSYS 高级接触问题
• 接触问题概述 • 在工程中会遇到大量的接触问题,如齿轮的啮合、法兰
联接、机电轴承接触、卡头与卡座、密封、板成形、冲 击等等。接触是典型的状态非线性问题,它是一种高度 非线性行为。接触例子如图1:
• 分析中常常需要确定两个或多个相互接触 物体的位移、接触区域的大小和接触面上 的应力分布。
触
• 3. 点-点接触单元用于模拟单点和另一个确定点 之间的接触。
• -建立模型时必须事先知道确切的接触位置;
• -多个点-点接触单元可以模拟两个具有多个单 元表面间的接触;
• ·每个表面的网格必须是相同的; • ·相对滑动必须很小; • ·只对小的转动响应有效。
• 例如: 点一点接触可以模拟一些面的接触。如地 基和土壤的接触
• 3、选取接触刚度的指导:
• Step 1.开始采用较小的刚度值 • Step 2.对前几个子步进行计算 • Step 3.检查穿透量和每一个子步中的平衡迭代次数
• ·在粗略的检查中,如以实际比例显示整个模型时就能观察到穿透, 则穿透可能太大了,需要提高刚度重新分析。
• ·如果收敛的迭代次数过多(或未收敛),降低刚度重新分析。 • 注意:罚刚度可以在载荷步间改变,并且可以在重启动中调整。 • 牢记:接触刚度是同时影响计算精度和收敛性的最重要的参数。如
ANSYS高级接触分析
图3-3
• 例:Seal.dat (图3-3) • Step 5. 建立接触面单元 • 设置接触单元属性、选择可变形体表面节点,并在可 变形体上建立接触单元(过程与在可变形体上建立目标单 元相同) Main Menu > Preprocessor > Modeling > Create > Elements > Surf/Contact > Surf to Surf (ESURF) • 这些接触单元与基体有同样的阶数(低阶或高阶)。 • 注意,在壳或梁单元上建立目标单元或接触单元时, 可以选择要在梁或壳单元的顶层还是底层建立单元。
例如: 超弹密封
• Step 3. 设置单元选项和实常数 • 接触对由实常数号来定义,接触单元和目标单元必须具有相同的实常数。 • Step 4. 建立目标单元(网格) • 此步中所采用的方法依赖于目标面是刚性的还是柔性的。 • 刚性目标面采用: • 直接生成 (E 命令) • 自动划分 (LMESH, AMEAH) • 可变形目标面采用 • Main Menu > Preprocessor > Modeling > Create > Elements > Surf/Contact > Surf to Surf(ESURF) • 对于直接生成刚性目标面,在建立目标单元之前需要要指定附加的单元属 性 TSHAP
• 2、接触刚度的选取
选定一个合适的接触刚度值需要一些经验。 • 对于面一面接触单元,接触刚度通常指定为基体单 元刚度的一个比例因子。 ◦ 开始估计时,选用 ◦ FKN = 1.0 大面积实体接触 ◦ FKN = 0.01-0.1 较柔软(弯曲占主导的部分) • 另外,也可以指定一个绝对刚度值,单位:(力/ 长度)/ 面积。 • 对于点一点(除 CONTA178)和点-面接触单元需 要为罚刚度 KN 输入绝对值: ◦ 初始估计时: ◦ 对于大变形: 0.1*E < KN < 1.0*E ◦ 对于弯曲: 0.01*E < KN < 0.1*E ◦ E 为弹性模量
ANSYS命令流学习笔记3-Solid单元的接触分析
!ANSYS命令流学习笔记3-Solid 单元的接触分析背景说明:两个半径均为100mm,的正交圆柱体发生正接触,作用在两圆柱接触体法线方向上的压力总和1000N,两圆柱体均为钢制分析两圆柱的接触情况。
分析思路:对模型进行对称约束,施加载荷后要进行节点自由度耦合,使加载面的所有节点能在承载后具有相同的Y方向位移,并且在变形后仍然保持为水平平面。
总结:1. 如果不进行节点耦合,结构变形会产生刚性位移,造成不收敛。
但是节点耦合造成了变形不符合实际情况,即忽略了施力面的泊松比造成的影响。
2. 由面网格生成体网格,几何只有面,有限元模型是体,需要定义网格单元,定义拉伸的目标单元,对拉伸选项做出设置。
3. 接触部分的细分十分重要,也是收敛的影响因素。
可以把下列文字,直接粘贴到TXT文档修改,更为方便。
! 本次学习重点:!1、命令流的复杂建模!重点学习下建模。
因为前处理尽可能用WB,所以这里也是规则模型。
还是太不好用了。
我选择狗带。
!2、接触的定义。
!Ansys 依据实常数识别接触对,接触对要具有同样的实常数定义。
而且命令流是根据选择的单元来定义接触单元。
GUI 操作,也不麻烦。
!3、分析步数的设定!非线性分析时步数设置,线性搜索,自动步长,定义子步。
!求解解析解finish/clearf=1000e=200r=100 ! 命令不区分大小写,参数也不区分大小写。
p=2.45*(f*e**2*((r+r)/(r*r))**2)**(1/3) !接触应力的解析解61.73,有限元法结果59.9MPa!正式分析命令流finish/clear/prep7 !进入前处理et,1,solid185et,2,mesh200keyopt,2,1,6 !查看help 中关于mesh200 和keyopt 的解释,KEYOPT, ITYPE, KNUM, VALU。
E 此处意思为定义了mesh(2)为4 nodes 的四边形单元(1,6)。