A Deep Learning Approach to Shape Optimization Problems for Flexoelectric Materials Using the Isogeometric Finite Element Method

A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization problems.To improve the fitting ability of the neural network,we use the idea of pre-training to determine the structure of the neural network and combine different optimizers for training.The isogeometric analysis-finite element method(IGA-FEM)is used to discretize the flexural theoretical formulas and obtain samples,which helps ANN to build a proxy model from the model shape to the target value.The effectiveness of the proposed method is verified through two numerical examples of parameter optimization and one numerical example of shape optimization.

Isogeometric Analysis of Hyperelastic Material Characteristics for Calcified Aortic Valve

This study explores the implementation of computed tomography(CT)reconstruction and simulation techniques for patient-specific valves,aiming to dissect the mechanical attributes of calcified valves within transcatheter heart valve replacement(TAVR)procedures.In order to facilitate this exploration,it derives pertinent formulas for 3D multi-material isogeometric hyperelastic analysis based on Hounsfield unit(HU)values,thereby unlocking foundational capabilities for isogeometric analysis in calcified aortic valves.A series of uniaxial and biaxial tensile tests is executed to obtain an accurate constitutive model for calcified active valves.To mitigate discretization errors,methodologies for reconstructing volumetric parametric models,integrating both geometric and material attributes,are introduced.Applying these analytical formulas,constitutive models,and precise analytical models to isogeometric analyses of calcified valves,the research ascertains their close alignment with experimental results through the close fit in displacement-stress curves,compellingly validating the accuracy and reliability of the method.This study presents a step-by-step approach to analyzing themechanical characteristics of patient-specific valves obtained fromCT images,holding significant clinical implications and assisting in the selection of treatment strategies and surgical intervention approaches in TAVR procedures.

A Hybrid Parallel Strategy for Isogeometric Topology Optimization via CPU/GPU Heterogeneous Computing

This paper aims to solve large-scale and complex isogeometric topology optimization problems that consumesignificant computational resources. A novel isogeometric topology optimization method with a hybrid parallelstrategy of CPU/GPU is proposed, while the hybrid parallel strategies for stiffness matrix assembly, equationsolving, sensitivity analysis, and design variable update are discussed in detail. To ensure the high efficiency ofCPU/GPU computing, a workload balancing strategy is presented for optimally distributing the workload betweenCPU and GPU. To illustrate the advantages of the proposedmethod, three benchmark examples are tested to verifythe hybrid parallel strategy in this paper. The results show that the efficiency of the hybrid method is faster thanserial CPU and parallel GPU, while the speedups can be up to two orders of magnitude.

Quantifying Uncertainty in Dielectric Solids’ Mechanical Properties Using Isogeometric Analysis and Conditional Generative Adversarial Networks

Accurate quantification of the uncertainty in the mechanical characteristics of dielectric solids is crucial for advancing their application in high-precision technological domains,necessitating the development of robust com-putational methods.This paper introduces a Conditional Generation Adversarial Network Isogeometric Analysis(CGAN-IGA)to assess the uncertainty of dielectric solids’mechanical characteristics.IGA is utilized for the precise computation of electric potentials in dielectric,piezoelectric,and flexoelectric materials,leveraging its advantage of integrating seamlessly with Computer-Aided Design(CAD)models to maintain exact geometrical fidelity.The CGAN method is highly efficient in generating models for piezoelectric and flexoelectric materials,specifically adapting to targeted design requirements and constraints.Then,the CGAN-IGA is adopted to calculate the electric potential of optimum models with different parameters to accelerate uncertainty quantification processes.The accuracy and feasibility of this method are verified through numerical experiments presented herein.

Application of Isogeometric Analysis Method in Three-Dimensional Gear Contact Analysis

