{"title":"3D formulation of mono-symmetrical composite beams with deformable connection","authors":"Yassir Wardi, Pisey Keo, Mohammed Hjiaj","doi":"10.1016/j.finel.2024.104163","DOIUrl":"https://doi.org/10.1016/j.finel.2024.104163","url":null,"abstract":"<div><p>This paper deals with a 3D linear formulation for mono-symmetric composite beams with deformable connection, taking into account non-uniform torsion. To simplify the development of the analytical solution, it is assumed that the warping of each layer of the composite section has no contribution on the stress resultants of each layer. Therefore, the warping function obtained with the classical St-Venant beam theory can be used for each subsection. As a result, the variables associated to both connection shearing plans become uncoupled. Using the virtual work principle, the governing equations are derived, and solved in closed-form. Based on the analytical expressions of the displacement fields, the exact stiffness matrix of the composite beam is computed. In addition, a displacement-based formulation is suggested. Appropriate polynomial interpolation functions are selected to circumvent slip-locking phenomenon. It has been shown that the slip-locking can be avoided by using quadratic shape function for axial displacement interpolations, by providing an additional middle node in each layer. Four examples are investigated in this paper. The prediction as well as the performance of the proposed direct stiffness method, are compared against an existing solution from the literature. In addition, slip-locking problem is addressed and the performance of the displacement-based method against the exact formulation is evaluated. The influence of warping effects on the composite beam response is assessed. Finally, a parametric study is conducted to evaluate the influence of connection rigidity and the coupling of the displacement fields on slip distributions.</p></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"237 ","pages":"Article 104163"},"PeriodicalIF":3.1,"publicationDate":"2024-04-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140640835","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":"Multiphysics modeling of magnetoelectric composite disks by a 2D axisymmetric finite element approach","authors":"S. Karimi , H. Talleb","doi":"10.1016/j.finel.2024.104169","DOIUrl":"https://doi.org/10.1016/j.finel.2024.104169","url":null,"abstract":"<div><p>A 2D axisymmetric finite element multiphysics model is proposed to study magnetoelectric composite disks. This modeling approach includes a nonlinear magneto-elastic model to replicate the behavior of magnetostrictive materials under static conditions. Additionally, it offers a harmonic regime resolution that considers frequency dependence, including the implicit inclusion of eddy currents in the formulation, as well as electrical load. To validate the model, simulation results regarding the dependence of the static magnetic field and frequency are presented and compared with experimental measurements from literature.</p></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"237 ","pages":"Article 104169"},"PeriodicalIF":3.1,"publicationDate":"2024-04-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140644550","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":"Stress–displacement stabilized finite element analysis of thin structures using solid-shell elements, Part I: On the need of interpolating the stresses","authors":"A. Aguirre , R. Codina , J. Baiges","doi":"10.1016/j.finel.2024.104168","DOIUrl":"https://doi.org/10.1016/j.finel.2024.104168","url":null,"abstract":"<div><p>This work studies the solid-shell finite element approach to approximate thin structures using a stabilized mixed displacement–stress formulation based on the Variational Multiscale framework. The work is divided in two parts. In Part I, the numerical locking effects inherent to the solid-shell approach are characterized using a variety of benchmark problems in the infinitesimal strain approximation. In Part II, the results are extended to formulate the mixed approach in finite strain hyperelastic problems. In the present work, the stabilized mixed displacement–stress formulation is proven to be adequate to deal with all kinds of numerical locking. Additionally, a more comprehensive analysis of each individual type of numerical locking, how it is triggered and how it is overcome is also provided. The numerical locking usually occurs when parasitic strains overtake the system of equations through specific components of the stress tensor. To properly analyze them, the direction of each component of the stress tensor has been defined with respect to the shell directors. Therefore, it becomes necessary to formulate the solid-shell problem in curvilinear coordinates, allowing to give mechanical meaning to the stress components (shear, twisting, membrane and thickness stresses) independently of the global frame of reference. The conditions in which numerical locking is triggered as well as the stress tensor component responsible of correcting the locking behavior have been identified individually by characterizing the numerical response of a set of different benchmark problems.</p></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"236 ","pages":"Article 104168"},"PeriodicalIF":3.1,"publicationDate":"2024-04-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0168874X24000623/pdfft?md5=7c545daa21b7e5817ed48636c9015b4a&pid=1-s2.0-S0168874X24000623-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140631635","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":"Multiscale topology optimization of structures by using isogeometrical level set approach","authors":"Masoud Aminzadeh, Seyed Mehdi Tavakkoli","doi":"10.1016/j.finel.2024.104167","DOIUrl":"https://doi.org/10.1016/j.finel.2024.104167","url":null,"abstract":"<div><p>This study aims to optimize topology of structures at macro and micro scales, simultaneously, by using a level set method in an isogeometric analysis (IGA) framework. To achieve this, equilibrium and homogenization equations in the model are solved by IGA method. The level set functions are defined over a grid in parameter space of associating b-splines of the IGA model. Therefore, control net of the model and level set grid are separated and there is no need to refine the control net for having smooth boundaries. Sensitivity analyses for both scales are performed to calculate the velocity of boundary points and the level set functions are updated by solving reaction-diffusion equations. Finally, several 2D and 3D examples with different geometry and boundary conditions are provided to show performance and efficiency of the method. Obtained results show good agreement with examples in literature in terms of both topology and final value of objective function. Also, by using IGA level set method, smooth boundaries are achieved in the final topology of micro and macro structures.</p></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"235 ","pages":"Article 104167"},"PeriodicalIF":3.1,"publicationDate":"2024-04-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140549469","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}
Tiannan Hu , Yaguang Wang , Hao Li , Minghao Yu , Kozo Furuta , Kazuhiro Izui , Shinji Nishiwaki
{"title":"Topology optimization of coated structures infilled with multiple materials","authors":"Tiannan Hu , Yaguang Wang , Hao Li , Minghao Yu , Kozo Furuta , Kazuhiro Izui , Shinji Nishiwaki","doi":"10.1016/j.finel.2024.104165","DOIUrl":"https://doi.org/10.1016/j.finel.2024.104165","url":null,"abstract":"<div><p>A novel topology optimization method for the design of coated structures infilled with multiple materials is proposed in this paper, where a novel material interpolation model for the topology description is developed based on the ordered SIMP scheme. With the introduction of two special Heaviside projections into the two-step filtering and projection procedure, the external coating and the substrate region can be well identified by using several modified design variables. Then, the material distribution of the multi-material infilling is obtained by multiplying the infill identification field with the piece-wisely projected design variables and optimized via the mathematical programming algorithm under the ordered SIMP framework. Using an eroded density field and its original field, the uniform thickness of the external coating can be well controlled. The proposed approach for optimizing coated structures with multi-phase infill materials is easy to implement due to its implementation relying on those frequently-used filtering and projection operations. Besides, without introducing any additional design variables, the method developed in this paper inherits the advantages of the ordered SIMP method and has great calculation efficiency and stable iteration performance. With the consideration of several issues such as different coating thicknesses and different design parameters, several 2D numerical examples are studied to demonstrate the effectiveness of the proposed approach, as well as a 3D example. The optimization results illustrate that the method developed in this paper is effective for the design of coated structures infilled with multiple materials and the advantages of considering multiple infill materials is also validated.</p></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"235 ","pages":"Article 104165"},"PeriodicalIF":3.1,"publicationDate":"2024-04-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140540676","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":"Finite elements procedure for evaluating temperature distribution in polyurethane-coated wheels","authors":"M. Palmieri , F. Cianetti , C. Braccesi","doi":"10.1016/j.finel.2024.104164","DOIUrl":"https://doi.org/10.1016/j.finel.2024.104164","url":null,"abstract":"<div><p>Polymeric materials find extensive applications across various engineering sectors. Among these, a particularly critical application for these materials is in the field of roller coasters. The wheels are typically made with an aluminum hub and a dense polyurethane coating, which, being in contact with the track, endures dynamic loads at high speeds. Due to the viscoelastic behavior typical of polymeric materials, these loads induce overheating of the coating leading to rapid degradation of the wheel. This results in machine downtime and a significant waste of time and money. In this manuscript, a methodology for finite element thermal-structural analysis has been developed. This method allows for the rapid evaluation of temperatures reached during operational cycles if compared to classical coupled-field thermal-structural analysis. The proposed methodology proves to be useful in selecting the appropriate type of wheels during the design phase requiring short computational time. The study first involved the development of the methodology, followed by validation through a comparison of analysis results with data obtained from experimental tests conducted by the manufacturer.</p></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"236 ","pages":"Article 104164"},"PeriodicalIF":3.1,"publicationDate":"2024-04-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140539787","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}
Chien H. Thai , P.T. Hung , H. Nguyen-Xuan , P. Phung-Van
{"title":"A free vibration analysis of carbon nanotube reinforced magneto-electro-elastic nanoplates using nonlocal strain gradient theory","authors":"Chien H. Thai , P.T. Hung , H. Nguyen-Xuan , P. Phung-Van","doi":"10.1016/j.finel.2024.104154","DOIUrl":"https://doi.org/10.1016/j.finel.2024.104154","url":null,"abstract":"<div><p>This study presents a combination approach of the higher-order shear deformation theory, nonlocal strain gradient theory (NSGT) and isogeometric analysis (IGA) for the free vibration of carbon nanotube-reinforced (CNT) magneto-electro-elastic (MEE) nanoplates. To account size-dependent effects at the nanoscale, the classical theory model is extended with two additional scale parameters. However, this extended model necessitates at least the third derivative of the approximation function, which is incompatible with the standard finite element method. So, IGA with NURBS offers higher-order continuity through its basis functions, making it well-suited for this size-dependent model. To simplify computations, a power-law scheme is employed to represent the material properties. Various distribution types of carbon nanotubes (CNTs) including UD, FG-X, FG-O and FG-V are incorporated to investigate their effects on mechanical behaviors of CNT-MEE nanoplates. The governing equations of motion are derived in their weak form using the principle of extended virtual displacement and then solved by isogeometric analysis (IGA). The impact of the magnetic, electric and elastic fields on the coupling behaviors of CNT-MEE nanoplates are studied. Specially, parametric studies are conducted to analyze the influence of geometrical parameters, CNT distributions, CNT volume fraction, matrix volume fraction, electric voltage, magnetic potential, nonlocal and strain gradient parameters on the natural frequencies of the CNT-MEE nanoplates. Comparisons between the results obtained using NSGT and the classical theory reveal significant findings. The natural frequencies calculated by NSGT exhibit dependence on the relative values of the nonlocal and strain gradient parameters.</p></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"236 ","pages":"Article 104154"},"PeriodicalIF":3.1,"publicationDate":"2024-04-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140344042","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}
A. Martínez-Martínez , E. Nadal , D. Bonete , O. Allix , J.J. Ródenas
{"title":"Patient-specific numerical simulation of the bone healing process including implant materials and gait conditions","authors":"A. Martínez-Martínez , E. Nadal , D. Bonete , O. Allix , J.J. Ródenas","doi":"10.1016/j.finel.2024.104153","DOIUrl":"https://doi.org/10.1016/j.finel.2024.104153","url":null,"abstract":"<div><p>Numerical simulations in biomechanics, particularly in bone healing, present a cost-effective option compared to experiments that demand prolonged observations with human or with animal models. However, to define in-silico simulations of the bone healing process requires considering multiple factors, such as the implant design and patient’s characteristics. As a result, the current challenge is integrating different numerical methodologies to simulate bone healing, aiming to facilitate the emergence of innovative clinical treatments and new implant designs.</p><p>In this paper, we present a patient-specific numerical methodology to simulate the bone healing process, able to consider patient’s load conditions and bone density distribution provided by CT-scans. The main novelty is the combination of the Cartesian grid Finite Element Method (cgFEM) with a bone callus healing model, complemented by a load-condition optimisation scheme to relate implant materials and load conditions while ensuring successful healing outcome.</p><p>This numerical methodology creates a finite element model based on the patient’s medical image, serving as a virtual testing tool for investigating the influence of implant materials on gait pattern requirements to ensure an optimal healing outcome. In practice, a personalised bone fracture model was employed to evaluate four distinct implant materials: two conventional materials (stainless steel and titanium) and two bioabsorbable candidates (polylactic acid plastic (PLA) and magnesium). The results offer personalised optimal load conditions for each studied material, showcasing the potential of in-silico studies in minimising uncertainties associated with exploring new clinical treatments.</p></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"236 ","pages":"Article 104153"},"PeriodicalIF":3.1,"publicationDate":"2024-04-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0168874X24000477/pdfft?md5=cb369afac55c556e3007c98e1e342ad3&pid=1-s2.0-S0168874X24000477-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140344043","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":"Geometric design and performance analysis of a foldcore sandwich acoustic metastructure for tunable low-frequency sound absorption","authors":"Yao Chen , Zerui Shao , Jialong Wei , Jian Feng , Pooya Sareh","doi":"10.1016/j.finel.2024.104150","DOIUrl":"https://doi.org/10.1016/j.finel.2024.104150","url":null,"abstract":"<div><p>Acoustic metamaterial structures have received extensive attention for sound and vibration engineering applications from the scientific community in recent years. However, the real-life application of conventional acoustic metamaterial structures is frequently limited by fixed frequency bands and increased structural thicknesses in low-frequency noise reduction. In this study, we introduce an origami-based acoustic metamaterial structure that consists of a Miura-ori foldcore, along with a perforated and an unperforated panel. The proposed Miura-ori foldcore sandwich acoustic metastructure (MOF-SAM) exhibits adjustable low-frequency sound absorption capacities due to the foldability of the origami foldcore. Moreover, we employ numerical methods to investigate the sound absorption properties of the MOF-SAM, quantified by the sound absorption coefficient. The results indicate that the structure has a single absorption peak which is superior to that of acoustic structures composed of conventional honeycomb cores. The dissipation of acoustic energy is due to the structural vibrations of the metastructure and the losses in the folding process of the origami foldcore. The numerical results of this study show that the proposed sound absorption mechanism enables tunable low-frequency sound absorption. The geometric design and periodicity of the origami unit fragments offer multiple distinct absorption peaks and thus tunable acoustic performance. These findings of this study are expected to inspire novel designs for next-generation acoustic devices.</p></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"235 ","pages":"Article 104150"},"PeriodicalIF":3.1,"publicationDate":"2024-03-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140308572","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":"Polytopal templates for semi-continuous vectorial finite elements of arbitrary order on triangulations and tetrahedralizations","authors":"Adam Sky, Ingo Muench","doi":"10.1016/j.finel.2024.104155","DOIUrl":"https://doi.org/10.1016/j.finel.2024.104155","url":null,"abstract":"<div><p>The Hilbert spaces <span><math><mrow><mi>H</mi><mrow><mo>(</mo><mi>curl</mi><mo>)</mo></mrow></mrow></math></span> and <span><math><mrow><mi>H</mi><mrow><mo>(</mo><mi>div</mi><mo>)</mo></mrow></mrow></math></span> are employed in various variational problems formulated in the context of the de Rham complex in order to guarantee well-posedness. Seeing as the well-posedness follows automatically from the continuous setting to the discrete setting in the presence of commuting interpolants as per Fortin’s criterion, the construction of conforming subspaces becomes a crucial step in the formulation of stable numerical schemes. This work aims to introduce a novel, simple method of directly constructing semi-continuous vectorial base functions on the reference element via template vectors associated with the geometric polytopes of the element and an underlying <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span>-conforming polynomial subspace. The base functions are then mapped from the reference element to the element in the physical domain via consistent Piola transformations. The method is defined in such a way, that the underlying <span><math><msup><mrow><mi>H</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span>-conforming subspace can be chosen independently, thus allowing for constructions of arbitrary polynomial order. We prove a linearly independent construction of Nédélec elements of the first and second type, Brezzi–Douglas–Marini elements, and Raviart–Thomas elements on triangulations and tetrahedralizations. The application of the method is demonstrated with two examples in the relaxed micromorphic model.</p></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":"236 ","pages":"Article 104155"},"PeriodicalIF":3.1,"publicationDate":"2024-03-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0168874X24000490/pdfft?md5=454fed7b42151a56aa00cd211ecadcd3&pid=1-s2.0-S0168874X24000490-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140296272","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}