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 transcathete...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 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...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.展开更多
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 parallelstr...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.展开更多
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 ba...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.展开更多
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 thre...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.展开更多
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 the...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.展开更多
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-compone...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.展开更多
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.T...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.展开更多
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 d...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(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 nonuni...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).展开更多
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 me...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.展开更多
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 th...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.展开更多
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 c...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.展开更多
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,an...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.展开更多
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...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.展开更多
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.A...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(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 compr...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.展开更多
A general and new explicit isogeometric topology optimisation approach with moving morphable voids(MMV)is proposed.In this approach,a novel multiresolution scheme with two distinct discretisation levels is developed t...A general and new explicit isogeometric topology optimisation approach with moving morphable voids(MMV)is proposed.In this approach,a novel multiresolution scheme with two distinct discretisation levels is developed to obtain high-resolution designs with a relatively low computational cost.Ersatz material model based on Greville abscissae collocation scheme is utilised to represent both the Young’s modulus of the material and the density field.Two benchmark examples are tested to illustrate the effectiveness of the proposed method.Numerical results show that high-resolution designs can be obtained with relatively low computational cost,and the optimisation can be significantly improved without introducing additional DOFs.展开更多
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 amo...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.展开更多
This paper presents a novel approximating method to construct highprecision single-patch representation of B-spline surface from a multi-patch representation for isogeometric applications.In isogeometric analysis,mult...This paper presents a novel approximating method to construct highprecision single-patch representation of B-spline surface from a multi-patch representation for isogeometric applications.In isogeometric analysis,multi-patch structure is not easy to achieve high continuity between neighboring patches which will reduce the advantage of isogeometric analysis in a sense.The proposed method can achieve high continuity at surface stitching region with low geometric error,and this technique exploits constructing the approximate surface with several control points are from original surfaces,which guarantees the local feature of the surface can be well-preserved with high precision.With the proposed approximating method,isogeometric analysis results using the new single-patch can be obtained efficiently compared with the original multi-patch structure.Several examples are presented to illustrate the effectiveness,accuracy and efficiency of the proposed method.展开更多
基金supported by the Natural Science Foundation of China(Project Nos.52075340 and 61972011)the Shanghai Special Research Project on Aging Population and Maternal and Child Health(Project No.2020YJZX0106).
文摘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.
基金supported by a Major Research Project in Higher Education Institutions in Henan Province,with Project Number 23A560015.
文摘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.
基金the National Key R&D Program of China(2020YFB1708300)the National Natural Science Foundation of China(52005192)the Project of Ministry of Industry and Information Technology(TC210804R-3).
文摘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.
基金supported by the National Key R&D Program of China(Grant Number 2020YFB1708300)China National Postdoctoral Program for Innovative Talents(Grant Number BX20220124)+1 种基金China Postdoctoral Science Foundation(Grant Number 2022M710055)the New Cornerstone Science Foundation through the XPLORER PRIZE,the Knowledge Innovation Program of Wuhan-Shuguang,the Young Top-Notch Talent Cultivation Program of Hubei Province and the Taihu Lake Innovation Fund for Future Technology(Grant Number HUST:2023-B-7).
文摘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.
基金support provided by the National Nature Science Foundation of China (Grant Nos.52075340,51875360)Project of Science and Technology Commission of Shanghai Municipality (No.19060502300).
文摘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.
基金supported by the National Natural Science Foundation of China(Grant Nos.62072148 and U22A2033)the National Key R&D Program of China(Grant Nos.2022YFB3303000 and 2020YFB1709402)+2 种基金the Zhejiang Provincial Science and Technology Program in China(Grant No.2021C01108)the NSFC-Zhejiang Joint Fund for the Integration of Industrialization and Informatization(Grant No.U1909210)the Fundamental Research Funds for the Provincial Universities of Zhejiang(Grant No.490 GK219909299001-028).
文摘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.
基金supported by the Fundamental Research Funds for the Cen-tral Universities(No.JUSRP12038)the Natural Science Foundation of Jiangsu Province(No.BK20200611)the National Natural Science Foundation of China(No.12102146).
文摘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.
基金sponsored by Natural Science Foundation of Henan under Grant No.222300420498.
文摘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.
基金Chao Zheng thanks the support from Sichuan Science and Technology Program[Grant No.2021JDRC0007].
文摘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.
基金supported by the National Natural Science Foundation of China under Grant Nos.61872316,62272406,61932018the National Key R&D Plan of China under Grant No.2020YFB1708900.
文摘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).
基金support fromthe National Natural Science Foundation of China (52079128).
文摘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.
基金supported by the National Key R&D Program of China (2020YFB1708300)the Project funded by the China Postdoctoral Science Foundation (2021M701310).
文摘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.
基金The research for this paper was supported by(1)the National Natural Science Foundation of China(Grants Nos.51708429,51708428)the Open Projects Foundation(Grant No.2017-04-GF)of State Key Laboratory for Health and Safety of Bridge Structures+1 种基金Wuhan Institute of Technology Science Found(Grant No.K201734)the science and technology projects of Wuhan Urban and Rural Construction Bureau(Grants Nos.201831,201919).
文摘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.
基金FF and LDL gratefully acknowledge the financial support of the German Research Foundation(DFG)within the DFG Priority Program SPP 1748“Reliable Simulation Techniques in Solid Mechanics”.AR has been partially supported by the MIUR-PRIN project XFAST-SIMS(No.20173C478 N).
文摘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.
基金Supported by National Key R&D Program of China(Grant No.2020YFB1708304)and Fundamental Research Funds for the Central Universities of the Huazhong University of Science and Technology(Grant No.5003123021)and the Program for HUST Academic Frontier Youth Team(Grant No.2017QYTD04).
文摘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.
基金National Natural Science Foundation of China(Grant Nos.51705158 and 51805174)the Fundamental Research Funds for the Central Universities(Grant Nos.2018MS45 and 2019MS059)。
文摘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.
基金This work was supported by National Natural Science Foundation of China(51705158)the Fundamental Research Funds for the Central Universities(2018MS45)Open Funds of National Engineering Research Center of Near-Net-Shape Forming for Metallic Materials(2018005).
文摘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.
基金National Natural Science Foundation of China under Grant Nos.51675525 and 11725211.
文摘A general and new explicit isogeometric topology optimisation approach with moving morphable voids(MMV)is proposed.In this approach,a novel multiresolution scheme with two distinct discretisation levels is developed to obtain high-resolution designs with a relatively low computational cost.Ersatz material model based on Greville abscissae collocation scheme is utilised to represent both the Young’s modulus of the material and the density field.Two benchmark examples are tested to illustrate the effectiveness of the proposed method.Numerical results show that high-resolution designs can be obtained with relatively low computational cost,and the optimisation can be significantly improved without introducing additional DOFs.
基金Hamid Ghasemi acknowledge the support of the Mechanical Engineering department at Arak University of Technology.Xiaoying Zhuang gratefully acknowledge the financial support by European Research Council for COTOFLEXI project(802205)Harold Park acknowledges the support of the Mechanical Engineering department at Boston University.Timon Rabczuk gratefully acknowledge financial support by the 2019 Foreign Experts Plan of Hebei Province.
文摘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.
基金This research was supported by the National Nature Science Foundation of China under Grant Nos.61602138,61772163 and 61761136010the NSFC-Zhejiang Joint Fund for the Integration of Industrialization and Informatization(Grant No.U1909210)Zhejiang Provincial Science and Technology Program in China(2018C01030).
文摘This paper presents a novel approximating method to construct highprecision single-patch representation of B-spline surface from a multi-patch representation for isogeometric applications.In isogeometric analysis,multi-patch structure is not easy to achieve high continuity between neighboring patches which will reduce the advantage of isogeometric analysis in a sense.The proposed method can achieve high continuity at surface stitching region with low geometric error,and this technique exploits constructing the approximate surface with several control points are from original surfaces,which guarantees the local feature of the surface can be well-preserved with high precision.With the proposed approximating method,isogeometric analysis results using the new single-patch can be obtained efficiently compared with the original multi-patch structure.Several examples are presented to illustrate the effectiveness,accuracy and efficiency of the proposed method.
