ANSYS拱坝建模命令流

ANSYS拱坝建模命令流

请问哪位朋友有ANSYS拱坝建模命令流,谢谢!!!!!!!!!!!

==============================

某砌石拱坝位于U型河谷中，坝高55.5m，为单曲等厚拱坝，顶宽5m，底宽16m，坝顶弧长115.65m，弧高比2.1。首先输入基本参数，自底向上生成坝体。复制内容到剪贴板

代码:!坐标系原点位于拱坝轴线与拱坝参考面的交点在坝基的投影，X轴指向右岸，Y轴指向下游，Z轴铅直向上

!单位:长度-M,力-KN,加速度-M/S,密度-KN/M3

!

!

FINI

/CLEAR,START

/FILNAME,ARCHDAM,1

/PLOPTS,DATE,0

/TRIAD,LBOT

/VIEW,1,1,1,1

/VUP,1,Z

*AFUN,DEG

Z_UP=53!上游水位

Z_DOWN=0!下游水位

Z_SAND=38!泥沙淤积水位

DENS_SAND=0.6!淤积泥沙浮容重

ANG_FRI=16!淤积泥沙内摩擦角

LAYER_NUM=8!定义拱层数

*DIM,ELEVATION,ARRAY,LAYER_NUM!定义拱圈高程

*DIM,T_ARCH,ARRAY,LAYER_NUM!定义拱厚

*DIM,ARCH_RAD,ARRAY,LAYER_NUM,2!定义圆弧半径，上下游

*DIM,RAD_CEN,ARRAY,LAYER_NUM!定义圆弧中心Y值

*DIM,ARCH_ANGLE,ARRAY,LAYER_NUM,2!定义圆弧左右中心角

ELEVATION(1)=55.5,47,39,29,24,19,9,0!拱圈高程

H_DAM=ELEVATION(1)-ELEVATION(LAYER_NUM)!坝高

T_ARCH(1)=5.00,6.68,8.27,10.25,11.24,12.23,14.22,16.00!拱厚

ARCH_RAD(1,1)=68.00,68.00,68.00,68.00,68.00,68.00,68.00,68.00!外半径

*VOPER,ARCH_RAD(1,2),ARCH_RAD(1,1),SUB,T_ARCH(1)!计算内半径

RAD_CEN(1)=68.00,68.00,68.00,68.00,68.00,68.00,68.00,68.00!圆心位置,Y坐标

ARCH_ANGLE(1,1)=-47.72,-46.36,-45.06,-43.36,-41.88,-40.38,-37.32,-34.48!左中心角ARCH_ANGLE(1,2)=49.73,47.13,44.68,41.57,40.44,39.29,36.92,35.72!右中心角

W_DAM1=-ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,1))!坝顶左弦长

W_DAM2=ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,2))!坝顶右弦长

LOCAL,11,1,0,RAD_CEN(1),ELEVATION(1),-90!局部柱面坐标,原点位于顶拱圆心,X轴为拱坝中心线,指向上游,Y轴为中心角

!为查值计算坝体温度，重新定义一套数组

*DIM,TEMPUP,TABLE,LAYER_NUM!上游面温度

*DIM,TEMPDOWN,TABLE,LAYER_NUM!下游面温度

*DIM,T2_ARCH,TABLE,LAYER_NUM!拱厚

*DIM,Y_RAD,TABLE,LAYER_NUM!圆心到坝轴线距离

*DIM,RADUP,TABLE,LAYER_NUM!外半径

TEMPUP(1)=-8.84,-7.30,-5.66,-4.94,-4.99,-5.17,-5.73,-6.25

TEMPUP(1,0)=0,6,14,24,29,34,44,53

TEMPDOWN(1)=-8.84,-8.21,-7.36,-6.32,-5.86,-5.45,-4.72,-6.71

TEMPDOWN(1,0)=0,6,14,24,29,34,44,53

T2_ARCH(1)=5.00,6.68,8.27,10.25,11.24,12.23,14.22,16.00

T2_ARCH(1,0)=0,6,14,24,29,34,44,53

Y_RAD(1)=68.00,68.00,68.00,68.00,68.00,68.00,68.00,68.00

Y_RAD(1,0)=0,6,14,24,29,34,44,53

RADUP(1)=68.00,68.00,68.00,68.00,68.00,68.00,68.00,68.00

RADUP(1,0)=0,6,14,24,29,34,44,53

CSYS,0

!/PNUM,KP,1

!/PNUM,LINE,1

/PREP7

KNN=0

*DO,II,1,LAYER_NUM!生成拱坝控制关键点,每层6个点,上下游个三个点

PX1=ARCH_RAD(II,1)*SIN(ARCH_ANGLE(II,1))

PY1=RAD_CEN(II)-ARCH_RAD(II,1)*COS(ARCH_ANGLE(II,1))

PX2=ARCH_RAD(II,1)*SIN(ARCH_ANGLE(II,2))

PY2=RAD_CEN(II

)-ARCH_RAD(II,1)*COS(ARCH_ANGLE(II,2))

K,KNN+1,PX1,PY1,ELEVATION(II)

K,KNN+2,0,RAD_CEN(II)-ARCH_RAD(II,1),ELEVATION(II)

K,KNN+3,PX2,PY2,ELEVATION(II)

PX3=ARCH_RAD(II,2)*SIN(ARCH_ANGLE(II,1))

PY3=RAD_CEN(II)-ARCH_RAD(II,2)*COS(ARCH_ANGLE(II,1))

PX4=ARCH_RAD(II,2)*SIN(ARCH_ANGLE(II,2))

PY4=RAD_CEN(II)-ARCH_RAD(II,2)*COS(ARCH_ANGLE(II,2))

K,KNN+4,PX3,PY3,ELEVATION(II)

K,KNN+5,0,RAD_CEN(II)-ARCH_RAD(II,2),ELEVATION(II)

K,KNN+6,PX4,PY4,ELEVATION(II)

KNN=6*II

*ENDDO

LARC,1,3,2!生成拱圈线

*REPEAT,LAYER_NUM*2,3,3,3

L,1,4

*REPEAT,LAYER_NUM,6,6

L,3,6

*REPEAT,LAYER_NUM,6,6

ASKIN,1,3,5,7,9,11,13,15!生成上游面

ASKIN,2,4,6,8,10,12,14,16!生成下游面

ASKIN,17,18,19,20,21,22,23,24!生成侧面

ASKIN,25,26,27,28,29,30,31,32

AL,1,17,2,25

AL,15,24,16,32

VA,6,1,3,2,4,5

==============================

拉伸坝体侧面、底面，并进行切割，生成坝基复制内容到剪贴板

代码:!拉伸生成坝基

ASEL,S,AREA,,4!选择右侧面

VEXT,ALL,,,4*H_DAM!拉伸成基岩

ASEL,S,AREA,,3!选择左侧面

VEXT,ALL,,,-4*H_DAM!拉伸成基岩

ALLSEL

ASEL,S,LOC,Z,ELEVATION(LAYER_NUM)-1,ELEVATION(LAYER_NUM)+1!选择底面VEXT,ALL,,,,,-NINT(H_DAM/10)*10

ASEL,S,AREA,,10,20,5

ASEL,A,AREA,,23,27,4

VEXT,ALL,,,,4*H_DAM

ASEL,S,AREA,,8,18,5

ASEL,A,AREA,,25,29,4

VEXT,ALL,,,,-4*H_DAM

ALLSEL,ALL

SAVE

WPROTA,,90

WPOFFS,,,NINT(H_DAM*1.5/10)*10

VSBW,ALL,,DELETE

VSEL,S,LOC,Y,-200,-H_DAM*1.4

VDEL,ALL,,,1

VSEL,ALL

WPOFFS,,,-NINT(H_DAM*3.5/10)*10

