各位虫友，帮忙看一下，我这个表面堆焊温度场的模拟有什么问题？我使用的是体生热率作为载荷加载的。代码如下：
FINISH
/CLEAR
/PREP7
/UNITS,SI

ET,1,SOLID70
PI=3.1415
LENGTH=0.03
I=160
U=22
V=0.007
Q=0.65*I*U

TIME_FREQUENCY=0.2
L_NUMBER=LENGTH/V/TIME_FREQUENCY

MPTEMP,1,20,100,200,300,500,700
MPTEMP,,900,1100,1300,1500,1700,1900
MPTEMP,,2000,2200,2400,2600

MPDATA,DENS,1,1,7800,7770,7750,7730,7790,7680
MPDATA,DENS,1,,7570,7460,7360,7030,6810,6620
MPDATA,DENS,1,,6540,6360,6170,6090

MPDATA,KXX,1,1,34.15,35.91,37.15,34.42,22.6,24.99
MPDATA,KXX,1,,27.32,29.84,32.17,33.45,36.72,40.19
MPDATA,KXX,1,,42.1,45.75,49.4,52.0

MPDATA,C,1,1,450,480,570,880,550,580
MPDATA,C,1,,610,640,680,12600,830,830
MPDATA,C,1,,830,830,830,820

MPDATA,ENTH,1,1,-6080,28910,83370,147720,364410,4.77E5
MPDATA,ENTH,1,,5.88E5,7.14E5,8.56E5,1.23E6,1.40E6,1.58E6
MPDATA,ENTH,1,,1.66E6,1.82E6,1.98E6,2.07E6

MPDATA,HF,1,1,1.336,6.68,13.36,20.4,38.7,71.15
MPDATA,HF,1,,110.04,170.6,180,264.4,270,280
MPDATA,HF,1,,290,310,340,360

BLOCK,0,0.06,0,LENGTH,0,0.015
BLOCK,0.06,0.08,0,LENGTH,0,0.015
BLOCK,0.08,0.14,0,LENGTH,0,0.015

BLOCK,0.066,0.074,0,LENGTH,0.015,0.018

VGLUE,ALL
AGLUE,ALL
LGLUE,ALL

NUMCMP,ALL

VSEL,S,,,2
VPLOT
ESIZE,0.002
TYPE,1
VSWEEP,ALL

VSEL,S,,,4,
VPLOT
ESIZE,0.002
TYPE,1
VSWEEP,ALL

VSEL,S,,,1,3,2
VPLOT
ESIZE,0.004
TYPE,1
VMESH,ALL

ALLSEL
*GET,ELEM_ACOUNT,ELEM,0,COUNT
VSEL,S,,,2
ESLV
*GET,AELEM,ELEM,0,COUNT
*DIM,ELEM_NUMBER,ARRAY,AELEM
*DIM,ELEM_CENTY,ARRAY,AELEM

*SET,A,1
*DO,K,1,ELEM_ACOUNT,1
        *IF,ESEL(K),EQ,1,THEN
                ELEM_NUMBER(A)=K
                *SET,A,A+1
        *ENDIF
*ENDDO

*DO,M,1,AELEM,1
        *GET,CENT_Y,ELEM,ELEM_NUMBER(M),CENT,Y
        ELEM_CENTY(M)=CENT_Y
*ENDDO

FINISH

/SOLU
ANTYPE,4
NROPT.FULL
NLGEOM,ON
OUTPR,ALL,ALL
OUTRES,ALL,ALL
TIMINT,OFF
KBC,1
AUTOTS,ON
NSUBST,4
TIME,0.01
ALLSEL
NSLV
D,ALL,TEMP,25
ALLSEL
SOLVE

DDELE,ALL,ALL
VSEL,S,,,2
ESLV
EKILL,ALL
TIME,0.02
ALLSEL
SOLVE

XC=0.070
ZC=0.015

TIME=0.02

QGEN=Q/0.008/0.002/0.007/TIME_FREQUENCY

*DO,I,1,L_NUMBER,1
        ALLSEL
        BFEDELE,ALL,ALL
        SFADELE,ALL,ALL
        YMIN=(I-1)*LENGTH/L_NUMBER
        YMAX=I*LENGTH/L_NUMBER
        *DO,N,1,AELEM,1
                *IF,ELEM_CENTY(N),LE,YMAX,THEN
                        *IF,ELEM_CENTY(N),GT,YMIN,THEN
                                BFE,ELEM_NUMBER(N),HGEN,1,QGEN
                                EALIVE,ELEM_NUMBER(N)
                        *ENDIF
                *ENDIF
        *ENDDO       

        ALLSEL
        ASEL,ALL
        SFA,ALL,,CONV,-1,25
        TIME=TIME+TIME_FREQUENCY
        TIMINT,ON
        TIME,TIME
        LNSRCH,ON
        NSUBST,5
        ALLSEL
        SOLVE
*ENDDO
计算出来的结果有点高


1.png 返回小木虫查看更多