分享

ANSYS常用命令合集(二) ANSYS ANSYS常用命令,ANSYS命令,CAE AN...

 牧轩 2012-03-03

2011年02月23日  中华工程师网  

2.6根据需要耦合某些节点自由度

    ◆cp, nset, lab,,node1,node2,……node17

    nset: 耦合组编号

    lab:  ux,uy,uz,rotx,roty,rotz

    node1-node17: 待耦合的节点号。如果某一节点号为负,则此节点从该耦合组中删去。如果node1=all,则所有选中节点加入该耦合组。

    注意:1,不同自由度类型将生成不同编号

                2,不可将同一自由度用于多套耦合组

    ◆CPINTF, LAB, TOLER       将相邻节点的指定自由度定义为耦合自由度

    LAB:UX,UY,UZ,ROTX,ROTY,ROTZ,ALL

    TOLER: 公差,缺省为0.0001

    说明:先选中欲耦合节点,再执行此命令

2.7定义单元表

    说明:1,单元表仅对选中单元起作用,使用单元表之前务必选择一种类型的单元

               2,单元表各行为选中各单元,各列为每单元的不同数据

    ◆ETABLE, LAB, ITEM, COMP    定义单元表,添加、删除单元表某列

      LAB:用户指定的列名(REFL, STAT, ERAS 为预定名称)

      ITEM: 数据标志(查各单元可输出项目)

      COMP: 数据分量标志

2.8存盘

    ◆save, fname, ext,dir, slab    存盘               

    fname : 文件名(最多32个字符)缺省为工作名

    ext:    扩展名(最多32个字符)缺省为db
   
    dir:    目录名(最多64个字符)缺省为当前

    slab:   “all” 存所有信息

            “model” 存模型信息

            “solv”   存模型信息和求解信息

3/solu

    ◆/solu           进入求解器

3.1加边界条件

    ◆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设置求解选项

    ◆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:从以有结果的任一步继续(用不着)

    ◆pred,sskey, --,lskey…..              在非线性分析中是否打开预测器

    sskey: off       不作预测(当有旋转自由度时或使用solid65时缺省为off)

    on       第一个子步后作预测(除非有旋转自由度时或使用solid65时缺省为on)

    -- :         未使用变量区

    lskey:  off         跨越荷载步时不作预测(缺省)

                 on         跨越荷载步时作预测(此时sskey必须同时on)

    注意:此命令的缺省值假定solcontrol为on

    ◆autots, key        是否使用自动时间步长

    key:on: 当solcontrol为on时缺省为on

           off: 当solcontrol为off时缺省为off

         1:  由程序选择(当solcontrol为on且不发生autots命令时在 .log文件中纪录“1”

    注意:当使用自动时间步长时,也会使用步长预测器和二分步长

    ◆NROPT, option,--,adptky    指定牛顿拉夫逊法求解的选项

    OPTION: AUTO:程序选择

      FULL:完全牛顿拉夫逊法

      MODI:修正的牛顿拉夫逊法

      INIT:使用初始刚阵

      UNSYM:完全牛顿拉夫逊法,且允许非对称刚阵

     ADPTKY:ON: 使用自适应下降因子

                      OFF:不使用自适应下降因子

    ◆NLGEOM,KEY

    KEY: OFF:不包括几何非线性(缺省)

              ON:包括几何非线性

    ◆ncnv, kstop, dlim, itlim, etlim, cplim                    终止分析选项

    kstop: 0 如果求解不收敛,也不终止分析

           1 如果求解不收敛,终止分析和程序(缺省)
 
           2如果求解不收敛,终止分析,但不终止程序

    dlim:最大位移限制,缺省为1.0e6

    itlim:  累积迭代次数限制,缺省为无穷多

    etlim:程序执行时间(秒)限制,缺省为无穷

    cplim:cpu时间(秒)限制,缺省为无穷

    ◆solcontrol ,key1, key2,key3,vtol     指定是否使用一些非线性求解缺省值

    key1: on 激活一些优化缺省值(缺省)

CNVTOL

Toler=0.5%

Minref=0.01(对力和弯矩)

NEQIT

最大迭代次数根据模型设定在15~26之间

ARCLEN

如用弧长法则用较ansys5.3更先进的方法

PRED

除非有rotx,y,zsolid65,否则打开

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:

    ◆outres, item, freq, cname              规定写入数据库的求解信息

    item: all  所有求解项

          basic   只写nsol, rsol, nload, strs

          nsol    节点自由度

          rsol    节点作用荷载

          nload   节点荷载和输入的应变荷载(?)

          strs     节点应力

    freq: 如果为n,则每n步(包括最后一步)写入一次

          none:  则在此荷载步中不写次项

          all:    每一步都写

          last:   只写最后一步(静力或瞬态时为缺省)

3.3定义载荷步

    ◆nsubst, nsbstp, nsbmx, nsbmn, carry     指定此荷载步的子步数

    nsbstp:  此荷载步的子步数

            如果自动时间步长使用autots,则此数定义第一子步的长度;如果solcontrol打开,且3D面-面接触单元使用,则缺省为1-20步;如果solcontrol打开,并无3D接触单元,则缺省为1子步;如果solcontrol关闭,则缺省为以前指定值;如以前未指定,则缺省为1)

    nsbmx, nsbmn:最多,最少子步数(如果自动时间步长打开)?

    ◆time, time               指定荷载步结束时间

    注意:第一步结束时间不可为“0”