Gears are pivotal in mechanical drives,and gear contact analysis is a typically difficult problem to solve.Emerging isogeometric analysis(IGA)methods have developed new ideas to solve this problem.In this paper,a threedimensional body parametric gear model of IGA is established,and a theoretical formula is derived to realize single-tooth contact analysis.Results were benchmarked against those obtained from commercial software utilizing the finite element analysis(FEA)method to validate the accuracy of our approach.Our findings indicate that the IGA-based contact algorithmsuccessfullymet theHertz contact test.When juxtaposed with the FEA approach,the IGAmethod demonstrated fewer node degrees of freedomand reduced computational units,all whilemaintaining comparable accuracy.Notably,the IGA method appeared to exhibit consistency in analysis accuracy irrespective of computational unit density,and also significantlymitigated non-physical oscillations in contact stress across the tooth width.This underscores the prowess of IGA in contact analysis.In conclusion,IGA emerges as a potent tool for addressing contact analysis challenges and holds significant promise for 3D gear modeling,simulation,and optimization of various mechanical components.

Full-Scale Isogeometric Topology Optimization of Cellular Structures Based on Kirchhoff-Love Shells

Cellular thin-shell structures are widely applied in ultralightweight designs due to their high bearing capacity and strength-to-weight ratio.In this paper,a full-scale isogeometric topology optimization(ITO)method based on Kirchhoff-Love shells for designing cellular tshin-shell structures with excellent damage tolerance ability is proposed.This method utilizes high-order continuous nonuniform rational B-splines(NURBS)as basis functions for Kirchhoff-Love shell elements.The geometric and analysis models of thin shells are unified by isogeometric analysis(IGA)to avoid geometric approximation error and improve computational accuracy.The topological configurations of thin-shell structures are described by constructing the effective density field on the controlmesh.Local volume constraints are imposed in the proximity of each control point to obtain bone-like cellular structures.To facilitate numerical implementation,the p-norm function is used to aggregate local volume constraints into an equivalent global constraint.Several numerical examples are provided to demonstrate the effectiveness of the proposed method.After simulation and comparative analysis,the results indicate that the cellular thin-shell structures optimized by the proposed method exhibit great load-carrying behavior and high damage robustness.

New Perspective to Isogeometric Analysis:Solving Isogeometric Analysis Problem by Fitting Load Function

Isogeometric analysis(IGA)is introduced to establish the direct link between computer-aided design and analysis.It is commonly implemented by Galerkin formulations(isogeometric Galerkin,IGA-G)through the use of nonuniform rational B-splines(NURBS)basis functions for geometric design and analysis.Another promising approach,isogeometric collocation(IGA-C),working directly with the strong form of the partial differential equation(PDE)over the physical domain defined by NURBS geometry,calculates the derivatives of the numerical solution at the chosen collocation points.In a typical IGA,the knot vector of the NURBS numerical solution is only determined by the physical domain.A new perspective on the IGAmethod is proposed in this study to improve the accuracy and convergence of the solution.Solving the PDE with IGA can be regarded as fitting the load function defined on the NURBS geometry(right-hand side)with derivatives of the NURBS numerical solution(left-hand side).Moreover,the design of the knot vector has a close relationship to theNURBS functions to be fitted in the area of data fitting in geometric design.Therefore,the detected feature points of the load function are integrated into the initial knot vector of the physical domainto construct thenewknot vector of thenumerical solution.Then,they are connected seamlessly with the IGA-C framework for its great potential combining the accuracy and smoothness merits with the computational efficiency,which we call isogeometric collocation by fitting load function(IGACL).In numerical experiments,we implement our method to solve 1D,2D,and 3D PDEs and demonstrate the improvement in accuracy by comparing it with the standard IGA-C method.We also verify the superiority in the accuracy of our knot selection scheme when employed in the IGA-G method,which we call isogeometric Galerkin by fitting load function(IGA-GL).

Parameterization Transfer for a Planar Computational Domain in Isogeometric Analysis

