版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

姜呃呃

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 80.4
散金: 10
帖子: 84
在线: 13.7小时
虫号: 6319137
注册: 2017-04-14
性别: GG
专业: 制造系统与自动化

[求助] 仿真慢已有1人参与

刚才投错区了，

ANSYS运算非常慢，请求帮助，程序如下方

回复此楼

» 猜你喜欢

拟解决的关键科学问题还要不要写已经有7人回复
请教限项目规定已经有3人回复
存款400万可以在学校里躺平吗已经有15人回复
Materials Today Chemistry审稿周期已经有6人回复
基金委咋了？2026年的指南还没有出来？已经有10人回复
基金申报已经有6人回复
推荐一本书已经有13人回复
国自然申请面上模板最新2026版出了吗？已经有17人回复
纳米粒子粒径的测量已经有8人回复
疑惑？已经有5人回复

1楼 2018-03-23 15:34:22

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

姜呃呃

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 80.4
散金: 10
帖子: 84
在线: 13.7小时
虫号: 6319137
注册: 2017-04-14
性别: GG
专业: 制造系统与自动化

mp,reft,1,10
mp,reft,2,10
mp,reft,3,10
v=0.03       !熔覆线速度
lx_number=125 !x方向分段数
ly_number=75  !y方向分段数
tinc=0.13    !单步时间长度
Q=2300       !激光功率
yita=0.03    !激光效率
R=0.002    !激光有效加热半径
pi=3.1415926
qq=Q*yita/(pi*R**2)  !一个分步值

/SOLU
ANTYPE,4
TRNOPT,FULL
LUMPM,0
OUTRES,ERASE
OUTRES,BASI,LAST
EQSLV,PCG,1E-6
KBC,0.0
TIMINT,0
TRNOPT,FULL
THOPT,FULL,0.05,64, , , ,0
EQSLV,PCG, ,2, ,DELE
MSAVE,0
TUNIF,10,
TREF,10,
OUTPR,BASIC,LAST,
OUTRES,ALL,LAST,
SOLCONTROL,OFF, ,NOPL
TIME,0.13
AUTOTS,0
NSUBST,1,0,0,0
KBC,0
TSRES,ERASE
TIMINT,1
TINTP,0.005, , ,1,0.5,0.2,
SOLVE

alls
esel,s,mat,,3
esel,a,mat,,2
estif,1.0e-10
ekill,all
alls
timint,on
asel,s,,,1,6,1
nsla,s,1
sf,all,conv,-2,10
alls
!toffst,273
tintp,0.005,,,1.0,0.5,0.2
t=0.26
time,t
!nsubst,2
outpr,basic,all
outres,basic,all
alls
solve

alls
!autots,on
!*do,i,1,lx_number
*do,i,1,2
!*do,j,1,ly_number
*do,j,1,4
sfedele,all,6,hflux
tsres,erase
xc=0.002+(i-1)*0.008
yc=0.002+(j-1)*0.004
zc=0.043
*do,enum,46876,75000
dist1=SQRT((centrx(enum)-xc)**2+(centry(enum)-yc)**2+(centrz(enum)-zc)**2)
*if,dist1,le,R,then
ealive,enum
eplot
qqq=qq
sfe,enum,6,hflux,,qqq
*endif
*enddo
alls
*do,enumm,1,46875
dist1=SQRT((centrx(enumm)-xc)**2+(centry(enumm)-yc)**2+(centrz(enumm)-0.042)**2)
*if,dist1,le,R,then
esel,s,,,enumm
sfdele,all,conv
*endif
*enddo
alls
*do,enumn,103126,112500
dist1=SQRT((centrx(enumn)-xc)**2+(centry(enumn)-yc)**2+(centrz(enumn)-zc)**2)
*if,dist1,le,R,then
ealive,enumn
eplot
esel,s,,,enumn
nsle
sf,all,conv,-2,10
*endif
*enddo
alls
esel,s,live
eplot
allsel,all
lnsrch,on
!nsubst,2
t=t+tinc
time,t
!*if,mod(j,5),ne,0,then
!outres,,none
!outpr,,none
!*else
outres,basic,last
outpr,basic,last
!*endif
solve
*enddo
alls

