Finite Elements in Analysis and Design最新文献

{"title":"A pure Stokes approach for coupling fluid flow with porous media flow","authors":"Modesar Shakoor,&nbsp;Chung Hae Park","doi":"10.1016/j.finel.2023.104106","DOIUrl":"https://doi.org/10.1016/j.finel.2023.104106","url":null,"abstract":"<div><p>Most numerical approaches for coupling fluid flow with porous media<span><span> flow rely either on Stokes equations in the fluid part of the domain and Darcy’s law in the porous part, or on Brinkman’s equation. In both cases, difficulties arise at the boundary between the two parts because the equations used in the porous part are not of Stokes type. In this paper, an alternative to Darcy’s law is proposed for modeling </span>flows in porous media<span>. This alternative relies on equations of Stokes type where the permeability tensor is replaced by force and stress derivative tensors. Numerical procedures are presented to compute these tensors from simulations at pore scale. Simulations in domains containing both fluid and porous parts are finally conducted simply assuming continuity of velocity and pressure and hence without imposing any condition at the boundary between the two parts. Results show that the proposed method is accurate and hence a promising alternative to Darcy’s law for problems involving both fluid and porous parts.</span></span></p></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":null,"pages":null},"PeriodicalIF":3.1,"publicationDate":"2024-01-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139406536","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":"Interplay of liquid particles and interphases on the macroscopic elastic response of Liquid-filled composites","authors":"J. Sadeghi,&nbsp;F. Kamarei,&nbsp;T. Goudarzi","doi":"10.1016/j.finel.2023.104102","DOIUrl":"https://doi.org/10.1016/j.finel.2023.104102","url":null,"abstract":"<div><p><span><span>This paper deals with providing the effective elastic response<span> of three-phase composites consist of a matrix filled with a random suspension of liquid-filled capsules firmly bonded to the matrix in the realm of small deformation theory. The capsules shell (interphases) and the matrix are considered to be elastic solids and the liquid is considered ideal. For this purpose, the solution for dilute concentrations of the interphases and the liquid particles are derived analytically. Then, the dilute solution is passed to an </span></span>iterative technique to generate the non-dilute version of the solution. The non-dilute response is available in terms of a system of two nonlinear-coupled-ODE-initial-value problems. The results are confronted with the results of a 3D full-field FE analysis and </span>experimental data<span> and a good agreement is observed. The solution, provided here for polydisperse microstructures, allows for considering different thicknesses for the interphases and in turn can possibly account for different physical phenomena like the size and the surface tension effects. A comprehensive study on the effects of interphases (capsules shell) thickness and mechanical properties as well as the volume concentration of the capsules and the liquid particles on the overall elastic properties of three-phase liquid-filled composites are carried out. We show that although the addition of liquid particles to the matrix has deteriorating effects on the overall elastic properties of the composites, proper selection of the interphases can compensate for these negative effects.</span></p></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":null,"pages":null},"PeriodicalIF":3.1,"publicationDate":"2023-12-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139050442","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":"NURBS-enhanced finite element method (NEFEM) on quadrilateral meshes","authors":"Mattia Montanari ,&nbsp;Gian Maria Santi ,&nbsp;Ruben Sevilla ,&nbsp;Liverani Alfredo ,&nbsp;Nik Petrinic","doi":"10.1016/j.finel.2023.104099","DOIUrl":"https://doi.org/10.1016/j.finel.2023.104099","url":null,"abstract":"<div><p>This paper formulates quadrilateral elements for the NURBS-enhanced finite element method (NEFEM). The objective is to extend the application of NEFEM to problems where the use of quadrilateral elements is preferred. By leveraging a mapping, between reference and physical spaces, that encapsulates the exact boundary representation of the domain, a tight integration with computer aided design (CAD) systems is achieved. The contribution of this work is an enhanced quadrilateral finite element that incorporates the exact CAD geometry purely from the boundary representation (B-rep) from CAD and without the need for a whole volume representation (V-rep) as a NURBS entity. Numerical examples involving heat transfer and linear elastic problems are used to numerically demonstrate the optimal convergence properties of the method under mesh refinement.</p></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":null,"pages":null},"PeriodicalIF":3.1,"publicationDate":"2023-12-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0168874X23001920/pdfft?md5=aca1828fdc0b83332c707503b2866d87&pid=1-s2.0-S0168874X23001920-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138839088","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":"Method of matched sections as a beam-like approach for plate analysis","authors":"Igor Orynyak,&nbsp;Kirill Danylenko","doi":"10.1016/j.finel.2023.104103","DOIUrl":"10.1016/j.finel.2023.104103","url":null,"abstract":"<div><p>A <em>new numerical method in application to the plate problem is suggested. It starts from consideration of the rectangular elements, each operating by 6 beam-like parameters: four bending parameters (displacement, angle of rotation, bending moment, transverse force) and two rotation parameters (angle of rotation and twisting moment). So, contrary to the classical FEM approach, the conjugation between the adjacent elements occurs between the adjacent sections rather than in polygon vertexes (nodes), and this gives the name to the method – method of matched sections, MMS. Technically, 6 left-side and 6 lower-side parameters are 12 inlet parameters, while 6 upper-side and right-side parameters are 12 outlet parameters; each element is defined by 24 parameters. Outlet parameters are related to inlet ones by 12 matrix relations, which are derived from the approximate solution of each differential equation (equilibrium, physical, and geometrical equations) of the plate theory. The matrix relation between inlet and outlet parameters is written in a form suitable for applying the transfer matrix method. The numerical examples for the thin and Mindlin plates show the high efficiency and accuracy of the method. In particular, the results for the Mindlin plate for minimal thickness give the same results as the thin plate (no shear locking); the method is insensitive to cases when, for adjacent elements, the ratio of dimensions or ratio of rigidities (elastic constants) differ by several orders of magnitude. The application of the method for the thermal as well as for the vibration problem is considered. The possible extension of the method to any curved geometry is discussed, too.</em></p></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":null,"pages":null},"PeriodicalIF":3.1,"publicationDate":"2023-12-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0168874X23001968/pdfft?md5=1e67aa56e8df556c59be3ffd9dcaec2e&pid=1-s2.0-S0168874X23001968-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138714224","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":"An efficient method for the finite element analysis of shell structures by placing feature-fitted local shell meshes in a global shell mesh","authors":"Thuan Ho-Nguyen-Tan ,&nbsp;Hyun-Gyu Kim","doi":"10.1016/j.finel.2023.104101","DOIUrl":"10.1016/j.finel.2023.104101","url":null,"abstract":"<div><p><span>This paper presents an efficient method for the finite element analysis<span><span> of shell structures using feature-fitted local shell </span>meshes<span> that are placed in a global shell mesh. Feature-fitted local shell meshes are independently constructed to accurately represent the geometric features of shell structures. Non-matching interfaces between global and local shell meshes are connected by interface shell elements with an arbitrary number of nodes on the element sides. Assumed covariant strains are used to alleviate the shear and membrane </span></span></span>locking phenomena in interface shell elements. Numerical results show that more accurate solutions near the geometric features of shell structures can be obtained from feature-fitted shell meshes compared to those from conventional shell meshes composed of mixed-typed shell elements.</p></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":null,"pages":null},"PeriodicalIF":3.1,"publicationDate":"2023-12-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138679177","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":"Local refinement for the modeling of composite beam based on the partition of the unity method","authors":"P. Vidal,&nbsp;L. Gallimard,&nbsp;O. Polit","doi":"10.1016/j.finel.2023.104100","DOIUrl":"https://doi.org/10.1016/j.finel.2023.104100","url":null,"abstract":"<div><p><span>In this work, a refined model is superimposed on a simple one only in a region of interest to improve the accuracy of the modeling of composite beam structures. The purpose is to concentrate the computational effort in a localized zone without loss of local precision. Outside the region of interest is modelized using a simple cheap model. The present approach is based on the Partition Unity Method (PUM) to ensure the continuity condition on the displacement field. Two refined models based on the Sinus model are used to demonstrate the capability of the method. Each of the kinematics needs a suitable Finite Element (FE) </span>approximation<span>. The method is assessed on numerical examples involving both homogeneous and laminated beams. The results are compared with reference solutions and computations using a mono-model. The computational cost is also evaluated to show the interest of the method. The results are very promising.</span></p></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":null,"pages":null},"PeriodicalIF":3.1,"publicationDate":"2023-12-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138548898","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":"Simulating 3D printing on hydrogel inks: A finite element framework for predicting mechanical properties and scaffold deformation","authors":"M.C.P. Vila Pouca ,&nbsp;M.R.G. Cerqueira ,&nbsp;J.P.S. Ferreira ,&nbsp;R. Darabi ,&nbsp;N.A.G. Ramião ,&nbsp;R. Sobreiro-Almeida ,&nbsp;A.P.G. Castro ,&nbsp;P.R. Fernandes ,&nbsp;J.F. Mano ,&nbsp;RM Natal Jorge ,&nbsp;M.P.L. Parente","doi":"10.1016/j.finel.2023.104098","DOIUrl":"10.1016/j.finel.2023.104098","url":null,"abstract":"<div><h3>Background</h3><p>Difficulties during the wound healing process may result in scarring, chronic wounds and sepsis. A common tissue engineering strategy to solve these problems rely on the development of 3D hydrogel scaffolds that mimic the structure, stiffness, and biological proprieties of the target tissue. One of the most effective biofabrication techniques to precisely control spatial deposition, architecture and porosity of hydrogels is 3D printing technology. However, final architectures of 3D printed structures can be compromised if the printing properties are not adequately selected.</p></div><div><h3>Purpose</h3><p>Our main goal was to create a numerical framework able to predict the deformations that arise due to the 3D printing process of hydrogel scaffolds. Our secondary goal was to analyze if the overall mechanical properties of the 3D printed scaffolds were affected by these deformations.</p></div><div><h3>Methods</h3><p>We applied finite element analysis using ABAQUS finite element software to develop our numerical framework. The finite elements were added in a time sequence, simulating the material deposition. The bulk material was experimentally characterized and represented numerically by the user-defined subroutine UMAT. We tested the simulation by ‘printing’ a 5.0 × 5.0 × 0.8 alginate ink at 5 and 10 mm/s. Afterwards, both the final 3D printed scaffolds and a theoretical non-deformed configuration were subjected to a uniaxial compression of 10 % of the initial height, and differences between their overall mechanical properties were analyzed.</p></div><div><h3>Results</h3><p>The numerical framework captured the bending between the scaffold filaments and the compression of the bottom layers. On average, the scaffold printed at 5 mm/s deformed ∼6 % more, compared to the scaffold printed at 10 mm/s. However, in terms of overall mechanical properties, both showed similar behavior. This behavior, however, was highly nonlinear and significantly different from the theoretical, non-deformed scaffold, particularly in a small strains’ regime.</p></div><div><h3>Conclusions</h3><p>A numerical framework that can be used as a preliminary tool to define the printing velocity, sequence and geometry, minimizing the deformations during the 3D printing process was developed. This framework can help to minimize experimentation and consequently, material waste. We also saw that these deformations should not be neglected when predicting the mechanical behavior using finite element analysis, particularly for small strains application.</p></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":null,"pages":null},"PeriodicalIF":3.1,"publicationDate":"2023-12-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0168874X23001919/pdfft?md5=418ed2e4f4849bc5fc4222bed4059804&pid=1-s2.0-S0168874X23001919-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138492088","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":"ANN strategies for the stress–strain analysis of metallic materials: Modeling, database, supervised learning, validation and performance analysis","authors":"P.G. Marques Flávio,&nbsp;L.R. Cabral Muniz,&nbsp;T. Doca","doi":"10.1016/j.finel.2023.104097","DOIUrl":"https://doi.org/10.1016/j.finel.2023.104097","url":null,"abstract":"<div><p>Artificial neural networks<span><span> (ANN) are developed and employed to characterize a wide range of metallic materials. Focus is given to the evaluation of stress–strain behavior via sphere-to-flat indentation. Each ANN is trained using a supervised machine learning<span> procedure comprised of two steps: (i) generation of a training dataset via calibrated finite element model, and (ii) validation using </span></span>experimental data retrieved from indentations tests. The developed frameworks aim to establish a fast and low-cost tool for the assessment of loading conditions in industrial applications. The best proposed solution is able to predict stress–strain behavior with a quasi-instantaneous response and errors of less than 3%. Moreover, outputs are attained with minute costs (processing and memory bandwidth) when compared to finite element simulations.</span></p></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":null,"pages":null},"PeriodicalIF":3.1,"publicationDate":"2023-11-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138472384","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":"Computation of the dynamic scalar response of large two-dimensional periodic and symmetric structures by the wave finite element method","authors":"D. Duhamel","doi":"10.1016/j.finel.2023.104096","DOIUrl":"10.1016/j.finel.2023.104096","url":null,"abstract":"<div><p><span><span><span>In the past, the study of periodic media mainly focused on one-dimensional periodic structures (meaning periodic along one direction), on the one hand to determine the dispersion curves linking the frequencies to the wavenumbers and on the other hand to obtain the response of a structure to an external excitation, both for bounded or unbounded structures. In the latter case, effective approaches have been obtained, based on methods such as the Wave </span>Finite Element (WFE). Two-dimensional periodic media are more complex to analyse but dispersion curves can be obtained rather easily as in the one-dimensional case. Obtaining the </span>steady state response of two-dimensional periodic structures to time-harmonic excitations is much more difficult than for one-dimensional media and the results mainly concern infinite media. This work is about this last case of the steady state response of finite two-dimensional periodic structures to time-harmonic excitations by limiting oneself to structures described by a </span>scalar variable<span><span> (acoustic, thermal, membrane behaviour) and having symmetries compared to two orthogonal planes parallel to the edges of a substructure. Using the WFE for a rectangular substructure and imposing the wavenumber in one direction, we can numerically calculate the wavenumbers and mode shapes associated with propagation in the perpendicular direction. By building solutions with </span>null forces on parallel boundaries, we can decouple the waves in the two directions parallel to the sides of the rectangle. The solution of each of these two problems is obtained by a fast Fourier transformation giving the amplitudes associated with the waves. By summing the contributions of all these waves we obtain the global solution for a two-dimensional periodic medium with a large number of substructures and a low computing time. Examples are given for the case of a two-dimensional membrane with many substructures and different types of heterogeneities.</span></p></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":null,"pages":null},"PeriodicalIF":3.1,"publicationDate":"2023-11-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138455753","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":"Fatigue response of open hole plates: A finite element simulation investigating the influence of dynamic and static cold expansion processes","authors":"Guo Zheng ,&nbsp;Zengqiang Cao ,&nbsp;Yuehaoxuan Wang ,&nbsp;Reza Talemi","doi":"10.1016/j.finel.2023.104085","DOIUrl":"10.1016/j.finel.2023.104085","url":null,"abstract":"<div><p><span><span><span>Cold Expansion (CE) techniques are extensively used in the aeronautical industry to enhance the fatigue life of open-hole plates. However, the availability of accurate </span>Finite Element<span> (FE) models to simulate the fatigue behavior of this process, particularly Dynamic Cold Expansion (DCE), is limited. This study introduces two novel methods for predicting the fatigue response of DCE and Static Cold Expansion (SCE) open-hole plates. The first method directly estimates the total fatigue life using </span></span>Continuum Damage Mechanics (CDM) and the Theory of Critical Distance (TCD). The second method separates the prediction of fatigue </span>crack initiation life<span><span> and propagation life, incorporating CDM, TCD, and the Extended Finite Element Method (XFEM). Moreover, FE models are developed to simulate residual stress, stress under external </span>cyclic loads<span>, and fatigue crack propagation<span> behavior for both DCE and SCE specimens. The proposed methods are evaluated, compared, and the mechanisms behind fatigue life enhancement and fatigue crack propagation<span> modes in CE specimens are discussed. It is found that the prediction accuracy is enhanced by considering stress distributions along the thickness direction and improving the Line Method (LM) in TCD through the introduction of a novel CE parameter. The results demonstrate that both methods achieve good predictive performance<span>, with an average error index within ±30%. Furthermore, it is observed that both DCE and SCE processes primarily improve the fatigue crack initiation life of open-hole plates, with the percentage of crack initiation fatigue<span> life increasing as the expansion size increases. The majority of the fatigue crack propagation life in CE specimens is concentrated in the initial stages of crack propagation. In addition, the effects of DCE and SCE processes on reducing the fatigue crack propagation rate are more pronounced along the thickness direction compared to the width direction, leading to distinct crack propagation modes between CE and non-cold expansion (NCE) specimens.</span></span></span></span></span></span></p></div>","PeriodicalId":56133,"journal":{"name":"Finite Elements in Analysis and Design","volume":null,"pages":null},"PeriodicalIF":3.1,"publicationDate":"2023-11-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138293996","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}