buddyd16 / pyMAOS

Code built from Matrix Analysis of Structures by Aslam Kassimali
BSD 3-Clause "New" or "Revised" License
4 stars 1 forks source link

Derivation PDF #1

Open EngDM opened 1 year ago

EngDM commented 1 year ago

Apologies in advance if this is missuse of the issues form.

I am trying to work my way through the linear load derivation, what exactly does the third page represent? Is this an equation for deflection at a point? Is there anything in that derivation would be deflection at point "x" for a linear load?

buddyd16 commented 1 year ago

Page 3 of that pdf is the simplified symbolic solutions for the integration constants.

Towards the bottom of page 1 you should see the word Deflection below that is the piecewise functions for deflection at a point for the first two bounds x=0->a, x=a->b. At the top of page 2 is the last bound of x=b->L.

On Sat, Aug 12, 2023 at 3:48 AM EngDM @.***> wrote:

Apologies in advance if this is missuse of the issues function.

I am trying to work my way through the linear load derivation, what exactly does the third page represent? Is this an equation for deflection at a point? Is there anything in that derivation would be deflection at point "x" for a linear load?

— Reply to this email directly, view it on GitHub https://github.com/buddyd16/pyMAOS/issues/1, or unsubscribe https://github.com/notifications/unsubscribe-auth/AC6XJSW47OYQVWHXZD7UNMTXU4YLVANCNFSM6AAAAAA3NYZOZM . You are receiving this because you are subscribed to this thread.Message ID: @.***>

EngDM commented 1 year ago

So given that there are C# constants in that equation, does this imply that there is no set piecewise function that could be implemented into a spreadsheet or something?  

buddyd16 commented 1 year ago

No, the equations shown in the derivation are the set piecewise functions. The formulas for constants C1-C12 are indicated on page 3.

a bit confused on what you are not understanding the derivation is laid out going from load functions to deflection functions via direct integration. Then initial boundary and compatibility conditions are defined. Then there are a set of 12 equations defined by applying the boundary/compatibility conditions these need to solved simultaneously for the integration constants c1-c12. Page 3 is the resulting solution for the integration constants, these then get plugged into the piecewise functions defined on page 1.

For the implementation in python see the class beginning on line 407 in the loading.py file.

EngDM commented 1 year ago

Ah I see, so on page 3 that is what C1-C12 are? I’m assuming broken up by line. Sorry, it’s not labeled and I’m unfamiliar with this derivation. Thanks

buddyd16 commented 1 year ago

“ Ah I see, so on page 3 that is what C1-C12 are? ”

yes, broken up by line.

No worries, the derivation itself is fairly straightforward (assuming you’re familiar with basic differential equations and how to apply boundary conditions) the equations just get huge for the linear load case so the hard part is making sure you have typed everything in correctly.

EngDM commented 1 year ago

Yea I see that looking at it further. It’s just massive since it’s a catch all for all UDL loads, sloped or constant over any portion of the beam. Anyways, I hope I haven’t been too much of a bother. Thanks,

buddyd16 commented 1 year ago

No bother at all, sometimes it’s a bit difficult to pin point how exactly to best help with the limited context of these messages. Happy to help if I am able.

EngDM commented 1 year ago

Hey again,

I'm in the middle of implementing the point load derivation, however I am noticing that it is only accurate when a>L/2 (or a>b). I've compared my values with 2 different online solvers, and they align when a>b.

I couldn't see anywhere in the derivation that said a has to be larger than b, however in the simplified equations online they are all listed only when xb.

buddyd16 commented 1 year ago

a bit confused there is no b variable in the point load derivation? the point load is fully located by a single position variable, a.

EngDM commented 1 year ago

Say b = a-L.

EngDM commented 1 year ago

image

Also, the bottom of the derivation appears to indicate that C6 should be positive?

image

buddyd16 commented 1 year ago

so no the initial equation is written with the assumption of a positive C6, the resulting solution shows that the value is actually (-P*a^3)/6.

I've run hand checks and checks against other software and am getting agreement with my functions?

EngDM commented 1 year ago

Ah that makes sense.

I'll do some more investigating, might have to re-enter everything. Given that it is only when a<x, it's probably an issue with the constants associated with it.

Just to clarify, it is for any variation of this loading and not just x<L/2 like some other published equations?

buddyd16 commented 1 year ago

My derivation is for the point load located at any value of a.

buddyd16 commented 1 year ago

my python class for the point load begins on line 318 here: https://github.com/buddyd16/pyMAOS/blob/main/pyMAOS/loading.py