焊接热源 选不中节点
最近在做焊接模拟,热源总是报错,请各位高手帮忙看下是哪儿的问题,小弟不胜感激。命令流如下,采用的是高斯热源,加载在面上的节点上/PREP7
K,1,0,500,1560,
K,2,750,500,1560,
K,3,750,-500,1560,
K,4,0,-500,1560,
K,5,0,200,1560,
K,6,0,-200,1560,
K,7,50,-200,1560,
K,8,50,200,1560,
LSTR, 1, 5
LSTR, 5, 8
LSTR, 8, 7
LSTR, 7, 6
LSTR, 6, 4
LSTR, 4, 3
LSTR, 3, 2
LSTR, 2, 1
FLST,2,8,4
FITEM,2,1
FITEM,2,8
FITEM,2,7
FITEM,2,6
FITEM,2,5
FITEM,2,3
FITEM,2,2
FITEM,2,4
AL,P51X
ET,1,PLANE82
FLST,5,3,4,ORDE,2
FITEM,5,2
FITEM,5,-4
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1,10, , , , , , ,1
FLST,5,5,4,ORDE,3
FITEM,5,1
FITEM,5,5
FITEM,5,-8
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1,50, , , , , , ,1
CM,_Y,AREA
ASEL, , , , 1
CM,_Y1,AREA
CMSEL,S,_Y
CMSEL,S,_Y1
AATT, , , 1, 0,
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
MSHAPE,0,2D
MSHKEY,0
CM,_Y,AREA
ASEL, , , , 1
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
AMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
ET,2,SOLID70
TYPE, 2
EXTOPT,ESIZE,3,0,
EXTOPT,ACLEAR,0
EXTOPT,ATTR,0,0,0
MAT,_Z2
REAL,_Z4
ESYS,0
FLST,2,1,5,ORDE,1
FITEM,2,1
VEXT,P51X, , ,0,0,-60,,,,
K,17,30,200,1480,
K,18,30,-200,1480,
K,19,0,-200,1480,
K,20,0,200,1480,
K,21,30,200,1430,
K,22,30,-200,1430,
K,23,0,-200,1430,
K,24,0,200,1430,
K,25,0,200,0,
K,26,0,-200,0,
K,27,30,-200,0,
K,28,30,200,0,
LSTR, 5, 6
LSTR, 20, 19
LSTR, 10, 20
LSTR, 19, 13
LSTR, 20, 17
LSTR, 17, 11
LSTR, 19, 18
LSTR, 18, 12
LSTR, 11, 12
LSTR, 20, 24
LSTR, 17, 21
LSTR, 24, 21
LSTR, 19, 23
LSTR, 18, 22
LSTR, 23, 22
LSTR, 21, 22
LSTR, 24, 23
LSTR, 24, 25
LSTR, 21, 28
LSTR, 23, 26
LSTR, 22, 27
LSTR, 25, 28
LSTR, 28, 27
LSTR, 25, 26
LSTR, 18, 17
LSTR, 26, 27
FLST,2,6,4
FITEM,2,25
FITEM,2,18
FITEM,2,27
FITEM,2,26
FITEM,2,28
FITEM,2,21
AL,P51X
FLST,2,4,4
FITEM,2,2
FITEM,2,25
FITEM,2,3
FITEM,2,4
AL,P51X
FLST,2,4,4
FITEM,2,29
FITEM,2,27
FITEM,2,10
FITEM,2,30
AL,P51X
FLST,2,4,4
FITEM,2,31
FITEM,2,28
FITEM,2,12
FITEM,2,32
AL,P51X
FLST,2,4,4
FITEM,2,29
FITEM,2,48
FITEM,2,26
FITEM,2,31
AL,P51X
FLST,2,4,4
FITEM,2,30
FITEM,2,11
FITEM,2,48
FITEM,2,32
AL,P51X
FLST,2,4,4
FITEM,2,29
FITEM,2,34
FITEM,2,35
FITEM,2,33
AL,P51X
FLST,2,4,4
FITEM,2,36
FITEM,2,38
FITEM,2,37
FITEM,2,31
AL,P51X
FLST,2,4,4
FITEM,2,35
FITEM,2,40
FITEM,2,39
FITEM,2,38
AL,P51X
FLST,2,4,4
FITEM,2,33
FITEM,2,40
FITEM,2,26
FITEM,2,36
AL,P51X
FLST,2,4,4
FITEM,2,34
FITEM,2,39
FITEM,2,48
FITEM,2,37
AL,P51X
FLST,2,4,4
FITEM,2,42
FITEM,2,41
FITEM,2,45
FITEM,2,35
AL,P51X
FLST,2,4,4
FITEM,2,43
FITEM,2,44
FITEM,2,49
FITEM,2,38
AL,P51X
FLST,2,4,4
FITEM,2,46
FITEM,2,47
FITEM,2,45
FITEM,2,49
AL,P51X
FLST,2,4,4
FITEM,2,41
FITEM,2,47
FITEM,2,43
FITEM,2,40
AL,P51X
FLST,2,4,4
FITEM,2,39
FITEM,2,44
FITEM,2,42
FITEM,2,46
AL,P51X
FLST,2,9,5,ORDE,4
FITEM,2,4
FITEM,2,-6
FITEM,2,11
FITEM,2,-16
VA,P51X
FLST,2,6,5,ORDE,3
FITEM,2,15
FITEM,2,17
FITEM,2,-21
VA,P51X
FLST,2,6,5,ORDE,3
FITEM,2,19
FITEM,2,22
FITEM,2,-26
VA,P51X
CM,_Y,VOLU
VSEL, , , , 2
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VSWEEP,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
FLST,5,4,4,ORDE,2
FITEM,5,41
FITEM,5,-44
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,30,3, , , ,1
FLST,5,2,4,ORDE,2
FITEM,5,45
FITEM,5,49
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , , , , , , ,1
FLST,5,2,4,ORDE,2
FITEM,5,45
FITEM,5,49
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1,30, , , , , , ,1
FLST,5,2,4,ORDE,2
FITEM,5,46
FITEM,5,-47
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1,50, , , , , , ,1
CM,_Y,VOLU
VSEL, , , , 4
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VSWEEP,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
MSHKEY,0
MSHAPE,1,3d
CM,_Y,VOLU
VSEL, , , , 3
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
NUMMRG,ALL, , , ,LOW
/PREP7
/units,si
et,,solid70
mptemp,1,20,250,500,750,1000,1500,1700,2500
mpDATA,ex,1,1,2.05e+11,1.87e+11,1.50e+11,7.00e+10,2.00e+10,1.90e+10,1.80e+10,1.20e+10
mpdata,kxx,1,1,50,47,40,27,30,35,140,142
mpdata,c,1,1,460,480,530,675,670,660,780,820
mpdata,dens,1,1,7820,7700,7610,7550,7490,7350,7300,7090
/SOL
*SET,r,0.045
*SET,remain,r/10
*SET,pi,3.1415926
*SET,v,0.01
*SET,i,180
*SET,u,25
*SET,y1,-0.2
*SET,t1,0.4/v
*SET,tinc1,0.01/v
*SET,ef,0.75
*SET,Q,ef*i*u
*SET,QM,3*Q/(PI*R**2)
*SET,q1,qm*(1-0.82902)*(pi*r**2)/3
*SET,q2,qm*(0.82902-0.47236)*(pi*r**2)/3
*SET,q3,qm*(0.47236-0.18498)*(pi*r**2)/3
*SET,q4,qm*(0.18498-0.04978)*(pi*r**2)/3
*SET,hgen1,q1/(0.002*pi*(r/4)**2)
*SET,hgen2,q2/(0.002*pi*((r/2)**2-(r/4)**2))
*SET,hgen3,q3/(0.002*pi*((3*r/4)**2-(r/2)**2))
*SET,hgen4,q4/(0.002*pi*(r**2-(3*r/4)**2))
ANTYPE,TRANS
TRNOPT,FULL
NROPT,AUTO,,
EQSLV,
SOLCONTROL,ON
AUTOTS,ON
KBC,1
OUTRES,ALL,5
TIME,0.001
DELTIM,0.001,0.001,0.001
TIMINT,OFF
TUNIF,25
SFA,1,,CONV,10,25
SFA,2,,CONV,10,25
SFA,3,,CONV,10,25
SFA,7,,CONV,10,25
SFA,8,,CONV,10,25
SFA,9,,CONV,10,25
SFA,10,,CONV,10,25
SFA,11,,CONV,10,25
SFA,12,,CONV,10,25
SFA,13,,CONV,10,25
SFA,14,,CONV,10,25
SFA,16,,CONV,10,25
SFA,17,,CONV,10,25
SFA,18,,CONV,10,25
SFA,20,,CONV,10,25
SFA,21,,CONV,10,25
SFA,22,,CONV,10,25
SFA,23,,CONV,10,25
SFA,24,,CONV,10,25
SFA,25,,CONV,10,25
SFA,26,,CONV,10,25
ALLSEL,ALL
SOLVE
TIMINT,ON
*DO,T,TINC1,T1-1,TINC1
TIME,T
*SET,Y1,Y1+0.01
AUTOTS,ON
AUTOTS,0
DELTIM,TINC1/2,TINC1/2,TINC1/2
KBC,1
OUTRES,ALL,5
TSRES,ERASE
WPSTYL,,,,,,,,1
CSYS,4
NWPAVE,2225
CSYS,4
WPOFF,,,Y1
CSWPLA,11,1,1,1,
ASEL,S,,,16
NSLA,R,1
NSEL,R,LOC,Y,0,R/4+REMAIN
NPLOT
BF,ALL,HGEN,HGEN1
ALLSEL,ALL
ASEL,S,,,16
NSLA,R,1
NSEL,R,LOC,Y,R/4-REMAIN,R/2+REMAIN
NPLOT
BF,ALL,HGEN,HGEN2
ALLSEL,ALL
ASEL,S,,,16
NSLA,R,1
NSEL,R,LOC,Y,R/2-REMAIN,3*R/4+REMAIN
NPLOT
BF,ALL,HGEN,HGEN3
ALLSEL,ALL
ASEL,S,,,16
NSLA,R,1
NSEL,R,LOC,Y,3*R/4-REMAIN,R+REMAIN
NPLOT
BF,ALL,HGEN,HGEN4
ALLSEL,ALL
SOLVE
TIMINT,ON
ASEL,S,,,16
NSLA,R,1
NSEL,R,LOC,Y,0,R/4+REMAIN
NPLOT
BF,ALL,HGEN,1E-10
ALLSEL,ALL
ASEL,S,,,16
NSLA,R,1
NSEL,R,LOC,Y,R/4-REMAIN,R/2+REMAIN
NPLOT
BF,ALL,HGEN,1E-10
ALLSEL,ALL
ASEL,S,,,16
NSLA,R,1
NSEL,R,LOC,Y,R/2-REMAIN,3*R/4+REMAIN
NPLOT
BF,ALL,HGEN,1E-10
ALLSEL,ALL
ASEL,S,,,16
NSLA,R,1
NSEL,R,LOC,Y,3*R/4-REMAIN,R+REMAIN
NPLOT
BF,ALL,HGEN,1E-10
ALLSEL,ALL
solve
*ENDDO
time,7200
nsubst,40
bfdele,all,hgen
solve
outres,all,5
FINISH 各位帮忙看看啊,自己检查了很多遍都没有查到错在哪了 我也不会,呵呵 SET,QM,3*Q/(PI*R**2)
*SET,q1,qm*(1-0.82902)*(pi*r**2)/3
*SET,q2,qm*(0.82902-0.47236)*(pi*r**2)/3
*SET,q3,qm*(0.47236-0.18498)*(pi*r**2)/3
*SET,q4,qm*(0.18498-0.04978)*(pi*r**2)/3
*SET,hgen1,q1/(0.002*pi*(r/4)**2)
*SET,hgen2,q2/(0.002*pi*((r/2)**2-(r/4)**2))
*SET,hgen3,q3/(0.002*pi*((3*r/4)**2-(r/2)**2))
*SET,hgen4,q4/(0.002*pi*(r**2-(3*r/4)**2))
楼主是不是在做多热源啊?? 4# tulongpower 不是多热源,热源要移动,要把先前一步加进的消除,设的那些是后便消除用的 自己发现了模型的参数有点问题,考完试改过来之后还是选不中节点 SET,QM,3*Q/(PI*R**2)
*SET,q1,qm*(1-0.82902)*(pi*r**2)/3
*SET,q2,qm*(0.82902-0.47236)*(pi*r**2)/3
*SET,q3,qm*(0.47236-0.18498)*(pi*r**2)/3
*SET,q4,qm*(0.18498-0.04978)*(pi*r**2)/3
*SET,hgen1,q1/(0.002*pi*(r/4)**2)
*SET,hgen2,q2/(0.002*pi*((r/2)**2-(r/4)**2))
*SET,hgen3,q3/(0.002*pi*((3*r/4)**2-(r/2)**2))
*SET,hgen4,q4/(0.002*pi*(r**2-(3*r/4)**2))
不是多热源,热源要移动,要把先前一步加进的消除,设的那些是后便消除用的 你这是热源离散呀 回答不对吧 要把热源运动的距离和时间建立起关系 楼主的问题解决了吗?
页:
[1]