VSBW,ALL,,DELETE

VSEL,S,LOC,Y,NINT(H_DAM*2/10)*10,NINT(H_DAM*2/10)*100

VDEL,ALL,,,1

VSEL,ALL

WPCSYS,,0

WPROTA,,,90

WPOFFS,,,NINT((W_DAM2+H_DAM)/10)*10

VSBW,ALL,,DELETE

VSEL,S,LOC,X,NINT((W_DAM2+H_DAM)/10)*10,NINT((W_DAM2+H_DAM)/10)*50 VDEL,ALL,,,1

VSEL,ALL

WPOFFS,,,-NINT((W_DAM2+W_DAM1+2*H_DAM)/10)*10

VSBW,ALL,,DELETE

VSEL,S,LOC,X,-NINT((W_DAM1+H_DAM)/10)*50,-NINT((W_DAM1+H_DAM)/10)*10 VDEL,ALL,,,1

WPCSYS,,0

VSEL,ALL

SAVE

==============================

用各控制高程切割整个模型，便于划分单元及计算成果和拱梁法比较复制内容到剪贴板代码:

*DO,II,LAYER_NUM,3,-1

WPOFFS,,,ELEVATION(II-1)-ELEVATION(II)

VSBW,ALL,,DELETE

*ENDDO

WPCSYS,,0

ALLSEL

WPROTA,,,90

VSBW,ALL,,DELETE

WPCSYS,,0

NUMMRG,ALL

NUMCMP,ALL

PX1=ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,1))

PX2=ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,2))

VSEL,S,LOC,X,PX1,PX2

VSEL,R,LOC,Z,ELEVATION(LAYER_NUM),ELEVATION(1)

CM,DAM,VOLU

VSEL,INVE

CM,BEDROCK,VOLU

ALLSEL

==============================

生成溢流堰复制内容到剪贴板

代码:

CMSEL,S,DAM

VSEL,R,LOC,Z,ELEVATION(1)-0.1,ELEVATION(2)+0.1

/VIEW,1,1,1,1

VPLOT

!生成溢流堰,为便于建模和施加荷载,与实际不太相符,根据圣维南原理,对大坝的总体应力不会有太大影响

W_WEIR=8!溢流堰单孔宽(坝轴线处)

NUM_WEIR=5!溢流堰孔数

W_WALL=2.5!溢流堰闸墩宽(坝

轴线处)

WPOFFS,,,ELEVATION(1)-ELEVATION(LAYER_NUM)

WPOFFS,,RAD_CEN(1)

WPROTA,,,90

*AFUN,RAD

PI=2*ACOS(0)

SITTA1=(W_WEIR*NUM_WEIR+W_WALL*(NUM_WEIR+1))/68*180/PI!溢流堰总中心角SITTA2=W_WEIR/ARCH_RAD(1,1)*180/PI!溢流堰单孔中心角

SITTA3=W_WALL/ARCH_RAD(1,1)*180/PI!溢流堰单个闸墩中心角

WPROTA,,-SITTA1/2

VSBW,ALL,,DELETE

*DO,II,1,NUM_WEIR

WPROTA,,SITTA3

VSBW,ALL,,DELETE

WPROTA,,SITTA2

VSBW,ALL,,DELETE

*ENDDO

WPROTA,,SITTA3

VSBW,ALL,,DELETE

SAVE

CMSEL,S,BEDROCK

VSEL,INVE

VSEL,R,LOC,Z,ELEVATION(1),ELEVATION(2)

CSYS,11

VSEL,R,LOC,Y,-SITTA1/2,SITTA1/2

CM,YLYY,VOLU

*DO,II,1,NUM_WEIR

CMSEL,S,YLYY

VSEL,R,LOC,Y,-SITTA1/2+II*SITTA3+(II-1)*SITTA2,-SITTA1/2+II*SITTA3+II*SITTA2

VDEL,ALL,,,1,,,1

*ENDDO

CSYS,0

CMSEL,S,BEDROCK

VSEL,INVE

CM,DAM,VOLU

CMSEL,S,DAM

VSEL,R,LOC,Z,ELEVATION(1),ELEVATION(2)

CSYS,11

VSEL,R,LOC,Y,-SITTA1/2,SITTA1/2,

CM,ZHADUN,VOLU

CSYS,0

WPCSYS,,0

VSEL,ALL

NUMMRG,ALL

NUMCMP,ALL

SAVE

==============================

划分单元复制内容到剪贴板

代码:!定义单元性质和材料

ET,1,SOLID95

MP,EX,1,1E7

MP,NUXY,1,0.25

MP,DENS,1,23!坝体容重

MP,ALPX,1,0.7E-5

MP,REFT,1,0

MP,EX,2,0.8E7

MP,NUXY,2,0.21

MP,ALPX,2,0.0

MP,REFT,2,0

!坝体单元划分

CMSEL,S,DAM

CMSEL,U,ZHADUN

ASLV,S

LSLA,S

*AFUN,DEG

PX1=ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,1))

PX2=ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,2))

PX3=ARCH_RAD(LAYER_NUM,2)*SIN(ARCH_ANGLE(LAYER_NUM,1)) PX4=ARCH_RAD(LAYER_NUM,2)*SIN(ARCH_ANGLE(LAYER_NUM,2)) LSEL,U,LOC,X,PX3,PX1

LSEL,U,LOC,X,PX4,PX2

LSEL,U,LOC,X,0

LSEL,U,LOC,Z,ELEVATION(1),ELEVATION(2) LESIZE,ALL,,,20,

CM,LTEMP1,LINE

CSYS,0

CMSEL,S,DAM

ASLV,S

LSLA,S

CMSEL,U,LTEMP1

CM,LTEMP2,LINE

LSEL,U,LENGTH,,T_ARCH(LAYER_NUM)+0.1,1000 LSEL,R,TAN1,Z

CSYS,11

LSEL,U,LOC,X,ARCH_RAD(1,1)

LSEL,U,LOC,X,ARCH_RAD(1,1)-T_ARCH(1) LSEL,U,LOC,X,ARCH_RAD(1,1)-T_ARCH(2) CSYS,0

CM,LTEMP3,LINE

LESIZE,ALL,,,5,

ALLSEL

CMSEL,S,LTEMP2

CMSEL,U,LTEMP3

LSEL,U,TAN1,Z

CM,LTEMP4,LINE

LESIZE,ALL,,,3,

CMSEL,S,DAM

CMSEL,U,ZHADUN

ASLV,S

LSLA,S

LSEL,R,LOC,Z,ELEVATION(1),ELEVATION(2) CSYS,11

LSEL,U,LOC,Y,-SITTA1/2-1,SITTA1/2+1 LSEL,U,LENGTH,,0,T_ARCH(1)+0.1

CSYS,0

CMSEL,U,LTEMP3

CMSEL,U,LTEMP4

LESIZE,ALL,,,10,

CM,LTEMP5,LINE

CMSEL,S,DAM

CMSEL,U,ZHADUN

ASLV,S

LSLA,S

LSEL,R,LOC,Z,ELEVATION(1),ELEVATION(2) CMSEL,U,LTEMP3

CMSEL,U,LTEMP4

CMSEL,U,LTEMP5

LSEL,U,LENGTH,,0,3

CM,LTEMP6,LINE

LESIZE,ALL,,,4,

CSYS,0

ALLSEL

MSHAPE,0

CMSEL,S,DAM

VSWEEP,ALL

SAVE

!坝基单元划分

CMSEL,S,BEDROCK

ASLV,S

LSLA,S

LSEL,U,TAN1,X

LSEL,U,LOC,X,PX2,PX1

CM,LTEMP7,LINE

LESIZE,ALL,,,5,,,,,1

CMSEL,S,BEDROCK

ASLV,S

