Nodal displacements A New Jacobi-based Iterative Method for the Classical Analysis of Structures

Traditionally, classical methods of structural analysis such as slope-deflection and moment distribution methods (Cross method) are used for primary analysis of structures and also controlling the results of computer programs. The main objective of this paper is to introduce a new method for classical computing and extending it to a matrix formulation. The proposed approach, named the “Slope Distribution Method (SDM)”, is based on a Jacobi iterative procedure, in which without forming the system of linear equations, structural displacement values are obtained. Also, to make the method applicable and to use it in computer softwares, the matrix formulation of the approach is developed, where there is no need for iterative procedures and the nodal rotations are obtained through solving only one matrix equation. The SDM is able to analyze frames with non-vertical columns and those with nodal vertical displacement. Whereas current analysis softwares have some elimination for the analysis of non-prismatic members, the proposed method can be applied to analyze structures with any non-prismatic member. The SDM process is also developed for the analysis of dual lateral load resisting systems (moment resisting frames with other lateral load resisting elements such as bracings and shear walls). The advantages of the method over previous ones and also, its accuracy and reliability are presented through the article.


INTRODUCTION
Multi-storey building frames may be considered the most widely used kind of structures, especially in urban and residential areas.Population growth and land scarcity increase the need of these types of structures.Substantial and rapid expansion of this necessity in the early decades of the twentieth century, led to the creation of different methods of frames analysis.On the other hand, because of the high indetermination degree of frames, their analysis with traditional force methods was very time consuming.(However, some researchers have conducted a number of studies on the analysis time reduction in the force method, recently (Kaveh, 1992(Kaveh, , 2006))).Based on these two reasons, engineers tried to use more applicable and less time consuming methods such as slope-deflection, cross and Kani methods (Cross, 1930;Kani, 1957).Slope-deflection is a manual method to analyze the beams and bending frames and was introduced in 1915 by George A. Maney (Maney, 1915).In this method, with formation of equations and applying nodal and shear force equilibrium conditions, rotate angle of nodes and members are calculated which are placed in corresponding relations to determine end moments of members.Slope-deflection was a revolutionary approach in comparison to the previous structural analysis methods.Instead of calculating static redundancies and using them to find structural deformations, nodal displacements are computed firstly and used to achieve support reactions.
The fundamental problem of the Slope-deflection method appears in structures with a high degree of indetermination, which leads to the formation and solving of linear equation systems; so a lot of efforts were done to find methods that do not have this problem.This method was used for one decade in many cases, until the moment distribution method was developed.In 1930, Cross proposed the new method and could eliminate this problem for non-sway structures (Cross, 1930).In the first step of MDM, the rigid connections of bending frames were assumed to be fixed, and the moments created by external loads on these connections were obtained.These moments are unbalanced at the joints of the original non-restrained structure, and in order to equilibrate the joints, the moments are distributed proportionally to corresponding members' stiffness.The procedure repeats until the unbalanced moments become negligible.The final moments at the joints of members are the sum of all distributed incremental moments (Volokh, 2002).The Cross method in the sway structures needs forming and solving algebraic equations with fewer numbers of unknowns.The Cross approach has easy interpretation and has been taught in different universities.This method could be used in simple programming of structural analysis, in which end moment of members is considered as the unknown.By using an iterative procedure, (Consecutive moment distribution and transmission among the members connected to rigid nodes), the system of linear equations created by the slope -deflection method is solved.Consequently, beams and frames are analyzed without solving any system of equation, directly.End shear force for each member is also obtained through static equilibrium.
In addition, the Gaspar Kani's method (Kani, 1957) which is used in structural analysis, uses iterative procedures to solve the slope deflection equations system (Behravesh and Kaveh, 1990).In the moment distribution method, the unknowns (end moments of structural members) are obtained through performing iteration on their changes (increments), whereas in the Kani method, the iteration procedure is applied on the unknowns themselves.
Latin American Journal of Solids and Structures 12 (2015) 2581-2617 Behravesh and Kaveh(Behravesh and Kaveh, 1990) explained the relationship between moment distribution and Kani methods and a numerical iterative procedure, and showed that the calculation trends in these two methods are similar to the Jacobi iteration procedure that has been used to solve the equations of classical displacement.The Jacobi iteration approach in both methods is converged, if the stiffness matrix is diagonally dominant.A study of Volokh (Volokh, 2002)also shows the correlation of MDM to the Jacobi iteration approach.
In this paper, a new analysis approach is proposed which is named "Slope Distribution Method" (SDM).The method is explained in manual formulation and is extended in matrix form.In this approach, there is no need to form and solve the system of algebraic equations which is its advantage compared to the slope-deflection approach.In comparison to the moment distribution method (Cross method), in the proposed method, the distribution and carry-over procedures are merged, and unlike the Cross and Kani approaches, instead of distributing and transmitting the moments at several members' ends attached to each structural rigid node, only the nodal rotations (slopes) are distributed.These properties decrease the analysis parameters and lead to analysis time reduction of the proposed method.As the numbers of unknowns depends on the number of nodes and not the number of members connected to each node, the unknowns are limited in comparison to the Cross and Kani methods.According to the above-mentioned advantages, it will be shown that the proposed approach is less time-consuming than well-known methods of slope-deflection, Cross and Kani.The SDM process is also a proper and low-cost procedure to analyze the structures with non-prismatic members.It should be noted that the current analysis softwares are not able to model every kind of non-prismatic member; whereas by defining the corresponding coefficients, the SDM is capable of analyzing these structures properly.
The paper is organized as follows: The following section describes the Jacobi scheme.Afterwards, a brief review of the slope -deflection method is explained.Thereafter, the proposed slope distribution method (SDM) will be illustrated.The relation between the new method and the Jacobi scheme is shown, and the relations of the proposed method for non-prismatic members are expanded.In the next part, the matrix formulation of the method is extended.The approach is followed with some numerical examples and is compared with other traditional methods to show its advantages.The ability of the approach in analysis of dual systems, frames with non-vertical columns and also frames with vertical nodal displacements are shown.The accuracy and reliability of the approach are shown through the paper, and finally, the paper ends with some conclusions that are useful for researchers.

