International Journal of Solids and Structures最新文献

{"title":"Wrinkling analysis of a stiff shallow film mounted on a cylindrically curved compliant substrate, Part I: Cylindrical mode","authors":"Alexis Kordolemis ,&nbsp;Antonios E. Giannakopoulos","doi":"10.1016/j.ijsolstr.2025.113382","DOIUrl":"10.1016/j.ijsolstr.2025.113382","url":null,"abstract":"<div><div>It is well known that when a material bilayer which consists of a stiff thin film perfectly bonded on a compliant substrate develops mechanical instabilities upon the application of in-plane compressive loads. This paper investigates in depth a particular surface instability, namely the one dimensional cylindrical mode, for the case where the initial configuration of the bilayer is cylindrically curved and is acted upon membrane biaxial strains. The investigation of the role of the curvature-induced anisotropy to the geometrical characteristics of the cylindrical wrinkling mode, i.e. the wavelength and the amplitude, is the main focus of this study. To this end, the film is modelled as a thin shallow shell and is analysed within the framework of the general shallow shell theory while the mechanical behaviour of the much softer substrate is investigated through linear elasticity theory for solids. Two cylindrical modes are presented analytically, namely along and across the ridge of the initial curvature, and it has been found that the latter is the favoured wrinkling mode due to the lower total energy of the bilayer in the buckled state. The perturbation method has been employed for the calculation of the displacement components in the compliant substrate and semi-analytical expressions for the wavelengths and the amplitudes are presented which may act as a useful design tool. The results of the present study are in good agreement compared to the results of other studies available in the open literature.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"317 ","pages":"Article 113382"},"PeriodicalIF":3.4,"publicationDate":"2025-04-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143867724","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":"Global stability formulation of a non-unified Miura origami-patterned slender tube","authors":"Yangqing Liu ,&nbsp;Xinrui Zhang ,&nbsp;Qin Yu ,&nbsp;Zhiqiang Liu ,&nbsp;Xinyu Cai ,&nbsp;Marco Meloni","doi":"10.1016/j.ijsolstr.2025.113406","DOIUrl":"10.1016/j.ijsolstr.2025.113406","url":null,"abstract":"<div><div>In structural engineering, axially compressed members with relatively high slenderness ratios are prone to global buckling. A promising solution to this issue is a tailored non-unified Miura-origami crease scheme, which has been proven to significantly mitigate global buckling in non-slender and slender tubes under axial compression. Building on previous work, this manuscript advances the research into stability evaluation of the patterned slender tube and proposes a quantification method for their global stability. First, the inertia moment of the patterned cross-section is theoretically derived, and the slenderness ratio of the tube is calculated with a modification coefficient of the critical load. Subsequently, a method for calculating the stability coefficient is obtained after parametric investigations on the effect of the geometry on the yield load by nonlinear finite element simulations. Finally, a quantification method for the global stability of the tube is proposed. The study shows that, in contrast to the typical Eulerian formula, the proposed method accounts for the reduction effect of the creases on the stability of the tube and effectively predicts its global stability. This method is expected to serve as a practical tool to promote and facilitate the application of origami-patterned tubes in structural engineering.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"317 ","pages":"Article 113406"},"PeriodicalIF":3.4,"publicationDate":"2025-04-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143876842","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":"Wrinkling analysis of a stiff shallow film mounted on a cylindrically curved compliant substrate, Part II: Checkerboard and hexagonal modes","authors":"Alexis Kordolemis ,&nbsp;Antonios E. Giannakopoulos","doi":"10.1016/j.ijsolstr.2025.113383","DOIUrl":"10.1016/j.ijsolstr.2025.113383","url":null,"abstract":"<div><div>A thin film mounted on a compliant substrate under biaxial compressive strains wrinkles in different patterns such as cylindrical, checkerboard, herringbone and hexagonal. In this paper, we provide a thorough analysis of the checkerboard and hexagonal wrinkling modes which have been well documented experimentally and numerically. Particular attention is paid to the role of the curvature-induced anisotropy in the geometrical characteristic of these particular surface wrinkling modes, <em>i.e</em>. the wavelengths in the two principal directions as well as the corresponding amplitudes. The film is assumed to be much stiffer than the substrate and the bilayer system is cylindrically curved and is acted upon biaxial compressive strains. The film is modelled as a shallow shell with finite rotations while the substrate is simulated as a linear three-dimensional elastic solid. Utilizing the minimization of the total energy of the system semi-analytical expressions for the critical values of the wavelengths and the corresponding amplitudes associated with the onset of the checkerboard and hexagonal mode are provided. The obtained results has been found to be in a very good agreement compared to experimental and numerical findings of other studies. It is shown that the presence of the initial curvature in the bilayer delays the critical strain and the wrinkling amplitudes significantly for both modes, compared to the flat system, and moreover explains the inward buckling of the hexagonal mode which has been observed experimentally.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"317 ","pages":"Article 113383"},"PeriodicalIF":3.4,"publicationDate":"2025-04-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143868217","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A physical construction of the Ogden-Hill equation for soft elastomeric networks","authors":"Ziyu Xing","doi":"10.1016/j.ijsolstr.2025.113393","DOIUrl":"10.1016/j.ijsolstr.2025.113393","url":null,"abstract":"<div><div>The theory of rubber elasticity stands as a cornerstone in the study of soft matter. Despite nearly a century of development, many models remain largely phenomenological in nature, lacking a firm physical grounding. Addressing this fundamental challenge, this study advances a strain energy function theory of molecular basis, drawing from the worm-like chain model and the tube-like entanglement model. Diverging from classical approaches like the freely-jointed chain (FJC) model, the proposed framework offers a nuanced analysis of semi-flexible chains based on their end-to-end distance distributions, thus providing a more comprehensive understanding of polymer mechanics. Moreover, the entanglement of chains is characterized through the assessment of tube potential energy. Through a defined set of parameters, the model adeptly predicts the large deformation behaviors of vulcanized rubber across various experimental conditions including uniaxial tension, uniaxial compression, pure shear, and equi-biaxial tension. Additionally, it offers analytical insights into phenomena such as uniaxial tension and the inflation of an ideal balloon based on a set of parameters. During its application, the model’s resemblance to the Ogden-Hill form was noted, prompting a comparative analysis with well-known equations in this form (e.g., Varga equation, Neo-Hookean equation, Mooney-Rivlin equation, Mullins-Tobin equation, classical Ogden-Hill equation), thereby elucidating the physical underpinnings of the strain energy function. The proposed model posits the strain energy function as comprising three distinct components, i.e., affine motion, semi-flexibility, and entanglement—each manifesting distinct mechanical characteristics denoted by invariants. Furthermore, comparative assessments against the Carroll model, Pucci-Saccomnadi model, full-chain FJC model and the semi-flexible worm-like chain (WLC) model underscore the advantages of the proposed framework. Notably, the model exhibits a capacity to accurately predict rubber stress responses under large deformations without encountering singularities, thus rendering it amenable to finite element analysis (FEA). Finally, the efficacy of the proposed constitutive models is corroborated through rigorous comparisons with experimental data drawn from a spectrum of literature sources encompassing vulcanized rubber, natural rubber, elastomeric hydrogel, supramolecular elastomeric networks, and highly entangled elastomeric networks.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"316 ","pages":"Article 113393"},"PeriodicalIF":3.4,"publicationDate":"2025-04-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143859812","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A geometrically nonlinear reduced-order method using hybrid-stress solid-shell formulations for 3D large deformation analysis of thin-walled structures","authors":"Ke Liang ,&nbsp;Jiaqi Mu ,&nbsp;Saullo G.P. Castro","doi":"10.1016/j.ijsolstr.2025.113385","DOIUrl":"10.1016/j.ijsolstr.2025.113385","url":null,"abstract":"<div><div>Post-buckling of thin-walled aeronautical structures induces failure modes related to skin-stiffener separation that require three-dimensional large deformation analyses for an accurate numerical prediction. Conventional finite elements are capable of solving such analyses, but with high associated computational costs, being therefore more utilized for the simulation of smaller structural components, or even restricted to perform virtual testing at coupon-level models. In this paper, a geometrically nonlinear reduced-order method using a hybrid-stress solid-shell formulation is proposed for large deformation analysis of thin-walled structures. Current reduced-order methods are mainly applicable to buckling problems, considering only the out-of-plane deformation. Furthermore, existing displacement-based reduced models involve a computationally expensive fourth-order tensor obtained with the higher-order strain energy variations. Here, a reduced-order model with only one degree of freedom is constructed for both in-plane and out-of-plane large deformation problems. It is shown that, with the hybrid-stress formulation, the constructional efficiency of the reduced system is largely improved by zeroing the fourth-order strain energy variation using the two-field Hellinger–Reissner variational principle, followed by a condensation of the stress terms that lead to a third-order approximation of the equilibrium equation. The nonlinear predictor solved by the reduced-order model can be corrected when its numerical accuracy is not satisfactory during the path-following analysis. A simple plate, a honeycomb cell with negative Poisson ratio, and a swept-back wing structure; are used as numerical examples to verify that the proposed method enables a superior path-following capability for the three-dimensional analysis of thin-walled structures undergoing large deflection, large rotation and large strains. Furthermore, an experimental validation of the proposed method is presented using a variable-thickness plate with mixed composite-metallic materials, undergoing large out-of-plane deflection.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"317 ","pages":"Article 113385"},"PeriodicalIF":3.4,"publicationDate":"2025-04-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143874734","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":"Characterization of residual stress distribution of film/substrate structure with elastic boundary: Theoretical modeling and experimental identification","authors":"Kunjie Sun ,&nbsp;Chengyi Shou ,&nbsp;Jubing Chen ,&nbsp;Chen Sun","doi":"10.1016/j.ijsolstr.2025.113379","DOIUrl":"10.1016/j.ijsolstr.2025.113379","url":null,"abstract":"<div><div>Film/substrate structures are extensively utilized in many fields due to their unique properties. Curvature-based techniques (CBTs) are employed to measure and analyze residual stresses within these structures. While existing theories have relaxed many limitations of CBTs, they primarily focus on the free boundaries of circular objects. To address this issue, this work presents an analytical method for evaluating residual stresses in film/substrate structures subjected to arbitrary misfit strains under elastic boundary conditions for the first time, and the proposed theory can be applied to objects of arbitrary shape. Based on linear elasticity theory, a system of coupled governing equations involving seven variables is established. The equations are decoupled and solved using Fourier series expansion and the method of constant variation. Random displacement fields with uniform, normal, and exponential distributions are utilized to examine the convergence of the results under various boundary stiffness conditions. Finite Element Method (FEM) is employed to validate the film stresses under different boundary conditions. Using Stoney formula as a rough estimate, an iterative algorithm has been developed to calculate equivalent boundary conditions for circular measurement boundaries. This algorithm enables dividing the object into multiple circular regions for individual measurement, facilitating the assessment of arbitrarily shaped objects. Finally, a monoscopic fringe reflectometry measurement system is developed to experimentally characterize the topography of specular surface. This system provides precise and efficient measurement of the topography of film/substrate structure. Utilizing this data, the full-field distribution of residual stresses is subsequently determined through an inversion process.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"316 ","pages":"Article 113379"},"PeriodicalIF":3.4,"publicationDate":"2025-04-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143859813","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-free temperature fields in thermoelastic quasicrystals","authors":"Viktoriya Pasternak ,&nbsp;Iaroslav M. Pasternak ,&nbsp;Heorhiy Sulym ,&nbsp;Ihor Hotsyk ,&nbsp;Roman Pasternak","doi":"10.1016/j.ijsolstr.2025.113390","DOIUrl":"10.1016/j.ijsolstr.2025.113390","url":null,"abstract":"<div><div>This paper studies the thermoelasticity of quasicrystal solids and identifies temperature distributions that do not induce thermal stresses. Unlike anisotropic crystalline solids, where a linear temperature distribution results in a stress-free state, such a distribution can generate thermal stresses in quasicrystal media. Thermal stresses arise due to incompatible thermal strain. Compatibility conditions for phason strain are derived, and stress-free linear temperature distributions are presented for different types of quasicrystals. Special attention is given to plane strain conditions, demonstrating the two-step mechanism by which thermal stresses develop in quasicrystals under plane strain. Using the Stroh formalism and a least-squares approach, the problem is solved for a finite quasicrystal solid, revealing a significantly nonlinear distribution of thermal stresses in response to certain linear temperature distributions (uniform heat flux).</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"316 ","pages":"Article 113390"},"PeriodicalIF":3.4,"publicationDate":"2025-04-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143850696","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":"Elastic wave propagation in magneto-active fibre composites","authors":"Harold Berjamin,&nbsp;Stephan Rudykh","doi":"10.1016/j.ijsolstr.2025.113373","DOIUrl":"10.1016/j.ijsolstr.2025.113373","url":null,"abstract":"<div><div>Fibre-reinforced elastomers are lightweight and strong materials that can sustain large deformations. When filled with magnetic particles, their effective mechanical response can be modified by an external magnetic field. In the present study, we propose an effective theory of fibre-reinforced composite, based on a neo-Hookean elastic response and a linear magnetic law in each phase. The theory is shown suitable to describe the motion of composite cylinders. Furthermore, it is found appropriate for the modelling of fibre-reinforced composites subjected to a permanent magnetic field aligned with the fibres. To reach this result, we use the incremental theory (‘small on large’), in combination with homogenisation theory and the Bloch–Floquet method. This way, we show that wave directivity is sensitive to the application of a permanent magnetic field, whereas the frequency range in which wave propagation is forbidden is not modified by such a load (the band gaps are invariant). In passing, we describe a method to deduce the total stress in the material based on the measurement of two wave speeds. Furthermore, we propose an effective energy function for the description of nonlinear composites made of Yeoh-type generalised neo-Hookean fibres within a neo-Hookean matrix.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"316 ","pages":"Article 113373"},"PeriodicalIF":3.4,"publicationDate":"2025-04-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143852067","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":"Eulerian rates of elastic incompatibilities for crystal plasticity applied to size-dependent hardening in finite bending","authors":"Lorenzo Bardella ,&nbsp;M.B. Rubin ,&nbsp;Andrea Panteghini","doi":"10.1016/j.ijsolstr.2025.113376","DOIUrl":"10.1016/j.ijsolstr.2025.113376","url":null,"abstract":"<div><div>By following the work in (Rubin and Bardella, 2024), this investigation develops measures of rates of elastic incompatibilities, denoted as <span><math><msub><mrow><mi>R</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub></math></span>, for crystal plasticity. This effort relies on Eulerian constitutive equations for finite-deformation anisotropic elastoplasticity governed by the evolution of microstructural material vectors. The rates <span><math><msub><mrow><mi>R</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub></math></span> depend on the crystallography as the latter enters the rate of plasticity <span><math><msub><mrow><mi>L</mi></mrow><mrow><mi>p</mi></mrow></msub></math></span> and the <span><math><msub><mrow><mi>R</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub></math></span> are obtained by evaluating the opposite of the current curl of <span><math><msub><mrow><mi>L</mi></mrow><mrow><mi>p</mi></mrow></msub></math></span> relative to the microstructural vectors. Each component of <span><math><msub><mrow><mi>R</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub></math></span> is invariant under superposed rigid body motions, such that it can be independently employed in the constitutive equations. In crystal plasticity, the adopted Eulerian framework allows for singling out in <span><math><msub><mrow><mi>R</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub></math></span> the contributions due to rates of densities of geometrically necessary dislocations and to the elastic distortion of the crystal lattice. In this work, <span><math><msub><mrow><mi>R</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub></math></span> are used to enhance the hardening, which is applied to the size-dependent material response of structurally thick circular sectors subjected to uniform large-deformation bending.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"316 ","pages":"Article 113376"},"PeriodicalIF":3.4,"publicationDate":"2025-04-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143842664","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":"Elastic–plastic crack-tip-opening-displacement-based description for surface, corner and embedded cracks tip stress field","authors":"Jianqiang Zhang ,&nbsp;Pengfei Cui ,&nbsp;Wanlin Guo","doi":"10.1016/j.ijsolstr.2025.113395","DOIUrl":"10.1016/j.ijsolstr.2025.113395","url":null,"abstract":"<div><div>Since part-through crack growth stages occupy most of crack growth life of engineering structures, it is essential to investigate the fracture parameters of part-through cracks. However, the complex three-dimensional (3D) stress states make it difficult to efficiently dominate the crack-tip fields. Here, the 3D elastic–plastic stress intensity factor <em>K_δ-Tz</em> is extended to dominate part-through cracks. Systematic 3D finite element (FE) analyses are conducted for typical part-through cracks (embedded, corner, and surface cracks) considering different elliptical ratios and hardening exponents. It is found that the predicted stress distributions by the <em>δ-T_z</em> solution agree well with 3D FE results. Additionally, the predictive performance of the <em>δ-T_z</em> solution improves with increasing hardening exponents. Across all experimental and numerical results, the variation of <em>J</em>-integral along the crack front line can reach 200%, while remaining within 21% for <em>K_δ-Tz</em>. These results demonstrate that <em>K_δ-Tz</em> can reduce geometric constraints effectively and be a more stable elastic–plastic fracture parameter for part-through cracks in engineering structures.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"317 ","pages":"Article 113395"},"PeriodicalIF":3.4,"publicationDate":"2025-04-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143867725","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}