声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1120|回复: 0

[综合讨论] 一条焊缝的命令流程序,请高手指教!

[复制链接]
发表于 2008-11-11 15:39 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
大家好
我编写了一条焊缝的焊接命令流,但是只是加热阶段的程序,冷却阶段的老是处错误,希望有高手指教下,本人将不胜感激!补充一点,我用的是直接法耦合的,下面是我焊接加热阶段的命令流程序,运行没有错误。(注意记得glue一下)
/PreP7
mPtemP,l,20,250,500,750,1000,1500,1700,2500
mpdate,kxx,l,l,50,47,40,27,30,35,140,142
mpdata,dens,l,l,7820,7700,7610,7550,7490,7350,7300,7090
mpdata,c,1,l,460,480,530,675,670,660,780,820
MPDATA,PRXY,l,l,0.28,0.29,0.31,0.35,0.4,0.45,0.48,0.5
mPdata,alPx,l,l,1.10E-05,1.22E-05,1.39E-05,l.48E-05,l.34E-05,l.33E-05,1.32E-05,1.3lE-05
mPdata,ex,l,l,2.05E+11,1.87E+11,1.50E+11,7.00E+10,2.00E+10,1.90E+10,1.80E+10,l.20E+10
MPTEMP,,,,,,,,
MPTEMP,l,20
MPTEMP,2,1533
MPTEMP,3,1590
MPTEMP,4,1670
MPDATA,ENTH,l,,0
MPDATA,ENTH,l,,7.5E+009
MPDATA,ENTH,l,,9.6E+009
MPDATA,ENTH,l,,1.1E+010
MPDATA,MURX,1,,1
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,250
MPTEMP,3,500
MPTEMP,4,750
MPTEMP,5,1000
MPTEMP,6,1500
MPTEMP,7,1700
MPTEMP,8,2500
MPDATA,HF,1,,1
MPDATA,HF,1,,3.5
MPDATA,HF,1,,5.2
MPDATA,HF,1,,10
MPDATA,HF,1,,15
MPDATA,HF,1,,30
MPDATA,HF,1,,31
MPDATA,HF,1,,35
!下面在各温度点下的屈服应力和屈服后的弹性模量
TB,BISO,l,6,2,
TBTEMP,20
TBDATA,,2.20E+08,2.05E+10,,,,
TBTEMP,500
TBDATA,,1.30E+08,l.50E+10,,,,
TBTEMP,750
TBDATA,,4.00E+07,7.00E+09,,,,
TBTEMP,1000
TBDATA,,2.50E+07,2.00E+09,,,,
TBTEMP,1700
TBDATA,,1.00E+06,1.80E+09,,,,
TBTEMP,2500
TBDATA,,1.00E+05,1.20E+09,,,,
SAVE
/prep7
/units,si
Rb=0.002
Length=0.12
width=0.03
height=0.003
et,l,solid5
et,2,solid98
block,0,length,0,width,0,height
block,0,length,0,3*Rb,0,height
block,0,length,0,6*Rb,0,height
block,0,length,0,12*Rb,0,height
VOVLAP,ALL
/view,1,1,1,1
!下面定义焊缝区的单元尺寸
Lesize,4,rb
Lesize,5,rb
Lesize,9,rb
Lesize,14,rb
Lesize,50,rb
Lesize,10,rb
Lesize,52,rb
Lesize,23,rb
Lesize,51,rb
Lesize,9,rb
Lesize,49,rb
Lesize,24,rb
!定义过渡区的单元尺寸
Lesize,31,2*rb
Lesize,26,2*rb
!定义远离焊缝区的单元尺寸
Lesize,43,4*rb
Lesize,38,4*rb
Lesize,2,4*rb
Lesize,7,4*rb
Lesize,53,4*rb
Lesize,54,4*rb
Lesize,55,4*rb
Lesize,56,4*rb
Lesize,11,1.5*rb
Lesize,12,1.5*rb
Lesize,47,1.5*rb
Lesize,48,1.5*rb
Save
记得glue一下
Type,1
Mshkey,1
Mshape,0
Vmesh,5
Vmesh,6
Type,2
Mshkey,0
Mshape,1
Vmesh,7
Vmesh,8
Save

/SOLU
DA,3,
/GO
DL,53,,ux,
DL,61,,UX,
DL,57,,UX,
DL,49,,UX,
DL,4,,uz,
DL,9,,uy,
/Solu
!时间步长设置
Antype,trans,new!
Nlgeom,on!
Timint,1,struct!
Timint,1,therm!
Timint,0,mag!
Timint,0,elect
TintP,0.005,,,-1,0.5,0.2
Tref, 25
Nropt,auto,,
Autos,on
Outres,,all
v=0.01
tinc=rb/v
!移动热源设置
U=25
I=180
ef =0.75
P=25*180*ef
r=0.006
qmax=p/(3.14159*r**2)
Max_time=length/rb+1
max_x=length/rb+1
max_y=width/rb+1
*dim,flux2,table,max_x,max_y,max_time,x,y,time
*do,k,1,max_time,1
*do,i,1,max_x,1
Flux2(i,0,k)=(i-1)*rb
*enddo
*do,j,1,max_y,1
Flux2(0,j,k)= (j-1)*rb
*enddo
*do,i,1,max_x,1
*do,j,1,max_y,1
!设置热源中心位置坐标
xcenter=v*(k-1)*tinc
ycenter=0
distance=sqrt(((i-1)*rb-xcenter)**2+((j-1)*rb-ycenter)**2)
*if,distance,le,r,then
Flux2(i,j,k)=qmax/EXP(3*DISTANCE**2/(R**2))
*else
Flux2(i,j,k)=0
*endif
*enddo
*enddo
Flux2(0,0,k)=(k-1)*tinc
Nsel,s,loc,z,height
Sf,all,hflux,%flux2%
Allsel,all
Nsel,s,loc,x,0
Sf,all,conv,10,20
Allsel,all
Nsel,s,loc,x,length
Sf,all,conv,10,20
Allsel,all
Nsel,s,loc,y,0
Sf,all,conv,10,20
Allsel,all
Nsel,s,loc,y,width
Sf,all,conv,10,20
Allsel,all
Nsel,s,loc,z,0
Sf,all,conv,10,20
Allsel,all
Time,k*rb/v
Nsubst,5
solve
*enddo

save
下面冷却阶段的命令流哪位高手能给指教一下,谢谢了!
回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-30 19:05 , Processed in 0.096118 second(s), 23 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表