International Journal of Solids and Structures最新文献

In this paper, we investigate wave propagation in cubic periodic architectured materials. We analyse three different types of unit cells, with distinct symmetries (centrosymmetric, non-centrosymmetric chiral and non-centrosymmetric achiral) with the aim of investigating the consequences of such symmetries on the elastodynamic behaviour of the architectured material. To this end, numerical simulations are performed on unit cells representative of the three types, to extract phase velocities and polarisations of waves along different directions. It is shown that some unconventional couplings between the different eigensolutions give rise to circular or elliptically polarised waves, associated with dispersive effects (acoustical activity). Subsequently, a theoretical analysis using a generalised equivalent continuum model (strain gradient elasticity) is performed to analyse these results and unveil the links between the symmetries of the architecture and the macroscopic elastodynamic behaviour. Indeed, it is shown that strain gradient elasticity is able to discriminate between the three symmetry classes, that are seen as equivalent by a classic continuum theory.

Peeling the front-side film from the flexible and ultra-thin wafer is a critical procedure for the fabrication of ultra-thin chips. For a successful peeling process, the following conditions are required simultaneously: the interface between the film and the wafer is debonded, the interface between the wafer and the substrate remains undelaminated, and the wafer stays intact. However, there are relatively few studies focusing on the underlying mechanism in this peeling process. Here, a theoretical model is developed to investigate the competing behavior of interface delamination and wafer cracking for the bilayer film-substrate system. Based on the constant-stress (Dugdale) cohesive law and Euler-Bernoulli beam theory, both the competing interface delamination criterion and the wafer cracking criterion are determined. The corresponding competing maps of interface delamination and wafer cracking are obtained, in which the interface delamination path and the wafer safety status can be predicted. The effect of several dimensionless parameters on the competing behavior of interface delamination and wafer cracking is examined systematically, including the property of the geometry, the material, and the interface of the bilayer film-substrate system. The theoretical model is validated by both finite element analysis (FEA) and experimental results. This research aims to provide some guidance for optimizing the peeling parameters and contribute to a higher success rate of peeling process.

This paper presents a topology optimization framework utilizing a deformation plasticity model to approximate the isotropic hardening von-Mises incremental elastoplasticity model under monotone proportional loading. One advantage of the model is that it is based on a yield surface allowing for precise matching to uniaxial elastoplastic isotropic hardening response. The deformation plasticity model and the incremental plasticity model coincides for proportional loading and since the deformation plasticity model is path-independent, the computational cost and implementation complexity reduce significantly compared to the conventional incremental elastoplasticity. To investigate the deformation plasticity model combined with topology optimization, we compare three common elastoplastic optimization objectives: stiffness, strain energy and plastic work. The possibility to limit the peak local plastic work while maximizing the strain energy is also investigated. The consistent analytical sensitivity analysis which only requires the terminal state is derived using adjoint method. Numerical examples demonstrate that the proportionality assumption is reasonable and the deformation plasticity model combined with topology optimization is a competitive alternative to cumbersome incremental elastoplasticity.

Additive manufacturing (AM) of high-temperature alloys through processes such as laser powder bed fusion (LPBF) has gained significant interest and is rapidly expanding due to its exceptional design freedom, which enables the fabrication of complex parts that contribute to the increased efficiency of aerospace and energy systems. The materials produced through this process exhibit unique microstructures and mechanical properties, which necessitate dedicated study and characterization. In this context, our research focuses on the experimental characterization of the isothermal cyclic viscoplastic mechanical response of Hastelloy X (HX) over the temperature range of 22 to 1000 °C and at various strain rates, addressing a current gap in the literature. Recognizing the need for material models that can accurately represent the cyclic mechanical response of LPBF HX across a broad temperature range, we developed a robust extension of the viscoplastic isotropic-kinematic hardening Chaboche model, intended for applications in the thermomechanical simulation of the LPBF process for the analysis of residual stress and distortion, as well as for assessing the mechanical integrity of LPBF components. The extension involves expressing the entire set of model parameters explicitly with analytical functions to account for their temperature dependence. Consequently, the model includes a relatively large number of parameters to represent the isotropic-kinematic hardening viscoplastic response of the alloy over a wide temperature range, and hence to overcome the endeavor of its systematic calibration, a dedicated calibration approach was introduced. The model ultimately demonstrated its capability to precisely represent the isothermal response of the alloy over the examined temperatures and strain rates. To evaluate the model’s predictiveness for non-isothermal conditions, out-of-phase thermomechanical cyclic experiments were also conducted as independent benchmark tests, where the model’s predictions were fairly consistent with the experimental results. As a part of this study, the derived material model has been integrated into the UMAT subroutine, complete with an analytical derivation of the consistent Jacobian matrix.

