24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者lk12将赠送您 5 个金币

lk12

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 90
帖子: 32
在线: 6.8小时
虫号: 24619712
注册: 2020-11-22
性别: GG
专业: 零件成形制造

[求助] ansys中使用表面单元加载对流载荷，求救如何使用

请教大佬表面效应单元的问题，我想在这个程序中增加一个表面效应的对流载荷，显示单元4未定义，esurf只能应用在线或者面上，请问怎么改呢？

FINISH
/CLEAR,START
/TITLE,Laser
/units,si
/VIEW,1,1,1,1
/prep7

!===定义单元
et,1,solid70

et,4,surf151

keyopt,4,4,1
keyopt,4,5,0
keyopt,4,8,4

!===定义材料
!===熔覆层为铝zl104
mp,dens,1,2630 !密度(kg/m^3)
mp,hf,1,100

mptemp,,,,,,,,
mptemp,1,25
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,1300

mpdata,c,1,,725 !比热容(J/(kg℃))
mpdata,c,1,,754
mpdata,c,1,,796
mpdata,c,1,,837
mpdata,c,1,,921
mpdata,c,1,,1011

mpdata,c,1,,1100
mpdata,c,1,,1130
mpdata,c,1,,1150
mpdata,c,1,,1170
mpdata,c,1,,1190
mpdata,c,1,,1210
mpdata,c,1,,1250

mpdata,kxx,1,,113 !热导率(W/(m℃))
mpdata,kxx,1,,144.9
mpdata,kxx,1,,149.1
mpdata,kxx,1,,149.1
mpdata,kxx,1,,144.9
mpdata,kxx,1,,145.1

mpdata,kxx,1,,145.1
mpdata,kxx,1,,145.1
mpdata,kxx,1,,145.1
mpdata,kxx,1,,145.1
mpdata,kxx,1,,145.1
mpdata,kxx,1,,145.1
mpdata,kxx,1,,145.1
mpdata,kxx,1,,145.1

!===熔覆层尺寸(m)
Width_Clad=2.5e-3 !熔覆层宽度
Length_Clad=30e-3
Height_Clad=0.35e-3

Layer=5

!==基材尺寸(m)
Width_Base=30e-3
Length_Base=50e-3
Height_Base=6e-3

!===建模
block,-Width_Clad/2,Width_Clad/2,0,Length_Clad,0,Height_Clad*Layer
block,-Width_Clad/2,Width_Clad/2,-10/1000,Length_Base-10/1000,0,-Height_Base/2
block,-Width_Clad/2,-Width_Base/2,-10/1000,Length_Base-10/1000,0,-Height_Base/2
block,Width_Clad/2,Width_Base/2,-10/1000,Length_Base-10/1000,0,-Height_Base/2
vglue,all
nummrg,all
numcmp,all

!===网格划分
type,1
mat,1
lesize,9,,,Layer
lesize,10,,,Layer
lesize,11,,,Layer
lesize,12,,,Layer
esize,2.2/1000/6
vmesh,1
type,1
mat,1
lesize,21,0.001,,,3
lesize,22,0.001,,,3
lesize,19,0.001,,,3
lesize,20,0.001,,,1/3
lesize,21,0.001,,,1/3
lesize,22,0.001,,,3
vsweep,4
type,1
mat,1
lesize,31,0.001,,,3
lesize,32,0.001,,,3
lesize,33,0.001,,,3
lesize,34,0.001,,,3
lesize,37,0.001,,,3
lesize,38,0.001,,,3
lesize,39,0.001,,,3
lesize,40,0.001,,,3
vsweep,2
vsweep,3

!=========激光参数
LaswerPower=1400 !激光功率(W)
Radius1=0.85/1000 !熔覆层半径(m)
Radius2=1.2/1000 !光斑半径(m)
Factor=0.35 !激光利用率
TruePower=LaswerPower*Factor

pi=3.14159265358979310
dt=0.00001 !小量
Init_Temp=20!初始温度
Coff_Conv=100 !对流换热系数
AirTemp=20
AirConvection=30
BaseTemp=20
BaseConvection=80

