24小时热门版块排行榜    

查看: 1304  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-05 6/300 2026-11-05 23:30 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 6/300 2026-11-05 23:19 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 6/300 2026-11-05 23:08 by qL1cU46YIdK0
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ttH3f1oA6ace 2026-11-04 12/600 2026-11-05 22:57 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 13/650 2026-11-05 22:55 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 9/450 2026-11-05 22:55 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 14/700 2026-11-05 22:46 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 ncTVBM6AOJHF 2026-11-04 13/650 2026-11-05 22:44 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-05 22:33 by qL1cU46YIdK0
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +7 ncTVBM6AOJHF 2026-11-04 20/1000 2026-11-05 22:33 by qL1cU46YIdK0
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +7 ncTVBM6AOJHF 2026-11-04 17/850 2026-11-05 22:24 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 PYZmcHr5kSaq 2026-11-05 5/250 2026-11-05 21:27 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-05 20:43 by qL1cU46YIdK0
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-04 7/350 2026-11-05 20:43 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 8/400 2026-11-05 20:21 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ncTVBM6AOJHF 2026-11-04 9/450 2026-11-05 20:10 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ttH3f1oA6ace 2026-11-04 6/300 2026-11-05 18:00 by My20PsvIExYO
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 11:10 by PYZmcHr5kSaq
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-05 10:39 by PYZmcHr5kSaq
[基金申请] 各位基础老师们,看完中央《从严治党》的决定后,对学术界更加有信心了吧 +5 iaeyuan 2026-11-01 5/250 2026-11-04 14:43 by jinchengzi
信息提示
请填处理意见