An official website of the United States government
Summary This study presents a thermo‐hydro‐mechanical (THM) model of unsaturated soils using isogeometric analysis (IGA). The framework employs Bézier extraction to connect IGA to the conventional finite element analysis (FEA), featuring the current study as one of the first attempts to develop an IGA‐FEA framework for solving THM problems in unsaturated soils. IGA offers higher levels of interelement continuity making it an attractive method for solving highly nonlinear problems. The governing equations of linear momentum, mass, and energy balance are coupled based on the averaging procedure within the hybrid mixture theory. The Drucker‐Prager yield surface is used to limit the modified effective stress where the model follows small strain, quasi‐static loading conditions. Temperature dependency of the surface tension is implemented in the soil‐water retention curve. Nonuniform rational B‐splines (NURBS) basis functions are used in the standard Galerkin method and weak formulations of the balance equations. Displacement, capillary pressure, gas pressure, and temperature are four independent quantities that are approximated by NURBS in spatial discretization. The framework is used to simulate strain localization in an undrained dense sand subjected to plane strain biaxial compression under different temperatures and displacement velocities. Results show that an increase in the displacement rate leads to reduction in the equivalent plastic strain while an increase in the temperature leads to an increase in the equivalent plastic strain. The findings suggest that the proposed IGA‐based framework offers a viable alternative for solving THM problems in unsaturated soils.
more »
« less
Enhancing Isogeometric Analysis With NURBS-Based Synthesis
Abstract Isogeometric analysis (IGA) is a computational technique that integrates computer-aided design (CAD) with finite element analysis (FEA) by employing the same basis functions for both geometry representation and solution approximation. While IGA offers numerous advantages, such as improved accuracy and efficiency, it also presents several challenges related to geometric modeling. Some of these challenges include accurately representing complex geometries with NURBS (Non-Uniform Rational B-Splines) or other basis functions used in IGA and generating high-quality meshes that conform to the complex geometry represented by NURBS curves/surfaces. This paper introduces an analytical framework to provide a more efficient and theoretically grounded method for generating curvilinear configurations and its analytical solution in IGA, bridging the gap between generated data and its physical representations. This innovative approach is distinguished by integrating the NURBS parameterization in curve generation and providing a corresponding framework to achieve a broader and more accurate explanation of meshes and properties, especially constructing new coordinates and calculating the physical displacements under these conditions. Our model enables the analytical understanding of complex curves from the UIUC airfoil and superformula datasets, demonstrating a deeper dive into simulations. This study signifies a pivotal juncture wherein machine-learning-based complex geometrical formulations are synergistically combined with actual isogeometric analysis.
more »
« less
- Award ID(s):
- 2245299
- PAR ID:
- 10593465
- Publisher / Repository:
- American Society of Mechanical Engineers
- Date Published:
- ISBN:
- 978-0-7918-8835-3
- Format(s):
- Medium: X
- Location:
- Washington, DC, USA
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Achieving stable stress solutions at large strains using the Material Point Method (MPM) is challenging due to the accumulation of errors associated with geometry discretization, cell-crossing noise, and volumetric locking. Several simplified attempts exist in the literature to mitigate these errors, including higher-order frameworks. However, the stability of the MPM solution in such frameworks has been limited to simple geometries and the single-phase formulation (i.e., neglecting pore fluid). Although never explored, multipatch isogeometric analysis offers desirable qualities to simulate complex geometries while mitigating errors in the MPM. The degree of required high-order spatial integration has also never been investigated to infer a minimum limit for the stability of the stress solution in MPM. This paper presents a general-purpose numerical framework for simulating stable stresses in porous media, capturing both near incompressibility and multiphase interactions. First, the numerical framework is presented considering Non-Uniform Rational B-splines (NURBS) to perform isogeometric analysis (IGA) in MPM. Additionally, a volumetric strain smoothing algorithm is used to alleviate errors associated with volumetric locking. Second, the manifestation of cell-crossing errors is assessed via a series of problems with orders ranging from linear to cubic interpolation functions. Third, the use of NURBS is investigated and verified for problems with circular geometries. Finally, multipatch analysis is deployed to simulate plane strain and 3D penetration in soils, considering nearly incompressible elastoplastic (total stress) analysis and fully-coupled hydro-mechanical (effective stress) analysis. The stability of the solution is also analyzed for different constitutive models. From the results, it can be concluded that the framework using cubic interpolation functions with strain smoothing is the most convenient, presenting stable stress solutions for a broad range of multiphase geotechnical applications.more » « less
-
Abstract The superior accuracy isogeometric analysis (IGA) brought to computations in fluid and solid mechanics has been yielding higher fidelity in computational aerodynamics. The increased accuracy we achieve with the IGA is in the flow solution, in representing the problem geometry, and, when we use the IGA basis functions also in time in a space–time (ST) framework, in representing the motion of solid surfaces. It is of course as part of a set of methods that the IGA has been very effective in computational aerodynamics, including complex-geometry aerodynamics. The set of methods we have been using can be categorized into those that serve as a core method, those that increase the accuracy, and those that widen the application range. The core methods are the residual-based variational multiscale (VMS), ST-VMS and arbitrary Lagrangian–Eulerian VMS methods. The IGA and ST-IGA are examples of the methods that increase the accuracy. The complex-geometry IGA mesh generation method is an example of the methods that widen the application range. The ST Topology Change method is another example of that. We provide an overview of these methods for IGA-based computational aerodynamics and present examples of the computations performed. In computational flow analysis with moving solid surfaces and contact between the solid surfaces, it is a challenge to represent the boundary layers with an accuracy attributed to moving-mesh methods and represent the contact without leaving a mesh protection gap.more » « less
-
Neurons exhibit complex geometry in their branched networks of neurites which is essential to the function of individual neuron but also brings challenges to transport a wide variety of essential materials throughout their neurite networks for their survival and function. While numerical methods like isogeometric analysis (IGA) have been used for modeling the material transport process via solving partial differential equations (PDEs), they require long computation time and huge computation resources to ensure accurate geometry representation and solution, thus limit their biomedical application. Here we present a graph neural network (GNN)-based deep learning model to learn the IGA-based material transport simulation and provide fast material concentration prediction within neurite networks of any topology. Given input boundary conditions and geometry configurations, the well-trained model can predict the dynamical concentration change during the transport process with an average error less than 10% and 120∼330 times faster compared to IGA simulations. The effectiveness of the proposed model is demonstrated within several complex neurite networks.more » « less
-
null (Ed.)Abstract Neurons exhibit complex geometry in their branched networks of neurites which is essential to the function of individual neuron but also brings challenges to transport a wide variety of essential materials throughout their neurite networks for their survival and function. While numerical methods like isogeometric analysis (IGA) have been used for modeling the material transport process via solving partial differential equations (PDEs), they require long computation time and huge computation resources to ensure accurate geometry representation and solution, thus limit their biomedical application. Here we present a graph neural network (GNN)-based deep learning model to learn the IGA-based material transport simulation and provide fast material concentration prediction within neurite networks of any topology. Given input boundary conditions and geometry configurations, the well-trained model can predict the dynamical concentration change during the transport process with an average error less than 10% and $$120 \sim 330$$ 120 ∼ 330 times faster compared to IGA simulations. The effectiveness of the proposed model is demonstrated within several complex neurite networks.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Conference Paper:
- https://doi.org/10.1115/DETC2024-142195
