Mechanics of Materials最新文献

{"title":"Thermal stress analysis in a graded thermoelectric film bonded to a homogeneous substrate","authors":"Ali Farhadian,&nbsp;Yadolah Alinia","doi":"10.1016/j.mechmat.2025.105256","DOIUrl":"10.1016/j.mechmat.2025.105256","url":null,"abstract":"<div><div>In this paper, the thermoelastic behavior of a functionally graded thermoelectric thin film attached to a homogeneous substrate is investigated. First assuming a one-dimensional temperature field, the steady-state temperature distribution is obtained using the equations governing the physics of thermoelectric material. Then, the integral equations are derived for the problem by combining the equilibrium and strain compatibility equations. The simple shearing deformation is adopted for the bonding layer since its tensile strength is much smaller than that of the film and the substrate. Finally, employing the Gauss-Chebyshev discretization method, the integral equation is solved for the interfacial stress distribution. A detailed parametric study is conducted to explore the effect of non-homogeneity parameters on the stress components for the film and the substrate. The numerical results indicate that proper adjustment of the non-homogeneity parameters for the functionally graded thermoelectric thin film can extend the service life and the mechanical stability of the device. Additionally, choosing a softer and/or thicker bonding layer can reduce the stress concentration near the film ends as well as the delamination edges. The interfacial shear stress distribution changes direction when the delamination length exceeds 40% of the film length. Depositing a thinner thermoelectric film on an elastic substrate can delay delamination failure driven by interfacial shear stress.</div></div>","PeriodicalId":18296,"journal":{"name":"Mechanics of Materials","volume":"203 ","pages":"Article 105256"},"PeriodicalIF":3.4,"publicationDate":"2025-01-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143161122","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":"Exploring mechanical damage in fascia: Experiments and advanced constitutive modeling approaches","authors":"Alejandro Aparici-Gil ,&nbsp;Marta M. Pérez ,&nbsp;Estefanía Peña","doi":"10.1016/j.mechmat.2025.105239","DOIUrl":"10.1016/j.mechmat.2025.105239","url":null,"abstract":"<div><div>Biological tissues exhibit complex structures that necessitate mechanical models incorporating details of their key components and the physical processes occurring within the material. Our objective is to enhance the understanding of damage mechanisms in fibered tissues through mechanical testing. This includes conducting uniaxial tensile tests on fascia beyond physiological stretch limits and developing two constitutive models to describe damage and rupture. These models integrate both phenomenological and microstructural perspectives.</div><div>Two perpendicular directions, corresponding to the two families of collagen fibers, were compared: the longitudinal direction, characterized by greater stiffness, and the transverse direction. The mean Cauchy rupture stress (<span><math><msub><mrow><mi>σ</mi></mrow><mrow><mi>r</mi></mrow></msub></math></span>) was reported as 16.67 for the longitudinal direction and 4.76 MPa for the transverse direction, with a significant difference observed between them (<span><math><mi>p</mi></math></span>-value <span><math><mo>&lt;</mo></math></span> 0.05). Similarly, a significant difference in stored strain energy was found between the two directions (<span><math><mi>p</mi></math></span>-value <span><math><mo>&lt;</mo></math></span> 0.05) between directions, being in longitudinal equal to 1.33 <span><math><mrow><mtext>N</mtext><mi>⋅</mi><mtext>mm</mtext><mo>/</mo><msup><mrow><mtext>mm</mtext></mrow><mrow><mn>3</mn></mrow></msup></mrow></math></span> and 0.49 in transversal one. However, rupture stretches (<span><math><msub><mrow><mi>λ</mi></mrow><mrow><mi>r</mi></mrow></msub></math></span>) did not exhibit a significant difference (<span><math><mi>p</mi></math></span>-value <span><math><mo>&gt;</mo></math></span> 0.05) with values of 1.17 and 1.22 for the longitudinal and transverse directions, respectively.</div><div>In this study, a hyperelastic constitutive model for fascia was modified to incorporate damage effects into the strain energy function. Additionally, an extended version of a microstructural damage model was developed to effectively replicate the experimental data. The proposed damage models successfully captured the stress–strain behavior and accurately represented the damage process. The coefficient of determination <span><math><msup><mrow><mi>R</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span> for the fitted data ranged from 0.616 to 0.973, except for Sample IV, which exhibited an <span><math><msup><mrow><mi>R</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span> value of 0.251 when using the phenomenological model. In all cases, the microstructural model provided a more accurate fit compared to the phenomenological model, with <span><math><msup><mrow><mi>R</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span> values ranging from 0.748 to 0.927.</div></div>","PeriodicalId":18296,"journal":{"name":"Mechanics of Materials","volume":"202 ","pages":"Article 105239"},"PeriodicalIF":3.4,"publicationDate":"2025-01-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143180273","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":"Anisotropy of the fracture toughness in β-HMX crystals: A computational study","authors":"Camilo A. Duarte ,&nbsp;Catalin R. Picu ,&nbsp;Vikas Tomar ,&nbsp;WaiChing Sun","doi":"10.1016/j.mechmat.2025.105242","DOIUrl":"10.1016/j.mechmat.2025.105242","url":null,"abstract":"<div><div>The reliability, stability, and secure handling of <span><math><mi>β</mi></math></span>-cyclo-tetramethylene tetranitramine (<span><math><mi>β</mi></math></span>-HMX) crystals, a high explosive (HE) material commonly used in polymer-bonded explosives (PBX), depend heavily on their fracture properties. Cracks in HE crystals are known to localize temperature or form hotspots due to interfacial friction and can also facilitate the propagation of chemical reactions, leading to early ignition and initiation. Hence, to develop safe and reliable HEs, it is essential to characterize the fracture toughness of <span><math><mi>β</mi></math></span>-HMX, which is believed to be highly anisotropic. Furthermore, it is important to understand the origin of fracture anisotropy in HMX crystals, which is hypothesized to depend not only on the surface energy of the fracture plane, as it occurs in brittle fracture, but also on plastic deformation due to crystallographic slip. For this purpose, we performed finite element simulations of single <span><math><mi>β</mi></math></span>-HMX crystals under Mode I deformation with an atomistic-informed crystal plasticity model. Fracture toughness is estimated computationally for crystals with cracks oriented in different directions using the J-integral method. Our results confirm that the fracture toughness of HMX is highly dependent on the crystal orientation, owing to both elastic and plastic anisotropy. Furthermore, we conclude that although brittle HMX crystals may not sustain extensive plastic deformation, the contribution of plasticity to the fracture toughness is not negligible, and the anisotropy of the plastic deformation should not be neglected.</div></div>","PeriodicalId":18296,"journal":{"name":"Mechanics of Materials","volume":"202 ","pages":"Article 105242"},"PeriodicalIF":3.4,"publicationDate":"2025-01-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143180275","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":"Micromechanics-based numerical analysis of failure in calcified abdominal aortic aneurysm","authors":"Jaynandan Kumar,&nbsp;Anshul Faye","doi":"10.1016/j.mechmat.2025.105241","DOIUrl":"10.1016/j.mechmat.2025.105241","url":null,"abstract":"<div><div>Abdominal aortic aneurysms (AAAs) are a critical medical concern characterized by the dilation of the abdominal aorta, with the potential for life-threatening rupture. Calcification of AAAs in varying amount is identified as one of the factors affecting their mechanical and failure behaviour. However, reasons behind the same are not clear. The current work presents a micro-mechanics based numerical method to analyse the effect of calcification on the rupture behaviour of aneurysmatic tissue. An anisotropic material model suitable for modelling biological tissues is used and it is calibrated against available experimental data under bi-axial loading. To model failure in tissues, an energy-based failure criterion is used and failure parameters are identified from available experimental data. Calcified tissues are modelled as a composite material with tissue as a matrix and calcium (Ca) particles as an inclusion. Multiple representative volume elements are generated and used for simulation to capture the effect of morphology and amount of calcification. Contact conditions between the tissue and Ca particles are also assumed for the investigation. Thus, failure envelopes of calcified tissues are generated under different conditions. Our findings reveal that calcification affects the aneurysm rupture significantly. Amount of calcification is more critical than its morphology. Highly calcified tissues fail at lower stretches and the failure initiates in the ligaments joining Ca particles. Failure location could be correlated with available experimental observations. With higher calcification, tissues also become more isotropic in nature. The study also emphasizes that stretch-based criterion is a better candidate for predicting the failure of the aneurysm than a stress-based criterion. Further, parameters for constitutive model and failure model are identified for homogenized calcified tissue with low calcification.</div></div>","PeriodicalId":18296,"journal":{"name":"Mechanics of Materials","volume":"202 ","pages":"Article 105241"},"PeriodicalIF":3.4,"publicationDate":"2025-01-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143179398","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":"Enhancing robustness in machine-learning-accelerated molecular dynamics: A multi-model nonparametric probabilistic approach","authors":"Ariana Quek ,&nbsp;Niuchang Ouyang ,&nbsp;Hung-Min Lin ,&nbsp;Olivier Delaire ,&nbsp;Johann Guilleminot","doi":"10.1016/j.mechmat.2024.105237","DOIUrl":"10.1016/j.mechmat.2024.105237","url":null,"abstract":"<div><div>In this work, we present a system-agnostic probabilistic framework to quantify model-form uncertainties in molecular dynamics (MD) simulations based on machine-learned (ML) interatomic potentials. Such uncertainties arise from the design and selection of ML potentials, as well as from training aspects pertaining to the definition of datasets and calibration strategies. Our approach relies on a stochastic reduced-order model (SROM) where the approximation space is expanded through the randomization of the projection basis. The construction of the underlying probability measure is achieved in the context of information theory, by leveraging the existence of multiple model candidates, corresponding to different ML potentials for instance. To assess the effectiveness of the proposed approach, the method is applied to capture model-form uncertainties in a sodium thiophosphate system, relevant to sodium-ion-state batteries. We demonstrate that the SROM accurately encodes model uncertainties from different ML potentials – including a Neuro-Evolution Potential (NEP) and a Moment Tensor Potential (MTP) – and can be used to propagate these uncertainties to macroscopic quantities of interest, such as ionic diffusivity. Additionally, we investigate the impact of augmenting the snapshot matrix with momenta, and of introducing a frequency-based split in the construction of the random projection matrix. Results indicate that including momenta improves the accuracy of the SROM, while frequency splitting enables stabilization around nominal responses during uncertainty propagation. The proposed enhancements contribute to more robust and stable predictions in MD simulations involving ML potentials.</div></div>","PeriodicalId":18296,"journal":{"name":"Mechanics of Materials","volume":"202 ","pages":"Article 105237"},"PeriodicalIF":3.4,"publicationDate":"2025-01-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143179942","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":"The bending-buckling coupled model for blistering behavior in anti-corrosion coatings","authors":"Liangji Ma ,&nbsp;Yin Yao ,&nbsp;Bo Zhang ,&nbsp;Zhilong Peng ,&nbsp;Shaohua Chen","doi":"10.1016/j.mechmat.2024.105238","DOIUrl":"10.1016/j.mechmat.2024.105238","url":null,"abstract":"<div><div>Anti-corrosion coatings are widely applied in marine engineering and marine equipment. Understanding their blistering failure mechanisms is vital for optimizing coating designs and extending their service life. This paper develops a bending-buckling coupled model and employs the Rayleigh-Ritz method to investigate the axisymmetric circular blister of coatings with initial deflections, with a special focus on the situation where the transverse load is opposite to the deflection. By incorporating the chemo-mechanical coupling, an analytical solution of the critical buckling load in terms of diffusion strain is well achieved, concisely explaining the impacts of transverse loads, initial deflections, and aspect ratios on the critical buckling load. The influence of these parameters on the post-buckling behavior of the coating is further discussed and the contour of the coating blister can be presented. The results should have theoretical guidance significance for predicting and analyzing the service behavior of anti-corrosion coatings.</div></div>","PeriodicalId":18296,"journal":{"name":"Mechanics of Materials","volume":"203 ","pages":"Article 105238"},"PeriodicalIF":3.4,"publicationDate":"2025-01-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143160604","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 additive Mori–Tanaka scheme for elastic–viscoplastic composites based on a modified tangent linearization","authors":"K. Kowalczyk-Gajewska ,&nbsp;S. Berbenni ,&nbsp;S. Mercier","doi":"10.1016/j.mechmat.2024.105191","DOIUrl":"10.1016/j.mechmat.2024.105191","url":null,"abstract":"<div><div>Mean-field modeling based on the Eshelby inclusion problem poses some difficulties when the non-linear Maxwell-type constitutive law is used for elasto–viscoplasticity. One difficulty is that this behavior involves different orders of time differentiation, which leads a long-term memory effect. One of the possible solutions to this problem is the additive interaction law. Generally, mean field models solely use the mean values of stress and strain fields per phase, while variational approaches consider the second moments of stresses and strains. It is seen that the latter approach improves model predictions allowing to account for stress fluctuation within the phases. However, the complexity of the variational formulations still makes them difficult to apply in the large scale finite element calculations and for non-proportional loadings. Thus, there is a need to include the second moments within homogenization models based on the additive interaction law. In the present study, the incorporation of the second moments of stresses into the formulation of the additive Mori–Tanaka model of two-phase elastic–viscoplastic material is discussed. A modified tangent linearization of the viscoplastic law is proposed, while the Hill–Mandel’s lemma is used to track the evolution of second moments of stresses. To study the model performance and efficiency, the results are compared to the full-field numerical calculations and predictions of other models available in the literature. Very good performance of the modified tangent linearization is demonstrated from these benchmarks for both monotonic and non monotonic loading responses.</div></div>","PeriodicalId":18296,"journal":{"name":"Mechanics of Materials","volume":"200 ","pages":"Article 105191"},"PeriodicalIF":3.4,"publicationDate":"2025-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143132827","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":"Identification for elastoplastic constitutive parameters of 316L stainless steel lattice structures using finite element model updating and integrated digital image correlation","authors":"Zhaozhen Huang,&nbsp;Caroline Antion,&nbsp;Franck Toussaint","doi":"10.1016/j.mechmat.2024.105232","DOIUrl":"10.1016/j.mechmat.2024.105232","url":null,"abstract":"<div><div>Lattice structures are widely considered for industrial applications owing to their excellent energy absorption and mechanical properties. In this work, octet-truss lattice structures are manufactured from 316L stainless steel powder by selective laser melting (SLM). The geometrical information of lattice structures is captured by SEM and X-ray tomography. It reveals that realistic dimensions of struts differ slightly from CAD-designed ones. The mechanical behaviors are investigated both experimentally and numerically. Quasi-static uni-axial compression experiments with 2D digital image correlation (DIC) technology are conducted to measure displacement/strain fields. Finite element analysis based on an elastic and anisotropic plastic constitutive model is used to simulate mechanical behaviors. To improve the predictive accuracy, a finite element model updating approach is implemented to identify constitutive parameters. The results show that numerical simulation with optimized parameters match well with experiments in aspect of force-displacement curve at elastic–plastic stage and displacement fields.</div></div>","PeriodicalId":18296,"journal":{"name":"Mechanics of Materials","volume":"202 ","pages":"Article 105232"},"PeriodicalIF":3.4,"publicationDate":"2024-12-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143179401","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":"Relating stiffness changes in porous materials to the evolution of pore space","authors":"Yulia Pronina ,&nbsp;Maria Narykova ,&nbsp;Mark Kachanov","doi":"10.1016/j.mechmat.2024.105236","DOIUrl":"10.1016/j.mechmat.2024.105236","url":null,"abstract":"<div><div>The work aims at relating stiffness changes in porous materials to the evolution of pore space geometry. After a brief review of the relevant micromechanics tools, we apply them to case studies on several metals. In particular, it is clarified, when porosity can or cannot be used as a single quantitative characteristic of the pore space in whose terms the effective stiffness is to be expressed, and when it must be changed to crack density. Namely, the use of porosity parameter is legitimate in cases of isotropic mixtures of pores having approximately equal shape factors, provided the shapes are not strongly oblate (aspect ratios larger than about 0.08). Considered examples show that, in cases of strongly oblate, crack-like pores, noticeable stiffness changes may occur at very low values of porosity; in such cases, the crack density parameter must be used. Besides predicting the effective stiffness in terms of proper characteristics of the pore space, the developed methodology allows monitoring the evolution of pore shapes based on stiffness changes and porosity data. In our analysis, pore geometries are modeled by spheroids of appropriate aspect ratios; they provide sufficient flexibility and allow quantitative modeling. The adequacy of such modeling is supported by agreement of the theoretical results with experimental data.</div></div>","PeriodicalId":18296,"journal":{"name":"Mechanics of Materials","volume":"202 ","pages":"Article 105236"},"PeriodicalIF":3.4,"publicationDate":"2024-12-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143179399","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":"Homogenised modelling of the electro-mechanical behaviour of a vascularised poroelastic composite representing the myocardium","authors":"Laura Miller ,&nbsp;Raimondo Penta","doi":"10.1016/j.mechmat.2024.105215","DOIUrl":"10.1016/j.mechmat.2024.105215","url":null,"abstract":"<div><div>We propose a novel model for a vascularised poroelastic composite representing the myocardium which incorporates both mechanical deformations and electrical conductivity. Our structure comprises a vascularised poroelastic extracellular matrix with an embedded elastic inclusions (representing the myocytes) and we consider the electrical conductance between these two solid compartments. There is a distinct lengthscale separation between the scale where we can visibly see the connected fluid compartment separated from the poroelastic matrix and the elastic myocyte and the overall size of the heart muscle. We therefore apply the asymptotic homogenisation technique to derive the new model. The effective governing equations that we obtain describe the behaviour of the myocardium in terms of the zero-th order stresses, current densities, relative fluid–solid velocities, pressures, electric potentials and elastic displacements. It effectively accounts for the fluid filling in the pores of the poroelastic matrix, flow in the vessels, the transport of fluid between the vessels and the matrix, and the elastic deformation and electrical conductance between the poroelastic matrix and the myocyte. This work paves the way towards a myocardium model that incorporates multiscale deformations and electrical conductivity whilst also considering the effects of the vascularisation and indeed the impact on mechanotransduction.</div></div>","PeriodicalId":18296,"journal":{"name":"Mechanics of Materials","volume":"202 ","pages":"Article 105215"},"PeriodicalIF":3.4,"publicationDate":"2024-12-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143179396","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}