diff options
Diffstat (limited to '3776/CH12/EX12.5/Ex12_5.sce')
| -rw-r--r-- | 3776/CH12/EX12.5/Ex12_5.sce | 27 |
1 files changed, 27 insertions, 0 deletions
diff --git a/3776/CH12/EX12.5/Ex12_5.sce b/3776/CH12/EX12.5/Ex12_5.sce new file mode 100644 index 000000000..c7bb91a72 --- /dev/null +++ b/3776/CH12/EX12.5/Ex12_5.sce @@ -0,0 +1,27 @@ +clear +//Given +//Virtual loading Without f_d +p_ab = -0.833 //lb The recorded virtual loading +p_bc = + 0.833 //lb The recorded virtual loading +F_ab = 2500 //lb +F_bc = -2500 //lb +l_ab = 60 //in - The length of the rod +l_bc = 60 //in - The length of the rod +A_ab = 0.15 //in2 the areaof ab +A_bc = 0.25 //in2 the areaof bc +E = 30*(10**6) //psi The youngs modulus of the material +//Part_a +e_a =p_ab*l_ab*F_ab/(A_ab*E) + p_bc*l_bc*F_bc/(A_bc*E) //lb-in the deflection +//With f_d +p_bd = 1 //lb The recorded virtual loading +F_bd = 1 //lb +l_bd = 40 //in - The length of the rod +A_bd = 0.1 //in2 the areaof ab +e_a_1 =p_ab*p_ab*l_ab/(A_ab*E) + p_bc*p_bc*l_bc/(A_bc*E) +p_bd*p_bd*l_bd/(A_bd*E) //lb-in the deflection +//Since the produced defelection should compensate the other one +x_d = e_a/e_a_1 +printf("\n The reaction force at D is %0.2f lb",-x_d) + +//Part - B +e_b = -x_d*l_bd/(A_bd*E ) //in - The deflection of nodal point B +printf("\n The deflection of nodal point B %0.4f in",e_b) |