The microscale effects of non-Fourier heat transfer are often ignored in studies concerning thermal shocks. This paper presents a one-dimensional physical model representing the composite structure of a flat plate coating and substrate. Model I considers the hyperbolic heat transfer of the coating and the parabolic heat transfer of the substrate. Additionally, appropriate boundary conditions are determined based on the heat transfer behavior at the interface. On this basis, a thermoelastic mechanics model of the coating and substrate is formulated. The model is discretized using the implicit difference method to acquire the numerical solution for the temperature field. Subsequently, the stress field is determined, and specific examples are provided. At the same time, mathematical model II of parabolic heat transfer for both the coating and substrate is established for comparative study. It is found that model I demonstrates delayed change, localized distribution, and fluctuation of thermal stress within the coating when taking into account the microscale effect of non-Fourier heat transfer, assuming identical initial conditions and thermal perturbations. In model I, the thermal stress at any position does not start from zero. Conversely, model II shows no fluctuation, and the thermal stress at any position starts to change from zero. After the generation of thermal stress of model I, it reaches the peak first, and the peak value is larger than that of model II. In the substrate, the thermal stress of model I is larger than that of model II, and the gradient of change is higher. At the interface, model I exhibits a “reflection effect”, where the stress value and the stress drop are larger than those of model II. The comparison shows that the thermal shock to model I is more complicated and intense. This study provides a useful reference for ensuring the reliability of coatings under extreme heat transfer environments.

In this study, the angle-preserving transformation method is employed to establish a propagation model for I/II composite lip-shaped cracks under tensile loading conditions. Based on Irwin's small-scale yielding equivalent hypothesis, a plastic propagation zone model is formulated for Ⅰ-Ⅱ composite lip-shaped cracks under tensile loading. This model provides expressions for the stress intensity factors (SIFs) of mode I and mode II at the tip of lip-shaped cracks within the plastic zone. Additionally, the stress distribution along the extension line of the lip-shaped crack tip is characterized. A tensile simulation model is developed, and the theoretical solution for stress distribution at the lip-shaped crack tip is compared with the elastoplastic and linear elastic simulation results. It is found that, based on Irwin's small-scale yielding equivalent hypothesis, the modified dimensions of lip-shaped cracks lead to increased crack sizes and greater equivalent SIFs. Geometric alterations in lip-shaped crack parameters also influence the plastic zone, with larger semi-lengths resulting in larger plastic zones under equivalent width-to-length ratios. Conversely, greater width-to-length ratios lead to smaller plastic zones under equivalent semi-lengths. Moreover, an increase in the inclination angle of the lip-shaped crack corresponds to a proportional increase in the plastic zone size. The plastic correction theory at the lip-shaped crack tip, founded on the basis of Irwin's small-scale yielding equivalent hypothesis, aligns well with plastic finite element simulations. As the inclination angle of the lip-shaped crack rises, stress levels at the crack tip diminish. On the one hand, this phenomenon arises from the transition from mode I crack extension to Ⅰ-Ⅱ composite crack extension, coupled with stress yielding at the concave region of the lip-shaped crack for larger inclination angles. On the other hand, this stress yielding serves to mitigate stress concentration at the crack tip, ultimately resulting in reduced stress levels at the crack tip.

Bounds on the mechanical properties provide fundamental guidelines for finding materials or structures with extreme mechanical performance. However, the bounds on some important mechanical properties, such as Young's modulus and tensile strength, remain unknown, while the search for target extreme materials from infinite potential materials of element combinations across the periodic table is challenging. It has long been questioned: have we approached the bounds on these mechanical properties? Is there a material that is stiffer or harder than diamond? To determine the bounds on the mechanical properties and find materials or structures with extreme mechanical performance, the key is to understand and quantify the structure-property relationship. Over the past decades, many attempts and achievements have been made to model the structure-property relationship, such as empirical/semiempirical formulas, first-principles calculations, machine learning, but these approaches often suffer from limitations in terms of accuracy, efficiency, universality, or interpretability. With the accumulation of knowledge and data, knowledge and data-driven understanding and modeling of structure-property relationships have shown immense potential. Recent studies within the knowledge and data-driven framework have led to the development of powerful theories for structure-property relationships. Based on these structure-property relationships, material properties can be predicted from structures, and conversely, structures can be designed for target material properties. Consequently, the bounds on some important mechanical properties have been determined, and numerous materials or structures with mechanical properties close to the theoretical bounds have been designed and fabricated. Our work provides an overview of the recent progress in these explorations of bounds on mechanical properties. First, we present the advances in knowledge and data-driven approaches for understanding and modeling structure-property relationships. Then, we review the determined bounds on mechanical properties and discovered materials or structures with extreme mechanical performance based on the knowledge and data-driven approaches. Finally, we discuss the challenges, opportunities, and some future directions in this field.