In this paper,we propose a parameterization transfer algorithm for planar domains bounded by B-spline curves,where the shapes of the planar domains are similar.The domain geometries are considered to be similar if their simplified skeletons have the same structures.One domain we call source domain,and it is parameterized using multi-patch B-spline surfaces.The resulting parameterization is C1 continuous in the regular region and G1 continuous around singular points regardless of whether the parameterization of the source domain is C1/G1 continuous or not.In this algorithm,boundary control points of the source domain are extracted from its parameterization as sequential points,and we establish a correspondence between sequential boundary control points of the source domain and the target boundary through discrete sampling and fitting.Transfer of the parametrization satisfies C1/G1 continuity under discrete harmonic mapping with continuous constraints.The new algorithm has a lower calculation cost than a decomposition-based parameterization with a high-quality parameterization result.We demonstrate that the result of the parameterization transfer in this paper can be applied in isogeometric analysis.Moreover,because of the consistency of the parameterization for the two models,this method can be applied in many other geometry processing algorithms,such as morphing and deformation.

Explicit Isogeometric Topology Optimization Method with Suitably Graded Truncated Hierarchical B-Spline

This work puts forward an explicit isogeometric topology optimization(ITO)method using moving morphable components(MMC),which takes the suitably graded truncated hierarchical B-Spline based isogeometric analysis as the solver of physical unknown(SGTHB-ITO-MMC).By applying properly basis graded constraints to the hierarchical mesh of truncated hierarchical B-splines(THB),the convergence and robustness of the SGTHB-ITOMMC are simultaneously improved and the tiny holes occurred in optimized structure are eliminated,due to the improved accuracy around the explicit structural boundaries.Moreover,an efficient computational method is developed for the topological description functions(TDF)ofMMC under the admissible hierarchicalmesh,which consists of reducing the dimensionality strategy for design space and the locally computing strategy for hierarchical mesh.We apply the above SGTHB-ITO-MMC with improved efficiency to a series of 2D and 3Dcompliance design problems.The numerical results show that the proposed SGTHB-ITO-MMC method outperforms the traditional THB-ITO-MMCmethod in terms of convergence rate and efficiency.Therefore,the proposed SGTHB-ITO-MMC is an effective way of solving topology optimization(TO)problems.

Multi-Patch Black-White Topology Optimization in Isogeometric Analysis

Topological optimization plays a guiding role in the conceptual design process.This paper conducts research on structural topology optimization algorithm within the framework of isogeometric analysis.For multi-component structures,the Nitsche’smethod is used to glue differentmeshes to performisogeometricmulti-patch analysis.The discrete variable topology optimization algorithm based on integer programming is adopted in order to obtain clear boundaries for topology optimization.The sensitivity filtering method based on the Helmholtz equation is employed for averaging of curved elements’sensitivities.In addition,a simple averaging method along coupling interfaces is proposed in order to ensure the material distribution across coupling areas is reasonably smooth.Finally,the performance of the algorithm is demonstrated by numerical examples,and the effectiveness of the algorithm is verified by comparing it with the results obtained by single-patch and ABAQUS cases.

Topology Optimization of Sound-Absorbing Materials for Two-Dimensional Acoustic Problems Using Isogeometric Boundary Element Method

In this work,an acoustic topology optimizationmethod for structural surface design covered by porous materials is proposed.The analysis of acoustic problems is performed using the isogeometric boundary elementmethod.Taking the element density of porousmaterials as the design variable,the volume of porousmaterials as the constraint,and the minimum sound pressure or maximum scattered sound power as the design goal,the topology optimization is carried out by solid isotropic material with penalization(SIMP)method.To get a limpid 0–1 distribution,a smoothing Heaviside-like function is proposed.To obtain the gradient value of the objective function,a sensitivity analysis method based on the adjoint variable method(AVM)is proposed.To find the optimal solution,the optimization problems are solved by the method of moving asymptotes(MMA)based on gradient information.Numerical examples verify the effectiveness of the proposed topology optimization method in the optimization process of two-dimensional acoustic problems.Furthermore,the optimal distribution of sound-absorbingmaterials is highly frequency-dependent and usually needs to be performed within a frequency band.

An Isogeometric Cloth Simulation Based on Fast Projection Method

