abaqus里面如何通过Fortran编译周期性的分段函数？

风平老涡

[Fortran] 纯文本查看 复制代码

real function p(t, t1, t2, p0)
implicit none
real, intent(in) :: t, t1, t2, p0
real :: r

r = mod(t, t1 + t2)
if (r <= t1 ) then
   p = p0
else
   p = 0
end if
return
end function p

风平老涡

中心思想就是把任意时间 t 转换成（T1+T2）的余数，然后根据余数判断P值。