MATLAB Answers

0

How can i solve this integration using newmark?

Asked by Leonardus Risky on 12 Aug 2018 at 14:11
Latest activity Commented on by Leonardus Risky on 13 Aug 2018 at 1:57

Hello everybody, I now create a code for solving the dynamic equation of motion in the time response using Newmark-Beta Method but it can't run because there's a problem using loop (for). you can use this reference for understanding the meaning of Newmark-beta Method. can anyone help me to solve this code? Thank You for your help.

  0 Comments

Sign in to comment.

1 Answer

Answer by Stephan Jung on 12 Aug 2018 at 16:41
Edited by Stephan Jung on 12 Aug 2018 at 16:43

Hi,

in your for loop is an end missing. I guess that this missing end should be positioned before your plot commands.

Then the following line throws an error:

   dis(i+1) = delF/Keff;

I think it should be:

dis(i+1) = delF(i)/Keff;

or:

dis(i+1) = delF(i+1)/Keff;

After fixing this there are some errors appearing due to different lenghts of x- and y-values whole plot is running, which is a result from your for loop and should be solved by you. Note that sbx has a length of 20001.

If you solved this then your plots will overwrite eachother, so a recommendation would be to use subplot to show all results in one figure.

I hope this helped, if not come back with more questions.

Best regards

Stephan

  1 Comment

thanks for your advise, i think i solved my problem with making new for loop. but anyway thank you

Best Regards

Leonardus

Sign in to comment.