The FEA model runs. The solver converges without warnings. The results are presented to the design team, the safety factor looks adequate, and the design is approved. Six months later, during testing or worse, during service, something fails in a way that the analysis did not predict. The investigation that follows invariably finds one or more of the same categories of error that appear in this article, committed during the analysis phase and not detected before the decision was made.
What makes FEA errors particularly dangerous is not that they are hard to understand once identified. It is that many of them produce results that look entirely plausible. The stress contour map has smooth gradients. The deformed shape looks reasonable. The solver did not report any errors or warnings. The peak stress is in a location that makes intuitive sense. The only problem is that the actual stress is three times higher, or the failure mode is entirely different, or the model is six times stiffer than reality because of an over-constraining boundary condition that was never questioned.
This article covers 16 specific FEA errors organized into six categories, each explained with the root cause that produces it, the magnitude and type of result error it causes, and the specific detection method that will catch it before it leads to a wrong engineering decision. The goal is not just to document mistakes but to give engineers the diagnostic toolkit to find and correct them systematically.
Error Overview: 16 Mistakes Mapped to Impact and Detection
The following table maps all 16 errors covered in this article to their typical result impact, detection method, and severity. Use it as a quick reference during model review, and refer to the detailed sections for each error category for the full technical explanation.
| Error Category | Specific Error | Typical Impact on Results | Detection Method | Severity |
| Unit system | Mixing mm and m, N and kN in same model | Factor of 1000 to 1,000,000 error on stresses | Reaction force check vs expected; dimensional sanity check | Critical – always catastrophic |
| Boundary conditions | Over-constraining with fixed wall instead of pin | 50-500% stiffness overestimation in bending | Deformation mode review; compare reaction moments vs applied | High – systematic error |
| Boundary conditions | Under-constraining – rigid body motion | Solver failure or near-singular matrix | Modal analysis with 6 zero-frequency modes expected | Critical – analysis is invalid |
| Boundary conditions | Artificial stiffness from enforced displacement on unintended DOF | Local stress artifacts near constraint; global stiffness wrong | Remove constraint and observe deformation change | High |
| Material properties | Wrong Young’s modulus (10x too high/low) | Displacements off by 10x; stress unchanged if load-controlled | Verify against published data; check material units | High |
| Material properties | Linear material used beyond yield | Predicted stress above Sy with no yielding shown | Check peak von Mises vs Sy; run with elastoplastic model | Very High – unsafe |
| Material properties | Unit system error in material (GPa vs MPa) | 1000x error on stress; displacement changes by 1000x | Sanity check displacement magnitude vs expected | Critical |
| Connection | Bonded instead of frictional contact | Artificially high force transfer; no sliding captured | Check contact pressure distribution; sliding in physical test? | High |
| Connection | Disconnected mesh nodes at part interface | Load not transferred; stress concentration at gap | Plot deformed shape; check force transfer through interface | Critical – load path wrong |
| Mesh | Coarse mesh at stress concentration | Peak stress underestimated by 50-90% | Mesh convergence study at concentration | High |
| Mesh | Linear tet in bending-dominated region | Bending stiffness 2-5x too high; stress wrong | Switch to quadratic tet; re-run and compare | High |
| Loads | Wrong load direction (global vs local coordinates) | Completely wrong deformation mode | Deformation mode check; verify against expected | Critical – completely wrong |
| Loads | Load magnitude in wrong unit (N vs kN) | 1000x error on all stresses | Equilibrium check; compare reaction to applied | Critical |
| Analysis type | Linear used for geometric nonlinear problem | Stiffness wrong; snap-through missed | Check displacement/dimension ratio; run NL and compare | High |
| Post-processing | Reading averaged instead of unaveraged stress | Peak stress artificially reduced by averaging | Switch to unaveraged; check gradient across element | High – masks failure risk |
| Post-processing | Von Mises instead of principal for brittle failure | Wrong failure criterion applied | Check failure mode; use max principal for brittle | High |