THE JACOBI SCHEME
The Jacobi scheme is a iterative method for solving the diagonally-dominant systems of equations (Golub and Van Loan, 1996;Young, 2013).In mathematics, a system of equations can be presented by a matrix format.If the magnitude of the diagonal entry in each row exceeds the sum of the magnitudes of the non-diagonal entries in that row, the matrix is called diagonally-dominant.In fact, matrix A is strictly diagonally dominant if (Behravesh and Kaveh, 1990;Datta, 2010).Solids and Structures 12 (2015) 2581-2617 Where, a ୧୨ denotes the entry in the ith row and jth column.For solving a set of equations in the matrix form Ax=B, the Jacobi relation is illustrated as (Golub and Van Loan, 1996;McCormac and Nelson, 1997;McGuire et al., 2000): In which, a ୧୨ and b ୧ denote the elements of coefficient matrix.

A BRIEF REVIEW OF SLOPE -DEFLECTION METHOD
As mentioned previously, in the suggested method, nodal displacements are considered the main unknowns.As nodal displacements are also the unknowns of the equations in the slope deflection method, this method will be reviewed briefly in this part.It should be noticed that unlike the slopedeflection method, the proposed method uses iterative procedures to solve the equations, and in general, it will not be necessary to solve the equations directly.The Slope-deflection method, as applied nowadays in analyzing the structures with rigid joints, was introduced by G.A.Maney in 1915.This method is based on evaluating the nodal rotations and displacements.After finding nodal rotations, end moment and shear of members can be calculated using the slopedeflection relations.These relations express end moment of a member based on rotations of nodes and elements, and also the external loads on the members.For element ij with constant bending rigidity, and length L (Figure 1), the relation is as follows (McCormac and Nelson, 1997;Norris et al., 1976).

