ANSYS建模apdl命令流实例应用
ANSYS中的APDL命令
结合自身经验,谈ANSYS中的APDL命令(二)发表时间:2009-5-10 作者: 倪欣来源: e-works关键字: ANSYS APDL 命令流在ANSYS中,命令流是由一条条ANSYS的命令组成的一个命令组合,这些命令按照一定顺序排布,能够完成一定的ANSYS功能,本文是作者结合自身经验所总结的一些命令。
1.1 /prep7(进入前处理)定义几何图形:关键点、线、面、体(1).csys,kcnkcn , 0 迪卡尔坐标系1 柱坐标2 球4 工作平面5 柱坐标系(以Y轴为轴心)n 已定义的局部坐标系(2).numstr, label, value 设置以下项目编号的开始nodeelemkplineareavolu注意:vclear, aclear, lclear, kclear 将自动设置节点、单元开始号为最高号,这时如需要自定义起始号,重发numstr(3).K, npt, x,y,z, 定义关键点Npt:关键点号,如果赋0,则分配给最小号(4).Kgen,itime,Np1,Np2,Ninc,Dx,Dy,Dz,kinc,noelem,imoveItime:拷贝份数Np1,Np2,Ninc:所选关键点Dx,Dy,Dz:偏移坐标Kinc:每份之间节点号增量noelem: “0” 如果附有节点及单元,则一起拷贝。
“1”不拷贝节点和单元imove:“0” 生成拷贝“1”移动原关键点至新位置,并保持号码,此时(itime,kinc,noelem)被忽略注意:MAT,REAL,TYPE 将一起拷贝,不是当前的MAT,REAL,TYPE(5).A, P1, P2, ……… P18 由关键点生成面(6).AL, L1,L2, ……,L10 由线生成面面的法向由L1按右手法则决定,如果L1为负号,则反向。
(线需在某一平面内坐标值固定的面内)(7).vsba, nv, na, sep0,keep1,keep2 用面分体(8).vdele, nv1, nv2, ninc, kswp 删除体kswp: 0 只删除体1 删除体及面、关键点(非公用)(9).vgen, itime, nv1, nv2, ninc, dx, dy, dz, kinc, noelem, imove 移动或拷贝体itime: 份数nv1, nv2, ninc:拷贝对象编号dx, dy, dz :位移增量kinc: 对应关键点号增量noelem,:0:同时拷贝节点及单元1:不拷贝节点及单元imove:0:拷贝体1:移动体(10).cm, cname, entity 定义组元,将几何元素分组形成组元cname: 由字母数字组成的组元名entity: 组元的类型(volu, area, line, kp, elem, node)(11).cmgrp, aname, cname1, ……,cname8 将组元分组形成组元集合aname: 组元集名称cname1……cname8: 已定义的组元或组元集名称1.2 定义几个所关心的节点,以备后处理时调用节点号。
ANSYS计算结果无难事,APDL经典命令让你的模型舞起来
ANSYS计算结果无难事,APDL经典命令让你的模型舞起来1、让你的ANSYS模型'舞'起来ANSYS计算结果的动画可采用ANTIME、ANMODE、ANCNTR、ANHARM等自动生成动画(AVI格式),使结果展示更加生动直观,相信使用ANSYS的都会制作。
然而,几何模型或有限元模型则无动画显示功能,有时为展示模型本身,会从多个角度截取图片。
那么,模型能否也可制作动画呢?答案是肯定的。
利用ANSYS的图形存储命令/SEG可以实现此功能,让你的模型动起来。
具体过程详见命令流中及其注释,动画上传总是失败,自己生成不要观看吧。
Finish$/clear$/prep7!简单的创建几何模型以减少篇幅blc4,0,0,4,2,5cyl4,2,4,1,,2,,4!关闭图例信息/plopts,info,off!以下开始制作模型动画!删除当前储存的图形/seg,dele/seg,multi,jhdh,1 !独立存储且不覆盖,文件名为jhdh/auto,1 !自动计算与图形区合适显示方式!正视/view,1,0,0,1$vplot!侧视/view,1,1$vplot!俯视/view,1,,1$vplot!D视图/view,1,1,1,1$vplot!循环36次,每次改变10度视角*do,i,1,36$/ang,1,10,ys,1$/replot$*enddo!关闭图形存储操作,保存为jhdh.avi文件/seg,off$/anfile,save,jhdh,avi其实比较简单,一旦进入模型动画制作过程,所有的xPLOT(x=KLAVNE)绘制的图形都将进入动画序列,按显示过程形成一部连续的动画。
2、用一个命令解决ANSYS数据列表分页早年初学ANSYS时,经常用到xLIST(如NLIST、ELIST、KLIST、LLIST、ALIST、VLIST等命令)和PRxSOL(如PRNSOL、PRESOL、PRRSOL、PRETAB、PRPATH)等列表命令,并希望将这些内容保存到TXT文件中,然后再导入EXCEL中处理。
ANSYS中的APDL命令总结
在ANSYS中,命令流是由一条条ANSYS的命令组成的一个命令组合,这些命令按照一定顺序排布,能够完成一定的ANSYS功能,这些功能一般来说通过菜单操作也能够实现(而那些命令流能够实现,菜单操作实现不了的单个命令比较少见)。
以下命令是结合我自身经验,和前辈们的一些经验而总结出来的,希望对大家有帮助。
(1).Lsel, type, item, comp, vmin, vmax, vinc, kswp 选择线type: s 从全部线中选一组线r 从当前选中线中选一组线a 再选一部线附加给当前选中组aunoneu(unselect)inve: 反向选择item: line 线号loc 坐标length 线长comp: x,y,zkswp: 0 只选线1 选择线及相关关键点、节点和单元(2).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”仅用绝对值(3).Esel, type, item, comp, vmin, vmax, vinc, kabs 选择一组单元type: S: 选择一组单元(缺省)R: 在当前组中再选一部分作为一组A: 为当前组附加单元U: 在当前组中不选一部分单元All: 选所有单元None: 全不选Inve: 反向选择当前组Stat: 显示当前选择状态Item:Elem: 单元号Type: 单元类型号Mat: 材料号Real: 实常数号Esys: 单元坐标系号(4). 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次项的系数(5). 定义DP材料:首先要定义EX和泊松比:MP,EX,MAT,……MP,NUXY,MAT,……定义DP材料单元表(这里不考虑温度):TB,DP,MAT进入单元表并编辑添加单元表:TBDATA,1,CTBDATA,2,ψTBDATA,3,……如定义:EX=1E8,NUXY=0.3,C=27,ψ=45的命令如下:MP,EX,1,1E8MP,NUXY,1,0.3TB,DP,1TBDATA,1,27TBDATA,2,45这里要注意的是,在前处理的最初,要将角度单位转化到“度”,即命令:*afun,deg(6). 根据需要耦合某些节点自由度cp, nset, lab,,node1,node2,……node17nset: 耦合组编号lab: ux,uy,uz,rotx,roty,rotz ,allnode1-node17: 待耦合的节点号。
ANSYS APDL命令流详解-3坐标系和工作平面
3.节点坐标系的旋转与修改
⑴ 将某些节点的坐标系旋转到与当前激活坐标系(简 称“当前坐标系”)方向一致 命令:NROTAT, NODE1, NODE2, NINC
其中NODE1、NODE2、NINC ---要旋转节点的起始号、末编号 (缺省为NODE1)及递增值(缺省值为1)。如NODE1=ALL则 其后参数将被忽略,NODE1也可为元件名。
⑸ 根据激活的坐标系定义局部坐标系 命令:CLOCAL, KCN, KCS, XL, YL, ZL, THXY, THYZ, THZX, PAR1,
PAR2
⑹ 删除局部坐标系 命令:CSDELE, KCN1, KCN2, KCINC 其中:
KCN1---为要删除的局部坐标系的起始编号,如果KCN1=ALL,则其后参 数将忽略。
例如:NMODIF,8,,,,15--- 修改节点8的节点坐标系方向,使之绕Z轴旋转 15°。
⑶在创建节点时直接定义其坐标系的旋转角度 命令:N, NODE, X, Y, Z, THXY, THYZ, THZX
例如:N,4,1,2,4,10,15,30---表示新建4号节点在当前坐标系中的坐标为 1,2,4,其节点坐标系绕Z,X,Y轴的角度分为10°、 15°和 30°。
KCN2---为要删除的局部坐标系的最终编号。 KCINC---为编号的递增数值,缺省为1。
例如:CSDELE,ALL---则删除了所有的局部坐标系。 CSDELE,11,15,2---则删除了11、13、15号局部坐标系。
⑺ 查看激活坐标系和局部坐标系 命令:CSLIST, KCN1, KCN2, KCINC 例如:CSLIST,ALL---则列表显示所有坐标系,并列出相关信息。
4. 单元坐标系的定义与修改
ANSYS_数据文件读写的APDL命令详解及实例
ANSYS 数据文件读写的APDL命令详解及实例作者:huright一 FORTRAN数据格式I格式(又叫整数格式)一般形式:Iw 或:Iw.m其中:w 一个数据占的位数宽度(又称“字段宽度”),m 需要输出的最少数字位数。
例1:(1)数字在指定的区域内向右端靠齐,如果数字位数比指定的字段宽度w小,则左边补以空格。
负数的符号也包含在字段宽度内。
(2)如果数字的位数超过了规定的字段宽度w,则不输出有效数据,而在该字段宽度范围内充满“*”符号。
(3)如果数字的位数超过了m,则按实际应输出的位数输出(但条件是不能超过w)。
m 不包括负号所占的一列。
F格式(又叫小数型格式)一般形式:Fw.dw 各数值占的总位数 d 输出数据的小数位数(小数点后的位数)。
例1:(1)数字在指定的区域内向右端靠齐,如果数字位数(含小数点和符号位)比指定的字段宽度w小,则左边补以空格;如果数字的位数超过了规定的字段宽度w,则不输出有效数据,而在该字段宽度范围内充满“*”符号。
(2)如果数据的小数位数比指定的小数位数d小,则在小数右边补0以凑足d位;如果小数位数大于d位,则输出时多于的小数位数按“四舍五入”规则舍去。
(3)假设b为数据整数部分的位数,则应使w≥b+d+1(小数点占一列),如果输出负数,则应保证w≥b+d+2(小数点和负号各占一列)。
(4)用F格式输出时应注意,由于难以事先确切估计出数据的大小,输出大的数时容易产生“宽度不够”的错误(由于w不够大),输出小的数时会出现丢掉有用数字的情况(由于d不够大而将后面的数字截去),这就是“大数印错,小数印丢”。
E格式(又叫指数型格式)一般形式:w.dw 各数值占的总位数,d 输出数据的小数位数(小数点后的位数)。
例1:(1)采取标准化的指数形式输出一个实数,d为以指数形式出现的数据的数字部分的小数位数。
(2)指数部分一般占4列,其中字母“E”和指数的符号各占一列,指数2列。
ANSYS_数据文件读写的APDL命令详解及实例
ANSYS 数据文件读写的APDL命令详解及实例作者:huright一 FORTRAN数据格式I格式(又叫整数格式)一般形式:Iw 或:Iw.m其中:w 一个数据占的位数宽度(又称“字段宽度”),m 需要输出的最少数字位数。
例1:(1)数字在指定的区域内向右端靠齐,如果数字位数比指定的字段宽度w小,则左边补以空格。
负数的符号也包含在字段宽度内。
(2)如果数字的位数超过了规定的字段宽度w,则不输出有效数据,而在该字段宽度范围内充满“*”符号。
(3)如果数字的位数超过了m,则按实际应输出的位数输出(但条件是不能超过w)。
m 不包括负号所占的一列。
F格式(又叫小数型格式)一般形式:Fw.dw 各数值占的总位数 d 输出数据的小数位数(小数点后的位数)。
例1:(1)数字在指定的区域内向右端靠齐,如果数字位数(含小数点和符号位)比指定的字段宽度w小,则左边补以空格;如果数字的位数超过了规定的字段宽度w,则不输出有效数据,而在该字段宽度范围内充满“*”符号。
(2)如果数据的小数位数比指定的小数位数d小,则在小数右边补0以凑足d位;如果小数位数大于d位,则输出时多于的小数位数按“四舍五入”规则舍去。
(3)假设b为数据整数部分的位数,则应使w≥b+d+1(小数点占一列),如果输出负数,则应保证w≥b+d+2(小数点和负号各占一列)。
(4)用F格式输出时应注意,由于难以事先确切估计出数据的大小,输出大的数时容易产生“宽度不够”的错误(由于w不够大),输出小的数时会出现丢掉有用数字的情况(由于d不够大而将后面的数字截去),这就是“大数印错,小数印丢”。
E格式(又叫指数型格式)一般形式:w.dw 各数值占的总位数,d 输出数据的小数位数(小数点后的位数)。
例1:(1)采取标准化的指数形式输出一个实数,d为以指数形式出现的数据的数字部分的小数位数。
(2)指数部分一般占4列,其中字母“E”和指数的符号各占一列,指数2列。
ANSYS APDL命令流详解-5几何模型的布尔运算
!进入前处理 !利用DO循环创建关键点 !如果I能被2整除则执行下面命令 !创建坐标为(2*i,4)的关键点 !否则(I不能被2整除) !创建坐标为(2*i,-4)的关键点 !结束IF语句 !结束循环语句 !利用循环创建线
l,2,19$l,1,20
LINL,all
!作线交运算,由于没有公共部分不能运算
PTOL=DEFA时,则恢复缺省设置; PTOL=STAT时,则列表输出当前设置。
2.3.2 交运算 Intersection
交运算就是由图素的共同部分形成一个新的图素,其运算 结果只保留两个或多个图素的重叠部分。
交运算分为公共相交和两两相交两种。公共相交就是仅保 留所有图素的重叠部分,即只生成一个图素,当图素很多时可 能不存在公共部分,这时布尔运算不能完成。两两相交是保留 任意两个图素的公共部分,有可能生成很多图素。
Value=1布尔操作失败时不产生一个警告信息。
Value=-1布尔操作失败时产生一个错误信息。
当Lab=VERSION时:Value=RV52(缺省)激活5.2版本兼容性选项
Value=RV51激活5.1版本兼容性选项。
该命令的全部缺省设置是操作失败产生一个警告信息,删除输入
图 素,不输出编号警告信息,使用5.2版本布尔兼容性选项。该命令可多
其中NX1~NX9为相交X的编号,NX1可以为P、ALL或元件名(其中X 表示L、A或V)。
4. 交运算的命令流示例
⑴ 线相交
任意创建一组线,分别作交运算和两两相交运算,命令流如下
/prep7 *do,i,1,20 *if,mod(i,2),eq,0,then k,i,2*i,4 *else k,i,2*i,-4 *endif *enddo *do,i,1,19$l,i,i+1$*enddo
ANSYS APDL命令流详解-13通用与时间历程后处理技术
5.1 通用后处理—图形显示结果
命令
功能
PLDISP 显示结构变形图
PLNSOL 显示节点结果图
PLESOL 显示单元结果图
PLVECT
PLCRACK
以矢量方式显示 结 果图
显示裂缝或压碎图
备注
结果内容很多,根据需要选用Item及 Comp 结果内容很多,根据需要选用Item及 Comp 如U、S等
Lab2,Lab3---用户定义的分矢量标识符,如Item为预定标识符,必须为空。 LabP---合成矢量标识符,缺省为Item。 Mode---显示方式控制。
如为空,则采用/DEVICE中的KEY参数指定的方式; 如MODE=RAST则采用光栅模式; 如MODE=VECT则采用矢量模式显示。 Loc---显示单元场结果的矢量位置。
5.1 通用后处理—结果输出控制选项
用于图形显示和列表显示,如导出结果的方式和显示比例等设置。
命令
功能
备注
AVPRIN
定义矢量和主轴的计算方法
用于计算主应力或主应变等时
AVRES
定义结果数据平均处理
仅适用于PowerGraphics模式
/EFACET
设置单元每边的分段数目
见2.4.1中
/DSCALE
第5章 通用与时间历程后处理技术
ANSYS两个后处理器: ★通用后处理器POST1:
查看整个模型在各个时间点上的结果. ★时间历程后处理器POST26:
查看整个模型上的某一点结果随时间的变化曲线。 后处理可在求解完后直接进入,也可在重新进入 ANSYS后读入文件进入后处理。
5.1 通用后处理—读入结果文件
ANSYS 分析的命令流(APDL语言)
ANSYS 分析的命令流(APDL语言)前面分析的过程都是基于ANSYS用户图形界面(GUI命令流)形式,GUI 形式非常的直观明了,分析过程中通常会采用此种的分析方式。
但是这种分析方式会有一个缺点,就是当我们操作失误时候,没有后退的功能,因此我们不得不重新操作,这样就会给分析带来很多的麻烦。
但如果以命令流的形式进行计算的话就能很容易的减少这些麻烦。
下面就是本次点机前处理和分析计算过程的命令流/UNITS,SI !采用国际单位制/TITLE,2D DIANJI Static Analysis !定义分析名称KEYW,MAGNOD,1*AFUN,DEG !指定角度单位为度/PREP7 !进入前处理器ET,1,PLANE53 !设定单元类型为plane53EMUNIT,MKS !电磁单位LOCAL,12,0,0,0,0,90 !定义局部坐标系12--17LOCAL,13,0,0,0,0,30LOCAL,14,0,0,0,0,330LOCAL,15,0,0,0,0,270LOCAL,16,0,0,0,0,210LOCAL,17,0,0,0,0,150MP,MURX,1,1.0 ! 定义第一种材料的相对磁导率MP,MURX,2,1.0 ! 定义第二种材料的相对磁导率TB,BH,3TBPT,,35.03,0.1 ! 第三种材料的B-H磁化特性TBPT,,46.97,0.2TBPT,,57.32,0.3TBPT,,66.08,0.4TBPT,,74.04,0.5TBPT,,82.01,0.6TBPT,,91.56,0.7TBPT,,103.5,0.8TBPT,,117.83,0.9TBPT,,136.15,1TBPT,,160.83,1.1TBPT,,203.03,1.2TBPT,,286.62,1.3TBPT,,461.78,1.4TBPT,,955.41,1.5TBPT,,2547.8,1.6!TBPLOT,BH,3 !绘制B-H曲线TB,BH,4TBPT,,130,0.1 ! 定义第四种材料的B-H磁化特性TBPT,,170,0.2TBPT,,197,0.3TBPT,,218,0.4TBPT,,250,0.5TBPT,,290,0.6TBPT,,338,0.7TBPT,,400,0.8TBPT,,472,0.9TBPT,,570,1TBPT,,682,1.1TBPT,,810,1.2TBPT,,970,1.3TBPT,,1600,1.4TBPT,,2520,1.5TBPT,,3520,1.6!TBPLOT,BH,4 !绘制B-H曲线MP,MURX,5,1.17 ! 定义第五种材料的相对磁导率MP,MGXX,5,4.4E+005 ! 定义矫顽力矢量为X方向分量MP,MURX,6,1.17 ! 定义第六种材料的相对磁导率MP,MGXX,6,4.4E+005 ! 定义矫顽力矢量为X方向分量MP,MURX,7,1.17 ! 定义第七种材料的相对磁导率MP,MGXX,7,4.4E+005 ! 定义矫顽力矢量为X方向分量MP,MURX,8,1.17 ! 定义第八种材料的相对磁导率MP,MGXX,8,4.4E+005 ! 定义矫顽力矢量为X方向分量MP,MURX,9,1.17 ! 定义第九种材料的相对磁导率MP,MGXX,9,4.4E+005 ! 定义矫顽力矢量为X方向分量MP,MURX,10,1.17 ! 定义第十种材料的相对磁导率MP,MGXX,10,4.4E+005 ! 定义矫顽力矢量为X方向分量RECTNG,-0.00275,0.00275,0.0485,0.062 ! 创建一个矩形RECTNG,-0.00025,0.00025,0.048,0.0485 ! 再创建一个矩形AADD,ALL ! 将两个矩形融合在一起,形成槽结构CSYS,1 !转动工作平面到柱坐标下AGEN,36,3,38,1,,10,0 ! 复制槽/PNUM,AREA,ALL ! 给所有槽编号PCIRC,0.048,0.074 ! 创建定子圆环AOVLAP,ALL !对定子圆环和槽进行交迭布尔操作使得在连接线上共节点PCIRC,0.037,0.047,78,102 ! 创建永磁体极CSYS,1 !转到柱坐标下AGEN,6,3,8,1,,60,0 ! 复制6个磁极PCIRC,0.023,0.047 ! 创建转子圆环AOVLAP,3,10,45,51,57,63,69 ! 对转子圆环与磁极进行交迭操作PCIRC,0.047,0.0475 !创建两层气隙PCIRC,0.0475,0.048AOVLAP,ALLAGLUE,ALL ! 将所有交接面黏在一起NUMCMP,AREA !压缩面编号ASEL,S,AREA,,1,35 !选中槽ASEL,A,AREA,,42AATT,1 !材料定义为1号材料ASEL,S,AREA,,45ASEL,A,AREA,,43 !选中气隙AATT,2 !材料定义为2号材料ASEL,S,AREA,,46 !选中定子铁芯AATT,3 !材料定义为3号材料ASEL,S,AREA,,44 !选中转子磁轭AATT,4 !材料定义为4号材料ASEL,S,AREA,,40 !选中第1-6个永磁体磁极AATT,5,,,12 !材料定义为5-10号材料ASEL,S,AREA,,39AATT,6,,,13ASEL,S,AREA,,38AATT,7,,,14ASEL,S,AREA,,41AATT,8,,,15ASEL,S,AREA,,37AATT,9,,,16ASEL,S,AREA,,36AATT,10,,,17/PNUM,MAT,1 !打开材料编号APLOT !重新显示ALLSEL,ALL !选择所有实体SMRTSIZE,1 !指定智能网格划分等级L AMESH,ALL !智能网格划分!ESEL,S,MAT,,1!CM,ARM,ELEM!FMAGBC,'ARM'!SA VE,DIANJI_2D_MESH.db !保存!SA VEFINISH !完成网格划分/SOLU !进入求解器ANTYPE,STATIC !选择静态磁场分析NROPT,AUTO !采用牛顿---拉夫森迭代方法!ASEL,S,AREA,,1,35 !选中槽!ASEL,A,AREA,,42!BFE,ALL,JS,1,,,15000000 !施加电流密度LOCAL,11,1 ! 定义柱坐标系FLST,2,8,4,ORDE,4FITEM,2,289FITEM,2,-292FITEM,2,561FITEM,2,-564DL,P51X, ,AZ,0, !加边界条件ALLSEL,ALLMAGSOLV!SA VE,EMAGE_2D_MESH.db !保存!SA VEFINISH/POST1PLF2D,27,0,10,1 !显示磁力线图FMAGSUM,'ARM' !对电磁力求和plvect,B,,,,VECT,ELEM,ON,0 !显示磁感应强度矢量plvect,H,,,,VECT,ELEM,ON,0 !显示磁场强度矢量/GRAPHICS,POWERA VRES,2PLNSOL,B,SUM !显示磁通密度等值云图flishESEL,S,MAT,,2CM,ARM,ELEM!生成一个组件FMAGSUM,'ARM'!以下命令流是保存图片时使用,使其背景为白,结合PlotCtrlsHard ﹥Copy命令可设置格式/GRAPHICS,POWER/RGB,INDEX,100,100,100, 0/RGB,INDEX, 80, 80, 80,13/RGB,INDEX, 60, 60, 60,14/RGB,INDEX, 0, 0, 0,15/REPLOT !背景反白显示/SHOW,JPEG,,0JPEG,QUAL,75,JPEG , ORIENT, HORIZJPEG , Color, 2JPEG , TMOD, 1/GFILE,600,EPLOT/SHOW,CLOSE/DEVICE,VECTOR,0。
APDL与命令流的使用
数组参数必须先定义后使用。
参数名不分大小写, 即, “RAD” 和 “Rad” 是相同的。所有参
数 内部都是以大写形式存储的。
中南大学
工程软件应用
2.1标量参数
• 标量参数的使用应注意:
– 必须以字母开头,只能包含字母、数字和下划线;
– 不能以下划线开头,以下划线开头/结尾的参数为系统隐
•
•
•
•
其中:Par——将被赋值的参数名;
Entity——将提取对象关键字(NODE,ELEM,KP……);
ENTNUM——提取对象的号码(若为0指全部实体);
Item1,IT1NUM——是指对象的项目名和项目号;
•
Item2,IT2NUM——同上。
中南大学
工程软件应用
菜单路径:
Utility Menu>Parameters>Get Scalar Data
– *GET提取命令。
– 内嵌提取函数。
*GET提取命令:
从特定对象(点、单元或面等)中提取数据库数据,并赋予某个用户命令
的Scalar参数,有两种途径:
– 菜单路径:Utility Menu>Parameters>Get Scalar Data。
– *GET命令的使用格式。
– *GET,Par,Entity,ENTNUM,Item1,IT1NUM,Item2,IT2NUM
-*CFOPEN,Fname,Ext,Dir,Loc
……
*CFWRITE,Command
方法3-通过菜单创建:
Utility Menu>Macro>Create Macro
方法4
-通过文本编辑器进行编辑。这是最常用的方法
ANSYS APDL命令流详解-9网格划分技术
面映射网格划分/体映射网格划分 扫掠生成体网格/单元有效性检查/网格修改
3.4 网格划分实例
基本模型的网格划分/复杂面模型的网格划分 复杂体模型的网格划分
创建几何模型后,必须生成有限元模型才能 分析计算,生成有限元模型的方法就是对几何 模型进行网格划分,网格划分主要过程包括三 个步骤: ⑴ 定义单元属性
3.1 定义单元属性 —材料属性
每一组材料属性有一个材料参考号,用于识别各个材料特
性组。一个模型中可有多种材料特性组。 1. 定义线性材料属性 命令:MP,Lab,MAT,C0,C1,C2,C3,C4 Lab---材料性能标识,其值可取:
EX:弹性模量(也可为EY、EZ)。 ALPX:线膨胀系数(也可为ALPY、ALPZ)。 PRXY:主泊松比(也可为PRYZ、PRXZ)。 NUXY:次泊松比(也可为NUYZ、NUXZ)。 GXY:剪切模量(也可为GYZ、GXZ)。 DAMP:用于阻尼的K矩阵乘子,即阻尼比。 DMPR:均质材料阻尼系数。 MU:摩擦系数 。 DENS:质量密度。
与下列两条命令等效:
et,1,beam4,,1
et,3,beam189,1
3.1 定义单元属性 —单元类型
3. 自由度集 命令:DOF,Lab1,Lab2,Lab3,Lab4,Lab5,Lab6,Lab7,Lab8,Lab9,Lab10
4.改变单元类型 命 令:ETCHG,Cnv
5. 单元类型的删除与列表 删除命令:ETDELE,ITYP1,ITYP2,INC 列表命令:ETLIST,ITYP1,ITYP2,INC
单元属性包括单元类型、实常数、材料特性、单元坐标系和截 面号等。
⑵ 定义网格控制选项
★对几何图素边界划分网格的大小和数目进行设置; ★没有固定的网格密度可供参考; ★可通过评估结果来评价网格的密度是否合理。
ANSYS APDL命令流详解-7几何建模技巧
⑷ 改变当前工作目录 GUI:Utility Menu>File>Change Directory 命令:/CWD, DIRPATH 其中DIRPATH为新工作目录的全路径名。当指定的新工作路 径不存在时,则不会改变路径,且给出错误信息。该命令可以 使用系统认可的任何目录,包括中文命名的目录。
⑸ 指定主标题 GUI:Utility Menu>File>Change Title 命令:/TITLE, Title 其中Title为主标题,最多72个字符,用%将参数或表达式括起 来也可进行替换。该主标题可显示在屏幕上的图形区,还可用 /STITLE指定子标题。
LOG和ERR文件,而不另外再新建文件; 如为1或ON则创建与工作文件名同名的LOG和ERR 文件,但原来的LOG和ERR文件不删除。
该命令执行后,其后所有要生成的文件都会以新的文件名命
名。工作文件名不能为中文,可以为数字、字母及特殊符号等
组成;一旦不接受所给出的工作文件名,ANSYS会直接采用缺 省工作文件名。
NP1,NP2,NINC---为进行缩放的关键点的编号范围。其中NP1也可为ALL或组件名。 NL1,NL2,NINC---为进行缩放的线的编号范围。其中NL1也可为ALL或组件名。 NA1,NA2,NINC---为进行缩放的线的编号范围。其中NA1也可为ALL或组件名。 NV1,NV2,NINC---为进行缩放的线的编号范围。其中NV1也可为ALL或组件名。 RX,RY,RZ---在当前坐标系下,关键点的X、Y、Z方向坐标值的比例因子。
该 命 令 退 出 ANSYS的 四 大 处 理 模 块 或 图 形 显 示 , 回 到 Begin层。在各个处理模块之间使用该命令可保证数据库的完 整,但数据库不会自动保存,需要执行Save命令或其它命令。
APDL命令流:将ansys分析结果输出为tecplot格式
APDL命令流:将ansys分析结果输出为tecplot格式/post1file,Jobname,rst !指明从哪⼀个结果⽂件中读取数据*get,NodeNum,NODE,0,COUNT !得到模型的所有节点数*get,nd,NODE,0,NUM,MIN !得到模型的最⼩节点编号*dim,nodes,array,NodeNum !定义⼀个存储节点的数组*dim,xyz,array,NodeNum,3 !定义⼀个存储节点坐标的数组*dim,NodeTemp,array,NodeNum !定义⼀个存储节点的数组set,4,1 !读⼊第四载荷步、第⼀个⼦步的结果*do,i,1,NodeNum,1nodes(i)=nd !将节点编号存储在nodes数组中xyz(i,1)=NX(nd) !在xyz(i,1)中存储节点的X坐标xyz(i,2)=NY(nd) !在xyz(i,1)中存储节点的Y坐标xyz(i,3)=NZ(nd) !在xyz(i,1)中存储节点的Z坐标*Get,NodeTemp(i),NODE,nodes(i),TEMPnd=NDNEXT(nd) !取得下⼀个节点编号*enddo*get,ElementNum,ELEM,0,COUNT !得到模型的所有单元数*get,el,ELEM,0,NUM,MIN !得到模型的最⼩单元编号*dim,elems,array,ElementNum !定义⼀个存储单元的数组*dim,ndlst,array,ElementNum,8 !定义⼀个存储节点坐标的数组*do,i,1,ElementNum,1elems(i)=el !将单元编号存储在elems数组中*do,j,1,8,1ndlst(i,j)=NELEM(el,j) !将单元对应的节点编号存储在ndlst数组中*enddoel=ELNEXT(el) !取得下⼀个单元编号*enddo*cfopen,mesh,dat !输出的⽂件名为:mesh.dat!写TecPlot⽂件的⽂件头*vwrite(’TITLE=Ansys Temperature Analysis’)!写变量名*vwrite(’VARIABLES=”X”,”Y”,”Z”,”Temp”‘)!写节点数和单元数*vwrite,NodeNum,ElementNum(’ZONE n=’f6.1,’ e=’f6.1,’ f=fepoint’,' et=brick’)!写节点的坐标和温度值*vwrite,xyz(1,1),xyz(1,2),xyz(1,3),NodeTemp(1)(4f12.6)!写⼀个空⾏*vwrite(’ ‘)!写各单元对应的节点号*vwrite,ndlst(1,1),ndlst(1,2),ndlst(1,3),ndlst(1,4),ndlst(1,5),ndlst(1,6),ndlst(1,7),ndlst(1,8)(8f6.0)!关闭⽂件*CFCLOSEfinish。
ANSYS 数据文件读写的APDL命令详解及实例
ANSYS 数据文件读写的APDL命令详解及实例源自机知网作者:huright一 FORTRAN数据格式I格式(又叫整数格式)一般形式:Iw 或:Iw.m其中:w 一个数据占的位数宽度(又称“字段宽度”),m 需要输出的最少数字位数。
例1:(1)数字在指定的区域内向右端靠齐,如果数字位数比指定的字段宽度w小,则左边补以空格。
负数的符号也包含在字段宽度内。
(2)如果数字的位数超过了规定的字段宽度w,则不输出有效数据,而在该字段宽度范围内充满“*”符号。
(3)如果数字的位数超过了m,则按实际应输出的位数输出(但条件是不能超过w)。
m 不包括负号所占的一列。
F格式(又叫小数型格式)一般形式:Fw.dw 各数值占的总位数 d 输出数据的小数位数(小数点后的位数)。
例1:(1)数字在指定的区域内向右端靠齐,如果数字位数(含小数点和符号位)比指定的字段宽度w小,则左边补以空格;如果数字的位数超过了规定的字段宽度w,则不输出有效数据,而在该字段宽度范围内充满“*”符号。
(2)如果数据的小数位数比指定的小数位数d小,则在小数右边补0以凑足d位;如果小数位数大于d位,则输出时多于的小数位数按“四舍五入”规则舍去。
(3)假设b为数据整数部分的位数,则应使w≥b+d+1(小数点占一列),如果输出负数,则应保证w≥b+d+2(小数点和负号各占一列)。
(4)用F格式输出时应注意,由于难以事先确切估计出数据的大小,输出大的数时容易产生“宽度不够”的错误(由于w不够大),输出小的数时会出现丢掉有用数字的情况(由于d不够大而将后面的数字截去),这就是“大数印错,小数印丢”。
E格式(又叫指数型格式)一般形式:w.dw 各数值占的总位数,d 输出数据的小数位数(小数点后的位数)。
例1:(1)采取标准化的指数形式输出一个实数,d为以指数形式出现的数据的数字部分的小数位数。
(2)指数部分一般占4列,其中字母“E”和指数的符号各占一列,指数2列。
ANSYS房屋框架结构模型命令流APDL
ANSYS房屋框架结构模型命令流APDL!********************提示********************************************!第一步:建立框架结构模型!第二步:框架结构施加重力荷载!第三步:框架结构施加活荷载!第四步:框架结构施加风荷载!第五步:框架结构荷载组合运算!第六步:框架结构模态分析!第七步:框架结构屈曲特征值分析!第八步:框架结构地震时程弹性分析!第九步:结束!采用地震波时间文件为:TIME.TXT!采用地震波文件分别为:AC_X.TXT/AC_Y.TXTFINISH/CLEAR/FILENAME,FRAME_SHEAR_WALL/TITLE, FRAME_SHEAR_WALL!************************************************************** ************ !************************第一步******************************************** !********************框架结构建模型***************************************/PREP7!采用单位为N/m/S 国际单位制!定义单元类型ET,1,BEAM4ET,2,SHELL63!定义实常数R,1,0.25,0.0052,0.0052,0.5,0.5R,2,0.2025,0.0034,0.0034,0.45,0.45 R,3,0.06,0.0002,0.00045,0.3,0.2 R,4,0.25,0.25,0.25,0.25R,5,0.1,0.1,0.1,0.1!定义材料参数MP,EX,1,3.0E10MP,NUXY,1,0.2MP,DENS,1,2500!关键点*DO,II,1,1K, 2+(II-1)*20, 6, 0, (II-1)*6K, 3+(II-1)*20, 12, 0, (II-1)*6K, 4+(II-1)*20, 17, 0, (II-1)*6K, 5+(II-1)*20, 0, 5, (II-1)*6K, 6+(II-1)*20, 6, 5, (II-1)*6K, 7+(II-1)*20, 12, 5, (II-1)*6K, 8+(II-1)*20, 17, 5, (II-1)*6K, 9+(II-1)*20, 0, 10, (II-1)*6K, 10+(II-1)*20, 6, 10, (II-1)*6K, 11+(II-1)*20, 12, 10, (II-1)*6 K, 12+(II-1)*20, 17, 10, (II-1)*6 K, 14+(II-1)*20, 6, 1.875, (II-1)*6 K, 16+(II-1)*20, 6, 3.125, (II-1)*6 K, 18+(II-1)*20, 6, 1.875, 3.9K, 20+(II-1)*20, 6, 3.125, 3.9*ENDDO*DO,II,1,10K, 1+II*20, 0, 0, (II-1)*3+6K, 2+II*20, 6, 0, (II-1)*3+6K, 3+II*20, 12, 0, (II-1)*3+6K, 4+II*20, 17, 0, (II-1)*3+6K, 5+II*20, 0, 5, (II-1)*3+6K, 6+II*20, 6, 5, (II-1)*3+6K, 7+II*20, 12, 5, (II-1)*3+6K, 8+II*20, 17, 5, (II-1)*3+6K, 9+II*20, 0, 10, (II-1)*3+6K, 10+II*20, 6, 10, (II-1)*3+6K, 11+II*20, 12, 10, (II-1)*3+6K, 12+II*20, 17, 10, (II-1)*3+6K, 13+II*20, 0, 1.875, (II-1)*3+6K, 14+II*20, 6, 1.875, (II-1)*3+6K, 15+II*20, 0, 3.125, (II-1)*3+6K, 16+II*20, 6, 3.125, (II-1)*3+6K, 17+II*20, 0, 1.875, (II-1)*3+6+1.5 K, 18+II*20, 6, 1.875, (II-1)*3+6+2 K, 19+II*20, 0, 3.125, (II-1)*3+6+1.5 K, 20+II*20, 6, 3.125, (II-1)*3+6+2 *ENDDO*DO,II,11,11K, 1+II*20, 0, 0, (II-1)*3+6K, 2+II*20, 6, 0, (II-1)*3+6K, 3+II*20, 12, 0, (II-1)*3+6K, 4+II*20, 17, 0, (II-1)*3+6K, 6+II*20, 6, 5, (II-1)*3+6K, 7+II*20, 12, 5, (II-1)*3+6K, 8+II*20, 17, 5, (II-1)*3+6K, 9+II*20, 0, 10, (II-1)*3+6K, 10+II*20, 6, 10, (II-1)*3+6K, 11+II*20, 12, 10, (II-1)*3+6K, 12+II*20, 17, 10, (II-1)*3+6K, 13+II*20, 0, 1.875, (II-1)*3+6K, 14+II*20, 6, 1.875, (II-1)*3+6K, 15+II*20, 0, 3.125, (II-1)*3+6K, 16+II*20, 6, 3.125, (II-1)*3+6K, 17+II*20, 0, 1.875, (II-1)*3+6+1.8 K, 18+II*20, 6, 1.875, (II-1)*3+6+2.5 K, 19+II*20, 0, 3.125, (II-1)*3+6+1.8 K, 20+II*20, 6, 3.125, (II-1)*3+6+2.5 *ENDDO*DO,II,12,12K, 1+II*20, 0, 0, 39.5K, 2+II*20, 6, 0, 39.5K, 3+II*20, 12, 0, 39.5K, 4+II*20, 17, 0, 39.5K, 5+II*20, 0, 5, 39.5K, 6+II*20, 6, 5, 39.5K, 7+II*20, 12, 5, 39.5K, 8+II*20, 17, 5, 39.5K, 9+II*20, 0, 10, 39.5K, 10+II*20, 6, 10, 39.5K, 11+II*20, 12, 10, 39.5K, 12+II*20, 17, 10, 39.5*ENDDO*DO,II,2,13L, (II-1)*20+1, (II-1)*20+2L, (II-1)*20+2, (II-1)*20+3L, (II-1)*20+3, (II-1)*20+4L, (II-1)*20+5, (II-1)*20+6L, (II-1)*20+6, (II-1)*20+7L, (II-1)*20+7, (II-1)*20+8L, (II-1)*20+9, (II-1)*20+10L, (II-1)*20+10, (II-1)*20+11L, (II-1)*20+11, (II-1)*20+12L, (II-1)*20+5, (II-1)*20+9L, (II-1)*20+6, (II-1)*20+10 L, (II-1)*20+3, (II-1)*20+7 L, (II-1)*20+4, (II-1)*20+8 L, (II-1)*20+8, (II-1)*20+12 *ENDDO*DO,II,1,12L, (II-1)*20+1, II*20+1L, (II-1)*20+2, II*20+2L, (II-1)*20+3, II*20+3L, (II-1)*20+4, II*20+4L, (II-1)*20+5, II*20+5L, (II-1)*20+6, II*20+6L, (II-1)*20+7, II*20+7L, (II-1)*20+8, II*20+8L, (II-1)*20+9, II*20+9L, (II-1)*20+10, II*20+10 L, (II-1)*20+11, II*20+11 L, (II-1)*20+12, II*20+12*ENDDO*DO,II,1,1L, (II-1)*20+14, (II-1)*20+18 L, (II-1)*20+16, (II-1)*20+20 *ENDDO*DO,II,2,12L, (II-1)*20+13, (II-1)*20+17 L, (II-1)*20+15, (II-1)*20+19 L, (II-1)*20+14, (II-1)*20+18 L, (II-1)*20+16, (II-1)*20+20 *ENDDOA, 1, 21, 25, 5*DO,II,2,12A, (II-1)*20+1, II*20+1, II*20+5, (II-1)*20+5, (II-1)*20+15, (II-1)*20+19, (II-1)*20+17, (II-1)*20+13 *ENDDO*DO,II,1,12A, (II-1)*20+2, II*20+2, II*20+6, (II-1)*20+6, (II-1)*20+16, (II-1)*20+20, (II-1)*20+18, (II-1)*20+14 *ENDDO*DO,II,1,12A, II*20+1, II*20+2, II*20+6, II*20+5A, II*20+2, II*20+3, II*20+7, II*20+6A, II*20+5, II*20+6, II*20+10,II*20+9A, II*20+6, II*20+7, II*20+11,II*20+10A, II*20+7, II*20+8, II*20+12,II*20+11*ENDDO!立柱网格划分LSEL,S,,,181,192 !底层立柱LATT,1,1,1LESIZE,ALL,,,10LMESH,ALLLSEL,S,,,192,324 !2~12层立柱LATT,1,2,1LESIZE,ALL,,,5LMESH,ALL!梁网格划分LSEL,S,,,1,180LATT,1,3,1,LESIZE,ALL,,,8LMESH,ALL!剪力墙网格划分LSEL,S,,,372 !外剪力墙两端、内剪力墙顶端LSEL,A,,,413LSEL,A,,,461LESIZE,ALL,,,8LSEL,s,,,418 !内剪力墙底端LSEL,A,,,420LESIZE,ALL,,,3LSEL,s,,,325 !内剪力墙底端内侧LSEL,A,,,326LESIZE,ALL,,,6LSEL,S,,,375 !剪力墙空洞顶部LSEL,A,,,379LSEL,A,,,383LSEL,A,,,387LSEL,A,,,391LSEL,A,,,395LSEL,A,,,399LSEL,A,,,403LSEL,A,,,407LSEL,A,,,415LSEL,A,,,419LSEL,A,,,423LSEL,A,,,427LSEL,A,,,431LSEL,A,,,435LSEL,A,,,439LSEL,A,,,443LSEL,A,,,447LSEL,A,,,451LSEL,A,,,455LSEL,A,,,459LSEL,A,,,463LESIZE,ALL,,,2!剪力墙网格划分ASEL,S,,,1,24AATT,1,4,2AMESH,ALL!楼板网格划分ASEL,S,,,25,96AATT,1,5,2AMESH,ALLNSEL,S,LOC,Z,0 !选取模型底端节点D,ALL,ALL !施加位移约束ALLSEL,ALL !重新选取所有节点/eshape,1.0/VIEW,1,1,1,1/ANG,1,270,XM,0/REPlotFINISH!************************************************************** ************ !************************第二步******************************************** !********************框架结构施加重力荷载*********************************ANTYPE,STATICNSEL,S,LOC,Z,0 !选取模型底端节点D,ALL,ALL !施加位移约束ALLSEL,ALL !重新选取所有节点ACEL,0,0,9.8SOLVEFINISH/POST1SET,FIRSTPLNSOL,U,Z,0,1FINISH!************************************************************** ************ !************************第三步******************************************** !*************框架结构施加楼面活荷载D=3KN/M^2***************************** !框架结构施加楼面活荷载D=3KN/M^2FINISH/SOLUANTYPE,STATICNSEL,S,LOC,Z,0 !选取模型底端节点D,ALL,ALL !施加位移约束ALLSEL,ALL !重新选取所有节点ACEL,0,0,9.8SOLVEFINISH/POST1SET,FIRSTPLNSOL,U,Z,0,1FINISH!************************************************************** ************ !************************第四步******************************************** !***********框架结构施加风荷载(基本风压=0.25KN/M^2,体形系数为1.0********** ANTYPE,STATIC*DIM,LOAD_1,ARRAY,12LOAD_1(1)=3.78,2.16,2.39,2.57,2.72,2.84,2.95,3.17,3.20,3.29, 3.39,3.51*DIM,LOAD_2,ARRAY,12LOAD_2(1)=7.56,4.32,4.78,5.14,5.44,5.68,5.90,6.34,6.40,6.58, 6.76,7.02*DIM,LOAD_3,ARRAY,12LOAD_3(1)=6.93,3.96,4.38,4.71,4.99,5.21,5.41,5.81,5.87,6.03, 6.20,6.44*DIM,LOAD_4,ARRAY,12LOAD_4(1)=3.15,1.80,1.99,2.14,2.27,2.37,2.46,2.64,2.67,2.74, 2.82,2.93*DIM,LOAD_A,ARRAY,12LOAD_A(1)=3.15,1.80,1.99,2.14,2.27,2.37,2.46,2.64,2.67,2.74, 2.82,2.93*DIM,LOAD_B,ARRAY,12LOAD_B(1)=6.30,3.60,3.98,4.28,4.54,4.74,4.92,5.28,5.34,5.48, 5.64,5.86*DIM,LOAD_C,ARRAY,12LOAD_C(1)=3.15,1.80,1.99,2.14,2.27,2.37,2.46,2.64,2.67,2.74, 2.82,2.93*DO,II,1,12FK,1+20*II,FY,LOAD_1(II)*ENDDO*DO,II,1,12FK,2+20*II,FY,LOAD_2(II)*ENDDO*DO,II,1,12FK,3+20*II,FY,LOAD_3(II)*ENDDO*DO,II,1,12FK,4+20*II,FY,LOAD_4(II)*ENDDOSOLVEFKDELE,ALL,ALL !第二次施加风荷载,删除第一次的。
ANSYS中的APDL命令
ANSYS中的APDL命令ANSYS中的APDL命令(一)(1).Lsel, type, item, comp, vmin, vmax, vinc, kswp 选择线 type: s 从全部线中选一组线r 从当前选中线中选一组线a 再选一部线附加给当前选中组aunoneu(unselect)inve: 反向选择item: line 线号loc 坐标length 线长comp: x,y,zkswp: 0 只选线1 选择线及相关关键点、节点和单元(2).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”仅用绝对值(3).Esel, type, item, comp, vmin, vmax, vinc, kabs 选择一组单元type: S: 选择一组单元(缺省)R: 在当前组中再选一部分作为一组A: 为当前组附加单元U: 在当前组中不选一部分单元All: 选所有单元None: 全不选Inve: 反向选择当前组Stat: 显示当前选择状态Item: Elem: 单元号Type: 单元类型号Mat: 材料号Real: 实常数号Esys: 单元坐标系号(4). 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次项的系数(5). 定义DP材料:首先要定义EX和泊松比:MP,EX,MAT,……MP,NUXY,MAT,……定义DP材料单元表(这里不考虑温度):TB,DP,MAT进入单元表并编辑添加单元表:TBDATA,1,CTBDATA,2,ψTBDATA,3,…… 如定义:EX=1E8,NUXY=0.3,C=27,ψ=45的命令如下:MP,EX,1,1E8MP,NUXY,1,0.3TB,DP,1TBDATA,1,27TBDATA,2,45这里要注意的是,在前处理的最初,要将角度单位转化到“度”,即命令:*afun,deg(6). 根据需要耦合某些节点自由度cp, nset, lab,,node1,node2,……node17nset: 耦合组编号lab: ux,uy,uz,rotx,roty,rotz ,allnode1-node17: 待耦合的节点号。
ANSYS-APDL的实例学习
ANSYS-APDL的实例学习10.1有限元模型的建立/TITLE,3D analysis on shield tunnel in Metro ! 确定分析标题/NOPR !菜单过滤设置/PMETH,OFF,0KEYW,PR_SET,1KEYW,PR_STRUC,1 !保留结构分析部分菜单/COM,/COM,Preferences for GUI filtering have been set to display: /COM,Structural1.材料、实常数和单元类型定义/clear !更新数据库/prep7 !进入前处理器et,1,solid45 !设置单元类型et,2,mesh200,6save !保持数据(2)定义模型中的材料参数。
!土体材料参数mp,ex,1,3.94e6 !地表层土弹性模量mp,prxy,1,0.35 !地表层土泊松比mp,dens,1,1828 !地表层土密度mp,ex,2,20.6e6 !盾构隧道所在地层参数mp,prxy,2,0.30mp,dens,2,2160mp,ex,3,500e6 !基岩地层参数mp,prxy,3,0.33mp,dens,3,2160!管片材料参数,管片衬砌按各向同性计算mp,ex,4,27.6e9 !管片衬砌弹性模量mp,prxy,4,0.2 !管片衬砌泊松比mp,dens,4,2500 !管片衬砌密度!注浆层,参数按水泥土取值mp,ex,5,1e9 !注浆层弹性模量mp,prxy,5,0.2 !注浆层泊松比mp,dens,5,2100 !注浆层密度save !保持数据第1章大型有限元软件ANSYS简介2 2.建立平面内模型并划分单元(1)在隧道中心线定义局部坐标,便于后来的实体选取。
local,11,0,0,0,0 !局部笛卡儿坐标local,12,1,0,0,0 !局部极坐标csys,11 !将当前坐标转换为局部坐标wpcsys,-1 !同时将工作平面转换到局坐标cyl4,,,,,2.7,90 !画部分圆半径为2.7cyl4,0,0,2.7,0,3,90 !画管片层部分圆cyl4,0,0,3,0,3.2,90 !画注浆层部分圆rectng,0,4.5,0,4.5 !画外边界矩形aovlap,all !做面递加nummrg,all !合并所有元素numcmp,all !压缩所有元素编号rectng,4.5,31.5,0,4.5 !画矩形面nummrg,all !合并所有元素numcmp,all !压缩所有元素编号save !保持数据(2)划分单元,如图10-1所示。
ANSYS建模apdl命令流实例应用
E,49*25-1+i+c,49*25+i+c *enddo *do,d,1,2*25,1 E,72*25-1+i+d,72*25+i+d *enddo !4 截面 type,1 mat,1 secnum,11 *do,a,1,3*25,1 E,9*25-1+i+a,9*25+i+a *enddo *do,b,1,4*25,1 E,26*25-1+i+b,26*25+i+b *enddo *do,c,1,3*25,1 E,51*25-1+i+c,51*25+i+c *enddo *do,d,1,4*25,1 E,68*25-1+i+d,68*25+i+d *enddo !3 截面 type,1 mat,1 secnum,10 *do,a,1,1*25,1 E,7*25-1+i+a,7*25+i+a *enddo *do,b,1,5*25,1 E,12*25-1+i+b,12*25+i+b *enddo *do,c,1,4*25,1 E,22*25-1+i+c,22*25+i+c *enddo *do,d,1,4*25,1 E,54*25-1+i+d,54*25+i+d *enddo *do,e,1,5*25,1 E,63*25-1+i+e,63*25+i+e *enddo !2 截面 type,1 mat,1 secnum,9 *do,a,1,1*25,1
/nerr,0 ! 将参数读入数据库 finish /clear /PREP7 ! 单元类型编号 et,1,beam44 et,2,shell63 et,3,link8 ! 材料类型编号 Q420/Q370/ 拉索 -1 ,混凝土 C40/C50-2 ,刚臂 -3 mp,ex,1,2.1e11 mp,prxy,1,0.3 mp,dens,1,7850 mp,alpx,1,1e-5
ANSYS中的APDL命令总结
在ANSYS中,命令流是由一条条ANSYS的命令组成的一个命令组合,这些命令按照一定顺序排布,能够完成一定的ANSYS功能,这些功能一般来说通过菜单操作也能够实现(而那些命令流能够实现,菜单操作实现不了的单个命令比较少见)。
以下命令是结合我自身经验,和前辈们的一些经验而总结出来的,希望对大家有帮助。
(1).Lsel, type, item, comp, vmin, vmax, vinc, kswp 选择线type: s 从全部线中选一组线r 从当前选中线中选一组线a 再选一部线附加给当前选中组aunoneu(unselect)inve: 反向选择item: line 线号loc 坐标length 线长comp: x,y,zkswp: 0 只选线1 选择线及相关关键点、节点和单元(2).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”仅用绝对值(3).Esel, type, item, comp, vmin, vmax, vinc, kabs 选择一组单元type: S: 选择一组单元(缺省)R: 在当前组中再选一部分作为一组A: 为当前组附加单元U: 在当前组中不选一部分单元All: 选所有单元None: 全不选Inve: 反向选择当前组Stat: 显示当前选择状态Item:Elem: 单元号Type: 单元类型号Mat: 材料号Real: 实常数号Esys: 单元坐标系号(4). 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次项的系数(5). 定义DP材料:首先要定义EX和泊松比:MP,EX,MAT,……MP,NUXY,MAT,……定义DP材料单元表(这里不考虑温度):TB,DP,MAT进入单元表并编辑添加单元表:TBDATA,1,CTBDATA,2,ψTBDATA,3,……如定义:EX=1E8,NUXY=0.3,C=27,ψ=45的命令如下:MP,EX,1,1E8MP,NUXY,1,0.3TB,DP,1TBDATA,1,27TBDATA,2,45这里要注意的是,在前处理的最初,要将角度单位转化到“度”,即命令:*afun,deg(6). 根据需要耦合某些节点自由度cp, nset, lab,,node1,node2,……node17nset: 耦合组编号lab: ux,uy,uz,rotx,roty,rotz ,allnode1-node17: 待耦合的节点号。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
/nerr,0 ! 将参数读入数据库 finish /clear /PREP7 ! 单元类型编号 et,1,beam44 et,2,shell63 et,3,link8 ! 材料类型编号 Q420/Q370/ 拉索 -1 ,混凝土 C40/C50-2 ,刚臂 -3 mp,ex,1,2.1e11 mp,prxy,1,0.3 mp,dens,1,7850 mp,alpx,1,1e-5
项目
混
主塔
凝 塔座
土
预应力钢筋
普通钢筋
索导管
钢料
垫板 检查设施
劲性骨架
锚具
波纹管 金属扩张网
材料规格 C50
C40
钢绞线 粗钢筋 Hrb335 20 号热轧无缝 Q235B Q235B Q235B 群锚 15-19 JL-32 粗钢筋锚具 100mm 45mm GM50075 型
单位 m3
m3
T T T T T T T 套 套 M M M2
! 中桁 *if,i,eq,94*25+2,then !7 截面 type,1 mat,1 secnum,14 *do,a,1,3*25,1 E,i-1+a,i+a *enddo *do,b,1,11*25,1 E,35*25-1+i+b,35*25+i+b *enddo *do,c,1,8*25,1 E,77*25-1+i+c,77*25+i+c *enddo !6 截面 type,1 mat,1 secnum,13 *do,a,1,3*25,1 E,i+3*25-1+a,i+3*25+a *enddo *do,b,1,3*25,1 E,32*25-1+i+b,32*25+i+b *enddo *do,c,1,3*25,1 E,46*25-1+i+c,46*25+i+c *enddo *do,d,1,1*25,1 E,i+76*25-1+d,i+76*25+d *enddo *do,e,1,7*25,1 E,87*25-1+i+e,87*25+i+e *enddo !5 截面 type,1 mat,1 secnum,12 *do,a,1,25,1 E,8*25-1+i+a,8*25+i+a *enddo *do,b,1,2*25,1 E,30*25-1+i+b,30*25+i+b *enddo *do,c,1,2*25,1
据不同的受力区段选用不同的杆件截面,在辅助墩附近的压重区梁段,腹杆采用箱型 截面杆件。主桁采用焊接杆件,整体节点。在节点外以高强度螺栓拼接的结构形式,
上下弦杆四面等强对接拼装。 H 型腹杆采用插入式连接。箱型腹杆采用四面与主桁节
点对拼的连接形式。主桁拼接采用 M30 高强螺栓。
表 1.2 安庆长江大桥主要工程数量表
name(19)='Sec_12' name(20)='Sec_13' name(21)='Sec_14' name(22)='Sec_15' name(23)='Sec_16' name(24)='Sec_x2'
*do,i,1,24,1 sectype,i,beam,mesh, secread,name(i),,,mesh *enddo
mp,ex,2,3.5e7 mp,prxy,2,0.16 mp,dens,2,2500 mp,alpx,2,1e-5
mp,ex,3,1e15 mp,prxy,3,1 mp,dens,3,0 mp,alpx,3,0
! 读取自定义截面并编号 *dim,name,char,24 name(1)='Sec_1a' name(2)='Sec_2a' name(3)='Sec_3a' name(4)='Sec_4a' name(5)='Sec_5a' name(6)='Sec_6a' name(7)='Sec_7a' name(8)='Sec_1' name(9)='Sec_2' name(10)='Sec_3' name(11)='Sec_4' name(12)='Sec_5' name(13)='Sec_6' name(14)='Sec_7' name(15)='Sec_8' name(16)='Sec_9' name(17)='Sec_10' name(18)='Sec_11'
结构构件 上弦杆、腹杆
下弦杆 主塔 塔座
斜拉索
表 2.3 模型实常数
结构部分 梁单元( beam44) 板壳单元( shell63 )
拉索单元( link8 )
实常数 R1 R2
R3,R4, R5, R6
截面性质 未定义参数 (由截面文件直接读入) 厚度 16mm (拉索面积) R,3,134.31e-4, R,4,115.84e-4, R,5,108.91e-4, R,6,100.000e-4,
E,49*25-1+i+c,49*25+i+c *enddo *do,d,1,2*25,1 E,72*25-1+i+d,72*25+i+d *enddo !4 截面 type,1 mat,1 secnum,11 *do,a,1,3*25,1 E,9*25-1+i+a,9*25+i+a *enddo *do,b,1,4*25,1 E,26*25-1+i+b,26*25+i+b *enddo *do,c,1,3*25,1 E,51*25-1+i+c,51*25+i+c *enddo *do,d,1,4*25,1 E,68*25-1+i+d,68*25+i+d *enddo !3 截面 type,1 mat,1 secnum,10 *do,a,1,1*25,1 E,7*25-1+i+a,7*25+i+a *enddo *do,b,1,5*25,1 E,12*25-1+i+b,12*25+i+b *enddo *do,c,1,4*25,1 E,22*25-1+i+c,22*25+i+c *enddo *do,d,1,4*25,1 E,54*25-1+i+d,54*25+i+d *enddo *do,e,1,5*25,1 E,63*25-1+i+e,63*25+i+e *enddo !2 截面 type,1 mat,1 secnum,9 *do,a,1,1*25,1
连续箱梁; 西引桥采用 15 孔梁长 32.6m 预应力混凝土简支梁及 2 孔梁长 24.6m 预应力
混凝土简支梁,其中宁安线采用箱梁,阜景线采用
T 梁。
主桥采用 103+188.5+580+217.5+159.5+117.5m 两塔钢桁斜拉桥方案,全长 1366m 。 主梁为三片主桁钢桁梁,桁间距 2x14m ,节间长 14.5m ,桁高 15m。主塔为钢筋混凝
E,6*25-1+i+a,6*25+i+a *enddo *do,b,1,2*25,1 E,17*25-1+i+b,17*25+i+b *enddo *do,c,1,1*25,1 E,21*25-1+i+c,21*25+i+c *enddo *do,d,1,1*25,1 E,58*25-1+i+d,58*25+i+d *enddo *do,e,1,2*25,1 E,61*25-1+i+e,61*25+i+e *enddo *do,f,1,2*25,1 E,74*25-1+i+f,74*25+i+f *enddo *do,g,1,1*25,1 E,85*25-1+i+g,85*25+i+g *enddo !1 截面 type,1 mat,1 secnum,8 *do,a,1,2*25,1 E,19*25-1+i+a,19*25+i+a *enddo *do,b,1,2*25,1 E,59*25-1+i+b,59*25+i+b *enddo *do,c,1,1*25,1 E,86*25-1+i+c,86*25+i+c *enddo
大桥全长 2996.8m ,其中主桥采用跨度为 101.5+188.5+580+217.5+159.5+116m 的
钢桁梁斜拉桥;非通航孔正桥采用 6 孔跨径 64m 预应力混凝土简支箱梁;东引桥采用
16 孔梁长 32.6m 预应力混凝土简支箱梁;跨大堤桥采用
48.9+86+48.8m 预应力混凝土
216 207
材料 Q370qD钢 Q420qD钢 C50混凝土 C40混凝土
斜拉索
表 2.2 安庆长江铁路大桥模型主要参数列表
弹性模量 (Pa) 2.1e11 2.1e11 3.5e7 3.5e7 2.1e11
泊松比 0.3 0.3 0.16 0.16 0.3