International Journal of Solids and Structures最新文献

{"title":"The global-local mechanical behaviors of multilayered structure and applications to superconducting coils","authors":"Sijian Wang ,&nbsp;Yuchen Han ,&nbsp;Huadong Yong ,&nbsp;Youhe Zhou","doi":"10.1016/j.ijsolstr.2025.113689","DOIUrl":"10.1016/j.ijsolstr.2025.113689","url":null,"abstract":"<div><div>Differently from the continuum medium, multilayered structures are discontinuous systems, which consist of a series of plates or shells stacked on each other. The global kinematic feature of the multilayered structure is determined by the local kinematics of its monolayers and the interfacial movement between them. At the interface, contact pressure works in the normal direction, and friction works in the tangential direction, which remarkably impacts the mechanical responses of multilayered structures. The multilayered structures have the following kinematic characteristics: interfacial slip is kinematically permissible and the contact gap can be negligible during deformation. Based on these features, the midplane displacements of each shell layer can be constructed as a continuous field. Layer-to-layer interactions are treated as internal forces, and the Coulomb friction law is incorporated as a material constitutive relationship. In kinematics, interfacial slip is taken into account in the global strain–displacement relationship, which leads to a correction term in the governing equation according to the dual relationship. In dynamics, the principle of minimum potential energy is utilized to derive the governing equation, with friction being accounted for through the mechanism of energy dissipation. Then, a continuum theory for multilayered structures is proposed in this paper and validated by comparisons with the discrete contact model and experiments. The proposed model is further extended to the mechanical study of high-temperature superconducting (HTS) magnets. The complexity of HTS magnets stems from the large number of contacts and the corresponding contact nonlinearity, which gives rise to computational inefficiencies and convergence problems. The proposed model could address these challenges and facilitate the mechanical analysis of HTS magnets.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"324 ","pages":"Article 113689"},"PeriodicalIF":3.8,"publicationDate":"2025-09-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145269770","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":"Large plastic deformation of voids in crystals","authors":"Jalal Smiri ,&nbsp;Joseph Paux ,&nbsp;Oğuz Umut Salman ,&nbsp;Ioan R. Ionescu","doi":"10.1016/j.ijsolstr.2025.113657","DOIUrl":"10.1016/j.ijsolstr.2025.113657","url":null,"abstract":"<div><div>The mechanisms of void growth and coalescence are key contributors to the ductile failure of crystalline materials. At the grain scale, single crystal plastic anisotropy induces large strain localization leading to complex shape evolutions. In this study, an Arbitrary Lagrangian–Eulerian (ALE) framework for a 2D Eulerian crystal plasticity model combined with dynamic remeshing is used to study the 2D shape evolution of cylindrical voids in single crystals. The large deformation and shape evolution of the voids under two types of loading are considered: (i) radial and (ii) uni-axial loadings. In both cases, the voids undergo complex shape evolutions induced by the interactions between slip bands, lattice rotations and large strain phenomena.</div><div>In case (i), the onset of the deformation revealed the formation of a complex fractal network of slip bands around the voids. Then, large deformations unearth an unexpected evolution of the slip bands network associated with significant lattice rotations, leading to a final hexagonal shape for the void. In case (ii), we obtain shear bands with very large accumulated plastic strain (<span><math><mrow><mo>&gt;</mo><mn>200</mn><mtext>%</mtext></mrow></math></span>) compared to the macroscopic engineering strains (<span><math><mrow><mo>&lt;</mo><mn>15</mn><mtext>%</mtext></mrow></math></span>). A high dependence between crystalline orientations, slip band localization and therefore shape evolution was observed, concluding in a high dependency between crystalline orientation and void shape elongation, which is of prime importance regarding coalescence of the voids, and thus to the formation of macro-cracks.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"324 ","pages":"Article 113657"},"PeriodicalIF":3.8,"publicationDate":"2025-09-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145268766","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":"Constitutive modeling of plastic anisotropy and strength differential effect of commercially pure titanium and its application to finite element analysis of circular cup drawing","authors":"Yuji Sato ,&nbsp;Honghao Wang ,&nbsp;Akira Yoshimura ,&nbsp;Jun Yanagimoto","doi":"10.1016/j.ijsolstr.2025.113686","DOIUrl":"10.1016/j.ijsolstr.2025.113686","url":null,"abstract":"<div><div>Owing to its strong corrosion resistance, moderate mechanical properties, and excellent biocompatibility, commercially pure titanium (CP Ti) has been widely applied in the fields of marine sciences, aerospace, and biomedicine. To reduce the trial-and-error production of CP Ti components and optimize the forming processes, an advanced constitutive model for accurate forming simulation of CP Ti was developed in this study based on the non-associated flow rule. The developed non-associated constitutive model can capture the orthotropic flow stress anisotropy, tension–compression strength differential (T-C SD) effect, and planar deformation anisotropy of CP Ti, and it was validated using material characterization data of a hot-rolled CP Ti plate. The developed constitutive model was implemented in Abaqus/Standard via the user subroutine UMAT and achieved an accurate finite element (FE) analysis of the circular cup drawing for the CP Ti sheet. All the material properties considered in the developed constitutive model affect the FE analysis results, which indicate the necessity of considering the orthotropic flow stress anisotropy, T-C SD effect, and in-plane deformation anisotropy during the sheet metal forming simulation of CP Ti.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"324 ","pages":"Article 113686"},"PeriodicalIF":3.8,"publicationDate":"2025-09-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145268765","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 constitutive and numerical model for P-formulation of flexoelectricity with type II strain gradient elasticity","authors":"Thuc Pham-Phu,&nbsp;Sergey Kozinov","doi":"10.1016/j.ijsolstr.2025.113668","DOIUrl":"10.1016/j.ijsolstr.2025.113668","url":null,"abstract":"<div><div>As the size of sensors, actuators, and electronic components is driven to the micro- and nanoscale ranges, flexoelectricity, a higher order electromechanical coupling effect, becomes more prominent. Conventional finite element framework requires elements of at least <span><math><msup><mrow><mi>C</mi></mrow><mrow><mn>1</mn></mrow></msup></math></span>-continuity to resolve strain gradient. This requirement can be circumvented by using mixed finite element method, which considers strain components as independent degrees of freedom. Most previous studies have employed the <span><math><mi>E</mi></math></span>-formulation because it is relatively straightforward to implement. However, developing a <span><math><mi>P</mi></math></span>-formulation is necessary, as it enables phase-field modeling of ferroelectric domain evolution (with electric polarization <span><math><mi>P</mi></math></span> serving as an order parameter) that is additionally influenced by strain gradients. Current research based on type II strain gradient elasticity presents a mixed finite element formulation capable of capturing the flexoelectric coupling effect using a kinematic constraint between relaxed strain and mechanical strain. The proposed mixed finite element is verified with a patch test and a thick cylinder problem, whose results are extensively discussed.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"324 ","pages":"Article 113668"},"PeriodicalIF":3.8,"publicationDate":"2025-09-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145268778","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":"Comment on “A critical note on the role of the contact stiffness in fracture of the multi-punch contact interface” by I. Argatov","authors":"Yongbin Wang ,&nbsp;Jinsheng Zhao ,&nbsp;Yuxiang He ,&nbsp;Mingshan Yang ,&nbsp;Jielei Chu ,&nbsp;Jianghong Yuan ,&nbsp;Weiqiu Chen ,&nbsp;Xiangyu Li","doi":"10.1016/j.ijsolstr.2025.113684","DOIUrl":"10.1016/j.ijsolstr.2025.113684","url":null,"abstract":"","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"324 ","pages":"Article 113684"},"PeriodicalIF":3.8,"publicationDate":"2025-09-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145218369","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":"Applications of a continuum model for brush gels","authors":"Yifei Ren,&nbsp;Prashant K. Purohit","doi":"10.1016/j.ijsolstr.2025.113664","DOIUrl":"10.1016/j.ijsolstr.2025.113664","url":null,"abstract":"<div><div>Brush gels, defined by their branched molecular architecture, exhibit unique mechanical adaptability and tunable properties. A continuum model is adopted to describe their mechanical behavior, governed by a parameter <span><math><mi>β</mi></math></span> related to the mean-squared end-to-end distance of the underlying polymer chain, which controls the transition from tensile-dilating to tensile-contracting behavior under uniaxial tension. This and other topological parameters of brush gels can be engineered to produce mechanical behavior like that of linear polymer gels and fibrous gels. Adjustments in branch topology of the polymer, such as grafting density and side-chain length, allow for tuning the stiffness while holding the shape of the force–stretch curve fixed. The rate-dependence of the force–stretch curve is more pronounced under compression due to osmotic effects. Rotational shear tests reveal the impact of chemical potential gradients driven by normal stress variations and it is shown that pre-compression can enhance network stiffness. Furthermore, brush gels exhibit increased resistance to drying-induced volume changes, attributed to their network architecture, which enables rapid equilibration. These findings underscore the potential of brush gels in applications requiring mechanical adaptability and stability under drying conditions.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"324 ","pages":"Article 113664"},"PeriodicalIF":3.8,"publicationDate":"2025-09-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145156739","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 two-step strategy for designing stretch-dominated metamaterials simultaneously satisfying desired elastic tensor and bandgap","authors":"Zijian Wang,&nbsp;Hua Deng","doi":"10.1016/j.ijsolstr.2025.113683","DOIUrl":"10.1016/j.ijsolstr.2025.113683","url":null,"abstract":"<div><div>In engineering, the design of bandgap mechanical metamaterials is generally required to satisfy the necessary static properties, which are macroscopically characterized by an elastic tensor. In this work, a design strategy is developed for stretch-dominated mechanical metamaterials that need to simultaneously satisfy the desired elastic tensors and bandgaps by modelling them as pin-jointed bar structures. Considering that the elastic tensor of a metamaterial is determined by the stiffness matrix of its unit cell, an analytical relationship between them is established using the Cauchy‒Born rule. Based on the analysis of the stiffness matrix, a sensitivity expression of the elastic tensor with respect to the element cross-sectional areas and the node coordinates of the unit cell is established, which can be used to design two types of unit cell parameters for the desired elastic tensor. This sensitivity matrix has a broad null space because the number of unit cell parameters to be designed is generally much larger than that of the desired elastic tensor components. Using the basis vectors of this null space, these design parameters can be modified to create the desired bandgap while keeping the elastic tensor components constant. Even if the elastic tensor components deviate during numerical iteration, they can be corrected by adjusting the design parameters in the null space of the bandgap sensitivity matrix. Thus, a two-step strategy is proposed to design the elastic tensor and the bandgap sequentially. The proposed method is verified by an example in which the desired elastic tensor components of a three-dimensional metamaterial are designed and two examples in which both the desired elastic tensor components and bandgaps are achieved. The elastic tensors of the obtained metamaterials are verified by comparing the deformations of the metamaterials and those of their equivalent continuum finite element models, and the bandgaps are also confirmed by the frequency–response curves.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"324 ","pages":"Article 113683"},"PeriodicalIF":3.8,"publicationDate":"2025-09-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145218368","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Finite axisymmetric deformations of liquid-amplified electrostatic zipping structures for shape displays","authors":"Shengjun Fan ,&nbsp;Fei Jia ,&nbsp;Jianglong Guo ,&nbsp;Yanju Liu","doi":"10.1016/j.ijsolstr.2025.113667","DOIUrl":"10.1016/j.ijsolstr.2025.113667","url":null,"abstract":"<div><div>Active shape morphing surfaces, that generate shape geometries through actuation, are an expanding field of soft robotics and have great potential for applications in haptic feedback, information displays, aerodynamics, and more. Liquid amplified electrostatic zipping actuation is a favorable candidate for morphing surfaces due to its fast response, low power consumption and considerable force output. Here, we propose an axisymmetric liquid amplified electrostatic zipping structure, containing a convex rigid electrode and an elastic film-electrode as morphing surfaces. A modeling method is developed for the deformation of the film-electrode, by combining the non-linear elasticity theory and energy minimum principle. We find numerically and experimentally that wrinkling may occur at the kink of the deformed film-electrode. A relaxed strain energy function is adopted to consider the tensed and wrinkled regions. This model allows us to better understand the impact of applied voltage, material, geometry, and even film pre-stretching on the quasi-static behavior of the film-electrode. The effects of non-dimensional voltage, rigid electrode radius and aspect ratio, and film pre-stretch are discussed. The experimental validation of this model is presented. Finally, we preliminarily demonstrate that the proposed shape morphing surface is used for information displays.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"324 ","pages":"Article 113667"},"PeriodicalIF":3.8,"publicationDate":"2025-09-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145156704","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":"Stiffness tailoring for enhanced buckling performance of isotropic panels","authors":"Sandesh Amgai,&nbsp;Paul Davidson","doi":"10.1016/j.ijsolstr.2025.113660","DOIUrl":"10.1016/j.ijsolstr.2025.113660","url":null,"abstract":"<div><div>This study investigates the buckling performance of isotropic panels with spatially varying stiffness distributions, achieved through three design strategies: thickness variation, modulus variation, and a combined approach. A semi-analytical formulation based on Classical Plate Theory and the Rayleigh–Ritz method was developed to compute critical buckling loads for arbitrary stiffness profiles. To explore the design space, multi-objective optimization was conducted using Fourier-based parameterization, producing Pareto fronts that capture the trade-off between buckling capacity and structural weight. Results show that tailored stiffness distributions can significantly improve buckling resistance—by up to 76%—compared to constant-stiffness designs of equal mass. Among the strategies, combined thickness–modulus variation outperformed individual approaches, and higher-order Fourier terms enabled finer control of stiffness gradients. The semi-analytical model was validated through finite element simulations, with discrepancies within 5%, confirming its accuracy. This work provides both a computational framework and mechanistic insights for designing lightweight, buckling-resistant structures via spatial stiffness tailoring.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"324 ","pages":"Article 113660"},"PeriodicalIF":3.8,"publicationDate":"2025-09-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145156705","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":"Forecast buckling and wrinkling analysis of Föppl von-Kármàn plates in functionally graded materials using Hermite-type approaches","authors":"Mohammed Rammane ,&nbsp;Said Mesmoudi ,&nbsp;Oussama Elmhaia ,&nbsp;Youssef Hilali ,&nbsp;Omar Askour ,&nbsp;Oussama Bourihane","doi":"10.1016/j.ijsolstr.2025.113650","DOIUrl":"10.1016/j.ijsolstr.2025.113650","url":null,"abstract":"<div><div>This work presents a numerical analysis of buckling and wrinkling phenomena in thin Functionally Graded (FG) isotropic and FG sandwich plates using the Föppl von-Kàrmàn theory. Three numerical approaches based on the use of the Asymptotic Numerical Method (ANM) with Hermite-type approximations under either strong or weak formulations are introduced. These approaches include the recently suggested Collocation Hermite-type Weighted Least Squares method (CHtWLS) coupled with ANM, the Finite Element Method (FEM) with ANM, and the Hermite-type Element Free Galerkin method (HEFG) with ANM. They are applied to study the influence of boundary conditions, plate geometry, and material distribution in FG sandwich plates on the stationary bifurcation behavior. The presented approaches employ path-following techniques to accurately and efficiently predict critical load factors and wrinkling shape modes. The use of different numerical methods adds versatility and robustness to the analysis, making the proposed CHtWLS with ANM particularly efficient and straightforward for nonlinear problems with Hermitian boundary conditions and a low degree of freedom, unlike FEM and HEFG with ANM, which are used under a weak form without the need for special treatment techniques. The buckling phenomenon is investigated in this study on square and rectangular plates using a biaxial compressive load, while the wrinkling phenomenon is investigated using a uniaxial tensile and compressive loads. The numerical findings underscore the significance of the FG power-law index as a crucial parameter influencing buckling behavior and critical loads under various boundary conditions. Notably, the study unveils an unusual wrinkling phenomenon under simply-supported boundaries, revealing a nonlinear bifurcation diagram featuring multiple critical points and folding behaviors in both square and rectangular plates.</div></div>","PeriodicalId":14311,"journal":{"name":"International Journal of Solids and Structures","volume":"324 ","pages":"Article 113650"},"PeriodicalIF":3.8,"publicationDate":"2025-09-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145156701","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}