text:
Finite element method in structural mechanics
From Wikipedia, the free encyclopedia
Jump to: navigation, search
This article needs attention from an expert on the subject. Please add a reason or a talk parameter to this template to explain the issue with the article. Consider associating this request with a WikiProject. (October 2010)
The 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 a set of appropriate finite elements interconnected at discrete points called nodes. Elements may have physical properties such as thickness, coefficient of thermal expansion, density, Young's modulus, shear modulus and Poisson's ratio.
Contents
[hide]
1 History
2 Element properties
2.1 Element interconnection and displacement
3 Practical considerations
4 Theoretical overview of FEM-Displacement Formulation: From elements, to system, to solution
5 Interpolation or shape functions
6 Internal virtual work in a typical element
6.1 Element matrices
6.2 Element virtual work in terms of system nodal displacements
7 System virtual work
8 Assembly of system matrices
9 See also
10 References
History[edit]
The origin of finite method can be traced to the matrix analysis of structures [1][2] where the concept of a displacement or stiffness matrix approach was introduced. Finite element concepts were developed based on engineering methods in 1950s. The finite element method obtained its real impetus in the 1960s and 1970s by John Argyris, and co-workers; at the University of Stuttgart, by Ray W. Clough; at the University of California, Berkeley, by Olgierd Zienkiewicz, and co-workers Ernest Hinton, Bruce Irons;[3] at the University of Swansea, by Philippe G. Ciarlet; at the University of Paris; at Cornell University, by Richard Gallagher and co-workers. The original works such as those by Argyris [4] and Clough [5] became the foundation for today’s finite element structural analysis methods. Earlier books such as by Zienkiewicz [6] and more recent books such as by Yang [7] give comprehensive summary of developments in finite-element structural analysis. Implementing the method in software is described in the classic text by Smith, Griffiths and Margetts.[8]
Element properties[edit]
Straight or curved one-dimensional elements with physical properties such as axial, bending, and torsional stiffnesses. This type of element is suitable for modeling cables, braces, trusses, beams, stiffeners, grids and frames. Straight elements usually have two nodes, one at each end, while curved elements will need at least three nodes including the end-nodes. The elements are positioned at the centroidal axis of the actual members.
Two-dimensional elements that resist only in-plane forces by membrane action (plane stress, plane strain), and plates that resist transverse loads by transverse shear and bending action (plates and shells). They may have a variety of shapes such as flat or curved triangles and quadrilaterals. Nodes are usually placed at the element corners, and if needed for higher accuracy, additional nodes can be placed along the element edges or even within the element. The elements are positioned at the mid-surface of the actual layer thickness.
Torus-shaped elements for axisymmetric problems such as membranes, thick plates, shells, and solids. The cross-section of the elements are similar to the previously described types: one-dimensional for thin plates and shells, and two-dimensional for solids, thick plates and shells.
Three-dimensional elements for modeling 3-D solids such as machine components, dams, embankments or soil masses. Common element shapes include tetrahedrals and hexahedrals. Nodes are placed at the vertexes and possibly in the element faces or within the element.
Element interconnection and displacement[edit]
The elements are interconnected only at the exterior nodes, and altogether they should cover the entire domain as accurately as possible. Nodes will have nodal (vector) displacements or degrees of freedom which may include translations, rotations, and for special applications, higher order derivatives of displacements. When the nodes displace, they will drag the elements along in a certain manner dictated by the element formulation. In other words, displacements of any points in the element will be interpolated from the nodal displacements, and this is the main reason for the approximate nature of the solution.
Practical considerations[edit]
From the application point of view, it is important to model the system such that:
Symmetry or anti-symmetry conditions are exploited in order to reduce the size of the model.
Displacement compatibility, including any required discontinuity, is ensured at the nodes, and preferably, along the element edges as well, particularly when adjacent elements are of different types, material or thickness. Compatibility of displacements of many nodes can usually be imposed via constraint relations.
Elements' behaviours must capture the dominant actions of the actual system, both locally and globally.
The element mesh should be sufficiently fine in order to produce acceptable accuracy. To assess accuracy, the mesh is refined until the important results shows little change. For higher accuracy, the aspect ratio of the elements should be as close to unity as possible, and smaller elements are used over the parts of higher stress gradient.
Proper support constraints are imposed with special attention paid to nodes on symmetry axes.
Large scale commercial software packages often provide facilities for generating the mesh, and the graphical display of input and output, which greatly facilitate the verification of both input data and interpretation of the results.
Theoretical overview of FEM-Displacement Formulation: From elements, to system, to solution[edit]
While the theory of FEM can be presented in different perspectives or emphases, its development for structural analysis follows the more traditional approach via the virtual work principle or the minimum total potential energy principle. The virtual work principle approach is more general as it is applicable to both linear and non-linear material behaviours. The virtual work method is an expression of conservation of energy: for conservative systems, the work added to the system by a set of applied forces is equal to the energy stored in the system in the form of strain energy of the structure's components.
The principle of virtual displacements for the structural system expresses the mathematical identity of external and internal virtual work:
External virtual work = ∫ V δ ϵ T σ d V ( 1 ) {\displaystyle {\mbox{External virtual work}}=\int _{V}\delta {\boldsymbol {\epsilon }}^{T}{\boldsymbol {\sigma }}\,dV\qquad \mathrm {(1)} }
In other words, the summation of the work done on the system by the set of external forces is equal to the work stored as strain energy in the elements that make up the system.
The virtual internal work in the right-hand-side of the above equation may be found by summing the virtual work done on the individual elements. The latter requires that force-displacement functions be used that describe the response for each individual element. Hence, the displacment of the structure is described by the response of individual (discrete) elements collectively. The equations are written only for the small domain of individual elements of the structure rather than a single equation that describes the response of the system as a whole (a continuum). The latter would result in an intractable problem, hence the utility of the finite element method. As shown in the subsequent sections, Eq.(1) leads to the following governing equilibrium equation for the system:
R = K r + R o ( 2 ) {\displaystyle \mathbf {R} =\mathbf {Kr} +\mathbf {R} ^{o}\qquad \qquad \qquad \mathrm {(2)} }
where
R {\displaystyle \mathbf {R} } = vector of nodal forces, representing external forces applied to the system's nodes.
K {\displaystyle \mathbf {K} } = system stiffness matrix, which is the collective effect of the individual elements' stiffness matrices : k e {\displaystyle \mathbf {k} ^{e}} .
r {\displaystyle \mathbf {r} } = vector of the system's nodal displacements.
R o {\displaystyle \mathbf {R} ^{o}} = vector of equivalent nodal forces, representing all external effects other than the nodal forces which are already included in the preceding nodal force vector R. These external effects may include distributed or concentrated surface forces, body forces, thermal effects, initial stresses and strains.
Once the supports' constraints are accounted for, the nodal displacements are found by solving the system of linear equations (2), symbolically:
r = K − 1 ( R − R o ) ( 3 ) {\displaystyle \mathbf {r} =\mathbf {K} ^{-1}(\mathbf {R} -\mathbf {R} ^{o})\qquad \qquad \qquad \mathrm {(3)} }
Subsequently, the strains and stresses in individual elements may be found as follows:
ϵ = B q ( 4 ) {\displaystyle \mathbf {\epsilon } =\mathbf {Bq} \qquad \qquad \qquad \qquad \mathrm {(4)} }
σ = E ( ϵ − ϵ o ) + σ o = E ( B q − ϵ o ) + σ o ( 5 ) {\displaystyle \mathbf {\sigma } =\mathbf {E} (\mathbf {\epsilon } -\mathbf {\epsilon } ^{o})+\mathbf {\sigma } ^{o}=\mathbf {E} (\mathbf {Bq} -\mathbf {\epsilon } ^{o})+\mathbf {\sigma } ^{o}\qquad \qquad \qquad \mathrm {(5)} }
where
q {\displaystyle \mathbf {q} } = vector of a nodal displacements--a subset of the system displacement vector r that pertains to the elements under consideration.
B {\displaystyle \mathbf {B} } = strain-displacement matrix that transforms nodal displacements q to strains at any point in the element.
E {\displaystyle \mathbf {E} } = elasticity matrix that transforms effective strains to stresses at any point in the element.
ϵ o {\displaystyle \mathbf {\epsilon } ^{o}} = vector of initial strains in the elements.
σ o {\displaystyle \mathbf {\sigma } ^{o}} = vector of initial stresses in the elements.
By applying the virtual work equation (1) to the system, we can establish the element matrices B {\displaystyle \mathbf {B} } , k e {\displaystyle \mathbf {k} ^{e}} as well as the technique of assembling the system matrices R o {\displaystyle \mathbf {R} ^{o}} and K {\displaystyle \mathbf {K} } . Other matrices such as ϵ o {\displaystyle \mathbf {\epsilon } ^{o}} , σ o {\displaystyle \mathbf {\sigma } ^{o}} , R {\displaystyle \mathbf {R} } and E {\displaystyle \mathbf {E} } are known values and can be directly set up from data input.
Interpolation or shape functions[edit]
Let q {\displaystyle \mathbf {q} } be the vector of nodal displacements of a typical element. The displacements at any other point of the element may be found by the use of interpolation functions as, symbolically:
u = N q ( 6 ) {\displaystyle \mathbf {u} =\mathbf {N} \mathbf {q} \qquad \qquad \qquad \mathrm {(6)} }
where
u {\displaystyle \mathbf {u} } = vector of displacements at any point {x,y,z} of the element.
N {\displaystyle \mathbf {N} } = matrix of shape functions serving as interpolation functions.
Equation (6) gives rise to other quantities of great interest:
Virtual displacements that are a function of virtual nodal displacements: δ u = N δ q ( 6 b ) {\displaystyle \delta \mathbf {u} =\mathbf {N} \delta \mathbf {q} \qquad \qquad \qquad \mathrm {(6b)} }
Strains in the elements that result from displacements of the element's nodes: ϵ = D u = D N q ( 7 ) {\displaystyle \mathbf {\epsilon } =\mathbf {Du} =\mathbf {DNq} \qquad \qquad \qquad \qquad \mathrm {(7)} }
where D {\displaystyle \mathbf {D} } = matrix of differential operators that convert displacements to strains using linear elasticity theory. Eq.(7) shows that matrix B in (4) is
B = D N ( 8 ) {\displaystyle \mathbf {B} =\mathbf {DN} \qquad \qquad \qquad \qquad \mathrm {(8)} }
Virtual strains consistent with element's virtual nodal displacements: δ ϵ = B δ q ( 9 ) {\displaystyle \delta {\boldsymbol {\epsilon }}=\mathbf {B} \delta \mathbf {q} \qquad \qquad \qquad \qquad \mathrm {(9)} }
Internal virtual work in a typical element[edit]
For a typical element of volume V e {\displaystyle V^{e}} , the internal virtual work due to virtual displacements is obtained by substitution of (5) and (9) into (1):
Internal virtual work = ∫ V e δ ϵ T σ d V e = δ q T ∫ V e B T { E ( B q − ϵ o ) + σ o } d V e ( 10 ) {\displaystyle {\mbox{Internal virtual work}}=\int _{V^{e}}\delta {\boldsymbol {\epsilon }}^{T}{\boldsymbol {\sigma }}\,dV^{e}=\delta \ \mathbf {q} ^{T}\int _{V^{e}}\mathbf {B} ^{T}{\big \{}\mathbf {E} (\mathbf {Bq} -\mathbf {\epsilon } ^{o})+\mathbf {\sigma } ^{o}{\big \}}\,dV^{e}\qquad \mathrm {(10)} }
Element matrices[edit]
Primarily for the convenience of reference, the following matrices pertaining to a typical elements may now be defined:
Element stiffness matrix K e = ∫ V e B T E B d V e ( 11 ) {\displaystyle \mathbf {K} ^{e}=\int _{V^{e}}\mathbf {B} ^{T}\mathbf {E} \mathbf {B} \,dV^{e}\qquad \mathrm {(11)} }
Equivalent element load vector Q o e = ∫ V e − B T ( E ϵ o − σ o ) d V e ( 12 ) {\displaystyle \mathbf {Q} ^{oe}=\int _{V^{e}}-\mathbf {B} ^{T}{\big (}\mathbf {E} \mathbf {\epsilon } ^{o}-\mathbf {\sigma } ^{o}{\big )}\,dV^{e}\qquad \mathrm {(12)} }
These matrices are usually evaluated numerically using Gaussian quadrature for numerical integration. Their use simplifies (10) to the following:
Internal virtual work = δ q T ( K e q + Q o e ) ( 13 ) {\displaystyle {\mbox{Internal virtual work}}=\delta \ \mathbf {q} ^{T}{\big (}\mathbf {K} ^{e}\mathbf {q} +\mathbf {Q} ^{oe}{\big )}\qquad \mathrm {(13)} }
Element virtual work in terms of system nodal displacements[edit]
Since the nodal displacement vector q is a subset of the system nodal displacements r (for compatibility with adjacent elements), we can replace q with r by expanding the size of the element matrices with new columns and rows of zeros:
Internal virtual work = δ r T ( K e r + Q o e ) ( 14 ) {\displaystyle {\mbox{Internal virtual work}}=\delta \ \mathbf {r} ^{T}{\big (}\mathbf {K} ^{e}\mathbf {r} +\mathbf {Q} ^{oe}{\big )}\qquad \mathrm {(14)} }
where, for simplicity, we use the same symbols for the element matrices, which now have expanded size as well as suitably rearranged rows and columns.
System virtual work[edit]
Summing the internal virtual work (14) for all elements gives the right-hand-side of (1):
System internal virtual work = ∑ e δ r T ( k e r + Q o e ) = δ r T ( ∑ e k e ) r + δ r T ∑ e Q o e ( 15 ) {\displaystyle {\mbox{System internal virtual work}}=\sum _{e}\delta \ \mathbf {r} ^{T}{\big (}\mathbf {k} ^{e}\mathbf {r} +\mathbf {Q} ^{oe}{\big )}=\delta \ \mathbf {r} ^{T}{\big (}\sum _{e}\mathbf {k} ^{e}{\big )}\mathbf {r} +\delta \ \mathbf {r} ^{T}\sum _{e}\mathbf {Q} ^{oe}\qquad \mathrm {(15)} }
Considering now the left-hand-side of (1), the system external virtual work consists of:
The work done by the nodal forces R: δ r T R ( 16 ) {\displaystyle \delta \ \mathbf {r} ^{T}\mathbf {R} \qquad \mathrm {(16)} }
The work done by external forces T e {\displaystyle \mathbf {T} ^{e}} on the part S e {\displaystyle \mathbf {S} ^{e}} of the elements' edges or surfaces, and by the body forces f e {\displaystyle \mathbf {f} ^{e}}
∑ e ∫ S e δ u T T e d S e + ∑ e ∫ V e δ u T f e d V e {\displaystyle \sum _{e}\int _{S^{e}}\delta \ \mathbf {u} ^{T}\mathbf {T} ^{e}\,dS^{e}+\sum _{e}\int _{V^{e}}\delta \ \mathbf {u} ^{T}\mathbf {f} ^{e}\,dV^{e}}
Substitution of (6b) gives:
δ q T ∑ e ∫ S e N T T e d S e + δ q T ∑ e ∫ V e N T f e d V e {\displaystyle \delta \ \mathbf {q} ^{T}\sum _{e}\int _{S^{e}}\mathbf {N} ^{T}\mathbf {T} ^{e}\,dS^{e}+\delta \ \mathbf {q} ^{T}\sum _{e}\int _{V^{e}}\mathbf {N} ^{T}\mathbf {f} ^{e}\,dV^{e}}
or − δ q T ∑ e ( Q t e + Q f e ) ( 17 a ) {\displaystyle -\delta \ \mathbf {q} ^{T}\sum _{e}(\mathbf {Q} ^{te}+\mathbf {Q} ^{fe})\qquad \mathrm {(17a)} }
where we have introduced additional element's matrices defined below:
Q t e = − ∫ S e N T T e d S e ( 18 a ) {\displaystyle \mathbf {Q} ^{te}=-\int _{S^{e}}\mathbf {N} ^{T}\mathbf {T} ^{e}\,dS^{e}\qquad \mathrm {(18a)} }
Q f e = − ∫ V e N T f e d V e ( 18 b ) {\displaystyle \mathbf {Q} ^{fe}=-\int _{V^{e}}\mathbf {N} ^{T}\mathbf {f} ^{e}\,dV^{e}\qquad \mathrm {(18b)} }
Again, numerical integration is convenient for their evaluation. A similar replacement of q in (17a) with r gives, after rearranging and expanding the vectors Q t e , Q f e {\displaystyle \mathbf {Q} ^{te},\mathbf {Q} ^{fe}} :
− δ r T ∑ e ( Q t e + Q f e ) ( 17 b ) {\displaystyle -\delta \ \mathbf {r} ^{T}\sum _{e}(\mathbf {Q} ^{te}+\mathbf {Q} ^{fe})\qquad \mathrm {(17b)} }
Assembly of system matrices[edit]
Adding (16), (17b) and equating the sum to (15) gives: δ r T R − δ r T ∑ e ( Q t e + Q f e ) = δ r T ( ∑ e k e ) r + δ r T ∑ e Q o e {\displaystyle \delta \ \mathbf {r} ^{T}\mathbf {R} -\delta \ \mathbf {r} ^{T}\sum _{e}(\mathbf {Q} ^{te}+\mathbf {Q} ^{fe})=\delta \ \mathbf {r} ^{T}{\big (}\sum _{e}\mathbf {k} ^{e}{\big )}\mathbf {r} +\delta \ \mathbf {r} ^{T}\sum _{e}\mathbf {Q} ^{oe}}
Since the virtual displacements δ r {\displaystyle \delta \ \mathbf {r} } are arbitrary, the preceding equality reduces to:
R = ( ∑ e k e ) r + ∑ e ( Q o e + Q t e + Q f e ) {\displaystyle \mathbf {R} ={\big (}\sum _{e}\mathbf {k} ^{e}{\big )}\mathbf {r} +\sum _{e}{\big (}\mathbf {Q} ^{oe}+\mathbf {Q} ^{te}+\mathbf {Q} ^{fe}{\big )}}
Comparison with (2) shows that:
The system stiffness matrix is obtained by summing the elements' stiffness matrices:
K = ∑ e k e {\displaystyle \mathbf {K} =\sum _{e}\mathbf {k} ^{e}}
The vector of equivalent nodal forces is obtained by summing the elements' load vectors:
R o = ∑ e ( Q o e + Q t e + Q f e ) {\displaystyle \mathbf {R} ^{o}=\sum _{e}{\big (}\mathbf {Q} ^{oe}+\mathbf {Q} ^{te}+\mathbf {Q} ^{fe}{\big )}}
In practice, the element matrices are neither expanded nor rearranged. Instead, the system stiffness matrix K {\displaystyle \mathbf {K} } is assembled by adding individual coefficients k i j e {\displaystyle {k}_{ij}^{e}} to K k l {\displaystyle {K}_{kl}} where the subscripts ij, kl mean that the element's nodal displacements q i e , q j e {\displaystyle {q}_{i}^{e},{q}_{j}^{e}} match respectively with the system's nodal displacements r k , r l {\displaystyle {r}_{k},{r}_{l}} . Similarly, R o {\displaystyle \mathbf {R} ^{o}} is assembled by adding individual coefficients Q i e {\displaystyle {Q}_{i}^{e}} to R k o {\displaystyle {R}_{k}^{o}} where q i e {\displaystyle {q}_{i}^{e}} matches r k {\displaystyle {r}_{k}} . This direct addition of k i j e {\displaystyle {k}_{ij}^{e}} into K k l {\displaystyle {K}_{kl}} gives the procedure the name Direct Stiffness Method.
See also[edit]
Finite element method
Flexibility method
Matrix stiffness method
Modal analysis using FEM
List of finite element software packages
Structural analysis
Virtual work
Interval finite element
References[edit]
Jump up ^ Matrix Analysis Of Framed Structures, 3rd Edition by Jr. William Weaver, James M. Gere, Springer-Verlag New York, LLC, ISBN 978-0-412-07861-3, 1966
Jump up ^ Theory of Matrix Structural Analysis, J. S. Przemieniecki, McGraw-Hill Book Company, New York, 1968
Jump up ^ Hinton, Ernest; Irons, Bruce (July 1968). "Least squares smoothing of experimental data using finite elements". Strain. 4: 24–27.
Jump up ^ Argyris, J.H and Kelsey, S. Energy theorems and Structural Analysis Butterworth Scientific publications, London
Jump up ^ Clough, R.W, “The Finite Element in Plane Stress Analysis.” Proceedings, 2nd ASCE Conference on Electronic Computations, Pittsburgh, Sep 1960
Jump up ^ The Finite Element Method for Solid and Structural Mechanics, Zienkiewicz O. C and Taylor R L ISBN 978-0-7506-6321-2, 1967, McGraw Hill, New York
Jump up ^ Finite Element Structural Analysis , T.Y Yang, Prentice-Hall, Inc, Englewood, NJ, 1986
Jump up ^ "Wiley: Programming the Finite Element Method, 5th Edition - I. M. Smith, D. V. Griffiths, L. Margetts". eu.wiley.com. Retrieved 2015-09-18.
Retrieved from "https://en.wikipedia.org/w/index.php?title=Finite_element_method_in_structural_mechanics&oldid=784585510"
Categories:
Finite element method
Numerical differential equations
Hidden categories:
Articles needing expert attention with no reason or talk parameter
Articles needing unspecified expert attention
Articles needing expert attention from October 2010
All articles needing expert attention
Navigation menu
Personal tools
Not logged in
Talk
Contributions
Create account
Log in
Namespaces
Article
Talk
Variants
Views
Read
Edit
View history
More
Search
Navigation
Main page
Contents
Featured content
Current events
Random article
Donate to Wikipedia
Wikipedia store
Interaction
Help
About Wikipedia
Community portal
Recent changes
Contact page
Tools
What links here
Related changes
Upload file
Special pages
Permanent link
Page information
Wikidata item
Cite this page
Print/export
Create a book
Download as PDF
Printable version
Languages
Español
Українська
Edit links
This page was last edited on 9 June 2017, at 02:23.
Text is available under the Creative Commons Attribution-ShareAlike License; additional terms may apply. By using this site, you agree to the Terms of Use and Privacy Policy. Wikipedia® is a registered trademark of the Wikimedia Foundation, Inc., a non-profit organization.
Privacy policy
About Wikipedia
Disclaimers
Contact Wikipedia
Developers
Cookie statement
Mobile view