‫ܯ‬
In these equations, the effects of axial and shear deformation are neglected, because the effect of these two parameters on the bending deformation is small.Using equations ( 3) and ( 4) and writing the equations of static equilibrium for the member ij, the end shear will be found as follows (Megson, 2005): Latin American Journal of Solids and Structures 12 (2015) 2581-2617 Ѳ ୧ and Ѳ ୨ are the rotation of node i and j respectively, and φ ୧୨ shows the sway rotation of member ij that is calculated through the below equation.FEM and FEV are the fixed end moment and shear. ) Δ:Relative displacement of the member ends.In Table 1, values of fixed end shear, moment and sign convention for different loading cases are shown (Mau, 2003): Table 1: Values of fixed end shear and moment (Mau, 2003).
Also, Equations (3) to ( 6) could be written in matrix formulation (Megson, 2005): By replacing the values of nodal displacements in the above matrix formulation, end shear force and moment of members are calculated.This matrix formation can be used in programming of the SDM matrix procedure.If a static equilibrium equation is written for a non-prismatic member, the shear force at the end of each member can be calculated.Given that the values of bending moments are defined according to nodal rotations, lateral rotation of members and also the external loads, the shears of members' ends are also obtained based on these three parameters: Latin American Journal of Solids and Structures 12 (2015) 2581-2617 where T ୧ and T ୨ are defined as follows: In the above equations, C ୧୨ and C ୨୧ are moment carry-over factors.Also, S ୧୨ and S ୨୧ are bending stiffness values of the two ends of the non-prismatic member and could be calculated by integration methods or using related handbooks (Association, 1958).

Introduction
SDM is an iterative method which is based on a series of computational cycles that are repeated until the results converge to a final value in each stage.By this repetition, simultaneous solving of algebraic equations is not necessary.This is a unique structural analysis method, in which the solution is obtained through an iterative process, without solving any equations to find problem parameters.SDM uses a Jacobi iterative scheme to find displacements in the system of equations produced by the slope-deflection method.

