The power battery pack of an electric vehicle is a complex system, whose efficient and high fidelity modeling is an urgent need for solving collision problems. This paper firstly proposes a hybrid reduced order battery pack model, with fine modeling in the central region and centralized mass in other parts through a simple connection contact. Compared with the refined full model, the maximum stress error is relatively small, and the computational efficiency is improved by about 3.5 times. Then, to consider the effect of the boundary further, a centroid module structural model with equivalent contact is proposed. Compared with the hybrid reduced order model, the collision displacement and the maximum stress of the centroid module structure model are much closer to those of the fine full model. Finally, an impact response analysis is conducted on the centroid module structural model and the fine full model, by comparing the behaviors of the collision center point and the impact object. It is found that both models have the similar trend in the response curves with small errors, which verifies the feasibility of the proposed model.

Damage and fracture are the main causes of structural failure, which have a significant impact on engineering safety. Crack propagation problem is also a fundamental scientific challenge that needs to be solved urgently. In this paper, the relevant theoretical basis for simulating damage and fracture, such as a fracture mechanics model, damage evolution model and numerical calculation methods, such as the finite element method, boundary element method and peridynamics theory, are introduced in the form of literature review. This paper also reviews the commonly used CAE software for structural damage and fracture analysis, including general-purpose finite element programs such as the damage and fracture analysis module that comes with ABAQUS, as well as specialized fracture analysis software, damage tolerance tools, fatigue life analysis tools, etc. The development status of some autonomous CAE software is also discussed. Finally, this paper analyzes some challenges faced by CAE software for damage and fracture simulation, and looks forward to the future development direction of domestic CAE software.

Accurately constructing the nonlinear hysteresis loop model at the bolt connection is crucial for the vibration reduction and safety performance evaluation of a satellite load-carrying structure. Traditional time-domain analysis methods of computational models require substantial time costs, and typical data-driven models struggle to construct high-precision hysteresis models. To address these challenges, a novel Residual Improvement Deep Learning Algorithm (RIDLA) is proposed for constructing the hysteresis loop model of displacement and force at the bolt connection. The algorithm fully leverages the capacity of Long Short-Term Memory (LSTM) neural networks to fit nonlinear relationships in time series. It adopts an innovative approach by creating a multi-level residual improvement deep learning model that iteratively refines predictions based on measured responses, resulting in highly accurate modeling of hysteresis at bolt connections. The performance of the RIDLA method is validated using experimental data from cyclic loading of a subcomponent of a satellite load carrying structure. The findings demonstrate that RIDLA achieves highly accurate predictions of the displacement and force hysteresis loop at the bolt connection. Additionally, the RIDLA method could be applied to predict the dynamic responses of other complex non-linear systems.

The non-isothermal complex flow caused by fluid impacting obstacles is very important to the industrial processes such as nuclear energy utilization. Through coupling various numerical techniques such as the density diffusive term, artificial viscous term, particle shifting technique, a stable and accurate non-isothermal smoothed particle hydrodynamics (SPH) scheme is established, and accurate simulation of non-isothermal complex flow caused by fluid impacting obstacles is realized. Based on the simulation for the non-isothermal flow past a heated cylinder, the non-isothermal dam break past single/multiple obstacles, it is demonstrated that: (1) the developed non-isothermal SPH scheme can not only compute a smooth pressure field and avoid the spurious oscillation of numerical solutions, but also predict accurately the temperature field and the key physical quantities; (2) this SPH scheme can also accurately show the interaction between the heat conduction process and the complex free-surface evolution, and has the capability to simulate non-isothermal complex flows past multiple obstacles.

The discontinuous Galerkin (DG) method has been widely adopted due to its excellent properties such as high accuracy and ease of parallelization. The adaptive mesh refinement (AMR) technique has been widely adopted to improve computational efficiency with much less computational cost compared with uniform global refinement to the same level with AMR. This paper combines the advantages of DG and AMR, and a new hybrid limiter is applied to the DG method on adaptive Cartesian grid based on p4est, an open-source library. The limiter exhibits advantages of high precision, compactness, robustness, and ease of implementation. The shock wave is captured with a shock indictor and the performance of the new hybrid limiter is compared with that of the total variational bounded (TVB) limiter in this paper. The result shows that the performance of the former is significantly better than that of the latter. A series of numerical examples for Euler equations and Navier-Stokes equations are used to verify the feasibility and efficiency of the proposed method. The results show that the new hybrid limiter performs very well in the AMRDG method, it has lower dissipation and great shock capture ability, and the computational efficiency is greatly improved while the accuracy is guaranteed.