values: mw-page-base noprint

values: mw-head-base noprint

attrs: id class style

values: firstHeading firstHeading en

values: siteSub noprint

values: jump-to-nav mw-jump

values:
plainlinks metadata ambox ambox-content
presentation

values:
/wiki/Coefficient_of_thermal_expansion
mw-redirect
Coefficient of thermal expansion

values:
/wiki/Ray_W._Clough
mw-redirect
Ray W. Clough

values:
/wiki/Pierre-and-Marie-Curie_University
mw-redirect
Pierre-and-Marie-Curie University

values:
/wiki/Stress_(physics)
mw-redirect
Stress (physics)

values:
/wiki/Strain_(materials_science)
mw-redirect
Strain (materials science)

values: Edit section: Element interconnection and displacement

values:
/wiki/Degrees_of_freedom_(engineering)
mw-redirect
Degrees of freedom (engineering)

values: Edit section: Theoretical overview of FEM-Displacement Formulation: From elements, to system, to solution

values: Minimum total potential energy principle

values: vertical-align: -2.338ex; width:44.449ex; height:5.843ex;

values: \mbox{External virtual work} = \int_{V}\delta\boldsymbol{\epsilon}^T \boldsymbol{\sigma} \, dV \qquad \mathrm{(1)}

values: vertical-align: -0.838ex; width:31.18ex; height:2.843ex;