The severity classification reflects the potential for the error to lead to an engineering decision that would be different if the model were correct. Critical errors produce results so wrong that no engineering decision made from them should be trusted. High errors produce systematically biased results that may lead to unconservative or over-conservative decisions. The most dangerous errors are those that produce plausible-looking results that do not trigger the analyst’s suspicion.
Error Category 1: Unit System Inconsistency
Unit system errors are the most catastrophically damaging FEA mistakes because they produce errors by factors of 1,000, 1,000,000, or more, in results that look entirely reasonable in magnitude because the analyst has no independent reference for what the correct answer should be. A unit system error is a silent multiplier that scales every result in the model by a constant factor without triggering any solver warning, any convergence issue, or any plausibility check that is not deliberately applied by the analyst.
How Unit System Errors Happen
FEA solvers do not have a built-in unit system. They process numbers. The solver does not know whether the number 210,000 you entered as Young’s modulus represents 210,000 MPa (correct for steel in MPa units) or 210,000 Pa (steel modulus 1,000,000 times too low) or 210,000 GPa (steel modulus 1,000 times too high). The solver accepts whatever numbers you provide and produces results in the same unit system those numbers imply. If you enter modulus in MPa, forces in N, and geometry in mm, the solver returns stresses in MPa and displacements in mm. If you mix these units, the results are in whatever undefined mixed unit system your inputs created.
The most common mixing error: geometry imported from a CAD system in millimeters, forces applied in kilonewtons (as copied from a load specification), and Young’s modulus entered in GPa (from a materials datasheet that uses GPa). The solver receives: geometry in mm, forces in kN, modulus in GPa. It computes stresses in GPa·kN/mm^2, which is not a standard unit, and the result is numerically somewhere between 10^3 and 10^6 times the correct stress value depending on the specific combination. The displacement result has the same problem. The contour plot still looks smooth and plausible because the color scale adjusts to whatever range the results happen to cover.
The Detection Method: Four Mandatory Unit Checks
- Consistent unit table before any model is started: write out your unit system explicitly before beginning. For SI: force in N, length in m, stress in Pa. For SI-mm: force in N, length in mm, stress in MPa. For US customary: force in lbf, length in in, stress in psi. Every input to the model must be in this system.
- Dimensional sanity check on displacements: run a quick estimate of the expected displacement before reviewing FEA results. A steel cantilever beam 200mm long, 10mm square, loaded with 10N at the tip should deflect approximately 5mm by beam theory. If your FEA shows 0.005mm or 5000mm, you have a unit error.
- Equilibrium check on reactions: the reaction forces at your boundary conditions must sum to the applied loads. If you applied 1000N and the solver reports a 1.0N reaction, your forces are in kN but you entered them expecting N.
- Stress sanity check: the peak stress should be in a physically plausible range. For steel with 10N applied to a 10mm square bar (cross-section area 100mm^2), the nominal stress is 0.1 MPa. If your FEA shows 100 GPa, there is a unit error in the material or load.
| Unit System Reference Card SI-mm (most common for mechanical engineering): Length: mm | Force: N | Mass: tonne (1000 kg) | Time: s | Stress: MPa (N/mm^2) | Modulus: MPa | Density: tonne/mm^3 (steel: 7.85e-9) | Thermal: mm, C, W/mm.C. SI (structures/civil): Length: m | Force: N or kN | Stress: Pa or kPa | Modulus: Pa or GPa. CRITICAL: never mix mm-geometry with GPa-modulus without explicitly converting. Steel modulus in SI-mm units = 210,000 MPa, NOT 210 GPa. |
Error Category 2: Boundary Condition Errors
Boundary condition errors are the most consequential modeling mistakes for structural accuracy. The boundary conditions define how the structure is supported and loaded, and an incorrect constraint fundamentally changes the structural problem being solved. No amount of mesh refinement or solver sophistication can correct for a boundary condition that does not represent the physical support condition. The mesh quality article in this series can be thought of as optimizing the numerical solution to a mathematical problem; boundary condition errors change the mathematical problem itself.
Over-Constraining: Adding Stiffness That Does Not Exist
The most common boundary condition error in structural FEA is over-constraining: applying more constraint than the physical support actually provides. The classic example is using a fixed wall (all six degrees of freedom constrained: three translations and three rotations) to represent a bolted connection where the bolts provide translational constraint but do not prevent rotation. The fixed constraint provides infinite rotational stiffness at the connection. The real bolted joint provides finite (and often low) rotational stiffness. The result: the FEA model is systematically stiffer than the real structure, producing lower deflections and potentially lower peak bending stresses than will exist in the physical part.
The quantitative impact of over-constraining depends on the structural geometry but can be substantial. For a cantilever beam of moderate slenderness ratio (length/height approximately 10:1), replacing a pin support (translational constraint, no rotational constraint) with a fixed support (full constraint) increases the structural stiffness by approximately 4 times and reduces the tip deflection by 4 times for the same applied load. The peak stress location also shifts from the constraint location to the mid-span in some load cases. A design decision made on fixed-support FEA results for a pin-supported real structure may be unconservative by this factor.
Under-Constraining: Rigid Body Motion
The opposite error, under-constraining, produces an ill-conditioned or singular stiffness matrix that the solver cannot invert. The symptom is either a solver error (negative pivot, singular matrix, rigid body motion detected) or, in some solvers with soft springs for stabilization enabled by default, extremely large displacements that indicate unconstrained motion. Under-constraining typically happens when a model is intended to be symmetric but the symmetry boundary conditions are not correctly applied, when a part is connected to adjacent structure through contact only and the contact has not yet engaged, or when the analyst is modeling a sub-structure and has not fully defined the interface with the surrounding structure.
Detection is straightforward: run a modal analysis (natural frequency extraction) on the model before the static analysis. An unconstrained model will show six rigid body modes at or near zero frequency (three translations and three rotations with zero elastic stiffness). Each unconstrained DOF in the model corresponds to one zero-frequency mode. The mode shapes of the zero-frequency modes directly identify which translations and rotations are unconstrained, pointing to the specific missing boundary conditions.
Artificial Stiffness from Constraint Location and Type
A subtler constraint error that produces locally incorrect results without making the global analysis obviously wrong is applying displacement constraints directly to mesh nodes that are on or adjacent to the region of engineering interest. A fully fixed node creates a stress concentration artifact at the constraint location that is entirely a modeling artifact, not a real physical feature. The reported stress at and near the fixed node is meaningless, and if the stress concentration from the constraint overlaps with the real stress concentration from the geometry, the analyst cannot easily separate the physical and artificial contributions to the peak stress.
The correct approach for any support condition that is not literally a fully fixed rigid wall is to apply constraints through remote points, rigid elements, or multi-point constraints (MPC) that distribute the constraint over a realistic contact area, or to replace the support with a spring element calibrated to the actual support stiffness. This separates the constraint mechanism from the geometry of interest and prevents the artificial stress concentration artifact at the constraint location.
| Boundary Condition Best Practice Before applying any constraint, ask: What physical mechanism provides this support in the real structure? A bolted joint provides translational stiffness and partial rotational stiffness, not infinite rotational stiffness. A bearing provides radial stiffness but typically not axial or moment stiffness unless specifically designed to do so. A weld provides all six DOF. Model the mechanism, not the assumption. When in doubt about the rotational stiffness: run the analysis twice, once with all six DOF fixed and once with only translational DOF fixed (pin-equivalent). The true answer lies between these bounds, and if they differ significantly, the rotational stiffness assumption is important and needs investigation. |

