Stiffness-Matrix-Free Power Sweep for Dynamic Simulation Stability
Find Innovative SolutionsGenerate Solutions
Solution Overview
Problem
Existing finite element analysis (FEA) methods for dynamic motion simulations are computationally expensive and require conservative estimates for the Courant Stability Limit, which can lead to inefficient time step calculations, especially for poorly shaped finite elements and nonlinear materials.
Innovation Solution
A stiffness-matrix-free power sweep method that iteratively computes the highest natural frequency of each finite element in every time step, reusing computationally intensive components and utilizing hypoelastic material moduli to determine the Courant Stability Limit, allowing for precise time interval calculations without deriving new stability limit equations for new element and material types.
Engineering Contradictions & Design Principles
Engineering Contradiction Analysis
1Reliability
If conservative least upper bounds formulas are used to estimate natural frequency for every element in every time step, then the Courant Stability Limit can be determined to ensure simulation convergence, but the computational expense increases significantly
Solution Approach 1:
The patent segments the computational domain into finite elements and applies different frequency estimation strategies to different element types. By categorizing elements into groups with similar characteristics, the method computes natural frequencies more efficiently while maintaining the Courant Stability Limit requirements for each segment, thereby reducing overall computational expense while preserving simulation convergence.
Solution Approach 2:
The patent changes the approach from using conservative least upper bounds formulas to using more accurate but computationally efficient frequency estimation methods. By adjusting the estimation parameters and formulas based on element type and deformation characteristics, the method achieves better computational efficiency while maintaining the reliability needed for simulation convergence.
2Measurement precision
If the time step is reduced to satisfy the Courant Stability Limit, then simulation accuracy is improved, but the total simulation time increases
Solution Approach 1:
The patent implements dynamic time step adjustment based on the instantaneous state of the simulation. By monitoring element deformation, material stiffness changes, and frequency estimates in real-time, the method dynamically adapts the time step size to maintain Courant Stability Limit compliance while maximizing simulation efficiency. This allows larger time steps when conditions permit, reducing total simulation time while preserving accuracy when needed.
3Measurement precision
If new stability limit equations are derived for each new finite element type, then accuracy for that element type is improved, but device complexity increases
Solution Approach 1:
The patent develops a universal framework for stability limit calculation that works across multiple finite element types without requiring separate equations for each element. By creating a generalized approach that accommodates different element geometries, material models, and deformation states within a single computational structure, the method maintains accuracy for various element types while significantly reducing software complexity and ease of implementation.
Data Source
AI summary
A method, apparatus, and system provide the ability to simulate dynamic motion for a computerized model (of finite mesh elements). An element diagonal lumped mass matrix of the mesh, an estimate of a highest element eigenvector and eigenvalue of the mesh, and a kinematic state of the model are computed. Processing iterates until exceeding a time duration. Incremental strain and stress tensors, and hypo-elastic material constants are computed. Within the time duration iteration, eigenvalues are converged, a power-sweep stress field is computed from the strain field using the material constants; divergence of the power-sweep stress field is computed using the current gradient operator; and a power-sweep estimate is computed. Upon convergence, the stability limit is determined and utilized as a time interval for simulating the dynamic motion.