请教各位一个边界的问题

simontintin

我想模拟一个拉伸的力学实验,一般做拉伸实验的时候都不是用恒定的力,而是给试样端面一个恒定的位移速度,比如0.05mm/min,请教各位这个条件该怎么加呢?谢谢

seeksilence

在bft.for中有一个子程序bound.for是修改边值的
      real*8 function bound(r,t,j)
      implicit real*8 (a-h,o-z)
      dimension r(2)
      bound=0.0
      加入判断，是否是动边界
      是的话，修改bound值
c     write(*,*) 'bound =',bound
      return
      end
  
r是坐标，t是时间，j对应自由度（看bft.for源程序）如
      if ((abs(r(1)-700)+abs(r(2)-290)).lt.1.0e-6) then
      bound=0.05*t        
      end if

simontintin

谢谢斑竹,我试着做做.