Acta Mechanica Sinica最新文献

This study examines the hemodynamics of non-Newtonian blood flow in stenosed arteries, focusing on the roles of stenosis severity, guidewire presence, and various non-Newtonian constitutive models. Computational simulations using the generalized power-law, Casson, and Carreau-Yasuda models are conducted for stenosis severities of 50%, 70%, 80%, and 90%. Results indicate that stenosis severity exerts the greatest influence on pressure and wall shear stress (WSS), with increasing severity leading to higher pressure drops and WSS maxima. Guidewire presence reduces recirculation zone lengths by nearly 60% across different severities and raises the trans-stenotic pressure drop up to 120%. While the choice of constitutive model has minimal impact on hemodynamics within the stenotic region, it becomes crucial in healthy vessels, where non-Newtonian effects are more pronounced. In cases with a guidewire, pressure gradients in the healthy region show up to 18.8% differences between non-Newtonian models. These findings highlight the dominant roles of stenosis severity and guidewire presence in shaping hemodynamics within stenotic regions while emphasizing the need for precise constitutive modeling to capture flow characteristics in healthy vascular segments.

As a clean energy source, microwaves have significant potential to enhance the efficiency of breaking hard rock for mechanical excavation in tunnelling engineering. In this study, the variable angle shear tests were carried out on granite specimens after microwave irradiation. Two-dimensional digital image correlation (2D-DIC) and acoustic emission (AE) technologies were used to monitor the mechanical behavior change during loading. The results indicate as the microwave heating time increased from 0 to 5 min and then to 10 min, the cohesion of the granite specimen decreased from 32.603 to 29.355 MPa and then to 24.666 MPa. But there is only a slight increase in the internal friction angle. As the shear angle and microwave irradiation duration increase, the cumulative AE counts and AE energy significantly decrease, indicating that microwaves can effectively reduce the energy needed in hard rock breaking. The change in AE b-value during loading can be divided into five stages: a rapid rise period, a stable fluctuation period, a rapid downward period, a small rising period, and a final downward period. Notably, the small rising period of the b-value between two rapid downward periods can be used as a precursor to failure. However, microwave irradiation can shorten or even eliminate this period. As observed through 2D-DIC, the strain evolution and crack propagation process of granite specimens are significantly influenced by both the shear angle and microwave irradiation. Decreasing shear angle and increasing heating time leads to a more dispersed distribution of the maximum principal strain concentration zone during loading, and the final failure pattern of the sample becomes more complex. Through the quantitative calculation of the normal force on the hob, it is found that microwave pretreatment can effectively enhance excavation rates. The maximum penetration increased to 15.67 mm/r for the specimens after 10 min irradiation, representing a 46.18% improvement compared to untreated specimens.

During the past decades, the numerical methods based on Navier-Stokes (N-S) equations and direct simulation Monte Carlo (DSMC) methods have been proven effective in simulating flows in the continuum and rarefied regimes, respectively. However, as single-scale methods, they face challenges in addressing common multi-scale problems, which are essential to simulate hypersonic flows around near-space vehicles and the flows in the micro-electro-mechanical systems. Hence, there is an urgent need for a method to predict multi-scale flows. In this work, a quantified model-competition (QMC) mechanism for diatomic multi-scale flows is derived from the integral solution of the Rykov model equations. This mechanism encapsulates both continuum and rarefied behaviors in a cell, weighted according to its local physical scale. By building upon the QMC mechanism, the N-S solver and DSMC solver are directly integrated within a cell to devise a simplified unified wave-particle (SUWP) method for diatomic gases. Specifically, the two-temperature equations considering the rotational energy are introduced into the kinetic inviscid flux scheme and the N-S solver. As to the particle part, the collisionless DSMC solver is utilized to describe the non-equilibrium phenomenon. The proposed SUWP method for diatomic gases undergoes validation across a series of cases, including zero-dimensional homogeneous gas relaxation, one-dimensional normal shock structure, two-dimensional flows around a flat plate and cylinder, and three-dimensional flows past a sphere and blunt cone. Additionally, the implementation details of multi-scale wave-particle methods analysis and discussion are also undertaken in this work.

Various physiological diseases and injuries, such as stroke or brain tissue concussion, are closely associated with the rupture of blood vessels, so that blood pressure plays an important role in their progression. Predicting the pressure-induced rupture of soft tissue can be challenging, however, due to the mechanical nonlinearity, tension-compression asymmetry, and softness-induced large deformations. A phase-field model (PFM) in the finite deformation regime is thus herein proposed to study the blood pressure-induced rupture, considering the tissue’s tension-compression asymmetry and nonlinearity. A staggered scheme is proposed to solve the coupling of the deformation and the phase field problems, and implemented in commercial software. With a proposed regularization of the phase field inside the pressure-induced rupture region, the model’s accuracy is first verified in terms of its capability to predict the crack opening displacement, and compare well with direct numerical simulation performed under the assumption of no crack propagation. Pressure-induced rupture in soft brain tissue with multiple initial tears (seen as crack-like defects) is subsequently examined in detail. This proposed PFM further incorporates blood transport during the rupture process, showing potential in predicting the realistic behavior of soft tissue.

