Flexibility method

﻿
Flexibility method

In structural engineering, the flexibility method is the "classical" consistent deformation method for computing member forces and displacements in structural systems. Its modern version formulated in terms of the members' flexibility matrices also has the name the matrix force method due to its use of member forces as the primary unknowns.

Member flexibility

Flexibility is the inverse of stiffness. For example, consider a spring that has "Q" and "q" as, respectively, its force and deformation:
* The spring stiffness relation is "Q = k q" where "k" is the spring stiffness.
* Its flexibility relation is "q = f Q", where "f" is the spring flexibility.
* Hence, "f" = 1/"k".

A typical member flexibility relation has the following general form:

:$mathbf\left\{q\right\}^m = mathbf\left\{f\right\}^m mathbf\left\{Q\right\}^m + mathbf\left\{q\right\}^\left\{om\right\} qquad qquad qquad mathrm\left\{\left(1\right)\right\}$where:"m" = member number "m".:$mathbf\left\{q\right\}^m$ = vector of member's characteristic deformations.:$mathbf\left\{f\right\}^m$ = member flexibility matrix which characterises the member's susceptibility to deform under forces.:$mathbf\left\{Q\right\}^m$ = vector of member's independent characteristic forces, which are unknown internal forces. These independent forces give rise to all member-end forces by member equilibrium.:$mathbf\left\{q\right\}^\left\{om\right\}$ = vector of member's characteristic deformations caused by external effects (such as known forces and temperature changes) applied to the isolated, disconnected member (i.e. with $mathbf\left\{Q\right\}^m = 0$).

For a system composed of many members interconnected at points called nodes, the members' flexibility relations can be put together into a single matrix equation, dropping the superscript m:

:$mathbf\left\{q\right\}_\left\{M imes 1\right\} = mathbf\left\{f\right\}_\left\{M imes M\right\} mathbf\left\{Q\right\}_\left\{M imes 1\right\} + mathbf\left\{q\right\}^\left\{o\right\}_\left\{M imes 1\right\} qquad qquad qquad mathrm\left\{\left(2\right)\right\}$

where "M" is the total number of members' characteristic deformations or forces in the system.

Unlike the matrix stiffness method, where the members' stiffness relations can be readily integrated via nodal equilibrium and compatibility conditions, the present flexibility form of equation (2) poses serious difficulty. With member forces $mathbf\left\{Q\right\}_\left\{M imes 1\right\}$ as the primary unknowns, the number of nodal equilibrium equations is insufficient for solution, in general--unless the system is statically determinate.

Nodal equilibrium equations

To resolve this difficulty, first we make use of the nodal equilibrium equations in order to reduce the number of independent unknown member forces. The nodal equilibrium equation for the system has the form:

:$mathbf\left\{R\right\}_\left\{N imes 1\right\} = mathbf\left\{b\right\}_\left\{N imes M\right\} mathbf\left\{Q\right\}_\left\{M imes 1\right\} + mathbf\left\{W\right\}_\left\{N imes 1\right\} qquad qquad qquad mathrm\left\{\left(3\right)\right\}$

where: $mathbf\left\{R\right\}_\left\{N imes 1\right\}$: Vector of nodal forces at all N degrees of freedom of the system.: $mathbf\left\{b\right\}_\left\{N imes M\right\}$: The resulting nodal equilibrium matrix: $mathbf\left\{W\right\}_\left\{N imes 1\right\}$: The vector of forces arising from loading on the members.

In the case of determinate systems, matrix b is square and the solution for Q can be found immediately from (3) provided that the system is stable.

The primary system

For statically indeterminate systems, "M > N", and hence, we can augment (3) with "I = M-N" equations of the form:

:

The vector X is the so-called vector of redundant forces and "I" is the degree of statical indeterminacy of the system. We usually choose "j", "k", ... , $alpha$, and such that $X_i$ is a support reaction or an internal member-end force. With suitable choices of redundant forces, the equation system (3) augmented by (4) can now be solved to obtain:

