[eliocpereira]

I would like to know how to impose an oscillating linear motion with a sudden stop to a body. I know how to impose an oscillating linear motion. I wrote the following lines in the dynamicMeshDict:

dynamicFvMesh   dynamicMotionSolverFvMesh;

motionSolver    solidBody;

solidBodyMotionFunction multiMotion;

oscillation
{
    solidBodyMotionFunction oscillatingLinearMotion; //Sinusoidal translation motion 
    oscillatingLinearMotionCoeffs  
    {
        amplitude     (0.02 0 0);  //amplitude vector of the sinusoidal translation [m]
        omega         21.9838;     //angular frequency of the sinusoidal translation [rad/s]
    }
}

Let us consider now a sudden stop at instant t=t_s. How could I express this in this code? For example, is it possible to use a step function at the amplitude entry, so that it goes zero when t=t_s?

[simrego]

Hi!


What if you try the timeActivatedFIleUpdate function object:
https://cpp.openfoam.org/v6/timeActi...8H_source.html
Just create 2 dynamicMeshDicts, and at the given time just copy the 2nd one.
I'm not sure if OF reads the dynamicMeshDict if it is modified. If not, the much easier way is to run the simulation to the time when the motion should stop, manually (through a script) change the dynMesDict file, and the end time of the simulation. And continue the run. I think this is not a professional move, but simple and working.

[eliocpereira]

I added this code to my controlDict file inside "functions":

    fileUpdate1                  //substitute a file at a given instant
    {
        type              timeActivatedFileUpdate;
        libs ("libutilityFunctionObjects.so");
        writeControl     timeStep;
        writeInterval    1;
        fileToUpdate      "$FOAM_CASE/constant/dynamicMeshDict"; //file to be substituted
        timeVsFile
        (
            (-1  "$FOAM_CASE/constant/dynamicMeshDict.1")        //substitution at t=0s
            (0.3 "$FOAM_CASE/constant/dynamicMeshDict.2")        //substituition at t=0.3s
        );
    }

I also placed my dynamicMeshDict.1 and dynamicMeshDict.2 files in the "constant" paste. I found out that the dynamicMeshDict file isn't re-read after the changes. Is there any option that needs to be added to the fileUpdate1 function so that OpenFOAM could read again the dynamicMeshDict file?

[eliocpereira]

I now tried to change the oscillatingLinearMotion.C and oscillatingLinearMotion.H files, so that they could do what I desired. I defined a new SBMFCoeffs variable: ts, which is the stopping instant of the oscillation. I changed the formula for the oscillation, using comparison operators:

const vector displacement = amplitude_*sin(omega_*t)*(t<ts_) + amplitude_*sin(omega_*ts_)*(t>=ts_);

You can see the new oscillatingLinearMotion.C and oscillatingLinearMotion.H files here ( and the dynamicMeshDict file as well):https://www.dropbox.com/sh/9mueknp5y...zBWamXcEa?dl=0

I discovered that the changes to the files actually didn't change anything. It felt like OpenFOAM used the old files instead of the new ones, although the old were substituted. I even tried a simple change: only putting a factor 2 multiplying the frequency in the sine function, so that I could expect a two times faster oscillation. Still, the results were the same as before, nothing changed. Why isn't OpenFOAM reading my changes?.

[eliocpereira]

I eventually made OpenFOAM recognize the changes that I did to the .C and .H files, by running the related makefile after. Basically I changed directory to where the "Make" folder lives, then I logged in as root and ran the makefile. To do this, I wrote the following commands at the terminal:

cd  /opt/openfoam6/src/dynamicMesh
sudo bash
wmake

[simrego]

Sorry for the late answer. So yes. That's why I told you in this case you can try to simply stop the simulation and the continua with the changed motion.
The next approach would be the new dynMesh solver, or a more better one if you use a coded one.
You should NOT modify the original code. You can do it, but next time you forget that, you use this motion solver, and you will have no idea what the hack is happening. :/
In this case it is better if you create a new library for yourself. (copy paste original, rename, rewrite, you can use.)