风平老涡 发表于 2020-8-24 22:13
中心思想就是把任意时间 t 转换成（T1+T2）的余数，然后根据余数判断P值。

嗯嗯，多谢了

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

风平老涡 发表于 2020-8-23 04:51
[mw_shl_code=fortran,true]real function p(t, t1, t2, p0)
implicit none
real, intent(in) :: t, t1, t2 ...

感谢回复，我从网上查了一下说您这个属于自定义函数的用法，而我这个是从abaqus仿真软件中定义一个pressure，再通过子程序定义这个pressure的幅值，所以在编写的时候还是稍微有点差别，我按照您这个简单又编写了一遍，在附件里面有，可以帮忙看一下我编写的对吗？

[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