24小时热门版块排行榜    

查看: 1410  |  回复: 0

焱涼

新虫 (初入文坛)

[求助] ansys模拟,关于时间求解,如何修改

/title,Weld Analysis by "Element Birth and Death"
/prep7
/unit,si
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
x1=0.05
y1=0.05
hlimit=2
D=2
H=0.015
account=100
numcmp,all
esel,all
vsel,s,volume,,1
eslv
*get,emax,elem,,num,max
*get,emin,elem,,num,min
*get,nse,elem,,count
*dim,ne,,nse
*dim,nex,,nse
*dim,ney,,nse
*dim,nez,,nse
*dim,mmz,,nse
*dim,mzz,,nse
mine=0
*do,i1,1,nse
esel,u,elem,mine
*get,nse1,elem,,count
ii=0
*do,i,emin,emax
*if,esel(i),eq,1,then
ii=ii+1
ne(ii)=i
*endif
*enddo
j=1
*do,i,1,nse1
*set,mmz(i),0
*set,mzz(i),0
*enddo
*do,i,1,nse1
*get,nex(i),elem,ne(i),cent,x
*get,ney(i),elem,ne(i),cent,y
*get,nez(i),elem,ne(i),cent,z
*enddo
*do,i,1,nse1
*if,nez(i),eq,h,then
*set,mmz(j),i
j=j+1
*enddo
*do,i3,1,account
*set,nem(i),0
*set,mzz(i3),mmz(i3)
*enddo
dt=1
t=0
esel,all
eplot
/auto,1
/replot
*do,i,1,nse
ekill,ne(i)
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
t=t+dt
time,t
*do,i3,1,ccount
k=mz(i3)
m=D*exp(-((nex(k)-x1)**2+(ney(k)-y1)**2))
nem(k)=nem(k)+m
*if,nem(k),ge,hlimit,then
ne(k),alive
nem(k)=nem(k)-hlimit
k1=k-1
nem(k1)=nem(k)
mz(i3)=k1
*endif
*enddo
solve
finish
/post1
重点在于!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!号后面的这个如何编程
建模部分未给出,是0.1*0.1*0.1的立方体,剖分的单元是0.01的立方体块
我要实现的是在在每个dt算一次m值,并且累加,满足大于hlimit的条件的单元就激活
现在关于求解这个部分完全不会。求帮忙
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 焱涼 的主题更新
信息提示
请填处理意见