本帖最后由 Flamsky 于 2009-5-31 11:01 编辑
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 |