Considering the issue that low frequency and broadband of elastic wave metamaterials cannot coexist, this paper realizes a low-frequency band gap by rotating asymmetric mechanical metamaterials, and further widens the low-frequency bandgap by introducing multiple orders. By utilizing the node rotation and ligament bending deformation characteristics of anti-chiral materials, the vibration body size and ligament stiffness of the anti-tetra chiral unit cell diagonal are gradually adjusted through ligament folding, and the multi-order asymmetric unit cell design is realized. The generation and change mechanism of elastic wave bandgap are explained by analyzing the resonance mode and transmission characteristics of the upper and lower bounds of the bandgap. The study shows that: in the asymmetric mode, the node rotation and ligament bending deformation characteristics of chiral materials are utilized to realize the rotation resonance of the mass block and open the bandgap, and the band gap is widened by the resonant superposition between two pairs of different mass blocks arranged alternately. Finally, the proposed asymmetric mechanical metamaterial is verified by experiments to have to demonstrate improved broadband and low-frequency vibration isolation performance.

The formation of porous media is influenced by a number of factors, including the deposition and fragmentation of particles, which result in the formation of interlayers with varying structures. These interlayers exert a significant influence on the mechanical behavior of porous media. This paper presents a systematic investigation into the influence of the inclination angle and thickness of the interlayer on the mechanical behavior of porous media, employing the discrete element method. The results demonstrate that the stress intensity of porous media containing interlayers is between those of the two homogeneous porous media and varies with changes in the inclination angle and thickness of the interlayers. The average coordination number between grains is found to be significantly affected by the thickness of the interlayer at the beginning of loading, but stabilized at the end of loading. The variation of the coordination number affects the distribution of strong and weak force chains, while the inclination angle and thickness of the interlayer determine the magnitude and direction of stress transfer in the force chains. Furthermore, the contact unit normal force and normal contact force are deflected with the increase of the inclination angle of the interlayer, demonstrating significant anisotropy. This study advances our understanding of the intricate mechanical behavior of porous media containing interlayers in strata, offering invaluable insights for optimization and practical application in geological engineering.

This paper presents an analytical method, namely interface stiffness transfer method, for evaluating the responses of multilayered elastic structures. Based on the Love function and general solutions, the stiffness matrix relationship of the displacement-stress state vectors is introduced to obtain the interface stiffness transfer matrix equation between adjacent layers, which satisfies an algebraic Riccati matrix equation. When the elastic layer is a half-space, an explicit solution is obtained directly for the interface stiffness matrix. The interface stiffness transfer matrix method starts from the bottom layer with a known stiffness, and then deals with one layer at a time until the uppermost layer is reached, obtaining the interface stiffness of the multilayered structure. Finally, by solving the symmetric equilibrium equations of the boundary conditions, the displacement-stress state vector of an arbitrary layer is obtained. This method keeps the advantages of the classical transfer matrix method, but naturally excludes its exponential growth terms. In particular, the proposed method is a powerful candidate for efficiently solving the algebraic Riccati equation for the optimal control problems. Numerical examples show the properties of the interface stiffness transfer method.

A surrogate model is a new research idea direction of aerodynamic data generation. A traditional surrogate model relies on a large number of high-precision simulation sample points and their responses values to ensure the accuracy of the model. A multi-credibility surrogate model can reduce the computational cost while maintaining a certain accuracy by integrating multi-layer high and low credibility models, which is of great significance for reducing a missile development cycle. In this paper, the influence of different numbers of low-confidence sample points on the multi-confidence model and the optimal ratio of high-low-confidence sample points are studied, and a multi-confidence sampling method suitable for aerodynamics data is proposed. It is applied to the construction of three multi-credibility surrogate models in the prediction of missile aerodynamic data, among which the Co-Kriging model has the best comprehensive prediction effect. The recommended ratio of high and low confidence sample size is between 1: 4 and 1: 3.

Accurately simulating the impact characteristics of a dam break is of paramount significance for the prediction and mitigation of dam-break flow disasters. The B-spline material point method (BSMPM), as an improved algorithm of the material point method (MPM), effectively enhances computational accuracy and improves convergence. However, the BSMPM solves the governing equations based on a tensor grid rather than an Eulerian background grid. Moreover, its interpolation shape functions have a larger influence domain. Consequently, when solving problems involving fluid-structure coupling and contact, issues such as premature contact, difficulty in capturing contact interfaces, and challenges in calculating contact forces arise. Within the same tensor grid space of the BSMPM, accurate capture of contact interfaces is achieved based on the relative velocities and unit outward normals of the same nodes; employing the Greville Abscissa enables precise contact of contacting objects, thereby avoiding premature or spurious contact; through the Lagrange multiplier method, interface contact forces are accurately determined, thus constructing a high-precision contact algorithm for fluid-structure strongly coupled problems, facilitating research on the simulation of dam-break fluid impact with rigid and elastic obstacles, and enabling a comparison with existing experimental or simulated results. The results demonstrate that the simulated impact loads and structural deformation evolution patterns correspond well with existing experimental/simulated results. For rigid obstacles, the peak impact pressure exhibits concave parabolic growth and positive correlation exponential function growth with increasing water level and dam-break slope, respectively. For elastic obstacles, the peak impact pressure decreases exponentially with the increase in the height of the probing point. The feasibility and effectiveness of simulating dam-break flow impact problems using the BSMPM contact algorithm are validated, providing a new perspective for simulating dam-break flow impact problems.