# Random and harmonic vibration analysis

**URL:** <https://prepomax.discourse.group/t/random-and-harmonic-vibration-analysis/52>\
**Category:** Feature Requests\
**Created:** [February 10, 2022, 5:03am UTC](https://prepomax.discourse.group/t/random-and-harmonic-vibration-analysis/52 "2022-02-10T05:03:31Z")\
**Posts on this page:** 12\
**Page:** 1

<div class="post-metadata">

**Author:** ![chumairo](https://avatars.discourse-cdn.com/v4/letter/c/848f3c/32.png) [@chumairo](https://prepomax.discourse.group/u/chumairo)\
**Post date:** [February 10, 2022, 5:03am UTC](https://prepomax.discourse.group/t/random-and-harmonic-vibration-analysis/52/1 "2022-02-10T05:03:31Z")

</div>

Dear All,

Is there any plan to add dynamic analysis to prepomax? I normally use FEA for modal analysis and propose structure modification to avoid resonance vibration problem. However, it would be nice if we could predict dynamic stress from measured vibration data as well in some cases which natural frequency could not be shifted away from excitation sources with sufficient margin. With predicted dynamic stress, it can be compared against S-N curve of material for fatigue failure analysis. Also, with proposed structure modification, dynamic stress and vibration can be predicted to see if any improvement compared to the existing ones.

With current version, can we edit keywords to do above analysis? Any guidance or suggestion would be highly appreciated.

Regards,

---

<div class="post-metadata">

**Author:** ![FEAnalyst](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/feanalyst/32/688_2.png) [@FEAnalyst](https://prepomax.discourse.group/u/FEAnalyst)\
**Post date:** [February 10, 2022, 9:11am UTC](https://prepomax.discourse.group/t/random-and-harmonic-vibration-analysis/52/2 "2022-02-10T09:11:34Z")

</div>

Unfortunately, random vibration analyses (as well as response spectrum ones) are not supported by CalculiX. When it comes to linear dynamics analyses, only transient modal dynamic and harmonic (steady-state) vibration simulations can be performed: [CalculiX: Overview of the finite element capabilities of CalculiX Version 2.18](http://www.dhondt.de/ov_calcu.htm)

In fact, steady-state dynamic analyses are one of the most requested features for PrePoMax. Currently, you could use keywords but it would be quite complicated. It will be somewhat easier if you apply load in form of a concentrated force acting at a reference point connected to the rest of the model with rigid body constraint.

 ![ssd](https://global.discourse-cdn.com/free1/uploads/prepomax/original/1X/2e5696b10eb87a18a4c3732866dbcfce14bc9955.png)

However, it seems that the results from this type of analysis can’t be displayed in PrePoMax (likely because of their nature - they are presented in the frequency domain:

 ![ssd 2](https://global.discourse-cdn.com/free1/uploads/prepomax/original/1X/2b48f79c8ee4dd980ad43cf42b8a90b32fab6615.png)

You would have to use CalculiX GraphiX to visualize such results.

---

<div class="post-metadata">

**Author:** ![chumairo](https://avatars.discourse-cdn.com/v4/letter/c/848f3c/32.png) [@chumairo](https://prepomax.discourse.group/u/chumairo)\
**Post date:** [February 10, 2022, 10:29am UTC](https://prepomax.discourse.group/t/random-and-harmonic-vibration-analysis/52/3 "2022-02-10T10:29:35Z")

</div>

Thank you a lot for your guidance. Do I have to add one more step of Frequency analysis before Steady state dynamic analysis?

---

<div class="post-metadata">

**Author:** ![FEAnalyst](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/feanalyst/32/688_2.png) [@FEAnalyst](https://prepomax.discourse.group/u/FEAnalyst)\
**Post date:** [February 10, 2022, 10:59am UTC](https://prepomax.discourse.group/t/random-and-harmonic-vibration-analysis/52/4 "2022-02-10T10:59:23Z")

</div>

Yes because it’s a mode-based procedure. A preceding Frequency step with storage enabled has to be included but this can be defined without the need to use keywords.

---

<div class="post-metadata">

**Author:** ![JakaP](https://avatars.discourse-cdn.com/v4/letter/j/da6949/32.png) [@JakaP](https://prepomax.discourse.group/u/JakaP)\
**Post date:** [September 25, 2022, 11:14am UTC](https://prepomax.discourse.group/t/random-and-harmonic-vibration-analysis/52/5 "2022-09-25T11:14:05Z")

</div>

Hi all,  
are there any plans to include steady state harmonic analysis in the PrePoMax in the near future?

---

<div class="post-metadata">

**Author:** ![lucas\_bueno](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/lucas_bueno/32/1168_2.png) [@lucas\_bueno](https://prepomax.discourse.group/u/lucas_bueno)\
**Post date:** [September 26, 2022, 12:34am UTC](https://prepomax.discourse.group/t/random-and-harmonic-vibration-analysis/52/6 "2022-09-26T00:34:08Z")

</div>

It would be a nice feature to have. I have been using the steady-state dynamics procedure for some time, and it would be nice to at least see the results in PrePoMax for more detailed post-processing.

---

<div class="post-metadata">

**Author:** ![Matej](https://avatars.discourse-cdn.com/v4/letter/m/a698b9/32.png) [@Matej](https://prepomax.discourse.group/u/Matej)\
**Post date:** [September 27, 2022, 12:54pm UTC](https://prepomax.discourse.group/t/random-and-harmonic-vibration-analysis/52/7 "2022-09-27T12:54:28Z")

</div>

The next step will be the support for explicit analysis.

If there are result files that are currently not supported by PrePoMax please let me know. Adding support for result files takes much less work than adding the UI for the preparation of the unsupported procedures.

---

<div class="post-metadata">

**Author:** ![FEAnalyst](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/feanalyst/32/688_2.png) [@FEAnalyst](https://prepomax.discourse.group/u/FEAnalyst)\
**Post date:** [September 27, 2022, 1:55pm UTC](https://prepomax.discourse.group/t/random-and-harmonic-vibration-analysis/52/8 "2022-09-27T13:55:14Z")

</div>

> [@Matej](#):
>
> If there are result files that are currently not supported by PrePoMax please let me know.

It seems that the results from SSD analyses can’t be postprocessed in PrePoMax, likely because they are in frequency domain:

> [@Random and harmonic vibration analysis](https://prepomax.discourse.group/t/random-and-harmonic-vibration-analysis/52/2):
>
> Unfortunately, random vibration analyses (as well as response spectrum ones) are not supported by CalculiX. When it comes to linear dynamics analyses, only transient modal dynamic and harmonic (steady-state) vibration simulations can be performed: [CalculiX: Overview of the finite element capabilities of CalculiX Version 2.18](http://www.dhondt.de/ov_calcu.htm) In fact, steady-state dynamic analyses are one of the most requested features for PrePoMax. Currently, you could use keywords but it would be quite complicated. It will be …

I attached the exemplary file:  
[SSD test.pmx](https://prepomax.discourse.group/uploads/short-url/fccY3GU2DSybCGj7bAQhqWc1eB5.pmx) (168.7 KB)

---

<div class="post-metadata">

**Author:** ![Matej](https://avatars.discourse-cdn.com/v4/letter/m/a698b9/32.png) [@Matej](https://prepomax.discourse.group/u/Matej)\
**Post date:** [October 1, 2022, 6:14pm UTC](https://prepomax.discourse.group/t/random-and-harmonic-vibration-analysis/52/9 "2022-10-01T18:14:53Z")

</div>

I have tried to make some progress on this topic. First of all, I could not get the analysis to finish from PrePoMax. Using ccx 2.18 and Pastix solver gave me an error. I changed the solver to Spooles, and it was a little better. The analysis terminates with the bottom monitor message, but actually, the .frd file is about 14 MB in size and contains only a part of the solution.

 ![image](https://global.discourse-cdn.com/free1/uploads/prepomax/original/1X/6bf7508c8be888764a4e75255ef90d3aef834d4c.png)

So I tried running Calculix 2.18 from the command line with the default solver and got the following error:

 ![image](https://global.discourse-cdn.com/free1/uploads/prepomax/original/1X/2dd4a254e2dbab76b0db7bdd4ec0c7c6869fffff.png)

Then I changed the version to 2.19. Run the analysis from PrePoMax. Using the default solver (Pastix) it now runs till the end. But not always. Sometimes the analysis terminates with an error. Running the same version from the command line also works most of the time but not always. When it does not work, I get the error about POSIX WinThreads.

Then I changed the version to 2.20. Run the analysis from PrePoMax. Using the default solver (Pastix) it now runs till the end. But only sometimes. Most times, the analysis terminates with an error. Running the same version from the command line also works sometimes. When it does not work, I get the error about POSIX WinThreads.

**Can someone confirm this?**

After successful analysis, the results could not be read due to NAN values in the imaginary fields. I have now changed the .frd reader to account for that. The results can now be read but are recognized and shown as ordinary static results.

 ![image](https://global.discourse-cdn.com/free1/uploads/prepomax/original/1X/0805c5e1b0d9912aa426bf7bc2039ee7df9b5a3d.png)  
 ![image](https://global.discourse-cdn.com/free1/uploads/prepomax/original/1X/5763142b1b6d4ec2740028bbd07de88b4aba679e.png)

---

<div class="post-metadata">

**Author:** ![FEAnalyst](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/feanalyst/32/688_2.png) [@FEAnalyst](https://prepomax.discourse.group/u/FEAnalyst)\
**Post date:** [October 1, 2022, 9:33pm UTC](https://prepomax.discourse.group/t/random-and-harmonic-vibration-analysis/52/10 "2022-10-01T21:33:09Z")

</div>

I also observe such inconsistent behavior when running this input file using the three mentioned versions of CalculiX. The analysis either terminates at some stage (usually the same as in your case - “Composing the steady state response from the eigenmodes” or during factorization) with no error messages (I don’t get that pop-up error regarding POSIX WinThreads) or is completed successfully. However, the difference is that I run it in parallel when submitting from the command line (I set the OMP\_NUM\_THREADS variable using cmd previously).

When submitting from PrePoMax, I sometimes get the “Job finished” message but often the output is the same as in your case ("Perturbation parameter is active. Process elapsed time: …). Again with parallelization.

What’s interesting, I tried with the beamdy8.inp file from ccx test examples and it runs successfully in most cases but sometimes fails during factorization (in all the 3 versions of CalculiX).

I can report it on the CalculiX forum.

---

<div class="post-metadata">

**Author:** ![Matej](https://avatars.discourse-cdn.com/v4/letter/m/a698b9/32.png) [@Matej](https://prepomax.discourse.group/u/Matej)\
**Post date:** [October 2, 2022, 2:26pm UTC](https://prepomax.discourse.group/t/random-and-harmonic-vibration-analysis/52/11 "2022-10-02T14:26:26Z")

</div>

I think it is worth reporting it if it is not only the problem with my computer.

---

<div class="post-metadata">

**Author:** ![FEAnalyst](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/feanalyst/32/688_2.png) [@FEAnalyst](https://prepomax.discourse.group/u/FEAnalyst)\
**Post date:** [October 3, 2022, 6:51pm UTC](https://prepomax.discourse.group/t/random-and-harmonic-vibration-analysis/52/12 "2022-10-03T18:51:17Z")

</div>

I’m trying to report it on the CalculiX forum but each time I click “Create Topic”, I get error 404. I will keep trying, maybe it’s some temporary issue with that forum.
