有限元命令整理.doc_咨信网zixin.com.cn

资源描述

U,TB,lab,,mat,ntemp,npts,tbopt,eosopt 定义非线性材料特性表 Lab: 材料特性表之种类 Bkin: 双线性随动强化 Bis 双线性等向强化 Mkin: 多线性随动强化(最多5个点) Mis 多线性等向强化（最多100个点） Dp: dp模型 Mat: 材料号 Ntemp: 数据的温度数对于bkin: ntemp缺省为6 mis ntemp缺省为1，最多20 bis ntemp缺省为6，最多为6 dp: ntemp, npts, tbopt 全用不上 Npts: 对某一给定温度数据的点数 u TBTEMP,temp,kmod 为材料表定义温度值 temp: 温度值 kmod: 缺省为定义一个新温度值如果是某一整数，则重新定义材料表中的温度值注意：此命令一发生，则后面的TBDATA和TBPT均指此温度，应该按升序若Kmod为crit, 且temp为空，则其后的tbdata数据为solid46,shell99,solid191中所述 mp,lab, mat, co, c1,…….c4 定义材料号及特性 lab: 待定义的特性项(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，C TBDATA，2，ψ /TBDATA，3，…… 如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下： MP，EX，1，1E8 MP，NUXY，1，0.3 TB，DP，1 TBDATA，1，27 TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,ded 定义非线性材料定义弹性模量： MPDATA,Lab,MAT,STLOC,C1,C2,C3,C4,C5,C6 MPTEMP,MPTEMP,STLOC,T1,T2,T3,T4,T5 激活塑性材料选项： TB,Lab,MAT,NTEMP,NPTS,TBOPT,EOSOPT 输入塑性材料数据： TETEMP,TEMP,KMOD TBDATA,STLOC,C1,C2,C3,C4,C5,C6 u *do, par, ival, fval, inc 定义一个do循环的开始 par: 循环控制变量 ival, fval, inc：起始值，终值，步长（正，负） u *enddo 定义一个do循环的结束 u *if,val1, oper, val2, base: 条件语句 val1, val2: 待比较的值（也可是字符，用引号括起来） oper: 逻辑操作（当实数比较时，误差为1e-10） eq, ne, lt, gt, le, ge, ablt, abgt base: 当oper结果为逻辑真时的行为 lable: 用户定义的行标志 stop: 将跳出anasys exit: 跳出当前的do循环 cycle: 跳至当前do循环的末尾 then: 构成if-then-else结构一个循环命令 *do, par, ival, fval, inc 定义一个do循环的开始 par: 循环控制变量 ival, fval, inc：分别为起始值，终值，步长（可正可负） …… *enddo 定义一个do循环的结束定义表、数组等真的很好用哦 dim, par, type, imax, jmax, kmax, var1, vae2, var3 定义数组 par: 数组名 type: array 数组，如同fortran,下标最小号为1，可以多达三维 char 字符串组（每个元素最多8个字符） table 表 imax,jmax, kmax 各维的最大下标号 var1,var2,var3 各维变量名，缺省为row,column,plane 在ANSYS帮助系统中关于*SET命令的注释下列出了ANSYS中可以使用的数学函数。所有这些数学函数均可以在ANSYS环境中使用，这些数学函数包括： ABS(X) 求绝对值 ACOS(X) 反余弦 ASIN(X) 反正弦 ATAN(X) 反正切 ATAN2(X,Y) 反正切, ArcTangent of (Y/X) , 可以考虑变量X,Y 的符号 COS(X) 求余弦 COSH(X) 双曲余弦 EXP(X) 指数函数 GDIS(X,Y) 求以X为均值,Y为标准差的高斯分布,在使用蒙地卡罗法研究随机荷载和随机材料参数时,可以用该函数处理计算结果 LOG(X) 自然对数 LOG10(X) 常用对数(以10为基) MOD(X,Y) 求 X/Y的余数. 如果 Y=0, 函数值为 0 NINT(X) 求最近的整数 RAND(X,Y) 取随机数,其中X 是下限, Y是上限 SIGN(X,Y) 取 X的绝对值并赋予Y的符号. Y>=0, 函数值为|X|, Y<0, 函数值为-|X|,. SIN(X) 正弦 SINH(X) 双曲正弦 SQRT(X) 平方根 TAN(X) 正切 TANH(X) 双曲正切 esel,s,mat,,1 选择材料号为1的单元 *get,emin,elem,,num,min 获得最小的单元号 *get,emax,elem,,num,max 获得最大的单元号 *DO,I,emin,emax 作循环 *GET,V1,ELEM,I,VOLU 获得单元的体积存到V1的变量中 V=V+V1 求和获得材料1的总体积 *enddo 1.*GET命令 *GET命令的使用格式为： *GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM 其中： Par是存储提取项的参数名； Entity是被提取项目的关键字，有效地关键字是NODE, ELEM, KP, LINE, 　　AREA, VOLU, PDS等； ENTNUM是实体的编号（若为０指全部实体）； Item1是指某个指定实体的项目名．例如，如果Entity是ELEM，那么Item1 　　要么是NUM（选择集中的最大或最小的单元编号），要么是COUNT　（选择集中的单元数目）．可以把*GET命令看成是对一种树型结构从上至下的路径搜索，即从一般到特殊的确定．单元生死法的使用在大多数静态和非线形瞬态分析小，都可以使用单元死活行为，与其他分析一样，分析过程包括建摸、加载并求解和查看结果3个主要步骤。 1．建立模型在PREP7中创建所有单元，包括那些在开始“死掉”，在以后的荷载少中被激活的单元。不能在求解过程中创建新的单元。 2．加载和求解 (1)指定分析类型。 (2)定义第—个荷载步。在结构分析中应激活大变形效应： ● 命令：NLGEOM,ON GUI： mainnMenu->preprocessor->Loads->Analysis Options Main Menu->Solution->Sol'n Controls Main Menu->Solution->AnalysisOptions 使用单元生死选项叫，应设置Newton-Raphson选项：命令：NROPT，Option，—，Adptky GUI： Main Menu->Preprocessor->Loads->AnalysisOptions Main Menu->Solution->Analysis Options 提示：打开自适应下降因子的全牛顿-拉普森选项通常会产生更好的结果。杀死所有要在后续荷载步“生”(激活)的单元：命令：EKILL，ELEM GUI： main Menu->Preprocessor->Loads->Other->Kill Elements Main Menu->Solution->Other->Kill Elements 重新定义刚度缩减因子： ● 命令：ESTIF,KMULT GUI： Main Menu->Preprocessor->Loads->Other->SfiffnessMult Main Menu->Solution->Other->StiffnessMult 注童：不与任何“生”的单元相连的结点将可能“漂移”，为了减少求解的方程数和避免病态条件，需要约束死的自由度。当单元变“生”时，必须删除这些人为约束。第一个荷载步的命令流示例： TIME,... !定义荷载步时间 NLGEOM，ON !打开大变形效应 NROPT,FULL !设定牛顿-拉普森选项 ESTIF,... !设定刚度端蔽因子(可选) ESEL,... !选择在本荷效步要杀死的单元 EKILL,... !杀死选择的单元 ESEL,S,LIVE !选择所有活的单元 NSLE,S !选择所有活的结点 NSEI,INVE !选择所有死的结点(不与任何活动单元相连的结点) 0,ALL,ALL,0 !约束所有死的结点(可选) NSELALL !选择所有结点 ESEL,ALL !选择所打单元 D,... !施加合适的约束 F,... !在“生”的结点上施加合适的自由度荷载 SF,... !施加合适的单元分布苘载 BF,... !施加合适的体荷载 SAVE !保存数怀库 SOLVE !求解 (3)定义后续荷载步。在后续荷载步中，用户可以杀死或激活单元。杀此单元：杀死所有要在后续荷载步“生”(激活)的单元： ' 命令：EKILL，ELEM ' GUI： Main Menu->preprocessor->Loads->Other->Kill Elements Main Menu->Solution->Other->Kill Elements 激活单元： ● 命令：EALIVE，ELEM ● GUI： Main Menu->Preprocessor->Loads->Other->Activate Elem Main Menu->Solution->Other Activate Elem 后续荷载步的命令流示例： TIME,... ESEL,... !选择要杀死的单元 EKILL,... !杀死所选择的单元 ESEL,... !选择要徼活的单元 EALIVE，... !激活所选择的单元 FDELE,... !删除死结点自由度荷载 D,... !约束死的结点自由度 F,... !在生的结点施加合适的结点自由度荷载 DDELE， !删除重新激活的结点自由度上的约束 SAVE SOLVE 注意：单元包含生死行为时，求解多个荷栽步不能使用荷载步文件法。荷载步文件不能记录单元的生死状态。 3．查看结果对包含单元生死的分析进行后处理时，方法基本按标准的后处理方法。需要指出的是， “死”的单元仍保持在模型中，在图形显示和列表输出时仍包括“死” 的单元。在使用PLNSOL 命令(GUI：Utility Menu->Plot->Results->Contour Plot->Nodal Solution)显示结点平均结果时， “死”的单元将会“污染”结果。在后处理时，建议使用ESEL命令移走“死”的单元。 4．通过ANSYS计算结果控制单元生死有些时候，无法明确指定单元“生”和”死”的时间点或空间位置。例如，热分析中要杀死已经熔化的单元，但事先并不知道单元的确切位置。单元的位置要以程序计算的温度为标准进行判断。可以使用下面的命令来判断和选择临界单元。将单元讣算结果存储到单元表：命令： ETABlE，Lab，Item，Comp ● GUI： Main Menu->General Postproc->Element Table->Define Table 使用选择命令选择临界单元： ● 命令： ESEl，Type，item，ComP，VMIN，VMAX，VINC，KABS ● GUI： Utility Menu->Select->Entities 通过训算结果控制单元生死命令流： /SOLU !进入求解器 ... !设置求解选项 SOLVE !开始求解 FINISH ! /POST1 !进入后处理器POST1 SET,... !存储总的等效匝变到单元表 ETABLE，STRAIN，EPTO，EQV !选择所有总的等效应变大于0．2的单元 ESEL,S，ETAB，STRAIN,O,20 FINISH ! /SOLU !垂新进入求解器 ANTYPE,REST !定义进行重启动分析 EKILL,ALL !杀死析选定的单元 ESEL,All !重新选择所有单元 ... !继续完成下面的求解分析

展开阅读全文