Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

Frequency domain analysis with two different frequencies

Please login with a confirmed email address before reporting spam

Hello,

I don't know if the title makes sense. Currently, with a frequency domain analysis I can evaluate steady dynamic response of a system with a set of frequencies for an harmonic load. Instead of sweeping through frequencies I would like to run the study but with this time, loads to which I have set different frequencies.

I tried running a frequency domain analysis with no harmonic load but with instead manual time varying loads (sin) and I get no result. It only produces constant stress and displacement.

In a nutshell I want to run frequency domain analysis in these two cases:

  • case #1 (already possible): one harmonic load through freqs [1,10,100] Hz

  • case #2 (target) : two harmonic load at once, with respectively freqs 9Hz and 11Hz.

Has anyone any idea on how to do that?

I must mention that a working example for two loads will need to be extended to up to 50 loads.

Thanks in advance for any tip


4 Replies Last Post 10.01.2019, 17:19 GMT-5
Henrik Sönnerlind COMSOL Employee

Please login with a confirmed email address before reporting spam

Posted: 6 years ago 10.01.2019, 03:00 GMT-5
Updated: 6 years ago 10.01.2019, 03:01 GMT-5

Hi Thomas,

Your question is quite interesting, and there is not a single short answer.

First, you must notice that the problem you describe may or may not be periodic at all depending on the selected frequencies. If the ratio between all pairs of frequencies is a rational number, then a periodic solution exists. The period may however be long when compared to the periods of the contributing loads.

In your case, the periods of the two loads are 1/9 s and 1/11 s. The whole process is however periodic with the period 1 s, so in order to evaluate the full solution you need to combine 9 cycles of the 9 Hz solution with 11 cycles of the 11 Hz solution. Clearly, this is nothing that can be done by a pure frequency domain approach.

What you need to do is to solve the harmonic problem for each of the constituent frequencies, and then combine the results in the time domain. An example of how that can be done (for two frequencies, but the principle is general) can be found in one of the studies in https://www.comsol.com/model/vibration-analysis-of-a-thick-beam-20301 .

If you are only interested in an upper bound to the response, you can skip looking at the result in time domain, and instead just sum the amplitudes of the result for all frequencies.

If the process is not is not truly periodic (say that you have the frequencies 1 Hz and sqrt(2) Hz), you can always get an acceptable accuracy by assuming periodicity, and for example setting the second frequency to 1.41 or 1.414. For each decimal you add, the period of the combined result becomes a factor of 10 larger.

Regards,
Henrik

-------------------
Henrik Sönnerlind
COMSOL
Hi Thomas, Your question is quite interesting, and there is not a single short answer. First, you must notice that the problem you describe may or may not be periodic at all depending on the selected frequencies. If the ratio between all pairs of frequencies is a rational number, then a periodic solution exists. The period may however be long when compared to the periods of the contributing loads. In your case, the periods of the two loads are 1/9 s and 1/11 s. The whole process is however periodic with the period 1 s, so in order to evaluate the full solution you need to combine 9 cycles of the 9 Hz solution with 11 cycles of the 11 Hz solution. Clearly, this is nothing that can be done by a pure frequency domain approach. What you need to do is to solve the harmonic problem for each of the constituent frequencies, and then combine the results in the time domain. An example of how that can be done (for two frequencies, but the principle is general) can be found in one of the studies in . If you are only interested in an upper bound to the response, you can skip looking at the result in time domain, and instead just sum the amplitudes of the result for all frequencies. If the process is not is not truly periodic (say that you have the frequencies 1 Hz and sqrt(2) Hz), you can always get an acceptable accuracy by assuming periodicity, and for example setting the second frequency to 1.41 or 1.414. For each decimal you add, the period of the combined result becomes a factor of 10 larger. Regards, Henrik

Please login with a confirmed email address before reporting spam

Posted: 6 years ago 10.01.2019, 09:27 GMT-5

Hi Henrik,

