This section reviews the basics of the finite element method. A simple example of a truss, constrained at one end and loaded at the other end as shown in Figure 12, is used to introduce some terms and conventions used in the rest of this document. The objective of the analysis is to find the displacement of the free end of the truss, the stress in the truss, and the reaction force at the constrained end of the truss.
The first step of any finite element simulation is to discretize the actual geometry of the structure using a collection of finite elements. Each finite element represents a discrete portion of the physical structure. The finite elements are joined by shared nodes. The collection of nodes and finite elements is called the mesh. The number of elements per unit of length, area, or volume in a mesh is referred to as the mesh density. In a stress analysis the displacements of the nodes are the fundamental variables that ABAQUS calculates. Once the nodal displacements are known, the stresses and strains in each finite element can be determined easily.
In this case the rod shown in Figure 12 will be modeled with two truss elements. In ABAQUS truss elements can carry axial loads only. The discretized model is shown in Figure 13 together with the node and element labels.
Free-body diagrams for each node in the model are shown in Figure 14. In general each node will carry an external load applied to the model, P, and internal loads, I, caused by stresses in the elements attached to that node. For a model to be in static equilibrium, the net force acting on each node must be zero; i.e., the internal and external loads at each node must balance each other. For node a this equilibrium equation can be obtained as follows.
Assuming that the change in length of the rod is small, the strain in element 1 is given by
Assuming that the material is elastic, the stress in the rod is given by the strain multiplied by the Young's modulus, E:
The axial force acting on the end node is equivalent to the stress in the rod multiplied by its cross-sectional area, A. Thus, a relationship between internal force, material properties, and displacements is obtained:
Equilibrium at node b must take into account the internal forces acting from both elements joined at that node. The internal force from element 1 is now acting in the opposite direction and so becomes negative. The resulting equation is
These equilibrium equations need to be solved simultaneously to obtain the displacements of all the nodes. This requirement is best achieved by matrix techniques; therefore, write the internal and external force contributions as matrices. If the properties and dimensions of the two elements are the same, the equilibrium equations can be simplified as follows:
In general, it may be that the element stiffnesses, the terms, are different from element to element; therefore, write the element stiffnesses as and for the two elements in the model. We are interested in obtaining the solution to the equilibrium equation in which the externally applied forces, P, are in equilibrium with the internally generated forces, I. When discussing this equation with reference to convergence and nonlinearity, we write it as
For the complete two-element, three-node structure we, therefore, modify the signs and rewrite the equilibrium equation as