!*do,p,1,ly_number
*do,p,1,4
sfedele,all,6,hflux
tsres,erase
xc=0.006+(i-1)*0.008
!yc=0.298-0.004*(p-1)
yc=0.014-0.004*(p-1)
zc=0.043
*do,enum,46876,75000
dist1=SQRT((centrx(enum)-xc)**2+(centry(enum)-yc)**2+(centrz(enum)-zc)**2)
*if,dist1,le,R,then
ealive,enum
eplot
qqq=qq
sfe,enum,6,hflux,,qqq
*endif
*enddo
*do,enumm,1,46875
dist1=SQRT((centrx(enumm)-xc)**2+(centry(enumm)-yc)**2+(centrz(enumm)-0.042)**2)
*if,dist1,le,R,then
esel,s,,,enumm
sfdele,all,conv
*endif
*enddo
alls
*do,enumn,103126,112500
dist1=SQRT((centrx(enumn)-xc)**2+(centry(enumn)-yc)**2+(centrz(enumn)-zc)**2)
*if,dist1,le,R,then
ealive,enumn
eplot
esel,s,,,enumn
nsle
sf,all,conv,-2,10
*endif
*enddo
alls
esel,s,live
eplot
allsel,all
outres,basic,last
outpr,basic,last
lnsrch,on
!nsubst,2
t=t+tinc
time,t
!*if,mod(p,5),ne,0,then
!outres,,none
!outpr,,none
!*else
outres,basic,last
outpr,basic,last
!*endif
solve
*enddo
*enddo
alls
outpr,basic,all
outres,basic,all
sfedele,all,6,hflux
t=t+300
time,t
solve

!第二层
alls
!*do,j,1,ly_number
*do,j,1,2
!*do,i,1,lx_number
*do,i,1,4
sfedele,all,6,hflux
tsres,erase
xc=0.002+(i-1)*0.004
yc=0.002+(j-1)*0.008
zc=0.044
*do,enum,75001,103125
dist1=SQRT((centrx(enum)-xc)**2+(centry(enum)-yc)**2+(centrz(enum)-zc)**2)
*if,dist1,le,R,then
ealive,enum
eplot
qqq=qq
sfe,enum,6,hflux,,qqq
*endif
*enddo
alls
*do,enumm,46876,75000
dist1=SQRT((centrx(enumm)-xc)**2+(centry(enumm)-yc)**2+(centrz(enumm)-0.043)**2)
*if,dist1,le,R,then
esel,s,,,enumm
sfdele,all,conv
*endif
*enddo
alls
*do,enumn,112501,121875
dist1=SQRT((centrx(enumn)-xc)**2+(centry(enumn)-yc)**2+(centrz(enumn)-zc)**2)
*if,dist1,le,R,then
ealive,enumn
eplot
esel,s,,,enumn
nsle
sf,all,conv,-2,10
*endif
*enddo
alls
esel,s,live
eplot
allsel,all
lnsrch,on
nsubst,1
t=t+tinc
time,t
!*if,mod(i,5),ne,0,then
!outres,,none
!outpr,,none
!*else
outres,basic,last
outpr,basic,last
!*endif
solve
*enddo
alls

!*do,k,1,lx_number
*do,k,1,4
sfedele,all,6,hflux
tsres,erase
!xc=0.48-0.004*(k-1)
xc=0.014-0.004*(k-1)
yc=0.006+(j-1)*0.008
zc=0.044
*do,enum,75001,103125
dist1=SQRT((centrx(enum)-xc)**2+(centry(enum)-yc)**2+(centrz(enum)-zc)**2)
*if,dist1,le,R,then
ealive,enum
eplot
qqq=qq
sfe,enum,6,hflux,,qqq
*endif
*enddo
*do,enumm,46876,75000
dist1=SQRT((centrx(enumm)-xc)**2+(centry(enumm)-yc)**2+(centrz(enumm)-0.043)**2)
*if,dist1,le,R,then
esel,s,,,enumm
sfdele,all,conv
*endif
*enddo
*do,enumn,112501,121875
dist1=SQRT((centrx(enumn)-xc)**2+(centry(enumn)-yc)**2+(centrz(enumn)-zc)**2)
*if,dist1,le,R,then
ealive,enumn
eplot
esel,s,,,enumn
nsle
sf,all,conv,-2,10
*endif
*enddo
alls
esel,s,live
eplot
allsel,all
outres,basic,last
outpr,basic,last
lnsrch,on
nsubst,1
time=time+tinc
time,time
!*if,mod(k,5),ne,0,then
! outres,,none
! outpr,,none
!*else
outres,basic,last
outpr,basic,last
!*endif
solve
*enddo
*enddo
alls
t=t+tinc
time,t
sfedele,all,6,hflux
solve
alls
outpr,basic,all
outres,basic,all
t=t+3600
time,t
solve
!save
!finish

