Get Help

Turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Community
- :
- PTC Mathcad
- :
- PTC Mathcad
- :
- Duhamel integral for a damped SDOF with mathcad

Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Highlighted

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

01-23-2018
01:37 PM

01-23-2018
01:37 PM

to all,

I am trying to set up the standard convulsion /Duhamel integral for a damped SDOF subject to an arbitrary input p(t) but can’t figure out the syntax required in MathCAD to make it work. Could someone help me on this one? (see attached mcad file)

Thanks

Regards

Solved! Go to Solution.

Labels:

1 ACCEPTED SOLUTION

Accepted Solutions

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

01-24-2018
10:07 AM

01-24-2018
10:07 AM

I think part of your problem stems from the definition of your impulse function p. You impulse frequency is much higher than the natural frequency of your system, and the pulse is over very fast. But the top illustration in your image shows an excitation imposed on a system with a higher frequency than the excitation.

Attached is an example with the same damped system, with a half-sine pulse with a frequency much lower than the natural frequency of your system. And the integral works (sort of.)

If you had an impulse defined by data points you would have to fit a function to the data before you integrated.

11 REPLIES 11

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

01-23-2018
02:07 PM

01-23-2018
02:07 PM

Re: Duhamel integral for a damped SDOF with mathcad

Mathcad REALLY likes units.

See if this helps

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

01-23-2018
02:13 PM

01-23-2018
02:13 PM

Re: Duhamel integral for a damped SDOF with mathcad

CONVOLUTION! please! You have to define x.0 and v.0

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

01-23-2018
03:29 PM

01-23-2018
03:29 PM

Re: Duhamel integral for a damped SDOF with mathcad

While sure Fred is right about units you see in his sheet that you have to strip out the units for the odesolve-block (one of the very few advatages of Prime is that you can use units here.).

But maybe we should explain why your attempt failed. The reason is that in your initial conditions you used x0 and v0, but forgot to define them.

EDIT: Sorry, Francesco already explained that. His reply did not show up until I posted mine.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

01-24-2018
02:53 AM

01-24-2018
02:53 AM

Re: Duhamel integral for a damped SDOF with mathcad

Thanks Fred, Francesco and Werner for the replies

Convolution indeed! Will call it Duhamel from now on

I didn't use unit because all the inputs are consistent (unit wise) and I did notice that the odesolve didn't like units. I also realized that I used m to define the mass

Thanks for pointing the error on the way I defined the Duhamel integral

Questions:

1. Why is the Duhamel integral which should give the response x(t) different to the x(t) obtained with the odesolve?

2. I am testing with a well defined p(t) function. How does one handle a function define with data points? given time step I suspect I need to define t[i

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

01-24-2018
08:45 AM

01-24-2018
08:45 AM

Re: Duhamel integral for a damped SDOF with mathcad

...Indeed, In all branches of engineering and physics, is used the term - convolution - ......

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

01-24-2018
10:07 AM

01-24-2018
10:07 AM

I think part of your problem stems from the definition of your impulse function p. You impulse frequency is much higher than the natural frequency of your system, and the pulse is over very fast. But the top illustration in your image shows an excitation imposed on a system with a higher frequency than the excitation.

Attached is an example with the same damped system, with a half-sine pulse with a frequency much lower than the natural frequency of your system. And the integral works (sort of.)

If you had an impulse defined by data points you would have to fit a function to the data before you integrated.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

01-24-2018
11:31 AM

01-24-2018
11:31 AM

Re: Duhamel integral for a damped SDOF with mathcad

Thanks. I just plugged numbers to make it work. Will need to have a look at the sheet provided at a slow pace.

I don't like setting up units and then having to remove them. Looks awful !

Haven't had a look a the detail yet but it's surprising that one gets different answer. Looks very good up to t=0.1 then after that the ODE solve is going "crazy" !

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

01-24-2018
12:20 PM

01-24-2018
12:20 PM

Re: Duhamel integral for a damped SDOF with mathcad

wrote:I don't like setting up units and then having to remove them. Looks awful !

Haven't had a look a the detail yet but it's surprising that one gets different answer. Looks very good up to t=0.1 then after that the ODE solve is going "crazy" !

I work a lot with actual hardware, so I need/like to have units to define things. But you have seen this problem set up without units and made to work. The issue is that the scaling has to be right, or the answers are wrong; if (like me) you are proscribing units in whatever comes in handy converting to a uniform system (aka "UnitsOf") is required for the correct answer.

As far a who goes crazy, the Duhamel calculation generated the spikes:

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Notify Moderator

01-25-2018
03:56 AM

01-25-2018
03:56 AM

Re: Duhamel integral for a damped SDOF with mathcad

Apologies. Indeed the Duhamel "goes" crazy. As a general interest: What is causing this "spiky" behaviour?