求教有关双椭球热源加载问题
对于薄板来说,一般加载高斯就行了,但是对于中厚板来说,加载高斯面热源,存在很大误差,加载双椭球热源是一个可行的办法,可是双椭球热源分为前后两部分,该怎么加呢? 希望能够讨论讨论,看过几个帖子都是关于高斯的!!!! 可以查看我以前的文章,有双拖求的加载的apdl网上有的高斯加载办法,两者实质是一样的,只有加载的时候公式不同,参数不一样,
如果没有弄懂高斯加载,那双拖求比较困难 学习一下,谢谢 /PREP7
!*
ET,1,SOLID70
!*
!*
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,600
MPTEMP,8,700
MPTEMP,9,800
MPTEMP,10,900
MPTEMP,11,1000
MPTEMP,12,1100
MPTEMP,13,1200
MPTEMP,14,1500
MPTEMP,15,1700
MPTEMP,16,2500
MPDATA,KXX,1,,50.7
MPDATA,KXX,1,,48.2
MPDATA,KXX,1,,46.5
MPDATA,KXX,1,,41.5
MPDATA,KXX,1,,38.1
MPDATA,KXX,1,,34.2
MPDATA,KXX,1,,30
MPDATA,KXX,1,,31.8
MPDATA,KXX,1,,26
MPDATA,KXX,1,,26.8
MPDATA,KXX,1,,26.8
MPDATA,KXX,1,,28.8
MPDATA,KXX,1,,29.2
MPDATA,KXX,1,,35
MPDATA,KXX,1,,14
MPDATA,KXX,1,,14.2
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,600
MPTEMP,8,700
MPTEMP,9,800
MPTEMP,10,900
MPTEMP,11,1000
MPTEMP,12,1500
MPTEMP,13,1700
MPTEMP,14,2500
MPDATA,C,1,,470
MPDATA,C,1,,490
MPDATA,C,1,,515
MPDATA,C,1,,550
MPDATA,C,1,,590
MPDATA,C,1,,645
MPDATA,C,1,,725
MPDATA,C,1,,860
MPDATA,C,1,,835
MPDATA,C,1,,680
MPDATA,C,1,,650
MPDATA,C,1,,660
MPDATA,C,1,,780
MPDATA,C,1,,820
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,7850
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,600
MPTEMP,8,700
MPTEMP,9,800
MPTEMP,10,1000
MPTEMP,11,1200
MPTEMP,12,1500
MPTEMP,13,2500
MPDATA,KXX,2,,11.9
MPDATA,KXX,2,,12.5
MPDATA,KXX,2,,14.3
MPDATA,KXX,2,,16.2
MPDATA,KXX,2,,17.3
MPDATA,KXX,2,,34.2
MPDATA,KXX,2,,32
MPDATA,KXX,2,,31.8
MPDATA,KXX,2,,26
MPDATA,KXX,2,,26.8
MPDATA,KXX,2,,28.8
MPDATA,KXX,2,,35
MPDATA,KXX,2,,14.2
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,2,,
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,500
MPTEMP,5,700
MPTEMP,6,1000
MPTEMP,7,1100
MPTEMP,8,1500
MPTEMP,9,1700
MPTEMP,10,2500
MPDE,C,2
MPDATA,C,2,,460
MPDATA,C,2,,490
MPDATA,C,2,,530
MPDATA,C,2,,570
MPDATA,C,2,,610
MPDATA,C,2,,645
MPDATA,C,2,,855
MPDATA,C,2,,670
MPDATA,C,2,,780
MPDATA,C,2,,820
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,2,,8350
BLOCK,0,-0.016,-0.016,0,0,0.3,
BLOCK,0,0.02,-0.016,0,0,0.3,
BLOCK,0.02,0.022,-0.016,0,0,0.3,
BLOCK,0.022,0.142,-0.016,0,0,0.3,
BLOCK,-0.016,0,0,0.02,0,0.3,
BLOCK,-0.016,0,0.02,0.022,0,0.3,
BLOCK,-0.016,0,0.022,0.15,0,0.3,
CYLIND,0.002,0,0,0.3,0,90,
FLST,2,8,6,ORDE,2
FITEM,2,1
FITEM,2,-8
VGLUE,P51X
FLST,5,7,6,ORDE,3
FITEM,5,1
FITEM,5,10
FITEM,5,-15
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 1, , 1, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
CM,_Y,VOLU
VSEL, , , , 9
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 2, , 1, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
/UI,MESH,OFF
ESIZE,0.002,0,
FLST,5,4,6,ORDE,4
FITEM,5,1
FITEM,5,9
FITEM,5,14
FITEM,5,-15
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VSWEEP,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
ESIZE,0.01,0,
FLST,5,2,6,ORDE,2
FITEM,5,11
FITEM,5,13
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VSWEEP,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
MSHKEY,0
MSHAPE,1,3d
FLST,5,2,6,ORDE,2
FITEM,5,10
FITEM,5,12
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
NSEL,S,EXT
FINISH
/SOL
FLST,2,11088,1,ORDE,18
FITEM,2,1
FITEM,2,-3303
FITEM,2,5539
FITEM,2,-8395
FITEM,2,9439
FITEM,2,-9587
FITEM,2,10631
FITEM,2,-13487
FITEM,2,39563
FITEM,2,-41066
FITEM,2,41096
FITEM,2,-41234
FITEM,2,41264
FITEM,2,-41406
FITEM,2,42074
FITEM,2,-42141
FITEM,2,42144
FITEM,2,-42211
/GO
!*
SF,P51X,CONV,200,8.9
*AFUN,DEG
*SET,ef,0.85
*SET,U,27
*SET,I,230
*set,pi,3.1415926
*SET,Q,ef*U*I
*SET,f1,0.6
*SET,f2,1.4
*SET,v,0.005
*SET,t1,60
*SET,step,0.002
*SET,tinc1,t1/150
*SET,i,1,150
*SET,XCENTER,0.141421356184E-02
*SET,YCENTER,0.141421356184E-02
*SET,ZCENTER,0
*SET,b11,0.01
*SET,b12,0.01
*SET,a1,0.01
*SET,c1,0.005
*SET,Qmax11,6*3**(1/2)*f1*q/(pi*pi**(1/2)*a1*b11*c1)
*SET,Qmax12,6*3**(1/2)*f2*q/(pi*pi**(1/2)*a1*b12*c1)
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNC_C5
*DEL,_FNC_C6
*DEL,_FNC_C7
*DEL,_FNCCSYS
*SET,_FNCNAME,'tfront1'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*DIM,_FNC_C4,,1
*DIM,_FNC_C5,,1
*DIM,_FNC_C6,,1
*DIM,_FNC_C7,,1
*SET,_FNC_C1(1),qmax11
*SET,_FNC_C2(1),xcenter
*SET,_FNC_C3(1),a1
*SET,_FNC_C4(1),ycenter
*SET,_FNC_C5(1),b11
*SET,_FNC_C6(1),zcenter
*SET,_FNC_C7(1),c1
*SET,_FNCCSYS,0
! /INPUT,tfront1.func,,,1
*DIM,%_FNCNAME%,TABLE,9,23,1,,,,%_FNCCSYS%
!
! Begin of equation: Qmax11*exp(-3((({x}-xcenter)/a1)^2+(({y}-ycenter)/b11)^2+
! (({z}-zcenter)/c1)^2))
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
*SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)%
*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%
*SET,%_FNCNAME%(8,0,1), %_FNC_C6(1)%
*SET,%_FNCNAME%(9,0,1), %_FNC_C7(1)%
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, 2, 2, 18
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 4, 19
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 2, 0, 0, -2
*SET,%_FNCNAME%(0,7,1), 0.0, -4, 0, 1, -2, 17, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 1, 3, 2, 20
*SET,%_FNCNAME%(0,9,1), 0.0, -2, 0, 1, -1, 4, 21
*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 2, 0, 0, -2
*SET,%_FNCNAME%(0,11,1), 0.0, -5, 0, 1, -2, 17, -1
*SET,%_FNCNAME%(0,12,1), 0.0, -1, 0, 1, -4, 1, -5
*SET,%_FNCNAME%(0,13,1), 0.0, -2, 0, 1, 4, 2, 22
*SET,%_FNCNAME%(0,14,1), 0.0, -4, 0, 1, -2, 4, 23
*SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 2, 0, 0, -4
*SET,%_FNCNAME%(0,16,1), 0.0, -5, 0, 1, -4, 17, -2
*SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 1, -1, 1, -5
*SET,%_FNCNAME%(0,18,1), 0.0, -1, 0, 3, 0, 0, -2
*SET,%_FNCNAME%(0,19,1), 0.0, -4, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 1, -3, 3, -4
*SET,%_FNCNAME%(0,21,1), 0.0, -1, 7, 1, -1, 0, 0
*SET,%_FNCNAME%(0,22,1), 0.0, -2, 0, 1, 17, 3, -1
*SET,%_FNCNAME%(0,23,1), 0.0, 99, 0, 1, -2, 0, 0
! End of equation: Qmax11*exp(-3((({x}-xcenter)/a1)^2+(({y}-ycenter)/b11)^2+
! (({z}-zcenter)/c1)^2))
!-->
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNC_C5
*DEL,_FNC_C6
*DEL,_FNC_C7
*DEL,_FNCCSYS
*SET,_FNCNAME,'tbehind1'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*DIM,_FNC_C4,,1
*DIM,_FNC_C5,,1
*DIM,_FNC_C6,,1
*DIM,_FNC_C7,,1
*SET,_FNC_C1(1),qmax12
*SET,_FNC_C2(1),xcenter
*SET,_FNC_C3(1),a1
*SET,_FNC_C4(1),ycenter
*SET,_FNC_C5(1),b12
*SET,_FNC_C6(1),zcenter
*SET,_FNC_C7(1),c1
*SET,_FNCCSYS,0
! /INPUT,tbehind1.func,,,1
*DIM,%_FNCNAME%,TABLE,9,23,1,,,,%_FNCCSYS%
!
! Begin of equation: Qmax12*exp(-3((({x}-xcenter)/a1)^2+(({y}-ycenter)/b12)^2+
! (({z}-zcenter)/c1)^2))
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
*SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)%
*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%
*SET,%_FNCNAME%(8,0,1), %_FNC_C6(1)%
*SET,%_FNCNAME%(9,0,1), %_FNC_C7(1)%
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, 2, 2, 18
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 4, 19
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 2, 0, 0, -2
*SET,%_FNCNAME%(0,7,1), 0.0, -4, 0, 1, -2, 17, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 1, 3, 2, 20
*SET,%_FNCNAME%(0,9,1), 0.0, -2, 0, 1, -1, 4, 21
*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 2, 0, 0, -2
*SET,%_FNCNAME%(0,11,1), 0.0, -5, 0, 1, -2, 17, -1
*SET,%_FNCNAME%(0,12,1), 0.0, -1, 0, 1, -4, 1, -5
*SET,%_FNCNAME%(0,13,1), 0.0, -2, 0, 1, 4, 2, 22
*SET,%_FNCNAME%(0,14,1), 0.0, -4, 0, 1, -2, 4, 23
*SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 2, 0, 0, -4
*SET,%_FNCNAME%(0,16,1), 0.0, -5, 0, 1, -4, 17, -2
*SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 1, -1, 1, -5
*SET,%_FNCNAME%(0,18,1), 0.0, -1, 0, 3, 0, 0, -2
*SET,%_FNCNAME%(0,19,1), 0.0, -4, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 1, -3, 3, -4
*SET,%_FNCNAME%(0,21,1), 0.0, -1, 7, 1, -1, 0, 0
*SET,%_FNCNAME%(0,22,1), 0.0, -2, 0, 1, 17, 3, -1
*SET,%_FNCNAME%(0,23,1), 0.0, 99, 0, 1, -2, 0, 0
! End of equation: Qmax12*exp(-3((({x}-xcenter)/a1)^2+(({y}-ycenter)/b12)^2+
! (({z}-zcenter)/c1)^2))
!-->
ANTYPE,4
TRNOPT,FULL
LUMPM,0
KBC,0
THOPT,FULL,0.05,64, , ,
EQSLV,PCG, ,1,
PRECISION,0
MSAVE,0
TOFFST,273, !
TUNIF,8.9,
*do,i,1,150,1 !
*SET,t,tinc1*i
time,t
DELTIM,TINC1,,,1 !
AUTOTS,on
outres,nsol,last, !
*SET,ZCENTER1,0.002*i
*SET,tfront1(4,0,1),0.141421356184E-02 !
*SET,tfront1(6,0,1),0.141421356184E-02
*SET,tfront1(8,0,1),ZCENTER1
LOCAL,11,2,0.141421356184E-02,0.141421356184E-02,ZCENTER1,,,,b11/a1,c1/a1,
ESEL,S,LIVE
nsle,r
NSEL,r,LOC,X,,a1
NSEL,R,LOC,Y,0,180
NSEL,R,LOC,Z,0,90
CM,NODEF,NODE
BF,NODEF,HGEN, %TFRONT1% !
CMDELE,NODEF
ALLSEL,ALL
*SET,tbehind1(4,0,1),0.141421356184E-02
*SET,tbehind1(6,0,1),0.141421356184E-02
*SET,tbehind1(8,0,1),ZCENTER1
LOCAL,12,2,0.141421356184E-02,0.141421356184E-02 ,ZCENTER1,,,,b12/a1,c1/a1,
ESEL,S,LIVE
nsle,r
NSEL,r,LOC,X,,a1
NSEL,R,LOC,Y,-180,0
NSEL,R,LOC,Z,0,90 !
CM,NODEB,NODE
BF,NODEB,HGEN, %TBEHIND1% !
CMDELE,NODEF
ALLSEL,ALL
SOLVE
BFDELE,all,HGEN
*ENDDO 这是我毕业设计时做的双椭球热源焊接问题,你看看吧 论坛里双椭球的例子还是很多的!
页:
[1]