A novel continuum-based fast projection scheme is proposed for cloth simulation.Cloth geometry is described by NURBS,and the dynamic response is modeled by a displacement-only Kirchhoff-Love shell element formulated directly on NURBS geometry.The fast projection method,which solves strain limiting as a constrained Lagrange problem,is extended to the continuum version.Numerical examples are studied to demonstrate the performance of the current scheme.The proposed approach can be applied to grids of arbitrary topology and can eliminate unrealistic over-stretching efficiently if compared to spring-based methodologies.

Isogeometric Analysis of Longitudinal Displacement of a Simplified Tunnel Model Based on Elastic Foundation Beam

Serious uneven settlement of the tunnel may directly cause safety problems.At this stage,the deformation of the tunnel is predicted and analyzed mainly by numerical simulation,while the commonly used finite element method(FEM)uses low-order continuous elements.Therefore,the accuracy of tunnel settlement prediction is not enough.In this paper,a method is proposed to study the vertical deformation of the tunnel by using the combination of isogeometric analysis(IGA)and Bézier extraction operator.Compared with the traditional IGA method,this method can be easily integrated into the existing FEM framework,and ensure the same accuracy.A numerical example of an elastic foundation beam subjected to uniformly distributed load and an engineering example of an equivalent elastic foundation beamof the tunnel are given.The results show that the solution of the IGA method is closer to the theoretical solution of the initial-parameter method than the FEM,and the accuracy and reliability of the proposedmodel are verified.Moreover,it not only provides some theoretical support for the longitudinal design of the tunnel,but also provides a new way for the application and popularization of IGA in tunnel engineering.

Interpolating Isogeometric Boundary Node Method and Isogeometric Boundary Element Method Based on Parameter Space

In this paper,general interpolating isogeometric boundary node method(IIBNM)and isogeometric boundary element method(IBEM)based on parameter space are proposed for 2D elasticity problems.In both methods,the integral cells and elements are defined in parameter space,which can reproduce the geometry exactly at all the stages.In IIBNM,the improved interpolating moving leastsquare method(IIMLS)is applied for field approximation and the shape functions have the delta function property.The Lagrangian basis functions are used for field approximation in IBEM.Thus,the boundary conditions can be imposed directly in both methods.The shape functions are defined in 1D parameter space and no curve length needs to be computed.Besides,most methods for the treatment of the singular integrals in the boundary element method can be applied in IIBNM and IBEM directly.Numerical examples have demonstrated the accuracy of the proposed methods.

Isogeometric Collocation:A Mixed Displacement-Pressure Method for Nearly Incompressible Elasticity Dedicated to Professor Karl Stark Pister for his 95th birthday

We investigate primal and mixed u−p isogeometric collocation methods for application to nearly-incompressible isotropic elasticity.The primal method employs Navier’s equations in terms of the displacement unknowns,and the mixed method employs both displacement and pressure unknowns.As benchmarks for what might be considered acceptable accuracy,we employ constant-pressure Abaqus finite elements that are widely used in engineering applications.As a basis of comparisons,we present results for compressible elasticity.All the methods were completely satisfactory for the compressible case.However,results for low-degree primal methods exhibited displacement locking and in general deteriorated in the nearly-incompressible case.The results for the mixed methods behaved very well for two of the problems we studied,achieving levels of accuracy very similar to those for the compressible case.The third problem,which we consider a“torture test”presented a more complex story for the mixed methods in the nearly-incompressible case.

A Comprehensive Review of Isogeometric Topology Optimization:Methods,Applications and Prospects

Topology Optimization(TO)is a powerful numerical technique to determine the optimal material layout in a design domain,which has accepted considerable developments in recent years.The classic Finite Element Method(FEM)is applied to compute the unknown structural responses in TO.However,several numerical deficiencies of the FEM significantly influence the effectiveness and efficiency of TO.In order to eliminate the negative influence of the FEM on TO,IsoGeometric Analysis(IGA)has become a promising alternative due to its unique feature that the Computer-Aided Design(CAD)model and Computer-Aided Engineering(CAE)model can be unified into a same mathematical model.In the paper,the main intention is to provide a comprehensive overview for the developments of Isogeometric Topology Optimization(ITO)in methods and applications.Finally,some prospects for the developments of ITO in the future are also presented.

