我做的双椭球热源焊接过程

Flamsky

ansys-apdl高级工程应用实例分析上面的东东，不过书上有错误哦
另外，结果温度场分布与书上的也不太一样，书上的温度曲线峰值为960，我得出的有将近1800。不知道是哪里错了，麻烦大侠们帮我看下。
附apdl如下
FINI
/clear  
/filname,ex-welding,1
/title,the thermal analysis of the welding by Flamsky
/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,250
MPTEMP,3,400
MPTEMP,4,500
MPTEMP,5,750
MPTEMP,6,800
MPTEMP,7,1000   
MPTEMP,8,1200   
MPTEMP,9,1500   
MPTEMP,10,1700  
MPTEMP,11,2500  
MPDATA,HF,1,,100
MPDATA,HF,1,,350
MPDATA,HF,1,,480
MPDATA,HF,1,,520
MPDATA,HF,1,,1000   
MPDATA,HF,1,,1050   
MPDATA,HF,1,,1500   
MPDATA,HF,1,,2000   
MPDATA,HF,1,,3000   
MPDATA,HF,1,,3100   
MPDATA,HF,1,,3500   
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.0  
MPDATA,KXX,2,,26.8  
MPDATA,KXX,2,,28.8  
MPDATA,KXX,2,,35.0  
MPDATA,KXX,2,,14.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  
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,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,400
MPTEMP,5,500
MPTEMP,6,1000   
MPTEMP,7,1100   
MPTEMP,8,1500   
MPTEMP,9,1700   
MPTEMP,10,2500  
MPDATA,HF,2,,120
MPDATA,HF,2,,210
MPDATA,HF,2,,300
MPDATA,HF,2,,460
MPDATA,HF,2,,540
MPDATA,HF,2,,1200   
MPDATA,HF,2,,1200   
MPDATA,HF,2,,2800   
MPDATA,HF,2,,3000   
MPDATA,HF,2,,3100   
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,2,,8350
MPWRITE,'a','mat',' '   
/REPLOT,RESIZE  
*SET,length1,0.004  
*SET,length2,0.008  
*SET,length3,0.03   
*SET,width1,0.02
*SET,width2,0.04
*SET,height,0.003   
*SET,dia,0.001  
*SET,bh,0.005
*SET,w,500  
*SET,step,0.002
*SET,v,0.004
!
ah=dia/2   
chf=2*ah   
chb=8*ah   
rf=2*chf/(chf+chb)
rb=2*chb/(chf+chb)
count=width1/step  
!
PARSAV,SCALAR,'ex5-4','pram',' '
BLOCK,,length1,,width1,,height,
BLOCK,length1,length2,,width2,,height,  
BLOCK,length2,length3,,width2,,height,  
VGLUE,ALL
SAVE
TYPE,   1   
MAT,       1
ESIZE,HEIGHT/8,0,   
VSWEEP,1
MAT,2
ESIZE,HEIGHT/2
VSWEEP,4
MSHAPE,1,3D
MSHKEY,0
VMESH,5
/SOL
TUNIF,293,  
*DO,I,0,COUNT
DISY=I*STEP
T=(DISY+STEP)/V
/SOL
!*
ANTYPE,4
!*
TRNOPT,FULL
LUMPM,0
TIME,t  
NROPT,FULL,,ON
AUTOTS,1
NSUBST,5, , ,1  
KBC,1   
TSRES,ERASE
TIMINT,1,THERM
TINTP,0.005, , ,1,0.5,0.2,  
OUTRES,BASIC,LAST,  
LNSRCH,1
!*
NSEL,S,EXT  
NSEL,U,LOC,X,0  
SF,ALL,CONV,200,293
!*设定双椭球坐标系
ALLSEL,ALL
LOCAL,11,2,0,DISY,0, , , ,3,8,  
NSEL,S,LOC,X,,dia
NSEL,R,LOC,Y,,90  
NSEL,R,LOC,Z,,90
CM,NODEF,NODE
!*  
LOCAL,12,0,0,DISY,0, , , ,1,1,  
*DEL,TALBENAME
*DEL,TALBEMTID
*DEL,TABLE_ROW1
*DEL,TABLE_ROW2
*DEL,TABLE_ROW3
*DEL,TABLE_ROW4
*DEL,TABLE_ROW5
*DEL,TABLECSYS
*DEL,TEMPF
!
*SET,TABLENAME,'TEMPF'
*DIM,TABLE_ROW1,,1
*DIM,TABLE_ROW2,,1
*DIM,TABLE_ROW3,,1
*DIM,TABLE_ROW4,,1
*DIM,TABLE_ROW5,,1
*SET,TABLE_ROW1(1),RF
*SET,TABLE_ROW2(1),W
*SET,TABLE_ROW3(1),CHF
*SET,TABLE_ROW4(1),AH
*SET,TABLE_ROW5(1),BH  
*SET,TABLECSYS,12
*DIM,%TABLENAME%,TABLE,7,38,1,,,,%TABLECSYS%   
*SET,%TABLENAME%(0,0,1),0.0,-999
*SET,%TABLENAME%(2,0,1),0.0
*SET,%TABLENAME%(3,0,1),%TABLE_ROW1(1)%
*SET,%TABLENAME%(4,0,1),%TABLE_ROW2(1)%
*SET,%TABLENAME%(5,0,1),%TABLE_ROW3(1)%
*SET,%TABLENAME%(6,0,1),%TABLE_ROW4(1)%
*SET,%TABLENAME%(7,0,1),%TABLE_ROW5(1)%
*SET,%TABLENAME%(0,1,1),1.0,-1,16,3,0,0,0   
*SET,%TABLENAME%(0,2,1) , 0.0,-2,0,6,0,0,-1
*SET,%TABLENAME%(0,3,1) , 0.0,-3,0,1,-2,3,-1   
*SET,%TABLENAME%(0,4,1) , 0.0 ,-1,0,1,-3,3,17
*SET,%TABLENAME%(0,5,1) , 0.0   ,-2,0,1,-1,3,18
*SET,%TABLENAME%(0,6,1) , 0.0   ,-1,0,0,0,0,0
*SET,%TABLENAME%(0,7,1) , 0.0   ,-3,0,1,0,0,-1
*SET,%TABLENAME%(0,8,1) , 0.0   ,-4,0,1,-1,2,-3
*SET,%TABLENAME%(0,9,1) , 0.0   ,-1,0,2,0,0,2
*SET,%TABLENAME%(0,10,1) , 0.0  ,-3,0,1,2,17,-1
*SET,%TABLENAME%(0,11,1) , 0.0  ,-1,0,2,0,0,19
*SET,%TABLENAME%(0,12,1) , 0.0  ,-5,,1,19,17,-1
*SET,%TABLENAME%(0,13,1) , 0.0  ,-1,,1,-3,4,-5
*SET,%TABLENAME%(0,14,1) , 0.0  ,-3,0,2,0,0,3
*SET,%TABLENAME%(0,15,1) , 0.0  ,-5,0,1,3,17,-3
*SET,%TABLENAME%(0,16,1) , 0.0  ,-3,0,2,0,0,20
*SET,%TABLENAME%(0,17,1) , 0.0  ,-6,0,1,20,17,-3
*SET,%TABLENAME%(0,18,1) , 0.0  ,-3,0,1,-5,4,-6
*SET,%TABLENAME%(0,19,1) , 0.0  ,-5,0,1,-1,1,-3
*SET,%TABLENAME%(0,20,1) , 0.0  ,-1,0,2,0,0,4
*SET,%TABLENAME%(0,21,1) , 0.0  ,-3,0,1,4,17,-1
*SET,%TABLENAME%(0,22,1) , 0.0  ,-1,0,2,0,0,21
*SET,%TABLENAME%(0,23,1) , 0.0  ,-6,0,1,21,17,-1
*SET,%TABLENAME%(0,24,1) , 0.0  ,-1,0,1,-3,4,-6
*SET,%TABLENAME%(0,25,1) , 0.0  ,-3,0,1,-5,1,-1
*SET,%TABLENAME%(0,26,1) , 0.0  ,-1,0,3,0,0,-3
*SET,%TABLENAME%(0,27,1) , 0.0  ,-5,0,1,-1,3,-3
*SET,%TABLENAME%(0,28,1) , 0.0  ,-1,0,1,-4,3,-5
*SET,%TABLENAME%(0,29,1) , 0.0  ,-1,7,1,-1,0,0
*SET,%TABLENAME%(0,30,1) , 0.0  ,-3,0,1,-2,3,-1
*SET,%TABLENAME%(0,31,1) , 0.0  ,-1,0,1,20,3,21
*SET,%TABLENAME%(0,32,1) , 0.0  ,-2,0,1,-1,3,19
*SET,%TABLENAME%(0,33,1) , 0.0  ,-1,0,3.14,0,0,-2
*SET,%TABLENAME%(0,34,1) , 0.0  ,-4,0,1,-2,3,-1
*SET,%TABLENAME%(0,35,1) , 0.0  ,-1,16,3.14,0,0,0
*SET,%TABLENAME%(0,36,1) , 0.0  ,-2,0,1,-4,3,-1
*SET,%TABLENAME%(0,37,1) , 0.0  ,-1,0,1,-3,4,-2
*SET,%TABLENAME%(0,38,1) , 0.0  ,99,0,1,-1,0,0
BF,NODEF,HGEN,%TEMPF%
CMDELE,NODEF
ALLSEL,ALL
!*
LOCAL,13,2,0,DISY,0, , , ,3,8,  
NSEL,S,LOC,X,,DIA   
NSEL,R,LOC,Y,-90,0  
NSEL,R,LOC,Z,0,90
CM,NODEB,NODE   
!*  
LOCAL,14,0,0,DISY,0, , , ,1,1,   
!
*DEL,TABLENAME
*DEL,TABLEMTID
*DEL,TABLE_ROW1
*DEL,TABLE_ROW2
*DEL,TABLE_ROW3
*DEL,TABLE_ROW4
*DEL,TABLE_ROW5
*DEL,TABLECSYS
*DEL,TEMPB
!
*SET,TABLENAME,'TEMPB'
*DIM,TABLE_ROW1,,1
*DIM,TABLE_ROW2,,1
*DIM,TABLE_ROW3,,1
*DIM,TABLE_ROW4,,1
*DIM,TABLE_ROW5,,1
*SET,TABLE_ROW1(1),RB
*SET,TABLE_ROW2(1),W
*SET,TABLE_ROW3(1),CHB
*SET,TABLE_ROW4(1),AH
*SET,TABLE_ROW5(1),BH  
*SET,TABLECSYS,14
*DIM,%TABLENAME%,TABLE,7,38,1,,,,%TABLECSYS%   
*SET,%TABLENAME%(0,0,1),0.0,-999
*SET,%TABLENAME%(2,0,1),0.0
*SET,%TABLENAME%(3,0,1),%TABLE_ROW1(1)%
*SET,%TABLENAME%(4,0,1),%TABLE_ROW2(1)%
*SET,%TABLENAME%(5,0,1),%TABLE_ROW3(1)%
*SET,%TABLENAME%(6,0,1),%TABLE_ROW4(1)%
*SET,%TABLENAME%(7,0,1),%TABLE_ROW5(1)%
*SET,%TABLENAME%(0,1,1),1.0,-1,16,3,0,0,0   
*SET,%TABLENAME%(0,2,1) , 0,-2,0,6,0,0,-1
*SET,%TABLENAME%(0,3,1) , 0,-3,0,1,-2,3,-1   
*SET,%TABLENAME%(0,4,1) , 0   ,-1,0,1,-3,3,17
*SET,%TABLENAME%(0,5,1) , 0   ,-2,0,1,-1,3,18
*SET,%TABLENAME%(0,6,1) , 0   ,-1,0,0,0,0,0
*SET,%TABLENAME%(0,7,1) , 0   ,-3,0,1,0,0,-1
*SET,%TABLENAME%(0,8,1) , 0   ,-4,0,1,-1,2,-3
*SET,%TABLENAME%(0,9,1) , 0   ,-1,0,2,0,0,2
*SET,%TABLENAME%(0,10,1) , 0  ,-3,0,1,2,17,-1
*SET,%TABLENAME%(0,11,1) , 0  ,-1,0,2,0,0,19
*SET,%TABLENAME%(0,12,1) , 0  ,-5,,1,19,17,-1
*SET,%TABLENAME%(0,13,1) , 0  ,-1,,1,-3,4,-5
*SET,%TABLENAME%(0,14,1) , 0  ,-3,0,2,0,0,3
*SET,%TABLENAME%(0,15,1) , 0  ,-5,0,1,3,17,-3
*SET,%TABLENAME%(0,16,1) , 0  ,-3,0,2,0,0,20
*SET,%TABLENAME%(0,17,1) , 0  ,-6,0,1,20,17,-3
*SET,%TABLENAME%(0,18,1) , 0  ,-3,0,1,-5,4,-6
*SET,%TABLENAME%(0,19,1) , 0  ,-5,0,1,-1,1,-3
*SET,%TABLENAME%(0,20,1) , 0  ,-1,0,2,0,0,4
*SET,%TABLENAME%(0,21,1) , 0  ,-3,0,1,4,17,-1
*SET,%TABLENAME%(0,22,1) , 0  ,-1,0,2,0,0,21
*SET,%TABLENAME%(0,23,1) , 0  ,-6,0,1,21,17,-1
*SET,%TABLENAME%(0,24,1) , 0  ,-1,0,1,-3,4,-6
*SET,%TABLENAME%(0,25,1) , 0  ,-3,0,1,-5,1,-1
*SET,%TABLENAME%(0,26,1) , 0  ,-1,0,3,0,0,-3
*SET,%TABLENAME%(0,27,1) , 0  ,-5,0,1,-1,3,-3
*SET,%TABLENAME%(0,28,1) , 0  ,-1,0,1,-4,3,-5
*SET,%TABLENAME%(0,29,1) , 0  ,-1,7,1,-1,0,0
*SET,%TABLENAME%(0,30,1) , 0  ,-3,0,1,-2,3,-1
*SET,%TABLENAME%(0,31,1) , 0  ,-1,0,1,20,3,21
*SET,%TABLENAME%(0,32,1) , 0  ,-2,0,1,-1,3,19
*SET,%TABLENAME%(0,33,1) , 0  ,-1,0,3.14,0,0,-2
*SET,%TABLENAME%(0,34,1) , 0  ,-4,0,1,-2,3,-1
*SET,%TABLENAME%(0,35,1) , 0  ,-1,16,3.14,0,0,0
*SET,%TABLENAME%(0,36,1) , 0  ,-2,0,1,-4,3,-1
*SET,%TABLENAME%(0,37,1) , 0  ,-1,0,1,-3,4,-2
*SET,%TABLENAME%(0,38,1) , 0  ,99,0,1,-1,0,0
BF,NODEB,HGEN, %TEMPB%   
CMDELE,NODEB
ALLSEL,ALL
SOLVE
BFDELE,ALL,HGEN
*ENDDO
SAVE
TIME,20
NSUBST,5,,,1
SOLVE

