Discussion Forum

Hello,

I'm running a large batch of time-dependent models. The simulation time should change for each model depending on the parameters I input, so right now I'm using a really long solution time and a stop condition that looks like:

if(t > 0.003 && mod1.focalPointDisplacement < 10^-7, -1, 1)

Such that when the domain point probe "focalPointDisplacement" falls below 10^-7, the simulation stops and I can carry on with the next one.

The problem is that for some inputs, *all* displacement will be below 10^-7. However, what I really want is to stop it after the values have "relaxed", as the focalPointDisplacement will increase to some maximum then drop down to (essentially) 0 (at which point I want to stop it). Ideally, I want to use a stop condition like

if(mod1.focalPointDisplacement < (0.05 * max(mod1.focalPointDisplacement)), -1, 1)

where "max(mod1.focalPointDisplacement)" is the maximum value recorded by that probe over time, up till the current solution time. Of couse, "max" doesn't work, and using a maxop doesn't seem to work either: maxop1(solid.disp) seems to always be exactly equal to mod1.focalPointDisplacement when they're both defined at the same location in space.

So, any way to grab the maximum of a probe over time and use it in a stop condition?

Thanks so much for any help you can offer!