Mrennui
Structural
- Oct 14, 2010
- 5
I will try and keep this brief
I'm trying to calculate the force on impact of a tile falling onto a glass roof using the work-energy method. I understand that there many many unknown variables and I have expalined this to the client and he requested simplified results anyhow.
I know the kinetic energy at impact 1/2mv². Then using work-energy the average force can be calculated from F=1/2mv²/d (later d = v). Where d is the displacement or deflection of the glass.
I then supposed that the maximum displacemt could be found using the engineers bending equation (1) M/I = E/R = -s/y
Where s is the stress, in this case the limiting stress(the example attached is for a steel beam). I supposed using max stress would give max deflection.
Knowning that (2) 1/R = d²v/dx² where v is the vertical displacement (page 324 Megson, I know this may not be strictly valid as glass can be subeject large deflections, but ignoring this)
Putting (2) into (1) Ey d²v/dx² = s
Integrating twice (I have attached my example)
Ey v = sx²/2 + s L/2 x
which gives a deflection of 41.9mm
I checked this supposing Mmax = smax Z
Then Mmax = WL/4 giving Wmax = 4Mmax/L
Finally v = WmaxL³/48EI
But this gives v as -27.9mm
I freely admit I'm not the best engineer or mathematician, can someone explain to me what I'm doing wrong? Is it flawed methodology, have I just got the number wrong? In the sheets attached I was troubled that v came out as positive. Any idea's are very much appreicated.
I'm trying to calculate the force on impact of a tile falling onto a glass roof using the work-energy method. I understand that there many many unknown variables and I have expalined this to the client and he requested simplified results anyhow.
I know the kinetic energy at impact 1/2mv². Then using work-energy the average force can be calculated from F=1/2mv²/d (later d = v). Where d is the displacement or deflection of the glass.
I then supposed that the maximum displacemt could be found using the engineers bending equation (1) M/I = E/R = -s/y
Where s is the stress, in this case the limiting stress(the example attached is for a steel beam). I supposed using max stress would give max deflection.
Knowning that (2) 1/R = d²v/dx² where v is the vertical displacement (page 324 Megson, I know this may not be strictly valid as glass can be subeject large deflections, but ignoring this)
Putting (2) into (1) Ey d²v/dx² = s
Integrating twice (I have attached my example)
Ey v = sx²/2 + s L/2 x
which gives a deflection of 41.9mm
I checked this supposing Mmax = smax Z
Then Mmax = WL/4 giving Wmax = 4Mmax/L
Finally v = WmaxL³/48EI
But this gives v as -27.9mm
I freely admit I'm not the best engineer or mathematician, can someone explain to me what I'm doing wrong? Is it flawed methodology, have I just got the number wrong? In the sheets attached I was troubled that v came out as positive. Any idea's are very much appreicated.