◆f, node, lab, value, value2, nend, ninc        在指定节点加集中荷载

    node:节点号

    lab: Fx,Fy,Fz,Mx,My,Mz

    value: 力大小

    value2: 力的第二个大小(如果有复数荷载)

    nend,ninc:在从node到nend的节点(增量为ninc)上施加同样的力

    注意:(1)节点力在节点坐标系中定义,其正负与节点坐标轴正向一致

    ◆sfa, area, lkey, lab, value, value2             在指定面上加荷载

    area: n 面号

    all 所有选中号

    lkey: 如果是体的面,忽略此项

    lab: pres

    value: 压力值

    ◆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 节点距离

    ◆lswrite, lsnum              将荷载与荷载选项写入荷载文件中

    lsnum :荷载步文件名的后缀,即荷载步数

    当 stat  列示当前步数

    init  重设为“1”

    缺省为当前步数加“1”

3.3.1注意

    1. 尽量加面载,不加集中力,以免奇异点

    2. 面的切向荷载必须借助面单元

3.4求解载荷步

    ◆lssolve, lsmin, lsmax, lsinc                读入并求解多个荷载步

    lsmin, lsmax, lsinc  :荷载步文件范围

4/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    显示变形结构和未变形结构的边缘

    *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模块

    ◆fsum, lab, item                            对单元之节点力和力矩求和

    lab: 空  在整体迪卡尔坐标系下求和

    rsys 在当前激活的rsys坐标系下求和

    item: 空 对所有选中单元(不包括接触元)求和

    cont: 仅对接触节点求和

◆PRSSOL, ITEM, COMP   打印BEAM188、BEAM189截面结果

item

comp

截面数据及分量标志

S

COMP

X,XZ,YZ应力分量

PRIN

S1,S2,S3主应力

SINT应力强度,SEQV等效应力

EPTO

COMP

总应变

PRIN

总主应变,应变强度,等效应变

EPPL

COMP

塑性应变分量

PRIN

主塑性应变,塑性应变强度,等效塑性应变

    说明:只有刚计算完还未退出ANSYS时可用,重新进入ANSYS时不可用

    ◆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

应力分量

123

主应力

Int,eqv

应力intensity,等效应力

epeo

x,y,z,xy,yz,xz

总位移分量

1,2,3

主应变

Int,eqv

应变intensity,等效应变

epel