values: \mathbf{R} = \mathbf{Kr} + \mathbf{R}^o \qquad \qquad \qquad \mathrm{(2)}

values: vertical-align: -0.338ex; width:2.014ex; height:2.176ex;

values: vertical-align: -0.338ex; width:2.104ex; height:2.176ex;

values: vertical-align: -0.338ex; width:2.427ex; height:2.343ex;

values: vertical-align: -0.338ex; width:1.113ex; height:1.676ex;

values: vertical-align: -0.338ex; width:3.051ex; height:2.343ex;

values: vertical-align: -0.838ex; width:35.357ex; height:3.176ex;

values: \mathbf{r} = \mathbf{K}^{-1} (\mathbf{R}-\mathbf{R}^o ) \qquad \qquad \qquad \mathrm{(3)}

values: vertical-align: -0.838ex; width:28.985ex; height:2.843ex;

values: \mathbf{\epsilon} = \mathbf{Bq} \qquad \qquad \qquad \qquad \mathrm{(4)}

values: vertical-align: -0.838ex; width:56.127ex; height:2.843ex;

values: \mathbf{\sigma} = \mathbf{E}(\mathbf{\epsilon} - \mathbf{\epsilon}^o)+\mathbf{\sigma}^o = \mathbf{E}(\mathbf{Bq} - \mathbf{\epsilon}^o)+\mathbf{\sigma}^o\qquad \qquad \qquad \mathrm{(5)}