Thanks for this tip ! I looked into the example you proposed and located the study that correspond well to what you were explaining and what I am looking for (overimposed frequency responses, time wise). It's called "periodic forced vibration" if I understand right, the 3rd one in 4.

One thing I do not get though, is through which steps the frequency responses are calculated and overimposed before extracting the time response with Inverse Discrete Fourier Transform. It seems there is only one harmonic load and the other one is not harmonic but has its phase shifting as a function of frequency. I was trying to get one of the two frequency only to be applied (as to perfectly understand the process) but it does not seem to work. Time wise, I'm only able to get two frequencies, or none.

In a nutsheel: inside Step 1 : Frequency Domain, only Edge Load 2 is enabled (the one with the phase) and in Step 2: Frequency to Time FFT all loads are disabled, but still I'm getting two frequencies in the time visualization.

Would you be please kind enough as to explain what I'm not getting here?

Thanks in advance

Hi Henrik, Thanks for this tip ! I looked into the example you proposed and located the study that correspond well to what you were explaining and what I am looking for (overimposed frequency responses, time wise). It's called "periodic forced vibration" if I understand right, the 3rd one in 4. One thing I do not get though, is through which steps the frequency responses are calculated and overimposed before extracting the time response with Inverse Discrete Fourier Transform. It seems there is only one harmonic load and the other one is not harmonic but has its phase shifting as a function of frequency. I was trying to get one of the two frequency only to be applied (as to perfectly understand the process) but it does not seem to work. Time wise, I'm only able to get two frequencies, or none. In a nutsheel: inside **Step 1 : Frequency Domain**, only Edge Load 2 is enabled (the one with the phase) and in **Step 2: Frequency to Time FFT** all loads are disabled, but still I'm getting two frequencies in the time visualization. Would you be please kind enough as to explain what I'm not getting here? Thanks in advance

Henrik Sönnerlind COMSOL Employee

Please login with a confirmed email address before reporting spam

Posted: 6 years ago 10.01.2019, 11:15 GMT-5

Hi,

If you look at the settings for Step1: Frequency Domain you will see that the only active load in that study step is Edge Load 2.

It is not tagged as Harmonic Perturbation since that study step is not of the perturbation type. In that study, an 'ordinary' load is considered as harmonic. If you were to change to Linear perturbation in the solver settings (see screenshot), then the solver would respond to the Harmonic Perturbation tag. You can compare with Modal Solver 1 in the previous study, which uses Edge Load 1.

Regards,
Henrik

-------------------
Henrik Sönnerlind
COMSOL
Hi, If you look at the settings for **Step1: Frequency Domain** you will see that the only active load in that study step is **Edge Load 2**. It is not tagged as **Harmonic Perturbation** since that study step is not of the perturbation type. In that study, an 'ordinary' load is considered as harmonic. If you were to change to **Linear perturbation** in the solver settings (see screenshot), then the solver would respond to the **Harmonic Perturbation** tag. You can compare with **Modal Solver 1** in the previous study, which uses **Edge Load 1**. Regards, Henrik


Dave Greve Certified Consultant

Please login with a confirmed email address before reporting spam

Posted: 6 years ago 10.01.2019, 17:19 GMT-5

In case #2 you want to calculate the response when there are two sources (forces) acting at the same time at different frequencies.

IF THE SYSTEM IS LINEAR you can do this by calculating the response at the two frequencies and summing them in the time domain.

If the system is not linear- the most obvious approach is to perform two time dependent calculations with sinusoids as forcing functions. However the calculated response will include both the forced response and an initial transient. The transient may take a considerable time to die away in order for you to get the steady state response (if that IS what you want).

In case #2 you want to calculate the response when there are two sources (forces) acting at the same time at different frequencies. IF THE SYSTEM IS LINEAR you can do this by calculating the response at the two frequencies and summing them in the time domain. If the system is not linear- the most obvious approach is to perform two time dependent calculations with sinusoids as forcing functions. However the calculated response will include both the forced response and an initial transient. The transient may take a considerable time to die away in order for you to get the steady state response (if that IS what you want).

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.