!===求解参数设置
StepDis=2.2/1000/6 !每个载荷步位移(m)
V=6/1000 !扫描速度(m/s)
StepTime=StepDis/V !单步时间长度(s)
TotalTime=(Length_Clad-2*Radius1)/V !总时长(s)
!StepNum=(Length_Clad-2*Radius1)/StepDis !载荷步数
StepNum=76 !载荷步数

!===对流
nsel,s,ext
nsel,u,loc,z,Height_Clad*Layer
SF,all,CONV,30,20
allsel,all

!===对流
nsel,s,ext
nsel,r,loc,z,Height_Clad*5

asel,s,loc,z,Height_Clad*5

esurf,all
ESEL,S,TYPE,,4

SFE,all,1,CONV,1,100
SFE,all,1,CONV,2,20

allsel,all

!===求解部分
/solu
antype,trans
nropt,full
tintpr,,,,1.0
trnopt,full !指定瞬态分析选项
timint,on !时间积分效应
autots,on !指定使用自动步长跟踪
deltim,,,,
nsubst,1
kbc,1
TUNIF,BaseTemp !初始温度20度
!toffst,273

!===杀死所有的熔覆层单元
VSEL,S,,,1 !选择熔覆层体
ESLV,S !选择依附于所选体上的单元
EKILL,ALL !杀死所选的单元(熔覆层)
ESEL,S,LIVE !选择活着的单元(基体)
EPLOT !刷新显示