values: vertical-align: -0.671ex; width:1.426ex; height:2.009ex;

values: vertical-align: -0.338ex; width:1.911ex; height:2.176ex;

values: vertical-align: -0.338ex; width:1.767ex; height:2.176ex;

values: vertical-align: -0.338ex; width:1.992ex; height:2.343ex;

values: vertical-align: -0.338ex; width:2.378ex; height:2.343ex;

values: Edit section: Interpolation or shape functions

values: vertical-align: -0.838ex; width:25.071ex; height:2.843ex;

values: \mathbf{u} = \mathbf{N} \mathbf{q} \qquad \qquad \qquad \mathrm{(6)}

values: vertical-align: -0.338ex; width:1.496ex; height:1.676ex;

values: vertical-align: -0.338ex; width:2.102ex; height:2.176ex;

values: vertical-align: -0.838ex; width:28.492ex; height:2.843ex;

values: \delta \mathbf{u} = \mathbf{N} \delta \mathbf{q} \qquad \qquad \qquad \mathrm{(6b)}

values: vertical-align: -0.838ex; width:37.9ex; height:2.843ex;

values: \mathbf{\epsilon} = \mathbf{Du} = \mathbf{DNq} \qquad \qquad \qquad \qquad \mathrm{(7)}

values: vertical-align: -0.338ex; width:2.06ex; height:2.176ex;

