International Journal for Numerical Methods in Engineering最新文献

{"title":"Surface Structured Quadrilateral Mesh Generation Based on Topology Consistent-Preserved Patch Segmentation","authors":"Haoxuan Zhang,&nbsp;Haisheng Li,&nbsp;Xiaoqun Wu,&nbsp;Nan Li","doi":"10.1002/nme.7644","DOIUrl":"https://doi.org/10.1002/nme.7644","url":null,"abstract":"<div>\u0000 \u0000 <p>Surface-structured mesh generation is an important part of the Computational Fluid Dynamics (CFD) preprocessing stage. The traditional method cannot automatically divide the 3D surface topology in complex structures. Thus, we propose a surface-structured quadrilateral mesh generation method based on topology-consistent-preserved patch segmentation. The core idea is to segment the complex 3D model into several simple parts according to mesh quality and map each part to the 2D parametric domain based on the conformal parameterization method. Then, we utilize pattern-based topology partitioning to divide the parametric domain into multiple quadrilateral subdomains, facilitating the generation of 2D structured quadrilateral meshes. By using the inverse mapping algorithm based on barycentric weights, the generated 2D structured mesh is inversely mapped back to the 3D space. Finally, we splice each part accurately according to the structured mesh distribution. Experimental results show that our proposed method can generate higher-quality structured quadrilateral meshes than previous methods without losing the mesh topology of the original model.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 1","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-12-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143119051","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Time-Domain Spectral BFS Plate Element With Lobatto Basis for Wave Propagation Analysis","authors":"Hela Ambati,&nbsp;Sascha Eisenträger,&nbsp;Santosh Kapuria","doi":"10.1002/nme.7617","DOIUrl":"https://doi.org/10.1002/nme.7617","url":null,"abstract":"<div>\u0000 \u0000 <p>A computationally efficient spectral Kirchhoff plate element is presented for time-domain analysis of wave propagation at high frequencies in thin isotropic plates. It employs a <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <msup>\u0000 <mrow>\u0000 <mi>C</mi>\u0000 </mrow>\u0000 <mrow>\u0000 <mn>1</mn>\u0000 </mrow>\u0000 </msup>\u0000 </mrow>\u0000 <annotation>$$ {C}^1 $$</annotation>\u0000 </semantics></math>-continuous spectral interpolation based on the modified bi-Hermite polynomials using the Gauss–Lobatto–Legendre (GLL) points as a basis with selective collocation of rotational and twisting degrees of freedom (DOFs) at element edge and corner nodes. The lowest order version of the proposed element reduces to the classical Bogner–Fox–Schmit (BFS) element for Kirchhoff plates. The GLL basis allows diagonalisation of the mass matrix using the nodal quadrature technique, which enhances the computational efficiency. The numerical properties of the proposed element are comprehensively evaluated, including the conditioning of the system matrices. Moreover, the effect of employing different numerical integration schemes and nodal sets is examined in both static and free vibration analyses. The effectiveness of the proposed element in wave propagation problems is evaluated by comparing its performance to the converged solutions achieved using the BFS element with a very fine mesh. Results demonstrate that the current element, without and even with mass matrix diagonalisation delivers exceptional accuracy while also exhibiting faster convergence and enhanced computational efficiency than the existing Kirchhoff plate elements.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 3","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-12-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143363005","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"An Electro-Elastic Coupling Model for Piezoelectric Composites Based on the Voronoi Cell Finite Element Method","authors":"Huan Li,&nbsp;Nan Yang,&nbsp;Ran Guo","doi":"10.1002/nme.7631","DOIUrl":"https://doi.org/10.1002/nme.7631","url":null,"abstract":"<div>\u0000 \u0000 <p>The Voronoi cell finite element method (VCFEM) has successfully characterized the linear elastic behavior of the composites. This study is dedicated to develop an electro-elastic coupling VCFEM model mimicking the fully-coupled electro-elastic behavior of piezoelectric composites. For fiber-reinforced piezoelectric composites considering interfacial cracks, the interface traction reciprocity, the interface charge density reciprocity on bonded interfaces and the interface traction-free, the interface charge density-free on debonded interfaces are comprised in the new assumed stress and electric displacement hybrid variational functional. The new variational functional is derived on the base of the element multifield energy functionals. Independent stress/electric displacement fields are respectively assumed within the two-phase material domain. Several numerical examples considering perfectly-bonded interface and partially cracked interface were used to demonstrate the accuracy of the proposed method by comparing the piezoelectric Voronoi element model results with those obtained by ABAQUS. Then this model is used to study the effect of several microscopic details, such as the property ratio of fiber to matrix, volume fraction, interfacial crack length and polarization direction on macroscopic equivalent physical and mechanical properties, as well as local stress/electric displacement fields. It is clear that the proposed model is suitable for analyzing piezoelectric composites containing many microstructures with bonded interface or debonded interface.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 1","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-12-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142861666","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"An Improved Spectral Element Differential Method in Solving Nonlinear Thermoelastic Coupling Problems With Discontinuous Interfaces","authors":"Jianning Zhao,&nbsp;Dong Wei,&nbsp;Yuxi Wang,&nbsp;Donghuan Liu","doi":"10.1002/nme.7645","DOIUrl":"https://doi.org/10.1002/nme.7645","url":null,"abstract":"<div>\u0000 \u0000 <p>In this paper, the spectral element differential method (SEDM) is improved to solve the nonlinear thermoelastic coupling problems with interface thermal resistance and interface gap in composite structures. The utilization of both Lobatto and Chebyshev node sets in SEDM significantly enhances solution efficiency by replacing integration with direct differential. Moreover, for strongly nonlinear problems caused by thermal radiation, unknown terms are incorporated into the stiffness matrix, and the relaxation iteration technique is also employed, the convergence has been improved compared to traditional methods. Importantly, the element format of the SEDM for 3D problems with discontinuous interfaces is given specifically in this paper, and element-by-element loop assembly of stiffness matrices is realized. Numerical examples confirm the effectiveness of the present method in efficiently and accurately solving 2D and 3D problems with discontinuous interfaces. The present method not only achieves faster convergence than the traditional finite element method, but also attains higher accuracy with fewer degrees of freedom and shorter computational time. Compared to the spectral element method (SEM), the proposed method significantly reduces the computation time of the stiffness matrix. Furthermore, by employing the coupled SEM-SEDM approach, computational efficiency is enhanced while maintaining high precision in sensitive regions.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 1","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-12-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142861667","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Novel Implicit Integration Algorithms With Identical Second-Order Accuracy and Flexible Dissipation Control for First-Order Transient Problems","authors":"Jinze Li,&nbsp;Naigang Cui,&nbsp;Yijun Zhu,&nbsp;Kaiping Yu","doi":"10.1002/nme.7639","DOIUrl":"https://doi.org/10.1002/nme.7639","url":null,"abstract":"<div>\u0000 \u0000 <p>For the solutions of first-order transient problems with optimal efficiency, conventional single-step implicit methodologies, unaided by additional post-processing techniques, encounter limitations in concurrently attaining identical second-order accuracy and controllable numerical dissipation. Addressing this challenge, the present study contributes not only a comprehensive analytical framework for formulating implicit integration algorithms but also leverages the auxiliary variable and the composite sub-step technique to propose two distinct types of implicit integration algorithms. Each type is characterized by self-initiation, unconditional stability, identical second-order accuracy, controllable numerical dissipation, and zero-order overshoots. Recognizing that single-step implicit methods can be conceptualized as composite single-sub-step ones, both algorithm types achieve identical effective stiffness matrices, thereby reducing the computational effort for solving linear problems. The utilization of auxiliary variables endows the proposed single-step method with numerical attributes akin to established counterparts, while achieving identical second-order accuracy. Conversely, the adoption of the composite sub-step technique in the proposed two-sub-step methods surpasses the published algorithms by significantly reducing the error constants. This superiority persists when enforcing the same sub-step sizes and sub-step dissipation levels. Numerical simulations affirm that the proposed methods consistently outperform existing alternatives without incurring additional computational expenses.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 3","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-12-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143362997","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Reduced Order Modelling of Fully Coupled Electro-Mechanical Systems Through Invariant Manifolds With Applications to Microstructures","authors":"Attilio Frangi,&nbsp;Alessio Colombo,&nbsp;Alessandra Vizzaccaro,&nbsp;Cyril Touzé","doi":"10.1002/nme.7641","DOIUrl":"https://doi.org/10.1002/nme.7641","url":null,"abstract":"<p>This article presents the first application of the direct parametrisation method for invariant manifolds to a fully coupled multiphysics problem involving the nonlinear vibrations of deformable structures subjected to an electrostatic field. The formulation proposed is intended for model order reduction of electrostatically actuated resonating Micro-Electro-Mechanical Systems (MEMS). The continuous problem is first rewritten in a manner that can be directly handled by the parametrisation method, which relies upon automated asymptotic expansions. A new mixed fully Lagrangian formulation is thus proposed, which contains only explicit polynomial nonlinearities, which is then discretised in the framework of finite element procedures. Validation is performed on the classical parallel plate configuration, where different formulations using either the general framework or an approximation of the electrostatic field due to the geometric configuration selected are compared. Reduced-order models along these formulations are also compared to full-order simulations operated with a time integration approach. Numerical results show a remarkable performance both in terms of accuracy and the wealth of nonlinear effects that can be accounted for. In particular, the transition from hardening to softening behaviour of the primary resonance while increasing the constant voltage component of the electric actuation is recovered. Secondary resonances leading to superharmonic and parametric resonances are also investigated with the reduced-order model.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 3","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-12-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.7641","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143362372","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Topology Optimization of Elastoplastic Structure Based on Shakedown Strength","authors":"Songhua Huang,&nbsp;Lele Zhang,&nbsp;Geng Chen,&nbsp;Yugong Xu,&nbsp;Min Chen,&nbsp;Zhiyuan Liu,&nbsp;Eng Gee Lim","doi":"10.1002/nme.7627","DOIUrl":"https://doi.org/10.1002/nme.7627","url":null,"abstract":"<div>\u0000 \u0000 <p>The traditional approach to structural lightweight optimization design, which is based on the elastic limit rule, often results in a structure that exhibits either weight redundancy or strength redundancy to some extent. This study introduces a novel integration of shakedown analysis with structural topology optimization, departing from the conventional elastic limit rule. Shakedown analysis identifies a non-failure external load region beyond the elastic limit but below the plastic limit, independent of loading history. The proposed method, for the first time, accounts for the influence of self-equilibrium residual stress at the element level, redefining effective and ineffective elements in topology optimization. Shakedown total stress replaces elastic equivalent stress, offering a comprehensive measure. Utilizing Melan's lower bound theorem, a gradient-based topology optimization framework for shakedown analysis is developed, ensuring structures stay within the elastic–plastic range, preventing excessive plastic deformation. The approach, employing the moving asymptotes method after adjoint sensitivity analysis of shakedown total stress, is applied to a three-dimensional L-shaped bracket. Even with a remarkable 50% reduction in weight, the maximum total shakedown stress of the bracket reveals that it only increases by a modest 17.20% from its initial value. Moreover, compared to traditional topology optimization methods based on either elastic stress or stiffness, the proposed method based on total shakedown stress leads to a higher shakedown limit. Specifically, the configuration designed using the total shakedown stress exhibited increases of 2.01% and 9.82% in the shakedown limit compared to those obtained using stiffness and equivalent elastic stress, respectively. This suggests that the proposed method can effectively balance the trade-off between shakedown strength and structural stiffness, achieving a 2.01% rise in shakedown strength with only a 2.24% compromise in structural stiffness. These findings highlight the method's effectiveness and potential, emphasizing the benefit of redefining effective and ineffective elements using shakedown stress in topology optimization.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 3","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-12-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143362897","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A Novel Discrete Element Method for Smooth Polyhedrons and Its Application to Modeling Flows of Concave-Shaped Particles","authors":"Siqiang Wang,&nbsp;Qingwei Xu,&nbsp;Dongfang Liang,&nbsp;Shunying Ji","doi":"10.1002/nme.7628","DOIUrl":"https://doi.org/10.1002/nme.7628","url":null,"abstract":"<div>\u0000 \u0000 <p>The smooth polyhedral model has been commonly used to construct non-spherical particles with smooth surfaces, whereas it is mainly constrained to numerical simulations involving concave-shaped particles. This constraint arises from the limitations imposed by the contact algorithm. In this study, the contact detection between smooth polyhedrons is simplified to that between dilated triangular elements, and a discrete element method for concave polyhedral particles with smooth surfaces is developed. Subsequently, an automatic mesh simplification algorithm is established to enhance the computational efficiency without compromising accuracy. In validating the smooth polyhedral model, the simulation results of a hexahedron colliding with a plane are found to agree favorably with the experimental results. Then, the elastic collisions between the convex and concave particles are analyzed, and the total kinetic energy before and after the particle collision remains unchanged. Furthermore, the influences of particle morphology on the packing fraction, flow fluctuation, flow rate, mixing rate, velocity distribution, and system energy in hoppers and rotating drums are analyzed, revealing the underlying flow characteristics of concave polyhedral granular materials with smooth surfaces.</p>\u0000 </div>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 3","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-12-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143362899","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A Finite Interface Element for Modeling Inclusion-Matrix Interaction in Magneto-Active Elastomers","authors":"Will Klausler,&nbsp;Michael Kaliske","doi":"10.1002/nme.7625","DOIUrl":"https://doi.org/10.1002/nme.7625","url":null,"abstract":"<p>Magneto-active elastomers (MAEs) are an emerging smart composite material consisting of a compliant elastomer matrix filled with stiff, magnetizable inclusions. Meso-scale models, which treat the matrix and inclusions as discrete, are currently formulated with the assumption of perfect bonding at the interface. We apply a model for cohesive damage at the material interface and introduce a treatment for the transmission of the magnetic field in the gap. Examples demonstrate the concept and the difference in structural response between samples treated with the conventional and presented approaches.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 3","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-12-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.7625","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143362877","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Convergence of Damped Polarization Schemes for the FFT-Based Computational Homogenization of Inelastic Media With Pores","authors":"Elodie Donval,&nbsp;Matti Schneider","doi":"10.1002/nme.7632","DOIUrl":"https://doi.org/10.1002/nme.7632","url":null,"abstract":"<p>Porous microstructures represent a challenge for the convergence of FFT-based computational homogenization methods. In this contribution, we show that the damped Eyre–Milton iteration is linearly convergent for a class of nonlinear composites with a regular set of pores, provided the damping factor is chosen between zero and unity. First, we show that an abstract fixed-point method with non-expansive fixed-point operator and non-trivial damping converges linearly, provided the associated residual mapping satisfies a monotonicity condition on a closed subspace. Then, we transfer this result to the framework of polarization schemes and conclude the linear convergence of the damped Eyre–Milton scheme for porous materials. We present general arguments which apply to a class of nonlinear composites and mixed stress-strain loadings, as well. We show that the best contraction estimate is reached for a damping factor of <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mn>1</mn>\u0000 <mo>/</mo>\u0000 <mn>2</mn>\u0000 </mrow>\u0000 <annotation>$$ 1/2 $$</annotation>\u0000 </semantics></math>, that is, for the polarization scheme of Michel–Moulinec–Suquet, and derive the corresponding optimal reference material. Our results generalize the recent work of Sab and co-workers who showed that an adaptively damped Eyre–Milton scheme leads to linear convergence for a class of linear composites with pores. Finally, we report on computational experiments which support our findings.</p>","PeriodicalId":13699,"journal":{"name":"International Journal for Numerical Methods in Engineering","volume":"126 3","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-12-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nme.7632","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143362876","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}