This study employs the direct simulation Monte Carlo method to investigate two-dimensional compressible decaying isotropic turbulence under high Mach number conditions, focusing on the effects of thermal non-equilibrium (TNE) and molecular thermal fluctuations. Simulations are performed for low-temperature cases involving rotational non-equilibrium, followed by high-temperature cases emphasizing vibrational non-equilibrium. The results demonstrate that the initial TNE state significantly impacts turbulence compressibility. Specifically, for initially rotationally hot cases, elevated translational temperatures strongly suppress turbulence compressibility, resulting in a slower decay of turbulent kinetic energy. These findings are also applicable to initially vibrationally hot cases, but the influence of TNE diminishes as the vibrational relaxation number Z_vib increases. Moreover, increasing Z_vib leads to a significant lag of vibrational temperature fluctuations relative to translational and rotational temperature fluctuations. Analysis of the turbulent energy and temperature spectra reveals that molecular thermal fluctuations dominate at length scales (i.e., crossover length scales) comparable to the turbulent dissipation length scale, causing the spectra to increase linearly with the wavenumber. For cases with initially rotationally or vibrationally hot conditions, the suppression of compressibility leads to a significant increase in the crossover length scale.

Metamaterials programmed with target rate-dependent mechanical properties are efficient platforms for realizing advanced functionalities. Yet, the loading rate-dependent mechanical property programming has received limited attention. Here, the “stair-building” strategy is employed in the rate domain by combining the bistability with viscoelasticity. An arbitrary target curve in the programmable space can be approximated by a “stair” built by two kinds of “bricks”. The “bricks” can be realized by a dual-bistable unit, constructed by two bistable structures in series. The dual-bistable unit can switch between two efficient stable phases without inducing changes in the global morphology. Such a unit exhibits N-shaped stress-strain curves at both efficient stable phases with different peak values, resulting in different heights of “bricks”. Moreover, the N-shaped curves have rate-dependent peak values, indicating that the heights of “bricks” change with loading rate. The “stair-building” strategy is realized by array-structured mechanical metamaterials based on dual-bistable units. Different stress-strain curves under various loading rates can be reprogrammed in the same piece of metamaterial by intentionally selecting the efficient stable phases of units. Besides, the rate effect of the metamaterial can also be tuned by reprogramming stress-strain curves under both low and high loading rates, respectively. This reprogrammable metamaterial is promising in smart vibration isolators and adaptive energy absorbers.

This study investigated the impact fragmentation behavior of projectiles with different yield strengths using 45 steel, 35CrMnSiA steel, and T12A steel as core materials. Through experimental and simulation analyses, the fracture mechanisms and damage characteristics of these materials under high strain rate impacts were examined. The results indicated that 45 steel and 35CrMnSiA steel exhibited noticeable ductility, with fracture surfaces showing prominent dimples and shear lips, experiencing significant plastic deformation and necking during impact. In contrast, T12A steel demonstrated brittle behavior, with fracture surfaces characterized by smooth, bright areas and fine radial streaks. In ballistic tests, T12A steel projectiles displayed distinct damage patterns compared to the other two materials due to their brittleness, resulting in shorter residual lengths. The force-time curves obtained from experiments and simulations showed that higher core strength shortened the time to reach peak force and reduced the overall contact duration. However, the high strength of T12A steel was not always associated with higher peak stresses, as its brittle fracture mode led to early instability. Fragment analysis revealed that smaller fragments (< 4 mm) primarily originated from the projectile head, while larger fragments came from the tail. Microscopic examination of the fragments revealed a mix of ductile and brittle fracture modes, with different particle sizes exhibiting distinct fracture surface features. Overall, the study provides insights into the impact performance and failure mechanisms of different core materials, highlighting the importance of material properties in determining projectile behavior under high strain rate conditions.

Dielectric elastomer (DE) is an electroactive polymer with the characteristics of high energy output, great flexibility, light-weight, mechanical compliance, and low cost, which are particularly suitable for DE energy generators. Energy harvesting efficiency is a key index to evaluate the performance of the energy generator, which depends on the structural configuration and the mechanical and dielectric properties of the DE material. This paper proposes a fractional viscoelastic polarization (FVP) model by combining the fractional viscoelasticity model and the polarization-based lumped parameter model. A dynamical model of a cone dielectric energy generator (CDEG) considering stretch-dependent electrostriction and nonlinear viscoelasticity is established. Additionally, a deep neural network (DNN) model is developed to explore the relationships between various parameters and the output energy of CDEGs to efficiently and accurately predict the energy output of CDEGs. Based on the DNN model, optimal parameter designs for CDEGs are obtained by using non-dominated sorting genetic algorithm II (NSGA-II). The experiments verified that the FVP model predicts accurately the output energy of CDEG and the established optimal design framework can accurately provide the optimal design parameters of CDEG, which offers deep insights for the design and fabrication of a high-efficiency dielectric energy generator.

The time spectral approach, a spectral method based on the Fourier series with an appropriate convergence speed, can be utilized for a time-varying problem like the flow around a pitching airfoil. This approach has the drawback of having a constant number of time intervals over the entire computational domain, which unnecessarily uses up more computer memory and central processing unit (CPU) time. By distributing time intervals in the computational domain optimally (proportional to the flow gradient), the adaptive time spectral approach can overcome the shortcoming of the time spectral method. In the current study, the adaptive time spectral method is added to an inviscid fluid flow solver. Also, in the airfoil with pitching motion, a grid known as an overset grid has been used, including two grids with an overlapping region. The results for the three cases (Cases 1, 2, and 5) of the NACA0012 pitching airfoil with different angles of attack studied by AGARD Institute, with Mach numbers 0.6, 0.6, and 0.755, respectively, showed that while having an acceptable solution accuracy, the amount of computer memory and CPU time is significantly reduced compared to the standard time spectral method.