Flamsky

没发过这种贴，图片顺序贴倒了

y1949b

恩，是热分析教材中的例子。有时间做做看

wangzhi82

看不懂:

xtfbnu

我被命令流震撼了，从今天起学习APDL编程

fangdian

太强劲了！

lilichunyan

为什么工件是两种材料属性啊，这模拟的什么情况呢？

Flamsky

7# lilichunyan
个人认为，应该只是增加一下例子的复杂程度而已。根据例子，两种材料的工件之间只是一般的导热吧，实际工程中这种情况可能比较难以见到。

acailaogu

不错，，有空搞搞，，，，，
温度场相对好弄，，，，

zhuwenfeng

我最近也被这本书上的这个例子困扰，
例子中的APDL程序中最关键的节点HGEN表格是如何形成这些数据的，“SET,%TABLENAME%(0,38,1) , 0  ,99,0,1,-1,0,0”，这些负数又如何理解，请高手说说一二吧

sd9366

没有人会吗？高手在那里？

王亮誉

本来我决定不会在社区回任何帖子了，但是看了你的帖子，我告诉自己这个帖子是一定要回的！这是百年难得一见的好贴啊！苍天有眼啊，让我在有生之年得以观得如此精彩绝伦的帖子！
楼主，你要继续努力啊！你是社区的希望啊！
我也许不能再看你的帖子了，这个回帖可能是我在这个社区回的最后的帖子，因为我可能会不能再上网了，但是我今天能回你的这个帖子，我已经感到非常满足了，不枉此生啊!!!


















折扣网  打折网   精品男装  时尚女装

tulongpower

请问是哪本书的啊，我需要参考下

tulongpower

请问是哪本书的啊，我需要参考下

tulongpower

请问是哪本书的啊，我需要参考下

sd9366

还没有人回答啊？

lixiaodong163

ansys-apdl高级工程应用实例分析与二次开发 14# tulongpower

lixiaodong163

用函数加载，然后再读取为表格形式，就是这样的格式了 10# zhuwenfeng

bqc

用函数加载，然后再读取为表格形式，就是这样的格式了 10# zhuwenfeng
lixiaodong163 发表于 2011-6-7 17:10

你好，这本书我有，例子我也做过。
我也明白是在命令流里的这些格式应该是界面操作生成的结果，
只是一直不明白在界面里应该这么操作呢，加载什么样的函数？
就是双椭球模型的两个公式吗？是公式生成的还是表格生成的啊？

还请赐教！谢谢

鼻涕虫害

热源模型和另外一个很像啊啊，可是都搞不明白~