LOCAL,11,0,0,0,0 !定义局部坐标系
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNCCSYS
*SET,_FNCNAME,'huanxing'
*SET,_FNCCSYS,11 !将方程坐标系设定为自己定义的局部坐标系
! /INPUT,.\function\huanxing.func,,,1
*DIM,%_FNCNAME%,TABLE,6,22,1,,,,%_FNCCSYS%
! Begin of equation: 1.1E8*EXP(-2*((sqrt({X}^2+{Y}^2)-0.00106)^2)/(0.00000025))
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*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, 2, 0, 0, -3
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 2, 0, 0, 2
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 2, 17, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,9,1), 0.0, -4, 0, 1, 3, 17, -1
*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 1, -3, 1, -4
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 16, 1, -1, 0, 0
*SET,%_FNCNAME%(0,12,1), 0.0, -3, 0, 0.00106, 0, 0, -1
*SET,%_FNCNAME%(0,13,1), 0.0, -4, 0, 1, -1, 2, -3
*SET,%_FNCNAME%(0,14,1), 0.0, -1, 0, 2, 0, 0, -4
*SET,%_FNCNAME%(0,15,1), 0.0, -3, 0, 1, -4, 17, -1
*SET,%_FNCNAME%(0,16,1), 0.0, -1, 0, 1, -2, 3, -3
*SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 0.00000025, 0, 0, -1
*SET,%_FNCNAME%(0,18,1), 0.0, -3, 0, 1, -1, 4, -2
*SET,%_FNCNAME%(0,19,1), 0.0, -1, 7, 1, -3, 0, 0
*SET,%_FNCNAME%(0,20,1), 0.0, -2, 0, 1.1E8, 0, 0, -1
*SET,%_FNCNAME%(0,21,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,22,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: 1.1E8*EXP(-2*((sqrt({X}^2+{Y}^2)-0.00106)^2)/(0.00000025))
!-->

!===第一层
time=0
XC=0
YC=3*StepDis

*DO,K,0,StepNum,1
LOCAL,11,0,XC,YC,0

t=0.45+StepTime*(k+1)
time,t

!===基于局部坐标系下复活第一层熔覆层
nsel,s,loc,x,-Width_Clad/2,Width_Clad/2
nsel,r,loc,z,0,Height_Clad
nsel,r,loc,y,-3*StepDis,3*StepDis
esln,s,1
ealive,all
allsel,all
ESEL,S,LIVE
EPLOT

!===第一层加载热源
nsel,s,loc,x,-Width_Clad/2,Width_Clad/2
nsel,r,loc,z,Height_Clad
nsel,r,loc,y,-3*StepDis,3*StepDis
SF,all,HFLUX,%huanxing%
ESEL,S,LIVE
/PSF,HFLUX,,2
/REPLOT
EPLOT
ALLSEL,ALL
OUTRES,ALL,ALL
SOLVE
YC=YC+StepDis
*ENDDO

!===第二层
CSYS,0
XC=0
YC=3*StepDis

*DO,k,0,StepNum,1
LOCAL,11,0,XC,YC,0

t=5.45+StepTime*(k+1)
time,t

!===基于局部坐标系下复活第二层熔覆层
nsel,s,loc,x,-Width_Clad/2,Width_Clad/2
nsel,r,loc,z,Height_Clad,Height_Clad*2
nsel,r,loc,y,-3*StepDis,3*StepDis
esln,s,1
ealive,all
allsel,all
ESEL,S,LIVE
EPLOT

!===第二层加载热源
nsel,s,loc,x,-Width_Clad/2,Width_Clad/2
nsel,r,loc,z,Height_Clad*2
nsel,r,loc,y,-3*StepDis,3*StepDis
SF,all,HFLUX,%huanxing%
ESEL,S,LIVE
/PSF,HFLUX,,2
/REPLOT
EPLOT
ALLSEL,ALL
OUTRES,ALL,ALL
SOLVE
YC=YC+StepDis
*ENDDO

!===第三层
CSYS,0
XC=0
YC=3*StepDis

*DO,k,0,StepNum,1
LOCAL,11,0,XC,YC,0

t=10.45+StepTime*(k+1)
time,t

!===基于局部坐标系下复活第三层熔覆层
nsel,s,loc,x,-Width_Clad/2,Width_Clad/2
nsel,r,loc,z,Height_Clad*2,Height_Clad*3
nsel,r,loc,y,-3*StepDis,3*StepDis
esln,s,1
ealive,all
allsel,all
ESEL,S,LIVE
EPLOT

!===第三层加载热源
nsel,s,loc,x,-Width_Clad/2,Width_Clad/2
nsel,r,loc,z,Height_Clad*3
nsel,r,loc,y,-3*StepDis,3*StepDis
SF,all,HFLUX,%huanxing%
ESEL,S,LIVE
/PSF,HFLUX,,2
/REPLOT
EPLOT
ALLSEL,ALL
OUTRES,ALL,ALL
SOLVE
YC=YC+StepDis
*ENDDO

!===第四层
CSYS,0
XC=0
YC=3*StepDis

*DO,k,0,StepNum,1
LOCAL,11,0,XC,YC,0

t=15.45+StepTime*(k+1)
time,t

!===基于局部坐标系下复活第四层熔覆层
nsel,s,loc,x,-Width_Clad/2,Width_Clad/2
nsel,r,loc,z,Height_Clad*3,Height_Clad*4
nsel,r,loc,y,-3*StepDis,3*StepDis
esln,s,1
ealive,all
allsel,all
ESEL,S,LIVE
EPLOT

!===第四层加载热源
nsel,s,loc,x,-Width_Clad/2,Width_Clad/2
nsel,r,loc,z,Height_Clad*4
nsel,r,loc,y,-3*StepDis,3*StepDis
SF,all,HFLUX,%huanxing%
ESEL,S,LIVE
/PSF,HFLUX,,2
/REPLOT
EPLOT
ALLSEL,ALL
OUTRES,ALL,ALL
SOLVE
YC=YC+StepDis
*ENDDO

!===第五层
CSYS,0
XC=0
YC=3*StepDis

*DO,k,0,StepNum,1
LOCAL,11,0,XC,YC,0

t=20.45+StepTime*(k+1)
time,t

!===基于局部坐标系下复活第五层熔覆层
nsel,s,loc,x,-Width_Clad/2,Width_Clad/2
nsel,r,loc,z,Height_Clad*4,Height_Clad*5
nsel,r,loc,y,-3*StepDis,3*StepDis
esln,s,1
ealive,all
allsel,all
ESEL,S,LIVE
EPLOT

!===第五层加载热源
nsel,s,loc,x,-Width_Clad/2,Width_Clad/2
nsel,r,loc,z,Height_Clad*5
nsel,r,loc,y,-3*StepDis,3*StepDis
SF,all,HFLUX,%huanxing%
ESEL,S,LIVE
/PSF,HFLUX,,2
/REPLOT
EPLOT
ALLSEL,ALL
OUTRES,ALL,ALL
SOLVE
YC=YC+StepDis
*ENDDO



FINISH
/POST1
PLNSOL,TEMP,,0

回复此楼