1.1.1 注意
1. 尽量加面载,不加集中力,以免奇异点 2. 面的切向荷载必须借助面单元 1.2 求解载荷步
u lssolve, lsmin, lsmax, lsinc 读入并求解多个荷载步 lsmin, lsmax, lsinc :荷载步文件范围 2 /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,uy,uz 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 变形轮廓和未变形边缘
1 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 列出要显示的变量 2 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图形文件 说明:没必要用此命令,需要的图形文件可计算后再输出 1 参数化设计语言 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句 ? 2 理论手册
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) 等效应力 等效应力
或
若
则有 (弹性状态下) |