LSLA,S

LSEL,U,TAN1,Y

PY4=RAD_CEN(LAYER_NUM)-ARCH_RAD(LAYER_NUM,2)*COS(ARCH_ANGLE(LAYER_NUM,2)) +20

LSEL,U,LOC,Y,-10,PY4

CM,LTEMP8,LINE

LESIZE,ALL,,,6,4,,,,1

CMSEL,S,BEDROCK

ASLV,S

LSLA,S

LSEL,U,TAN1,Z

LSEL,R,LOC,Z,ELEVATION(LAYER_NUM),ELEVATION(LAYER_NUM)-200

LESIZE,ALL,,,5,5,,,,1

CMSEL,S,BEDROCK

ASLV,S

LSLA,S

LSEL,U,TAN1,Z

LSEL,R,LOC,Z,ELEVATION(1),ELEVATION(LAYER_NUM)

LESIZE,ALL,,,3,

CMSEL,S,BEDROCK

ASLV,S

LSLA,S

LSEL,R,LOC,Z,-50,-1000

LSEL,R,LOC,Y,-10,PY4

CSYS,11

LSEL,U,LOC,X,ARCH_RAD(1,1)

LSEL,U,LOC,X,ARCH_RAD(1,1)-T_ARCH(LAYER_NUM)

CSYS,0

LSEL,R,LENGTH,,T_ARCH(LAYER_NUM)-5,T_ARCH(LAYER_NUM)+3,

LESIZE,ALL,,,5,

CSYS,0

CMSEL,S,BEDROCK

ASLV,S

LSLA,S

LSEL,R,LOC,Z,ELEVATION(LAYER_NUM)-50,ELEVATION(LAYER_NUM)-1000

LSEL,R,LOC,X,PX2,PX1

LESIZE,ALL,,,20,

CMSEL,S,BEDROCK

MAT,2

VMESH,ALL

!单元划分完毕

ALLSEL

NUMMRG,ALL

NUMCMP,ALL

SAVE

==============================

定义边界条件，顺河向约束Y向位移，横河向约束X向位移，底面约束竖向Z位移复制内容到剪贴板

代码:

FINI

/SOLU

!边界位移条件

*GET,NXMAX,NODE,,MXLOC,X

*GET,NXMIN,NODE,,MNLOC,X

*GET,NYMAX,NODE,,MXLOC,Y

*GET,NYMIN,NODE,,MNLOC,Y

*GET,NZMIN,NODE,,MNLOC,Z

NSEL,S,LOC,X,NXMAX

NSEL,A,LOC,X,NXMIN

D,ALL,UX,0!约束上下游边界X向位移

NSEL,S,LOC,Y,NYMAX

NSEL,A,LOC,Y,NYMIN

D,ALL,UY,0!约束左右边界Y向位移

NSEL,S,LOC,Z,NZMIN

D,ALL,ALL,0!约束底面边界Z向位移

NSEL,ALL

==============================

施加水荷载，淤沙压力，坝体自重复制内容到剪贴板代码:!大坝荷载

ALLSEL

CMSEL,S,DAM

CMSEL,U,ZHADUN

ASLV,S

CSYS,11

ASEL,R,LOC,X,ARCH_RAD(1,1)-0.1,ARCH_RAD(1,1)+0.1 NSLA,S,1

CM,N_DAMUP,NODE!定义大坝上游面节点CSYS,0

ALLSEL

CMSEL,S,DAM

ASLV,S

ASEL,R,EXT

LSLA,S

LSEL,R,LOC,X,0

LSEL,U,LOC,Y,-0.1,0.1

LSEL,U,LOC,Z,ELEVATION(2)

LSEL,U,LOC,Z,ELEVATION(LAYER_NUM)

ASLL,S

ASEL,U,LOC,X,0

LSLA,S

NSLA,S,1

CM,N_DAMDOWN,NODE!定义大坝下游面节点ALLSEL

CMSEL,S,YLYY

ASLV,S

CSYS,11

ASEL,R,LOC,X,ARCH_RAD(1,1)-0.1,ARCH_RAD(1,1)+0.1 ASEL,U,,,291,326,326-291

NSLA,S,1

CM,N_ZHADUN1,NODE!定义溢流坝中墩节点ASEL,S,,,291,326,326-291

NSLA,S,1

CM,N_ZHADUN2,NODE!定义溢流坝边墩节点CSYS,0

*IF,Z_UP,NE,ELEVATION(LAYER_NUM),THEN

CMSEL,S,N_DAMUP

NSEL,R,LOC,Z,0,Z_UP

/PSF,PRES,NORM,2,0.1

SFGRAD,PRES,0,Z,0,-9.81

SF,ALL,PRESS,Z_UP*9.81!施加大坝上游面水压力

*ENDIF

*IF,Z_UP,GT,ELEVATION(2),THEN

CMSEL,S,N_ZHADUN1

NSEL,R,LOC,Z,ELEVATION(2),Z_UP

SFGRAD,PRES,0,Z,0,-(W_WEIR+W_WALL)/W_WALL*9.81

SF,ALL,PRESS,Z_UP*(W_WEIR+W_WALL)/W_WALL*9.81!施加中墩水压力,包含闸门传递的水压力

CMSEL,S,N_ZHADUN2

NSEL,R,LOC,Z,ELEVATION(2),Z_UP

SFGRAD,PRES,0,Z,0,-(W_WEIR/2+W_WALL)/W_WALL*9.81

SF,ALL,PRESS,Z_UP*(W_WEIR/2+W_WALL)/W_WALL*9.81!施加边墩水压力,包含闸门传递的水压力

*ENDIF

*IF,Z_SAND,GT,ELEVATION(LAYER_NUM),THEN

SFCUM,PRES,ADD!荷载为累加方式

CMSEL,S,N_DAMUP

NSEL,R,LOC,Z,0,Z_SAND

*AFUN,DEG

SAND_GRADS=DENS_SAND*(TA

N(45-ANG_FRI/2))**2*9.81!计算淤沙压力梯度

SFGRAD,PRES,0,Z,0,-SAND_GRADS

SF,ALL,PRESS,Z_SAND*SAND_GRADS!施加淤沙压力

SFCUM,PRES,REPL!荷载为替代方式

*ENDIF

ALLSEL

*IF,Z_DOWN,NE,ELEVATION(LAYER_NUM),THEN

CMSEL,S,N_DAMDOWN

NSEL,R,LOC,Z,0,Z_DOWN

SFGRAD,PRES,0,Z,0,-9.81

SF,ALL,PRESS,Z_DOWN*9.81!施加大坝下游面水压力

*ENDIF

ALLSEL

ACEL,0,0,1!施加大坝自重荷载

SAVE

==============================

施加温度体荷载复制内容到剪贴板

代码:!施加大坝温度荷载

CMSEL,S,DAM

NSLV,S,1

*GET,NCOUNT_DAM,NODE,,COUNT

*GET,NMIN_DAM,NODE,,NUM,MIN

TUNIF,0,

TREF,0,

NTT=NMIN_DAM

*DO,II,1,NCOUNT_DAM

Z_NTT=NZ(NTT)!Z

Z_TEMP=ELEVATION(1)-Z_NTT!到坝基高度

TL_ARCH=T2_ARCH(Z_TEMP)!拱厚

CSYS,11

X_NTT=NX(NTT)!X:距顶拱圆心距离

CSYS,0

RAD0=Y_RAD(Z_TEMP)!圆心到坝轴线距离

RAD1=RADUP(Z_TEMP)!外半径

LXUP=RAD_CEN(1)-RAD0+RAD1-X_NTT

TTT1=TEMPUP(Z_TEMP)

TTT2=TEMPDOWN(Z_TEMP)

TTT=TTT1+(TTT2-TTT1)/TL_ARCH*LXUP