:$mathbf\left\{Q\right\}_\left\{M imes 1\right\} = mathbf\left\{B\right\}_R mathbf\left\{R\right\}_\left\{N imes 1\right\} + mathbf\left\{B\right\}_X mathbf\left\{X\right\}_\left\{I imes 1\right\} + mathbf\left\{Q\right\}_\left\{v cdot M imes 1\right\} qquad qquad qquad mathrm\left\{\left(5\right)\right\}$

Substitution into (2) gives:

:$mathbf\left\{q\right\}_\left\{M imes 1\right\} = mathbf\left\{f\right\}_\left\{M imes M\right\}$

Big( mathbf{B}_R mathbf{R}_{N imes 1} + mathbf{B}_X mathbf{X}_{I imes 1} + mathbf{Q}_{v cdot M imes 1} Big)

+ mathbf{q}^{o}_{M imes 1} qquad qquad qquad mathrm{(6)}

Equations (5) and (6) are the solution for the "primary system" which is the original system that has been rendered statically determinate by cuts that expose the redundant forces $mathbf\left\{X\right\}$. Equation (5) effectively reduces the set of unknown forces to $mathbf\left\{X\right\}$.

Compatibility equation and solution

Next, we need to set up $I$ compatibility equations in order to find $mathbf\left\{X\right\}$. The compatibility equations restore the required continuity at the cut sections by setting the relative displacements $mathbf\left\{r\right\}_\left\{X\right\}$ at the redundants X to zero. That is, using the unit dummy force method:

:$mathbf\left\{r\right\}_\left\{X\right\} = mathbf\left\{B\right\}_X^T mathbf\left\{q\right\} = mathbf\left\{B\right\}_X^T Big \left[ mathbf\left\{f\right\} Big\left( mathbf\left\{B\right\}_R mathbf\left\{R\right\} + mathbf\left\{B\right\}_X mathbf\left\{X\right\} + mathbf\left\{Q\right\}_v Big\right) + mathbf\left\{q\right\}^\left\{o\right\} Big\right] = 0 qquad qquad qquad mathrm\left\{\left(7a\right)\right\}$

:or $mathbf\left\{r\right\}_\left\{X\right\} = mathbf\left\{F\right\}_\left\{XX\right\} mathbf\left\{X\right\} + mathbf\left\{r\right\}^o_X = 0 qquad qquad qquad mathrm\left\{\left(7b\right)\right\}$where: $mathbf\left\{F\right\}_\left\{XX\right\} = mathbf\left\{B\right\}_X^T mathbf\left\{f\right\} mathbf\left\{B\right\}_X$

:$mathbf\left\{r\right\}^o_X = mathbf\left\{B\right\}_X^T Big \left[ mathbf\left\{f\right\} Big\left( mathbf\left\{B\right\}_R mathbf\left\{R\right\} + mathbf\left\{Q\right\}_v Big\right) + mathbf\left\{q\right\}^\left\{o\right\} Big\right]$ Equation (7b) can be solved for X, and the member forces are next found from (5) while the nodal displacements can be found by

:$mathbf\left\{r\right\}_\left\{R\right\} = mathbf\left\{B\right\}_R^T mathbf\left\{q\right\} = mathbf\left\{F\right\}_\left\{RR\right\} mathbf\left\{R\right\} + mathbf\left\{r\right\}^o_R$

where: $mathbf\left\{F\right\}_\left\{RR\right\} = mathbf\left\{B\right\}_R^T mathbf\left\{f\right\} mathbf\left\{B\right\}_R$ is the "system flexibility matrix".

:$mathbf\left\{r\right\}^o_R = mathbf\left\{B\right\}_R^T Big \left[ mathbf\left\{f\right\} Big\left( mathbf\left\{B\right\}_X mathbf\left\{X\right\} + mathbf\left\{Q\right\}_v Big\right) + mathbf\left\{q\right\}^\left\{o\right\} Big\right]$

