[转帖]一个BRACKET的接触分析
这两天做了一个接触分析，是一个支架用螺栓固定在一立板上的大变形分析，采用双线性随动强化材料模型，考虑螺栓与支架间的接触以及支架与立板间的接触关系，三个螺栓分别施加5000N的预紧力，在支架的平板上施加1MPa的面力。为了考虑螺栓的预紧与锁定，将螺栓预紧力的施加和边界约束条件的施加放在同一个载荷步计算，在第二个载荷步才施加平板上的面力。命令流文件如下：/PREP7
ET,1,SOLID45
EX,1,2.1E11
PRXY,1,0.3
TB,BISO,1,1,2,
TBTEMP,0
TBDATA,,235E6,0.4E11,,,,
!建立初始几何模型
BLOCK,-0.030,0.030,-0.060,0.040,,0.005,
BLOCK,-0.03,0.03,0,0.005,,0.1,
FLST,2,2,6,ORDE,2
FITEM,2,1
FITEM,2,-2
VOVLAP,P51X
wpoff,0.010,,0.005
CSYS,4
K, ,,,,
K, ,,,0.045,
K, ,,-0.030,0,
FLST,2,3,3
FITEM,2,21
FITEM,2,22
FITEM,2,23
A,P51X
FLST,2,1,5,ORDE,1
FITEM,2,1
VEXT,P51X, , ,-0.005,0,0,,,,
CSYS,0
FLST,3,1,6,ORDE,1
FITEM,3,1
VSYMM,X,P51X, , , ,0,0
FLST,3,1,6,ORDE,1
FITEM,3,1
VGEN, ,P51X, , ,0.005, , , , ,1
FLST,3,1,6,ORDE,1
FITEM,3,2
VGEN, ,P51X, , ,-0.005, , , , ,1
!用工作平面分割几何体
KWPAVE, 31
FLST,2,6,6,ORDE,2
FITEM,2,1
FITEM,2,-6
V***W,P51X,,DELETE !将V***W修改为V S B W(中间无空格）
wpro,,,-90
FLST,2,7,6,ORDE,4
FITEM,2,1
FITEM,2,-5
FITEM,2,7
FITEM,2,-8
V***W,P51X,,DELETE !将V***W修改为V S B W(中间无空格）
/REPLO
KWPAVE, 28
FLST,2,12,6,ORDE,5
FITEM,2,1
FITEM,2,-2
FITEM,2,6
FITEM,2,9
FITEM,2,-17
V***W,P51X,,DELETE !将V***W修改为V S B W(中间无空格）
KWPAVE, 25
FLST,2,17,6,ORDE,8
FITEM,2,1
FITEM,2,-5
FITEM,2,7
FITEM,2,-10
FITEM,2,13
FITEM,2,-14
FITEM,2,17
FITEM,2,-22
V***W,P51X,,DELETE !将V***W修改为V S B W(中间无空格）
KWPAVE, 22
FLST,2,22,6,ORDE,8
FITEM,2,1
FITEM,2,-8
FITEM,2,11
FITEM,2,-12
FITEM,2,15
FITEM,2,-16
FITEM,2,18
FITEM,2,-27
V***W,P51X,,DELETE !将V***W修改为V S B W(中间无空格）
/REPLO
KWPAVE, 29
wpro,,-90
FLST,2,27,6,ORDE,8
FITEM,2,1
FITEM,2,-5
FITEM,2,7
FITEM,2,-10
FITEM,2,13
FITEM,2,-14
FITEM,2,17
FITEM,2,-32
V***W,P51X,,DELETE !将V***W修改为V S B W(中间无空格）
!工作平面移到螺栓孔的位置
FLST,2,2,3,ORDE,2
FITEM,2,8
FITEM,2,54
KWPAVE,P51X
wpoff,,,-0.0175
wpro,,,90.00
CSYS,4
CYLIND,0.004, ,-0.005,0.001,0,360,
V***V, 4, 3 !将V***V修改为V S B V(中间无空格）
CYLIND,0.003, ,-0.010,,0,360,
CYLIND,0.005, ,0,0.003,0,360,
FLST,2,2,3,ORDE,2
FITEM,2,102
FITEM,2,-103
KWPAVE,P51X
CYLIND,0.005, ,-0.003,,0,360,
CSYS,0
VDELE, 10, , ,1
!映射到另一边去
FLST,3,4,6,ORDE,4
FITEM,3,3
FITEM,3,-4
FITEM,3,9
FITEM,3,18
VSYMM,X,P51X, , , ,0,0
!工作平面移到下一个螺栓孔的位置
FLST,2,2,3,ORDE,2
FITEM,2,26
FITEM,2,32
KWPAVE,P51X
WPOFF,0.015
!产生下面的螺栓孔和螺栓
CSYS,4
CYLIND,0.004, ,-0.005,,0,360,
V***V, 34, 39 !将V***V修改为V S B V(中间无空格）
CYLIND,0.003, ,-0.010,,0,360,
CYLIND,0.005, ,0,0.003,0,360,
wpoff,,,-0.010
CYLIND,0.005, ,-0.003,0,0,360,
CSYS,0
!MERGE支架主体结构的KP
NUMMRG,KP, , , ,LOW
!把螺母螺栓GLUE起来
FLST,2,9,6,ORDE,9
FITEM,2,3
FITEM,2,-4
FITEM,2,10
FITEM,2,18
FITEM,2,23
FITEM,2,34
FITEM,2,38
FITEM,2,-39
FITEM,2,41
VGLUE,P51X
!移开螺栓
FLST,3,9,6,ORDE,5
FITEM,3,3
FITEM,3,10
FITEM,3,34
FITEM,3,42
FITEM,3,-47
VGEN, ,P51X, , , , ,-0.1, , ,1
!划分螺母螺栓单元
MSHAPE,1,3D
MSHKEY,0
!*
FLST,5,9,6,ORDE,5
FITEM,5,3
FITEM,5,10
FITEM,5,34
FITEM,5,42
FITEM,5,-47
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!划分支架
FLST,5,32,6,ORDE,13
FITEM,5,1
FITEM,5,-2
FITEM,5,5
FITEM,5,-9
FITEM,5,11
FITEM,5,-17
FITEM,5,19
FITEM,5,-22
FITEM,5,24
FITEM,5,-33
FITEM,5,35
FITEM,5,-37
FITEM,5,40
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VSWEEP,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL !由支架板映射生成固定板实体和网格
FLST,3,20,6,ORDE,16
FITEM,3,4
FITEM,3,-5
FITEM,3,8
FITEM,3,10
FITEM,3,-12
FITEM,3,14
FITEM,3,-15
FITEM,3,17
FITEM,3,-18
FITEM,3,21
FITEM,3,-22
FITEM,3,25
FITEM,3,-27
FITEM,3,30
FITEM,3,32
FITEM,3,-35
VSYMM,Z,P51X, , , ,0,0
!将螺栓螺母移到指定位置
FLST,3,9,6,ORDE,5
FITEM,3,3
FITEM,3,9
FITEM,3,31
FITEM,3,36
FITEM,3,-41
FLST,3,9,6,ORDE,5
FITEM,3,3
FITEM,3,9
FITEM,3,31
FITEM,3,36
FITEM,3,-41
VGEN, ,P51X, , , , ,0.1, , ,1
!建立接触对
/COM, CONTACT PAIR CREATION - START
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,
MAT,1
MP,EMIS,1,7.88860905221e-031
R,3
REAL,3
ET,2,170
ET,3,174
R,3,,,1.0,0.1,0,-0.0001
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,3,4,0
KEYOPT,3,5,0
KEYOPT,3,7,0
KEYOPT,3,8,0
KEYOPT,3,9,0
KEYOPT,3,10,1
KEYOPT,3,11,0
KEYOPT,3,12,0
KEYOPT,3,2,0
KEYOPT,2,5,0
! Generate the target surface
ASEL,S,,,132
CM,_TARGET,AREA
TYPE,2
NSLA,S,1
ESLN,S,0
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,42
CM,_CONTACT,AREA
TYPE,3
NSLA,S,1
ESLN,S,0
ESURF
ALLS
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
ESEL,R,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
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
/COM, CONTACT PAIR CREATION - END
VPLOT
/COM, CONTACT PAIR CREATION - START
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
MAT,1
MP,EMIS,1,7.88860905221e-031
R,4
REAL,4
ET,4,170
ET,5,174
R,4,,,1.0,0.1,0,-0.0001
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,5,4,0
KEYOPT,5,5,0
KEYOPT,5,7,0
KEYOPT,5,8,0
KEYOPT,5,9,0
KEYOPT,5,10,1
KEYOPT,5,11,0
KEYOPT,5,12,0
KEYOPT,5,2,0
KEYOPT,4,5,0
! Generate the target surface
ASEL,S,,,91
CM,_TARGET,AREA
TYPE,4
NSLA,S,1
ESLN,S,0
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,194
CM,_CONTACT,AREA
TYPE,5
NSLA,S,1
ESLN,S,0
ESURF
ALLS
ESEL,ALL
ESEL,S,TYPE,,4
ESEL,A,TYPE,,5
ESEL,R,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
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
/COM, CONTACT PAIR CREATION - END
VPLOT
/COM, CONTACT PAIR CREATION - START
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
MAT,1
MP,EMIS,1,7.88860905221e-031
R,5
REAL,5
ET,6,170
ET,7,174
R,5,,,1.0,0.1,0,-0.0001
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,7,4,0
KEYOPT,7,5,0
KEYOPT,7,7,0
KEYOPT,7,8,0
KEYOPT,7,9,0
KEYOPT,7,10,1
KEYOPT,7,11,0
KEYOPT,7,12,0
KEYOPT,7,2,0
KEYOPT,6,5,0
! Generate the target surface
ASEL,S,,,154
CM,_TARGET,AREA
TYPE,6
NSLA,S,1
ESLN,S,0
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,159
CM,_CONTACT,AREA
TYPE,7
NSLA,S,1
ESLN,S,0
ESURF
ALLS
ESEL,ALL
ESEL,S,TYPE,,6
ESEL,A,TYPE,,7
ESEL,R,REAL,,5
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,6
ESEL,A,TYPE,,7
ESEL,R,REAL,,5
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
/COM, CONTACT PAIR CREATION - END
VPLOT /COM, CONTACT PAIR CREATION - START
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
MAT,1
MP,EMIS,1,7.88860905221e-031
R,6
REAL,6
ET,8,170
ET,9,174
R,6,,,1.0,0.1,0,-0.0001
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,9,4,0
KEYOPT,9,5,0
KEYOPT,9,7,0
KEYOPT,9,8,0
KEYOPT,9,9,0
KEYOPT,9,10,1
KEYOPT,9,11,0
KEYOPT,9,12,0
KEYOPT,9,2,0
KEYOPT,8,5,0
! Generate the target surface
ASEL,S,,,142
CM,_TARGET,AREA
TYPE,8
NSLA,S,1
ESLN,S,0
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,244
CM,_CONTACT,AREA
TYPE,9
NSLA,S,1
ESLN,S,0
ESURF
ALLS
ESEL,ALL
ESEL,S,TYPE,,8
ESEL,A,TYPE,,9
ESEL,R,REAL,,6
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
VPLOT
*CREATE,cwzplot,mac
/COM,
/COM,PLOT CONTACT PAIR(S)
~eui,'::apdl::noprint 1'
~eui,'::apdl::nooutput 1'
!*
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,,6
ESEL,R,ENAME,,169,175
NSLE
KSLN,S
LSLK,S,1
ASLL,S,1
/PNUM,TYPE,1
/NUM,1
/PSYMB,ESYS,0
EPLOT
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
!*
/PSYMB,ESYS,0
/NUM,0
/PNUM,TYPE,0
/PNUM,REAL,0
/mrep,cwzplot
~eui,'::apdl::nooutput 0'
~eui,'::apdl::noprint 0'
*END
cwzplot
VPLOT
*CREATE,cwzplot,mac
/COM,
/COM,PLOT CONTACT PAIR(S)
~eui,'::apdl::noprint 1'
~eui,'::apdl::nooutput 1'
!*
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,,5
ESEL,R,ENAME,,169,175
NSLE
KSLN,S
LSLK,S,1
ASLL,S,1
/PNUM,TYPE,1
/NUM,1
/PSYMB,ESYS,0
EPLOT
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
!*
/PSYMB,ESYS,0
/NUM,0
/PNUM,TYPE,0
/PNUM,REAL,0
/mrep,cwzplot
~eui,'::apdl::nooutput 0'
~eui,'::apdl::noprint 0'
*END
cwzplot
VPLOT
*CREATE,cwzplot,mac
/COM,
/COM,PLOT CONTACT PAIR(S)
~eui,'::apdl::noprint 1'
~eui,'::apdl::nooutput 1'
!*
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,,5
ESEL,R,ENAME,,169,175
NSLE
KSLN,S
LSLK,S,1
ASLL,S,1
/PNUM,TYPE,1
/NUM,1
/PSYMB,ESYS,0
EPLOT
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
!*
/PSYMB,ESYS,0
/NUM,0
/PNUM,TYPE,0
/PNUM,REAL,0
/mrep,cwzplot
~eui,'::apdl::nooutput 0'
~eui,'::apdl::noprint 0'
*END cwzplot
/COM, CONTACT PAIR CREATION - START
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
MAT,1
MP,EMIS,1,7.88860905221e-031
R,7
REAL,7
ET,10,170
ET,11,174
R,7,,,1.0,0.1,0,-0.0001
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,11,4,0
KEYOPT,11,5,0
KEYOPT,11,7,0
KEYOPT,11,8,0
KEYOPT,11,9,0
KEYOPT,11,10,1
KEYOPT,11,11,0
KEYOPT,11,12,0
KEYOPT,11,2,0
KEYOPT,10,5,0
! Generate the target surface
ASEL,S,,,160
CM,_TARGET,AREA
TYPE,10
NSLA,S,1
ESLN,S,0
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,272
CM,_CONTACT,AREA
TYPE,11
NSLA,S,1
ESLN,S,0
ESURF
ALLS
ESEL,ALL
ESEL,S,TYPE,,10
ESEL,A,TYPE,,11
ESEL,R,REAL,,7
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,10
ESEL,A,TYPE,,11
ESEL,R,REAL,,7
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
/COM, CONTACT PAIR CREATION - END
/COM, CONTACT PAIR CREATION - START
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
MAT,1
MP,EMIS,1,7.88860905221e-031
R,8
REAL,8
ET,12,170
ET,13,174
R,8,,,1.0,0.1,0,-0.0001
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,13,4,0
KEYOPT,13,5,0
KEYOPT,13,7,0
KEYOPT,13,8,0
KEYOPT,13,9,0
KEYOPT,13,10,1
KEYOPT,13,11,0
KEYOPT,13,12,0
KEYOPT,13,2,0
KEYOPT,12,5,0
! Generate the target surface
ASEL,S,,,161
CM,_TARGET,AREA
TYPE,12
NSLA,S,1
ESLN,S,0
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,173
CM,_CONTACT,AREA
TYPE,13
NSLA,S,1
ESLN,S,0
ESURF
ALLS
ESEL,ALL
ESEL,S,TYPE,,12
ESEL,A,TYPE,,13
ESEL,R,REAL,,8
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,12
ESEL,A,TYPE,,13
ESEL,R,REAL,,8
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
/COM, CONTACT PAIR CREATION - END
FLST,5,20,5,ORDE,20
FITEM,5,24
FITEM,5,-25
FITEM,5,43
FITEM,5,52
FITEM,5,59
FITEM,5,67
FITEM,5,82
FITEM,5,85
FITEM,5,96
FITEM,5,103
FITEM,5,106
FITEM,5,123
FITEM,5,126
FITEM,5,136
FITEM,5,141
FITEM,5,148
FITEM,5,150
FITEM,5,153
FITEM,5,163
FITEM,5,172
ASEL,R, , ,P51X
CM,a1,AREA
FLST,5,23,6,ORDE,5
FITEM,5,36
FITEM,5,38
FITEM,5,40
FITEM,5,42
FITEM,5,-61
VSEL,S, , ,P51X
FLST,5,3,6,ORDE,3
FITEM,5,36
FITEM,5,38
FITEM,5,40
VSEL,U, , ,P51X
ASLV,S
FLST,5,20,5,ORDE,20
FITEM,5,181
FITEM,5,191
FITEM,5,197
FITEM,5,199
FITEM,5,208
FITEM,5,212
FITEM,5,215
FITEM,5,221
FITEM,5,228
FITEM,5,230
FITEM,5,237
FITEM,5,239
FITEM,5,247
FITEM,5,252
FITEM,5,-253
FITEM,5,258
FITEM,5,261
FITEM,5,264
FITEM,5,269
FITEM,5,273
ASEL,R, , ,P51X
CM,a2,AREA ALLS
CMSEL,S,A1
CMSEL,S,A2
/REPLO
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
MAT,1
MP,EMIS,1,7.88860905221e-031
R,9
REAL,9
ET,14,170
ET,15,174
R,9,,,1.0,0.1,0,-0.0001
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,15,4,0
KEYOPT,15,5,0
KEYOPT,15,7,0
KEYOPT,15,8,0
KEYOPT,15,9,0
KEYOPT,15,10,1
KEYOPT,15,11,0
KEYOPT,15,12,0
KEYOPT,15,2,0
KEYOPT,14,5,0
! Generate the target surface
ASEL,S,,,24
ASEL,A,,,25
ASEL,A,,,43
ASEL,A,,,52
ASEL,A,,,59
ASEL,A,,,67
ASEL,A,,,82
ASEL,A,,,85
ASEL,A,,,96
ASEL,A,,,103
ASEL,A,,,106
ASEL,A,,,123
ASEL,A,,,126
ASEL,A,,,136
ASEL,A,,,141
ASEL,A,,,148
ASEL,A,,,150
ASEL,A,,,153
ASEL,A,,,163
ASEL,A,,,172
CM,_TARGET,AREA
TYPE,14
NSLA,S,1
ESLN,S,0
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,181
ASEL,A,,,191
ASEL,A,,,197
ASEL,A,,,199
ASEL,A,,,208
ASEL,A,,,212
ASEL,A,,,215
ASEL,A,,,221
ASEL,A,,,228
ASEL,A,,,230
ASEL,A,,,237
ASEL,A,,,239
ASEL,A,,,247
ASEL,A,,,252
ASEL,A,,,253
ASEL,A,,,258
ASEL,A,,,261
ASEL,A,,,264
ASEL,A,,,269
ASEL,A,,,273
CM,_CONTACT,AREA
TYPE,15
NSLA,S,1
ESLN,S,0
ESURF
ALLS
ESEL,ALL
ESEL,S,TYPE,,14
ESEL,A,TYPE,,15
ESEL,R,REAL,,9
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,14
ESEL,A,TYPE,,15
ESEL,R,REAL,,9
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
ALLS
!接触对产生完毕
!产生预紧单元
PSMESH,1,pre1, ,VOLU,3, 0,Z, , , , , ,
PSMESH,2,pre2, ,VOLU,9, 0,Z, , , , , ,
PSMESH,3,pre3, ,VOLU,31, 0,Z, , , , , ,
FINI
!施加边界条件
!产生预紧力
/SOLU
*DEL,_z
*DIM,_z,ARRAY,1,1,1
*SET,_z(1) ,1
*DO,_z1,1,1
SLOAD,_z(_z1),9,LOCK,FORC, 5000,1,2
*ENDDO
*DEL,_z
*DIM,_z,ARRAY,1,1,1
*SET,_z(1) ,2
*DO,_z1,1,1
SLOAD,_z(_z1),9,LOCK,FORC, 5000,1,2
*ENDDO
*DEL,_z
*DIM,_z,ARRAY,1,1,1
*SET,_z(1) ,3
*DO,_z1,1,1
SLOAD,_z(_z1),9,LOCK,FORC, 5000,1,2
*ENDDO
!位移约束条件
FLST,2,30,4,ORDE,24
FITEM,2,217
FITEM,2,-220
FITEM,2,255
FITEM,2,-258
FITEM,2,295
FITEM,2,-298
FITEM,2,315
FITEM,2,324
FITEM,2,-325
FITEM,2,337
FITEM,2,-338
FITEM,2,348
FITEM,2,356
FITEM,2,-357
FITEM,2,372
FITEM,2,375
FITEM,2,384
FITEM,2,390
FITEM,2,403
FITEM,2,419
FITEM,2,-420
FITEM,2,428
FITEM,2,446
FITEM,2,457
DL,P51X, ,ALL,
FLST,5,5,6,ORDE,5
FITEM,5,5
FITEM,5,11
FITEM,5,15
FITEM,5,33
FITEM,5,35
VSEL,S, , ,P51X
ASLV,R
LSLA,R
FLST,2,5,4,ORDE,5
FITEM,2,2
FITEM,2,71
FITEM,2,117
FITEM,2,141
FITEM,2,168
DL,P51X, ,ALL,
ALLS
!求解控制
ANTYPE,0
NLGEOM,1
NSUBST,20,28,14
OUTRES,ERASE
OUTRES,BASI,1
LNSRCH,1
PRED,0
TIME,1
!写载荷步文件
LSWRITE,1,
FLST,2,10,5,ORDE,10
FITEM,2,37
FITEM,2,72
FITEM,2,76
FITEM,2,89
FITEM,2,92
FITEM,2,110
FITEM,2,113
FITEM,2,117
FITEM,2,130
FITEM,2,133
SFA,P51X,1,PRES,1e6
LSWRITE,2,
!求解
LSSOLVE,1,2
计算的结果如下：
（为了更清楚的反映变形的趋势，SCALE FACTOR 值设为50。）
整体应力分布
整体应变图
螺栓的应力分布
螺栓与支架及螺栓与立板间的接触应力（支架与立板间的接触应力没有将其显示，值很小）
请问：V***W,P51X,,DELETE !将V***W修改为V S B W(中间无空格）中的命令是什么意思啊，干什么用的啊，谢谢！运行时也运行不了。 要通过文件导入运行。。。。
页:
[1]