The behavior of many engineering structures, excited from external stresses, is governed by differential equations. We can solve most of these differential equations analytically as long as we oversimplify or idealize the structure in consideration. But unfortunately, the design geometry, material properties, and loadings of these structures are complex in reality and idealization may not give the correct response due to external loads, hence the need to break the structural problem into subproblems, get the solution of each individual subproblem and assemble the various solutions into a global solution can be achieved using a numerical technique called finite element method.
The finite element method(FEM) is a numerical method for solving differential equations of boundary value problems(also known as field problems). The field normally represents a physical structure constrained by boundary conditions.
It is worth remembering that studying or analyzing an engineering structure or phenomenon with a finite element method is what is referred to as finite element analysis.
The finite element method (FEM) is a computational technique used to solve boundary value problems in engineering. Boundary value problems are also called field problems. The field is the domain of interest and most often represents a physical structure. The field variables are the dependent variables of interest governed by the differential equation. The boundary conditions are the specified values of the field variables (or related variables such as derivatives) on the boundaries of the field. The point in the finite element in which field variables are calculated is called a node. The values of the field variable computed at the nodes are used to approximate the values at non-nodal points by interpolation of the nodal values.
These are functions that describe the variations of the field variables(dependent variables) in the finite element. They are normally polynomial forms of the independent variables, derived to satisfy certain required conditions at the nodes.
The stiffness matrix encodes the characteristic of a finite element. Considering engineering structure finite elements, the stiffness matrix will contain information relating to the geometry and material behavior of the structure. This is like the data representing the structural system. This information can indicate the resistance of the element to deformation when subjected to loading. Such deformation may include axial, bending, shear, and torsional response.
1d element
2d element
3. 3d element
In this stage, we define the following:
In this stage, we:
We are going to analyze a cantilever beam of length 5m. It has fixed support at point A and it's free at the other end. A uniformly distributed load of 10kN/m is acting throughout the span.
This is easily implemented using a Python finite element library called anastruct.
from anastruct import SystemElements
ss = SystemElements()
ss.add_element(location=[[0, 0], [5, 0]])
ss.add_support_fixed(node_id=1)
ss.q_load(element_id=1, q=-10)
ss.show_structure()
ss.solve()
ss.show_shear_force()
ss.show_bending_moment()
ss.show_displacement()
The ability of an engineer to have the skills of finite element methods in his/her toolbox will help in analyzing complex structures that analytical methods can not easily solve due to assumptions and oversimplification. And with computers becoming more powerful, the practice of utilizing finite element methods will continue to increase.