A stable deformation mode is highly desired for mechanical metamaterials, especially when coupled with a negative Poisson’s ratio. However, such metamaterials often face challenges in terms of scalability toward large deformation or strain. In response, we propose a multi-step hierarchical auxetic metamaterial design paradigm, incorporating a series of incrementally scaled-down structures with same scale factor $\alpha$ into a re-entrant framework. This design enables instability regulation and multi-step deformation capabilities while preserving auxetic behavior, even under significant strain. Such multi-step metamaterials exhibit excellent properties, including tailored multi-phase compression modulus and strength, along with an enhanced energy absorption capacity that is as large as 2.1 times that of the original auxetic metamaterial. Experiments and simulations demonstrate that the deformation mechanism and compression response of the proposed multi-step auxetics are strongly influenced by the reduction factor and the order of the inner structure. A particularly intriguing observation is that the incorporation of embedded microstructures can restore stable deformation, even in the presence of significant initial instability, particularly with a reduction factor of $1/5$. At high relative density, its specific energy absorption stands out favorably compared to other configurations, highlighting the success of the recoverable buckling mechanism. This work paves the way for designing multi-step mechanical metamaterials for use in impact resistance and body protection.

To understand the influence of friction on the shear-slip behavior of heterogeneous brittle composites, a novel cohesive interlayer model that can effectively capture the friction effect was proposed based on the classical Park-Paulino-Roesler model. Meanwhile, the unified potential energy function governing the interface tangential and normal behaviors was introduced to realize the mechanical interaction between Mode I fracture and Mode II fracture, and a smooth friction growth function was added in the elastic deformation stage for calculating the accurate contact pressure and friction force. Furthermore, the capability of the proposed model in addressing unloading and reloading was improved, and the fracture energy can vary accordingly during cyclic loading. To verify the effectiveness of the proposed model, it was examined by modelling the shear behavior of a masonry wallette. The results show that the relative error of the proposed model is 14.92% which is much lower than those of the other three pre-existing models when calculating the displacement corresponding to peak shear stress. Meanwhile, in terms of peak shear stress and initial displacement at residual stage, the relative errors of the proposed model are only 1.82% and 5.04%, respectively, indicating the high accuracy. Besides, the tangent stiffness determined by the second-order integration of the potential energy function is also continuous and smooth, which ensures the effective convergence of the proposed cohesive model.

In this study, we investigate the effect of Van der Waals and Casimir forces on the mathematical model of nano-electromechanical systems (NEMS) such as nano-beam actuators that contain cantilever and double cantilever beams. The singular nonlinear boundary value problem governing the beam-type actuators, including geometric nonlinearity is solved by using an intelligent strength of feedforward artificial neural networks (ANNs) and hybridization of optimization algorithms such as arithmetic optimization algorithm (AOA) and active set algorithm (ASA). The proposed ANN-AOA-AS algorithm is employed to quantify the effect of changes in applied voltage, dispersion forces, geometric nonlinearity parameters, and initial axial strain on the deflection of the beam. Furthermore, to validate the results obtained by the proposed algorithm, statistical analyses are conducted to compare the approximate solutions with state-of-the-art methodologies available in the latest literature. In addition, performance indicators are defined such as mean square error (MSE), Nash–Sutcliffe efficiency (NSE), mean absolute deviations (MAD), root mean square error (RMSE), and Error in Nash–Sutcliffe efficiency (ENSE) to study the accuracy and efficiency of the solutions. The results show that these indicators’ mean percentage values lie around $10^{-4}$ to $10^{-6}$ which reflects the perfect modeling of the approximate solutions.