SDM equations
From static equilibrium, if external moment (M ୣ୶ ) is applied at the node B, the rotation at this node will take place until moment equilibrium is achieved, ∑ M = M ୣ୶ (Figure 2).‫ܯ‬ + ‫ܯ‬ େ = M ௫ (13) Using the equations from the Slope-Deflection Method given by Laursen (Laursen, 1988) and Kassimali (Kassimali, 1995), we have: L is the length of the member, EI is the cross-section rigidity, and B Ѳ is the rotation at B (Kassimali, 1995;Laursen, 1988).Placing Equations ( 14) and (15) into equation ( 13) and solving for B Ѳ , we have: In general, initial rotation of rigid node "i" has a direct relationship with the external load and the lateral rotation angle of members connected to the node.This dependency is shown in equation ( 17): In the mentioned equation, N is the number of elements connected to node i, is the relative bending stiffness of members connected to node i and K ୧୨ = ( ୍ ) ୧୨ is bending stiffness of each of these members.
In SDM, to calculate the slope of node i, the recursive series is defined as in equation ( 18): Where n is the number of analysis stage and θ ୧ () is the rotation value of node i under primary loads that are obtained through Equation ( 17).Δθ ୧ (க) is the difference between the connection slope at node i, in two successive stages, due to incremental unbalanced moments.The value of Δθ ୧ (க) will decrease during the analysis procedure and finally tends toward zero.In this condition, the connection rotation converges to its actual value, and the system will be balanced.
As discussed before, in the SDM method, incremental unbalanced moments at members' ends decrease to achieve nodes' equilibrium.In each step, for moment equilibrium condition at each node, it is necessary that: ∆M ୧ Repeating the slope distribution approach, equilibrium is established in connection when n→ ∞.
is the sum of incrementally unbalanced moments at node i and ∑ K ୧୨ ୨ୀ ୨ୀଵ is the sum of the flexural stiffness of members connected to the ith node.n is the number of the analysis cycle and N is the number of members connected to node i. Equation ( 14) is obtained through a simple static equilibrium on node i.The equation is the same as the Jacobi iterative procedure for solving the system of linear equations which has been described in previous parts.Using equation ( 19), the incremental unbalanced rotations caused by unbalanced moments, and consequently, the nodal rotations are calculated.
Inserting equations ( 3), ( 4) and ( 18) in equation ( 19), the fundamental equations of the SDM are consequently achieved: To compute the values of lateral rotation of the member ij,( φ ୧୨ (୬) ), a shear equilibrium equation should be used in each storey of the structure ( Figure 3).Through writing static equilibrium in each storey of the structure (i.e.Sth storey), the auxiliary equation is obtained to analyze structure without forming any equation system: In which, V ୗ is the shear of the Sth storey resulted from lateral loads such as earthquake or wind and V ୨ is the shear of the jth column of that storey and obtained through relations ( 5) and ( 6).The φ ୱ (୬) is calculated as follows: m is the number of columns in the Sth storey and i and j are the two ends of column ij; φ ୱ () is the initial rotation of member and D ୧୨ is the relative shear stiffness of the storey columns which is defined through equations ( 24) and ( 25): FEV is the fixed end shear of columns in the Sth storey; μ ୧୨ is shear stiffness of column ij, and D ୧୨ is the shear stiffness ratio of columns in the storey of interest.
The SDM method can be extended for the condition where the beams jointed to a node have a lateral displacement, and also when the φ values of structure members in a storey are not the same.As the lateral displacements of structural members (∆) in each storey are related by geometric relationships, the equations can be expressed based on the lateral rotation of a storey member (φ ୗ ) that is considered as the base and independent one.So, the above equations will be corrected as follows: Solids and Structures 12 (2015) 2581-2617   ‫ܦ‬ In the above equations, N is the number of elements connected to node i, m is the number of columns in the storey "S" and n is the step number of the recursive process.
As the lateral independent rotation of each storey of structure (φ ୱ (୬) ) is related to nodal rotations (Ѳ ୧ (୬) ), we can place equation ( 23) into equation ( 27) and simplify the analysis stages: For 2 ≤ n we have: Where, ቁ ω ୧୨ and R ୧୨ is the relative bending stiffness of members connected to node i which was introduced before.Also, the value of the carry-over factor (δ ୧୨ ) is defined as: The current equation performs as a recurrence relation.The rotation value of each node will be updated in each stage and finally converges to its accurate values.Closing the values of nodal rotations in two consecutive steps by desired accuracy is considered as the end of the SDM process.If a structure has non-prismatic members, then the above equations will be changed as follows: Latin American Journal of Solids and Structures 12 (2015) 2581-2617 For n≥ 2 we have: As the bending stiffness coefficients and carry-over factor are different at two ends of storey columns in non-prismatic structures, the shear stiffness ratio (D ୰ ) in node "r" of the column element should be calculated separately as follows: T ୰ is the shear stiffness of node "r" of the non-prismatic column which is obtained through equations ( 11) and ( 13), and ∑ T ୰ ୫ ଵ is the total shear stiffness of all nodes of columns (nodes of the two ends of the columns) in the storey of interest (i.e.storey "s") of the structure.Also, θ ୰ is the rotation of node "r" of the non-prismatic column.
When the non-prismatic structure has no lateral displacement in the storey "S"(φ ୱ = 0), by simplification, the below relations are conducted: For a special situation, prismatic structure without any lateral displacement in the storey "S" is replaced in the above equations.