BF,NTT,TEMP,TTT

*IF,II,LT,NCOUNT_DAM,THEN

NTT=NDNEXT(NTT)

*ENDIF

*ENDDO

==============================

计算求解，采用PCG求解器

ALLSEL

SAVE

EQSLV,PCG,1E-8

SOLVE

SAVE

FINI

最大变形2cm，位于拱冠梁顶，与拱梁分载法一致。

ANSYS APDL命令流学习参数化建模

第一天 目标:熟悉ANSYS基本关键字的含义 k --> Keypoints 关键点 l --> Lines 线 a --> Area 面 v --> V olumes 体 e --> Elements 单元 n --> Nodes 节点 cm --> component 组元 et --> element type 单元类型 mp --> material property 材料属性 r --> real constant 实常数 d --> DOF constraint 约束 f --> Force Load 集中力 sf --> Surface Force on nodes 表面载荷 bf --> Body Force on Nodes 体载荷 ic --> Initial Conditions 初始条件 第二天 目标:了解命令流的整体结构,掌握每个模块的标识 !文件说明段 /BATCH /TITILE,test analysis !定义工作标题/FILENAME,test !定义工作文件名 /PREP7 !进入前处理模块标识!定义单元,材料属性,实常数段 ET,1,SHELL63 !指定单元类型 ET,2,SOLID45 !指定体单元 MP,EX,1,2E8 !指定弹性模量 MP,PRXY,1,0.3 !输入泊松比 MP,DENS,1,7.8E3 !输入材料密度 R,1,0.001 !指定壳单元实常数-厚度...... !建立模型 K,1,0,0,, !定义关键点 K,2,50,0,, K,3,50,10,, K,4,10,10,, K,5,10,50,, K,6,0,50,, A,1,2,3,4,5,6, !由关键点生成面...... !划分网格 ESIZE,1,0, AMESH,1 ...... FINISH !前处理结束标识

ANSYS命令流实例

/PREP7 !进入前处理 ANTYPE,STATIC !设置分析类型为静力结构分析 PSTRES,ON !用于后面的模态分析中考虑预应力(该开关不影响静力分析) ET,1,LINK10 !选取单元类型1(单向杆单元) KEYOPT,1,3,0 !设置仅承受拉应力，KEYOPT(3)=0 R,1,306796E-8,543248E-8 !设置实常数，包括绳索截面积(306796E-8)，初始应变(543248E-8) MP,EX,1,30E6 !定义材料的弹性模量(1号材料) MP,DENS,1,73E-5 !定义材料的密度(1号材料) N,1 ! 定义第1号节点 N,14,100 ! 定义第14号节点 FILL ! 均分填满第2号至第13号节点 E,1,2 !由节点1及节点2生成单元 EGEN,13,1,1 !依序复制生成13个单元 D,ALL,ALL ! 对所有节点施加固定约束 FINISH ! 前处理结束 /SOLU ! 进入求解模块，求解预应力引起的应力状态 SOLVE ! 求解 FINISH ! 退出求解模块 /POST1 ! 进入一般的后处理 ETABLE,STRS,LS,1 !针对LINK10单元,建立单元列表STRS,通过LS及特征号1来获得单元的轴向应力 *GET,STRSS,ELEM,13,ETAB,STRS !针对单元列表STRS, 提取13号单元的应力 FINISH ! 后处理结束 /POST26 ! 进入时间历程后处理，处理支反力 RFORCE,2,1,F,X !将1号节点上的x方向支反力提取，并存储到2号变量中 STORE ! 存储 *GET,FORCE,V ARI,2,EXTREM,VMAX !将2号变量的最大值赋给参数FORCE /SOLU ! 再次进入求解模块，模态分析 ANTYPE,MODAL ! 模态分析 MODOPT,SUBSP,3 ! 选择子空间迭代法，求3阶模态 MXPAND,3 ! 设定3阶模态扩展 PSTRES,ON ! 用于在模态分析中考虑预应力(还需在前面的静力分析中也同时打开) DDELE,2,UX,13 ! 删除从2号节点到13号节点上的UX约束 DDELE,2,UY,13 !删除从2号节点到13号节点上的UY约束 SOLVE !求解 *GET,FREQ1,MODE,1,FREQ ! 提取第1阶模态共振频率，并赋值给参数FREQ1 *GET,FREQ2,MODE,2,FREQ ! 提取第2阶模态共振频率，并赋值给参数FREQ2 *GET,FREQ3,MODE,3,FREQ ! 提取第3阶模态共振频率，并赋值给参数FREQ3 *STATUS !列出所有参数的实际内容

ansys旋转经典命令流

1 旋转摩擦 (1) 2. 电磁三d命令流实例（论坛看到） (11) 3. 帮助感应加热例子induction heating of a solid cylinder billet (15) 4. 感应加热温度场的数值模拟（论文）inducheat30命令流 (19) 5. 如何施加恒定的角速度？Simwe仿真论坛 (24) 6. 旋转一个已经生成好的物体 (27) 7. 产生这样的磁力线 (28) 8. 旋转摩擦生热简单例子（二维旋转） (32) 8.1. 原版 (32) 8.2. 部分gui操作 (35) 9. VM229 Input Listing (39) 10 轴承---耦合+接触分析 (47) 11. 板的冲压仿真 (52) 1 旋转摩擦 FINISH /FILNAME,Exercise24 !定义隐式热分析文件名 /PREP7 !进入前处理器 ET,1,SOLID5 !选择单元类型 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,7800 !定义材料1的密度 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,C,1,,460 !定义材料1的比热 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,KXX,1,,66.6 !定义材料1的热传导系数 MPTEMP,,,,,,,, MPTEMP,1,0 UIMP,1,REFT,,,30 !定义材料1的热膨胀系数的参考温度 MPDATA,ALPX,1,,1.06e-5 !定义材料1的热膨胀系数MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,206e9 !定义材料1的弹性模量 MPDATA,PRXY,1,,0.3 !定义材料1的泊松比 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,2,,8900 !定义材料2的密度 MPTEMP,,,,,,,, MPTEMP,1,0

Ansys常见命令流

Ansys命令流 第一天 目标:熟悉ANSYS基本关键字的含义 k --> Keypoints 关键点 l --> Lines 线 a --> Area 面 v --> Volumes 体 e --> Elements 单元 n --> Nodes 节点 cm --> component 组元 et --> element type 单元类型 mp --> material property 材料属性 r --> real constant 实常数 d --> DOF constraint 约束 f --> Force Load 集中力 sf --> Surface load on nodes 表面载荷 bf --> Body Force on Nodes 体载荷 ic --> Initial Conditions 初始条件 第二天 目标:了解命令流的整体结构,掌握每个模块的标识 !文件说明段 /BATCH /TITILE,test analysis !定义工作标题 /FILENAME,test !定义工作文件名 /PREP7 !进入前处理模块标识 !定义单元,材料属性,实常数段 ET,1,SHELL63 !指定单元类型 ET,2,SOLID45 !指定体单元 MP,EX,1,2E8 !指定弹性模量 MP,PRXY,1,0.3 !输入泊松比 MP,DENS,1,7.8E3 !输入材料密度 R,1,0.001 !指定壳单元实常数-厚度...... !建立模型 K,1,0,0,, !定义关键点 K,2,50,0,,

