L=0.1
L1=0.04
H1=0.02
H=0.005
DERTA=0.95*H
L2=L1-2*DERTA
H2=H1
/PREP7
K,1,0,0,0
K,2,L,0,0
K,3,L,-H1,0
K,4,L+L1,-H1,0
K,5,L+L1,0,0
K,6,2*L+L1,0,0
K,7,0,H2+H,0
K,8,L+DERTA,H2+H,0
K,9,L+DERTA,H,0
K,10,L+DERTA+L2,H,0
K,11,L+DERTA+L2,H+H2,0
K,12,2*L+L1,H+H2,0
LSTR, 1, 2
LSTR, 2, 3
LSTR, 3, 4
LSTR, 4, 5
LSTR, 5, 6
LSTR, 7, 8
LSTR, 8, 9
LSTR, 9, 10
LSTR, 10, 11
LSTR, 11, 12
LFILLT,6,7,0.004, ,
!*
LFILLT,9,10,0.004, ,
!*
LFILLT,7,8,0.004, ,
!*
LFILLT,8,9,0.004, ,
LFILLT,1,2,0.004, ,
!*
LFILLT,4,5,0.004, ,
!*
LFILLT,4,3,0.007, ,
!*
LFILLT,3,2,0.007, ,
RECTNG,0.02,0.22,0,0.005,
ET,1,PLANE82
!*
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,1E11
MPDATA,PRXY,1,,0.3
TB,BISO,1,1,2,
TBTEMP,0
TBDATA,,2E8,2E9,,,,
ESIZE,0.005/3,0,
MSHAPE,0,2D
MSHKEY,1
!*
CM,_Y,AREA
ASEL, , , , 1
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
!*
AMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,0.3
MAT,1
R,3
REAL,3
ET,2,169
ET,3,172
KEYOPT,3,9,0
KEYOPT,3,10,2
R,3,
RMORE,
RMORE,,0
RMORE,0
! Generate the target surface
LSEL,S,,,6
LSEL,A,,,7
LSEL,A,,,8
LSEL,A,,,9
LSEL,A,,,10
LSEL,A,,,11
LSEL,A,,,12
LSEL,A,,,13
LSEL,A,,,14
CM,_TARGET,LINE
TYPE,2
LATT,-1,3,2,-1
TYPE,2
LMESH,ALL
! Create a pilot node
KSEL,S,,,7
KATT,-1,3,2,-1
KMESH,7
! Generate the contact surface
LSEL,S,,,21
CM,_CONTACT,LINE
TYPE,3
NSLL,S,1
ESLN,S,0
ESURF
*SET,_REALID,3
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
ESEL,R,REAL,,3
LSEL,S,REAL,,3
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
ESEL,R,REAL,,3
LSEL,S,REAL,,3
CMSEL,A,_NODECM
CMDEL,_NODECM
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM
CMSEL,S,_KPCM
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM
CMSEL,S,_AREACM
CMDEL,_AREACM
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
!*
!*
!*
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,0.3
MAT,1
R,4
REAL,4
ET,4,169
ET,5,172
KEYOPT,5,9,0
KEYOPT,5,10,2
R,4,
RMORE,
RMORE,,0
RMORE,0
! Generate the target surface
LSEL,S,,,1
LSEL,A,,,2
LSEL,A,,,3
LSEL,A,,,4
LSEL,A,,,5
LSEL,A,,,15
LSEL,A,,,16
LSEL,A,,,17
LSEL,A,,,18
CM,_TARGET,LINE
TYPE,4
LATT,-1,4,4,-1
TYPE,4
LMESH,ALL
! Create a pilot node
KSEL,S,,,1
KATT,-1,4,4,-1
KMESH,1
! Generate the contact surface
LSEL,S,,,19
CM,_CONTACT,LINE
TYPE,5
NSLL,S,1
ESLN,S,0
ESURF
*SET,_REALID,4
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,4
ESEL,A,TYPE,,5
ESEL,R,REAL,,4
LSEL,S,REAL,,4
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
! Reverse target normals
FLST,5,9,4,ORDE,4
FITEM,5,1
FITEM,5,-5
FITEM,5,15
FITEM,5,-18
CM,_Y,LINE
LSEL, , , ,P51X
CM,_YEL,ELEM
CM,_YND,NODE
NSLL,S,1
ESLN,S,1
ESEL,R,REAL,,_REALID
ESURF,,REVERSE
CMSEL,S,_Y
CMSEL,S,_YEL
CMSEL,S,_YND
CMDELE,_Y
CMDELE,_YEL
CMDELE,_YND
/REPLOT
!*
ESEL,ALL
ESEL,S,TYPE,,4
ESEL,A,TYPE,,5
ESEL,R,REAL,,4
LSEL,S,REAL,,4
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,4
ESEL,A,TYPE,,5
ESEL,R,REAL,,4
LSEL,S,REAL,,4
CMSEL,A,_NODECM
CMDEL,_NODECM
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM
CMSEL,S,_KPCM
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM
CMSEL,S,_AREACM
CMDEL,_AREACM
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
CM,_CWZ_EL,ELEM
CM,_CWZ_ND,NODE
CM,_CWZ_KP,KP
CM,_CWZ_LN,LINE
CM,_CWZ_AR,AREA
CM,_CWZ_VL,VOLU
ESEL,NONE
ESEL,A,REAL,,4
ESEL,R,ENAME,,169,170
NSLE
*GET,_z1,ELEM,,NUM,MAX
KSLN,S
LSLK,S,1
ASLL,S,1
*CREATE,cwzplot,mac
~eui,'::apdl::noprint 1'
~eui,'::apdl::nooutput 1'
/PNUM,REAL,1
/NUM,1
/PSYMB,ESYS,0
EPLOT
/PSYMB,ESYS,0
/NUM,0
/PNUM,TYPE,0
/PNUM,REAL,0
/mrep,cwzplot
~eui,'::apdl::nooutput 0'
~eui,'::apdl::noprint 0'
*END
cwzplot
*SET,_REALID,4
FLST,5,2,2,ORDE,2
FITEM,5,498
FITEM,5,-499
CM,_ELMCM,ELEM
ESEL,S, , ,P51X
ESURF,,REVERSE
/REPLOT
CMSEL,S,_ELMCM
CMDELE,_ELMCM
!*
CMSEL,S,_CWZ_EL
CMDEL,_CWZ_EL
CMSEL,S,_CWZ_ND
CMDEL,_CWZ_ND
CMSEL,S,_CWZ_KP
CMDEL,_CWZ_KP
CMSEL,S,_CWZ_LN
CMDEL,_CWZ_LN
CMSEL,S,_CWZ_AR
CMDEL,_CWZ_AR
CMSEL,S,_CWZ_VL
CMDEL,_CWZ_VL
!*
/MREP,EPLOT
NROPT,UNSYM
FINISH
/SOL
EPLOT
ANTYPE,0
NLGEOM,1
!*
FLST,2,1,3,ORDE,1
FITEM,2,1
!*
/GO
DK,P51X, , , ,0,ALL, , , , , ,
FLST,2,1,3,ORDE,1
FITEM,2,7
!*
/GO
DK,P51X, , , ,0,UX,ROTZ, , , , ,
FLST,2,1,3,ORDE,1
FITEM,2,7
!*
/GO
DK,P51X, ,-1.1*h2, ,0,UY, , , , , ,
!*
OUTRES,ALL,ALL,
!*
TIME,1
AUTOTS,-1
DELTIM,0.0006,0.0005,0.002,1
KBC,0
!*
TSRES,ERASE
lswrite,1
FLST,2,1,3,ORDE,1
FITEM,2,7
!*
/GO
DK,P51X, ,0.2*h2, ,0,UY, , , , , ,
!*
OUTRES,ALL,ALL,
!*
TIME,1
AUTOTS,-1
DELTIM,0.001,0.0005,0.002,1
KBC,0
!*
TSRES,ERASE
lswrite,2,
lssolve,1,2,1
[ 本帖最后由 liangdiansheng 于 2008-11-9 14:56 编辑 ] |