values: Strain-displacement relations (page does not exist)

values: vertical-align: -0.838ex; width:30.766ex; height:2.843ex;

values: \mathbf{B} = \mathbf{DN} \qquad \qquad \qquad \qquad \mathrm{(8)}

values: vertical-align: -0.838ex; width:31.282ex; height:2.843ex;

values: \delta \boldsymbol{\epsilon} = \mathbf{B} \delta \mathbf{q} \qquad \qquad \qquad \qquad \mathrm{(9)}

values: Edit section: Internal virtual work in a typical element

values: vertical-align: -0.338ex; width:2.933ex; height:2.343ex;

values: vertical-align: -2.338ex; width:86.593ex; height:5.843ex;

values: \mbox{Internal virtual work} = \int_{V^e}\delta\boldsymbol{\epsilon}^T \boldsymbol{\sigma} \, dV^e = \delta\ \mathbf{q}^T \int_{V^e} \mathbf{B}^T \big\{\mathbf{E}(\mathbf{Bq} - \mathbf{\epsilon}^o)+\mathbf{\sigma}^o\big\} \, dV^e \qquad \mathrm{(10)}

values: vertical-align: -2.338ex; width:30.651ex; height:5.843ex;

values: vertical-align: -2.338ex; width:40.648ex; height:5.843ex;