K,3,50,10,, K,4,10,10,, K,5,10,50,, K,6,0,50,, A,1,2,3,4,5,6, !由关键点生成面 ...... !划分网格 ESIZE,1,0, AMESH,1 ...... FINISH !前处理结束标识 /SOLU !进入求解模块标识 !施加约束和载荷 DL,5,,ALL SFL,3,PRES,1000 SFL,2,PRES,1000 ...... SOLVE !求解标识 FINISH !求解模块结束标识 /POST1 !进入通用后处理器标识 ...... /POST26 !进入时间历程后处理器 …… /EXIT,SAVE !退出并存盘 以下是日志文件中常出现的一些命令的标识说明,希望能给大家在整理LOG文件时有所帮助 /ANGLE !指定绕轴旋转视图 /DIST !说明对视图进行缩放 /DEVICE !设置图例的显示,如:风格,字体等 /REPLOT !重新显示当前图例 /RESET !恢复缺省的图形设置 /VIEW !设置观察方向 /ZOOM !对图形显示窗口的某一区域进行缩放

ANSYS命令流学习笔记10-利用APDL在WorkBench中进行非线性屈曲分析

!ANSYS命令流学习笔记10-利用APDL在WorkBench中进行非线性屈曲分析 !学习重点： !1、强化非线性屈曲知识 首先了解屈曲问题。在理想化情况下,当F < Fcr时, 结构处于稳定平衡状态，若引入一个小的侧向扰动力，然后卸载, 结构将返回到它的初始位置。当F > Fcr时, 结构处于不稳定平衡状态, 任何扰动力将引起坍塌。当F = Fcr时,结构处于中性平衡状态，把这个力定义为临界载荷。在实际结构中, 几何缺陷的存在或力的扰动将决定载荷路径的方向。在实际结构中, 很难达到临界载荷，因为扰动和非线性行为, 低于临界载荷时结构通常变得不稳定。 要理解非线性屈曲分析，首先要了解特征值屈曲。特征值屈曲分析预测一个理想线弹性结构的理论屈曲强度，缺陷和非线性行为阻止大多数实际结构达到理想的弹性屈曲强度，特征值屈曲一般产生非保守解, 使用时应谨慎。 !理论解，根据Euler公式。其中μ取决于固定方式。 !有限元方法， 已知在特征值屈曲问题： 求解，即可得到临界载荷 而非线性屈曲问题： 其中为结构初始刚度，为有缺陷的结构刚度，为位移矩阵，为载荷矩阵。 非线性屈曲分析时考虑结构平衡受扰动（初始缺陷、载荷扰动）的非线性静力分析，该分析时一直加载到结构极限承载状态的全过程分析，分析中可以综合考虑材料塑性、几何非线性、接触、大变形。非线性屈曲比特征值屈曲更精确，因此推荐用于设计或结构的评价。 !2、熟悉WB中非线性屈曲分析流程 (1) 前处理，施加单元载荷，进行预应力静力分析。 (2) 基于预应力静力分析，指定分析类型为特征值屈曲分析，完成特征值屈曲分析。 (3) 在APDL模块将一阶特征屈曲模态位移乘以适当系数，将此变形后的形状当做非线性分析的初始模型。

ansys实例命令流-弹塑性分析命令流

/FILNAME,Elastic-Plasitc,1 /TITLE, Elastic-Plasitc Analysis !前处理。 /PREP7 !**定义梁单元189。 ET,1,BEAM189 !定义单元。 !**梁截面1。 SECTYPE, 1, BEAM, HREC, , 0 !定义梁截面。SECOFFSET, CENT SECDATA,50,100,6,6,6,6,0,0,0,0 !定义梁截面完成。 !**定义材料。 MPTEMP,,,,,,,, !定义弹塑性材料模型。MPTEMP,1,0 MPDATA,EX,1,,2.05e5 MPDATA,PRXY,1,,0.3 TB,BISO,1,1,2, TBTEMP,0 TBDATA,,150,18600,,,, !定义弹塑性材料模型。!**建立几何模型。 K,1, , , , K,2 ,900, K,3 ,,50 LSTR, 1, 2 !**网格划分。 FLST,5,1,4,ORDE,1 !定义网格密度。FITEM,5,1 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,50, , , , ,1 !定义网格密度完成。CM,_Y,LINE !网格划分。 LSEL, , , , 1 CM,_Y1,LINE CMSEL,S,_Y CMSEL,S,_Y1 LATT,1, ,1, , 3, ,1 CMSEL,S,_Y CMDELE,_Y CMDELE,_Y1 LMESH, 1 !网格划分完成。 !施加载荷及求解。 FINISH /SOL

!**施加约束。 FLST,2,1,3,ORDE,1 !施加约束。FITEM,2,1 /GO DK,P51X, , , ,0,UX,UY,UZ,ROTX, , , FLST,2,1,3,ORDE,1 FITEM,2,2 /GO DK,P51X, , , ,0,UY,UZ,ROTX, , , , !施加约束完成。 !**加载。 FLST,2,50,2,ORDE,2 FITEM,2,1 FITEM,2,-50 SFBEAM,P51X,1,PRES,100, , , , , , LSWRITE,1, !定义载荷步1完成。FLST,2,50,2,ORDE,2 !定义载荷步2。FITEM,2,1 FITEM,2,-50 SFEDELE,P51X,1,PRES LSWRITE,2, !定义载荷步2完成。!设定求解步并求解。 LSSOLVE,1,2,1,

个人总结ansys命令流

Q235 属性：弹性模量E=2.1e5 N/mm2 密度=7.85e-6kg/mm3 泊松比=0.3 mp,ex,1,2.1e5 mp,prxy,1,0.3 mp,dens,1,7.85e-6 1,ksymm 镜像点 2,arsym 镜像面 3,kgen 复制点 4.adele删除面 6，kdist,k1,k2 测量两关键点的距离 7,adele,a,,,1 删除area and below 8,创建圆柱面： circle 创建圆 然后创建直线 然(轴线) 利用拉伸命令创建圆柱面creat__areas__by Lines adrag 线拉伸成面modeling>operate>extrude>lines>>along lines VDRAG 面拉伸成体modeling>operate>extrude>areas>>along lines ！创建空心圆柱体 这个命令 CYLIND, RAD1, RAD2, Z1, Z2, THETA1, THETA2 Main Menu>Preprocessor>Modeling>Create>Volumes>Cylinder>By Dimensions Main Menu>Preprocessor>Trefftz Domain>TZ Geometry>Create>Volume>Cylinder>By Dimensions 9,aptn 分割面 10,asbw 用工作平面切割面 11.wpoffs 12.wprota

https://www.wendangku.net/doc/c515931099.html,ng 过圆外一点做圆的切线(0°或180°) 14，nummrg 将重复的点消除 15，asba 面减去面 16,两个圆柱面的相贯线作法：做出两个相穿的圆柱面，利用APTN命令 17，选择面，不选择一部分面 asel,u,loc,z,kz(735) 18.在工作平面上生成一个矩形面 RECTING,X1,X2,Y1,Y2 X1，X2——矩形在工作平面X方向坐标值的变化范围 Y1,Y2——矩形在工作平面Y方向坐标值的变化范围 18,圆阵列 建立工作平面与圆柱的横截面平行，在工作平面情况下建立局部坐标系（柱坐标系），然后利用agen命令复制。 19，转换成局部柱坐标系 20,kfill 在两个关键点之间生成一个或多个关键点 21.网格划分 aatt,1,14,1, !aatt,mat,real,type,esys,secn aesize,all,1000 !aesize,anum,size, 单元尺寸 mshape,0,2d !mshape,key,dimension 指定划分单元形状amesh,all k,1,24000,33000,2230 k,2,24000,33000,-2230 k,3,-24000,33000,-2230 k,4,-24000,33000,2230 kfill,2,3,23,5,1,1 kfill,1,4,23,28,1,1 *do,i,5,26 l,i,i+1 *enddo

ANSYS-结构稳态(静力)分析之经典实例-命令流格式

