Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Load Combination Bending Moment Diagram #121

Closed
egyangelo opened this issue Jan 21, 2023 · 2 comments · Fixed by #133
Closed

Load Combination Bending Moment Diagram #121

egyangelo opened this issue Jan 21, 2023 · 2 comments · Fixed by #133
Labels

Comments

@egyangelo
Copy link

egyangelo commented Jan 21, 2023

I noticed that bending moment diagram for load combination is not presented accurate and shown as linear value, while it should be 2nd degree curve, although the values are correct.

Bending Moment for Dead Load
image

Bending Moment for Live Load
image

Bending Moment for LoadCombination (D+L)
image

@smith120bh smith120bh added the bug label Jan 24, 2023
@smith120bh
Copy link
Collaborator

Thanks, @egyangelo ! I'll take a look at this shortly (probably this weekend). Looks like the load combination is adding the nodal bending moments, but not the elements' bending moments.

@smith120bh
Copy link
Collaborator

@egyangelo : Apologies for this taking a bit longer than planned! It is fixed in the just-released version 1.4.1 of anastruct. Please let me know if you still have any issues or if you find any other problems, or have any other suggestions for anastruct!

I tested that this works now using the following code:

system = SystemElements()

system.add_element([[0, 0], [2, 0]])
system.add_element([[2, 0], [4, 2]])
system.add_element([[4, 2], [6, 2]])
system.add_support_hinged([1, 4])

lc_dist = LoadCase("dist")
lc_dist.q_load(q=-1, element_id=[1, 2, 3])
system.apply_load_case(lc_dist)
system.remove_loads()

lc_point = LoadCase("point")
lc_point.point_load(node_id=[2, 3], Fy=-2)

combination = LoadCombination("combo")
combination.add_load_case(lc_dist, factor=2)
combination.add_load_case(lc_point, factor=1)
results = combination.solve(system)
results["combination"].show_bending_moment()

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

Successfully merging a pull request may close this issue.

2 participants