cancel
Showing results for 
Search instead for 
Did you mean: 
cancel
Showing results for 
Search instead for 
Did you mean: 

Community Tip - Learn all about the Community Ranking System, a fun gamification element of the PTC Community. X

Bending Moment and Shear Force Diagram Issues

MikeArmstrong
5-Regular Member

Bending Moment and Shear Force Diagram Issues

I have started drafting some material for student engineers relating to bending moments and shear forces. Initially,  the SF and BM were calculated long hand and the plotted graphs look correct. However, when integration and differentiation is used the graphs are not correct. I am using M15 with Windows 10.

 

Thanks,

 

Mike

1 ACCEPTED SOLUTION

Accepted Solutions

I broke the definite integral into three parts and the defect disappears, look here:

MA answer 2.jpg

View solution in original post

16 REPLIES 16

Hi, 

I will go on so:

MA answer.jpg

 

MikeArmstrong
5-Regular Member
(To:-MFra-)

Hi,

I have tried this and still get the same error. I have just modified the functions and swapped the 'otherwise' statements to 'if' as shown below:

This has corrected the differentiation error when plotting the shear force diagrams, but there is still an error with the bending moment (MM(x)) plot when integrating V(x).

Thanks,

Mike

MikeArmstrong
5-Regular Member
(To:-MFra-)

Nice work and the graphs are look correct. However, it still isn't clear why the bending moment formula in my worksheet (MM(x)) is returning a zero value for MM(0.9599*m) and above.

 

Mike

Hi,

Bending Moment & Shear Diagrams

Have you checked the Mathcad Roark's ebook? MM & SF diagrams are plotted for some typical examples

MikeArmstrong
5-Regular Member
(To:JBlackhole)

Thanks for the response. Yes I have these but wanted to draft the examples from first principles. The only issue I have is the discontinuity with the bending moment when integrating the shear force.

 

Mike 

I broke the definite integral into three parts and the defect disappears, look here:

MA answer 2.jpg

MikeArmstrong
5-Regular Member
(To:-MFra-)

Good work and thanks for the effort.

 

Could you please post the worksheet so I can investigate further?

 

Mike

Voilà:

MikeArmstrong
5-Regular Member
(To:-MFra-)

Thank you.

That is a very clever workaround and something I have not seen before. I actually don’t think I have seen the Phi fuction (not sure if that’s the correct name) used ever.

Is the issue in my example due to a discontinuity in the shear force assessment?

Mike

1) Welcome back, Mike!



Is the issue in my example due to a discontinuity in the shear force assessment?

2) No, its the splitting of the integral in three parts (which shouldn't be necessary) which does the trick. It works OK with the way you defined V(x), too. Have no idea whats going on an why the result suddenly is 0. Tweaking TOL did not change anything.

B.png

BTW, using the Heaviside function Phi makes the definition a bit less clear but has the advantage that you can evaluate a function defined that way symbolically, too, if necessary.

MikeArmstrong
5-Regular Member
(To:Werner_E)

Thanks Werner and nice example. I have to agree that the Heaviside function Phi makes it a little less clear, but certainly a good tool to have.

The behaviour of the integration suddenly going to zero is certainly strange, especially as differentiating the bending moment for shear force works.

 

Clipboard01.jpg

Mike


The behaviour of the integration suddenly going to zero is certainly strange, especially as differentiating the bending moment for shear force works.

 

Yes, and the Mathcad help only succinctly recommends to split the integration into parts which Francesco did and could make it work.

B.png

I tried decreasing TOL and/or changing the method of integration, but in vain.

BTW, you should replace the expression "4wa-wx" in your definition of the bending moment by "w*(L-x)" as otherwise it would yield wrong results if you change "a" from 0.25m to some other value.

MikeArmstrong
5-Regular Member
(To:Werner_E)

Noted and well spotted. I was leaving it as shown to display how Mathcad simplifies equations in respect to the textbook examples.

Mike

F(x) is a special function often called Heaviside step function or simply unitary step. Its derivative is the Dirac pulse.  For example, you can plot functions like the one here depicted:

staircase.jpg

 For the other question, see the answers of the other community members.

MikeArmstrong
5-Regular Member
(To:-MFra-)

Thanks for the explanation and help, much appreciated.

 

Mike

Top Tags