ANSYS 结构稳态(静力)分析之经典实例-命令流格式.txt两人之间的感情就像织毛衣，建立 的时候一针一线，小心而漫长，拆除的时候只要轻轻一拉。。。。/FILNAME,Allen-wrench,1 ! Jobname to use for all subsequent files /TITLE,Static analysis of an Allen wrench /UNITS,SI ! Reminder that the SI system of units is used /SHOW ! Specify graphics driver for interactive run; for batch ! run plots are written to pm02.grph ! Define parameters for future use EXX=2.07E11 ! Young's modulus (2.07E11 Pa = 30E6 psi) W_HEX=.01 ! Width of hex across flats (.01m=.39in) *AFUN,DEG ! Units for angular parametric functions定义弧度单位 W_FLAT=W_HEX*TAN(30) ! Width of flat L_SHANK=.075 ! Length of shank (short end) (.075m=3.0in) L_HANDLE=.2 ! Length of handle (long end) (.2m=7.9 in) BENDRAD=.01 ! Bend radius of Allen wrench (.01m=.39 in) L_ELEM=.0075 ! Element length (.0075 m = .30 in) NO_D_HEX=2 ! Number of divisions on hex flat TOL=25E-6 ! Tolerance for selecting nodes (25e-6 m = .001 in) /PREP7 ET,1,SOLID45 ! 3维实体结构单元；Eight-node brick element ET,2,PLANE42 ! 2维平面结构；Four-node quadrilateral (for area mesh) MP,EX,1,EXX ! Young's modulus for material 1；杨氏模量 MP,PRXY,1,0.3 ! Poisson's ratio for material 1；泊松比 RPOLY,6,W_FLAT ! Hexagonal area创建规则的多边形 K,7 ! Keypoint at (0,0,0) K,8,,,-L_SHANK ! Keypoint at shank-handle intersection K,9,,L_HANDLE,-L_SHANK ! Keypoint at end of handle L,4,1 ! Line through middle of hex shape L,7,8 ! Line along middle of shank L,8,9 ! Line along handle LFILLT,8,9,BENDRAD ! Line along bend radius between shank and handle! 产生 一个倒角圆，并生成三个点 /VIEW,,1,1,1 ! Isometric view in window 1 /ANGLE,,90,XM ! Rotates model 90 degrees about X! 不用累积的旋转 /TRIAD,ltop /PNUM,LINE,1 ! Line numbers turned on LPLOT

ansys命令流

第一天目标： 熟悉ANSYS基本关键字的含义k --> Keypoints关键点l --> Lines线a --> Area 面v --> Volumes体e --> Elements单元n --> Nodes节点cm --> component组元et --> element type单元类型mp --> material property材料属性r --> real constant实常数d --> DOF constraint约束f --> Force Load集中力sf --> Surface load on nodes 表面载荷bf --> Body Force on Nodes体载荷ic --> Initial Conditions初始条件第二天目标： 了解命令流的整体结构,掌握每个模块的标识!文件说明段/BATCH/TILE,test analysis!定义工作标题/FILENAME,test!定义工作文件名/PREP7!进入前处理模块标识!定义单元,材料属性,实常数段ET,1,SHELL63!指定单元类型ET,2,SOLID45!指定体单元MP,EX,1,2E8!指定弹性模量MP,PRXY,1, 0.3!输入泊松比MP,DENS,1, 7.8E3!输入材料密度R,1, 0.001!指定壳单元实常数-厚度......!建立模型K,1,0,0,,!定义关键点 K,2,50,0,,K,3,50,10,,K,4,10,10,,K,5,10,50,,K,6,0,50,,A,1,2,3,4,5,6,!由关键点生成面......!划分网格ESIZE,1,0,AMESH, 1......FINISH!前处理结束标识/SOLU!进入求解模块标识!施加约束和载荷DL,5,,ALLSFL,3,PRES,1000SFL,2,PRES, 1000......SOLVE!求解标识FINISH!求解模块结束标识/POST1!进入通用后处理器标识....../POST26!进入时间历程后处理器……/EXIT,SAVE!退出并存盘以下是日志文件中常出现的一些命令的标识说明,希望能给大家在整理LOG文件时有所帮助/ANGLE!指定绕轴旋转视图/DIST!说明对视图进行缩放/DEVICE!设置图例的显示,如： 风格,字体等/REPLOT!重新显示当前图例/RESET!恢复缺省的图形设置/VIEW!设置观察方向/ZOOM!对图形显示窗口的某一区域进行缩放第三天生成关键点和线部分 1.生成关键点K,关键点编号,X坐标,Y坐标,Z坐标例：

ANSYS APDL命令流建模及模态分析实例相关内容

本文介绍了轮毂的ANSYS APDL命令流建模及模态分析实例相关内容。 ANSYS命令流及注释 五个辐条的轮毂 ! !初始化ANSYS环境 ! FINISH /CLEAR !清空内存 /FILNAM,WHEEL5 !文件名 /TITILE,WHEEL5 PARAMETER MODELING !工作名 ! !定义几何尺寸参数 ! R1=180 R2=157 R3=75 R4=75 R5=30 R6=28 R7=20 R8=90 R9=60 S_HOLE=5 TH1=48 TH2=23 TH3=11 TH4=180 TH5=40 TH6=45 TH7=105

TH8=25 TH9=15 TH10=25 TH11=13 /VIEW,1,1,1,1 !改变视图/ANG,1 /PNUM,LINE,1 /PNUM,AREA,1 /PNUM,VOLU,1 /NUMBER,1 ! !关键点 ! /PREP7 k,1,r5,r7,0 k,2,r4-ky(1),ky(1),0 k,3,r4,0,0 k,4,r1,0,0 k,5,kx(4),th5-th9,0 k,6,r1-th8,ky(5),0 k,7,kx(6),th4/2,0 k,8,kx(7)+th11,ky(7)+th10,0 k,9,kx(8),th4-th3,0 k,10,kx(4),ky(9),0 k,11,kx(4),th4,0 k,12,r2,ky(11),0 k,13,kx(12),ky(8),0 k,14,kx(7)-th3,ky(7),0 k,15,kx(14),th5,0 k,16,r3+r6,ky(15),0

ansys实例命令流-谱分析命令流

/FILNAME, Beam,1 ！定义工作文件名。 /TITLE, Beam Analysis ！定义工作标题。/PREP7 !定义单元。 ET,1,BEAM188 !定义材料属性。 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,2.1e5 MPDATA,PRXY,1,,0.3 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,7.9e-6 ! 定义杆件截面■200。 SECTYPE, 1, BEAM, RECT, , 0 SECOFFSET, CENT SECDATA,10,10,0,0,0,0,0,0,0,0 !建立几何模型。 K,1, ,, , K,2,350,, , !生成立柱。 LSTR, 1, 2 !以上完成几何模型。 !以下进行网格划分。 FLST,5,1,4,ORDE,1 FITEM,5,1 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,50, , , , ,1 !定义单元大小。!分配、划分平板结构。 LMESH, 1 !分析类型施加载荷并求解。 FINISH /SOLUTION ANTYPE,2 ！定义分析类型及求解设置。MSAVE,0 ！模态提取方法。

MODOPT,LANB,10 EQSLV,SPAR MXPAND,10, , ,0 ！模态扩展设置。 LUMPM,0 PSTRES,0 MODOPT,LANB,10,0,0, ,OFF MXPAND,10,0,0,1,0.001, !施加约束。 FLST,2,2,3,ORDE,2 FITEM,2,1 FITEM,2,-2 /GO DK,P51X, , , ,0,ALL, , , , , , !求解。 FINISH /SOL /STATUS,SOLU SOLVE !定义谱分析。 FINISH /SOLUTION ANTYPE,8 SPOPT,PSD,10,1 PSDUNIT,1,DISP,386.4, PSDFRQ,1, ,13.8,40,50.6,73,120 !定义谱—频率表。PSDFRQ,1, ,134,178,233, , PSDV AL,1,1,4,0.6,3,5 PSDV AL,1,6,2,6, , FLST,2,2,1,ORDE,2 !施加谱。 FITEM,2,1 FITEM,2,-2 D,P51X,UX,1.0 PFACT,1,BASE, !计算PSD激励参与系数。PSDRES,DISP,REL !设置输出选项。PSDRES,VELO,OFF PSDRES,ACEL,OFF