An artificial intelligence machine learning-based design framework is proposed to design lattice-based metamaterials with hexagonal symmetry that deliver wide band gaps at user-desired frequency ranges between 0 and 1000 kHz. The design approach starts by selecting a traditional, easy-to-manufacture parent lattice-based material that does not necessarily exhibit wide or functional band gaps. Subsequently, the parent lattice is transformed into a band-gap-rich lattice by superposing periodic triangular-shaped perturbations (i.e., zigzag-sine-based curvatures) with controllable frequencies and magnitudes on its ligaments. Finally, the frequency and magnitude parameters needed to deliver a specific band gap between 0 and 1000 kHz are determined using a hybrid intelligent framework, developed based on an Adaptive Neuro-Fuzzy Inference Systems (ANFIS). The ANFIS network integrates fuzzy logic expert models and artificial neural networks’ machine learning capabilities. Such a hybrid network is known for its ability to model strongly nonlinear and complex data. The data used in training the ANFIS models is generated using parametric finite element-based simulations where band gaps corresponding to a wide range of perturbation frequencies and magnitudes are computationally determined. The parametric study showed a nonlinear and complex topology-band gap characteristic relation; however, the Adaptive Neuro-Fuzzy Inference System (ANFIS) proved capable of modeling the observed complex topology-band gap behavior efficiently. The accuracy of the ANFIS models exceeded 99 % in several design ranges (i.e., perturbation parameters ranges). These were designated as high-accuracy design regions and were highlighted in the proposed design approach. Using multiple case studies with different band gap requirements, the ANFIS-based design framework proved effective in delivering customized lattice-based metamaterials with user-defined band gap frequencies.

This paper extends the reduced Yld2004 (rYld2004) function to present the anisotropic hardening behavior for body-centered cubic and face-centered cubic metals under the proportional loading conditions based on neural network. The parameters of the rYld2004 anisotropic hardening model (AH_rYld2004) are determined by the uniaxial tensile yield stresses along ${\text{0}}^{°}$, ${\text{15}}^{°}$, ${\text{30}}^{°}$, ${\text{45}}^{°}$, ${\text{60}}^{°}$, ${\text{75}}^{°}$ and ${\text{90}}^{°}$ from the rolling direction as well as equibiaxial tension. The evolution of anisotropic parameters are described by the back propagation neural network optimized by ant colony optimization algorithm. The predicted data by AH_rYld2004 and some common anisotropic models are compared with the experimental results to verify the precision of the AH_rYld2004 in characterizing anisotropic hardening. The comparison proves that the AH_rYld2004 precisely characterize the anisotropic evolution with increasing plastic deformation for AA 3003-O and QP980. Simultaneously, the AH_rYld2004 function based on neural network is used to accurately simulate of circular cup deep drawing for AA 3003-O and uniaxial tension for QP980. The results indicate that the AH_rYld2004 model is capable to accurately represent the plastic anisotropic evolution for uniaxial tension along seven loading directions and equibiaxial tension.

We investigate the propagation of harmonic flexural waves in periodic two-phase phononic multi-supported continuous beams whose elementary cells are designed according to the quasicrystalline standard Fibonacci substitution rule. The resulting dynamic frequency spectra are studied with the aid of a trace-map formalism which provides a geometrical interpretation of the recursive rule governing traces of the relevant transmission matrices: the traces of three consecutive elementary cells can be represented as a point on the surface defined by an invariant function of the square root of the circular frequency, and the recursivity implies the description of a discrete orbit on the surface. In analogy with the companion axial problem, we show that, for specific layouts of the elementary cell (the canonical configurations), the orbits are almost periodic. Likewise, for the same layouts, the stop-/pass-band diagrams along the frequency domain are almost periodic. Several periodic orbits exist and each corresponds to a self-similar portion of the dynamic spectra whose scaling law can be investigated by linearising the trace map in the neighbourhood of the orbit. The obtained results provide a new piece of theory to better understand the dynamic behaviour of two-phase flexural periodic waveguides whose elementary cell is obtained from quasicrystalline generation rules.