程序
/PREP7
ET,1,SOLID186
MPTEMP,,,,,,,, !钛
MPTEMP,1,0
MPDATA,EX,1,,1.15e11
MPDATA,PRXY,1,,0.37
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,4.5e3
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,MU,1,,0.05
TB,BISO,1,1,2, !双线性等向强化材料模型
TBTEMP,0
TBDATA,1,1000E6,4.2E10,,,,
MPTEMP,,,,,,,, !金刚石
MPTEMP,1,0
MPDATA,EX,2,,1.1e12
MPDATA,PRXY,2,,0.07
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,2,,3.52e3
/PNUM,AREA,1
/PNUM,VOLU,1
CYLIND,0.5e-3,0,0,0.6e-3,0,90,
/ANG,1,30,YS,1
/REP,FAST
/ANG,1,30,YS,1
/REP,FAST
/ANG,1,30,YS,1
/REP,FAST
CYLIND,0.28e-3,0,0,0.2e-3,0,90,
VSBV, 1, 2
CONE,0.28e-3,0.3e-3,0.2e-3,0.3923e-3,0,90,
VSBV, 3, 1
CONE,0.3e-3,0.4e-3,0.3923e-3,0.6e-3,0,90,
VSBV, 2, 1
CYLIND,0.28e-3,0,-0.5e-3,0.2e-3,0,90,
CONE,0.28e-3,0.3e-3,0.2e-3,0.3923e-3,0,90,
CONE,0.3e-3,0.3e-3,0.3923e-3,1e-3,0,90,
VADD,1,2,4
SMRT,4
VSEL, , , , 3
VSWEEP,3
VSEL, , , , 5
VSWEEP,5
R,3,,,0.1
ET,2,TARGE170
ET,3,CONTA174
ASEL,S,,,8 !创建目标面
ASEL,A,,,10
ASEL,A,,,13
REAL,3
MAT,2
TYPE,2
NSLA,S,1
ESLN,S,0
ESURF
ALLSEL,ALL
ASEL,S,,,4 !创建接触面
ASEL,A,,,24
ASEL,A,,,27
REAL,3
MAT,1
TYPE,3
NSLA,S,1
ESLN,S,0
ESURF
ALLSEL,ALL
EPLOT
REAL,3 !创建目标面的控制点
MAT,2
TYPE,2
KPLOT
KMESH, 9
ALLSEL,ALL
eplot
FINISH
/SOLU
DA,16,SYMM
DA,17,SYMM
DA,5,SYMM
DA,6,SYMM
DA,25,SYMM
DA,26,SYMM
DA,28,SYMM
DA,29,SYMM
ANTYPE,4
TRNOPT,FULL
tintp,,0.25,0.5,0.5
outpr,all,all
OUTRES,ALL,ALL
AUTOTS,ON
pred,on
LNSRCH,on
KBC,1
nlgeom,on
ff=20945
pi=3.1415926
t=1/ff
n=30
*do,i,1,n*20
time,t/n*i
degree=sqrt(2)*5/180*3.1415926*sin(2*pi*ff*t/n*i)
zong=sqrt(2)*(70e-9)*sin(2*pi*ff*t/n*i)
DK,9,rotz,degree
DK,9,uz,zong
allsel,all
da,1,uz,-1.5873e-5
allsel,all
allsel
nsubt,10,10000,1
solve
*enddo
finish
|