/prep7
etchg,tts
mptemp,1,0,20,200,600,1200,2000
mpdata,ex,1,1,2.060e11,2.056e11,2.019e11,1.936e11,1.183e11,1.648e11
mpdata,ex,2,1,2.080e11,2.076e11,2.038e11,1.955e11,1.830e11,1.664e11
mpdata,ex,3,1,2.080e11,2.076e11,2.038e11,1.955e11,1.830e11,1.664e11
mpdata,alpx,1,1,11e-6,11.115e-6,12.15e-6,14.45e-6,17.9e-6,22.5e-6
mptemp,1,100,200,400,600,800,1000
mpdata,alpx,2,1,11e-6,12.1e-6,12.65e-6,13.6e-6,14.7e-6,15.21e-6
mpdata,alpx,3,1,11e-6,12.1e-6,12.65e-6,13.6e-6,14.7e-6,15.21e-6
mp,prxy,1,1,0.33
mp,prxy,2,1,0.33
mp,prxy,3,1,0.33
mp,reft,1,10
mp,reft,2,10
mp,reft,2,10
tb,bkin,1,4
tbtemp,0,1
tbdata,1,359.25e6,2.060e11
tbtemp,20,2
tbdata,1,350.28e6,2.056e11
tbtemp,200,3
tbdata,1,269.55e6,2.019e11
tbtemp,600,4
tbdata,1,90.15e6,1.936e11
!tbcopy,bkin,1,2
lx_number=250 !x方向分段数
ly_number=75  !y方向分段数
tinc=0.2    !单步时间长度

/solu
antype,trans
outpr,basic,last
outres,basic,last
kbc,0
autots,on
alls
nsel,s,,,189507
nsel,a,,,189582
d,all,ux,0
d,all,uy,0
time,0.01
alls
solve
esel,s,mat,,2
esel,a,mat,,3
ekill,all
allsel,all
time=0.02
time,time
solve

alls
esel,s,mat,,2
*do,j,1,lx_number
!j=1
*do,i,1,ly_number
tsres,erase
xc=0.002+(j-1)*0.004
yc=0.004*(i-1)
zc=0.043
*do,enum,37501,75000
dist1=SQRT((centrx(enum)-xc)**2+(centry(enum)-yc)**2+(centrz(enum)-zc)**2)
*if,dist1,le,R,then
ealive,enum
*endif
*enddo
esel,s,live
eplot
allsel,all
time=time+tinc-0.01
time,time
ldread,temp,,,time,,,rth
!nsubst,5,20,1
lnsrch,on
!nlgeom,on
solve
time=time+0.01
time,time
bfdele,all,temp
solve
*enddo

*do,i,1,ly_number
tsres,erase
xc=0.002+(j-1)*0.004
yc=0.3-0.004*(i-1)
zc=0.043
*do,enum,37501,75000
dist1=SQRT((centrx(enum)-xc)**2+(centry(enum)-yc)**2+(centrz(enum)-zc)**2)
*if,dist1,le,R,then
ealive,enum
*endif
*enddo
esel,s,live
eplot
allsel,all
time=time+tinc-0.01
time,time
ldread,temp,,,time,,,rth
!nsubst,5,20,1
lnsrch,on
!nlgeom,on
solve
time=time+0.01
time,time
bfdele,all,temp
solve
*enddo
*enddo
time=time+3600
time,time
ldread,temp,,,time,,,rth  !(注意这里只读到最后的时间点？此处读取载荷步会不会好一点)
!nsubst,40
solve
save

time=time+600
time,time
ldread,temp,,,time,,,rth
!nsubst,20
solve

alls
esel,s,mat,,3
*do,j,1,lx_number
!j=1
*do,i,1,ly_number
tsres,erase
xc=0.002+(j-1)*0.004
yc=0.004*(i-1)
zc=0.044
*do,enum,1,37500
dist1=SQRT((centrx(enum)-xc)**2+(centry(enum)-yc)**2+(centrz(enum)-zc)**2)
*if,dist1,le,R,then
ealive,enum
*endif
*enddo
esel,s,live
eplot
allsel,all
time=time+tinc-0.01
time,time
ldread,temp,,,time,,,rth
!nsubst,5,20,1
lnsrch,on
!nlgeom,on
solve
time=time+0.01
time,time
bfdele,all,temp
solve
*enddo

