Real-Time Elastic Body Simulation via Position-Based Dynamics

Resolve Bottlenecks,
Find Innovative Solutions
Generate Solutions

Solution Overview

Problem

Conventional methods for simulating the dynamics of elastic bodies, particularly in real-time scenarios, face challenges in achieving realistic and stable results due to issues with stiffness and energy balance, especially in conditions of large deformations and non-linear effects.

Innovation Solution

The introduction of analytic equations and algorithms for a position-based dynamics method, XPBD, which models elastic bodies as three-dimensional finite elements, allowing for efficient real-time simulation by computing energy and deformation using non-linear interpolation functions and material parameters, and incorporating FEM-based constraint functions to handle strain energy and compliance.

Engineering Contradictions & Design Principles

VSEngineering Contradiction Analysis

1Measurement precision

If conventional FEM methods are used for accurate simulation, then simulation accuracy is improved, but real-time performance deteriorates

Engineering Contradiction:
Improvesimulation accuracyVSAvoidreal-time performance
Core Design Contradiction:
Measurement precisionVSProductivity

Solution Approach 1:

The elastic body is segmented into discrete finite elements (tetrahedra, hexahedra, etc.) with defined nodes and interior points. This segmentation allows the complex continuous problem to be divided into manageable discrete elements that can be processed efficiently in real-time while maintaining accuracy through proper element formulation.

Inventive Principle:
Principle #1Segmentation

Solution Approach 2:

The patent transforms the traditional FEM approach by changing key parameters: using a position-based formulation with direct position updates instead of iterative force-based solving, employing analytic solutions for element constraints, and using a simplified time integration scheme. These parameter changes enable real-time performance while preserving accuracy.

Inventive Principle:
Principle #35Parameter changes

2Productivity

If position-based dynamics is used for real-time performance, then computational speed is improved, but simulation realism deteriorates in large deformations

Engineering Contradiction:
Improvecomputational speedVSAvoidsimulation realism
Core Design Contradiction:
ProductivityVSMeasurement precision

Solution Approach 1:

The patent substitutes the traditional force-based mechanical system with a position-based system. Instead of computing forces, accelerations, and velocities through iterative solvers, the method directly computes position corrections based on constraint violations. This substitution maintains computational speed while improving realism through analytic constraint satisfaction.

Inventive Principle:
Principle #28Mechanics substitution (Replace mechanical system)

Solution Approach 2:

The patent performs preliminary computation of element constraints and their gradients before the main simulation loop. By pre-computing the constraint functions and their derivatives for each finite element, the method enables fast real-time updates without sacrificing accuracy in large deformation scenarios.

Inventive Principle:
Principle #10Preliminary action

3Manufacturing precision

If iteration count is increased for constraint satisfaction, then constraint accuracy is improved, but computational efficiency deteriorates

Engineering Contradiction:
Improveconstraint accuracyVSAvoidcomputational efficiency
Core Design Contradiction:
Manufacturing precisionVSProductivity

Solution Approach 1:

The constraint functions and their gradients are pre-computed for each finite element before the simulation begins. This preliminary action eliminates the need for repeated iterative computations during the simulation, achieving high constraint accuracy in a single pass while maintaining computational efficiency.

Inventive Principle:
Principle #10Preliminary action

Solution Approach 2:

The patent uses analytic copies of the constraint satisfaction equations that can be evaluated directly without iteration. By deriving closed-form solutions for the constraint corrections based on pre-computed gradients, the method achieves exact constraint satisfaction in a single computational step rather than requiring multiple iterations.

Inventive Principle:
Principle #26Copying

4Stability of the object's composition

If time step is decreased for stability, then simulation stability is improved, but computational cost increases

Engineering Contradiction:
Improvesimulation stabilityVSAvoidcomputational cost
Core Design Contradiction:
Stability of the object's compositionVSLoss of time

Solution Approach 1:

The patent changes the time integration parameter by using a simplified explicit scheme with direct position updates. This parameter change allows for larger stable time steps compared to traditional implicit FEM methods, reducing computational cost while maintaining simulation stability through the position-based formulation's inherent numerical properties.

Inventive Principle:
Principle #35Parameter changes

Data Source

PatentUS20230061175A1Real-Time Simulation of Elastic Body
Publication Date: 2023.03.02 MIMYK MEDICAL SIMULATIONS PTE LTD
  • US20230061175A1 patent drawing
  • US20230061175A1 patent drawing
  • US20230061175A1 patent drawing

AI summary

An elastic body is modelled as a plurality of models of particles and models of geometric elements. Each particle model has a position attribute. Each geometric element has a boundary defined by two or more of the particles as nodes of the geometric element, nodes being points shared with neighbouring geometric elements. The interior area or volume of the geometric elements is modelled via non-linear interpolation functions that use the positions of the particles of the respective geometric elements as inputs to compute position and/or strain of interior points of the geometric elements. Energy is summed over the interior region of the element, the energy computation based on (a) position and/or stress of the geometric element computed by the non-linear interpolation functions and/or (b) non-linear material parameters that depend on position and/or strain at the interior points of the geometric elements. New positions of the particles are computed based minimizing total energy of the element.