找回密码
 注册
Simdroid-非首页
查看: 292|回复: 6

求教有关双椭球热源加载问题

[复制链接]
发表于 2006-12-8 16:17:07 | 显示全部楼层 |阅读模式 来自 安徽合肥
对于薄板来说,一般加载高斯就行了,但是对于中厚板来说,加载高斯面热源,存在很大误差,加载双椭球热源是一个可行的办法,可是双椭球热源分为前后两部分,该怎么加呢?
 楼主| 发表于 2006-12-8 16:26:33 | 显示全部楼层 来自 安徽合肥
Simdroid开发平台
希望能够讨论讨论,看过几个帖子都是关于高斯的!!!!

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
发表于 2006-12-21 21:41:49 | 显示全部楼层 来自 广东深圳
可以查看我以前的文章,有双拖求的加载的apdl
网上有的高斯加载办法,两者实质是一样的,只有加载的时候公式不同,参数不一样,
如果没有弄懂高斯加载,那双拖求比较困难
发表于 2009-11-26 09:41:38 | 显示全部楼层 来自 湖北武汉
学习一下,谢谢
回复 不支持

使用道具 举报

发表于 2009-12-5 10:02:48 | 显示全部楼层 来自 黑龙江哈尔滨
/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
回复 不支持

使用道具 举报

发表于 2009-12-5 10:03:32 | 显示全部楼层 来自 黑龙江哈尔滨
这是我毕业设计时做的双椭球热源焊接问题,你看看吧
回复 不支持

使用道具 举报

发表于 2010-4-13 14:57:51 | 显示全部楼层 来自 吉林长春
论坛里双椭球的例子还是很多的!
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Simapps系列直播

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-9-20 19:38 , Processed in 0.059637 second(s), 13 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表