JesseBonanno / IndeterminateBeam

A solver for 1D indeterminate beams
MIT License
69 stars 20 forks source link

Maximum Bending Moment / Shear Force #35

Open jujus98 opened 1 year ago

jujus98 commented 1 year ago

Hello,

first of all, thank you for this wonderful package!

I want to calculate the maximum bending moment and shear force for a specific section of the beam. Therefore I tried to use the following code:

# Generate x-coordinates with a step size of 1
x_coords = list(range(6))
# Calculate the maximum bending moment for the section of x coordinates
max_bending_moment = beam.get_bending_moment(*x_coord ,return_max=True) 

Unfortunately, the program always returns the absolute maximum of the entire beam. Can someone help me with this issue? I would be very grateful!

Best regards Julius

jujus98 commented 1 year ago

I found a solution, but I thought there may be a more elegant one:

# Generate x-coordinates with a step size of 0.1
x_coords = [i * 0.1 for i in range(int(front_fuselage / 0.1) + 1)]
# Calculate maximum bending moment for each x-coordinate
max_bending_moments_ff = beam.get_bending_moment(*x_coords)
# Calculate the overall max bending moment for this section
overall_max_bending_moment_ff = abs(max(max_bending_moments_ff, key=abs))
JesseBonanno commented 1 year ago

Hi Julius,

I think your solution is the same way i would go about it. You could put the beam.get_bending_moment in the list comprehension but thats essentially the same thing.

I just read through the documentation and for the beam.get_bending_moment function return_max = True isnt doing what you think it is. In reality the get_bending_moment function works as a wrapper this parameter shouldnt have been made as a transparent option for the user, and it will always return the max for the beam. I think i should rewrite this function to be similar to how you expected it to work.

Cheers, Jesse

jujus98 commented 1 year ago

Hi Jesse,

thank you very much for your fast answer!

Best regards Julius

jujus98 commented 1 year ago

Hi Jesse,

I hope that you are doing well!

Is there an option to remove the deflection graph from the output?

Best regards and thank you very much Julius