24小时热门版块排行榜    

查看: 1156  |  回复: 2

gaochangqing

银虫 (初入文坛)

[求助] 复合材料烧结ansys仿真APDL命令流例子已有1人参与

求助复合材料烧结ansys仿真APDL命令流例子
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

姜呃呃

新虫 (小有名气)

引用回帖:
1楼: Originally posted by gaochangqing at 2017-11-07 10:43:34
求助复合材料烧结ansys仿真APDL命令流例子

可以交流下我也在做ansys仿真 用apdl

发自小木虫Android客户端
2楼2017-12-14 19:43:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

姜呃呃

新虫 (小有名气)

【答案】应助回帖

我在simwe里面找到的
finish
/clear
/units,si
/prep7

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

!===定义材料
!===实验粉末为316L不锈钢
mp,dens,1,7980 !密度(kg/m^3)

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,1430
mptemp,15,2000

mpdata,c,1,,477 !热导率(W/(m℃))
mpdata,c,1,,496
mpdata,c,1,,515
mpdata,c,1,,525
mpdata,c,1,,550
mpdata,c,1,,577
mpdata,c,1,,582
mpdata,c,1,,611
mpdata,c,1,,640
mpdata,c,1,,669
mpdata,c,1,,675
mpdata,c,1,,711
mpdata,c,1,,739
mpdata,c,1,,760
mpdata,c,1,,812

mpdata,kxx,1,,12.6 !比热容(J/(kg℃))
mpdata,kxx,1,,14.7
mpdata,kxx,1,,16.3
mpdata,kxx,1,,18.0
mpdata,kxx,1,,19.8
mpdata,kxx,1,,20.8
mpdata,kxx,1,,22.6
mpdata,kxx,1,,23.9
mpdata,kxx,1,,25.5
mpdata,kxx,1,,26.4
mpdata,kxx,1,,27.5
mpdata,kxx,1,,28.6
mpdata,kxx,1,,29.7
mpdata,kxx,1,,31.7
mpdata,kxx,1,,42


!===基材为304不锈钢
mp,dens,1,7930 !密度(kg/m^3)

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,1430
mptemp,15,2000

mpdata,c,2,,477 !热导率(W/(m℃))
mpdata,c,2,,496
mpdata,c,2,,515
mpdata,c,2,,525
mpdata,c,2,,550
mpdata,c,2,,577
mpdata,c,2,,582
mpdata,c,2,,611
mpdata,c,2,,640
mpdata,c,2,,669
mpdata,c,2,,675
mpdata,c,2,,711
mpdata,c,2,,739
mpdata,c,2,,760
mpdata,c,2,,812

mpdata,kxx,2,,12.6 !比热容(J/(kg℃))
mpdata,kxx,2,,14.7
mpdata,kxx,2,,16.3
mpdata,kxx,2,,18.0
mpdata,kxx,2,,19.8
mpdata,kxx,2,,20.8
mpdata,kxx,2,,22.6
mpdata,kxx,2,,23.9
mpdata,kxx,2,,25.5
mpdata,kxx,2,,26.4
mpdata,kxx,2,,27.5
mpdata,kxx,2,,28.6
mpdata,kxx,2,,29.7
mpdata,kxx,2,,31.7
mpdata,kxx,2,,42

!===熔覆层尺寸(m)
lx1=2.2/1000
ly1=30/1000
lz1=0.6/1000

!==基材尺寸(m)
lx2=10/1000
ly2=30/1000
lz2=8/1000

!===建模
block,-lx1/2,lx1/2,0,ly1,0,lz1
block,-lx2/2,lx2/2,0,ly2,0,-lz2/4
vglue,all
nummrg,all
numcmp,all

!===网格划分
type,1
mat,1
esize,0.0005
vmesh,1

type,1
mat,2
lesize,24,0.0005,,,5
lesize,26,0.0005,,,5
vsweep,2


!=====激光参数:功率P=630w、送粉量Mp=8g/min、速度V=6mm/s、离焦量ZI=-3mm
!qqq=(3*kk*P/(pi*r1**2))exp(-3*((Xc-x)**2+(Yc-y))**2/r1**2)

P=630 !激光功率(W)                          
r1=1.1/1000 !光斑半径(m)
kk=0.3 !激光功率利用率
pi=3.1415926
qq=3*kk*P/(pi*r1**2) !一个分步值

lsize=0.0005 !(m)
v=6/1000 !速度(m/s)
tinc=lsize/v !单步时间长度(s)
tmax=ly1/v !总时长(s)
StepNum=ly1/lsize

!=======后处理
/solu
ANTYPE,TRANS           !声明分析类型
trnopt,full            !瞬态分析完全法
nropt,full         !完全牛顿=拉普森选项
TIMINT,ON              !打开瞬态效应
TINTP,,,,1         !指定瞬态积分参数
AUTOTS,on           !使用自动时间步长
kbc,1 !载荷阶跃变化
tunif,20 !初始温度20度
tref,20

!===生死单元
esel,s,elem,,1,600,1  
ekill,all  
esel,s,live  
allsel,all
!===
*do,k,1,StepNum,1
tsres,erase
time,tinc*k
nsubst,60
xc=0
yc=lsize*(k-1) !熔覆中心点
*do,Enum,1,600,1
    dist1=SQRT((centrx(Enum)-xc)**2+(centry(Enum)-yc)**2)
     *if,dist1,le,r1,then
      ealive,Enum     
     *endif
   *enddo
   *do,Enum,1,600,1  
    dist1=SQRT((centrx(Enum)-xc)**2+(centry(Enum)-yc)**2)
    qqq=qq*exp(-3*dist1**2/r1**2)
     *if,dist1,le,r1,then
      sfe,Enum,6,hflux,,qqq
     *endif
   *enddo
esel,s,live  
eplot
allsel,all
outres,all,all
solve
sfedele,all,,hflux
*enddo
3楼2017-12-14 19:46:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gaochangqing 的主题更新
信息提示
请填处理意见