values: \mathbf{Q}^{oe} = \int_{V^e} - \mathbf{B}^T \big( \mathbf{E}\mathbf{\epsilon}^o - \mathbf{\sigma}^o\big ) \, dV^e \qquad \mathrm{(12)}

values: vertical-align: -1.171ex; width:50.975ex; height:3.509ex;

values: Edit section: Element virtual work in terms of system nodal displacements

values: vertical-align: -1.171ex; width:50.347ex; height:3.509ex;

values: vertical-align: -3.005ex; width:93.232ex; height:5.509ex;

values: \mbox{System internal virtual work} = \sum_{e} \delta\ \mathbf{r}^T \big( \mathbf{k}^e \mathbf{r} + \mathbf{Q}^{oe} \big)  = \delta\ \mathbf{r}^T \big( \sum_{e} \mathbf{k}^e \big)\mathbf{r} + \delta\ \mathbf{r}^T \sum_{e} \mathbf{Q}^{oe}    \qquad \mathrm{(15)}

values: vertical-align: -0.838ex; width:14.984ex; height:3.176ex;

values: \delta\ \mathbf{r}^T \mathbf{R} \qquad \mathrm{(16)}

values: vertical-align: -0.338ex; width:2.875ex; height:2.343ex;

values: vertical-align: -0.338ex; width:2.502ex; height:2.343ex;

