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,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: ◆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 | 应力分量 | 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 | 塑性应变分量 |
◆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: 节点号 ◆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: 节点号 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) 等效应力 等效应力 或 若 则有 (弹性状态下)
|