An official website of the United States government
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
This content will become publicly available on February 1, 2027
Enhanced MPM framework with multipatch isogeometric analysis for geotechnical applications
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
- PAR ID:
- 10656655
- Publisher / Repository:
- Elsevier
- Date Published:
- Journal Name:
- Computers and Geotechnics
- Volume:
- 190
- Issue:
- C
- ISSN:
- 0266-352X
- Page Range / eLocation ID:
- 107743
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Abstract Continuous‐assumed‐strain (CAS) elements were recently introduced (Casquero and Golestanian. Comput Methods Appl Mech Eng. 2022; 399:115354.) to remove the membrane locking present in quadratic ‐continuous NURBS‐based discretizations of linear plane curved Kirchhoff rods. In this work, we generalize CAS elements to remove shear and membrane locking from quadratic NURBS‐based discretizations of linear plane curved Timoshenko rods. CAS elements are an assumed strain treatment that interpolates the shear and membrane strains at the knots using linear Lagrange polynomials. Consequently, the inter‐element continuity of the shear and membrane strains is maintained. The numerical experiments considered in this work show that CAS elements excise the spurious oscillations in shear and membrane forces caused by shear and membrane locking. Furthermore, when using CAS elements with either full or reduced integration, the convergence of displacements, rotations, and stress resultants is independent of the slenderness ratio up to while the convergence is highly dependent on the slenderness ratio when using NURBS elements. We apply the locking treatment of CAS elements to quadratic ‐continuous NURBS and the resulting element type is named discontinuous‐assumed‐strain (DAS) elements. Comparisons among CAS and DAS elements show that once locking is properly removed, continuity across element boundaries results in higher accuracy than continuity across element boundaries. Lastly, CAS elements result in a simple numerical scheme that does not add any significant computational burden in comparison with the locking‐prone NURBS‐based discretization of the Galerkin method.more » « less
-
Abstract Soft materials are of major interest for biomechanics applications due to their high deformability and susceptibility to experience damage events under different loading scenarios. The present study is concerned with modelling damage evolution processes in these nonlinear materials whose structural responses are prone to locking when low-order kinematic interpolation is employed in the context of nonlinear Finite Element schemes. For this reason, a pair of gradient-enhanced continuum damage schemes are proposed with the aim of tackling mechanical failure problems in applications that exhibit shear and volumetric locking. In particular, we present the consistent formulation and the assessment of the corresponding performance of (i) a mixed displacement-enhanced assumed strain employing a total Lagrangian formulation, and (ii) a three-field mixed displacement-pressure-Jacobian formulation. The novel and formulations are consistently derived and numerically implemented, providing a satisfactory agreement with respect to built-in elements handling the treatment of shear and volumetric locking, respectively, in conjunction to the modelling damage phenomena via the use of a penalty-based gradient-enhanced formulation. This performance is examined via several numerical applications. Furthermore, the final example justifies the need for a formulation combining both mixed FE approaches to simulate problems encompassing both locking issues (shear and volumetric locking), which can be performed using a combination of the and herein proposed.more » « less
-
Abstract This article presents a novel derivation for the governing equations of geometrically curved and twisted three-dimensional Timoshenko beams. The kinematic model of the beam was derived rigorously by adopting a parametric description of the axis of the beam, using the local Frenet–Serret reference system, and introducing the constraint of the beam cross ection planarity into the classical, first-order strain versus displacement relations for Cauchy’s continua. The resulting beam kinematic model includes a multiplicative term consisting of the inverse of the Jacobian of the beam axis curve. This term is not included in classical beam formulations available in the literature; its contribution vanishes exactly for straight beams and is negligible only for curved and twisted beams with slender geometry. Furthermore, to simplify the description of complex beam geometries, the governing equations were derived with reference to a generic position of the beam axis within the beam cross section. Finally, this study pursued the numerical implementation of the curved beam formulation within the conceptual framework of isogeometric analysis, which allows the exact description of the beam geometry. This avoids stress locking issues and the corresponding convergence problems encountered when classical straight beam finite elements are used to discretize the geometry of curved and twisted beams. Finally, this article presents the solution of several numerical examples to demonstrate the accuracy and effectiveness of the proposed theoretical formulation and numerical implementation.more » « less
