Skip to main content
F
fkcom1-Visitor
1-Visitor
October 19, 2018
Question

Newmark Integration

  • October 19, 2018
  • 1 reply
  • 4756 views

Dear helpers,

 

I am currently trying to solve the reaction of a laminated plate by using the Newmark-Integration-Scheme based on the First Order Shear Deformation Theory (the force acting on the plate is sinusoidal in time). Therefor I wrote a loop which calculates the needed Wmn for every next time step dependent on values of the current timestep.

 

My problem is that this calculation needs a lot of time at the moment (especially when I try to plot the last equation) so I can't see if my calculations are correct. 

 

My question is if there is any mistake I did by programming the loop because this is the first time I was programming in mathcad. 

 

I am thankful for every advice.

Thanks,

Felix

 

 

1 reply

-
-MFra-21-Topaz II
21-Topaz II
October 19, 2018

Halo,

In den Mathcad-Programmen können die Werte der rechten Seite nicht den Funktionen von zwei Variablen auf der linken Seite zugeordnet werden. Stattdessen kann dies getan werden, wenn auf der linken Seite ein Element der Matrix ist. Siehe Fotos:

fkcom.jpg

    F
    fkcom1-VisitorAuthor
    1-Visitor
    November 8, 2018

    Hello,

     

    thank you very much for your reply and sorry for my late answer.

    I understand what you mentioned and I tried it out (see attached file). Now the loop should work better.

    Unfortunately I'm still not getting the results that I want. Is there another mistake you can see in the file?

     

    Thank you very much in advance!

     

    deutsch:

    Besten Dank für deine Antwort und sorry für meine späte Rückmeldung.

    Deinen Vorschlag habe ich verstanden und umgesetzt (s. angehängte Datei). Ich denke, die Schleife sollte jetzt funktionieren. 

    Leider bekomme ich immer noch nicht die Ergebnisse, die ich gern hätte. Ist da noch etwas, was ich übersehen haben könnte? Meine Mathcad-Kenntnisse sind nicht die besten, vielleicht habe ich noch einen ähnlichen Fehler drin wie der bereits von dir erwähnte.

     

    Besten Dank vorab.

     

    -
    -MFra-21-Topaz II
    21-Topaz II
    November 8, 2018

    Hi fkcom,
    at the beginning of the file it is good to define the origin of the indexes, in this case one. Then I defined the intervals for the x and y. It seems that it works, even if the times are very long when you choose high-order matrices in the "createmash" arguments. (10x10 is fast enough while 20x20 takes etwa 30 seconds)

    fkcom1.jpg

      Terms & ConditionsAccessibility statement