24小时热门版块排行榜    

查看: 1295  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

gaochangqing

银虫 (初入文坛)

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

求助复合材料烧结ansys仿真APDL命令流例子
回复此楼
已阅   回复此楼   关注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的回帖
查看全部 3 个回答

姜呃呃

新虫 (小有名气)

引用回帖:
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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[职场人生] 40+中年人失业后,还能找到好工作吗? +10 xqy2009 2026-10-13 12/600 2026-10-15 19:24 by lpwang
[有机交流] 用NBS溴化时溶剂除了CCl4外还有什么 +3 18834002860 2026-10-12 6/300 2026-10-15 14:30 by 长沙摩尔信息科�
[基金申请] 申请科研项目纯靠实力能成功吗? 200+7 David1982 2026-10-13 11/550 2026-10-15 10:30 by zyqchem
[硕博家园] 极限拉扯,如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
信息提示
请填处理意见