*do,i,1,ly_number
tsres,erase
xc=0.002+(j-1)*0.004
yc=0.3-0.004*(i-1)
zc=0.044
*do,enum,1,37500
dist1=SQRT((centrx(enum)-xc)**2+(centry(enum)-yc)**2+(centrz(enum)-zc)**2)
*if,dist1,le,R,then
ealive,enum
*endif
*enddo
esel,s,live
eplot
allsel,all
time=time+tinc-0.01
time,time
ldread,temp,,,time,,,rth
!nsubst,5,20,1
lnsrch,on
!nlgeom,on
solve
time=time+0.01
time,time
bfdele,all,temp
solve
*enddo
*enddo
time=time+3600
time,time
ldread,temp,,,time,,,rth
!nsubst,40
solve
save
finish

赞一下

回复此楼

2楼2018-03-23 15:36:09

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

姜呃呃

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 80.4
散金: 10
帖子: 84
在线: 13.7小时
虫号: 6319137
注册: 2017-04-14
性别: GG
专业: 制造系统与自动化

!模型500mm×300mm×40mm，单元十三万左右，熔覆层大小4mm×4mm×0.2mm
!一层熔覆层单元75×126个
!定义材料参数

mp,reft,1,10
mp,reft,2,10
mp,reft,3,10
v=0.03       !熔覆线速度
lx_number=125 !x方向分段数
ly_number=75  !y方向分段数
tinc=0.13    !单步时间长度
Q=2300       !激光功率
yita=0.03    !激光效率
R=0.002    !激光有效加热半径
pi=3.1415926
qq=Q*yita/(pi*R**2)  !一个分步值

!施加初始温度
/SOLU
ANTYPE,4
TRNOPT,FULL
LUMPM,0
OUTRES,ERASE
OUTRES,BASI,LAST
EQSLV,PCG,1E-6
KBC,0.0
TIMINT,0
TRNOPT,FULL
THOPT,FULL,0.05,64, , , ,0
EQSLV,PCG, ,2, ,DELE
MSAVE,0
TUNIF,10,
TREF,10,
OUTPR,BASIC,LAST,
OUTRES,ALL,LAST,
SOLCONTROL,OFF, ,NOPL
TIME,0.13
AUTOTS,0
NSUBST,1,0,0,0
KBC,0
TSRES,ERASE
TIMINT,1
TINTP,0.005, , ,1,0.5,0.2,
SOLVE

!施加基体六个面的对流换热系数
alls
esel,s,mat,,3
esel,a,mat,,2
estif,1.0e-10
ekill,all
alls
timint,on
asel,s,,,1,6,1
nsla,s,1
sf,all,conv,-2,10
alls
!toffst,273 !不知道这个要不要加上
tintp,0.005,,,1.0,0.5,0.2
t=0.26
time,t
!nsubst,2
outpr,basic,all
outres,basic,all
alls
solve

!对基体上部熔覆层施加热流密度
alls
*do,i,1,lx_number
*do,j,1,ly_number
sfedele,all,6,hflux
tsres,erase
xc=0.002+(i-1)*0.008  !初始热源位置
yc=0.002+(j-1)*0.004
zc=0.043
*do,enum,46876,75000  !选择激活的单元
dist1=SQRT((centrx(enum)-xc)**2+(centry(enum)-yc)**2+(centrz(enum)-zc)**2)
*if,dist1,le,R,then
ealive,enum
eplot
qqq=qq
sfe,enum,6,hflux,,qqq
*endif
*enddo
alls
*do,enumm,1,46875    !选择激活熔覆层单元下面的基体单元，删除其对流
dist1=SQRT((centrx(enumm)-xc)**2+(centry(enumm)-yc)**2+(centrz(enumm)-0.042)**2)
*if,dist1,le,R,then
esel,s,,,enumm
sfdele,all,conv
*endif
*enddo
alls
*do,enumn,103126,112500  !选择相应的表面效应单元、激活、施加对流
dist1=SQRT((centrx(enumn)-xc)**2+(centry(enumn)-yc)**2+(centrz(enumn)-zc)**2)
*if,dist1,le,R,then
ealive,enumn
eplot
esel,s,,,enumn
nsle
sf,all,conv,-2,10
*endif
*enddo
alls
esel,s,live
eplot
allsel,all
lnsrch,on
!nsubst,2    !不知如何选择，选择的多了运算非常慢
t=t+tinc
time,t
*if,mod(j,15),ne,0,then !之前试过如果结果全部输出，rth文件巨大，硬盘放不下，所以选择步数输出
!outres,,none
!outpr,,none
*else
outres,basic,last
outpr,basic,last
*endif
solve
*enddo
alls