Shell structures are widely used in engineering, especially in aerospace and civil engineering. As a result, the study on dynamic behaviors of shell structures is crucial for engineering applications. Over the years, shell theory has undergone continuous improvement and development, leading to analytical solutions for specific shell structures. However, solving analytical solutions for complex shell structures with intricate shapes becomes highly challenging and even unattainable. Therefore, numerical methods like the finite element method (FEM) and the meshless method are employed for further solutions. The meshless method is a powerful complement to FEM, relying solely on nodal information for the formation of shape functions and enabling easy construction of higher-order smooth approximations. Consequently, it naturally holds an advantage in the numerical analysis of plate shell structures. Based on the 3D continuous shell theory and the moving least-squares (MLS) approximation, a meshless model for arbitrary shells is established in this paper. The MLS approximation is used not only for geometric surface interpolation, but also for displacement field approximation. The meshless equation governing the forced vibration of arbitrary shells is derived under Hamilton's principle and solved using the time-domain implicit Newmark method. The full transformation method is used to impose the essential boundary conditions. The code for the proposed method is developed in the MATLAB platform and used to compute several representative shell examples, obtaining the first ten natural frequencies of each shell type and the time history deflection at the center point under different pulse loads, considering both with and without damping. The calculated results are compared with ABAQUS solutions to verify the effectiveness and accuracy of the presented method. When using the meshless method to solve the forced vibration of arbitrary shell structures, it does not rely on grid partitioning. This enables effective analysis of various shell structures with different shapes, showcasing its strong adaptability and vast potential for applications.

Mechanical metamaterials (or meta-structures) exhibit extraordinary physical and mechanical properties due to their unique microstructural designs. By combining the design ideas of mechanical metamaterials with intelligent and flexible materials (IFMs), it is possible to create intelligent flexible mechanical metamaterials (IFMMs) with self-sensing and self-actuating capabilities. This paper takes conventional mechanical metamaterials as a starting point and analyzes the fundamental design ideas, deformation mechanisms, and mechanical properties of IFMMs. According to recent research progress on IFMMs, this novel metamaterial is categorized as mechanical metamaterial based on shape memory polymers (SMPs), hydrogel, magnetoactive soft materials, and dielectric elastomers, with a particular focus on the first two types. On the basis of our previous work, we present a general approach that utilizes analytical methods and numerical simulations to analyze the mechanical properties of negative Poisson's ratio, negative expansion, and multi-stable metamaterials under the assumptions of small deformation and large deformation with multi-field coupling, respectively. In the case of small deformation, the use of beam theory and energy methods proves to be essential for obtaining fundamental mechanical parameters of the materials. Moreover, accurate constitutive models and numerical implementation under large deformation and multi-field coupling offer the possibility to analyze more complex deformations and structures. In addition, the preparation and performance testing of IFMMs remain crucial. Advanced manufacturing techniques have introduced new opportunities for the preparation of IFMMs, and currently, various methods are available to effectively prepare these materials. The performance testing of IFMMs includes both experiments applicable to traditional materials and specialized experiments only for IFMs. Finally, this paper concludes by highlighting some key issues and potential trends of IFMMs. These challenges primarily revolve around material properties, fabrication methods, mechanical models, and structural designs. This review may bring beneficial inspiration for the future development of IFMMs.

During the manufacturing process of a rough workpiece, the non-uniformity of material mechanical properties can result in the generation of residual stresses within the workpiece, which can lead to structural failure. During the cutting removal process of the workpiece, the residual stresses gradually release, causing deformation. In this study, the "birth-death element" technique of finite element analysis was used to simulate the cutting removal process of the material. This process was then transformed into the release of residual stresses. A novel analytical method combining radial basis function interpolation with geometric and physical equations was proposed based on plate shell theory, small-deformation theory, theory of elasticity, and the superposition principle. The method aimed to invert the residual stress field and calculate the resulting deformation. This study was divided into two parts. In the first part, the influence of stress equilibrium equations was neglected, and the radial basis function interpolation method was used to invert the release of residual stresses in thin plates according to the initial residual stress field and the residual stress field after material removal. Next, the stresses were substituted into the physical equations to calculate the strain. The strain was then substituted into the geometric equations, and the plane displacement was calculated by strain integration from geometric equations. Based on the plate shell theory equations, the buckling deformation was inverted according to the plane displacement. In the second part, the reverse process of the first part was performed. Firstly, the buckling deformation caused by material removal in the thin plate was inverted. Then, the buckling deformation was substituted into the plate shell theory equations to invert the release of residual stress and the reconstructed residual stress field. The results demonstrated the reversibility of these two processes. Furthermore, the analytical solutions showed high agreement with the finite element solutions. This suggested that the analytical method proposed in this paper is applicable to thin-plate structures under elastic conditions and can be used to estimate residual stress distribution and predict deformation in the thin-plate cutting removal process.