ANSYS的基本使用

2ANSYS的基本使用；2.1ANSYS环境简介；ANSYS有两种模式：一种是交互模式（Inter；运行该程序一般采用Interactive进入，这；进入系统后会有6个窗口，提供使用者与软件之间的交；各窗口的功能如下：；1.应用命令菜单（UtilityMenu）：包含；设定（WorkPlane）、参数化设计（Para；及辅助说明（Help）等；2.主菜单（M 2 ANSYS 的基本使用 2.1 ANSYS环境简介 ANSYS有两种模式：一种是交互模式（Interactive Mode），另一个是非交互模式（Batch Mode）。交互模式是初学者和大多数使用者所采用，包括建模、保存文件、打印图形及结果分析等，一般无特别原因皆用交互模式。但若分析的问题要很长时间，如一、两天等，可把分析问题的命令做成文件，利用它的非交互模式进行分析。 运行该程序一般采用 Interactive 进入，这样可以定义工作名称，并且存放到指定的工作目录中。若使用 Run Interactive Now 进入还需使用命令定义工作文件名或使用默认的文件名，使用该方式进入一般是为恢复上一次中断的分析。所以在开始分析一个问题时，建议使用 Interactive 进入交互模式。 进入系统后会有6个窗口，提供使用者与软件之间的交流，凭借这6个窗口可以非常容易的输入命令、检查模型的的建立、观察分析结果及图形输出与打印。整个窗口系统称为GUI(Graphical User Interface).如图2-1所示。 各窗口的功能如下： 1. 应用命令菜单（Utility Menu）：包含各种应用命令，如文件控制（Fi le）、对象选择（Select）、资料列式（List）、图形显示（Pplot）、图形控制（PlotCtrls）、工作界面

几个ansys经典实例(长见识)

平面问题斜支座的处理 如图5-7所示，为一个带斜支座的平面应力结构，其中位置2及3处为固定约束，位置4处为一个45o的斜支座，试用一个4节点矩形单元分析该结构的位移场。 (a)平面结构(b)有限元分析模型 图5-7 带斜支座的平面结构 基于ANSYS平台，分别采用约束方程以及局部坐标系的斜支座约束这两种方式来进行处理。 (7) 模型加约束 左边施加X,Y方向的位移约束 ANSYS Main Menu: Solution →Define Loads →Apply →-Structural→Displacement On Nodes →选取2，3号节点→OK →Lab2: All DOF(施加X,Y方向的位移约束) →OK 以下提供两种方法处理斜支座问题，使用时选择一种方法。 ?采用约束方程来处理斜支座 ANSYS Main Menu：Preprocessor →Coupling/ Ceqn →Constraint Eqn ：Const :0, NODE1:4, Lab1: UX,C1:1,NODE2:4,Lab2:UY,C2:1→OK 或者?采用斜支座的局部坐标来施加位移约束 ANSYS Utility Menu：WorkPlane →Local Coordinate System →Create local system →At specified LOC + →单击图形中的任意一点→OK →XC、YC、ZC分别设定为2，0，0,THXY：45 →OK ANSYS Main Menu：Preprocessor →modeling →Move / Modify →Rotate Node CS →To active CS → 选择4号节点 ANSYS Main Menu：Solution →Define Loads →Apply →Structural →Displacement On Nodes →选取4号节点→OK →选择Lab2:UY(施加Y方向的位移约束) →OK 命令流； !---方法1 begin----以下的一条命令为采用约束方程的方式对斜支座进行处理 CE,1,0,4,UX,1,4,UY,-1 !建立约束方程(No.1): 0=node4_UX*1+node_UY*(-1) !---方法1 end --- !--- 方法2 begin --以下三条命令为定义局部坐标系，进行旋转，施加位移约束 !local,11,0,2,0,0,45 !在4号节点建立局部坐标系 !nrotat, 4 !将4号节点坐标系旋转为与局部坐标系相同 !D,4,UY !在局部坐标下添加位移约束 !--- 方法2 end

ansys命令流最全详细介绍二

三 生成关键点和线部分 1.生成关键点 K,关键点编号,X坐标,Y坐标,Z坐标 例:K,1,0,0,0 2.在激活坐标系生成直线 LSTR,关键点P1,关键点P2 例LSTR,1,2 3.在两个关键点之间连线 L,关键点P1,关键点P2 例L,1,2 注:此命令会随当前的激活坐标系不同而生成直线或弧线 4.由三个关键点生成弧线 LARC,关键点P1,关键点P2,关键点PC,半径RAD 例LARC,1,3,2,0.05 注:关键点PC是用来控制弧线的凹向 5.通过圆心半径生成圆弧

CIRCLE,关键点圆心,半径RAD,,,,圆弧段数NSEG 例:CIRCLE,1,0.05,,,,4 6.通过关键点生成样条线 BSPLIN,关键点P1,关键点P2,关键点P3,关键点P4,关键点P5,关键点P6 例:BSPLIN,1,2,3,4,5,6 7.生成倒角线 LFILLT,线NL1,线NL2,倒角半径RAD 例LFILLT,1,2,0.005 8.通过关键点生成面 A,关键点P1,关键点P2,关键点P3,关键点P4,关键点P5,关键点P6,P7,P8... 例:A,1,2,3,4 9.通过线生成面 AL,线L1,线L2,线L3,线L4,线L5,线L6,线L7,线L8,线L9,线L10 例:AL,5,6,7,8 10.通过线的滑移生成面

ASKIN,线NL1,线NL2,线NL3,线NL4,线NL5,线NL6,线NL7,线NL8,线NL9 例:ASKIN,1,4,5,6,7,8 注:线1为滑移的导向线 四 目标:掌握常用的实体-面的生成 生成矩形面 1.通过矩形角上定位点生成面 BLC4,定位点X方向坐标XCORNER,定位点Y方向坐标YCORNER,矩形宽度WIDTH,矩形高度HEIGHT,矩形深度DEPTH 例:BLC4,0,0,5,3,0 2.通过矩形中心定位点生成面 BLC5,定位点X方向坐标XCENTER,定位点Y方向坐标YCENTER,矩形宽度WIDTH,矩形高度HEIGHT,矩形深度DEPTH 注:与上条命令的不同就在于矩形的定位点不一样 例:BLC5,2.5,1.5,5,3,0 3.通过在工作平面定义矩形X.Y坐标生成面 RECTNG,矩形左边界X坐标X1,矩形右边界X坐标X2,矩形下边界Y

ANSYS建模apdl命令流实例应用

