本文为 ANSYS 有限元 分析 系列博客中所采用的法兰分析实例的系列命令流。
FINISH ! 退出当前处理器
/CLEAR,ALL ! 清除所有
/PREP7 ! 进人前处理器
*AFUN,RAD ! 指定角度单位为弧度(非必要步骤)
pi = ACOS(-1) ! 获取圆周率pi的数值(非必要步骤)
*AFUN,DEG ! 将角度单位切换为度(非必要步骤)
/com,************* 1. 创建坐标系 *************
CSYS,0 ! 激活总体直角坐标系统
xc = 500 $ yc = 0 $ zc = 0 ! 局部坐标系原点位置(相对应当前激活坐标系)
thxy = 0 $ thyz = 0 $ thzx = 90 ! 局部坐标轴旋转角度度(相对应当前激活坐标系)
LOCAL,15,0,xc,yc,zc,thxy,thyz,thzx ! 创建局部直角坐标系15
CSYS,15 ! 激活局部坐标系15
xc = 1000 $ yc = 0 $ zc = 0 ! 局部坐标系原点位置(相对应当前激活坐标系)
thxy = 0 $ thyz = 0 $ thzx = 0 ! 局部坐标轴旋转角度度(相对应当前激活坐标系)
CLOCAL,16,1,xc,yc,zc,thxy,thyz,thzx ! 根据激活的局部坐标系15定义新的局部柱坐标系16
/com,************** 2. 创建几何 **************
CSYS,16 ! 激活局部柱坐标系16
WPCSYS,-1 ! 根据当前坐标系定义工作平面
d0 = 40 ! 螺栓孔直径为40
K,1,d0/2,0 $ K,2,d0/2,45 $ K,3,d0/2,90 $ K,4,d0/2,135 $ K,5,d0/2,180
L,1,2 ! 由两个关键点创建一条线(可以是直线或是弧线取决于坐标系)
*REPEAT,4,1,1 ! 上一条命令重复执行4次,始末关键点编号增量为1。
CSYS,15 ! 激活局部坐标系15
CLOCAL,17,1 ! 根据激活的局部坐标系15定义新的局部柱坐标系17
WPCSYS,-1, ! 根据当前坐标系定义工作平面
width = 88 ! 环带宽度为88
num = 55 ! 螺栓总数
Ds = 2000 ! 螺栓分度圆直径
K,6,Ds/2+width/2 $ K,7,Ds/2+width/2,180/num $ K,8,Ds/2,180/num
K,9,Ds/2-width/2,180/num $ K,10,Ds/2-width/2,0
L,6,7
*REPEAT,4,1,1
L,1,6
*REPEAT,5,1,1
/PNUM,KP,1 ! 显示关键点编号
/PNUM,LINE,1 ! 显示线编号
/REPLOT ! Replot
GPLOT ! Multi-Plots
AL,1,10,5,9 ! 创建平面,平面由边界线围成,边界线的编号分别为1、10、5、9。
*REPEAT,4,1,1,1,1 ! 上一条命令重复执行4次,各边界线编号增量均为1。
/PNUM,KP,0 ! 关闭关键点编号显示
/PNUM,LINE,0 ! 关闭线编号显示
/PNUM,AREA,1 ! 显示面的编号
/REPLOT ! Replot
GPLOT ! Multi-Plots
/com,************** 3. 定义材性 **************
ALLSEL,ALL ! *GET前最好全选,以免出现意向不到的错误。
*GET,mtmax,MAT,,NUM,MAX ! 获取当前定义的最大单元类型编号并将其赋予变量mtmax。
MPTEMP,,,,,,,, ! 为材料属性定义温度表
MPTEMP,1,0 ! 为材料属性定义温度表
MPDATA,EX,mtmax+3,,2e5 ! 定义mtmax+1号材料属性弹性模量
MPDATA,PRXY,mtmax+3,,0.28 ! 定义mtmax+1号材料属性泊松比
MPDATA,DENS,mtmax+3,,2700 ! 定义mtmax+1号材料属性密度
/com,************** 4. 定义单元 **************
! MESH200单元
*GET,etmax,ETYP,,NUM,MAX ! 获取当前定义的最大单元类型编号并将其赋予变量etmax (etmax=0)。
ET,etmax+2,MESH200 ! 定义单元类型为Mesh200,其单元识别号为etmax+2。
KEYOPT,etmax+2,1,7 ! 设置单元etmax+2的关键选项,即令 Keyoption(1)=7。
! SOLID186单元
ET,etmax+3,SOLID186 ! 定义单元类型为SOLID186,其单元识别号为etmax+3。
KEYOPT,etmax+3,2,1 ! 设置单元关键选项,将单元类型参考号为etmax+1的单元设置为完全积分单元。
/com,************** 5. 设置种子 **************
LESIZE,9, , ,2, , , , ,0 ! 设置组成面的边(线)的网格大小,编号为9的线被分成2段。
*REPEAT,5,1 ! 上一条命令重复执行5次
LESIZE,1, , ,2, , , , ,0 ! 设置线的网格大小,编号为1的线被分成2段。
*REPEAT,4,1 ! 上一条命令重复执行4次
LESIZE,8, , ,2, , , , ,0 ! 设置线的网格大小,编号为8的线被分成2段。
/com,************** 6. 创建网格 **************
MAT,mtmax+3 ! 给随后生成的单元激活一个材料号
TYPE,etmax+2 ! 给随后生成的单元激活一个单元类型号
MSHAPE,0,2D ! 指定划分单元的形状,生成四边形单元。
MSHKEY,1 ! 采用映射网格划分方式划分网格
AMESH,ALL ! 将所有面划分网格
/PNUM,ELEM,1 ! 显示单元编号
/REPLOT ! Replot
EPLOT ! Elements plots
/VIEW,1,-1 ! Left View
/ANG,1
/REP,FAST
MAT,mtmax+3 ! 给随后生成的单元激活一个材料号
TYPE,etmax+3 ! 给随后生成的单元激活一个单元类型号
! 执行拉伸
EXTOPT,ESIZE,8,1, ! EXTOPT:由面单元生成体单元的控制选项;
! 8表示在体生成或体扫略方向上单元分割数量为8;
! 1表示在体生成或体扫略方向上的间隔率为1(默认)。
EXTOPT,ACLEAR,1 ! 体单元网格生成后清除面单元网格即MESH200。
tf = 40 ! 法兰厚度为40
VEXT,ALL, , ,0,0,tf ! 通过给的偏移量由面生体
EPLOT ! Elements plots
/REPLOT ! Replot
/com,************** 7. 镜像及阵列 **************
CSYS,15 ! 激活局部直角坐标系15
VSYMM,Y,ALL, , , ,0,0 ! 对所有几何体进行镜像操作,镜像平面为xoz平面。
CSYS,17 ! 激活局部柱坐标系17
num = 55 ! 螺栓总数
VGEN,4,ALL, , , ,360/num, , ,0 ! 复制体,重复生成次数4,复制全部体。
! 当前激活坐标系下,关键点坐标值的偏移量为(--,360/num,0)
! 单个偏转角度为360/num
/PNUM,ELEM,0 ! 关闭单元编号显示
/REPLOT ! Replot
EPLOT ! Elements plots
NUMMRG,ALL ! 对实体进行合并
NUMCMP,ALL ! 压缩所定义项的编号
ALLSEL,ALL,ELEM ! Selects all elements.
CM,TopFlangeElems,ELEM ! 由所选实体生成一个组件
! 元件名称:TopFlangeElems;元件中的数据类型:单元。
CSYS,16 ! 激活局部柱坐标系16
CLOCAL,19,0 ! 创建局部直角坐标系19
VSYMM,Z,ALL, , , ,0,0 ! 对所有几何体进行镜像操作,镜像平面为局部坐标系19的xoY平面。
ALLSEL,ALL,ELEM ! Selects all elements.
CMSEL,U,TopFlangeElems ! 上一选择集中除去组件TopFlangeElems内的单元
CM,BotFlangeElems,ELEM ! 由所选实体生成一个组件
! 元件名称:BotFlangeElems;元件中的数据类型:单元。
/com,************ 8. 输出cdb文件 *************
ALLSEL,ALL ! 全选 select all entities
CDWRITE,DB,'PartFlange','cdb',,'','' ! 当前模型另存为 PartFlange.cdb 文件。

法兰模型
FINISH ! 退出当前处理器
/CLEAR,ALL ! 清除所有
/PREP7 ! 进人前处理器
*AFUN,RAD ! 指定角度单位为弧度(非必要步骤)
pi = ACOS(-1) ! 获取圆周率pi的数值(非必要步骤)
*AFUN,DEG ! 将角度单位切换为度(非必要步骤)
/com,************* 1. 创建坐标系 *************
CSYS,0 ! 激活总体直角坐标系统
xc = 500 $ yc = 0 $ zc = 0 ! 局部坐标系原点位置(相对应当前激活坐标系)
thxy = 0 $ thyz = 0 $ thzx = 90 ! 局部坐标轴旋转角度度(相对应当前激活坐标系)
LOCAL,15,0,xc,yc,zc,thxy,thyz,thzx ! 创建局部直角坐标系15
CSYS,15 ! 激活局部坐标系15
xc = 1000 $ yc = 0 $ zc = 0 ! 局部坐标系原点位置(相对应当前激活坐标系
免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删