Skip to main content
O
OLUBUKOLAADEBOY10-Marble
10-Marble
August 11, 2017
Question

Wilson Theta Method with Matrix

  • August 11, 2017
  • 1 reply
  • 2820 views

Good Evening Mathcad Community,

 

About a week ago the following people have been help me with programming from translating Matlab into Mathcad on Wilson Theta Method, Werner, FM and Jamie. Now that the problem is solved, I'm trying to read in matrices Mass matrix, Stiffness matrix, Damping matrix and an external force vector into the function but get comfused with the n by n matrix into the function. Would any member of the community please see what I'm doing wrong and help me the code is below the function. I have followed Werner's work on previous worksheet. If I get this to work my goal is to solve the problem below the function on dynamic analysis for foundations. Please attached.

1 reply

W
Werner_E25-Diamond I
25-Diamond I
August 11, 2017

I don't have the time to dig into it in more depth, so just a loose collection of observations.

 

  • In Matlab there seems to be a function with four input matrices and four output vectors. I would implement as a function in Mathcad, too.
  • You initialize u,an,v as square matrices but you need and treat them later as simple vectors - that causes syntax error
  • The section "% Initial conditions" from the Matlab program seems to be missing in your implementation
  • The expression (R-R) in the Matlab code looks strange. This would just be a vector with all zeros. So apart from that R is used in the (missing) initial conditions it has no further influence on the calculation.

Good luck!

 

O
OLUBUKOLAADEBOY10-MarbleAuthor
10-Marble
August 11, 2017

Good Morning Werner,

 

You are correct, I missed the initial conditions. I have modified the function to include the initial condition. Also, I do not know how the author came about the (theta*(R-R)) in the equation. I checking my numerical analysis text on the topic to find out. My main problem is do I need two for loop or one for the matrix?. Please help.

 

Best Regards,

 

Boyer

W
Werner_E25-Diamond I
25-Diamond I
August 11, 2017

As you see you get an error because the various array dimensions don't match.

You define u for example as a 301 x 301 matrix. The Maltlab code defines the corresponding disp as a 2 x 301 matrix and seems to anly use the first row of it. Not sure whats going on here. Guess u etc. can be just vectors.

Terms & ConditionsAccessibility statement