Sign Convention
Sign convention in the proposed SDM method is similar to other manual analysis methods in which clockwise direction for θ,φ, FEM and M is considered as positive.According to this convention, the forces and linear displacements are taken as positive when the element rotates clockwise.It is in accordance with the positive directions used by some other researchers(e.g.(Timoshenko and Goodier, 1987).Also, the positive direction for shear force is the direction in which the force rotates the element clockwise (Figure 4).

MATRIX FORMULATION OF SDM
Manual approach of SDM should be repeated in a consecutive process until the final values of structure node displacements are obtained; this process is time-consuming and boring in structures with high degrees of indeterminacy.In fact, in matrix formulation of SDM, the initial information such as(φ () , θ () , δ and ω) will be calculated manually by the user and a computer will be used to perform the analysis process.Therefore, introducing an accurate solution without any repetition procedure and forming and solving only one matrix equation to give the unknowns, is the main purpose of the SDM matrix procedure.
If the equations ( 23), ( 25) and ( 26) are defined by matrix formulation, then we will have: If we have: And by inserting equation ( 46) in (47) we have: {θ} is a vector including the last rotations in structure nodes and has a rank of (p×1) and ൛θ () ൟ is also a column vector including the initial rotations of structure nodes due to external loads with a rank of (p×1); [ω] is a square matrix with a rank of (p×p) which includes the slope distribution factors of the members, and p is the number of rigid nodes and structure supports, and [Ι] is a unique matrix with an equal rank with [ω].The matrices of ሾω ሿ and ሾδሿ have the same rank with ሾωሿ and their entries are calculated just in two end nodes of the columns of each structure storey.They are replaced in corresponding matrixes to estimate the lateral displacement of structure (φ).
In non-prismatic structures, due to the difference between bending stiffness coefficients and carryover factors in two ends of storey columns, we have δ ୧୨ ≠ δ ୨୧ and they should be calculated separately for each column node.In a multi-storey structure, ሾω ሿand ሾδሿ matrices should be formed separately for each storey of the structure and the matrix ሾτሿ for all of the structure is calculated through the superposition principle by using the partial matrix of ሾτሿ in each storey of the structure.Therefore: In which, Ns is the number of stories in the structure.Each row in matrix ሾω ሿ is corresponded to the rigid node of structure storey, and the value of ሾω ሿ ୧ in node ί is equal to a total of (1 ) for the members connected to node ί on the storey "S"; which is placed in each entry of the ith row of matrix ሾω ሿ.Therefore: Indeed, the values of relative shear stiffness (D ୰ ), calculated through equation ( 38) are placed in matrix D, and matrix ሾδሿ is obtained through equation ( 50): It should be noted that during the process of placing the entries in the global matrix of the structure, the zero value is considered for entries corresponding to nodes which do not exist in the given storey.If the structure includes prismatic members, equations ( 49) and ( 50) are changed as follows: By calculating the above parameters, the rotations in structural nodes can be calculated.To this end, equation ( 27) can be changed to a matrix format.Equations ( 53) and (54) show the 2nd step of calculating ΔѲ: Latin American Journal of Solids and Structures 12 (2015) 2581-2617 In this way, through expanding relationship ( 27), other values of ൛ΔѲ (୬) ൟ can be obtained for n≥ 2: Then, through placing equations ( 44), ( 45) and ( 55) in ( 18), the values of ൛θ (୬) ൟ will be calculated in the last cycle.Therefore, it can be concluded that: If we define: ሾZሿ = InverseሼሾIሿ − ሾVሿሽ , ሾVሿ = ሾωሿ + ሾτሿ and ሾUሿ = ሾIሿ + ሾωሿ, then we will have: Final values of nodal rotations can be calculated through equation ( 56) by forming and solving only one matrix equation.End moments and shears of members can be evaluated through placing the nodal rotations in slope-deflection equations.The above equation can be used in computer programming of SDM.It should be noted that the slope distribution factors in matrix formulations of SDM are calculated for all members connected to the structure nodes, whereas the supports will also be considered as nodes.At first, to compute the bending stiffness, carry-over factors of members and fixed-end moments for non-prismatic members, the "hand book of frame constants"( (Association, 1958)) could be used.Then, the matrix of slope distribution factor, [ω], will be formed for a structure.This matrix is a square one with a rank of (p×p), where p is the number of structure nodes.The entry ω ୧୨ in matrix [ω] is equal to the slope distribution factor of member ij, which connects the nodes i and j.For those nodes which are not connected to each other by any member, ω ୧୨ = 0; for original diameter entries, ω ୧୧ = 0 and for fixed and two-roller supports in which R ୧୨ = 0, ω ୧୨ is also considered as zero.For hinge and roller supports, ω ୧୨ = − C ୧୨ , because the value of R ୧୨ is equal to unity (see equation ( 40)).The value of initial rotation for a rigid node of structure under external loading is obtained by equations ( 26) and (32).For fixed and two-roller supports, θ ୧ () = 0 and for hinged support and roller support, it is computed through equation ( 26).
It is worth noting that the initial rotation value in matrix formulation of SDM is computed for a hinged support; So, the bending stiffness, carry-over factor and the fixed-end moments obtained through "handbooks of frame constants" ( (Association, 1958)) will be used in the equations without any changes.If the structure has no lateral displacement (φ ୧୨ = 0), equation ( 56) will be simplified as follows:

NUMERICAL EXAMPLES
6.1 In the first example, the proposed procedure of SDM is applied on the continuous beam of A and D are fixed.Other specifications are shown in Figure 5.The cross section of beam is Ushape (UPN 80).EI is assumed to be constant through the length of the beam.Table 3 shows the steps of Jacobi-based SDM for analysis of the continuous beam of example 6.1:  Below, the parametric calculation of some steps of the SDM is shown: Using slope deflection relations, the values of end moments and shears for the beam are calculated as below: To apply the matrix formulation on the above example, [ω] and ሼθ ሽ are formed for the beam of Figure 5, by using the parameter values in Table 2: By applying equation (58) and using the calculated [ω] andሼθ ሽ, nodal rotation matrix is obtained as follow: As could be seen, the results of the manual and matrix formulation are in good agreement.The above beam is also analyzed by Sap2000 software.The results are shown in Figures.6 and 7.As can be seen, the difference between the results of the SDM with those of these diagrams is negligible and is less than 1% .It seems this small difference is because of considering shear deformation in software calculations.6.2In this example, bending frame with lateral displacement and different columns heights, will analyze by classic and matrix procedure of SDM and results will be compared with those of Sap2000 software.Characteristics, loading and support conditions are shown in Figure 8 and the flexural rigidity of frame members is constant.If it is supposed that φ is considered as the independent φ of the storey, and we have: φ ୱ =φ , then regarding geometric relationships between the members of the frame, it can be obtained that: At first, primarily needed parameters are calculated.These values are presented in Table 4.
Now, slope distribution relationships between rigid nodes are simply produced through computing needed parameters.The parametric form of the slope distribution procedure for calculating nodal rotations of Figure 8 is shown in Table 5 (for the first three steps).Similarly, by repeating this procedure, other step values of nodal rotations' changes are calculated.The results of all stages are shown in Figure 9.
߱ ∆ߠ (ଵ) + ߱ େ ቂߜ ∆ߠ (ଵ) + ߜ ∆ߠ (ଵ) ቃ ߱ ∆ߠ (ଵ) + ߱ ቂߜ ∆ߠ (ଵ) + ߜ ∆ߠ (ଵ) ቃ ∆ߠ (ଶ)   Table 5: Parametric form of SDM procedure for calculating nodal rotation of example 6.2.Now, for the matrix procedure of SDM, by calculating the required parameters, the corresponding matrices can be formed.Then, through SDM matrix procedure, frame node rotations will be defined: As can be observed, the calculated nodes rotation values through manual procedure of SDM after the sixth stages are equal to the results from their matrix procedure, till two decimals value; this shows the proper accuracy and convergence of the proposed method.Now, the lateral rotation of members, end moment and end shear values of members can be computed using equation ( 23) and slope deflection relations: To compare the results obtained through the proposed method with those of computer programs; the above frame is modeled and analyzed through Sap2000 software.In Figure 10, the shear and moment diagrams of frame members are shown as the two outputs of software:  In such a problem, the stiffness of braces in the storey should be included in calculations.What analyzing relations reveals is, while designing of the two given bracing members is performed based on allowable compression force, both members will participate in the lateral resisting system; so, the total lateral stiffness of the braces is equal to total lateral stiffness of the two members (Rezaee Pajand and Aftabi Suny, 2010; Zalka, 2002): Where, E is elasticity modulus, A is cross-section area, L, is the length of bracing member and Ө is the angle between the bracing member axis and horizontal direction.Regarding that the above system is a dual bending-bracing frame system, it is needed to consider the stiffness of these two systems during the computations processes, simultaneously.According to Figure 3, the storey shear (V ୱ ) is equal to the sum of the column's shear forces and the horizontal components of the bracing members forces in the desired storey: In which, m is the number of column members in the storey and V ୴ is obtained by the below equation: In this equation, the sum of lateral stiffness of the bracing system in the Sth storey is defined by K ୴ .
According to Figure 4 and based on the main assumption of the SDM method, i.e. neglecting axial deformation of frame members, relative displacements of all Sth storey columns are the same and equal to ∆ ୱ .
All bending frame equations are used in dual frames, but member shear stiffness will be changed.Therefore equations ( 28) and ( 29) are modified as follows: If the structural members are prismatic, the above equations will be simplified as below: Now, to solve the above problem, stiffness and carry-over coefficients of non-prismatic columns are computed at first, using the tables of "handbook of frame constants" (Association, 1958), as follows: Latin American Journal of Solids and Structures 12 (2015) 2581-2617 The frame beam is prismatic and has stiffness coefficient and carry-over factors as below: Now, by defining the required data, the rotation coefficients are calculated:    Through comparing manual and matrix procedures of SDM, it can be shown that the value of frame nodal rotation after the sixth step using a manual procedure is equal to its matrix results till two decimal values which indicate the accuracy and good convergence speed of the proposed method.
After calculating the nodal displacements, members' lateral rotations and bending moments are resulted by equation ( 34 In the above example, if the bracing system is eccentric, the only difference is the value of the stiffness of bracing members.Regarding studies (Rezaee Pajand and Aftabi Suny, 2010;Zalka, 2002), it is just enough that parameter "A" in equation ( 59) is replaced by "A ୣ " which is defined as follows: In the above equation, β is the proportion of height to length of braced span and α is the ratio of distance between the brace connection place to its neighbor column to the length of braced span.
Latin American Journal of Solids and Structures 12 (2015) 2581-2617 These two factors can be seen in Figure 13.In a situation where shear walls are used in the structure, it is just enough to define wall lateral stiffness through the following equation to analyze the structure by the SDM method: Where, E is elasticity modulus for concrete shear wall and I is the moment inertia of shear wall and are calculated by following equation: f ୡ is compressive strength of concrete based on N/mm2 .b and L ୵ are the thickness and the length of the wall, respectively.Coefficient γ in equation ( 67) shows the effect of shear deformation on the wall stiffness and obtained by the following equation: where, H is the height of the wall.
7 BENDING FRAME WITH NON-VERTICAL COLUMNS.

7.1
The SDM method can be used for frames at which the beams at storey level have lateral rotation or the columns are non-vertical.Through writing horizontal equilibrium equations for shear force, the horizontal components of the axial forces of non-vertical columns are participated in the equilibrium.Therefore, to ease the process, moment equilibrium equation around a virtual point resulted from the intersection of members or the continuation (point I in Figure 14), are used to remove the effect of column axial force in equations.If the concentrated load of P is applied on the connection "B" in the frame of Figure 14, the parameters of the SDM approach can be defined for the frame through writing moment equilibrium for point "I". Figure 14: An example of a frame with non-vertical columns and Lateral displacements of its members (Kaveh, 2012;Megson, 2005).
If φ ୱ =φ , then from equation ( 72), we have the following relations: In this example, the required parameters are calculated firstly and the structure is analyzed by the two manual and matrix forms of the proposed method.
Now, through calculating required parameters, the slope distribution process will be repeated between rigid connections.Parametric form of slope distribution procedure for calculating nodal rotation of the frame of Figure 15 is similar to Table 5.In  For analysis of the frame by matrix formulation of the presented approach, the corresponding matrices should be formed at first: By forming the initial rotation matrix and using equation ( 56), the final rotation value of the nodes are calculated: By using equation ( 23) and applying slope-deflection equations, the end-moments and shears of the members are calculated.
The comparison of the results with those of (Kaveh, 2012), shows the accuracy of the proposed method.

8.1
If a frame column does not continue to the foundation level, due to vertical displacements, the frame beams will have lateral rotation.This is a limitation of the Kani method which is not able to analyze a frame with vertical displacement, and structure columns should be continued to the foundation level.But through the SDM method, these kinds of structures can be analyzed.An example of this structure is shown in Figure 16.load and horizontal force resulted from an earthquake.The frame geometrical characteristics and loading are represented in Figure 18.Finally, the obtained results will be compared with those of SAP2000 software.Bending rigidity values of all frame members are constant and equal.Considering the geometrical relation between structural members, we have: The required coefficients and factors are calculated through simple relations of structural analysis (Tables 10 and 11).For analysis of the frame by matrix procedure, the corresponding matrices are formed firstly based on Tables 10 and 11. process are merged together.This characteristic and considering the fact that the SDM is repeated only on rigid nodes and not on the moment values of all the members connected to rigid nodes, analysis time consuming and steps are reduced, comparatively.Regarding the later characteristics, analysis time consuming and the number of stages is also less than the Kani method; and by obtaining nodal rotation and lateral rotation of the storey, end moment and shear values of members can be calculated, simultaneously.In another work, the steps of the manual formulation were shown by a geometric progression, which led to matrix formulation.The main advantage of the matrix formulation is using only a matrix equation for the unknowns that could be used in computer programming.To show the accuracy of the Matrix formulation, some examples were solved.The results also show the accuracy and effectiveness of the proposed matrix formulation of SDM.
Another advantage of the proposed method is that the approach does not contain limitations as those in the Kani method and is capable of analyzing frames with non-vertical columns and also, those with vertical displacement.In the proposed method, with some modification, lateral stiffness of bracing members could be applied in SDM equations, so the proposed method is able to analyze bracing and dual systems.The SDM process is also capable of analyzing structures containing nonprismatic members.It should be noted that the current analysis software are not usually able to model every kind of non-prismatic member; while by specifying some basic parameters, the SDM approach is able to simply analyze structures with these kinds of members.The results of this study could be used in structural engineering calculations and the method could be expanded for other specific structures.

Figure 1 :
Figure 1: Deformation curve of the elastic beam.

Figure 3 :
Figure 3: The horizontal Shear of the frame in the Sth Storey (Rezaee Pajand and Aftabi Suny, 2010).

Figure 4 :
Figure 4: Positive sign convention for force and displacement.
from reference(Lopes et al., 2011).This beam has 4 supports in A, B, C and D. The supports Latin American Journal ofSolids and Structures 12 (2015) 2581-2617

Figure 6 :Figure 7 :
Figure 6: Shear diagram (KN) for the beam of the example 6.1.
Fixed-end moment of members, slope distribution and carry-over factors for the frame of Figure8.Latin American Journal ofSolids and Structures 12 (2015) 2581-2617

Figure 9 :
Figure9: SDM procedure for analysis of the moment frame of example 6.2.

Figure 10 :
Figure 10: The shear diagram (in Ton) and moment diagram (in Ton-m) of the frame of example 6.2

Figure 11 :
Figure 11: The one bay moment resisting frame with non-prismatic columns and x-bracing system(Rezaee Pajand and Aftabi Suny, 2010)
moment of the members, slope distribution and carry-over factors of the frame of Figure15.

Figure 16 :
Figure 16: A moment resisting frame with vertical displacements of the nodes.

Figure 18 :
Figure 18: The bending frame of Figure 17.

Figure 19 :
Figure 19: Classical procedure of SDM for analysis of the frame of example 2.8.

Table 3
The steps of Jacobi-based SDM for analysis of the continuous beam of example 6.1.

Table 6 :
Parametric form of SDM procedure for nodal rotation of example 6.3.

Table 7 :
The results of SDM procedure for nodal rotation of example 6.3.
Table 9, the calculation process of the rotations for nodes B and C are presented:

Table 9 :
Calculation process of the rotation values for nodes B and C of the frame of Figure 15.

Table 10 :
Fixed-end moment of members, slop distribution and carry-over factors for example 8.2.