x,y,z,xy,yz,xz

弹性应变分量

123

弹性主应变

Int,eqv

弹性intensity,弹性等效应变

eppl

x,y,z,xy,yz,xz

塑性应变分量

    ◆PRNSOL, item, comp             打印选中节点结果

    item: 项目(见上表)

    comp: 分量

    ◆PRETAB, LAB1, LAB2, ……LAB9  沿线单元长度方向绘单元表数据

    LABn : 空:   所有ETABLE命令指定的列名

    列名: 任何ETABLE命令指定的列名

    ◆PLLS, LABI, LABJ, FACT, KUND  沿线单元长度方向绘单元表数据

    LABI:节点I的单元表列名

    LABJ:节点J的单元表列名

    FACT: 显示比例,缺省为1

    kund: 0 不显示未变形的结构

          1 变形和未变形重叠

          2 变形轮廓和未变形边缘

5/post26 (时间历程后处理)

    ◆nsol, nvar, node, item, comp,name 在时间历程后处理器中定义节点变量的序号

    nvar:变量号(从2到nv(根据numvar定义))

    node: 节点号

item

comp

u

x, y,z

rot

x, y,z

    ◆ESOL, NVAR, ELEM, NODE, ITEM, COMP, NAME 将结果存入变量

    NVAR: 变量号,2以上

    ELEM: 单元号

    NODE: 该单元的节点号,决定存储该单元的哪个量,如果空,则给出平均值

    ITEM:

    COMP:

    NAME: 8字符的变量名, 缺省为ITEM加COMP

    ◆rforce, nvar, node, item, comp, name   指定待存储的节点力数据

    nvar: 变量号

    node: 节点号

item

comp

F

x, y.z

M

x, y,z

    name: 给此变量一个名称,8个字符

◆add, ir, ia,ib,ic,name,--,--,facta, factb, factc  将ia,ib,ic变量相加赋给ir变量

    ir, ia,ib,ic:变量号

    name:     变量的名称

    ◆/grid, key

    key:  “0” 或“off”  无网络

          “1”或“on”  xy网络

          “2”或“x”   只有x线

          “3”或“y”    只有y线

    ◆xvar, n

    n: “0”或“1”  将x轴作为时间轴

       “n”        将x轴表示变量“n”

       “-1”     ?

    ◆/axlab, axis, lab                          定义轴线的标志

    axis: “x”或“y”

    lab:   标志,可长达30个字符

    ◆plvar, nvar, nvar2, ……,nvar10       画出要显示的变量(作为纵坐标)

    ◆prvar, nvar1, ……,nvar6            列出要显示的变量

6PLOTCONTROL菜单命令

    ◆pbc, ilem, ……,key, min, max, abs     在显示屏上显示符号及数值  

    item: u     所加的位移约束

          rot   所加的转角约束

    key:  0    不显示符号

            1       显示符号

           2 显示符号及数值

    ◆/SHOW, FNAME, EXT, VECT, NCPL  确定图形显示的设备及其他参数

    FNAME: X11:屏幕

    文件名:各图形将生成一系列图形文件

    JPEG: 各图形将生成一系列JPEG图形文件

    说明:没必要用此命令,需要的图形文件可计算后再输出

7参数化设计语言

    ◆*do, par, ival, fval, inc             定义一个do循环的开始

    par: 循环控制变量

    ival, fval, inc:起始值,终值,步长(正,负)   
  
    ◆*enddo                         定义一个do循环的结束

    ◆*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) 等效应力

    等效应力  
 
    或   
      
    若
         
    则有        (弹性状态下)

    本站是提供个人知识管理的网络存储空间,所有内容均由用户发布,不代表本站观点。请注意甄别内容中的联系方式、诱导购买等信息,谨防诈骗。如发现有害或侵权内容,请点击一键举报。
    转藏 分享 献花(0

    0条评论

    发表

    请遵守用户 评论公约

    类似文章 更多