values: vertical-align: -0.338ex; width:2.141ex; height:2.343ex;

values: vertical-align: -3.005ex; width:41.01ex; height:6.509ex;

values: \sum_{e} \int_{S^e} \delta\ \mathbf{u}^T \mathbf{T}^e \, dS^e +   \sum_{e} \int_{V^e} \delta\ \mathbf{u}^T \mathbf{f}^e \, dV^e

values: vertical-align: -3.005ex; width:48.645ex; height:6.509ex;

values: \delta\ \mathbf{q}^T \sum_{e} \int_{S^e} \mathbf{N}^T \mathbf{T}^e \, dS^e +   \delta\ \mathbf{q}^T \sum_{e} \int_{V^e} \mathbf{N}^T \mathbf{f}^e \, dV^e

values: vertical-align: -3.005ex; width:32.271ex; height:5.509ex;

values: -\delta\ \mathbf{q}^T \sum_{e} (\mathbf{Q}^{te} +  \mathbf{Q}^{fe}) \qquad \mathrm{(17a)}

values: vertical-align: -2.338ex; width:33.257ex; height:5.843ex;

values: \mathbf{Q}^{te} =  -\int_{S^e} \mathbf{N}^T \mathbf{T}^e \, dS^e  \qquad \mathrm{(18a)}

values: vertical-align: -2.338ex; width:33.638ex; height:5.843ex;

values: \mathbf{Q}^{fe} =  -\int_{V^e} \mathbf{N}^T \mathbf{f}^e \, dV^e  \qquad \mathrm{(18b)}

values: vertical-align: -0.671ex; width:8.605ex; height:3.176ex;

values: vertical-align: -3.005ex; width:32.086ex; height:5.509ex;

values: -\delta\ \mathbf{r}^T \sum_{e} (\mathbf{Q}^{te} +  \mathbf{Q}^{fe}) \qquad \mathrm{(17b)}

values: Edit section: Assembly of system matrices

values: vertical-align: -3.005ex; width:61.213ex; height:5.509ex;

