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

you99

我想定义一个图中那种的载荷函数，需要循环的周期很大，abaqus里面自带的amplitude功能实现不了，想用vuamp子程序定义这个函数，如果只是个普通的分段函数的话我直接就用if then... else...语句了但这个是周期性的，不知道该怎么实现，恳请大佬指点一二

附上abaqus中的子程序    vuamp02.for (1.59 KB, 下载次数: 7)

2020-8-24 16:24 上传

点击文件名下载附件

you99

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

嗯嗯，多谢了

风平老涡

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

you99

风平老涡 发表于 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