Latest ArticlesInterface consistency and error convergence are central issues in concurrent multiscale computational methods, particularly critical for atomistic-to-continuum coupling models. However, existing theoretical studies remain limited and are mostly confined to one-dimensional settings. This work focuses on the multiresolution molecular mechanics (MMM) approach and systematically investigates the impact of various energy sampling schemes on interface consistency and error convergence. Two-dimensional square and triangular lattice models containing both atomistic and coarse-grained regions are constructed under bilinear element interpolation. The results show that interface secondary sampling schemes can significantly improve consistency in the interfacial region, with the scheme incorporating all neighboring layers achieving the best performance. Error analysis reveals that discretization error dominates the total error, and increasing the number of secondary sampling points effectively reduces the sampling error, particularly under tensile loading conditions. Moreover, both lattice types exhibit consistent error convergence behavior, demonstrating high generality of the method to different structures. This study highlights the advantages of energy sampling strategies in improving interface treatment and convergence behavior in MMM, providing theoretical support for the development of high-accuracy multiscale computational mechanics methods.
To reveal the instability failure mechanism and energy evolution law of the rock mass with multiple cracks, a numerical model of red sandstone was established by using PFC2D. The mesoscopic parameters of the numerical model were calibrated based on the results of uniaxial compression tests and Brazilian splitting tests of intact red sandstone specimens. On this basis, the particle flow simulation tests of red sandstone with multiple cracks were carried out. The results show that with the increase of λ, the peak strength of the multiple-cracked red sandstone gradually decreases when the inclination angle remains unchanged; with the gradual increase of the inclination angle α, the peak strength of the multiple-cracked red sandstone gradually increases when the short-long axis ratio λ remains unchanged; the failure mode shows a diagonal tensile-shear failure, with tensile as the main and shear as the auxiliary, and the failure and instability mode of the specimens are all along the extension direction of the prefabricated cracks; the failure mode of the multiple-cracked red sandstone is jointly affected by the short-long axis ratio and the inclination angle; before the peak strength, the rock mainly shows the energy accumulation characteristics; at the peak strength, the total energy of the red sandstone is mainly elastic energy and supplemented by dissipated energy; at the instability failure, the total energy of the red sandstone is mainly dissipated energy and supplemented by elastic energy; the energy storage capacity and failure difficulty of the red sandstone change with the variation of the multiple-crack inclination angle α and the working conditions, which can be used as a reference for rock breaking operations; and a high-precision damage constitutive model of red sandstone based on different height-diameter ratios λ and different joint inclination angles α was established.
This study focuses on the integrated manufacturing process of sensors in unbonded flexible pipes for deep-sea mining engineering, investigating the sensor integration process with compensation reinforcement layers and their mechanical performance. A simulation model for the winding process of the compensation reinforcement layer in flexible pipes was established, systematically analyzing the influence patterns of winding tension and winding angles on sensor integration processes. Additionally, the effects of sensor quantity and winding angles on the tensile performance of integrated compensation reinforcement layers were explored. The results demonstrate that winding tension significantly affects sensor elongation rates, while winding angles predominantly influence stress distribution in the lining layer. The tensile performance of integrated compensation reinforcement layers shows minimal sensitivity to sensor quantity, and variations in winding angles within a specific range exhibit limited impact on anti-tensile properties. These findings reveal the correlation between sensor integration parameters and structural performance of flexible pipes, providing theoretical guidance for optimizing sensor integration processes in deep-sea flexible pipes.
The free vibration analysis of membranes is of significant importance in engineering structures, especially in the design and optimization of membrane structures. This paper presents a new type of triangular element, aiming to improve the computational accuracy in free vibration analysis of membranes. Traditional 3-node triangular elements in membrane vibration analysis typically rely on polynomial shape functions, but this method often suffers from insufficient accuracy in complex vibration modes and high-order frequencies. To address this issue, this paper constructs a 10-node triangular element with the shape function incorporating trigonometric functions. The proposed 10-node triangular element consists of three corner nodes, two points of trisection for every edge, and the centroid node with its shape functions derived using the area coordinate method. The stiffness matrix and mass matrix are derived, and the frequencies and modes for free vibration of the membrane are computed, thereby the dynamic characteristics can be studied. To evaluate the effectiveness of this element, several typical examples are chosen, including the free vibration analysis of rectangular membrane and triangular membrane. By comparing with theoretical solutions and the 3-node element calculations in Ansys, the obtained results show that the 10-node triangular element can approximate the theoretical solutions with few computational elements. And the precision of the presented 10-node element is similar with that of the standard 10-node triangular element. The high precision of the proposed element is demonstrated in analysis of free vibration of membrane, which has the potential of further research and promotion.
A study is conducted on the nonlinear dynamics of inclined flow pipes reinforced with graphene composite materials conveying fluid under two-phase flow conditions. Based on the von Karman nonlinear strain-displacement relationship and Hamilton's principle, the dynamic equations for inclined two-phase flow pipelines are established. The nonlinear dynamic model was solved using Galerkin method and fourth-order Runge-Kutta method to analyze the influence of the distribution pattern, weight fraction, and gas volume fraction of graphene platelets on the natural frequency and nonlinear dynamics of pipes conveying fluid. The experimental results show that the vibration amplitude of graphene reinforced pipes with V-shaped distribution is the smallest, followed by X-shaped distribution and A-shaped distribution. In addition, increasing the gas volume fraction can help alleviate fluid-induced vibration phenomena in pipelines. The above conclusion provides a theoretical basis for the application of graphene-reinforced inclined two-phase flow pipelines in practical engineering.
This study systematically investigates the influence of the wall slip length (Ls) on the statistical properties and flow structures of turbulent channel flow. The Navier slip boundary condition is applied at the boundary of turbulent channel flow, and direct numerical simulation (DNS) is employed to numerically explore the evolution of turbulence for Ls ranging from 0 (no-slip) to 0.1. The results reveal that as Ls increases, the viscous damping effect at the wall is substantially reduced, resulting in an overall elevation of the mean velocity profile. Within the viscous sublayer, the mean velocity increment exhibits a linear relationship with Ls, satisfying the relation . In the near-wall region, the turbulence fluctuation intensity demonstrates an enhanced dependence on Ls, with the intensification of Q2 (ejection) events leading to an elevated Reynolds stress peak that shifts closer to the wall. Analysis of wall-attached low-speed streaks indicates that, for a dimensionless wall-normal structure scale
, both their number and volume increase significantly with rising Ls. Furthermore, it is found that the effects of the wall slip condition are confined to the near-wall region, while the outer inertia-dominated region continues to follow the scaling laws of no-slip wall turbulence.
Optimization of the natural frequencies of exponential functionally graded plates is a critical issue in the engineering field, playing a vital role in enhancing the dynamic performance of plate structures. In response to this challenge, this paper presents an innovative optimization approach that synergistically integrates the smoothed finite element method (SFEM) with surrogate models, aiming to address the problem with both high efficiency and precision. Based on the first order shear deformation theory, a SFEM is established for free vibration analyses of functionally graded plates. During the computation of the system stiffness matrix, gradient smoothing operations are applied to bending strains within smoothed domains, effectively improving computational accuracy. To overcome the shear locking phenomenon, different interpolation forms are adopted to treat bending strains and shear strains separately. For natural frequency optimization, a series of sample points are selected, and their corresponding natural frequencies are calculated using the SFEM. Subsequently, a surrogate model is established to map the relationship between the gradient index and the natural frequencies. The golden section method is employed to determine the optimal gradient index that achieves preset natural frequency targets. Numerical examples demonstrate that the surrogate model based on piecewise cubic Hermite interpolation exhibits high computational accuracy. Moreover, the surrogate model-based optimization significantly reduces the number of SFEM frequency calculations required, substantially enhancing optimization efficiency. This approach provides an efficient and practical method for optimizing natural frequencies of exponential functionally graded plates.
For the lightweight design problem of composite grid-stiffened sandwich cylindrical shells, this paper proposes an optimization strategy that couples the differential quadrature method (DQM) with intelligent algorithms. First, the buckling governing equations are established based on the smeared stiffener method and energy principle, which are efficiently solved using DQM and validated through finite element analysis. Subsequently, an artificial neural network (ANN) surrogate model and genetic algorithm (GA) are employed to achieve structural optimization. The results demonstrate that the critical buckling loads obtained via DQM agree well with finite element results, confirming the accuracy of DQM in analyzing the buckling of sandwich cylindrical shells. The ANN-based surrogate model exhibits high reliability in predicting the critical buckling loads of grid-stiffened sandwich shells. Moreover, the genetic algorithm, combined with theoretical results, efficiently yields lightweight design parameters. Case studies show that the optimized structure not only achieves significant weight reduction but also exhibits a substantial increase in critical buckling load.
This paper investigates the nonlinear vibration isolation problem of simply supported beam bridge structures under displacement excitation, employing the Incremental Harmonic Balance Method (IHBM) to derive an approximate analytical solution for the nonlinear vibration response of the beam. The research focuses on a coupled system consisting of a quasi-zero stiffness (QZS) isolator, constructed using a three-spring system, and a simply supported beam bridge. Based on the classical Euler-Bernoulli beam theory, the governing equations of motion under displacement excitation are established and systematically solved using the IHBM, with the entire analytical process thoroughly derived. The study transforms the final solution into a linear matrix equation using generalized coordinates, achieving a procedural and standardized computational process. To validate the reliability of the IHBM approximate analytical results, the study compares the IHBM computational results with numerical solutions obtained using the fourth-order Runge-Kutta method (ODE45). The results demonstrate that the IHBM method exhibits significant advantages in computational stability and result completeness. Additionally, through parametric analysis, the study explores the influence of key isolator parameters on system amplitude, further confirming the effectiveness and engineering practicality of the IHBM in nonlinear vibration isolation research. The research outcomes provide new theoretical foundations and methodological references for the nonlinear vibration isolation analysis of simply supported beam structures, offering important guidance for engineering practice.
As a critical step in flexible electronics packaging, the ultra-thin chip peeling process plays a vital role in ensuring high-yield manufacturing. This study focused on mechanical behavior differences between two peeling methods: roller-stretching and needle-ejecting. A theoretical model of the "chip-adhesive-substrate" laminated structure was established and validated by finite element simulation. A dual-criteria safety criterion was proposed to quantify process safety based on the competing relationship between interfacial fracture energy of the adhesive layer and surface cracking stress of the chip layer, which overcame the limitations of traditional methods for quantitatively evaluating the safety of the peeling process. Results demonstrated that the needle-ejecting procedure outperforms roller-stretching in terms of the safety of ultra-thin chip peeling. The roller-stretching process only has high engineering application prospect for peeling large-sized and thick chips from soft and thick substrates. Furthermore, an innovative stretching-ejecting combination technology is proposed, introducing the concept of synergistic matching to achieve chip stress neutralization and fracture mode optimization. This research provides theoretical insights into non-destructive ultra-thin chip peeling technology, and delivers practical guidance for advancing high-yield flexible microelectronics packaging.