values: \delta\ \mathbf{r}^T \mathbf{R} -\delta\ \mathbf{r}^T \sum_{e} (\mathbf{Q}^{te} +  \mathbf{Q}^{fe}) =  \delta\ \mathbf{r}^T \big( \sum_{e} \mathbf{k}^e \big)\mathbf{r} + \delta\ \mathbf{r}^T \sum_{e} \mathbf{Q}^{oe}

values: vertical-align: -0.338ex; width:2.752ex; height:2.343ex;

values: vertical-align: -3.005ex; width:40.798ex; height:5.509ex;

values: \mathbf{R} = \big( \sum_{e} \mathbf{k}^e \big)\mathbf{r} + \sum_{e} \big( \mathbf{Q}^{oe} + \mathbf{Q}^{te} +  \mathbf{Q}^{fe} \big)

values: vertical-align: -3.005ex; width:11.393ex; height:5.509ex;

values: vertical-align: -3.005ex; width:29.154ex; height:5.509ex;

values: \mathbf{R}^o = \sum_{e} \big( \mathbf{Q}^{oe} + \mathbf{Q}^{te} +  \mathbf{Q}^{fe} \big)

values: vertical-align: -1.171ex; width:2.714ex; height:3.176ex;

values: vertical-align: -0.671ex; width:3.577ex; height:2.509ex;

values: vertical-align: -1.338ex; width:5.236ex; height:3.176ex;

values: vertical-align: -0.671ex; width:4.989ex; height:2.009ex;

values: vertical-align: -0.838ex; width:2.855ex; height:2.843ex;

values: vertical-align: -1.005ex; width:2.871ex; height:3.009ex;

values: vertical-align: -1.005ex; width:2.096ex; height:2.843ex;

values: vertical-align: -0.671ex; width:2.155ex; height:2.009ex;

values:
/wiki/Matrix_stiffness_method
mw-redirect
Matrix stiffness method

values: List of finite element software packages

values: International Standard Book Number (identifier)

values:
reflist mw-references-wrap
list-style-type: decimal;

values: Category:Articles needing expert attention with no reason or talk parameter

values: Category:Articles needing unspecified expert attention

values: Category:Articles needing expert attention from October 2010

values: Category:All articles needing expert attention

values: Discussion about edits from this IP address [ctrl-alt-n]

values: A list of edits made from this IP address [ctrl-alt-y]

values: You are encouraged to create an account and log in; however, it is not mandatory

values: You're encouraged to log in; however, it's not mandatory. [ctrl-alt-o]

values: Discussion about the content page [ctrl-alt-t]

values: p-namespaces navigation vectorTabs p-namespaces-label

values:
p-variants
navigation
vectorMenu emptyPortlet
p-variants-label

values: Past revisions of this page [ctrl-alt-h]

values: p-views navigation vectorTabs p-views-label

values:
p-cactions
navigation
vectorMenu emptyPortlet
p-cactions-label

values:
mw-wiki-logo
/wiki/Main_Page
Visit the main page

values:
/wiki/Main_Page
Visit the main page [ctrl-alt-z]
z

values: Featured content – the best of Wikipedia

values: Find background information on current events

values:
/wiki/Special:Random
Load a random article [ctrl-alt-x]
x

values: portal navigation p-navigation p-navigation-label

values: Guidance on how to use and edit Wikipedia

values: About the project, what you can do, where to find things

values: A list of recent changes in the wiki [ctrl-alt-r]

values: portal navigation p-interaction p-interaction-label

values: List of all English Wikipedia pages containing links to this page [ctrl-alt-j]

values: Recent changes in pages linked from this page [ctrl-alt-k]

values:
/wiki/Wikipedia:File_Upload_Wizard
Upload files [ctrl-alt-u]
u

values: A list of all special pages [ctrl-alt-q]

values: Permanent link to this revision of the page

values: Link to connected data repository item [ctrl-alt-g]

values: Printable version of this page [ctrl-alt-p]

values:
uls-settings-trigger
Language settings
0
button
true

values: Método de los elementos finitos en la mecánica estructural – Spanish

values: Метод скінченних елементів у механіці конструкцій – Ukrainian

values: /static/images/wikimedia-button-1.5x.png 1.5x, /static/images/wikimedia-button-2x.png 2x

values: /static/images/poweredby_mediawiki_132x47.png 1.5x, /static/images/poweredby_mediawiki_176x62.png 2x

values: mediawiki ltr sitedir-ltr mw-hide-empty-elt ns-0 ns-subject page-Finite_element_method_in_structural_mechanics rootpage-Finite_element_method_in_structural_mechanics skin-vector action-view