查看文章 |
ANSYS命令流中文说明 （2）
2008-05-27 15:46
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 合并相同位置的item label: 要合并的项目 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 再选一部线附加给当前选中组 au none u(unselect) inve: 反向选择 item: line 线号 loc 坐标 length 线长 comp: x,y,z kswp: 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,……node17 nset: 耦合组编号 lab: ux,uy,uz,rotx,roty,rotz node1-node17: 待耦合的节点号。如果某一节点号为负，则此节点从该耦合组中删去。如果node1=all,则所有选中节点加入该耦合组。 注意：1，不同自由度类型将生成不同编号 2，不可将同一自由度用于多套耦合组 u CPINTF, LAB, TOLER 将相邻节点的指定自由度定义为耦合自由度 LAB：UX,UY,UZ,ROTX,ROTY,ROTZ,ALL TOLER: 公差，缺省为0.0001 说明：先选中欲耦合节点，再执行此命令 2.7 定义单元表 说明：1，单元表仅对选中单元起作用，使用单元表之前务必选择一种类型的单元 2，单元表各行为选中各单元，各列为每单元的不同数据 u ETABLE, LAB, ITEM, COMP 定义单元表，添加、删除单元表某列 LAB:用户指定的列名（REFL, STAT, ERAS 为预定名称） ITEM: 数据标志（查各单元可输出项目） COMP: 数据分量标志 2.8 存盘 u save, fname, ext,dir, slab 存盘 fname : 文件名（最多32个字符）缺省为工作名 ext: 扩展名（最多32个字符）缺省为db dir: 目录名（最多64个字符）缺省为当前 slab: “all” 存所有信息 “model” 存模型信息 “solv” 存模型信息和求解信息 3 /solu u /solu 进入求解器 3.1 加边界条件 u D, node, lab, value, value2, nend, ninc, lab2, lab3, ……lab6 定义节点位移约束 Node : 预加位移约束的节点号，如果为all,则所有选中节点全加约束，此时忽略nend和ninc. Lab: ux,uy,uz,rotx,roty,rotz,all Value,value2: 自由度的数值（缺省为0） Nend, ninc: 节点范围为：node-nend，编号间隔为ninc Lab2-lab6: 将lab2-lab6以同样数值施加给所选节点。 注意：在节点坐标系中讨论 3.2 设置求解选项 u antype, status, ldstep, substep, action antype: static or 1 静力分析 buckle or 2 屈曲分析 modal or 3 模态分析 trans or 4 瞬态分析 status: new 重新分析（缺省），以后各项将忽略 rest 再分析，仅对static,full transion 有效 ldstep: 指定从哪个荷载步开始继续分析，缺省为最大的，runn数（指分析点的最后一步） substep: 指定从哪个子步开始继续分析。缺省为本目录中，runn文件中最高的子步数 action, continue: 继续分析指定的ldstep,substep 说明：继续以前的分析（因某种原因中断）有两种类型 singleframe restart: 从停止点继续 需要文件：jobname.db 必须在初始求解后马上存盘 jobname.emat 单元矩阵 jobname.esav 或 .osav : 如果.esav坏了，将.osav改为.esav results file: 不必要，但如果有，后继分析的结果也将很好地附加到它后面 注意：如果初始分析生成了.rdb, .ldhi, 或rnnn 文件。必须删除再做后继分析 步骤： （1）进入anasys 以同样工作名 （2）进入求解器，并恢复数据库 （3）antype, rest （4）指定附加的荷载 （5）指定是否使用现有的矩阵（jobname.trl）（缺省重新生成） kuse: 1 用现有矩阵 （6）求解 multiframe restart:从以有结果的任一步继续（用不着） u pred,sskey, --,lskey….. 在非线性分析中是否打开预测器 sskey: off 不作预测（当有旋转自由度时或使用solid65时缺省为off） on 第一个子步后作预测（除非有旋转自由度时或使用solid65时缺省为on） -- ： 未使用变量区 lskey: off 跨越荷载步时不作预测（缺省） on 跨越荷载步时作预测（此时sskey必须同时on） 注意：此命令的缺省值假定solcontrol为on u autots, key 是否使用自动时间步长 key:on: 当solcontrol为on时缺省为on off: 当solcontrol为off时缺省为off 1: 由程序选择（当solcontrol为on且不发生autots命令时在 .log文件中纪录“1” 注意：当使用自动时间步长时，也会使用步长预测器和二分步长 u NROPT, option,--,adptky 指定牛顿拉夫逊法求解的选项 OPTION: AUTO:程序选择 FULL:完全牛顿拉夫逊法 MODI:修正的牛顿拉夫逊法 INIT：使用初始刚阵 UNSYM：完全牛顿拉夫逊法，且允许非对称刚阵 ADPTKY:ON: 使用自适应下降因子 OFF：不使用自适应下降因子 u NLGEOM，KEY KEY: OFF:不包括几何非线性（缺省） ON：包括几何非线性 u ncnv, kstop, dlim, itlim, etlim, cplim 终止分析选项 kstop: 0 如果求解不收敛，也不终止分析 1 如果求解不收敛，终止分析和程序（缺省） 2如果求解不收敛，终止分析,但不终止程序 dlim：最大位移限制，缺省为1.0e6 itlim: 累积迭代次数限制，缺省为无穷多 etlim：程序执行时间（秒）限制，缺省为无穷 cplim：cpu时间（秒）限制，缺省为无穷 u solcontrol ,key1, key2,key3,vtol 指定是否使用一些非线性求解缺省值 key1: on 激活一些优化缺省值（缺省） CNVTOL Toler=0.5%Minref=0.01(对力和弯矩) NEQIT 最大迭代次数根据模型设定在15~26之间 ARCLEN 如用弧长法则用较ansys5.3更先进的方法 PRED 除非有rotx,y,z或solid65，否则打开 LNSRCH 当有接触时自动打开 CUTCONTROL Plslimit=15%, npoint=13 SSTIF 当NLGEOM,on时则打开 NROPT,adaptkey 关闭（除非：摩擦接触存在；单元12,26,48,49,52存在；当塑性存在且有单元20,23,24,60存在） AUTOS 由程序选择 off 不使用这些缺省值 key2: on 检查接触状态（此时key1为on） 此时时间步会以单元的接触状态（据keyopt(7)的假定）为基础 当keyopt(2)=on 时，保证时间步足够小 key3: 应力荷载刚化控制，尽量使用缺省值 空：缺省，对某些单元包括应力荷载刚化，对某些不包括（查） nopl:对任何单元不包括应力刚化 incp:对某些单元包括应力荷载刚化（查） vtol： u outres, item, freq, cname 规定写入数据库的求解信息 item: all 所有求解项 basic 只写nsol, rsol, nload, strs nsol 节点自由度 rsol 节点作用荷载 nload 节点荷载和输入的应变荷载（？） strs 节点应力 freq: 如果为n,则每n步（包括最后一步）写入一次 none: 则在此荷载步中不写次项 all: 每一步都写 last: 只写最后一步（静力或瞬态时为缺省） 3.3 定义载荷步 u nsubst, nsbstp, nsbmx, nsbmn, carry 指定此荷载步的子步数 nsbstp: 此荷载步的子步数 如果自动时间步长使用autots,则此数定义第一子步的长度；如果solcontrol打开，且3D面-面接触单元使用，则缺省为1-20步；如果solcontrol打开，并无3D接触单元，则缺省为1子步；如果solcontrol关闭，则缺省为以前指定值；如以前未指定，则缺省为1） nsbmx, nsbmn：最多，最少子步数（如果自动时间步长打开）？ u time, time 指定荷载步结束时间 注意：第一步结束时间不可为“0” u f, node, lab, value, value2, nend, ninc 在指定节点加集中荷载 node:节点号 lab: Fx,Fy,Fz,Mx,My,Mz value: 力大小 value2: 力的第二个大小（如果有复数荷载） nend,ninc：在从node到nend的节点（增量为ninc）上施加同样的力 注意：（1）节点力在节点坐标系中定义，其正负与节点坐标轴正向一致 u sfa, area, lkey, lab, value, value2 在指定面上加荷载 area: n 面号 all 所有选中号 lkey: 如果是体的面，忽略此项 lab: pres value: 压力值 u SFBEAM, ELEM, LKEY, LAB, VALI, VALJ, VAL2I, VAL2J, IOFFST, JOFFST 对梁单元施加线荷载 ELEM: 单元号，可以为ALL,即选中单元 LKEY: 面载类型号，见单元介绍。对于BEAM188，1为竖向；2为横向；3为切向 VALI,VALJ: I, J节点处压力值 VAL2I,VAL2J: 暂时无用 IOFFST, JOFFST: 线载距离I, J 节点距离 u lswrite, lsnum 将荷载与荷载选项写入荷载文件中 lsnum ：荷载步文件名的后缀，即荷载步数 当 stat 列示当前步数 init 重设为“1” 缺省为当前步数加“1” 3.3.1 注意 1. 尽量加面载，不加集中力，以免奇异点 2. 面的切向荷载必须借助面单元 3.4 求解载荷步 u lssolve, lsmin, lsmax, lsinc 读入并求解多个荷载步 lsmin, lsmax, lsinc ：荷载步文件范围 4 /post1（通用后处理） u set, lstep, sbstep, fact, king, time, angle, nset 设定从结果文件读入的数据 lstep :荷载步数 sbstep：子步数，缺省为最后一步 time： 时间点（如果弧长法则不用） nset： data set number u dscale, wn, dmult 显示变形比例 wn: 窗口号（或all）,缺省为1 dmult, 0或auto : 自动将最大变形图画为构件长的5% u pldisp, kund 显示变形的结构 kund： 0 仅显示变形后的结构 1 显示变形前和变形后的结构 2 显示变形结构和未变形结构的边缘 u *get, par, node, n, u, x(y,z) 获得节点n的x(y,z)位移给参数par 等价于函数 ux(n),uy(n),uz(z) node(x,y,z): 获得(x,y,z)节点号 arnode(x,y,z)：获得和节点n相连的面 注意：此命令也可用于/solu模块 u fsum, lab, item 对单元之节点力和力矩求和 lab: 空 在整体迪卡尔坐标系下求和 rsys 在当前激活的rsys坐标系下求和 item: 空 对所有选中单元（不包括接触元）求和 cont: 仅对接触节点求和 u PRSSOL, ITEM, COMP 打印BEAM188、BEAM189截面结果 说明：只有刚计算完还未退出ANSYS时可用，重新进入ANSYS时不可用 item comp 截面数据及分量标志 S COMP X,XZ,YZ应力分量 PRIN S1,S2,S3主应力SINT应力强度，SEQV等效应力 EPTO COMP 总应变 PRIN 总主应变，应变强度，等效应变 EPPL COMP 塑性应变分量 PRIN 主塑性应变，塑性应变强度，等效塑性应变 u plnsol, item, comp, kund, fact 画节点结果为连续的轮廓线 item: 项目（见下表） comp: 分量 kund: 0 不显示未变形的结构 1 变形和未变形重叠 2 变形轮廓和未变形边缘 fact: 对于接触的2D显示的比例系数，缺省为1 item comp discription u x,y,z,sum 位移 rot x,y,z,sum 转角 s x,y,z,xy,yz,xz 应力分量 1，2，3 主应力 Int,eqv 应力intensity,等效应力 epeo x,y,z,xy,yz,xz 总位移分量 1,2,3 主应变 Int,eqv 应变intensity,等效应变 epel x,y,z,xy,yz,xz 弹性应变分量 1，2，3 弹性主应变 Int,eqv 弹性intensity,弹性等效应变 eppl x,y,z,xy,yz,xz 塑性应变分量 u PRNSOL, item, comp 打印选中节点结果 item: 项目（见上表） comp: 分量 u PRETAB, LAB1, LAB2, ……LAB9 沿线单元长度方向绘单元表数据 LABn : 空： 所有ETABLE命令指定的列名 列名： 任何ETABLE命令指定的列名 u PLLS, LABI, LABJ, FACT, KUND 沿线单元长度方向绘单元表数据 LABI:节点I的单元表列名 LABJ:节点J的单元表列名 FACT: 显示比例，缺省为1 kund: 0 不显示未变形的结构 1 变形和未变形重叠 2 变形轮廓和未变形边缘 5 /post26 （时间历程后处理） u nsol, nvar, node, item, comp,name 在时间历程后处理器中定义节点变量的序号 nvar：变量号（从2到nv（根据numvar定义）） node: 节点号 item comp u x, y,z rot x, y,z u ESOL, NVAR, ELEM, NODE, ITEM, COMP, NAME 将结果存入变量 NVAR: 变量号，2以上 ELEM: 单元号 NODE: 该单元的节点号，决定存储该单元的哪个量，如果空，则给出平均值 ITEM: COMP: NAME: 8字符的变量名， 缺省为ITEM加COMP u rforce, nvar, node, item, comp, name 指定待存储的节点力数据 nvar: 变量号 node: 节点号 item comp F x, y.z M x, y,z name: 给此变量一个名称，8个字符 u add, ir, ia,ib,ic,name,--,--,facta, factb, factc 将ia,ib,ic变量相加赋给ir变量 ir, ia,ib,ic：变量号 name: 变量的名称 u /grid, key key: “0” 或“off” 无网络 “1”或“on” xy网络 “2”或“x” 只有x线 “3”或“y” 只有y线 u xvar, n n: “0”或“1” 将x轴作为时间轴 “n” 将x轴表示变量“n” “-1” ？ u /axlab, axis, lab 定义轴线的标志 axis: “x”或“y” lab: 标志，可长达30个字符 u plvar, nvar, nvar2, ……,nvar10 画出要显示的变量（作为纵坐标） u prvar, nvar1, ……,nvar6 列出要显示的变量 6 PLOTCONTROL菜单命令 u pbc, ilem, ……,key, min, max, abs 在显示屏上显示符号及数值 item: u 所加的位移约束 rot 所加的转角约束 key: 0 不显示符号 1 显示符号 2 显示符号及数值 u /SHOW, FNAME, EXT, VECT, NCPL 确定图形显示的设备及其他参数 FNAME: X11:屏幕 文件名：各图形将生成一系列图形文件 JPEG: 各图形将生成一系列JPEG图形文件 说明：没必要用此命令，需要的图形文件可计算后再输出 7 参数化设计语言 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,if循环至label句 ？ 8 理论手册 1.方程组解法：(1)直接解法；(2)迭代解法 (1) 直接解法：a.稀疏矩阵法；b. 波前解法 a. 稀疏矩阵法：占内存大，但运算次数少；通过变换刚度矩阵的顺序使得非零元素最少 b. 波前解法： 占内存小 波前是指在还没有一个单元被解完的时候激活的方程数？ (2) 迭代解法：JCG法；PCG法；ICCG法 JCG法：可解实数、对称、非对称矩阵 PCG法：高效求解各种矩阵（包括病态），但仅解实、对称矩阵 ICCG法：类似JCG,但更强 2. 应变密度，等效应变，应力密度，等效应力 （1）应变密度（strain intensity） 应变密度 是三个主应变 （2）等效应变 有效泊松比 ：用户由avprin 命令设定；0（如果不设定） （3）应力密度(stress intensity) 应力密度 (4) 等效应力 等效应力 或 若 则有 （弹性状态下） EGEN,ITIME, NINC, IEL1, IEL2, IEINC, MINC, TINC, RINC, CINC, SINC, DX, DY, DZ 单元复制命令是将一组单元在现有坐标下复制到其他位置， 但条件是必须先建立节点，节点之间的号码要有所关联。 ITIME:复制次数，包括自己本身。 NINC: 每次复制元素时，相对应节点号码的增加量。 IEL1,IEL2,IEINC:选取复制的元素，即哪些元素要复制。 MINC:每次复制元素时，相对应材料号码的增加量。 TINC:每次复制元素时，类型号的增加量。 RINC:每次复制元素时，实常数表号的增加量。 CINC:每次复制元素时，单元坐标号的增加量。 SINC:每次复制元素时，截面ID号的增加量。 DX, DY, DZ:每次复制时在现有坐标系统下，节点的几何位置的改变量。 mshape,key,dimension 指定网格化分时单元形状 key: 0 四边形(2D), 六面体(3D) 1 三角形(2D), 四面体(3D) dimension: 2D 二维 3D 三维 定义局部坐标： LOCAL,KCN,KCS,XC,YC,ZC,THXY,THYZ,THZX,PAR1,PAR2 KCN:坐标系统代号，大于10的任何一个号码都可以。 KCS:局部坐标系统的属性。 KCS=0 卡式坐标；KCS=1 圆柱坐标；KCS=2 球面坐标；KCS=3 自定义坐标；KCS=4 工作平面坐标；KCS=5 全局初始坐标。 XC,YC,ZC:局域坐标与整体坐标系统原点的关系。 THXY,THYZ,THZX:局域坐标与整体坐标系统X、Y、Z轴的关系。 声明单位： /UNITS,LABEL LABEL=SI （公制，米、千克、秒） LABEL=CSG （公制，厘米、克、秒） LABEL=BFT （英制，长度=ft英尺） LABEL=BIN （英制，长度=in英寸） 定义节点： N,NODE,X,Y,Z,THXY,THYZ,THZX NODE：欲建立节点的号码； X,Y,Z：节点在目前坐标系统下的坐标位置。 注意：若在圆柱坐标系统下x,y,z对应r,θ,z；在球面系统下对应r,θ,?。 定义节点的集中力： F,NODE,Lab,VALUE,VALUE2,NEND,NINC NODE:节点号码。 Lab:外力的形式。 Lab=FX,FY,FZ,MX,MY,MZ(结构力学的方向、力矩方向) =HEAT(热学的热流量) =AMP,CHRG(电学的电流、载荷) =FLUX(磁学的磁通量) VALUE:外力的大小。 NODE,NEND,NINC:选取施力节点的范围和关联，故在建立节点时应先规划节点的号码，以方便整个程序的编辑。 定义作用于元素的分布力： SFE,ELEM,LKEY,Lab,KVAL,VAL1,VAL2,VAL3,VAL4 ELEM:元素号码。 LKEY:建立元素后，依节点顺序，该分布力定义施加边或面的号码 Lab:力的形式。 Lab=PRES 结构压力 =CONV热学的对流 =HFLUX热学的热流率 VAL1~VAL4:相对应作用于元素边及面上节点的值。 例如：分布力位于编号为1的3d元素、第六个面，作用于此面的四个边上的力分别为：10，20，30，40。 SEF,1,6,PRES,,10,20,30,40 ansys中关于文件读取,保存,及退出程序的命令: /Filname,fname,key 指定新的工作文件名 fname:文件名及路径,默认为先前设置的工作路径 key: 0 使用已有的log和error文件 1 使用新的log和error,但不删除旧的. /Title,tile 指定一个标题 /Exit,slab,Fname,Ext,--, 退出程序 Slab: model, 仅保存模型数据文件(默认) solu 保存模型及求解数据 all, 保存所有的数据文件 nosave, 不保存任何数据文件 /Input,Fname,Ext,--,LIne,log 读入数据文件 Fname,文件名及目录路径,默认为先前设置的工作目录 Ext, 文件扩展名 后面的几个参数一般可以不考虑. (注): 用此命令时,文件名及目录路径都必须为英文,不能含有中文字符. /Pbc,item,--,key,min,max,abs 在显示屏上显示符号及数值 item: u, 所加的位移约束 rot, 所加的转角约束 temp 所加的温度荷载 F 所加的集中力荷载 cp 耦合节点显示 ce 所加的约束方程 acel 所加的重力加速度 all 显示所有的符号及数值 key : 0 不显示符号 1 显示符号 2 显示符号及数值 [以上只列出了一些常用的item,详细的可参考帮助文档] /plopts,vers,0 不在屏幕上显示ansys标记 1. wpoffs,xoff,yoff,zoff 移动工作平面 xoff－x方向移动的距离 yoff－y方向移动的距离 zoff－z方向移动的距离 2.csys,4 激活该局部坐标系 3.wprota，thxy，thyz，thzx 旋转工作平面 thxy－绕z轴旋转 thyz－绕x轴旋转 thzx－绕y轴旋转 4.改变划分网格后的单元 首先：esel,Type, Item, Comp, VMIN, VMAX, VINC, KABS type中有 s－选择新的单元 r－在所选中的单元中再次选单元 a－再选别的单元 u－在所选的单元中除掉某些单元 all－选中所有单元 none－不选 inve－反选刚才没有被选中的所有单元 stat－显示当前单元的情况 其中 Item, Comp一般系统默认 VMIN－选中单元的最小号 VMAX－选中单元的最大号 VINC－单元号间的间隔 KABS: 0---核对号的选取 1----取绝对值 如：esel 其次： emodif，IEL, STLOC, I1, I2, I3, I4, I5, I6, I7, I8 改变选中的单元类型为所需要的类型 5. 显示所有单元元素： /eshape,SCALE scale : 0--一般地显示面、体单元元素（系统默认） 1－－显示所有的元素 如：/eshape,1 6. eplot，all 可以看到所有单元 7.lfillt，NL1, NL2, RAD, PCENT 对两相交的线进行倒圆 NL1－第一条线号 NL2－第二条线号 RAD－圆角半径 PCENT－是否生成关键点，一般为默认 如：lfillt，1，2，0.5 D, NODE, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 -- 定义节点的自由度约束. NODE，节点编号， Lab，自由度编号，如X向，Y向等 VALUE，约束点位移，实部，VALUE2，如果位移为复数，则为虚部 NEND, NINC, ，定义的终止节点编号和节点编号增量 Lab2, Lab3, Lab4, Lab5, 该部分节点的其他自由度编号。同lab 2.LATT, MAT, REAL, TYPE, --, KB, KE, SECNUM -- 定义线的属性，有限元划分用 MAT, REAL, TYPE分别为材料，实常数，单元类型编号， KB，KE，定义截面的方向关键点，如beam18x系列，默认两关键点一致 SECNUM 定义的截面的编号 3.PLDISP, KUND 显示结构变形图 kund：0 只显示变形后的结构图 1 显示变形后的结构图＋变形前的结构图 2 显示变形后的结构图＋变形前的结构边界图 4.FLIST, NODE1, NODE2, NINC -- 列表节点力荷载 所列表的节点范围是：从节点编号node1到node2，以NINC的节点增加数 5.DLIST, NODE1, NODE2, NINC -- 列表节点约束. 所列表的节点范围是：从节点编号node1到node2，以NINC的节点增加数 6.LSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP -- 选择一组线的子集 Type 定义选择集的类型 可以为 s－选择一个新的子集，默认如此 r－从当前选择子集中选择一部分作为新的子集 a－选择一个新的子集附加到当前选择集上 inve－觉得有时比较重要，对当前子集取数学上集合的逆操作 all－选择全部的线 还有u，none，stat等选项 Item ，comp 一般取item comp line（材料mat 单元类型type 实常数R） 对应量的编号 loc坐标位置 x，y，z VMIN, VMAX, VINC，根据Item ，comp取的量，而与之对应的量的数值范围;起始量的数值，终止量的数值，量的增加数值 KSWP 0 仅选择线 1 选择线外还将与线有关的属性比如关键点，单元，节点等一起选中 在ansys下的ls-dyna中编的程序里写入 edwrite,both 可生成d3plot文件，这样可在“独立”的ls-dyna中读入该文件。这是我的经验。 wpcsys,-1,0 将工作平面与总体笛卡尔系对齐 csys,1 将激活坐标系转到总体柱坐标系 antype,static 定义分析类型为静力分析 /post1中的几个命令: set, lstep, sbstep, fact, king, time, angle, nset 设定从结果文件读入的数据 lstep :荷载步数 sbstep：子步数，缺省为最后一步 time： 时间点（如果弧长法则不用） nset： data set number dscale, wn, dmult 显示变形比例 wn: 窗口号（或all）,缺省为1 dmult, 0或auto : 自动将最大变形图画为构件长的5% pldisp, kund 显示变形的结构 kund： 0 仅显示变形后的结构 1 显示变形前和变形后的结构 2 显示变形结构和未变形结构的边缘 PRETAB, LAB1, LAB2, ……LAB9 沿线单元长度方向绘单元表数据 LABn : 空： 所有ETABLE命令指定的列名 列名： 任何ETABLE命令指定的列名 PLLS, LABI, LABJ, FACT, KUND 沿线单元长度方向绘单元表数据 LABI:节点I的单元表列名 LABJ:节点J的单元表列名 FACT: 显示比例，缺省为1 kund: 0 不显示未变形的结构 1 变形和未变形重叠 2 变形轮廓和未变形边缘 etable, lab,item,comp 将单元的某项结果制作成表格,以供pretable命令输出, lab: 字段名称,自己指定 item: 结果的顶目名称,在每个单元的说明中有(在单元说明表中冒号左边的 comp, 结果项目名称的分量,在单元说明表中冒号右边的 比如将plane42单元的x应力分量制成表 etable,sx,x,x LACAL,KCN,KCS,XC,YC,ZC,THXY,THYZ,THZX,PAR1,PAR2 定义区域坐标系统,该命令执行后,ANSYS坐标系统自动更改为新建立的坐标系统,故可以定义许多区域坐标系统,以辅助有限元模型的建立。 KCN:该区域坐标系统的确定代号,大于10的任何一个号码都可以。 KCS:该区域坐标系统的属性。0,1,2分别代表卡式坐标,圆柱坐标,球面坐标。 XC,YC,ZC:该区域坐标系统与整体坐标系统原点的关系。 THXY,THYZ,THZX:该区域坐标系统与整体系统X,Y,Z轴的关系 claer,nl1,nl2,lmesh 就是将后面的直线网格化之后的节点和元素都删除 但是共享节点依然存在 mshkey,key 声明是使用自由化网格（key=0） 对应网格（key=1） 或者是混合网格（key=2） 后面两种我因为是新手，所以不大会用，一般都用自由网格~~ 关于工作平面： KWPAVE, P1, P2, P3, P4, P5, P6, P7, P8, P9 把工作平面的中心移动到以上几点的平均点 最多9 如果只选一点，那么就是把工作平面的中心移动到此点 WPOFF, XOFF, YOFF, ZOFF 移动工作平面，注意xoff，yoff，zoff是相对当前点的移动量 而不是整体坐标 WPROT, THXY, THYZ, THZX 旋转工作平面 和上面的一样，是相对当前的工作平面选择一个角度，默认设置是角度为单位 wpstyl 关闭工作平面显示 Nummrg,label,toler, Gtoler,action,switch 合并相同位置的item label: 要合并的项目 node: 节点， Elem,单元，kp: 关键点（也合并线，面及点） mat: 材料，type: 单元类型，Real: 实常数 cp：耦合项，CE：约束项，CE: 约束方程，All：所有项 toler: 公差 Gtoler：实体公差 Action: sele 仅选择不合并 空 合并 switch: 较低号还是较高号被保留（low, high） 注意：可以先选择一部分项目，再执行合并。如果多次发生合并命令，一定要先合并节点，再合并关键点。合并节点后，实体荷载不能转化到单元，此时可合并关键点解决问题。 我也感觉和Glue效果一样，但是它有独到的好处的。 numcmp是压缩编号，对计算没有影响的。 一个条件命令 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循环的结束 一个网格划分命令 用desize定义单元尺寸时单元划分应遵守的级别 高：lesize kesize esize desize 用smartzing定义单元尺寸时单元划分应遵守的级别 高：lesize kesize smartsize 定义表、数组等真的很好用哦 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 把一个矩阵的一列加起来的方法 提取当前选择集中的结点总数存入变量aaa1； 提取当前选择集中的结点的最小结点号存入变量aaa2； 定义aaa1×2数组aaa3； 开始循环： aaa3数组的第一列存储结点号； aaa3数组的第二列存储Sx； 下一个结点号存入变量aaa2； 循环结束。 /post1 *get,aaa1,node,0,count *get,aaa2,node,0,num,min *dim,aaa3,array,aaa1,2 *do,i,1,aaa1 aaa3(i,1)=aaa2 *get,aaa3(i,2),node,aaa2,s,x aaa2=ndnext(aaa2) *enddo L, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2：在两个关键点之间定义一条线。 功能：在当前激活坐标系统下，在两个指定关键点之间生成直线或曲线。 P1,P2：线的起点和终点。 NDIV：这条线的单元划分数。一般不用，指定单元划分数推荐用LESIZE。这里需要说明一下：如果你的模型相对规则，为了得到高质量的网格，不妨在划线的时候指定单元划分数，这样，既方便又能按照自己的意愿来分网。 SPACE:间隔比。通常不用，指定间隔比推荐使用命令LESIZE。 说明： 线的形状由激活坐标系决定，直角坐标系中将产生一条直线，柱坐标系中，随关键的坐标不同可能产生直线，圆弧线或螺旋线。 KGEN, ITIME, NP1, NP2, NINC, DX, DY, DZ, KINC, NOELEM, IMOVE：通过一组关键点生成额外的关键点。 ITIME：生成操作总共执行的次数，如果要生成额外的点，该值必须大于1。 NP1,NP2,NINC：被生成的那组关键点的编号为NP1至NP2,编号增量为NINC（缺省为1）。 DX,DY,DZ：关键点在激活坐标系下的位置增量（柱坐标系和球面坐标系下要注意坐标的变换）。 KINC:生成的点集与原始点集之间的增量值。如果是0，则指定为最低可用关键点编号。 NOELEM：指定是否单元和节点也随之生成。0，生成；1，不生成。 IMOVE：指定关键点是否被移除或重新定义。0，按照ITIME要求生成额外关键点；1，移除原始关键点到新的位置，保持编号不变（ITIME,KINC,NOELEM被忽略）。 ANSYS的命令流里经常看到一些相对比较固定的代码，这些代码组合在一起构成ANSYS 的一个操作。比如，通过旋转命令将面生成体的操作，命令流如下： TYPE,2 !指定生成体的单元类型 EXTOPT,ESIZE,18,0, !指定单元划分数 EXTOPT,ACLEAR,1 !清除面网格 VROTAT,ALL,,,,,,15,16, !绕关键点15，16构成的轴将所划的面网格旋转360 这一组代码在GUI方式下通常由一个或几个对话框组成，还有的命令有时需要同时使用，否则就会出错。例如：想在自己定义的坐标系下编辑有限元模型，需要定义坐标系并激活，然后将所有节点移到当前坐标系中。命令流如下： WPRO,,,90 !将工作平面绕Y轴旋转90度 CSWPLA,11,1,1,1, !在工作平面原点创建柱坐标系，并激活 NROTAT,ALL !将所有节点旋转到激活坐标系 为了能充分发挥命令流的优势，建议各位把自己常用的代码贡献出来，这样，不仅熟悉了命令，更重要的是，掌握了一种通过命令流来实现的ANSYS操作。 |
最近读者：