大家来看看是不是气体保护焊

19850913

/filn,welding
/config,nres,1000000
/prep7
/units,si
et,1,plane55
et,2,solid70
mptemp,1,20,200,400,800,1200,1500
mpdata,kxx,1,1,13.31,19.47,25.07,27.53,29.76,31.95
mpdata,c,1,1,470,550,634,676,719,765
mpdata,dens,1,1,7966,7814,7630,7535,7436,7320
mpdata,enth,1,1,7.8e8,1.67e9,2.75e9,4.34e9,5.26e9,6.27e9
*dim,cnvtab,table,6,1,1,temp
cnvtab(1,0)=20.0,200.0,400.0,800.0,1200.0,1500.0
cnvtab(1,1)=18.85,33.44,55.95,88.92,134.6,273.4
*afun,deg
width=0.3               
height=0.012               
long=0.3
bf=3
k,1,0,0,0
k,2,long,0,0
k,3,long,width/2,0
k,4,0,width/2,0
a,1,2,3,4
wpoffs,0,0.021,0
wprota,0,90,0
asbw,1
lesize,1,,,75,,,,,1
lesize,3,,,75,,,,,1
lesize,7,,,75,,,,,1
lesize,5,,,7,,,,,1
lesize,6,,,7,,,,,1
lesize,8,,,5,0.3,,,,1
lesize,9,,,5,0.3,,,,1
mshkey,0
amesh,2
amesh,3
/view,1,1,1,1
/replot
extopt,aclear,1
extopt,attr,0
extopt,esiz,bf
mat,1
type,2
vext,2,3,,,,height
nummrg,all                     ! 合并所有实体
numcmp,all                     ! 压缩编号
/image,save,vmesh,jpeg         !保存图片为jpeg格式
save
wprota,0,-90,0
finish
/solu
antype,4                     !声明分析类型
tunif,20
kbc,1                        !   指定载荷是递增载荷或阶跃载荷0 递增 1阶跃
timint,1                    !   打开瞬态效应 时间积分
autots,0                  !   指定是否使用自动时间步长跟踪或载荷步跟踪
asel,s,area,,1,3
asel,a,area,,5
asel,a,area,,7,11
sfa,all,1,conv,%cnvtab%,20
nsel,s,loc,y,0,0.006
nsel,r,loc,z,0.008,0.012
nplot
esln,r,1
eplot
ekill,all
save
!*********************************************焊接参数
v=0.3/75!焊接速度4mm/s   
u=30  !电压
i=240 !电流
y=0.8 !热利用率
a1=0.000050 !焊缝面积
hgen1=y*u*i/(a1*v) !焊缝热生成率
!*********************************************
xx=0
*do,tm,1,75
  time,tm
  deltim,,,,1
  tsres,erase
  nsel,s,loc,y,0,0.006
  nsel,r,loc,z,0.008,0.012  
  nsel,r,loc,x,xx,xx+0.004
  esln,r,1
  ealive,all
  eplot
  bfe,all,hgen,1,hgen1
  allsel,all
  
  solve
  allsel,all
  
  bfe,all,hgen,1,0
  outres,all,all
  xx=xx+0.004
*enddo
tm=tm+100
time,tm
nsubst,20
outres,all,all
solve
tm=tm+200
time,tm
nsubst,20
outres,all,all
solve
tm=tm+500
time,tm
nsubst,25
outres,all,all
solve
save
finish
!残余应力的计算
/filn,structure
/prep7
etchg,tts
mptemp,1,20,200,400,800,1200,1500
mpdata,ex,1,1,2.09e11,1.95e11,1.72e11,1.33e11,0.84e11,0.45e11
mpdata,nuxy,1,1,0.267,0.322,0.262,0.229,0.223,0.223            !   泊松比
mpdata,alpx,1,1,15.24e-6,17.44e-6,18.83e-6,19.38e-6,19.95e-6,20.7e-6 !   线膨胀系数
   
tb,bkin,1,6             !经典双线性随动强化参数（材料1）
tbtemp,20,1  
tbdata,1,3.15e8,2.09e10
tbtemp,200,2
tbdata,1,2.78e8,1.95e10
tbtemp,400,3
tbdata,1,2.25e8,1.72e10
tbtemp,800,4
tbdata,1,1.53e8,1.33e10
tbtemp,1200,5  
tbdata,1,9.3e7,0.84e10
tbtemp,1500,6   
tbdata,1,3.3e7,0.45e10
finish
/solu
tref,20
nlgeom,on
nropt,full
lnsrch,on
neqit,10000
cnvtol,f,,0.01,2
kbc,0
outres,all,12
allsel,all
da,4,symm
d,P51X, , , , , ,UX,UZ, , , ,
D,P51X, , , , , ,UZ, , , , ,

mnsub1=12
*do,mt,1,300,1
  ldread,temp,,,mt,,'welding',rth
  nsubst,mnsub1
  time,mt
  solve
*enddo
*do,mt,301,500,2
  ldread,temp,,,mt,,'welding',rth
  nsubst,mnsub1
  time,mt
  solve
*enddo
*do,mt,501,875,2
  ldread,temp,,,mt,,'welding',rth
  nsubst,mnsub1
  time,mt
  solve
*enddo
save
finish

19850913

怎么没有人来顶一下