ANSYS高手请进！急！！！！！
浏览：164468 回答：1
我做的是一个关于焊接温度场的模拟！一块板上堆焊一个筒体，！程序如下：
/CONFIG, NRES, 3700
/TITLE,Thermal/Stress analysis for SFF zig-zag wall
/PREP7
/UNITS, SI
/dev,font,1,clean,medium,r,14,,,,,
!------------------------------------------------------
!---------------------- Parameters
CW=0.001
!width of the bead （圆柱体的厚度）
CV=0.005
!scanning speed
CL=0.010
!
圆柱体的外径 length of the cladding - wall
LAYER=60
!the number of layerS
WH=0.0195
!the height of the wall
CH=WH/LAYER
!the layer hight
PW=182
!laser power [W] 9% of full power
EFF=0.375
!laser power efficiency
SSET=1600
!cooling time after deposition [s] or number of steps
WPARA=16
!parameter which defines the width of the substrate
SH=0.006
!height of the substrate
SW=CW*WPARA
!width of the substrate
SL=0.040
!length of the substrate
!---------------------- Model geometry
! Keypoints of the substrate
K,1,-SW/2,0,0
K,2,SW/2,0,0
K,3,SW/2,SL,0
K,4,-SW/2,SL,0
! Base surface of the substrate
A,1,2,3,4
! 圆柱体的基面
CYL4,0,SL/2,CL/2, ,CL/2-CW
AOVLAP,ALL
!overlaps the two areas
!---------------------Choose element type - 4 node brick element
ET,1,PLANE55
!---------------------- Material thermo-physical properties
! Material of the wall - H13 tool steel - 9 point model
MPTEMP,1,273,423,500,580,644,1600,1673,1730,1800,3500
MPDATA,KXX,1,1,17.3,22.7,23.7,24.5,25,29.5,29.7,29.8,29.8,31.2
MPDATA,DENS,1,1,7760,7650,7600,7580,7550,7200,7150,7000,6900,5900
MPDATA,C,1,1,460,486,1000,1000,538,757,1400,1400,800,800
!---------------------------------- Mesh the area
! Buildup
AMESH,2
ESEL,ALL
EREF,ALL, , ,1,0,1,1
! Substrate
LESIZE,1,CV/4
! element size in y direction defined by dividing scanning speed by 4
LESIZE,2,CV/4
LESIZE,3,CV/4
LESIZE,4,CV/4
MSHAPE,0,2d
MSHKEY,0
AMESH,4
AMESH,3
!--------------------------------- Iso view
EPLOT
/VIEW,1,-1,1,1
/VUP,1,Z
/AUTO,1
/REPLOT,FAST
!---------------------- Choose element - 8 node brick element ------
ET,2,SOLID70
!---------------------- Extrude the model ------------------------
! Substrate
TYPE,2
EXTOPT,ESIZE,9,3.5
!division in vertical directon
EXTOPT,ACLEAR,1
VEXT,3,,,0,0,-SH
!extrude volume bellow the wall
EXTOPT,ACLEAR,0
VEXT,2,,,0,0,-SH
!extrude the rest of the substrate
EXTOPT,ACLEAR,0
VEXT,4,,,0,0,-SH
!extrude the rest of the substrate
! Buildup - wall
TYPE,2
EXTOPT,ESIZE,LAYER
!division in vertical direction = # of layers
EXTOPT,ACLEAR,1
VEXT,2,,,0,0,WH
!extrude wall volume
!---------------------- Iso view ---------------------------------
EPLOT
/VIEW,1,-1,1,1
/VUP,1,Z
/REPLOT,FAST
!
!---------------------------- THERMAL ANALYSIS -------------
/SOLU
ANTYPE,TRANS
! transient analysis
NROPT,FULL,,OFF
! You must explicitly set the Newton-Raphson option
ESTIF,1E-20
! Sets non-default reduction factor (optional)
NEQIT,500
! max number of iterations
!-------- METAL DEPOSITION
!--------------------------- Disable the elements of the wall--------------
NSEL,S,LOC,Z,0,0.031 !select nodes between 0.1mm to 30 mm (high enough to select
! all elements) in Z direction
ESLN,S,1
!select elements with all their nodes in the set of selected nodes
EKILL,ALL
!disable all selected elements
ALLSEL,ALL
!---------------------------Initial condition----------------------------------
ESEL,S,LIVE
NSLE,S,ALL
TUNIF,298
!all initial temprature are 298 K
NSEL,INVE
D,ALL,TEMP,298
!constrains on all inactive dofs
ALLSEL,ALL
! Convection coefficients on the substrate
! top surface
VSEL,S,,,1
NSLV,S,1
NSEL,R,LOC,Z,0
SF,ALL,CONV,77,298
! side surfaces
ASEL,S,AREA,,19,22
NSLA,S,1
SF,ALL,CONV,77,298
!convection at bottom surface of the substrate
VSEL,S,,,1,2
ESLV,S
NSLE,S
NSEL,R,EXT
NSEL,R,LOC,Z,-SH
SF,ALL,CONV,20,298
ALLSEL,ALL
T_POOL=1850
!initial temperature of the molten pool - to calculate convection
!on molten pool surface
T_INC=0.2
!time increment
FLST,2,1,8
FITEM,2,0,0.2E-01,0
WPAVE,P51X
CSYS,1
NX=SL/2-CW/2
!initial x position
NY=0
YINC=7.5*T_INC/0.2
!y increment
NSTEP=360/(7.5*T_INC/0.2)
!number of steps per layer
T=0.2
!initial time
NZ=CH
!initial z position
HEAT1=PW*EFF/(CW*CH*(YINC))
!heat rate for one step (absorbed laser energy per !volume for one time step)
PIC=0
!set counter to zero
*DO,H,1,LAYER,1
*IF,H,GT,10,AND,H,LE,36,THEN
!increase convection on the molten pool for
!the next layer
T_POOL=T_POOL*1.0255
*ELSEIF,H,GT,36,AND,H,LE,45,THEN
T_POOL=T_POOL*1.015
*ELSEIF,H,GT,45
T_POOL=T_POOL*1.005
*ENDIF
*DO,L,1,NSTEP,1
! loop for deposition
PIC=PIC+1
! increase counter
TIME,T
! time step
! increase power in the first 10 layers to account for thin first several
! beads (high conduction through substrate)
! we generally use higher power at the beginning of the deposition to get
! uniform wall thickness
*IF,H,EQ,1,THEN
HEAT=HEAT1*1.8
!first layer
*ELSEIF,H,EQ,2
HEAT=HEAT1*1.44
!second layer
*ELSEIF,H,EQ,3
HEAT=HEAT1*1.3
!third layer
*ELSEIF,H,EQ,4
HEAT=HEAT1*1.2
!fourth layer
*ELSEIF,H,EQ,5
HEAT=HEAT1*1.15
*ELSEIF,H,EQ,6
HEAT=HEAT1*1.1
*ELSEIF,H,EQ,7
HEAT=HEAT1*1.08
*ELSEIF,H,EQ,8
HEAT=HEAT1*1.06
*ELSEIF,H,EQ,9
HEAT=HEAT1*1.05
*ELSEIF,H,EQ,10
HEAT=HEAT1*1.025
*ELSE
HEAT=HEAT1
*ENDIF
ALLSEL,ALL
!--- reactivate elements for the current step
NSEL,S,NODE,LOC,X,NX-CW/2,NX+CW/2
!select nodes between X range
NSEL,R,LOC,Y,NY,NY+YINC
!select nodes between Y range
NSEL,R,LOC,Z,NZ-CH,NZ
!select nodes between Z range
ESLN,S,1
!----select elements whose nodes are all in the range
! of selected nodes
EALIVE,ALL
!----reactivate the elements
CM,AG%PIC%,ELEM
!----create group of reactivated elements
NSLE,S,ALL
DDELE,ALL,TEMP
! delete constraints on reactivated elements
BFE,ALL,HGEN,,HEAT
! apply heat generation on elements
ALLSEL,ALL
SOLVE
OUTRES,ALL,ALL
! results output
! remove load
ALLSEL,ALL
BFEDELE,ALL,HGEN
NY=NY+YINC
! increase y coordinate
T=T+T_INC
! increase time
!GROUP COMPONENT - for animation purpose
!ESEL,S,LIVE
!CM,AG%PIC%,ELEM
*ENDDO
NZ=NZ+CH
! increase z coordinate - next layer
FLST,2,1,8
! 移动坐标系
FITEM,2,0.2E-01,0,NZ
WPAVE,P51X
*ENDDO
ALLSEL,ALL
!
!------- COOLING AFTER DEPOSITION
!
! convection on the substrate - variant with area selection
VSEL,S, , ,3
ASLV,S
SF,ALL,CONV,20,298
ALLSEL,ALL
T_INCC=1
! time increment 1 [s]
*DO,LC,1,SSET,1
TIME,T
! start time
!DELTIM,T_INC
! You must explicitly set the Newton-Raphson option
NROPT,FULL,,OFF
! Sets non-default reduction factor (optional)
ESTIF,1E-20
SOLVE
OUTRES,ALL,ALL
T=T+T_INCC
! increase time
*ENDDO
SAVE
FINISH
!
! ------- SAVE RESULTS FOR STRUCTURAL ANALYSIS
!
/POST26
ALLSEL,ALL
SSET1=SSET+LAYER*NSTEP
! number of steps deposition + cooling
FILE,Sff_test,RTH
! file name
*DEL,TIM
! delete TIM parameter if exists
*DIM,TIM,ARRAY,SSET1
! define TIM array
NSOL,2,1,TEMP
! read temperatures in nodes
VGET,TIM(1),1
! move data into array parameter vector
FINISH
提示选不到单元，我用到了柱坐标系！！！请高手帮忙解决一下！！！
模型和要的得到的结果在附件中！！！！
/CONFIG, NRES, 3700
/TITLE,Thermal/Stress analysis for SFF zig-zag wall
/PREP7
/UNITS, SI
/dev,font,1,clean,medium,r,14,,,,,
!------------------------------------------------------
!---------------------- Parameters
CW=0.001
!width of the bead （圆柱体的厚度）
CV=0.005
!scanning speed
CL=0.010
!
圆柱体的外径 length of the cladding - wall
LAYER=60
!the number of layerS
WH=0.0195
!the height of the wall
CH=WH/LAYER
!the layer hight
PW=182
!laser power [W] 9% of full power
EFF=0.375
!laser power efficiency
SSET=1600
!cooling time after deposition [s] or number of steps
WPARA=16
!parameter which defines the width of the substrate
SH=0.006
!height of the substrate
SW=CW*WPARA
!width of the substrate
SL=0.040
!length of the substrate
!---------------------- Model geometry
! Keypoints of the substrate
K,1,-SW/2,0,0
K,2,SW/2,0,0
K,3,SW/2,SL,0
K,4,-SW/2,SL,0
! Base surface of the substrate
A,1,2,3,4
! 圆柱体的基面
CYL4,0,SL/2,CL/2, ,CL/2-CW
AOVLAP,ALL
!overlaps the two areas
!---------------------Choose element type - 4 node brick element
ET,1,PLANE55
!---------------------- Material thermo-physical properties
! Material of the wall - H13 tool steel - 9 point model
MPTEMP,1,273,423,500,580,644,1600,1673,1730,1800,3500
MPDATA,KXX,1,1,17.3,22.7,23.7,24.5,25,29.5,29.7,29.8,29.8,31.2
MPDATA,DENS,1,1,7760,7650,7600,7580,7550,7200,7150,7000,6900,5900
MPDATA,C,1,1,460,486,1000,1000,538,757,1400,1400,800,800
!---------------------------------- Mesh the area
! Buildup
AMESH,2
ESEL,ALL
EREF,ALL, , ,1,0,1,1
! Substrate
LESIZE,1,CV/4
! element size in y direction defined by dividing scanning speed by 4
LESIZE,2,CV/4
LESIZE,3,CV/4
LESIZE,4,CV/4
MSHAPE,0,2d
MSHKEY,0
AMESH,4
AMESH,3
!--------------------------------- Iso view
EPLOT
/VIEW,1,-1,1,1
/VUP,1,Z
/AUTO,1
/REPLOT,FAST
!---------------------- Choose element - 8 node brick element ------
ET,2,SOLID70
!---------------------- Extrude the model ------------------------
! Substrate
TYPE,2
EXTOPT,ESIZE,9,3.5
!division in vertical directon
EXTOPT,ACLEAR,1
VEXT,3,,,0,0,-SH
!extrude volume bellow the wall
EXTOPT,ACLEAR,0
VEXT,2,,,0,0,-SH
!extrude the rest of the substrate
EXTOPT,ACLEAR,0
VEXT,4,,,0,0,-SH
!extrude the rest of the substrate
! Buildup - wall
TYPE,2
EXTOPT,ESIZE,LAYER
!division in vertical direction = # of layers
EXTOPT,ACLEAR,1
VEXT,2,,,0,0,WH
!extrude wall volume
!---------------------- Iso view ---------------------------------
EPLOT
/VIEW,1,-1,1,1
/VUP,1,Z
/REPLOT,FAST
!
!---------------------------- THERMAL ANALYSIS -------------
/SOLU
ANTYPE,TRANS
! transient analysis
NROPT,FULL,,OFF
! You must explicitly set the Newton-Raphson option
ESTIF,1E-20
! Sets non-default reduction factor (optional)
NEQIT,500
! max number of iterations
!-------- METAL DEPOSITION
!--------------------------- Disable the elements of the wall--------------
NSEL,S,LOC,Z,0,0.031 !select nodes between 0.1mm to 30 mm (high enough to select
! all elements) in Z direction
ESLN,S,1
!select elements with all their nodes in the set of selected nodes
EKILL,ALL
!disable all selected elements
ALLSEL,ALL
!---------------------------Initial condition----------------------------------
ESEL,S,LIVE
NSLE,S,ALL
TUNIF,298
!all initial temprature are 298 K
NSEL,INVE
D,ALL,TEMP,298
!constrains on all inactive dofs
ALLSEL,ALL
! Convection coefficients on the substrate
! top surface
VSEL,S,,,1
NSLV,S,1
NSEL,R,LOC,Z,0
SF,ALL,CONV,77,298
! side surfaces
ASEL,S,AREA,,19,22
NSLA,S,1
SF,ALL,CONV,77,298
!convection at bottom surface of the substrate
VSEL,S,,,1,2
ESLV,S
NSLE,S
NSEL,R,EXT
NSEL,R,LOC,Z,-SH
SF,ALL,CONV,20,298
ALLSEL,ALL
T_POOL=1850
!initial temperature of the molten pool - to calculate convection
!on molten pool surface
T_INC=0.2
!time increment
FLST,2,1,8
FITEM,2,0,0.2E-01,0
WPAVE,P51X
CSYS,1
NX=SL/2-CW/2
!initial x position
NY=0
YINC=7.5*T_INC/0.2
!y increment
NSTEP=360/(7.5*T_INC/0.2)
!number of steps per layer
T=0.2
!initial time
NZ=CH
!initial z position
HEAT1=PW*EFF/(CW*CH*(YINC))
!heat rate for one step (absorbed laser energy per !volume for one time step)
PIC=0
!set counter to zero
*DO,H,1,LAYER,1
*IF,H,GT,10,AND,H,LE,36,THEN
!increase convection on the molten pool for
!the next layer
T_POOL=T_POOL*1.0255
*ELSEIF,H,GT,36,AND,H,LE,45,THEN
T_POOL=T_POOL*1.015
*ELSEIF,H,GT,45
T_POOL=T_POOL*1.005
*ENDIF
*DO,L,1,NSTEP,1
! loop for deposition
PIC=PIC+1
! increase counter
TIME,T
! time step
! increase power in the first 10 layers to account for thin first several
! beads (high conduction through substrate)
! we generally use higher power at the beginning of the deposition to get
! uniform wall thickness
*IF,H,EQ,1,THEN
HEAT=HEAT1*1.8
!first layer
*ELSEIF,H,EQ,2
HEAT=HEAT1*1.44
!second layer
*ELSEIF,H,EQ,3
HEAT=HEAT1*1.3
!third layer
*ELSEIF,H,EQ,4
HEAT=HEAT1*1.2
!fourth layer
*ELSEIF,H,EQ,5
HEAT=HEAT1*1.15
*ELSEIF,H,EQ,6
HEAT=HEAT1*1.1
*ELSEIF,H,EQ,7
HEAT=HEAT1*1.08
*ELSEIF,H,EQ,8
HEAT=HEAT1*1.06
*ELSEIF,H,EQ,9
HEAT=HEAT1*1.05
*ELSEIF,H,EQ,10
HEAT=HEAT1*1.025
*ELSE
HEAT=HEAT1
*ENDIF
ALLSEL,ALL
!--- reactivate elements for the current step
NSEL,S,NODE,LOC,X,NX-CW/2,NX+CW/2
!select nodes between X range
NSEL,R,LOC,Y,NY,NY+YINC
!select nodes between Y range
NSEL,R,LOC,Z,NZ-CH,NZ
!select nodes between Z range
ESLN,S,1
!----select elements whose nodes are all in the range
! of selected nodes
EALIVE,ALL
!----reactivate the elements
CM,AG%PIC%,ELEM
!----create group of reactivated elements
NSLE,S,ALL
DDELE,ALL,TEMP
! delete constraints on reactivated elements
BFE,ALL,HGEN,,HEAT
! apply heat generation on elements
ALLSEL,ALL
SOLVE
OUTRES,ALL,ALL
! results output
! remove load
ALLSEL,ALL
BFEDELE,ALL,HGEN
NY=NY+YINC
! increase y coordinate
T=T+T_INC
! increase time
!GROUP COMPONENT - for animation purpose
!ESEL,S,LIVE
!CM,AG%PIC%,ELEM
*ENDDO
NZ=NZ+CH
! increase z coordinate - next layer
FLST,2,1,8
! 移动坐标系
FITEM,2,0.2E-01,0,NZ
WPAVE,P51X
*ENDDO
ALLSEL,ALL
!
!------- COOLING AFTER DEPOSITION
!
! convection on the substrate - variant with area selection
VSEL,S, , ,3
ASLV,S
SF,ALL,CONV,20,298
ALLSEL,ALL
T_INCC=1
! time increment 1 [s]
*DO,LC,1,SSET,1
TIME,T
! start time
!DELTIM,T_INC
! You must explicitly set the Newton-Raphson option
NROPT,FULL,,OFF
! Sets non-default reduction factor (optional)
ESTIF,1E-20
SOLVE
OUTRES,ALL,ALL
T=T+T_INCC
! increase time
*ENDDO
SAVE
FINISH
!
! ------- SAVE RESULTS FOR STRUCTURAL ANALYSIS
!
/POST26
ALLSEL,ALL
SSET1=SSET+LAYER*NSTEP
! number of steps deposition + cooling
FILE,Sff_test,RTH
! file name
*DEL,TIM
! delete TIM parameter if exists
*DIM,TIM,ARRAY,SSET1
! define TIM array
NSOL,2,1,TEMP
! read temperatures in nodes
VGET,TIM(1),1
! move data into array parameter vector
FINISH
提示选不到单元，我用到了柱坐标系！！！请高手帮忙解决一下！！！
模型和要的得到的结果在附件中！！！！
