Every engineer running FEA faces this decision. The model is built, the boundary conditions are defined, and the load case is ready. Now: linear static analysis, which runs in minutes and converges every time, or nonlinear analysis, which may take hours, requires careful setup, and can fail to converge in ways that are genuinely difficult to diagnose? The choice is consequential in both directions: choose linear when the problem is actually nonlinear and you risk results that are not just inaccurate but potentially dangerously non-conservative. Choose nonlinear when linear would have been adequate and you pay a computational and setup cost that serves no engineering purpose.
The answer is not simply to always use nonlinear FEA because it is more realistic. Nonlinear analysis is not more accurate than linear analysis when the physical behavior is genuinely linear. It is slower, more complex to set up, harder to troubleshoot when it fails, and requires more engineering judgment to interpret correctly. Linear analysis, when its assumptions hold, provides exactly correct results with a fraction of the computational effort. The engineering skill is knowing which assumptions hold and when they do not.
This article provides that knowledge in specific, actionable terms. It covers the three types of nonlinearity that drive the decision, the quantitative thresholds that define when each type becomes significant, the computational cost multipliers engineers need to know when making this decision under schedule pressure, the 12-scenario decision table that covers the situations most engineers encounter, how to use a linear pre-check as a diagnostic tool before committing to nonlinear analysis, and the convergence troubleshooting framework for when nonlinear analyses fail.
The Three Linear FEA Assumptions: Understanding What You Are Committing To
Linear FEA rests on three simultaneous assumptions. When all three hold, linear analysis is not just computationally convenient: it is the theoretically correct analytical approach. When any one of the three fails, the linear model produces results that may range from slightly over-conservative to catastrophically wrong depending on the degree and type of nonlinearity present. Understanding exactly what each assumption states, why it breaks down, and what the consequences of its failure are is the foundation of making the linear-versus-nonlinear decision correctly.
Assumption 1: Small Displacements and Small Strains (Geometric Linearity)
Linear FEA assumes that the deformations of the structure are small enough that: the structure’s original geometry adequately represents its deformed geometry for equilibrium calculations, the stiffness of the structure does not change as it deforms, and the strain-displacement relationships remain linear. These three sub-assumptions together are called geometric linearity or the small-displacement assumption.
When this assumption holds, the stiffness matrix [K] computed at the undeformed configuration is the same stiffness matrix that applies throughout the loading, and the equilibrium equation [K]{u} = {F} is solved once to give the complete displacement field. When it fails, the stiffness matrix changes as the structure deforms, and the equilibrium equation must be re-solved at each increment of load with an updated stiffness matrix.
The practical threshold for geometric nonlinearity is generally: if the maximum deflection exceeds 5 to 10 percent of the smallest characteristic structural dimension (the minimum cross-sectional dimension, the plate thickness, the beam height), or if the maximum strain exceeds 1 to 5 percent, geometric nonlinearity is likely to produce results that differ meaningfully from linear analysis. Below these thresholds, geometric nonlinearity effects are typically less than 5 percent of the linear result, which is within the accuracy range of FEA for most engineering purposes.
For stocky, compact structures loaded in their primary stiffness direction, this assumption is rarely violated at service load levels. For slender structures, flexible membranes, snap-fit features, and any structure where the loading direction changes relative to the deformed shape, this assumption should be examined explicitly before committing to a linear analysis.
| Geometric Nonlinearity Quick Test Step 1: Run a linear static analysis. Step 2: Check the maximum displacement. Step 3: Divide by the smallest characteristic dimension (beam height, plate thickness, shell radius). Step 4: If the ratio exceeds 0.05 (5%), re-run with geometric nonlinearity enabled (large displacement option in your solver). Compare results. If they differ by more than 5%, the nonlinear result is more accurate and you should use it. If they are within 5%, linear is adequate for geometric effects. |
Assumption 2: Linear Elastic Material Behavior (Material Linearity)
Linear FEA assumes that the relationship between stress and strain is linear throughout the structure for the entire loading history. This is Hooke’s Law: stress equals Young’s modulus times strain, and the material returns to its original shape when the load is removed. The elastic modulus E is constant, and there is no plastic deformation, no creep, no viscosity, and no damage accumulation.
This assumption breaks down when the stress anywhere in the structure reaches or exceeds the material’s yield strength. Beyond yield, the stress-strain relationship is no longer linear: the material work-hardens (in the case of most metals) or softens, and the stiffness changes. More critically, the strain at a point is no longer uniquely determined by the current stress: it depends on the loading history, which linear analysis has no mechanism to capture.
For steel components at service load levels, material linearity is often an excellent assumption: structural steel has a well-defined linear elastic region up to approximately 250 MPa (for mild steel) or 690 MPa (for high-strength steel). If the peak von Mises stress in a linear analysis remains below approximately 70 to 80 percent of the yield strength, material nonlinearity effects are negligible. Above this threshold, and certainly above yield, material nonlinearity must be addressed explicitly.
Elastomers, rubber, biological tissues, and polymer foams are material nonlinear by nature: their stress-strain behavior is nonlinear even at small strains because these materials are inherently non-Hookean. These materials require hyperelastic material models (Mooney-Rivlin, Ogden, Neo-Hookean) that capture the nonlinear stress-strain behavior from the beginning of loading, not just after a yield point is reached.
Assumption 3: Linear Boundary Conditions (Contact and Constraint Linearity)
Linear FEA assumes that the boundary conditions, the points and surfaces where loads are applied and where motion is constrained, remain constant and fixed throughout the analysis. All contacts between parts are either fully bonded (no separation, no sliding) or fully free (no interaction). Loads are applied at fixed points in fixed directions. Constraints do not change as the structure deforms.
This assumption fails whenever contact between surfaces changes during loading: two surfaces that start in contact may separate under certain load conditions, two surfaces that start separated may come into contact, and sliding contact introduces friction-dependent tangential forces at the interface. It also fails when applied forces change direction as the structure deforms (follower forces), or when boundary conditions are load-dependent.
Contact nonlinearity is computationally the most challenging of the three types because it introduces discontinuous changes in the stiffness matrix: when a node goes from not-in-contact to in-contact, the number of active constraints in the model changes suddenly, causing numerical instability in the iterative solution process. This is why bolted joint analyses, press-fit analyses, rubber seal compression, and bearing-race interfaces are all inherently nonlinear: the contact state is part of the solution, not a known input.
The Three Nonlinearity Types: Engineering Examples and Thresholds
Understanding the three assumption types in the abstract is necessary but not sufficient. The engineering value comes from being able to recognize, in a specific design and loading scenario, which type of nonlinearity is present and how significant it is likely to be.
Geometric Nonlinearity: When Structure Changes Its Own Stiffness
The most common scenario for geometric nonlinearity is a slender structure under lateral load. A fishing rod bent by the weight of a fish, a snap-fit clip at its maximum deflection, a thin-walled tube under internal pressure combined with bending, a membrane structure under applied pressure: in all of these cases, the deformed shape of the structure is significantly different from the original shape, and the equilibrium of forces must be assessed in the deformed configuration.
There are two specific structural behaviors that cannot be captured at all by linear analysis: snap-through buckling and membrane stiffening. Both are purely geometric phenomena that depend on the nonlinear relationship between deformation and stiffness.
Snap-through buckling occurs when a shallow arch or dome structure under increasing central load reaches a critical point where the arch suddenly snaps through to an inverted configuration. Linear analysis finds a smooth, monotonically increasing displacement response up to and including the snap-through point. It completely misses the snap-through itself because snap-through is a geometric instability that depends on the changing stiffness of the deforming arch, which linear analysis assumes to be constant.
Membrane stiffening (also called geometric stiffening or stress stiffening) occurs when a flexible membrane develops significant in-plane tensile stress under transverse loading, and that tensile stress contributes to the membrane’s resistance to further transverse deflection. A stretched cable or a drum skin becomes stiffer as it is loaded because the tensile load in the membrane stiffens it against transverse forces. Linear analysis, which does not account for the stiffness contributed by membrane stresses, significantly underpredicts the load-carrying capacity of these structures.
Material Nonlinearity: Beyond Hooke’s Law
Material nonlinearity matters most when the design intent is to allow controlled plastic deformation, when the structure is loaded to failure, or when the material is inherently nonlinear (rubber, elastomers, some polymers, biological tissues). For these applications, linear analysis produces results that are not just inaccurate but actively misleading because the predicted stresses exceed the yield strength without the analysis having any mechanism to redistribute the excess load through plastic flow.
In metal forming and manufacturing simulations, the entire purpose of the process is to cause plastic deformation: the blank is plastically formed into the desired shape. Linear analysis is fundamentally inapplicable. The material model must include both the elastic and plastic regions of the stress-strain curve, and the analysis must track the accumulated plastic strain as the blank progressively deforms into the die.
For structural integrity and fitness-for-service assessments of pressure vessels and piping under overload conditions, material nonlinearity is required by ASME Section VIII Division 2 when using the Direct Route to Design by Analysis. The limit load calculation, which determines the load at which the structure reaches plastic collapse, explicitly requires an elastoplastic material model with no hardening (elastic-perfectly plastic) to identify the load at which the structure can no longer maintain equilibrium without unlimited plastic flow.
Contact Nonlinearity: When Interfaces Define the Solution
Contact nonlinearity is present in virtually every assembled mechanical structure: bolted connections, press-fit interfaces, bearing contacts, snap-fit joints, rubber seals against housing surfaces, gear tooth contacts. Any interface where the contact pressure, contact area, and sliding behavior are part of the structural response rather than known inputs to the analysis requires contact nonlinearity.
The specific challenges of contact nonlinearity extend beyond the three-body problem of contact mechanics. Friction at contact interfaces introduces load-path dependence: the contact forces depend not just on the current load but on the sequence of loading, because frictional interfaces lock in shear forces that persist after the load is removed. This makes contact analysis with friction inherently history-dependent and requires incremental loading even when the final load state is static.
The Decision Framework: 12 Engineering Scenarios
The following table maps 12 common engineering analysis scenarios to the correct analysis type, the primary indicator to check, and the solver approach required. Use this as a starting reference for any new analysis before reviewing the detailed sections below.
| Situation | Linear FEA Valid? | Nonlinear Type Needed | Primary Indicator | Solver Approach |
| Steel bracket under service load, stress < 0.7 Sy | Yes | None | Von Mises stress well below yield | Linear static |
| Rubber seal under compression | No | Material + geometric | Elastomeric material, large strain | Nonlinear static, hyperelastic material model |
| Snap-fit clip at maximum deflection | No | Geometric (large deformation) | Deflection > 5% of characteristic length | Nonlinear static, large displacement |
| Bolted joint assembly with preload | No | Contact nonlinearity | Parts interact and may separate | Nonlinear static with contact |
| Sheet metal forming simulation | No | Geometric + material + contact | Plastic yielding, large strain, die contact | Nonlinear with plasticity and contact |
| Post-buckling structural response | No | Geometric nonlinearity | Stiffness changes after buckling load | Nonlinear static or Riks method |
| Crash simulation (impact < 100ms) | No | Geometric + material + contact + dynamic | High strain rate, dynamic inertia effects | Explicit dynamic nonlinear |
| Vibration mode shapes of a stiff structure | Yes (modal) | None if stress < 0.5 Sy | Frequency and mode shape extraction | Linear modal analysis |
| Pressure vessel below design pressure (ASME Div 1) | Yes (by code) | None for code compliance | Code-mandated linear elastic stress basis | Linear static, stress categorization |
| Pressure vessel above 2/3 yield (limit load check) | No | Material nonlinearity | Plastic collapse assessment required | Nonlinear with elastoplastic material |
| Thin shell buckling (Euler column type) | Linear buckling first | Geometric (post-buckling optional) | Linear eigenvalue buckling, then verify | Linear buckling eigenvalue + optional NL |
| Biological soft tissue under load | No | Geometric + material + contact | Viscoelastic, large deformation | Nonlinear with viscoelastic or hyperelastic |
Reading this table: the green-highlighted scenarios in the ‘Linear FEA Valid?’ column are cases where linear analysis is justified by engineering assessment. The red-highlighted scenarios require nonlinear analysis. Note that code-mandated linear analysis (ASME Div 1 pressure vessels) appears as a green case even though the actual physical behavior may include some nonlinearity, because the design code’s safety factors are calibrated for linear elastic stress analysis and a nonlinear analysis used in this code context would require a different evaluation methodology (Div 2 Direct Route).
Computational Cost of Nonlinear FEA: What You Are Actually Paying For
The decision between linear and nonlinear analysis is not purely a matter of accuracy. It is an engineering decision that includes computational cost, setup time, result interpretation complexity, and convergence risk. Understanding the computational cost multipliers of nonlinear analysis helps engineers and engineering managers make this decision with realistic expectations about the time and resources required.
| Analysis Type | Solve Time (vs Linear) | Memory (vs Linear) | Iteration Method | Convergence Risk |
| Linear static | 1x (baseline) | 1x (baseline) | Direct solver, single pass | None – always converges |
| Linear buckling (eigenvalue) | 2x-5x | 1.5x-2x | Eigenvalue extraction (Lanczos) | Low – eigenvalue extraction is robust |
| Nonlinear static (geometric NL only) | 3x-15x | 1.5x-3x | Newton-Raphson incremental | Medium – diverges at snap-through points |
| Nonlinear static (material + geometric) | 5x-30x | 2x-5x | Newton-Raphson with arc-length | Medium-High – plasticity causes slow convergence |
| Nonlinear static (contact) | 5x-50x | 2x-4x | Lagrange multiplier or penalty | High – contact opening/closing causes instability |
| Fully nonlinear (all three types) | 10x-100x | 3x-8x | Newton-Raphson or explicit | Very High – requires expert setup and monitoring |
| Implicit nonlinear dynamic | 20x-200x | 4x-10x | Newmark-Beta or HHT-alpha | High – time step must meet stability requirements |
| Explicit nonlinear dynamic | 5x-50x (per cycle) | 2x-4x | Central difference (conditionally stable) | Low per step but needs very small time steps |
Why Nonlinear FEA Costs So Much More
The fundamental reason nonlinear analysis is so much more expensive than linear analysis is the need for incremental, iterative solution. A linear analysis solves one matrix equation: [K]{u} = {F}. The stiffness matrix K is computed once, factored once, and the displacement vector u is computed in one back-substitution operation. The cost is dominated by the matrix factorization, which scales roughly as O(n^1.5) for sparse matrices, where n is the number of degrees of freedom.
A nonlinear analysis solves this equation many times. The load is applied in increments (typically 10 to 100 load steps), and within each load step, the solution is iterated using the Newton-Raphson method until the residual (the difference between internal and external forces) falls below a specified tolerance. Each Newton-Raphson iteration requires recomputing the tangent stiffness matrix (which has changed because the geometry or material state has changed), re-factoring it, and performing another back-substitution. A nonlinear analysis with 50 load steps and 5 Newton-Raphson iterations per step requires 250 matrix factorizations compared to the single factorization in a linear analysis.
Contact Nonlinearity: The Most Computationally Expensive Case
Contact problems are computationally the most expensive because the active contact set can change from one Newton-Raphson iteration to the next. A node that was not in contact in one iteration may come into contact in the next, changing the number of active constraints and requiring the stiffness matrix to be modified. This causes chattering: nodes oscillating between in-contact and not-in-contact states without converging, which the solver must detect and address by stabilizing the contact behavior. Different FEA solvers (Ansys, Abaqus, NASTRAN, LS-DYNA) implement contact stabilization through different algorithms (augmented Lagrangian, penalty method, Lagrange multiplier) with different stability and accuracy tradeoffs.
Explicit vs Implicit Solvers: The Dynamic Nonlinear Choice
For dynamic nonlinear problems, engineers choose between implicit solvers (Newmark-Beta, HHT-alpha, backward difference) and explicit solvers (central difference). This choice has profound consequences for both computational cost and the types of problems each can handle.
Implicit solvers are unconditionally stable: they can use large time steps without numerical instability. They require matrix factorization at each time step but can advance through time efficiently for quasi-static or low-frequency dynamic problems. They are used for slow events: forming processes, material testing, structural response to slowly applied loads, and seismic analysis of buildings.
Explicit solvers are conditionally stable: they require a time step smaller than the Courant-Friedrichs-Lewy (CFL) condition, which is approximately equal to the element’s smallest dimension divided by the wave speed in the material. For steel (wave speed approximately 5,000 m/s), a 1mm element requires a time step smaller than 0.0000002 seconds. This very small time step makes explicit solvers suitable for fast events (crash simulations, blast loading, high-speed impact) where the event itself happens in milliseconds and the time step is naturally small relative to the event duration.
Read more on Scan-to-CAD vs Manual Modeling: Cost Comparison
The Linear Pre-Check: A Diagnostic Tool Before Committing to Nonlinear
One of the most time-efficient workflows in FEA practice is using a linear analysis as a diagnostic tool to determine whether nonlinear analysis is actually necessary. This pre-check workflow is standard practice among experienced FEA analysts but rarely documented explicitly in educational content. It avoids the expensive, time-consuming, and potentially non-converging nonlinear analysis when the structure actually behaves linearly under the given loading.
The Linear Pre-Check Workflow
- Run linear static analysis: Set up the model with linear material properties, no contact (bond all interfaces initially), and small displacement assumptions. Solve and extract results.
- Check maximum displacement ratio: Compute max displacement / characteristic dimension. If > 0.05, geometric nonlinearity is potentially significant.
- Check von Mises stress vs yield: Identify the maximum von Mises stress. If > 0.7 x Sy at any location, material nonlinearity is potentially significant. If > Sy, it is definitely significant.
- Check contact interfaces: Identify all interfaces where contact is assumed. Any interface where parts may separate or slide under load requires contact nonlinearity.
- Make the decision: If none of the above thresholds are exceeded, linear analysis is adequate. If any are exceeded, proceed to nonlinear with the specific nonlinearity type identified by the pre-check results.
- Run nonlinear and compare: For validation purposes, compare linear and nonlinear results on the first nonlinear analysis for a given part family. If the results agree within 5 to 10 percent, linear may be adequate for future similar analyses of that part family.
| The Pre-Check Rule of Thumb Linear is likely adequate when: max displacement / characteristic dimension < 5%, max von Mises < 70% Sy, no changing contact conditions. Nonlinear is required when: any of these thresholds is exceeded, the material is a rubber or elastomer, the loading involves dynamic impact or fast transients, or the structural response involves buckling or snap-through. When in doubt: run both for one representative load case. The additional time is almost always justified by the confidence it provides in the analysis strategy. |
Nonlinear FEA Convergence Troubleshooting
The single most common practical challenge in nonlinear FEA is convergence failure: the Newton-Raphson iteration does not converge within the maximum number of iterations at a specific load increment. Convergence failure is not random: it occurs at specific load levels for specific physical reasons, and understanding those reasons is the key to resolving it.
What Convergence Failure Actually Means
Convergence failure means that the solver could not find a displacement increment that reduces the residual (force imbalance) below the convergence tolerance within the maximum allowed number of iterations. This happens for three distinct reasons: the physical system is genuinely unstable at this load level (snap-through, plastic collapse), the numerical model has setup problems (poorly constrained model, inconsistent contact definitions, material model extrapolated beyond its calibrated range), or the increment size is too large for the solver to converge within the iteration limit even though the physical behavior is stable.
Convergence Troubleshooting Checklist
Nonlinear FEA Convergence Troubleshooting |
Code-Mandated Analysis Types: When the Standard Decides for You
In certain regulated engineering domains, the choice between linear and nonlinear FEA is not made by the analyst based on engineering judgment: it is specified by the applicable design code. Understanding which codes mandate which analysis types is essential for engineers working in pressure vessels, nuclear, aerospace, and civil structural applications.
ASME Section VIII: Two Division Philosophy
ASME Boiler and Pressure Vessel Code Section VIII governs the design of pressure vessels in two divisions with fundamentally different analytical philosophies. Division 1 uses linear elastic stress analysis combined with stress categorization (primary, secondary, peak stresses) and prescriptive safety factors. The analysis methodology is built on linear FEA results, and the evaluation criteria are specific to linear elastic results. Using a nonlinear analysis under Division 1 methodology would produce results that cannot be directly compared to the Division 1 acceptance criteria without additional interpretation.
Division 2 offers an alternative using the Direct Route to Design by Analysis, which explicitly permits and in some cases requires nonlinear analysis. The Direct Route includes limit load analysis (elastoplastic material, load factor approach) and elastic-plastic analysis (full nonlinear analysis at factored loads) as assessment methods. An analyst using Division 2 Direct Route is explicitly expected to perform nonlinear analysis for overload and progressive plastic deformation assessments.
Aerospace and Structural Codes
Aerospace structural analysis for primary structure typically follows linear elastic analysis with knockdown factors for buckling, combined loads, and material variability. The Federal Aviation Administration (FAA) and European Union Aviation Safety Agency (EASA) structural substantiation requirements are calibrated for linear analysis methods, and nonlinear analysis may require additional justification and validation before acceptance. However, nonlinear FEA is widely used in aerospace for detailed stress analysis of fastened joints, bearing analysis, and composite damage progression, where linear analysis is insufficient to capture the physical behavior.
Frequently Asked Questions
Q: When should I use nonlinear FEA instead of linear FEA?
Use nonlinear FEA when any of the three linear assumptions are violated: when maximum displacement exceeds 5 to 10 percent of the smallest characteristic structural dimension (geometric nonlinearity), when any location in the structure reaches or exceeds the material yield strength under service loads (material nonlinearity), or when contact conditions between parts change during loading, parts may separate or slide, or frictional forces are significant at interfaces (contact nonlinearity). As a practical shortcut: run a linear analysis first, check maximum displacement ratio and von Mises stress, and escalate to nonlinear only if the thresholds are exceeded.
Q: What is the difference between geometric and material nonlinearity in FEA?
Geometric nonlinearity (also called large displacement nonlinearity) occurs when the structure’s deformation is large enough that its stiffness changes as it deforms. The equilibrium equations must be evaluated in the deformed configuration rather than the original. Material nonlinearity occurs when the stress-strain relationship is not linear, which happens when the material yields plastically (metals above yield), when the material is inherently nonlinear even at small strains (rubber, elastomers, biological tissues), or when time-dependent effects like creep are present. Both types can occur simultaneously, and many real-world problems require both to be addressed together.
Q: How much slower is nonlinear FEA compared to linear FEA?
Nonlinear FEA is typically 5 to 50 times slower than an equivalent linear analysis for geometric and material nonlinearity alone, and 10 to 100 times slower when contact nonlinearity is included. The increase is due to the need for incremental loading (typically 10 to 100 load steps) and iterative solution within each step using the Newton-Raphson method, each iteration requiring a full matrix refactorization. Fully nonlinear analyses with all three nonlinearity types (geometry, material, contact) can require hundreds of matrix factorizations compared to the single factorization in a linear analysis.
Q: What is contact nonlinearity in FEA and when does it apply?
Contact nonlinearity occurs whenever the contact state between surfaces in an assembly changes during loading: two surfaces that are initially in contact may separate, initially separated surfaces may come into contact, or sliding occurs at a frictional interface. Any assembled mechanical structure where parts interact has potential contact nonlinearity: bolted joints, press-fit interfaces, bearing contacts, snap-fit joints, rubber seals against housings, and gear tooth contacts. Contact is computationally the most challenging nonlinearity type because the active contact set changes discontinuously during the solution, causing numerical instability that requires specialized contact algorithms to manage.
Q: What is the Newton-Raphson method in nonlinear FEA?
Newton-Raphson (N-R) is the iterative solution algorithm used in most implicit nonlinear FEA solvers. At each load increment, N-R iterates to find a displacement increment that reduces the residual (the difference between internal forces generated by the current displacement state and the externally applied forces) below a specified convergence tolerance. Each iteration requires computing the tangent stiffness matrix (which has changed from the previous iteration due to geometric or material state changes), factoring it, and computing a new displacement estimate. Convergence is achieved when the residual norm falls below the tolerance, typically 0.1 to 1 percent of the applied force norm. N-R convergence can fail if the load step is too large or if the physical system is genuinely unstable.
Q: Can I use linear FEA for buckling analysis?
Linear eigenvalue buckling analysis (Euler buckling) is valid for predicting the critical load at which a structure first loses stability, but it provides no information about the post-buckling behavior. If the structure is expected to carry load beyond the buckling point (shell buckling in aerospace structures, for example), or if imperfections significantly affect the buckling load (which they do for most real structures), nonlinear geometric analysis including initial geometric imperfections is required. The standard approach is to use linear eigenvalue buckling to identify the critical load and mode shape, then use the mode shape as an initial geometric imperfection for a nonlinear post-buckling analysis to determine the actual structural response.
Q: What does convergence failure in nonlinear FEA mean?
Convergence failure in nonlinear FEA means that the Newton-Raphson iterative solver could not find a displacement increment that reduces the force residual below the convergence tolerance within the maximum allowed number of iterations. This can indicate physical instability (the structure has reached a snap-through point, plastic collapse limit, or contact instability), model setup problems (unconstrained rigid body motion, poorly defined contact, material model extrapolated beyond its calibrated range), or numerical issues (increment size too large, mesh quality problems in high-strain regions). Troubleshooting begins by identifying the load level at which convergence fails and examining the deformed shape at the last converged increment to determine whether the failure is physical or numerical.
Conclusion:
The correct answer to when to use linear versus nonlinear FEA is not a preference for one approach over the other. It is a question answered by assessing three specific physical conditions against three specific thresholds, and choosing the analysis type that correctly captures the physical behavior that will govern the structural response under the loading conditions of interest.
Linear FEA is not a compromise or a shortcut when it is appropriate: it is the theoretically correct analysis for structures where all three linearity assumptions hold. The computational efficiency, the certainty of convergence, and the straightforward interpretation of results are properties of a correctly applied analysis, not concessions made to save time. Nonlinear FEA is not universally more accurate than linear FEA: it is more accurate for problems where the physical nonlinearity is present and significant, and it is equivalent to or marginally different from linear FEA for problems where it is not.
The engineering skill is the assessment: using the pre-check workflow to identify which thresholds are exceeded, understanding which type of nonlinearity dominates, selecting the appropriate nonlinear analysis type and solver approach, and knowing how to diagnose convergence problems when they arise. These skills are what separate an FEA practitioner who gets correct results efficiently from one who either misses real nonlinear behavior or spends excessive resources on nonlinear analyses for inherently linear problems.
Continue building your simulation and engineering analysis knowledge with our guides on CAD modeling best practices, design for assembly, multi-body modeling techniques, and master models for complex engineering programs.