How can i solve this integration using newmark?

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.

Answer 1

[Stephan]

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

[Leonardus Risky]

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

Best Regards

Leonardus