内容发布更新时间 : 2024/5/5 9:05:14星期一 下面是文章的全部内容请认真阅读。

6.3.3 命令流模式

/PREP7 ！定义标题

/TITLE, TRANSVERSE SHEAR STRESSES IN A CANTILEVER BEAM ANTYPE,STATIC

ET,1,SHELL99,,,,,2,4 ! 8-NODE 壳单元 PRINTOUT

R,1,4,1 ! 四层轴对称壳 RMORE

RMORE,1,,0.5,1,,0.5 ! 等厚度 MP,EX,1,30E6 MP,NUXY,1,0

/COM --- INPUT FAILURE STRESSES FOR MATERIAL #1 --- /COM --- COMPRESSION VALUES ARE LEFT TO DEFAULT --- TB,FAIL,1 TBTEMP,,CRIT

TBDATA,1,0,0,1 ! 使用TSAI-WU破坏准则 TBTEMP,0 ！定义破坏应力 TBDATA,10,25000 TBDATA,12,3000 TBDATA,14,5000 TBDATA,16,500 N,1 N,3,,1 FILL

NGEN,11,3,1,3,1,1 E,1,7,9,3,4,8,6,2 EGEN,5,6,-1

NSEL,S,LOC,X ！选择节点 D,ALL,ALL NSEL,S,LOC,X,10

CP,1,UZ,ALL ！耦合节点 NSEL, R, LOC, Y

F,ALL,FZ,10000 ！施加集中力 NSEL,ALL OUTPR,,1 FINISH /SOLU SOLVE FINISH ！后处理 /POST1

ETABLE,NX,SMISC,7 ETABLE,FC,NMISC,1 ETABLE,FCMX,NMISC,2 ETABLE,FCLN,NMISC,3 ETABLE,ILMX,NMISC,4

ETABLE,ILLN,NMISC,5 PRETAB,NX,ILLN,ILMX PRETAB,FC,FCLN,FCMX ETABLE,SXZ,S,XZ ETABLE,ILSXZ,SMISC,11 ！获取单元结果数据

*GET,SIGXZ1,ELEM,4,ETAB,SXZ *GET,SIGXZ2,ELEM,1,ETAB,ILSXZ *GET,SIGXZ3,ELEM,1,ETAB,ILMX *GET,FC3,ELEM,1,ETAB,FCMX *DIM,LABEL,CHAR,4,2 *DIM,VALUE,,4,3

LABEL(1,1) = 'SIGXZ,ps','SIGXZ,ps','SIGXZ,ps','FC3MAX (' LABEL(1,2) = 'i(Z=H/2)','i(Z=H/4)','i(Z= 0 )','FCMX) ' *VFILL,VALUE(1,1),DATA,0,5625,7500,225

*VFILL,VALUE(1,2),DATA,SIGXZ1,SIGXZ2,SIGXZ3,FC3

*VFILL,VALUE(1,3),DATA,0,ABS(SIGXZ2/5625),ABS(SIGXZ3/7500),ABS(FC3/225) /OUT,beam,vrt FINISH

7.3.3 命令流模式

/PREP7

/TITLE, NATURAL FREQUENCY OF A PIEZOELECTRIC TRANSDUCER ET,1,SOLID5,3 ！选择单元 ！设定材料属性 MP,DENS,3,7500 MP,PERX,3,804.6 MP,PERZ,3,659.7 ！定义表格 TB,PIEZ,3 TBDATA,16,10.5 TBDATA,14,10.5 TBDATA,3,-4.1 TBDATA,6,-4.1 TBDATA,9,14.1 TB,ANEL,3

TBDATA,1,13.2E10,7.1E10,7.3E10 TBDATA,7,13.2E10,7.3E10 TBDATA,12,11.5E10 TBDATA,16,3.0E10 TBDATA,19,2.6E10 TBDATA,21,2.6E10 ！建模 L=10E-3 W=10E-3 H=20E-3 K,1 K,2,L

K,3,L,W K,4,,W

KGEN,2,1,4,1,,,H L,1,5

LESIZE,1,,,4 ！设定单元尺寸 ESIZE,,2

MSHK,1 MSHA,0,3D V,1,2,3,4,5,6,7,8

VMESH,1 ！划分单元 NSEL,S,LOC,X,0 ！施加约束

DSYM,SYMM,X NSEL,S,LOC,Y,0 DSYM,SYMM,Y NSEL,ALL FINISH /SOLU ANTYP,MODAL

MODOPT,LANB,10,50000,150000 ！设置提取模态方式 MXPAND,10 ！扩展模态 /COM, NSEL,S,LOC,Z,0 D,ALL,VOLT,0 NSEL,S,LOC,Z,H

D,ALL,VOLT,0 NSEL,ALL SOLVE

*GET,F1,MODE,2,FREQ *GET,F2,MODE,5,FREQ FINISH /POST1 /VIEW,1,1,1,1 /VUP,,Z /TYPE,,4

*DO,I,1,10 SET,LS,I

PLDISP,1 ！显示模态 *ENDDO FINISH /SOLU ANTYPE,MODAL

MODOPT,LANB,10,50000,150000 MXPAND,10

/COM, NSEL,S,LOC,Z,H

DDELE,ALL,VOLT ！删除VOLT 约束 CP,1,VOLT,ALL ！耦合VOLT 约束 NSEL,ALL