I worked on base line correction before. So, I know what you're referring to. Could you give specific examples where this tool could be useful? In seismic, we typically get accelerograms from PEER or KiKNet databases that are baseline corrected.

The largest constants of integration, which are the root cause of the phenomena we are discussing here, are observed in roughly monotonic vibration signals, such as floor vibrations due to loads from pedestrians or machine operation, or vibrations observed in locomotives, which is extremely relevant for transportation of nuclear fuel. I actually didn't know PEER serves baseline corrected data, but even if they didn't you're right, seismic ground motions typically need not be baseline corrected, since the many different frequencies tend to work together to cancel out the drift constants. I've applied my baseline correction code to ground motions from PEER, and the results were almost identical to those obtained from Newmark integration alone.

However, it is my understanding that there are cases where earthquake ground motions exhibit approximately monotonic behavior, such as near-fault ground motions or the effect when the propagation of ground ripples tend to concentrate at point (I forget what you call this phenomena). If a user happened to have a raw accelerogram of these types of ground motions, they will need to perform baseline corrections. Also, let's say a user simply wants to run an analysis with a single sinusoidal function, like Asin(fpit). Any boundary where this displacement were prescribed would drift far off into space unless the correct initial conditions are set, but even then there is still some drift, since the double integral results in a C1 + C2t drift that cannot be removed with a scalar initial condition. It is sometimes the case that shake-table experiments use analytical functions like these to control the actuators, and so if a user wanted to model such an experiment with MASTODON, it might be helpful to apply this correction.

I suppose the reason why my baseline correction applied to accelerations from PEER did nothing is because they were probably already corrected 😆

Refer to this article: ftp://140.115.123.63/share/%B0%EA%A9m%AAF%B6%D5%B8%EA%AE%C6/08_References/Strong%20motion%20record%20processing%20for%20the%20PEER%20center.pdf

PEER motions are baseline corrected :)

For near fault motions, you would, in some cases, expect a nonzero displacement after the earthquake ends. This non-zero displacement is called fling step and it is due to nonlinear motion of the earth.

And so they are! This is a helpful reference, thank you, Som.

Right, if a permanent displacement is the expectation, than these corrections are not useful and, in fact, harmful if not used with care.

Hi @crswong888, thank you for the issue, and I'm glad you and @somu15 started this interesting discussion. This is definitely an important feature for seismic analysis and we need to discuss a good place for it. I'll let @somu15 figure this one out. @somu15 I'd start by asking @aeslaughter how to do this.

@cbolisetti Sure, I will check with Andrew.

@cbolisetti Thank you, Chandu, and Som as well.

@somu15 I will work with you as much as it takes to do this job correctly. Let me know when I can open a PR. Note that I have a functioning prototype that you could easily try out and I provided a link in my original message.

@crswong888 Please consider a few things:

• Add this tool in your MASTODON fork.
• In src/vectorprocessors and include/vectorpostprocessors create the .C and .h files, respectively.
• Create documentation in doc/content/source/vectorpostprocessors.

Before doing any the above, I am not clear what are the inputs and outputs, if your tools has to be used after or before a simulation. Please clarify these.

No problem. I will do this and report back when I get the chance.

@somu15 okay I added the items you requested in the latest commit to my MASTODON fork. Obviously, this commit needs a lot of fine-tuning, but the state of this branch currently serves the purpose of making clear to you how this vpp object works. I would like to add some more documentation to the theory manual describing the mathematical formulation behind this object, if you think this is appropriate. I would also like to eventually add more tests with varying combinations of the input parameters.

Regarding your question about whether this vpp should be used before or after the simulation, the answer is: both. It depends on what you are trying to accomplish. Please note that, if the user wanted to use the corrected data as a boundary condition, say with PresetAcceleration, getting it to work properly takes a bit of hackery in the input file. It is actually currently not possible in the master versions of the MOOSE Framework, but I have several PR's and Issues open to address this problem. Ultimately, I would like to create an action system that simplifies this for the user, but I figured I would deal with this in a future MASTODON PR, after all of the basic requirements in the Framework have been addressed.