*do,p,1,ly_number       !如上，只不过是光源往回走
sfedele,all,6,hflux
tsres,erase
xc=0.006+(i-1)*0.008
yc=0.298-0.004*(p-1)
zc=0.043
*do,enum,46876,75000
dist1=SQRT((centrx(enum)-xc)**2+(centry(enum)-yc)**2+(centrz(enum)-zc)**2)
*if,dist1,le,R,then
ealive,enum
eplot
qqq=qq
sfe,enum,6,hflux,,qqq
*endif
*enddo
*do,enumm,1,46875
dist1=SQRT((centrx(enumm)-xc)**2+(centry(enumm)-yc)**2+(centrz(enumm)-0.042)**2)
*if,dist1,le,R,then
esel,s,,,enumm
sfdele,all,conv
*endif
*enddo
alls
*do,enumn,103126,112500
dist1=SQRT((centrx(enumn)-xc)**2+(centry(enumn)-yc)**2+(centrz(enumn)-zc)**2)
*if,dist1,le,R,then
ealive,enumn
eplot
esel,s,,,enumn
nsle
sf,all,conv,-2,10
*endif
*enddo
alls
esel,s,live
eplot
allsel,all
outres,basic,last
outpr,basic,last
lnsrch,on
!nsubst,2
t=t+tinc
time,t
*if,mod(p,15),ne,0,then
outres,,none
outpr,,none
*else
outres,basic,last
outpr,basic,last
*endif
solve
*enddo
*enddo
alls                   !冷却阶段
outpr,basic,all
outres,basic,all
sfedele,all,6,hflux
t=t+300
time,t
solve

赞一下

回复此楼

3楼2018-03-23 15:38:39

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

姜呃呃

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 80.4
散金: 10
帖子: 84
在线: 13.7小时
虫号: 6319137
注册: 2017-04-14
性别: GG
专业: 制造系统与自动化

三楼为第一层

发自小木虫Android客户端

回复此楼

4楼2018-03-25 08:44:53

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

匿名

用户注销 (正式写手)

应助: 0 (幼儿园)
金币: 581.1
散金: 60
红花: 6
帖子: 716
在线: 70.3小时
虫号: 0
注册: 2015-05-18
专业: 水工结构和材料及施工

本帖仅楼主可见

5楼2018-03-25 09:10:16

已阅申请仿真EPI 回复此楼编辑查看我的主页

匿名

用户注销 (正式写手)

应助: 0 (幼儿园)
金币: 581.1
散金: 60
红花: 6
帖子: 716
在线: 70.3小时
虫号: 0
注册: 2015-05-18
专业: 水工结构和材料及施工

本帖仅楼主可见

6楼2018-03-25 09:21:37

已阅申请仿真EPI 回复此楼编辑查看我的主页

geng?1015

新虫 (正式写手)

应助: 0 (幼儿园)
金币: 65.9
散金: 1300
红花: 4
沙发: 3
帖子: 595
在线: 46.3小时
虫号: 8669674
注册: 2018-05-05
性别: MM
专业: 制造系统与自动化

厉害

发自小木虫IOS客户端

回复此楼

7楼2018-05-17 16:37:21

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

liyaobo0509

专家顾问 (著名写手)

专家经验: +353
应助: 214 (大学生)
金币: 8142.8
散金: 288
红花: 64
沙发: 4
帖子: 1838
在线: 480.5小时
虫号: 911339
注册: 2009-11-24
性别: GG
专业: 机械工程
管辖: 仿真模拟

【答案】应助回帖

1.通过 ANSYS Mechanical APDL Launch 的 ANSYS Batch 环境读入 APDL代码，有更快的执行速度；

2.设置 HPC（High Performance Computing Setup）采用并行计算提高CPU利用率；

可以先采用大点的单元尺寸试算一下，看看模型有没问题~

赞一下

回复此楼

8楼2018-05-17 17:23:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主姜呃呃的主题更新

返回列表

24小时热门版块排行榜

[求助] 仿真慢 已有1人参与

» 猜你喜欢

【答案】应助回帖

[求助] 仿真慢已有1人参与