Engineering Have I solved this structural engineering equation correctly?

  • Thread starter Thread starter Tygra
  • Start date Start date
AI Thread Summary
The discussion focuses on a structural engineering equation for calculating rotations at rigid joints using an iterative procedure. The user has provided a specific formula and MATLAB code to compute these rotations based on given parameters such as horizontal storey shear and storey height. The initial values for rotations are derived from a formula in a 1979 engineering book, and the user seeks validation of their calculations and code. They emphasize that the last rotation value is set to zero due to a fixed support condition. Assistance is requested to confirm the correctness of the approach and calculations.
Tygra
Messages
55
Reaction score
8
Homework Statement
A structural engineering iterative equations to find rotations at joints
Relevant Equations
$$ x_i = \frac {Q_ih_i + Q_{i+1}h_{i+1}}{4K} + \frac {C}{K}x_{i-1} + \frac {C}{K}x_{i+1} $$
Hi all,

I have a structural engineering book from 1979. I am trying to follow it as best as I can. I have come to a formula that calculates the rotations in radians at the rigid joint that requires an iterative procedure. This equation comes in the form of:

$$ x_i = \frac {Q_ih_i + Q_{i+1}h_{i+1}}{4K} + \frac {C}{K}x_{i-1} + \frac {C}{K}x_{i+1} $$

Where:

## Q ## is the horizontal storey shear

## h ## is the storey height

## K = (6G_i + C_i + C_{i+1}) ##

## G = \frac {I_g}{h} ##

## C = \frac {I_c}{L} ##

## L ## is the girder length

## x ## is the rotations for the ## ith ## storey

The book says to use initial values x by using this formula:

$$ x_i = \frac {Q_ih_i + Q_{i+1}h_{i+1}}{24G_i} $$

After you can compute the rotations for each level use the iterative formula to improve the values of the rotations.

Now, I have used MATALB to do this. Here is my code:

Matlab:
x1(1) = (Q*h/4)/(24*G1)
x2(1) = (2*Q*h/4 + Q*h/4)/(24*G1)
x3(1) = (3*Q*h/4 + 2*Q*h/4)/(24*G1)
x4(1) = (4*Q*h/4 + 3*Q*h/4)/(24*G1)
x5(1) = (5*Q*h/4 + 4*Q*h/4)/(24*G1)
x6(1) = (6*Q*h/4 + 5*Q*h/4)/(24*G1)
x7(1) = (7*Q*h/4 + 6*Q*h/4)/(24*G1 + 2*C)
x8(1) = 0

for i = 1:10
     x1(i+1) = Q*h/(4*K) + C/K*x2(i);
     x2(i+1) = (Q*h + 2*Q*h)/(4*K) + C/K*x1(i) + C/K*x3(i);
     x3(i+1) = (2*Q*h + 3*Q*h)/(4*K) + C/K*x2(i) + C/K*x4(i);
     x4(i+1) = (3*Q*h + 4*Q*h)/(4*K) + C/K*x3(i) + C/K*x5(i);
     x5(i+1) = (4*Q*h + 5*Q*h)/(4*K) + C/K*x4(i) + C/K*x6(i);
     x6(i+1) = (5*Q*h + 6*Q*h)/(4*K) + C/K*x5(i) + C/K*x7(i);
     x7(i+1) = (6*Q*h + 7*Q*h)/(4*K) + C/K*x6(i);
     x8(i+1) = 0;
end

So from this is get my rotations at each storey. I haven't posted my results. I merely hope for someone to check the above workings and let me know if I have done this correctly.

For a visual aid, here is a typical structure of seven storeys with 8 rotations. Note, x8 equals zero because it is clamped against rotation because I am using a fixed support.
7 Storey Rigid Frame.webp


I am hoping someone can help?

Many thanks
 
Physics news on Phys.org
Looking at your code, I think your Q and h variables need to have indexed values...Q(1) h(1), Q(i+1), etc. Otherwise, it looks ok.

It seems that you have loaded it backwards in the code or the image. As you have written it, a +1 height in your code has it going a higher floor instead of lower floor like you have in the image.
 
Thread 'PLL - How to find all the gains of a PI corrector and fix Ki ? MATLAB'
Hello everyone, it's not really for my homework (because I'm not at school anymore) but I'm new and I don't know how to start a new forum. I've never done automatic before and I need help. I have a 2nd ordre PLL on MATLAB, with an NCO and a PI corrector; the PLL works via the Mid/End algorithm (that's part is ok). I need to find the link between the theoretical formulas and their values in the code to find the value of the gain Ki, and possibly that of Kp if it's not correct. Unfortunately...
Thread 'Have I solved this structural engineering equation correctly?'
Hi all, I have a structural engineering book from 1979. I am trying to follow it as best as I can. I have come to a formula that calculates the rotations in radians at the rigid joint that requires an iterative procedure. This equation comes in the form of: $$ x_i = \frac {Q_ih_i + Q_{i+1}h_{i+1}}{4K} + \frac {C}{K}x_{i-1} + \frac {C}{K}x_{i+1} $$ Where: ## Q ## is the horizontal storey shear ## h ## is the storey height ## K = (6G_i + C_i + C_{i+1}) ## ## G = \frac {I_g}{h} ## ## C...

Similar threads

Replies
2
Views
323
Replies
0
Views
1K
Replies
5
Views
3K
Replies
2
Views
2K
Replies
0
Views
1K
Replies
1
Views
1K
Back
Top