Tygra_1983
Student
- Oct 8, 2021
- 121
Hi guys,
Back again!
I am attempting to solve a pinned support gable frame using the unit load method.
I am using BeamGuru software, but I want to check the results to esure the program is working properly.
I have replaced the pinned supporta at A with a roller and it looks like this:
The height to eaves is 8 metres, the height to the apex is 9.58 metres. The frame is 30 metres long and the roof is at an angle of 6 degrees.
I have been using Octave to perform the integrations:
For HA, from these calculations I am getting 111.90 kN. The software is giving me 98.6 kN. So, Im a bit off. I think I could be getting mixed up with signs, but I have checked and checked, and cannot see why I am getting the incorrect result.
Many Thanks!
Back again!
I am attempting to solve a pinned support gable frame using the unit load method.
I am using BeamGuru software, but I want to check the results to esure the program is working properly.
I have replaced the pinned supporta at A with a roller and it looks like this:
The height to eaves is 8 metres, the height to the apex is 9.58 metres. The frame is 30 metres long and the roof is at an angle of 6 degrees.
I have been using Octave to perform the integrations:
Code:
syms x E I
% Real structure
Mxr1 = (196.787*cosd(6) - 0.587*sind(6)).*x - 13.14.*x.^2/2 % left rafter
Mxr2 = (197.413*cosd(6) + 0.587*sind(6) - 1.174*sind(6)).*x + 9.392 - 13.14.*x.^2/2; % Right rafter
Mxc2 = -1.174*x % right-hand column
% Virtual structure
mxr1 = -sind(6)*x + 8 % Left rafter
mxr2 = sind(6)*x + 8 % Right rafter
mxc2 = -x % Right-hand column
% Integrations over excah member length for real and virtual structures
Ix1 = 1/(E*I)*int(Mxr1*mxr1,15.083,0)
Ix2 = 1/(E*I)*int(Mxr2*mxr2,15.083,0)
Ix3 = 1/(E*I)*int(Mxc2*mxc2,8,0)
ix1 = 1/(E*I)*int(mxr1.^2,15.083,0)
ix2 = 1/(E*I)*int(mxr2.^2,15.083,0)
ix3 = 1/(E*I)*int(mxc2.^2,8,0)
% calculate the horizontal reaction at A
HA = (Ix1+Ix2+Ix3)/(ix1+ix2+ix3)
For HA, from these calculations I am getting 111.90 kN. The software is giving me 98.6 kN. So, Im a bit off. I think I could be getting mixed up with signs, but I have checked and checked, and cannot see why I am getting the incorrect result.
Many Thanks!