Abstract This paper presents NURBS-OT (non-uniform rational B-splines—optimal transport), a new approach in the field of computer graphics and computer-aided design (CAD)/computer-aided manufacturing (CAM) for modeling complex free-form designs like aerodynamic and hydrodynamic structures, traditionally shaped by parametric curves such as Bézier, B-spline, and NURBS. Unlike prior models that used generative adversarial networks (GANs) involving large and complex parameter sets, our approach leverages a much lighter (0.37M versus 5.05M of BézierGAN), theoretically robust method by blending optimal transport with NURBS. This integration facilitates a more efficient generation of curvilinear designs. The efficacy of NURBS-OT has been validated through extensive testing on the University of Illinois Urbana-Champaign (UIUC) airfoil and superformula datasets, where it showed enhanced performance on various metrics. This demonstrates its ability to produce precise, realistic, and esthetically coherent designs, marking a significant advancement by merging classical geometrical techniques with modern deep learning.
more »
« less
Generalizations of non-uniform rational B-splines via decoupling of the weights: theory, software and applications
We introduce a new class of curves and surfaces by exploring multiple variations of non-uniform rational B-splines. These variations which are referred to as generalized non-uniform rational B-splines (GNURBS) serve as an alternative interactive shape design tool, and provide improved approximation abilities in certain applications. GNURBS are obtained by decoupling the weights associated with control points along different physical coordinates. This unexplored idea brings the possibility of treating the weights as additional degrees of freedoms. It will be seen that this proposed concept effectively improves the capability of NURBS, and circumvents its deficiencies in special applications. Further, it is proven that these new representations are merely disguised forms of classic NURBS, guaranteeing a strong theoretical foundation, and facilitating their utilization. A few numerical examples are presented which demonstrate superior approximation results of GNURBS compared to NURBS in both cases of smooth and non-smooth fields. Finally, in order to better demonstrate the behavior and abilities of GNURBS in comparison to NURBS, an interactive MATLAB toolbox has been developed and introduced.
more »
« less
- Award ID(s):
- 1661597
- PAR ID:
- 10167227
- Date Published:
- Journal Name:
- Engineering with computers
- Volume:
- 177
- ISSN:
- 0177-0667
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Heart valve disease (HVD), a significant cardiovascular complication, is one of the leading global causes of morbidity and mortality. Treatment for HVD often involves medical devices such as bioprosthetic valves. However, the design and optimization of these devices require a thorough understanding of their biomechanical and hemodynamic interactions with patient-specific anatomical structures. Parametric procedural geometry has become a powerful tool in enhancing the efficiency and accuracy of design optimization for such devices, allowing researchers to systematically explore a wide range of possible configurations. In this work, we present a robust framework for parametric and procedural modeling of stented bioprosthetic heart valves and patient-specific aortic geometries. The framework employs non-uniform rational B-splines (NURBS)-based geometric parameterization, enabling precise control over key anatomical and design variables. By enabling a modular and expandable workflow, the framework supports iterative optimization of valve designs to achieve improved hemodynamic performance and durability. We demonstrate its applicability through simulations on bioprosthetic aortic valves, highlighting the impact of geometric parameters on valve function and their potential for personalized device design. By coupling parametric geometry with computational tools, this framework offers researchers and engineers a streamlined pathway toward innovative and patient-specific cardiovascular solutions.more » « less
-
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
-
We present a theoretical and empirical study of the gradient dynamics of overparam- eterized shallow ReLU networks with one-dimensional input, solving least-squares interpolation. We show that the gradient dynamics of such networks are determined by the gradient flow in a non-redundant parameterization of the network function. We examine the principal qualitative features of this gradient flow. In particular, we determine conditions for two learning regimes: kernel and adaptive, which depend both on the relative magnitude of initialization of weights in different layers and the asymptotic behavior of initialization coefficients in the limit of large network widths. We show that learning in the kernel regime yields smooth interpolants, minimizing curvature, and reduces to cubic splines for uniform initializations. Learning in the adaptive regime favors instead linear splines, where knots cluster adaptively at the sample points.more » « less
An official website of the United States government

