Continue to Site

Eng-Tips is the largest engineering community on the Internet

Intelligent Work Forums for Engineering Professionals

  • Congratulations SE2607 on being selected by the Eng-Tips community for having the most helpful posts in the forums last week. Way to Go!

Defining matrix variable to get numerical solution to the Winkler foundation differential equation

Status
Not open for further replies.

gmprox23

Civil/Environmental
Aug 12, 2022
3
Following the solution to the differential equation for the deflection curve of a beam on an elastic foundation, am trying to get the bending moment M and shearing force V based on the 2nd- and 3rd-order derivatives.

I can come up with values based on the derivative expression BUT am unable to assign the results to matrix variables M and V for plotting the results to show variation along the length of the beam.

I tried to define matrices with matching size, BUT it won't evaluate the differential equation.

Should I use "for loop" instead? Or am I missing something? Can't seem to figure it out.

Thank you for the help.


1222_Winkler_foundation_PTC_query.pdf_vsnlyy.jpg

 
Replies continue below

Recommended for you

First of all, yes that looks as though it should work, but since it doesn't I'd guess the vectorising operator is confusing things, I don't think Mathcad does Matlab styl structures. If you use an range index for lambda (and possibly M) it might work.

Cheers

Greg Locock


New here? Try reading these, they might help FAQ731-376
 
Hi Greg,

I tried to implement a 'for loop' and defined the range index for the v (=distance x) and bending moment M, BUT it still won't assign the values to the matrix variable M (pls see below).

Not sure why the matrix variable M = zero, unlike the results of the 2nd order derivatire of the differential equation.

1222b_Winkler_foundation_PTC_query__vectorized_M_ufrpv2.jpg
1222c_Winkler_foundation_PTC_query__vectorized_M__shearing_force_3rd_order_derivative_flskdb.jpg
 
Sorry, beyond my pay grade, I'll leave it to someone with a current version of mathcad to comment. It still looks to me as though you are trying to assign multiple values to each row of M.

Cheers

Greg Locock


New here? Try reading these, they might help FAQ731-376
 
I’ve done a few beam on elastic foundation problems in the past. I’ve found the use of displacement, gradient, moment and shear equations (differentials of the displacement equation) as functions of an independent variable x along the beam, incorporating the unknown constants in each of the equations, a good way to solve the problem. I use a solve block (Given and Find) with usually defined gradient, moment and / or shear equations set to known boundary conditions at defined values of x, to solve for the unknown constants of the general solution equation.
 
Status
Not open for further replies.

Part and Inventory Search

Sponsor