Hierarchical Shape Functions for MPM Geotechnical Stability
Find Innovative SolutionsGenerate Solutions
Solution Overview
Problem
The Material Point Method (MPM) in geotechnical engineering faces instability due to the use of higher-order Lagrangian shape functions, which lead to badly conditioned system matrices and unstable numerical solutions, particularly when material points only partially cover elements, and previous solutions like removing negative parts or using B-splines have compromised accuracy and efficiency.
Innovation Solution
The implementation of hierarchical shape functions, which are a composition of lower-order and higher-order polynomials, ensures strict non-negativity, preventing internal cancellations and maintaining a well-conditioned system matrix, along with compensation factors to address partition-of-unit and Kronecker delta property issues during interpolation and boundary condition calculations.
Engineering Contradictions & Design Principles
Engineering Contradiction Analysis
1Measurement precision
If higher-order Lagrangian shape functions are used in MPM, then the accuracy of numerical integration is improved, but the system matrix becomes badly conditioned or singular, leading to unstable numerical solutions
Solution Approach 1:
The patent changes the functional form of shape functions from Lagrangian to hierarchical B-spline basis functions. This parameter change in the mathematical representation maintains non-negativity while providing higher-order continuity, thereby improving integration accuracy without compromising solution stability.
Solution Approach 2:
The patent employs composite basis functions that combine the advantages of different function types. The hierarchical B-spline basis functions integrate the non-negativity property of piecewise polynomials with the smoothness of higher-order continuous functions, creating a composite mathematical tool that resolves the contradiction between accuracy and stability.
2Reliability
If negative parts of Lagrangian shape functions are removed, then the system matrix conditioning is improved, but the accuracy of numerical integration over an element is significantly reduced
Solution Approach 1:
Instead of removing negative parts, the patent changes the fundamental parameter of the shape function definition to hierarchical B-spline basis functions. These functions are inherently non-negative while maintaining higher-order polynomial characteristics, thus improving matrix conditioning without sacrificing integration accuracy.
3Reliability
If B-splines are used as basis functions, then the stability of numerical solutions is improved, but the computational efficiency is reduced and memory consumption increases
Solution Approach 1:
The patent applies segmentation by using hierarchical B-spline basis functions that are defined locally over element domains. This segmentation approach limits the coupling of shape functions to only adjacent elements, creating a sparse system matrix that reduces computational complexity and memory requirements compared to global B-spline formulations.
Solution Approach 2:
The hierarchical B-spline basis functions exhibit local quality by having compact support within and around adjacent elements. This local property ensures that only shape functions within an element and its direct neighbors are coupled, maintaining computational efficiency while providing the stability benefits of B-splines.
Data Source
AI summary
In one embodiment, material points are received that cover at least a portion of an element of a background mesh that describes a continuum of soil, rock and/or groundwater. MPM-based geotechnical analysis and simulation is conducted at least in part by performing a numerical integration over the material points to produce a system matrix and right-hand side vector. The numerical integration applies hierarchical shape functions to the material points. The MPM-based geotechnical analysis and simulation also may subtract out contributions of any lower-order polynomials from higher-order polynomials of the hierarchical shape functions when interpolating one or more state variables for the material points to the background mesh. The MPM-based geotechnical analysis and simulation also may subtract out contributions any lower-order polynomials from higher-order polynomials of the hierarchical shape functions when calculating one or more boundary conditions for the material points.