In this paper, the filter function in the ICM method and the penalty function in the variable density method are both referred as the mapping functions. Different forms of mapping functions have a significant impact on the convergence efficiency of structural topology optimization. Therefore, it is necessary to study how to construct a suitable mapping function for the optimization model. Aimed at this problem, how to construct and select a mapping function in the establishment of the structural topology optimization model is studied, and the influence of different mapping functions on the convergence efficiency of structural topology optimization is discussed. An approach is proposed to construct a mapping function to achieve high-efficiency convergence in structural topology optimization. Five common forms of mapping functions are also given. An optimization model and a solution algorithm matching the mapping function with highly efficient convergence (MFHEC) are proposed. Firstly, the convergence rates of the filter function and the quasi-filter function of the same form of mapping functions are compared. Then the convergence rates of the fast filter function of different forms of mapping functions are compared. Taking the structural topology optimization problem of minimizing structural volume under displacement constraints as an example, the ICM method is adopted to establish the models and solve the problems. The higher convergence efficiency of MFHEC is verified by the results of numerical comparison. The results show that the fast filter function has a faster convergence rate than other functions in the same form of mapping functions. Compared with five different forms of mapping functions, the filter function of power function form has the fastest convergence efficiency. Finally, it should be emphasized that the conclusions of the mapping function studied in this paper are equally applicable to the filter function of the ICM method and the penalty function of the variable density method. The proposed method for constructing MFHEC is very useful for improving the efficiency of the ICM method and the variable density method.

Defects play a crucial role in understanding the physical and mechanical behavior of materials. In this study, the fracture problem of an infinite one-dimensional hexagonal piezoelectric quasicrystal material matrix containing secondary asymmetric straight cracks with lip-shaped pores is investigated. A defect mechanics model of secondary asymmetric cracks with lip-shaped pores is constructed for the first time. Utilizing conformal transformation technology, a conformal transformation formula from an infinite region containing secondary asymmetric cracks at the lip on the physical plane to the outer region of the unit circle is built. Using the complex variable method, analytical expressions for the field intensity factor and energy release rate at the crack tip are obtained. Under given conditions, these analytical results can be simplified into solutions for other defect models, such as secondary single cracks at the lip and secondary symmetric cracks at the lip. At the same time, they can also degenerate into the solutions of classical Griffith cracks and lip cracks without secondary cracks. Numerical examples reveal the effects of defect size, particularly the lip height and crack length, on the field intensity factor and energy release rate. The results show that increasing the length of both sides of the crack promotes crack propagation, while increasing the height of the lip inhibits crack propagation. These findings are consistent with the conclusions drawn from theoretical analysis. When the length of the secondary crack on one side of the lip is zero, as the height of the lip increases, the stress intensity factor and energy release rate at the crack tip on the other side first increase to a peak and then gradually decrease, eventually stabilizing at a constant level. The research results presented in this paper can contribute to the development of a theoretical framework for material fracture mechanics and provide technical assistance for nondestructive testing, reliability design, and optimization of piezoelectric quasicrystal material equipment and components.

During penetration, the temperature of the projectile will sharply rise due to the large amount of heat generated by the sliding friction between the projectile and the target. High temperature can soften the projectile, potentially change its shape, penetration mechanism, and further affect the penetration ability. In order to study the temperature rise of the projectile during high-speed penetration, a two-part temperature rise calculation model of the penetration is established. In the first part, the heat flux data set of the projectile at different positions in the process of penetration is obtained according to rigid body dynamics theory and the friction heat generation mechanism. The second part takes the heat flux data set as the boundary condition and calculates the temperature distribution of the projectile based on heat conduction theory and the finite difference algorithm. The stability of the model is discussed from the two aspects of time step and projectile mesh, and appropriate values are selected. The calculation model of temperature rise is used to study the heat flux and temperature distribution of the projectile, and the factors influencing the temperature rise of the projectile are discussed and analyzed. The results show that the temperature rise is obvious during high-speed penetration, but it only lasts for a very short time, and the high temperature is mainly distributed near the surface of the projectile head. The position of the highest surface temperature outside the projectile during penetration is related to the shape of the projectile. During the penetration time, the ratio of the heat conduction distance to the radius of the projectile decreases with the increase of the size of the projectile. The research results are useful for the design and material selection of high-speed penetration projectiles.

Electrical connectors widely used in electrical and electronic devices suffer from severe contact failure problems which determine the reliable service of these devices or even cause their destruction. The insertion and withdrawal forces are significant factors to evaluate the performance and quality of electrical connectors. The wear behavior due to the repeated insertion and withdrawal processes has a non-negligible influence on the forces. Therefore, the forces and wear behavior are investigated in this paper. A simulation model is developed to study and evaluate the performance of electrical connectors. The finite element method (FEM) is used to simulate the insertion and withdrawal processes and the wear profile of an electrical connector in the application. According to the development of the FEM model of the electrical connector, the equivalent insertion and withdrawal forces with and without the consideration of wear behavior are calculated. Frictional behavior is used to study the insertion and withdrawal forces without considering wear. In cooperation with the adaptive mesh technique of the commercial FEM software ABAQUS, the wear model based on frictional dissipation energy is applied to predict the wear morphology. The characteristics of the insertion and withdrawal procedure are analyzed using the forces vs. time curves. Furthermore, the wear profiles of the contact surface are obtained at different times of the insertion and withdrawal processes, and the effect of wear on the insertion and withdrawal forces of the electrical connector is discussed. The critical number of insertions and withdrawals leading to a larger wear depth is found. The size parameter of the contact component is modified to validate the conclusion of wear influence. It is worth emphasizing that this study is significant to the failure problem of electrical connectors and provides a guideline for the production and application of electrical connectors.