Supports' movements taking place at the redundants can be included in the right-hand-side of equation (7), while supports' movements at other places must be included in $mathbf\left\{r\right\}^o_X$ and $mathbf\left\{r\right\}^o_R$ as well. For educational purposes, the above computations are best done using the computer with a suitable software package such as [http://www.mathworks.com/ MATLAB] or the freeware [http://groups.yahoo.com/group/cmap/ CMAP]

Advantages and disadvantages

While the choice of redundants in (4) appears to be arbitrary and troublesome for automatic computation, this objection can be overcome by proceeding from (3) directly to (5) using a modified Gauss-Jordan elimination process. This is a robust procedure that automatically selects a good set of redundant forces to ensure numerical stability.

It is apparent from the above process that the matrix stiffness method is easier to comprehend and to implement for automatic computation. It is also easier to extend for advanced applications such as non-linear analysis, stability, vibrations, etc. For these reasons, the matrix stiffness method is the method of choice for use in general purpose structural analysis software packages. On the other hand, for linear systems with low degree of statical indeterminacy, the flexibility method has the advantage of being computationally less intensive. This advantage, however, is a moot point as personal computers are widely available and more powerful. The main redeeming factor in learning this method nowadays is its educational value in imparting the concepts of equilibrium and compatibility in addition to its historical value. In contrast, the procedure of the direct stiffness method is so mechanical that it risks being used without much understanding of the structural behaviors.

The upper arguments were valid up to the late 1990s. However, recent advances in numerical computing have shown a come back of the force method, especially in the case of nonlinear systems. New frameworks have been developed that allow "exact" formulations irrespectively of the type or nature of the system nonlinearities. The main advantages of the flexibility method is that the result error is independent of the discretization of the model and that it is indeed a very fast method. For instance, the elastic-plastic solution of a continuous beam using the force method requires only 4 beam elements whereas a commercial "stiffness based" FEM code requires 500 elements in order to give results with the same accuracy. To conclude, one can say that in the case where the solution of the problem requires recursive evaluations of the force field like in the case of structural optimization or system identification, the efficiency of the flexibility method is indisputable.

ee also

*Finite element method in structural mechanics
*Structural analysis

External links

[http://www.public.iastate.edu/~fanous/ce332/force/homepage.html Consistent Deformations - Force Method]

Wikimedia Foundation. 2010.

Look at other dictionaries:

• Method of conditional probabilities — In mathematics and computer science, the probabilistic method is used to prove the existence of mathematical objects with desired combinatorial properties. The proofs are probabilistic they work by showing that a random object, chosen from some… …   Wikipedia

• Direct stiffness method — As one of the methods of structural analysis, the direct stiffness method (DSM), also known as the displacement method or matrix stiffness method, is particularly suited for computer automated analysis of complex structures including the… …   Wikipedia

• Matrix stiffness method — In structural engineering, the matrix stiffness method (or simply stiffness method , also known as Direct stiffness method) is a matrix method that makes use of the members stiffness relations for computing member forces and displacements in… …   Wikipedia

• Matrix method — The matrix method is a structural analysis method used as a fundamental principle in many applications in civil engineering. The method is carried out, using either a stiffness matrix or a flexibility matrix. The flexibility method is not… …   Wikipedia

• Finite element method in structural mechanics — Finite element method (FEM) is a powerful technique originally developed for numerical solution of complex problems in structural mechanics, and it remains the method of choice for complex systems. In the FEM, the structural system is modeled by… …   Wikipedia

• Labour market flexibility — refers to the speed with which labour markets adapt to fluctuations and changes in society, the economy or production.DefinitionIn the past, the most common definition of labour market flexibility was the neo liberal definition. This entailed the …   Wikipedia

• Discrete element method — A discrete element method (DEM), also called a distinct element method is any of family of numerical methods for computing the motion of a large number of particles of micrometre scale size and above. Though DEM is very closely related to… …   Wikipedia

• Regulatory Flexibility Act — Infobox U.S. legislation name= Regulatory Flexibility Act fullname= The Regulatory Flexibility Act acronym= RFA / the Reg Flex Act enacted by= effective date= 1 January 1981 public law url= http://www.archives.gov/federal register/laws/regulatory …   Wikipedia

• Critical path method — PERT chart for a project with five milestones (10 through 50) and six activities (A through F). The project has two critical paths: activities B and C, or A, D, and F – giving a minimum project time of 7 months with fast tracking. Activity E is… …   Wikipedia

• Product software implementation method — A product software implementation method is a systematically structured approach to effectively integrate a software based service or component into the workflow of an organizational structure or an individual end user.This entry focuses on the… …   Wikipedia