Structural Design Optimization Using Isogeometric Analysis: A Comprehensive Review

Isogeometric analysis(IGA),an approach that integrates CAE into conventional CAD design tools,has been used in structural optimization for 10 years,with plenty of excellent research results.This paper provides a comprehensive review on isogeometric shape and topology optimization,with a brief coverage of size optimization.For isogeometric shape optimization,attention is focused on the parametrization methods,mesh updating schemes and shape sensitivity analyses.Some interesting observations,e.g.the popularity of using direct(differential)method for shape sensitivity analysis and the possibility of developing a large scale,seamlessly integrated analysis-design platform,are discussed in the framework of isogeometric shape optimization.For isogeometric topology optimization(ITO),we discuss different types of ITOs,e.g.ITO using SIMP(Solid Isotropic Material with Penalization)method,ITO using level set method,ITO using moving morphable com-ponents(MMC),ITO with phase field model,etc.,their technical details and applications such as the spline filter,multi-resolution approach,multi-material problems and stress con-strained problems.In addition to the review in the last 10 years,the current developmental trend of isogeometric structural optimization is discussed.

Data-Driven Structural Design Optimization for Petal-Shaped Auxetics Using Isogeometric Analysis

Focusing on the structural optimization of auxetic materials using data-driven methods,a back-propagation neural network(BPNN)based design framework is developed for petal-shaped auxetics using isogeometric analysis.Adopting a NURBSbased parametric modelling scheme with a small number of design variables,the highly nonlinear relation between the input geometry variables and the effective material properties is obtained using BPNN-based fitting method,and demonstrated in this work to give high accuracy and efficiency.Such BPNN-based fitting functions also enable an easy analytical sensitivity analysis,in contrast to the generally complex procedures of typical shape and size sensitivity approaches.

Isogeometric analysis of free-form Timoshenko curved beams including the nonlinear effects of large deformations

In the present paper, the isogeometric analysis（IGA） of free-form planar curved beams is formulated based on the nonlinear Timoshenko beam theory to investigate the large deformation of beams with variable curvature. Based on the isoparametric concept, the shape functions of the field variables（displacement and rotation） in a finite element analysis are considered to be the same as the non-uniform rational basis spline（NURBS） basis functions defining the geometry. The validity of the presented formulation is tested in five case studies covering a wide range of engineering curved structures including from straight and constant curvature to variable curvature beams. The nonlinear deformation results obtained by the presented method are compared to well-established benchmark examples and also compared to the results of linear and nonlinear finite element analyses. As the nonlinear load-deflection behavior of Timoshenko beams is the main topic of this article, the results strongly show the applicability of the IGA method to the large deformation analysis of free-form curved beams. Finally, it is interesting to notice that, until very recently, the large deformations analysis of free-form Timoshenko curved beams has not been considered in IGA by researchers.

Three-Dimensional Isogeometric Analysis of Flexoelectricity with MATLAB Implementation

Flexoelectricity is a general electromechanical phenomenon where the electric polarization exhibits a linear dependency to the gradient of mechanical strain and vice versa.The truncated pyramid compression test is among the most common setups to estimate the flexoelectric effect.We present a three-dimensional isogeometric formulation of flexoelectricity with its MATLAB implementation for a truncated pyramid setup.Besides educational purposes,this paper presents a precise computational model to illustrate how the localization of strain gradients around pyramidal boundary shapes contributes in generation of electrical energy.The MATLAB code is supposed to help learners in the Isogeometric Analysis and Finite Elements Methods community to learn how to solve a fully coupled problem,which requires higher order approximations,numerically.The complete MATLAB code which is available as source code distributed under a BSD-style license,is provided in the part of Supplementary Materials of the paper.