我做的是平板对接激光焊温度场模拟。加载过程中,先在长方体的五个面(SURF152表面效应单元)上加了对流,又在体(SOLID70单元)上加了生热率。但是结果只有一个未加对流载荷的面上加上了生热率,其余五个面温度是室温。百思不得其解,请高手指教。以下是我加的对流和生热率
ET,3,SURF152
KEYOPT,3,4,1
KEYOPT,3,5,1
KEYOPT,3,6,0
KEYOPT,3,8,1
KEYOPT,3,9,1
R,3,1,5.67E-14
TYPE,3
REAL,3
MAT,3
MPTEMP,,,,,,,,
MPTEMP,1,20,300,600,750,900,1200
MPTEMP,7,1500,1800,2000
MP,HF,MAT,1,6,50,120,180,200,250
MP,HF,MAT,7,378,700,850
ALLSEL
ASEL,S,LOC,X,0,0
ASEL,R,LOC,Z,0,PLATE_HEI
ASEL,R,LOC,Y,0,PLATE_WID
CM,AREA1,AREA
SFA,AREA1,,CONV,-3,20
ALLSEL
ASEL,S,LOC,X,PLATE_LEN,PLATE_LEN
ASEL,R,LOC,Z,0,PLATE_HEI
ASEL,R,LOC,Y,0,PLATE_WID
CM,AREA2,AREA
SFA,AREA2,,CONV,-3,20
ALLSEL
ASEL,S,LOC,Y,PLATE_WID,PLATE_WID
ASEL,R,LOC,Z,0,PLATE_HEI
ASEL,R,LOC,X,0,PLATE_LEN
CM,AREA3,AREA
SFA,AREA2,,CONV,-3,20
ALLSEL
ASEL,S,LOC,Z,0,0
ASEL,R,LOC,X,0,PLATE_LEN
ASEL,R,LOC,Y,0,PLATE_WID
CM,AREA4,AREA
SFA,AREA4,,CONV,-3,20
ALLSEL
ASEL,S,LOC,Z,PLATE_HEI,PLATE_HEI
ASEL,R,LOC,X,0,PLATE_LEN
ASEL,R,LOC,Y,0,PLATE_WID
CM,AREA5,AREA
SFA,AREA5,,CONV,-3,20
ALLSEL
SAVE
LO_COORD=0
RECY_NODE=1
LOAD_STEP=PLATE_LEN/WELD_GRID
*DO,CALCURECY,1,LOAD_STEP,1
LO_COORD=CALCURECY*WELD_GRID
LOCAL,20,1,LO_COORD,0,0,,,,1,1
ALLSEL,ALL
NSEL,ALL
NROTAT,ALL
ALLSEL,ALL
NSEL,S,LOC,Z,0,0
NSEL,R,LOC,X,0,1
ESLN,S
ESEL,R,TYPE,,2
CM,N01,ELEM
SAVE
ALLSEL,ALL
NSEL,S,LOC,Z,0.5,0.5
NSEL,R,LOC,X,0,0.7
ESLN,S
ESEL,R,TYPE,,2
CM,N02,ELEM
SAVE
ALLSEL,ALL
NSEL,S,LOC,Z,1,1
NSEL,R,LOC,X,0,0.4
ESLN,S
ESEL,R,TYPE,,2
CM,N03,ELEM
SAVE
ALLSEL,ALL
CMSEL,S,N01,ELEM
CMSEL,A,N02,ELEM
CMSEL,A,N03,ELEM
CM,N_%RECY_NODE%,ELEM
SAVE
CSDELE,20,,,
SAVE
CSYS,0
ALLSEL,ALL
NSEL,ALL
NROTAT,ALL
ALLSEL,ALL
RECY_NODE=RECY_NODE+1
*ENDDO
/SOLU
ANTYPE,TRANS,NEW
TRNOPT,FULL
TIMINT,ON
TINTP,,,,1
LUMPM,0
LNSRCH,1
NROPT,FULL,,
EQSLV,ICCG,1.0E-5
NEQIT,50
AUTOTS,1
KBC,1
OUTRES,BASIC,LAST
NSUBST,5,10,1
TUNIF,25
SAVE
ALLSEL
RECY_NODE=1
*DO,CALCURECY,1,LOAD_STEP,1
CSYS,0
ALLSEL,ALL
CMSEL,S,N_%RECY_NODE%,ELEM
BFE,ALL,HGEN,,QVOLUME
TIME,CALCURECY*WELDGDTM
ALLSEL,ALL
SOLVE
SAVE
ALLSEL,ALL
CMSEL,S,N_%RECY_NODE%,ELEM
BFEDELE,ALL,HGEN
ALLSEL
SAVE
RECY_NODE=RECY_NODE+1
*ENDDO |