Castigliano's second theorem states that the partial derivative of the total strain energy U with respect to a force P applied at a point is equal to the displacement δ at that point in the direction of the force. δ=∂P∂U Since support B is a support, the vertical displacement at B is zero. Therefore, we can say: ∂RB∂U=0, where RB is the reaction at support B. First, we introduce a redundant reaction RB at support B. Now, we can find the reactions at supports A and C. Sum of moments about C=0 (taking clockwise moments positive): RA∗10+RB∗4−10∗6−2∗4∗2=0 10RA+4RB−60−16=0 10RA=76−4RB RA=7.6−0.4RB Sum of vertical forces = 0 (taking upward forces positive):
RA+RB+RC−10−2∗4=0 RA+RB+RC=18 RC=18−RA−RB=18−(7.6−0.4RB)−RB=10.4−0.6RB Now, let us define the bending moment equations for the two spans: AB and BC. Span AB (0≤x≤6) measured from A: M1=RAx=(7.6−0.4RB)x for 0≤x≤4 M1=RAx−10(x−4)=(7.6−0.4RB)x−10x+40 for 4≤x≤6 Span BC (0≤x≤4) measured from C: M2=RCx−2x2/2=(10.4−0.6RB)x−x2 The total strain energy due to bending is given by:
U=∫0L2EIM2dx According to Castigliano's theorem:
∂RB∂U=∫0LEIM∂RB∂Mdx=0 Since EI is constant, we have: ∫0LM∂RB∂Mdx=0 ∫06M1∂RB∂M1dx+∫04M2∂RB∂M2dx=0 Now, let's compute the required partial derivatives:
∂RB∂M1=−0.4x for 0≤x≤4 ∂RB∂M1=−0.4x for 4≤x≤6 ∂RB∂M2=−0.6x Thus,
∫04((7.6−0.4RB)x)(−0.4x)dx+∫46(((7.6−0.4RB)x−10x+40)(−0.4x)dx+∫04((10.4−0.6RB)x−x2)(−0.6x)dx=0 −∫04(3.04x2−0.16RBx2)dx−∫46(3.04x2−0.16RBx2−4x2+16x)(−0.4x)dx−∫04(6.24x2−0.36RBx2−0.6x3)dx=0 −∫04(3.04x2−0.16RBx2)dx−∫46(−1.216x3+0.064RBx3+1.6x3−6.4x2)dx−∫04(6.24x2−0.36RBx2−0.6x3)dx=0 −∫04(3.04x2−0.16RBx2)dx−∫46(0.384x3+0.064RBx3−6.4x2)dx−∫04(6.24x2−0.36RBx2−0.6x3)dx=0 −[3.04x3/3−0.16RBx3/3]04−[0.384x4/4+0.064RBx4/4−6.4x3/3]46−[6.24x3/3−0.36RBx3/3−0.6x4/4]04=0 −[3.04(64)/3−0.16RB(64)/3]−[0.384(324−256)+0.064RB(324−256)−6.4(216−64)/3]−[6.24(64)/3−0.36RB(64)/3−0.6(256)/4]=0 −65.17+3.413RB−[0.384(68)+0.064RB(68)−6.4(152)/3]−[133.12−7.68RB−38.4]=0 −65.17+3.413RB−[26.112+4.352RB−324.27]−[94.72−7.68RB]=0 −65.17+3.413RB−26.112−4.352RB+324.27−94.72+7.68RB=0 138.298+6.741RB=0 RB=−138.298/6.741=−20.516 However, the calculation is too long and tedious and thus prone to error. Instead of this approach, let us analyze the effect of removing support B. If we remove the support, the deflection is d=48EIPL3+384EI5qL4. We then equate this displacement to the displacement at B due to the reaction at B, which is upwards. This simplifies the calculations significantly. Another simplified approach to estimate the reaction at B is to consider that the problem is roughly equal to one where the entire beam is equally loaded with 2kN/m. So, totalload=2kN/m∗10m=20kN. Then, RB=totalload/2=10kN. The 10kN load would then introduce additional downward reaction at B, which may be roughly 5-10 kN as well. Without EI value, we cannot obtain an accurate numerical solution.