大桥全长2996.8m，其中主桥采用跨度为101.5+188.5+580+217.5+159.5+116m的钢桁梁斜拉桥；非通航孔正桥采用6孔跨径64m预应力混凝土简支箱梁；东引桥采用16孔梁长32.6m预应力混凝土简支箱梁；跨大堤桥采用48.9+86+48.8m预应力混凝土连续箱梁；西引桥采用15孔梁长32.6m预应力混凝土简支梁及2孔梁长24.6m预应力混凝土简支梁，其中宁安线采用箱梁，阜景线采用T梁。 主桥采用103+188.5+580+217.5+159.5+117.5m两塔钢桁斜拉桥方案，全长1366m。主梁为三片主桁钢桁梁，桁间距2x14m，节间长14.5m，桁高15m。主塔为钢筋混凝土结构，塔顶高程+204.00m，塔底高程-6.00m，斜拉索为空间三索面，立面上每塔两侧共18对索，全桥216根斜拉索。所有桥墩上均设竖向和横向约束，4#塔与主梁之间设纵向水平约束，3#塔与梁间使用带限位功能的粘滞阻尼器。主梁为”N”字型桁式，横向采用三片桁结构，主桁的横向中心距各为14m，桁高15m，节间距14.5m[2]。 结构构造 主桥采用两塔钢桁斜拉桥方案，主梁为三片主桁钢桁梁，主桁上下弦杆均为箱型截面，上弦杆内高1000mm，内宽1200mm，板厚20～48mm。下弦杆内高1400mm，宽1200mm，板厚20～56mm。下弦杆顶板向桁内侧加宽700mm与整体桥面板焊接。腹杆主要采用H型截面。H型杆件宽1200mm，高720和760mm，板厚20～48mm。根据不同的受力区段选用不同的杆件截面，在辅助墩附近的压重区梁段，腹杆采用箱型截面杆件。主桁采用焊接杆件，整体节点。在节点外以高强度螺栓拼接的结构形式，上下弦杆四面等强对接拼装。H型腹杆采用插入式连接。箱型腹杆采用四面与主桁节点对拼的连接形式。主桁拼接采用M30高强螺栓。

ANSYS带轮建模命令流

Finish$/clear$/prep7 Dd=200$fai=38$b=13$z=4 S=14$ha=3 Hf=9$e=15$f=10 Ks=8$r1=0.5$r2=1.0 R3=1.5$dta=6 C1=2$c2=2 Pd=25$d0=24 D1=1.9*d0 S1=1.5*s S2=0.5*s L=2*d0 Da=dd+2*ha Ub=(z-1)*e+2*f Rk1=dd/2-hf-dta-0.5*(ub-s)/pd-s2 Rk2=d1/2+0.5*(l-s)/pd+s1 Rk=(rk1-rk2)/2 Dk=rk1+rk2 *afun,deg$y0=hf+ha B0=b-2*tan(fai/2)*y0 Local,12,0,-ub/2,dd/2-hf K,,0,y0$k,,f-b/2,y0$k,,f-b0/2 *do,i,1,z-1 X0=f+(i-1)*e$k,,x0+b0/2 K,,x0+b/2,y0$k,,x0+e-b/2,yo K,,x0+e-b0/2$*enddo K,,ub-f+b0/2$k,,ub-f+b/2,y0 K,,ub,y0 *get,kp1,kp,0,num,max *do,i,1,kp1-1$l,i,i+1$*enddo

*get,l1,line,0,num,max *do,i,1,z$j=4*i Lfillt,j,j+1,r1$Lfillt,j+1,j+2,r2 Lfillt,j+3,j+4,r2$*enddo Lfillt,3,4,r2 Csdele,12$ksll,s Ksel,inve$kdele,all Allsel Numcmp,all Cm,l1cm,line *get,kp1,kp,0,num,max Y0=dd/2-hf-dta$k,,-ub/2,y0+c1-c1/pd$k,,-ub/2+c1,y0-c1/pd K,,-s/2,y0-0.5*(ub-s)/pd$k,,-s/2,d1/2+0.5*(l-s)/pd K,,-l/2,d1/2$k,,-l/2,d0/2+c2$k,,-l/2+c2,d0/2$l,1,kp1+1 *do,i,kp1+1,kp1+6$l,i,i+1$*enddo Lsel,s,loc,y,y0-c1/pd,d1/2+0.5*(l-s)/pd$*get,l1,line,0,num,min L2=lsnext(l1)$lfillt,l1,l2,r3$lsel,all$lfillt,l2,l2+1,r3 Cmsel,u,l1cm Lsymm,x,all$ksel,s,loc,y,d0/2$*get,kp1,kp,0,num,min Kp2=kpnext(kp1)$l,kp1,kp2$allsel Nummrg,all$numcmp,all Al,all *get,kp1,kp,0,num,max K,kp1+10,-ub/2$k,kp1+20,ub/2 Vrotat,all,,,,,,kp1+10,kp1+20,,ks Kdele,kp1+10,kp1+20,10 !chouchou

ansys实例命令流-实体梁分析命令流

/FILNAME,SolidBeam ,1 !定义工作文件名。/TITLE,SolidBeam Analysis !定义工作标题。/PREP7 ET,1,SOLID95 !定义材料属性。 MPTEMP,,,,,,,, !定义材料属性。MPTEMP,1,0 MPDATA,EX,1,,2.06e5 MPDATA,PRXY,1,,0.3 !建立几何模型 K,1,,,, K,2,450,,, K,3,450,-55,, K,4,,-100,, FLST,2,4,3 FITEM,2,1 FITEM,2,2 FITEM,2,3 FITEM,2,4 A,P51X VOFFST,1,45, , !网格划分。 FLST,5,4,4,ORDE,2 FITEM,5,9 FITEM,5,-12 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,3, , , , ,1 !* FLST,5,4,4,ORDE,4

FITEM,5,2 FITEM,5,4 FITEM,5,6 FITEM,5,8 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,6, , , , ,1 !* FLST,5,4,4,ORDE,4 FITEM,5,1 FITEM,5,3 FITEM,5,5 FITEM,5,7 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,30, , , , ,1 !* CM,_Y,VOLU VSEL, , , , 1 CM,_Y1,VOLU CHKMSH,'VOLU' CMSEL,S,_Y !* VSWEEP,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 !加载。

ANSYS_结构稳态(静力)分析之经典实例-命令流格式

/FILNAME,Allen-wrench,1 ! Jobname to use for all subsequent files /TITLE,Static analysis of an Allen wrench /UNITS,SI ! Reminder that the SI system of units is used /SHOW ! Specify graphics driver for interactive run; for batch ! run plots are written to pm02.grph ! Define parameters for future use EXX=2.07E11 ! Young's modulus (2.07E11 Pa = 30E6 psi) W_HEX=.01 ! Width of hex across flats (.01m=.39in) *AFUN,DEG ! Units for angular parametric functions定义弧度单位 W_FLAT=W_HEX*TAN(30) ! Width of flat L_SHANK=.075 ! Length of shank (short end) (.075m=3.0in) L_HANDLE=.2 ! Length of handle (long end) (.2m=7.9 in) BENDRAD=.01 ! Bend radius of Allen wrench (.01m=.39 in) L_ELEM=.0075 ! Element length (.0075 m = .30 in) NO_D_HEX=2 ! Number of divisions on hex flat TOL=25E-6 ! Tolerance for selecting nodes (25e-6 m = .001 in) /PREP7 ET,1,SOLID45 ! 3维实体结构单元；Eight-node brick element ET,2,PLANE42 ! 2维平面结构；Four-node quadrilateral (for area mesh) MP,EX,1,EXX ! Young's modulus for material 1；杨氏模量 MP,PRXY,1,0.3 ! Poisson's ratio for material 1；泊松比 RPOLY,6,W_FLAT ! Hexagonal area创建规则的多边形 K,7 ! Keypoint at (0,0,0) K,8,,,-L_SHANK ! Keypoint at shank-handle intersection K,9,,L_HANDLE,-L_SHANK ! Keypoint at end of handle L,4,1 ! Line through middle of hex shape L,7,8 ! Line along middle of shank L,8,9 ! Line along handle LFILLT,8,9,BENDRAD ! Line along bend radius between shank and handle! 产生一个倒角圆，并生成三个点 /VIEW,,1,1,1 ! Isometric view in window 1 /ANGLE,,90,XM ! Rotates model 90 degrees about X! 不用累积的旋转 /TRIAD,ltop /PNUM,LINE,1 ! Line numbers turned on LPLOT ! Line numbers off !