Error Category 3: Material Property Errors
Material property errors are the category that most consistently surprises analysts because FEA software makes entering material properties feel authoritative: you type a number, the software accepts it, and the model runs. The software has no way to know whether the number you entered is correct for your material, in your unit system, at the relevant temperature, in the relevant manufacturing condition. Material property errors are therefore entirely the analyst’s responsibility to detect, and they can remain undetected through mesh convergence studies, equilibrium checks, and deformation mode reviews because these validation steps do not require correct material properties to pass.
Wrong Young’s Modulus: The Invisible Stiffness Error
An incorrect Young’s modulus shifts every displacement and stress in a load-controlled analysis by the ratio of the wrong modulus to the correct modulus. Using 210,000 Pa instead of 210,000 MPa for steel (a factor of 1,000,000 error from GPa-to-Pa confusion) produces displacements 1,000,000 times too large and reactions and stresses that appear incorrect because the structure is effectively compliant rubber rather than steel. This extreme case would be immediately obvious from a displacement sanity check.
More insidious: using 70,000 MPa (aluminum modulus) instead of 210,000 MPa (steel modulus) in a steel component, perhaps because the analyst copied material data from a previous aluminum project. The model runs correctly in every formal sense. Displacements are 3x too high (aluminum is 3x less stiff). Stresses, for a load-controlled analysis, are unchanged (stress = load/area, independent of modulus for statically determinate structures) but deflections affect the load path in statically indeterminate structures. The deformed shape looks qualitatively correct. Only a quantitative check against a hand calculation for the specific deflection catches this error.
Linear Material Beyond Yield: The Most Dangerous Material Error
Using a linear elastic material model in a situation where the true material behavior is elastoplastic is the single most dangerous material error in structural FEA because it produces results that appear to show an adequate safety factor when the real structure has actually yielded and may be near plastic collapse. A linear elastic model reports that the stress at a location is 450 MPa on a steel with a yield strength of 250 MPa. This is a physically impossible result: the real material cannot sustain 450 MPa elastically. But the solver has no knowledge of the yield strength and reports the linear elastic result without comment.
The analyst who reads 450 MPa from a linear elastic model on a 250 MPa yield-strength steel must recognize that the result cannot be the true stress. The real stress is bounded by the yield strength (in the absence of strain hardening), and the real strain is much larger than the linear analysis predicts because the material is yielding and absorbing energy that the linear model assumes is being stored elastically. For a safety assessment, any linear elastic FEA result exceeding the yield strength must be flagged and either investigated with an elastoplastic analysis or assessed using plasticity correction methods (Neuber’s rule, ESED method) specifically developed for this situation.
Temperature-Dependent Properties at Wrong Temperature
Many materials show significant changes in mechanical properties with temperature: Young’s modulus decreases, yield strength decreases, creep rate increases. An FEA analysis using room-temperature material properties for a component operating at 400 to 600 degrees Celsius, a typical turbine blade or exhaust manifold operating condition, may overestimate stiffness by 20 to 40 percent and overestimate yield strength by 50 to 70 percent compared to the actual elevated-temperature properties. The resulting safety factor is fictitious: it reflects the room-temperature material, not the material at operating temperature.
Detection requires knowing the operating temperature and verifying that the material properties in the model match published data at that temperature, not at room temperature. For thermal-structural coupled analyses, the temperature field must be correctly computed and the temperature-dependent material properties must be defined as functions of temperature in the material model, not as single values at one temperature.
Error Category 4: Connection and Contact Errors
In multi-body assemblies, the connection between parts is a modeling decision with direct consequences for load transfer, stress distribution, and overall structural stiffness. A wrong contact assumption is often worse than no contact at all, because it creates a plausible-looking stress distribution that silently transfers load in the wrong way. The most common FEA software defaults to bonded contact for all contact pairs, which means that unless the analyst explicitly changes the contact type, every touching surface in the model is assumed to be rigidly glued to every other touching surface.
Bonded Contact: When It Applies and When It Does Not
Bonded contact treats two surfaces as if they are welded or adhesively bonded with no possibility of separation, sliding, or relative displacement. It is appropriate for: welded connections (if the weld is not being analyzed for integrity), adhesive bonds that will not be stressed beyond their elastic limit, and press-fit interfaces where the contact pressure is sufficient to prevent any relative motion. It is not appropriate for: bolted connections (which can open under tension and slide under shear if the friction is overcome), bearing contacts (which can separate), snap-fit connections (which can disengage), and any interface where the contact state might change during loading.
The specific error from using bonded contact where frictionless or frictional contact is correct: the bonded interface transmits tensile force across the contact faces, which is physically impossible for surfaces that are merely in contact without adhesion. This results in artificially high load transfer across what should be a compression-only interface, changing the stress distribution in both parts and potentially masking a separation condition that would create a stress concentration in the physical assembly that the bonded model never shows.
Disconnected Mesh Nodes: The Silent Load Path Failure
In an assembly model where parts are meshed independently and then positioned in contact, it is possible for mesh nodes at the interface to be very close to each other but not actually connected. If the solver does not detect a contact pair between the surfaces (because the contact definition was omitted or the contact detection tolerance is not set wide enough to find the gap), the model treats the two surfaces as if they are in free air. No load transfers between them.
This error is particularly insidious because the model runs without error, the convergence is good, and the deformed shape may look physically reasonable on the non-loaded side. The error appears as one part moving through another in the deformed shape (interpenetration) or as a complete absence of stress in the part that should be receiving load through the contact interface. Detection: plot the deformed shape with the actual scale factor (1:1, not exaggerated) and look for interpenetration. Check that the force transferred through every interface equals the expected proportion of the applied load.
Error Category 5: Load Application Errors
Load application errors cause the analysis to solve the wrong physical problem. Unlike material errors, which affect the magnitude of results while leaving the qualitative pattern correct, load direction errors can produce a completely different deformation mode and stress distribution from the physically correct solution. A structure loaded in the Z-direction that is analyzed with the load in the X-direction (perhaps because global coordinate axes were confused with local component axes) will show maximum stress in the wrong member, maximum deflection in the wrong direction, and completely wrong reactions.
Load Direction Errors: Global vs Local Coordinates
The most common load direction error is applying a load in the global coordinate system direction when it should be in a local coordinate system direction, or vice versa. Gravity, for example, acts in the global Y (or Z, depending on the model orientation) direction. A component that is angled at 30 degrees to horizontal has gravity acting along that angle in global coordinates, but if the analyst applies gravity as a vertical downward force in global coordinates and the model is oriented with the component vertical in the model space, the load is applied in the wrong direction relative to the component geometry.
The detection method is the deformation mode check: review the deformed shape and ask whether the structure deflects in the direction you would expect given the applied load direction. If the deformation is perpendicular to the expected direction or shows a mode that does not match the loading, a load direction error is the likely cause. For models with distributed pressure loads, plot the load direction vectors as arrows on the model surface and verify they are pointing in the correct direction relative to the geometry.
Pressure Load on Wrong Surface Orientation
Surface pressure loads in FEA are applied normal to the surface element face. The direction of this normal depends on the element face orientation in the mesh. If the element normals are incorrectly oriented (pointing inward instead of outward on the external surface of a pressure vessel, for example), the pressure load is applied inward, which collapses the vessel instead of pressurizing it. This error produces a deformed shape that is immediately recognizable as wrong (inward deflection instead of outward), but only if the analyst inspects the deformed shape with a physically relevant scale factor.
Diagnostic for pressure direction errors: always plot load direction vectors before solving any pressure-loaded model. Most FEA pre-processors allow pressure direction vectors to be displayed as arrows on the mesh surface. Verify that all arrows point in the correct direction (outward for internal pressure, inward for external, always normal to the surface and in the direction the load actually acts) before running the analysis.
Error Category 6: Post-Processing and Interpretation Errors
Post-processing errors are distinct from modeling errors in one important way: the FEA model and its solution are correct, but the results are misread, misinterpreted, or presented in a way that produces wrong engineering conclusions These errors happen after the solver has finished and the results look plausible. They are entirely in the analyst’s hands and require engineering knowledge to avoid.
Averaged vs Unaveraged Stress: The Most Common Interpretation Error
When FEA computes stresses, it computes them at the integration points inside each element, then extrapolates them to the element nodes. At every node shared by multiple elements, there are therefore multiple stress values: one from each adjacent element. These values are generally different because the stress field is discontinuous across element boundaries in FEA. The FEA post-processor can either average these values (producing a smooth, continuous stress contour that artificially suppresses the discontinuity) or display them unaveraged (showing the individual element values with their true discontinuity, which is a measure of the error in the solution).
The error: reporting averaged stress at a stress concentration when unaveraged is more appropriate. Averaging smooths out the peak by mixing the high stress in the high-gradient element with the lower stress in the adjacent coarser element. The reported peak is lower than the true peak by an amount that depends on the element size at the concentration and the severity of the gradient. For a stress concentration with a physically real gradient, the averaged stress underestimates the peak. For a mesh convergence check, the difference between averaged and unaveraged is a direct measure of the mesh quality at that location: a large difference signals that the mesh is too coarse to accurately capture the gradient.
Wrong Stress Quantity for the Failure Mode
Different failure modes require different stress quantities. Using the wrong one can produce a safety factor that is wrong by more than the failure mode factor itself. The key distinctions:
- Von Mises (equivalent stress): appropriate for ductile metal yielding under multiaxial stress states. The von Mises criterion predicts yielding when the distortional strain energy equals the yield strain energy. It is correct for ASME Section VIII pressure vessel analysis (Division 2 uses von Mises as the basis for the Mises yield criterion) and for most structural steel assessments.
- Maximum principal stress: appropriate for brittle fracture assessment and for fatigue in materials where tensile cracks are the initiation mechanism. For a cast iron component (low ductility), the maximum principal stress governs failure, not von Mises. Reporting von Mises for a brittle material can give a safety factor that is 20 to 40 percent non-conservative for biaxial stress states where the principal stresses are not equal.
- Maximum shear stress (Tresca): used for ductile yielding assessment in some codes (ASME Division 1 uses a modified Tresca criterion). The Tresca criterion is more conservative than von Mises by a factor of up to 15 percent for equibiaxial stress states.
- Normal stress perpendicular to weld: for weld fatigue assessment (IIW recommendations, BS 7608), the relevant stress is typically the hot spot stress or the structural stress normal to the expected crack plane, not the von Mises stress at the weld toe.
Scale Factor Errors in Deformed Shape Interpretation
Deformed shape plots in FEA are typically displayed with an exaggerated scale factor (common values: 10x, 100x, 1000x) to make the deformation visible to the human eye when the actual deformation is small compared to the model dimensions. An exaggerated deformed shape is a useful visualization tool, but it cannot be used to assess the magnitude of deformation, the presence of interpenetration, or whether the deformation mode is physically reasonable
The specific error: using a highly exaggerated scale factor to assess contact behavior in an assembly. Two surfaces that appear to separate by a large gap in a 1000x scale factor plot may in fact overlap by 0.001mm in reality, which is physically impossible (interpenetration) and indicates either a mesh contact issue or an overly compliant model. Always switch to 1:1 true scale when assessing whether contact surfaces are behaving physically, whether parts are interpenetrating, or whether the actual displacement magnitude is acceptable.
The Pre-Analysis Checklist: Preventing These Errors Before They Propagate
The majority of the 16 errors in this article are preventable by a systematic pre-analysis and post-analysis review process. The following checklist covers the most critical checks at each stage of the FEA workflow.
| What is the difference between FEA verification and validation? Verification confirms that the numerical solver correctly implements the mathematical model (solving the equations correctly). Validation confirms that the mathematical model represents the physical system accurately (solving the correct equations). In practice: the software vendor is responsible for code verification; the analyst is responsible for solution verification (mesh convergence) and model validation (comparison to analytical solutions or experimental data) on every analysis. |
FEA Error Prevention Checklist |
Frequently Asked Questions
Q: What are the most common FEA mistakes that lead to wrong results?
The most consequential FEA errors fall into six categories: unit system inconsistency (mixing mm with GPa, or N with kN, which creates factors-of-thousands errors in all results), boundary condition errors (over-constraining with fixed supports instead of pins adds artificial stiffness; under-constraining causes rigid body motion), material property errors (wrong modulus, using linear material beyond yield), connection errors (bonded contact where separation or sliding should occur), load direction errors (global vs local coordinate confusion), and post-processing errors (averaged stress suppressing real peak, wrong stress quantity for the failure mode). All of these can produce results that look plausible while being systematically wrong.
Q: How do I detect a unit system error in my FEA model?
Perform a displacement sanity check: before reviewing any stress results, estimate the expected peak deflection using a hand calculation or analytical formula and compare it to the FEA result. A factor-of-1000 discrepancy indicates a unit error in the material modulus or applied force. Also check the reaction forces: they must sum to the applied loads. If you applied 1000 N and the reactions sum to 1.0 N, your forces were entered as kN when the model expects N. Prevent unit errors by writing out your unit system explicitly before building the model and verifying every material property input against its published source with explicit unit confirmation.
Q: What is the difference between averaged and unaveraged stress in FEA?
FEA computes stress at integration points inside each element, then extrapolates to the nodes at element corners. Where multiple elements share a node, each element produces a different stress value at that node because the stress field is discontinuous across element boundaries in finite element analysis. Averaged stress combines these multiple values into a single value at each node, producing a smooth contour. Unaveraged stress shows the individual element values without combining them. The difference between averaged and unaveraged at a location is a mesh quality indicator: a large difference signals that the mesh is too coarse to accurately resolve the stress gradient there. For peak stress reporting at stress concentrations, unaveraged stress is more conservative and more meaningful.
Q: Why does my FEA show stress above the material yield strength?
If a linear elastic material model is used, the FEA solver has no knowledge of the yield strength and will report stresses above yield without any warning. Linear elastic FEA can report any stress value regardless of whether it is physically achievable. Any linear elastic FEA result that exceeds the material yield strength is physically impossible as reported: the real material would have yielded and redistributed the stress. This does not mean the structure is safe – it means the model does not capture the real behavior. Options: run an elastoplastic nonlinear analysis to capture the post-yield behavior, or apply a plasticity correction method (Neuber’s rule) to estimate the true strain from the linear elastic stress result.
Q: How do I know if my FEA boundary conditions are correct?
Run a modal analysis before the static analysis. An unconstrained model will show 6 rigid body modes at near-zero frequency. Each zero-frequency mode represents one missing constraint, and the mode shape shows which translational or rotational direction is unconstrained. For over-constraining: run the analysis twice with different constraint types at the same location (fully fixed vs pin equivalent) and compare results. If the results differ by more than 10-20%, the rotational constraint assumption is significant and must be investigated. Also review the deformed shape: if the structure does not deform in the direction you expect given the applied loads, the boundary conditions are likely wrong.
Q: What is the most dangerous FEA error an engineer can make?
Using a linear elastic material model in a situation where the material is actually yielding under the applied loads. This produces a stress result above yield strength that the analyst may not recognize as physically impossible. The engineer then calculates a safety factor by dividing the reported stress by the yield strength, which gives a safety factor less than 1.0 (indicating imminent failure) but may interpret it as requiring a ‘redesign rather than immediate concern.
The real danger is when the analyst accepts the linear result, perhaps rounding down the peak to an averaged value, and arrives at a result just above yield that looks marginally safe. The actual behavior may involve significant plastic strain, fatigue initiation, and potential progressive failure that the linear model has no mechanism to predict.
Conclusion:
The 16 errors in this article share a common characteristic: every single one of them is predictable, recognizable in pattern, and preventable with the right pre-analysis and post-analysis discipline. They are not random artifacts of software complexity or numerical noise. They are the result of specific modeling decisions that do not correctly represent the physics of the problem, applied in ways that the FEA software cannot detect and cannot warn against.
The engineer’s defense against these errors is not just technical knowledge, though that is necessary. It is the intellectual discipline of questioning every result against an independent reference before accepting it. The equilibrium check is fast and catches load and unit errors. The hand calculation comparison catches magnitude errors. The deformation mode review catches direction and constraint errors. The averaged-vs-unaveraged comparison catches post-processing errors. None of these checks requires additional simulation runs. They require five to fifteen minutes of thoughtful review that transforms a result from an unverified number into a credible engineering evidence.
The engineers who consistently produce reliable FEA results are not those who never make any of these mistakes. They are the ones who have built systematic review habits that catch these mistakes before they propagate into engineering decisions. The checklist in this article is a starting point for building those habits. Apply it to your next analysis. The mistakes it prevents are not hypothetical.
Continue building your FEA competency with our guides on mesh quality and FEA accuracy, when to use linear vs nonlinear FEA, and the validation methods that confirm your results are physically correct.



























