Computational Mechanics最新文献_第8页

Raster approach to modelling the failure of arbitrarily inclined interfaces with structured meshes 用结构网格对任意倾斜界面的失效进行建模的栅格方法

IF 4.1 2区工程技术

Computational Mechanics Pub Date : 2024-03-11 DOI: 10.1007/s00466-024-02456-6

{"title":"Raster approach to modelling the failure of arbitrarily inclined interfaces with structured meshes","authors":"","doi":"10.1007/s00466-024-02456-6","DOIUrl":"https://doi.org/10.1007/s00466-024-02456-6","url":null,"abstract":"<h3>Abstract</h3> This paper presents an approach to evaluate the failure of arbitrarily inclined interfaces using FE models with structured spatial discretization, providing accurate prediction of crack propagation along paths known a priori that are not constrained to the element boundaries. The combination of algorithms for the generation of structured discretization of representative polycrystalline microstructures with novel cohesive element formulations allow modelling the failure of complex topologies along rasterised boundaries, with noticeably higher computational efficiency and comparable accuracy. Two formulations of raster cohesive elements are presented, adopting either elastic-brittle or Tvergaard–Hutchinson traction separation laws. The formulations proposed are first validated comparing the failure of the interface within bi-crystal structures discretised using hexahedral elements either within a structured mesh (i.e. with rasterised boundaries) or an unstructured mesh (i.e. with planar boundary). Subsequently, the effectiveness of the formulations is demonstrated comparing the inter-granular crack propagation within complex polycrystalline microstructures. The behaviour of the novel cohesive element formulation in structured meshes consisting of regular hexahedral elements is in excellent agreement with the deformation and failure of classic cohesive element formulations placed along the planar boundaries of unstructured meshes consisting of tetrahedral elements. The higher computational cost of the raster cohesive elements is more than compensated by the increase in computational efficiency of structured meshes when compared to unstructured meshes, leading to a reduction of the simulation time of up to over 200 times for the simulations presented in the paper, thus allowing the simulation of large domains.","PeriodicalId":55248,"journal":{"name":"Computational Mechanics","volume":"34 1","pages":""},"PeriodicalIF":4.1,"publicationDate":"2024-03-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140098681","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

Velocity-based space-time FEMs for solid dynamics problem: generalized framework for linear basis functions in time 用于固体动力学问题的基于速度的时空有限元：时间线性基函数的广义框架

IF 4.1 2区工程技术

Computational Mechanics Pub Date : 2024-03-11 DOI: 10.1007/s00466-024-02461-9

Vikas Sharma, Kazunori Fujisawa, Yuki Kuroda

{"title":"Velocity-based space-time FEMs for solid dynamics problem: generalized framework for linear basis functions in time","authors":"Vikas Sharma, Kazunori Fujisawa, Yuki Kuroda","doi":"10.1007/s00466-024-02461-9","DOIUrl":"https://doi.org/10.1007/s00466-024-02461-9","url":null,"abstract":"Time discontinuous Galerkin space-time finite element method (ST/FEM) can be used for developing arbitrary high-order accurate and unconditionally stable time integration schemes for elastodynamics problems. The existing ST/FEMs can be classified as the single-field and two-field ST/FEM: in the former method, either displacement or velocity, is independent and discontinuous in time. In contrast, in the latter method, both displacement and velocity fields are independent and discontinuous in time. Both methods have third-order accuracy for linear interpolation in time, higher than typical time integration schemes used in semi-discretized. However, these methods currently lack a unified computational framework, so each method requires a separate implementation. Therefore, the main goal of the present study is to develop a generalized computational framework that can facilitate the derivation and implementation of the existing linear-in-time ST/FEMs in a unified manner. This framework is developed by realizing that existing methods differ through the treatments of displacement-velocity relationships, which can be unified through displacement functions. In addition, by employing this framework, a new ST/FEM, which is designated as LC v-ST/FEM, is derived from the linear combination of displacement functions of single-field and two-field ST/FEMs. LC v-ST/FEM contains a user-defined parameter (alpha in [0,1]), which can be used for controlling the high-frequency dissipation characteristics. From finite difference analysis and numerical solutions of benchmark problems, it is demonstrated that the proposed method is the third order accurate in time, unconditionally stable, and contains negligible numerical dispersion error for all (0 le alpha le 1). Moreover, for (alpha ne 0), the method can attenuate the spurious high-frequency components from the velocity and displacement fields.","PeriodicalId":55248,"journal":{"name":"Computational Mechanics","volume":"30 1","pages":""},"PeriodicalIF":4.1,"publicationDate":"2024-03-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140098716","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

Plane crack problems within strain gradient elasticity and mixed finite element implementation 应变梯度弹性中的平面裂缝问题及混合有限元实现

IF 4.1 2区工程技术

Computational Mechanics Pub Date : 2024-03-08 DOI: 10.1007/s00466-024-02451-x

Aleksandr Yu Chirkov, Lidiia Nazarenko, Holm Altenbach

{"title":"Plane crack problems within strain gradient elasticity and mixed finite element implementation","authors":"Aleksandr Yu Chirkov, Lidiia Nazarenko, Holm Altenbach","doi":"10.1007/s00466-024-02451-x","DOIUrl":"https://doi.org/10.1007/s00466-024-02451-x","url":null,"abstract":"An alternative approach is proposed and applied to solve boundary value problems within the strain gradient elasticity theory. A mixed variation formulation of the finite element method (FEM) based on the concept of the Galerkin method is used. To construct finite-dimensional subspaces separate approximations of displacements, deformations, stresses, and their gradients are implemented by choosing the different sets of piecewise polynomial basis functions, interrelated by the stability condition of the mixed FEM approximation. This significantly simplifies the pre-requirement for approximating functions to belong to class C1 and allows one to use the simplest triangular finite elements with a linear approximation of displacements under uniform or near-uniform triangulation conditions. Global unknowns in a discrete problem are nodal displacements, while the strains and stresses and their gradients are treated as local unknowns. The conditions of existence, uniqueness, and continuous dependence of the solution on the problem’s initial data are formulated for discrete equations of mixed FEM. These are solved by a modified iteration procedure, where the global stiffness matrix for classical elasticity problems is treated as a preconditioning matrix with fictitious elastic moduli. This avoids the need to form a global stiffness matrix for the problem of strain gradient elasticity since it is enough to calculate only the residual vector in the current approximation. A set of modeling plane crack problems is solved. The obtained solutions agree with the results available in the relevant literature. Good convergence is achieved by refining the mesh for all scale parameters. All three problems under study exhibit specific qualitative features characterizing strain gradient solutions namely crack stiffness increase with length scale parameter and cusp-like closure effect.","PeriodicalId":55248,"journal":{"name":"Computational Mechanics","volume":"44 1","pages":""},"PeriodicalIF":4.1,"publicationDate":"2024-03-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140075579","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

Phase field modelling and simulation of damage occurring in human vertebra after screws fixation procedure 螺钉固定术后人体脊椎损伤的相场建模与模拟

IF 4.1 2区工程技术

Computational Mechanics Pub Date : 2024-03-07 DOI: 10.1007/s00466-024-02450-y

Deison Preve, Pietro Lenarda, Daniele Bianchi, Alessio Gizzi

{"title":"Phase field modelling and simulation of damage occurring in human vertebra after screws fixation procedure","authors":"Deison Preve, Pietro Lenarda, Daniele Bianchi, Alessio Gizzi","doi":"10.1007/s00466-024-02450-y","DOIUrl":"https://doi.org/10.1007/s00466-024-02450-y","url":null,"abstract":"The present endeavour numerically exploits the use of a phase-field model to simulate and investigate fracture patterns, deformation mechanisms, damage, and mechanical responses in a human vertebra after the incision of pedicle screws under compressive regimes. Moreover, the proposed phase field framework can elucidate scenarios where different damage patterns, such as crack nucleation sites and crack trajectories, play a role after the spine fusion procedure, considering several simulated physiological movements of the vertebral body. Spatially heterogeneous elastic properties and phase field parameters have been computationally derived from bone density estimation. A convergence analysis has been conducted for the vertebra-screws model, considering several mesh refinements, which has demonstrated good agreement with the existing literature on this topic. Consequently, by assuming different angles for the insertion of the pedicle screws and taking into account a few vertebral motion loading regimes, a plethora of numerical results characterizing the damage occurring within the vertebral model has been derived. Overall, the phase field results confirm and enrich the current literature, shed light on the medical community, which will be useful in enhancing clinical interventions and reducing post-surgery bone failure and screw loosening. The proposed computational approach also investigates the effects in terms of fracture and mechanical behaviour of the vertebral-screws body within different metastatic lesions opening towards major life threatening scenarios.","PeriodicalId":55248,"journal":{"name":"Computational Mechanics","volume":"137 1","pages":""},"PeriodicalIF":4.1,"publicationDate":"2024-03-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140053763","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

A hierarchic isogeometric hyperelastic solid-shell 分层等几何超弹性固壳

IF 4.1 2区工程技术

Computational Mechanics Pub Date : 2024-03-05 DOI: 10.1007/s00466-024-02452-w

Leonardo Leonetti, Hugo M. Verhelst

引用次数: 0

Image-based inverse characterization of in-situ microscopic composite properties 基于图像的原位微观复合材料性能反向表征

IF 4.1 2区工程技术

Computational Mechanics Pub Date : 2024-03-04 DOI: 10.1007/s00466-024-02454-8

Zimu Su, Nelson Carvalho, Michael W. Czabaj, Caglar Oskay

{"title":"Image-based inverse characterization of in-situ microscopic composite properties","authors":"Zimu Su, Nelson Carvalho, Michael W. Czabaj, Caglar Oskay","doi":"10.1007/s00466-024-02454-8","DOIUrl":"https://doi.org/10.1007/s00466-024-02454-8","url":null,"abstract":"An inverse characterization approach to identify the in-situ elastic properties of composite constituent materials is developed. The approach relies on displacement measurements available from image-based measurement techniques such as digital image correlation and template matching. An optimization problem is formulated, where the parameters of an assumed functional form describing spatially variable material properties are obtained by minimizing the discrepancies between noisy displacement measurements and the corresponding simulated values. The proposed formulation is analyzed from a statistical inference theory standpoint. It is shown that the approach exhibits estimation consistency, i.e. given noisy input data the identified material properties converge to the true material properties as the number of available measurements increases. The performance of the proposed approach is evaluated by a series of virtual characterizations that mimic physical characterization tests in which fiber centroid displacements are obtained through fiber template matching. The virtual characterizations demonstrate that the effect of measurement noise in identifying the in-situ constituent properties can be mitigated by selecting a sufficiently large measurement dataset. The numerical studies also show that, given a rich measurement dataset, the proposed approach is able to describe increasingly complex spatial variation of properties.","PeriodicalId":55248,"journal":{"name":"Computational Mechanics","volume":"71 1","pages":""},"PeriodicalIF":4.1,"publicationDate":"2024-03-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140033013","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

Computationally efficient stress reconstruction from full-field strain measurements 通过全场应变测量进行高效计算的应力重构

IF 4.1 2区工程技术

Computational Mechanics Pub Date : 2024-03-02 DOI: 10.1007/s00466-024-02458-4

Miroslav Halilovič, Bojan Starman, Sam Coppieters

{"title":"Computationally efficient stress reconstruction from full-field strain measurements","authors":"Miroslav Halilovič, Bojan Starman, Sam Coppieters","doi":"10.1007/s00466-024-02458-4","DOIUrl":"https://doi.org/10.1007/s00466-024-02458-4","url":null,"abstract":"Stress reconstruction based on experimentally acquired full-field strain measurements is computationally expensive when using conventional implicit stress integration algorithms. The computational burden associated with repetitive stress reconstruction is particularly relevant when inversely characterizing plastic material behaviour via inverse methods, like the nonlinear Virtual Fields Method (VFM). Spatial and temporal down-sampling of the available full-field strain data is often used to mitigate the computational effort. However, for metals subjected to non-linear strain paths, temporal down-sampling of the strain fields leads to erroneous stress states biasing the identification accuracy of the inverse method. Hence, a significant speedup factor of the stress integration algorithm is required to fully exploit the experimental data acquired by Digital Image Correlation (DIC). To this end, we propose an explicit stress integration algorithm that is independent on the number of images (i.e. strain fields) taken into account in the stress reconstruction. Theoretically, the proposed method eliminates the need for spatial and temporal down-sampling of the experimental full-field data used in the nonlinear VFM. Finally, the proposed algorithm is also beneficial in the emerging field of real-time DIC applications.","PeriodicalId":55248,"journal":{"name":"Computational Mechanics","volume":"79 1","pages":""},"PeriodicalIF":4.1,"publicationDate":"2024-03-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140017408","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

Fluid-evolving landform interaction by a surface-tracking method 利用地表跟踪法研究流体与演变地貌之间的相互作用

IF 4.1 2区工程技术

Computational Mechanics Pub Date : 2024-03-02 DOI: 10.1007/s00466-024-02464-6

Jorge Molina, Pablo Ortiz, Rafael Bravo

{"title":"Fluid-evolving landform interaction by a surface-tracking method","authors":"Jorge Molina, Pablo Ortiz, Rafael Bravo","doi":"10.1007/s00466-024-02464-6","DOIUrl":"https://doi.org/10.1007/s00466-024-02464-6","url":null,"abstract":"This paper introduces a continuous finite element model to simulate fluid flow-bedform interaction problems. The approach utilizes a non-oscillatory finite element algorithm to compute the fluid dynamics by solving the complete Navier–Stokes equations. Additionally, it addresses the evolution of the fluid–bedform interface as a consequence of spatially non-balanced sediment fluxes through the solution of a conservation equation for the erodible layer thickness. A sign preservation algorithm is particularly relevant for landform tracking because a positive definite thickness of the erodible sediment layer is essential to model the interaction between evolving cohesionless sediment layers and rigid beds. The fluid/terrain interface is explicitly captured through a surface tracking methodology. First, new nodes fitting the interface are incorporated into the finite element mesh; then, elements beneath this interface are deactivated, while intersected elements are restructured to get a mesh composed exclusively of tetrahedral elements. Numerical experiments demonstrate capabilities of the method by exploring relevant problems related with civil engineering, such as the evolution of trenches and the scour of a submerged pile.","PeriodicalId":55248,"journal":{"name":"Computational Mechanics","volume":"47 1","pages":""},"PeriodicalIF":4.1,"publicationDate":"2024-03-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140017544","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

A new crack-tip element for the logarithmic stress-singularity of Mode-III cracks in spring interfaces 用于弹簧界面中模式 III 裂纹的对数应力-奇异性的新裂纹尖端元素

IF 4.1 2区工程技术

Computational Mechanics Pub Date : 2024-02-24 DOI: 10.1007/s00466-024-02448-6

V. Mantič, A. Vázquez-Sánchez, M. Romero-Laborda, M. Muñoz-Reja, S. Jiménez-Alfaro, L. Távara

{"title":"A new crack-tip element for the logarithmic stress-singularity of Mode-III cracks in spring interfaces","authors":"V. Mantič, A. Vázquez-Sánchez, M. Romero-Laborda, M. Muñoz-Reja, S. Jiménez-Alfaro, L. Távara","doi":"10.1007/s00466-024-02448-6","DOIUrl":"https://doi.org/10.1007/s00466-024-02448-6","url":null,"abstract":"A new crack-tip finite element able to improve the accuracy of Finite Element Method (FEM) solutions for cracks growing along the Winkler-type spring interfaces between linear elastic adherents is proposed. The spring model for interface fracture, sometimes called Linear-Elastic (perfectly) Brittle Interface Model (LEBIM), can be used, e.g., to analyse fracture of adhesive joints with a thin adhesive layer. Recently an analytical expression for the asymptotic elastic solution with logarithmic stress-singularity at the interface crack tip considering spring-like interface behaviour under fracture Mode III was deduced by some of the authors. Based on this asymptotic solution, a special 5-node triangular crack-tip finite element is developed. The generated special singular shape functions reproduce the radial behaviour of the first main term and shadow terms of the asymptotic solution. This special element implemented in a FEM code written in Matlab has successfully passed various patch tests with spring boundary conditions. The new element allows to model cracks in spring interfaces without the need of using excessively refined FEM meshes, which is one of the current disadvantages in the use of LEBIM when stiff spring interfaces are considered. Numerical tests carried out by h-refinement of uniform meshes show that the new singular element consistently provides significantly more accurate results than the standard finite elements, especially for stiff interfaces, which could be relevant for practical applications minimizing computational costs. The new element can also be used to solve other problems with logarithmic stress-singularities.\u0000","PeriodicalId":55248,"journal":{"name":"Computational Mechanics","volume":"42 1","pages":""},"PeriodicalIF":4.1,"publicationDate":"2024-02-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139947362","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

A maximum-entropy length-orientation closure for short-fiber reinforced composites 短纤维增强复合材料的最大熵长度-取向封闭

IF 4.1 2区工程技术

Computational Mechanics Pub Date : 2024-02-24 DOI: 10.1007/s00466-024-02447-7

Alok Mehta, Matti Schneider

引用次数: 0