Author Center
Editorial
Member of CSTAM
Dynamics>>More
• Featured Article
Chinese Journal of Theoretical and Applied Mechanics was founded in 1957, is a comprehensive bimonthly academic periodical, administrated by the Chinese Academy of Sciences(CAS), sponsored by the Chinese Society of Theoretical and Applied Mechanics and undertaken by the Institute of Mechanics, CAS. The current Editor-in-Chief is academician Wei Yueguang. Previous editor-in-chiefs in chronological order: Prof. H. S. Tsien, Prof. Y. H. Kuo, Prof. Zheng Zhemin, Prof. Lin Tongji, Prof. Huang Kezhi, Prof. Wu Chengkang, Prof. Wang Ziqiang, Prof. Yang Wei, Prof. Zheng Quanshui....... More>>
 Author:Chen Yunmin,Ma Pengcheng,Tang Yao Abstract: Numerical simulation and physical simulation are two main methods to analyze the settlement and stability of soil mass. As the mathematical equations of the soil stress-strain relationships, constitutive models are the foundations of numerical simulation. Soil is a type of granular materials, leading to three essential character. . . More>>
• 18 September 2020, Volume 52 Issue 5
Research Review
STATE OF THE ART OF CONSTITUTIVE RELATIONS OF HYPERELASTIC MATERIALS1)
Peng Xiangfeng, Li Luxian
2020, 52(5):  1221-1234.  doi:10.6052/0459-1879-20-189
Abstract ( 14 )   PDF (305KB) ( 8 )
Hyperelastic materials are commonly used in practical engineering with the prominent feature that a very large deformation may be produced under a force but the initial state can be completely recovered when the force is removed. Hyperelastic materials are typically nonlinear elastic ones, whose behaviors are in general characterized by their strain energy functions. For several decades, a lot of mathematical models and physical models have been proposed to study their constitutive relations through constructing the form of energy functions. However, a complete constitutive relation suitable for varied deformation modes and the entire deformation range is still the significant issue to expect in this field. This paper summarizes and analyzes the latest research status of constitutive relations of hyperelastic materials from three perspectives: (1) volume change modes including incompressible and compressible ones; (2) deformation modes such as uniaxial tension, shearing, biaxial tension and combined stretch and shear; (3) the entire range of deformation including small deformation, moderate deformation and large deformation. The latest progresses indicate that, in order to comprehensively describe experimental data of a given hyperelastic material and to apply it in practical problems, it is necessary to establish a complete constitutive relationship of compressible hyperelastic materials, which is suitable for varied deformation modes and the entire range of deformation. The authors suggest an implementation procedure for establishing the complete constitutive relationship of an actual hyperelastic material and an approach to construct the strain energy function of a compressible material.
Theme Articles on Thermal Stresses
Hu Keqiang, Gao Cunfa, Zhong Zheng, Chen Zengtao
2020, 52(5):  1235-1244.  doi:10.6052/0459-1879-20-127
Abstract ( 44 )   PDF (407KB) ( 28 )
FRICTIONAL TEMPERATURE ANALYSIS OF TWO-DIMENSIONAL ELASTO-PLASTIC WHEEL-RAIL SLIDING CONTACT WITH TEMPERATURE-DEPENDENT MATERIAL PROPERTIES1)
Fu Peilin, Ding Li, Zhao Jizhong, Zhang Xu, Kan Qianhua, Wang Ping
2020, 52(5):  1245-1254.  doi:10.6052/0459-1879-20-122
Abstract ( 67 )   PDF (569KB) ( 20 )
The accurate prediction of frictional temperature during wheel-rail sliding contact is of great significance to the studies for wear and fatigue performance of wheel-rail system. Current analytical or semi-analytical models of wheel-rail frictional temperature usually employ the elliptical distribution of contact pressure in the Hertz elastic contact theory and a single temperature-dependent material property, which differs from the actual heat-transfer state in many wheel-rail sliding contact conditions. Therefore, introducing the plastic correction of contact pressure and the temperature dependence of various properties in the calculation model of wheel-rail sliding temperature rise simultaneously could be great helpful to improve the accuracy of prediction result. Based on the elasto-plastic contact theory, considering the temperature dependence of thermal conductivity, specific heat capacity and friction coefficient simultaneously, the integration of thermal conductivity with respect to temperature is set as the quantity that needs to be solved by applying the Kirchhoff transformation method, and the nonlinear Fourier heat conduction equation is transformed into a corresponding simple partial differential equation with single variable coefficient, a unified implicit difference scheme with arbitrary form of temperature dependence is derived, and the influences of convection coefficient, vertical load, creepage and train speed on the temperature rise over rail surface are discussed, respectively. Results show that the convection coefficient has little effect on the temperature rise in the high-speed condition; the increase of creepage and train speed can enhance the friction power, and thus causes the increase of frictional temperature; and the maximum temperature rise also increases approximately linearly with the increasing vertical load. In addition, considering the temperature dependence of various thermophysical properties in the calculation model of temperature rise induced by the wheel-rail sliding contact can effectively avoid the overestimation of temperature rise, and the temperature-dependent friction coefficient has a more remarkable effect on the prediction result than thermal conductivity and specific heat capacity.
A GENERALIZED THERMOELASTIC DIFFUSION PROBLEM OF THIN PLATE HEATED BY THE ULTRASHORT LASER PULSES1)
Li Yan, He Tianhu, Tian Xiaogeng
2020, 52(5):  1255-1266.  doi:10.6052/0459-1879-20-118
Abstract ( 3 )   PDF (944KB) ( 2 )
In recent years, ultrashort laser pulses are widely used in the fields of ultra-precision machining, optical storage and microelectronic manufacture due to the advantages of high power density, short duration and high machining accuracy. In the manuscript, the memory-dependent nonlocal generalized thermoelastic diffusion theory is established based on the L-S generalized thermoelastic diffusion theory as well as considering the memory-dependent effect and spatial nonlocal effect. The theory can accurately predict the thermoelastic diffusion responses of structures whose geometry size is equivalent to its internal characteristic scale. The control equations of the theory are derived, and the solution of the control equations are obtained based on the Laplace integral transformation. As a numerical example, the transient thermoelastic diffusion responses of a semi-infinite thin plate subjected to a non-Gaussian laser pulse and a chemical shock are studied. The variation of the temperature, chemical potential, displacement, stresses and concentration with different nonlocal parameters, thermal time delay factors and diffusion time delay factors are obtained. The results show that heat conduction has significant effect on mass transfer, while mass transfer has little effect on heat conduction; nonlocal parameter has significant influence on displacement and stress, but little effect on temperature, chemical potential and concentration. The establishment of this theory and the solution method are aimed at accurately predicting the transient responses of the heat and mass under the impact of mechanical loading, heat and chemical potential.
A GENERALIZED PIEZOELECTRIC-THERMOELASTIC THEORY WITH STRAIN RATE AND ITS APPLICATION1)
Li Jiwei, He Tianhu
2020, 52(5):  1267-1276.  doi:10.6052/0459-1879-20-120
Abstract ( 3 )   PDF (1163KB) ( 1 )
The deformation of a large number of materials in engineering is between elasticity and viscosity, which exhibits both features of elastic solid and viscous fluid, that is, viscoelasticity. Viscoelasticity causes many mechanical relaxation phenomena, such as strain relaxation and hysteresis loss. In the investigation of transient response of multi-field problems subjected to thermal loading, it is especially important to take the phenomena such as thermal relaxation and the strain relaxation into consideration to accurately describe their transient response. For the transient response of the generalized piezoelectric-thermoelastic problems, although the generalized piezoelectric-thermoelastic model was established by taking into account the thermal relaxation, no strain relaxation is included so far. In present paper, by considering the strain relaxation in the process of deformation, a generalized piezoelectric-thermoelastic theory is theoretically established by extending Chandrasekharaiah's theory through taking strain rate into consideration. By means of the thermodynamic laws, the theory is formulated and the corresponding state equations and governing equations are obtained. In constitutive equation, a term of the product of a strain relaxation time and the strain rate is introduced, meanwhile, thermal relaxation time factors are included in constitutive equation and energy equation respectively. Subsequently, this theory is applied to investigating the dynamic response of a one-dimensional piezoelectric-thermoelastic problem subjected to a moving heat source. The Laplace transform and its numerical inverse transform are used to solve the problem, and the transient response under different strain relaxation time and heat source moving speed is obtained, that is, the distribution law of dimensionless temperature, displacement, stress and electric potential. The effect of strain rate on each physical quantity was investigated, and the results were presented in graphical form. The results show that the strain rate has a significant effect on the distributions of temperature, displacement, stress and electric potential.
THERMOELASTIC WAVE PROPAGATION IN PLATES: AN IMPROVED LEGENDRE POLYNOMIAL APPROACH1)
Wang Xianhui, Li Fanglin, Liu Yujian, Chen Huitao, Yu Jiangong
2020, 52(5):  1277-1285.  doi:10.6052/0459-1879-20-124
Abstract ( 61 )   PDF (1075KB) ( 25 )
In recent years, the research of thermoelastic coupled wave has greatly promoted the development of high temperature online detection and laser ultrasonic technology. For its small attenuation, long propagation distance and wide signal coverage, ultrasonic guided wave has become one of the rapid development directions in the field of nondestructive testing. However, the development of guided wave high temperature on-line detection and laser ultrasonic guided wave technology is slow. The key lies in the difficulty in solving the coupled thermoelastic wave equation and the difficulty in studying the propagation and attenuation characteristics. As an effective method, Legendre polynomial approach has been widely used to solve the problem of guided wave propagation since 1999. But there are two shortcomings in this method, which limit its further development and application. Two defects are: (1) Due to the Legendre polynomial and its derivative in integral kernel function, the integrals in the solution process leads to low calculation efficiency; (2) Only the thermoelastic guided wave propagation with isothermal boundary conditions can be treated. In order to solve these two defects, an improved Legendre polynomial method is proposed to solve the fractional thermoelastic guided waves in plates. The analytical integral instead of numerical integration in the available conventional Legendre polynomial approach, which greatly improves the calculation efficiency. A new treatment of the adiabatic boundary condition for the Legendre polynomial is developed by introducing the temperature gradient expansion based on the rectangular window function. Compared with the available data shows the validity of the improved method. Comparison with the CPU time between two approaches indicates the higher efficiency of the presented approach. Finally, the phase velocity dispersion curves, attenuation curves, the stress, displacement and temperature distributions for a plate with different fractional orders are analysed. The fractional order has weak influence on the elastic mode velocity, but it has considerable influence on the elastic mode attenuation.
EXACT THERMO-ELECTRO-ELASTIC SOLUTION OF FUNCTIONALLY GRADED MULTILAYERED ONE-DIMENSIONAL QUASICRYSTAL CYLINDRICAL SHELLS1)
Li Yang, Qin Qinghua, Zhang Liangliang, Gao Yang
2020, 52(5):  1286-1294.  doi:10.6052/0459-1879-20-126
Abstract ( 63 )   PDF (399KB) ( 113 )
Layered structures made of two and more materials with different properties can meet the needs of industrial development. However, the abrupt change of material properties at the interface of the laminated structures can easily cause some interface problems, such as stress concentration, interface cracks, and interface delamination phenomena. Functionally graded materials refer to utilize a continuously changing component gradient instead of the original sudden change interface, which can eliminate or weaken the abrupt change of the physical properties and then increase the bonding strength of the layered structures. In this paper, the research object is the functionally graded multilayered one-dimensional quasicrystal cylindrical shells. By virtue of the pseudo-Stroh formalism and the propagator matrix method, we establish the layered one-dimensional quasicrystal cylindrical shells model with the material parameters following the power-law type distribution along its radius direction, and obtain the exact thermo-electro-elastic solution of the functionally graded layered one-dimensional quasicrystal cylindrical shells with simply supported boundary condition. Numerical examples are carried out to investigate the influences of the exponential factor on temperature, electric, phason and phonon fields of the functionally graded layered one-dimensional quasicrystal cylindrical shells subjected to both inner and outer surfaces temperature variations, especially the effects on physical quantities at the inner and outer surfaces of the layered one-dimensional quasicrystal cylindrical shells. The obtained results indicate that: the exponential factor can change the distribution characteristic of material parameters, which can cause a significant influence on the physical quantities in the temperature, electric, phason, and phonon fields; by increasing the exponential factor, the deformation at the internal surface induced by temperature stimuli is reduced and the strength of the layered one-dimensional quasicrystal cylindrical shells is improved. The results obtained in this paper can provide a reliable theoretical basis for the design and manufacture of functionally graded layered one-dimensional quasicrystal cylindrical shells.
ANALYSIS OF THE PIEZOTRONIC EFFECT OF A PIEZOELETRIC SEMICONDUCTOR FIBER UNDER MUTIPLE LOCAL TEMPERATURE LOADINGS1)
Cheng Ruoran, Zhang Chunli
2020, 52(5):  1295-1303.  doi:10.6052/0459-1879-20-128
Abstract ( 39 )   PDF (1724KB) ( 14 )
The polarization potential resulted from a temperature change through a series coupling effects can significantly change physical fields in piezoelectric semiconductor (PS) structures, which has been found important engineering applications in wearable electronics and temperature-related semiconductor electron devices. Using the one-dimensional multi-field coupling model for thermo-piezoelectric semiconductors, we study the effect of multiple local temperature changes on the behaviors of PS fibers. Based on the linearized current constitutive relations, we obtain the analytical solutions. For instance, a PS fiber under two local temperature changes is numerically studied. The effect of local temperature changes on the distributions of the displacement, potential, electric displacement, polarization and carrier are examined. For large temperature changes, we conduct a nonlinear analysis by COMSOL using the nonlinear current model. The numerical results show that the potential barrier and well produced by the temperature changes through a series of coupling effects depend on the magnitude of the temperature change and the thermal load position. This provides a useful theoretical guidance in the designing of devices.
Fluid Mechanics
DEEP-LEARNING-BASED INTELLIGENT FORCE MEASUREMENT SYSTEM USING IN A SHOCK TUNNEL1)
Wang Yunpeng, Yang Ruixin, Nie Shaojun, Jiang Zonglin
2020, 52(5):  1304-1313.  doi:10.6052/0459-1879-20-190
Abstract ( 28 )   PDF (5705KB) ( 8 )
Aerodynamic force measurement in high-enthalpy flow is very important for the design and optimization of hypersonic vehicles. Currently, impulse facilities are used for generating high-temperature and high-pressure driving gas to simulate the high-enthalpy flow with hypersonic flight-conditions, such as a shock tunnel. However, when force tests are conducted in an impulse facility, the inertial force has a large influence on the measuring results, which creates low-frequency vibrations of the test model and its motion cannot be addressed through digital filtering. In the case of a few milliseconds of test time, the structural design of the six-component balance is greatly challenged. Therefore, dynamic calibration becomes very important for improving the precision and accuracy of force measurement during short-duration. A new method, deep-learning-based single-vector dynamic self-calibration of the force measurement system, and intelligent force measurement system are proposed for obtaining high-accurate aerodynamic force in impulse facilities. One of the main features of this dynamic calibration method is the calibration of the overall force measurement system, not just the balance. Applying this method, the calibrated force measurement system is the wind tunnel test object, which ensures the consistency of calibration and application. In the evaluation, the test verification has achieved relatively ideal results, the large-scale low-frequency vibration interference has been basically eliminated, and the accuracy and reliability of the force measurement in impulse facility have been greatly improved.
A RESIDUAL-BASED UNRESOLVED-SCALE FINITE ELEMENT MODELLING FOR IMPLICT LARGE EDDY SIMULATION1)
Chen Linfeng
2020, 52(5):  1314-1322.  doi:10.6052/0459-1879-20-055
Abstract ( 55 )   PDF (6002KB) ( 27 )
In consistence with large and small scales in turbulent flows, shape function space can be divided into resolved and unresolved scale spaces in a frame of finite element method. Introducing the same decomposition of the weighting function space, the variational formulations of Navier-Stokes equations can be divided into two systems of equations: resolved- and unresolved-scale equations. Generally, only the resolved-scale equation is computed, and the unresolved scales are modeled. Based on the unresolved-scale equations, an approximate residual-based unresolved-scale modeling is proposed in the present study. The large-scale equations are then computed by substituting the unresolved-scale modeling. The method is called residual-based large eddy simulation, in which unlike in the classical LES a filtering for Navier-Stokes equations is needed, multiscale decomposition is instead used. Numerical simulations of a turbulent channel flow are implemented with in-house codes of the residual-based large eddy simulation. The results show that, with a low number of elements, the mean streamwise velocity obtained using the present method is in agreement with the DNS data in the inner layer, and it is slightly overpredicted in the outer layer. Underprediction of the Reynolds stress by the present method causes a reduction of turbulence intensity transportation from the streamwise direction to the normal direction. Isosurfaces of the streamwise velocity reveals its capability of capturing the large-eddy structures. Meanwhile, low-speed streaks can be clearly observed in the sublayer near the wall.
MODELLING AND CALCULATION OF THE TURBULENT BOUNDARY LAYER ON A ROTATING CYLINDER SURFACE WITH STRONG SUCTION1)
Ievgen Mochalin, Lin Jingwen, Cai Jiancheng, Volodymyr Brazhenko, E Shiju
2020, 52(5):  1323-1333.  doi:10.6052/0459-1879-20-032
Abstract ( 33 )   PDF (1732KB) ( 6 )
A simple, fast and adequate approach to the turbulent boundary layer calculation on the surface of a rotating permeable cylinder has been elaborated for the case of a strong suction through the cylinder surface. Firstly, the rotational gap flow between two concentric cylinders was analyzed theoretically; the outer cylinder is stationary, and the inner porous cylinder is rotating with suction condition. Based on the fact that the stationary outer cylinder does not influence the flow near the rotating inner one, it can be treated as a boundary layer on the surface of rotating permeable cylinder, and an analytical expression for the circumferential velocity distribution is obtained. Secondly, the Cebeci-Smith two-layer algebraic turbulence model has been adjusted to account for centrifugal force field (streamlines curvature), wall suction, and low-Reynolds-number effects. Analytical corrections and empirical coefficients are used to tune the model for the specific conditions of coupled influence of the factors mentioned above. The calibration database was used which has been obtained by detailed numerical simulation based on the Reynolds stress turbulent model. The numerical simulation approach has been comprehensively verified in the known study for the specific flow conditions under consideration. Finally, the solution algorithm based on generalized Cebeci-Smith two-layer algebraic turbulence model was offered to solve the boundary layer flow over the rotating porous cylinder surface. The algorithm is suited for the situation of flow uniformity in the azimuthal and axial directions that required a special iterative procedure to be elaborated. The results of the algebraic turbulence model with different combinations of the rotational speed and the suction velocity agree well with the simulation results of the Reynolds stress turbulent model. It is demonstrated that the method developed reproduces also the laminar boundary layer at the same initial conditions when the detailed numerical simulation predicts the stable laminar flow in the inner cylinder boundary layer.
THE INFLUENCE OF GLOBAL-DIRECTION STENCIL ON GRADIENT AND HIGH-ORDER DERIVATIVES RECONSTRUCTION OF UNSTRUCTURED FINITE VOLUME METHODS1)
Kong Lingfa, Dong Yidao, Liu Wei
2020, 52(5):  1334-1349.  doi:10.6052/0459-1879-20-093
Abstract ( 10 )   PDF (6125KB) ( 2 )
The accuracy of unstructured finite volume methods is greatly influenced by different stencils. In previous work, based on the existing problems of local-direction stencil, we explored a more concise global-direction stencil selection method for the second-order unstructured finite volume solver, and stencil cells selected by this novel stencil selection method are always along the boundary normal and circumferential directions even on grids with high aspect ratio. As a result, the variation of flowfield is effectively captured, and flow anisotropy are well reflected. In addition, the novel method is topology-independent, since global directions are determined by the flowfield, while the local directions are strongly coupled with the grid. Therefore, the complex process of advancing front as well as local directions estimation are completely avoided in the novel stencil selection method, and the phenomenon that stencil cells deviate from the boundary normal vector is effectively eliminated on high-aspect-ratio triangular grids. What's more, a better computational accuracy and lower truncation errors on the second-order accurate finite volume solver are obtained by the employment of global-direction stencil. In order to further test the effectiveness of global-direction stencil on high-order unstructured finite volume methods, we will preliminarily utilize this stencil to test the effect of gradient and high-order derivatives reconstruction. After verification, computational errors of global-direction stencil are lower than that of local-direction stencil, and also lower than that of commonly used vertex-neighbor stencil on different grid types. Besides, errors of variable and derivatives at the Gauss point obtained by global-direction stencil are also the lowest among three methods we tested. Therefore, the global-direction stencil is well performed on gradient as well as high-order derivatives reconstruction, and it is feasible to extend this novel stencil selection method to high-order unstructured finite volume methods.
PIV MEASUREMENT OF CLOSE IMPINGING JET ON FLAT PLATE1)
Hu Jianjun, Zhu Qing, Wang Meid, Jin Yaolan, Wang Simin, Kong Xiangdong
2020, 52(5):  1350-1361.  doi:10.6052/0459-1879-20-142
Abstract ( 20 )   PDF (15552KB) ( 7 )
The time-resolved particle image velocimetry (TR-PIV) method was used to measure the flow field of a jet impinging on a flat plate in two orthogonal plane at very close distance. The effects of impinging distance and Reynolds number on the flow characteristics and the vortex topology in the clearance were analyzed. The results show that there are three kinds of typical vortex structure in the clearance, namely, double vortex ring mode, single vortex ring mode and complete entrainment mode. However, under the condition of turbulent state with large flow rate, the jet may break through the vortex ring and form random high-speed outflow. The occurrence of each flow pattern is mainly related to the strength of wall constraints. The energy transportation and flow loss of three typical flow patterns are investigated by vorticity analysis. The results show that the energy of the jet is transmitted outward through the vortex-ring mode at close impingement. In the double vortex ring mode, the two vortex rings have opposite rotational direction. Due to the constraint of the end face, both vortex rings are strictly constrained within the end face of the jet nozzle. The strength of the primary vortex ring is significantly greater than that of the secondary vortex ring. Finally, the proper orthogonal decomposition (POD) method is used to analyze the flow mode and energy distribution of the impinging jet. The first ten modal analysis of single and double vortex mode show that the energy fluctuation occurs in a paired pattern at a lower order, which indicates that both the primary and secondary vortex rings have good symmetry. Meanwhile, in the double vortex ring mode, the primary vortex ring is the dominant large-scale flow structure. The first three modes of the complete entrainment mode show that the energy of the jet is concentrated in the upstream of the jet and decreases sharply with the turbulent diffusion.
ARCJET ABLATION EXPERIMENT TO SIMULATE THE CHELYABINSK ASTEROID ENTRY1)
Luo Yue, Wang Lei, Dang Leining, Liu Jinbo, Zhang Jun, Liu Sen
2020, 52(5):  1362-1370.  doi:10.6052/0459-1879-20-081
Abstract ( 25 )   PDF (21285KB) ( 10 )
Ablation is one of the most important phenomenon when an asteroid enters the earth atmosphere at hypervelocity, which largely determines the mass loss, flight trajectory, and even radiation characteristics of the asteroid. To research the typical ablation process of asteroids when entering the earth atmosphere, experiments were conducted in an arc heater to simulate the typical conditions (velocity: 6 km/s, height: 17 km, diameter: 1 m) of Chelyabinsk asteroid event. The blunt-shaped test samples with the head radius of 20mm were made by carbon steel and basalt. In this work, the ablation process of test samples were clearly recorded, in which the melt flow of two different materials and the spallation of fragments as well as vaporization of basalt were observed. The evolutions of emission spectroscopy, recession profile and surface temperature profile during the whole process were acquired. The results indicate that the ablation phenomenon and the mechanism of mass loss of two materials are obviously different: Under the impact of the high-temperature flow, the carbon steel was sputtering into mass of tiny droplets which were washed away by the flow rapidly, while the mass loss of basalt were the shear flow of molten matter with small amount of massive spalling and evaporation spraying. All samples were exposed for four seconds in the plasma stream, mass loss and recession of the basalt and carbon steel were [37.9 g, 7.3 mm] and [72.7 g, 13.1 mm] respectively. The estimated effective enthalpy of ablation was 2.6 MJ/kg, the component measured by emission spectroscopy conforms to the electron microscopy (EDS) scanning.
RESEARCH ON IMPROVEMENTS OF LRN TURBULENCE MODEL BASED ON FLOW AROUND AUTOMOBILE BODY1)
Qiang Guanglin, Yang Yi, Chen Zhen, Gu Zhengqi, Zhang Yong
2020, 52(5):  1371-1382.  doi:10.6052/0459-1879-20-095
Abstract ( 64 )   PDF (9792KB) ( 23 )
The flow around automobiles was modularized into typical local flows in this paper. Through analyzing the characteristics of typical local flows, it is verified that the capture ability of turbulence model to transition is the key to accurately simulate the flow around automobiles. The paper optimized the solutions of steady-state and transient-state problems by analyzing the separation and transition mechanism of the flow, promoted the prediction ability of turbulence model for transition and improved the accuracy of turbulence model for automobile flow field simulation. For the steady-state solution of the flow around automobiles, by introducing the streamline curvature factor and the response threshold into the low Reynolds number (LRN) $k$-$\varepsilon$ model proposed by Jones and Lauder, a modified low Reynolds number turbulence model (modified LRN $k$-$\varepsilon$) which can predict transition more accurately was obtained. This model alleviated the problems of the original model's over-dependence on the turbulent dissipation rate and the insufficient prediction of the total stress development. For the transient-state solution, by analyzing the characteristics of the RNAS(Reynolds-averaged Navier-Stokes equations)/LES(large eddy simulation) mixed turbulence model, introducing the constrained large eddy simulation (CLES) method and the modified LRN $k$-$\varepsilon$ turbulence model proposed in this paper, a transition LRN CLES model that can accurately predict the transition was proposed. These improved models were applied to the simulation of the external flow field and buffeting noise of a real automobile model respectively. Computations were carried out using the ANSYS Fluent solver. The calculation results were compared with the simulation results of the commonly turbulence models, HD-2 wind tunnel test results and real vehicle road test results, it show that the improved turbulence models can more accurately simulate the steady-state and transient-state characteristics of the complex real automobiles, which provides a reliable theoretical basis and effective numerical solution method for the study of automotive aerodynamic.
Solid Mechanics
ANLYTICAL SOLUTION OF THERMOELASTIC DAMPING IN RECTANGULAR MINDLIN MICRO PLATES1)
Ma Hangkong, Zhou Chenyang, Li Shirong
2020, 52(5):  1383-1393.  doi:10.6052/0459-1879-20-175
Abstract ( 4 )   PDF (5474KB) ( 5 )
Analytical solution of thermoelastic damping for rectangular Mindlin micro plate with the four edges simply supported is presented for the first time. Based on the Mindlin plate theory and the one-way coupled heat conduct theory, governing differential equations for thermo-elastically coupled free vibration of the micro plate are formulated. Ignoring the in-plane variation of the temperature gradient, analytical solution of temperature field in terms of the kinematic parameters is obtained under the adiabatic boundary conditions at the top and bottom surfaces. Furthermore, the equations of structural vibration including the thermal bending moment are transformed into a fourth-order partial differential equation only in terms of the deflection. By using the mathematical similarity between the eigenvalue problems under the simply supported boundary conditions, analytical solution of the complex natural frequency of the Mindlin plate is expressed in terms of the frequency of isothermal Kirchhoff plate. Then the inverse quality factor which represents the level of the thermoelastic damping is obtained. Finally, the effects of the shear deformation, the material and the geometry parameters on the thermoelastic damping are examined in detail by the numerical results. The numerical results show that thermoelastic damping estimated by Mindlin plate theory is less than that by Kirchhoff plate theory. The difference between the values evaluated by the two plate theories becomes very significant near the critical thickness. Moreover, along with the increase of the side-to-thickness aspect ratio,the maximum of the thermoelastic damping in micro Mindlin plate increase monotonically, however, that of micro Kirchhoff plate keeps constant.
ANALYSIS OF 3-D NOTCHED/CRACKED STRUCTURES BY USING EXTENDED BOUNDARY ELEMENT METHOD1)
Li Cong, Niu Zhongrong, Hu Zongjun, Hu Bin
2020, 52(5):  1394-1408.  doi:10.6052/0459-1879-20-129
Abstract ( 38 )   PDF (13036KB) ( 53 )
According to the theory of linear elasticity, the conventional numerical methods are difficult to calculate the singular stress fields of three dimensional V-notched/cracked structures because of the stress singularity in the V-notch/crack tip region. In this paper, the extended boundary element method (XBEM) is first proposed to calculate the whole displacement and stress fields of three dimensional V-notch/crack structures. Firstly, the three dimensional V-notched/cracked structure is divided into two parts, which are a small sectoral column around the notch/crack tip and the outer region without the tip sectorial column. The displacement and stress components in the small sector column are expressed as the asymptotic series expansions with respect to the radial coordinate from the tip. The stress singular orders and the associated displacement and stress eigen-functions in the tip region are determined by the interpolating matrix method. The amplitude coefficients in the asymptotic series expansions are taken as the basic unknowns. Secondly, the boundary element method is used to analyze the three dimensional V-notched/cracked structure removed the small sector column. Hence, the whole displacement and stress fields of both the tip region and outer region are obtained by combining the boundary element analysis and the asymptotic series expansions of the displacement and stress fields in the notch/crack tip region, where the XBEM has the characteristics of the semi-analytic approach. The XBEM is suitable for the displacement and stress analysis of the three dimensional V-notched/cracked structures, and its solution can accurately describe the displacement and stress fields from the notch/crack tip to the whole region of the V-notched/cracked structures. Finally, two typical examples are given to demonstrate the effectiveness and accuracy of the extended boundary element method.
A MULTIAXIAL HIGH CYCLE FATIGUE LIFE PREDICTION MODEL CONSIDERING THE EFFECT OF MEAN STRESS FOR TENSION-TORSION LOADINGS WITH SAME FREQUENCY1)
Li Jing, Tong Xiaolong, Yang Shuo, Qiu Yuanying
2020, 52(5):  1409-1421.  doi:10.6052/0459-1879-20-066
Abstract ( 6 )   PDF (1022KB) ( 1 )
The influence of normal mean stress on fatigue life prediction has been well reflected for most of the critical plane-based high cycle fatigue life prediction models, whereas the effect of shear mean stress isn't well considered in these models. It is found that the fatigue life of 7075-T651 aluminum alloy is substantially reduced due to the existing of the shear mean stress by analyzing the experimental data of this aluminum alloy, which is similar to the effect of tensile mean stress. Therefore, nonconservative predictions maybe obtained under the loading paths with substantial shear mean stresses for these life prediction models ignoring the effect of mean shear stress. In order to estimate the fatigue life better, a new critical plane-based multiaxial high cycle fatigue life prediction model is proposed to take into account the effects of both normal and shear mean stresses. In the proposed model, the strain-based Fatemi-Socie criterion is first extended to the high cycle fatigue life prediction. And then a stress-based Fatemi-Socie criterion is developed. The shear and normal Walker factors are introduced in the developed criterion to consider the effects of shear and normal mean stresses, respectively. Both the shear and normal Walker factors vary from 0 to 1, which reflects the sensitivity of the material to shear and normal mean stresses. Procedures to determine the damage parameters acting on the critical plane and to calculate the constants contained in the proposed model are all presented. The proposed model is valid for the metallic materials with the ratio ${0.5<\tau _{-1} } / {\sigma _{-1} }<0.8$. Comparisons between test results of 5 kinds of metallic materials and model predictions under 12 types of loading paths with different mean stress levels showed that the proposed model presents relatively accurate predictions. Most of the predictions are fell within a life factor of $\pm$3.
FUNDAMENTAL FREQUENCY MAXIMIZATION DESIGN FOR CONTINUOUS FIBER-REINFORCED COMPOSITE STRUCTURES1)
Cheng Changzheng, Bian Guangyao, Wang Xuan, Long Kai, Li Jingchuang, Wu Qiaoguo
2020, 52(5):  1422-1430.  doi:10.6052/0459-1879-20-083
Abstract ( 88 )   PDF (3732KB) ( 109 )
Compared with traditional metal-materials, fiber-reinforced composite materials have better performance in many aspects such as strength, stiffness, and fracture resistance. At present, fiber-reinforced composite materials have been widely used in automotive, aerospace, and other industrial fields. This paper proposes a topology optimization method for solving the fundamental frequency maximization of undamped free vibration of continuous fiber-reinforced composite structures. To achieve the simultaneous optimization of the structural topological configuration and the fiber angle. A dynamic topological optimization model is established with the permitted material usage as the constraint and the structure's first-order eigenvalue as the objective function. The model includes density design variables that characterizes the topological configuration of the structure and angular design variables that characterizes the fiber orientation. The analytical sensitivity formulas of the objective function of eigenvalue with respect to density design variables and angle design variables are derived in detail, and the method of moving asymptotes (MMA) is used to solve the optimization problem. Finally, three numerical examples are performed to verify the effectiveness of the proposed method, which includes a static optimization example with the stiffness maximization as the goal and two dynamic optimization examples with the first-order eigenvalue maximization as the goal. The results show that the proposed method can achieve a stable iterative history and fast convergence, and can effectively improve the structural frequency while achieving the integrated optimization of the structural topological configuration and the fiber angle.
AN INTERVAL FINITE ELEMENT METHOD BASED ON THE NEUMANN SERIES EXPANSION1)
Wu Pengge, Ni Bingyu, Jiang Chao
2020, 52(5):  1431-1442.  doi:10.6052/0459-1879-20-152
Abstract ( 25 )   PDF (1126KB) ( 25 )
Uncertainty is common in the practical engineering. The interval finite element method is an interval method which introduces the numerical computational method of finite element to structural uncertainty analysis. The aim of the interval finite element analysis is to obtain the upper and lower response bounds of the structure with interval uncertain parameters, where solving the interval finite element equilibrium equations is the key issue. But the solution of interval linear equations belongs to a class of NP-hard problems which are often difficult to solve. This paper classifies and defines a type of linearly decomposable interval finite element problems, which exist commonly in practical engineering. To solve this type of problems, an interval finite element method based on Neumann series is proposed. It is named as the linearly decomposable interval finite element problem if the stiffness matrix in the interval finite element analysis formulation can be expressed as a linear superposition of a set of independent interval variables when the interval uncertain parameter is expressed as a linear superposition form of the independent interval variables. For this kind of problems, the inverse of the stiffness matrix can be represented by its Neumann series expansion. Thus the explicit expressions of structural responses with interval variables can be then obtained, with which the upper and lower bounds of the structural response can be solved efficiently. Finally, two numerical examples show the effectiveness and accuracy of the proposed method.
DYNAMIC MECHANICAL BEHAVIOR AND PENETRATION PERFORMANCE OF WFeNiMo HIGH-ENTROPY ALLOY1)
Chen Haihu, Zhang Xianfeng, Xiong Wei, Liu Chuang, Wei Haiyang, Wang Haiying, Dai Lanhong
2020, 52(5):  1443-1453.  doi:10.6052/0459-1879-20-166
Abstract ( 42 )   PDF (18906KB) ( 45 )
In order to investigate the deformation behavior and penetration performance of WFeNiMo high-entropy alloy under different strain rates, the static mechanical properties of the high-entropy alloy was tested by universal material testing machine and the dynamic mechanical properties of the high-entropy alloy was tested by the SHPB (split Hopkinson pressure bar). The micro mechanism of deformation characteristics of the alloy under different strain rates was also discussed. Based on the ballistic gun test platform, the fragments penetration performance of the high-entropy alloy and the typical tungsten alloy (93W-4.9Ni-2.1Fe, wt%) to the finite thickness steel target was studied. The relationship between the penetration process of the two kinds of alloy fragments and the target damage characteristics, the energy consumption of penetration and the impact velocity was analyzed. The results show that the yield strength and strain rate of the high-entropy alloy and tungsten alloy present a positively correlation. The yield strength of the high-entropy alloy is higher than the tungsten alloy under the same strain rate. With the increase of strain rate of deformation, the high-entropy alloy develops from the brittle fracture, quasi-cleavage with the mixing of tough and brittle characters to the fracture deformation mode with adhesive characteristics. The high-entropy alloy has a strong local adiabatic deformation ability and high shear sensitivity when the fragments penetrate into the thin steel targets. The energy consumption of the high-entropy alloy fragments penetrating into the target is lower than the tungsten alloy fragments under the same impact velocity. The high--entropy alloy has excellent mechanical properties and superior performance in the penetration ability. In addition to the traditional shear plug effect, there is a certain energy release characteristic when the thin target is impacted at high speed by the high-entropy alloy fragments and it has a good application prospect in the field of the preformed fragments.
Dynamics, Vibration and Control
EQUIVALENT DYNAMICS MODELING AND ANALYSIS OF LIQUID-FILLED SPACECRAFT WITH FUEL CONSUMPTION1)
Liu Feng, Yue Baozeng, Ma Bole, Shen Yunfeng
2020, 52(5):  1454-1464.  doi:10.6052/0459-1879-20-027
Abstract ( 4 )   PDF (1310KB) ( 3 )
The liquid in the tank of in-orbit spacecraft may experiences different types of motion, mainly including the overall rigid motion of liquid respect to the tank, the lateral sloshing of the liquid free surface, and the liquid rotation starting which will maybe eventually graduate into the rotary sloshing of the liquid free surface and/or the spinning motion of liquid, etc. The composite 3DOF-rigid-pendulum sloshing model is able to describe all of these main motion types of liquid, and it has been validated to be effective for the analysis of the liquid sloshing dynamics during the rotation-starting period. In this paper, the composite rigid-pendulum model is developed to investigate the large motion coupling dynamics of a variable-mass liquid-filled spacecraft with nonlinear liquid slosh and fuel consumption, by taking account of the variation of equivalent model parameters over liquid-fill ratio. Based on the Lagrangian equations, the orbital-attitude-slosh coupling dynamics model of a liquid-filled spacecraft is equivalently established by using the composite rigid-pendulum model to replace the liquid in a spherical tank. Then, simulations of a three-axis stabilized large angle attitude maneuver and a zero-impulse orbital maneuver of spacecraft are given for the coupling dynamics response analysis of the liquid-filled spacecraft system. Simulation results show that the liquid motion with respect to the tank will cause the position offset of spacecraft, and the orbital velocity of spacecraft will not converge to zero when the spacecraft executes zero-impulse maneuvers with fuel consumption; and that the liquid is likely to experience violent and complicated sloshing motion, which may cause the instability in the rigid-body motion of spacecraft, during the orbital maneuvering when the tank is eccentrically installed in the rigid hub.
DYNAMICS MODELING AND EXPERIMENT OF A FLEXIBLE CAPTURING MECHANISM IN A SPACE MANIPULATOR1)
Li Haiquan, Liang Jianxun, Wu Shuang, Liu Qian, Zhang Wenming
2020, 52(5):  1465-1474.  doi:10.6052/0459-1879-20-106
Abstract ( 8 )   PDF (1180KB) ( 3 )
An end-effector with flexible capturing mechanisms, which can accomplish on-orbit capturing operations with large tolerance, is a vital component of a large-scale space manipulator. Dynamics modeling and theoretical analysis of the flexible capturing mechanisms are very important for on-orbit servicing task simulation and prediction. In this paper, a dynamics model of a flexible capturing mechanism with three cables in a space end-effector is developed firstly. The three-dimensional absolute nodal coordinate formulation (ANCF) is used to create nonlinear finite elements of flexible cables. Both bending and longitudinal deformation of the cables are considered, furthermore, contact between the flexible cables and the rigid target is analyzed by introducing an intermediate cylinder reference coordinate system. Then, an experiment with passive spring suspension is built to validate the proposed model and signals of both motions and forces of the target are collected and compared with the simulation results. The comparison shows that the values of the simulation match well with the experimental measuring ones. The presented model could be used as supplements for the two-dimension planar air-bearing experiment and could be used for capturing task simulations of large-scale space manipulators. At last, capturing simulations of two representative on-orbit operations are conducted by co-simulation with the proposed model and a dynamics model of a large-scale space manipulator. One of the operations is the soft capturing process of an inspection task on the spacecraft surface, the other is the soft capturing of a floating target. The main difference between the two simulations is that the target in the first simulation is fixed on the base of the manipulator. Results of these simulations show that the soft capturing process can be accomplished on the prescribed condition.
FIXED POINT CHAOS AND FOLD/FOLD BURSTING OF A CLASS OF DUFFING SYSTEMS AND THE MECHANISM ANALYSIS1)
Chen Yani, Meng Wenjing, Qian Youhua
2020, 52(5):  1475-1484.  doi:10.6052/0459-1879-20-098
Abstract ( 30 )   PDF (20033KB) ( 36 )
In this paper, a class of bistable Duffing type system with two slow variables under new materials is explored. The system is simulated by time history diagram, phase diagram and bifurcation diagram, then the dynamic mechanism of the system under different parameters is analyzed theoretically. Firstly, this manuscript describes that when the amplitude parameter value is greater than 1, the system may exhibit fixed point chaos and explains the reason of fixed point chaos. Secondly, this manuscript introduces the phenomenon of Fold/Fold bursting in parameter space which is caused by the movement of the system from one side of the saddle-node surface to the other side. We also call it saddle-node bursting. In fact, when the system passes through the saddle-node surface, the number of equilibrium points changes. Then this manuscript uses the path of longitudinal parabolic to explain the mechanism of Fold/Fold bursting. And it is found that regardless of the value of constant coefficient term and amplitude, as long as a certain relationship is satisfied, there will always be Fold/Fold bursting. Next this manuscript uses the linear path to discuss the influence of newly added constant coefficient term. It is found that the position where the path intersects the saddle-node surface will affect the symmetry of the bursting, and the span of the path will affect the magnitude of the bursting oscillation. Finally, this manuscript uses the multiple inflection curve path to discuss the phenomenon when two incentive terms have specific relation. When $n=3$, the change of the constant coefficient term will make the system show Fold/Fold bursting with different times, and the maximum can reach triple bursting. Moreover, it is found that if you can find a path that can be divided into $n$ segments, and each segment will have an intersection with the saddle surface, then $n$ times Fold/Fold bursting will occur.
ANALYSIS OF FRICTION INDUCED VIBRATION UNDER THE ACTIVE DISTURBANCE REJECTION CONTROL FRAMEWORK1)
Piao Minnan, Wang Ying, Zhou Yajing, Sun Mingwei, Zhang Xinhua, Chen Zengqiang
2020, 52(5):  1485-1497.  doi:10.6052/0459-1879-20-149
Abstract ( 22 )   PDF (419KB) ( 35 )
Active disturbance rejection control(ADRC) is a practical control method with a two-degree-of-freedom structure. Due to its capability of handling multifarious disturbances in a straightforward and effective manner, ADRC has been successfully applied to many mechanical systems. However, the limit cycle vibration may be induced when employing the ADRC for mechanical systems with friction. At present, there is no precise analysis work about the friction induced vibration under the ADRC framework. Therefore, this paper investigates this problem by using the analysis tools of nonlinear dynamic systems. First, two representative friction models, static switch model and dynamic LuGre model, respectively, are considered, and active disturbance rejection controllers of different orders are designed for a class of second-order motion systems. Equivalent forms of the controllers are obtained and their relationships with the proportional-integral-derivative(PID) controller are revealed. Then, the limit cycle is calculated by using the shooting method combined with the pseudo arc-length continuation approach. Based on the Floquet theory, the stability, occurrence and type of bifurcation of the limit cycle can be determined. In addition, the local stability of the equilibrium points is analyzed based on the Jacobian matrix and approximate numerical method. Finally, the effects of the model and parameter of friction, the order and parameters of the ADRC on the limit cycle are investigated by numerical calculations. As shown by the calculation results, the parameter $\beta$, which determines the negative slope of the Stribeck effect, has a significant effect. When $\beta>1$, closed-loop systems with these two friction models have the same characteristics. Cyclic fold bifurcation(CFB) of the limit cycle occurs and the set of equilibrium points is locally stable. However, characteristics of these two closed-loop systems are totally different when $\beta<1$. As for the ADRC order, it is found that the order does not affect the conclusions in terms of the existence and stability of the limit cycle, and the stability of the set of equilibrium points. Moreover, low-order ADRC has superior performance in tackling the conflict between the friction compensation and stability robustness. These results can provide some guidelines on the understanding of practical phenomena, selection of the ADRC order, and parameter tuning.
VIBRO-IMPACT DYNAMICS OF PIPE CONVEYING FLUID SUBJECTED TO RIGID CLEARANCE CONSTRAINT1)
Wang Yikun, Wang Lin, Ni Qiao, Yang Mo, Liu Dezheng, Qin Tao
2020, 52(5):  1498-1508.  doi:10.6052/0459-1879-20-137
Abstract ( 33 )   PDF (16691KB) ( 38 )
The vibro-impact dynamics due to loose constraints have become one of the key scientific problems in the dynamical system of pipes conveying fluid. The impact force modeled by smoothed nonlinear springs varies continuously with time and displacement of the vibrating pipe, which cannot exactly capture the non-smooth characteristics of the saltation of state vectors for the pipes before and after an impact. In this paper, a non-smooth mathematical model of simply supported pipes conveying pulsating fluid, subjected to a rigid constraint somewhere along the pipe length is established, with consideration of the effect of the values of clearance and coefficient of restitution of the constraint. Especially, the periodic and aperiodic oscillations are investigated under various pulsating frequencies of the internal fluid. The transition matrices of the displacement and velocity of all nodes along the pipe before and after impact were derived based on a Galerkin's approach. The nonlinear equations of motion are solved via a fourth-order Runge-Kutta method, by applying the impact boundary conditions. Results show that the pipe is capable of displaying interesting vibro-impact behaviors in the presence of the rigid clearance constraint with the variation of pulsating frequency of the flowing fluid. With a clearance close to the maximum displacement of the pipe without constraint, periodic vibro-impact behaviors are observed with multiple impacts. The vibration velocities before the pipe impacts on the edge of the rigid clearance constraint decrease to zero gradually with the displacement invariant, which is called a dynamical stick-slip motion, also known as a typical non-smooth phenomenon. By decreasing the value of coefficient of restitution, the responses of the pipe may change from periodic vibrations to chaotic ones. This work provides an attractive strategy for further understanding of the nonlinear impact dynamics of pipes conveying fluid subjected to rigid clearance constraint based on non-smooth theories.
Biomechanics, Engineering and Interdiscipliary Mechanics
ANALYTICAL STUDY ON GROUND VIBRATION INDUCED BY MOVING VEHICLE1)
Wang Li'an, Zhao Jianchang, Wang Zuowei
2020, 52(5):  1509-1518.  doi:10.6052/0459-1879-20-033
Abstract ( 3 )   PDF (14782KB) ( 1 )
Based on vehicle dynamics, the coupling vibration model between vehicle and ground is established, and the contact model of elastic roller is adopted to reflect the tyre tolerance. At the same time, the longitudinal and vertical forces between the wheel and the ground are considered, the system dynamic control equation is constructed, and the analytical solution of the surface vibration displacement is obtained by using Fourier and Laplace integral transformation. In the numerical example, the inverse discrete Fourier transform and Crump's method are used to do the numerical inversion, and the time domain solution of the surface vibration displacement is obtained, the influence of the parameters of the surface vibration displacement is analyzed. The results show that the surface irregularity has the most significant influence on the wheel-earth interaction, and the more uneven the ground, the greater the wheel-earth interaction and the greater the surface vibration displacement. The influence of vehicle speed on the wheel-ground force is limited, but it has a great influence on the excitation frequency of load. When vehicle speed increases, the excitation frequency increases, and the surface vibration displacement increases accordingly. At a low speed, tire inclusivity has a certain effect on the wheel-ground force and surface vibration. With the increase of tire inflation pressure, the wheel-ground force and surface vibration displacement increase, but with the increase of speed, this effect will gradually disappear.
ANISOTROPIC STRENGTH CRITERION BASED ON T CRITERION AND THE TRANSFORMATION STRESS METHOD1)
Wan Zheng, Meng Da
2020, 52(5):  1519-1537.  doi:10.6052/0459-1879-20-138}
Abstract ( 5 )   PDF (1817KB) ( 1 )
Geotechnical materials are usually distributed horizontally in layers, which can be regarded as transverse isotropic materials. Transverse isotropy has a significant effect on the deformation and strength of geotechnical materials. Based on the proposed t-strength criterion, the t-strength criterion is proposed on the physical concept of the existence of an effective slip plane in the isotropic element body, and the failure condition is fulfilled when the ratio of the principal shear stress to the principal normal stress on the space effective slip plane reaches a certain threshold. Effective slip plane in space and physical deposition surface, based on the position relationship between the two surfaces in the space, the idea is proposed that two surfaces angle as a characterization of transverse isotropy can influence the shear strength parameters, and assume that when the angle of value, the greater the strength of the anisotropy of the contribution degree is, the greater the corresponding stress ratio strength value is, on the contrary, the corresponding smaller stress ratio strength value is. Based on the above ideas, the two surfaces angle is regarded as a parameter to construct anisotropic stress ratio strength formula, and use the stress strength formula for correction, a new anisotropic criterion has been put forward based on t strength criterion, finally transverse isotropic t criterion formula of strength criterion is proposed. On the basis of the above criteria, considering transverse isotropic stress space into isotropic stress space, based on the rule of anisotropic t, transformation equation is derived based on the transverse isotropic strength t criterion. It can be used to convert the traditional p and q as the variables of the isotropic constitutive model into the transverse isotropic three-dimensional constitutive model. By predicting the strength of rock and soil materials and the stress-strain relationship test data under true triaxial condition, the validity and applicability of the proposed transverse isotropic t criterion and its stress transformation formula are verified.
STUDY ON DEBONDING FAILURE MODEL OF POLYUREA-BASED COATING WITH HIGH VELOCITY WATER FLOW BASED ON COHESIVE ZONE MODEL
Li Bingqi, Zhang Zhenyu, Li Bin, Liu Xiaonan, Yang Xuhui
2020, 52(5):  1538-1546.  doi:10.6052/0459-1879-20-064
Abstract ( 22 )   PDF (3002KB) ( 32 )
The surface of the flood discharge building is usually sprayed with polyurea-based coating to improve the impact resistance. However, the research on the debonding failure mechanism of the anti-wear polyurea-based coating under the action of high-speed water flow velocity is still blank. Based on the flow pattern of high-velocity water flow, the mechanical model of high-velocity water flow to flood-discharge building was determined, and the load on flood discharge building caused by high-velocity water flow mainly includes drag force, impact force, fluctuating force and lifting force. The cohesive zone model was used to characterize the debonding failure process of the interface between polyurea-based coating and flood-reducing building, and the debonding failure model of polyurea-coating with high-speed water flow was established, and the finite element formal equations, constitutive relationship, damage initiation principle, evolution principle and contact and collision model of the model are given. The relationship between stress-displacement in the process of debonding failure was obtained by the debonding failure tests, and the variation law between peeling failure stress and inclination angle of interface was obtained. According to the peeling failure test, the parameters of debonding failure model were obtained, and the model was also verified. The test results were in good agreement with the model calculation results, which provides a theoretical basis for the anti-shock and wear-resistant design of flood discharge building.
 NUMERICAL STUDY ON DYNAMIC STALL FLOW CONTROL FOR WIND TURBINE AIRFOIL USING PLASMA ACTUATOR null Abstract ( 8 ) PDF (0KB) (4) Accepted：2020-09-15
 ADVANCES IN KEY MECHANICAL PARAMETERS FOR RELIABILITY ASSESSMENT OF HIGH-SPEED TRAIN BEARINGS null Abstract ( 11 ) PDF (0KB) (3) Accepted：2020-09-15
 Research on Load Test Spectrum of EMU Car Bogies Based on Damage Consistency null Abstract ( 7 ) PDF (0KB) (3) Accepted：2020-09-15
 DYNAMICS ANALYSIS OF STOCHASTIC SPATIAL FLEXIBLE MULTIBODY SYSTEM Abstract ( 13 ) PDF (656KB) (8) Accepted：2020-09-15
 Stability and mode evolution of a cantilevered fluid-conveying pipe with lumped mass Hao-Ran HaoYi Hu-Liang HuDAI Lin Wang Q. Ni Abstract ( 4 ) PDF (1067KB) (4) Accepted：2020-09-15
•  2020 No.4 No.3 No.2 No.1 2019 No.6 No.5 No.4 No.3 No.2 No.1 2018 No.6 No.5 No.4 No.3 No.2 No.1 2017 No.6 No.5 No.4 No.3 No.2 No.1 2016 No.6 No.5 No.4 No.3 No.2 No.1 2015 No.6 No.5 No.4 No.3 No.2 No.1 2014 No.6 No.5 No.4 No.3 No.2 No.1 2013 No.6 No.5 No.4 No.3 No.2 No.1 2012 No.6 No.5 No.4 No.3 No.2 No.1 2011 No.6 No.5 No.4 No.3 No.2 No.1
•  ADVANCES IN MULTIPHASE SEEPAGE CHARACTERISTICS OF NATURAL GAS HYDRATE SEDIMENTS 1) Cai Jianchao, Xia Yuxuan, Xu Sai, Tian Haitao Chinese Journal of Theoretical and Applied Mechanics. 2020 Vol. 52 (1): 208-223 doi: 10.6052/0459-1879-19-362 Abstract (1647)   HTML (88)    PDF (899KB) (676)
 ADVANCES IN DYNAMIC MODELING AND OPTIMIZATION OF FLEXIBLE MULTIBODY SYSTEMS1) Sun Jialiang,Tian Qiang,Hu Haiyan Chinese Journal of Theoretical and Applied Mechanics. 2019 Vol. 51 (6): 1565-1586 doi: 10.6052/0459-1879-19-212 Abstract (509)   HTML (48)    PDF (7422KB) (543)
 PROGRESS IN CONSTITUTIVE BEHAVIOR OF ALUMINUM FOAM 1) Jia Ran, Zhao Guiping Chinese Journal of Theoretical and Applied Mechanics. 2020 Vol. 52 (3): 603-622 doi: 10.6052/0459-1879-20-020 Abstract (440)   HTML (49)    PDF (17934KB) (527)
 MECHNICAL PROPERTIES AND BEHAVIORS OF HIGH ENTROPY ALLOYS 1) Li Jianguo,Huang Ruirui,Zhang Qian,Li Xiaoyan Chinese Journal of Theoretical and Applied Mechanics. 2020 Vol. 52 (2): 333-359 doi: 10.6052/0459-1879-20-009 Abstract (347)   HTML (29)    PDF (85411KB) (265)
 DOES STRUCTURE DETERMINE PROPERTY IN AMORPHOUS SOLIDS? 1) Wang Yunjiang, Wei Dan, Han Dong, Yang Jie, Jiang Mingqiang, Dai Lanhong Chinese Journal of Theoretical and Applied Mechanics. 2020 Vol. 52 (2): 303-317 doi: 10.6052/0459-1879-19-368 Abstract (345)   HTML (29)    PDF (4835KB) (330)
 CONSTITUTIVE MODELS AND HYPERGRAVITY PHYSICAL SIMULATION OF SOILS$^{\bf 1)}$ Chen Yunmin,Ma Pengcheng,Tang Yao Chinese Journal of Theoretical and Applied Mechanics. 2020 Vol. 52 (4): 901-915 doi: 10.6052/0459-1879-20-059 Abstract (331)   HTML (28)    PDF (13172KB) (419)
 THEORY, MODELLING AND COMPUTATION OF NONLINEAR OCEAN INTERNAL WAVES1) Wang Zhan, Zhu Yuke Chinese Journal of Theoretical and Applied Mechanics. 2019 Vol. 51 (6): 1589-1604 doi: 10.6052/0459-1879-19-326 Abstract (315)   HTML (17)    PDF (23102KB) (280)
 International and national academic activities organized or sponsored by CSTAM in 2020 Chinese Journal of Theoretical and Applied Mechanics. 2020 Vol. 52 (1): 292-300 Abstract (282)      PDF (186KB) (196)
 THE BOUNCE DYNAMICS OF A RISING SINGLE BUBBLE NEAR A VERTICAL WALL 1) Sun Jiao, Zhou Wei, Cai Runze, Chen Wenyi Chinese Journal of Theoretical and Applied Mechanics. 2020 Vol. 52 (1): 1-11 doi: 10.6052/0459-1879-19-228 Abstract (275)   HTML (34)    PDF (17585KB) (222)
 CRACKING SIMULATION BASED ON A NONLOCAL MACRO-MESO-SCALE DAMAGE MODEL1) Lu Guangda, Chen Jianbing Chinese Journal of Theoretical and Applied Mechanics. 2020 Vol. 52 (3): 749-762 doi: 10.6052/0459-1879-19-319 Abstract (238)   HTML (4)    PDF (9560KB) (126)
 WIND-TUNNEL EXPERIMENT ON SAND DEPOSITION MECHANISM AND OPTIMAL MEASURES OF WIND-BREAK WALL ALONG RAILWAY IN STRONG WIND AREA1) Xin Guowei, Huang Ning, Zhang Jie Chinese Journal of Theoretical and Applied Mechanics. 2020 Vol. 52 (3): 635-644 doi: 10.6052/0459-1879-20-045 Abstract (232)   HTML (9)    PDF (14395KB) (81)
 DUALITY RELATIONS OF BEAMS IN NATURAL VIBRATIONS 1) Hu Haiyan Chinese Journal of Theoretical and Applied Mechanics. 2020 Vol. 52 (1): 139-149 doi: 10.6052/0459-1879-20-019 Abstract (229)   HTML (13)
 THE MECHANISM OF JETTING BEHAVIORS OF AN OSCILLATING BUBBLE1) Li Shuai, Zhang Aman, Han Rui Chinese Journal of Theoretical and Applied Mechanics. 2019 Vol. 51 (6): 1666-1681 doi: 10.6052/0459-1879-19-219 Abstract (227)   HTML (14)    PDF (38087KB) (139)
 RESERCH PROGRESS ON INTERACTION BETWEEN A SINGLE CAVITATION AND FREE SURFACE1) Guo Wenlu,Li Hongchen,Wang Jingzhu,Wang Yiwei,Huang Chenguang Chinese Journal of Theoretical and Applied Mechanics. 2019 Vol. 51 (6): 1682-1698 doi: 10.6052/0459-1879-19-328 Abstract (224)   HTML (14)    PDF (33133KB) (181)
 AN IMPROVED FIFTH-ORDER WENO-Z+ SCHEME 1) Luo Xin, Wu Songping Chinese Journal of Theoretical and Applied Mechanics. 2019 Vol. 51 (6): 1927-1939 doi: 10.6052/0459-1879-19-249 Abstract (224)   HTML (8)    PDF (3779KB) (105)
 DYNAMICS MODELING, SIMULATION, AND CONTROL OF ROBOTS WITH FLEXIBLE JOINTS AND FLEXIBLE LINKS$^{\bf 1)}$ Fang Wuyi,Guo Xian,Li Liang,Zhang Dingguo Chinese Journal of Theoretical and Applied Mechanics. 2020 Vol. 52 (4): 965-974 doi: 10.6052/0459-1879-20-067 Abstract (213)   HTML (12)    PDF (2010KB) (107)
 ACTIVE-CONTROL-ORIENTED DYNAMIC MODELLING FOR ON-ORBIT ASSEMBLY SPACE STRUCTURE1) Wang Enmei, Wu Shunan, Wu Zhigang Chinese Journal of Theoretical and Applied Mechanics. 2020 Vol. 52 (3): 805-816 doi: 10.6052/0459-1879-19-375 Abstract (212)   HTML (5)    PDF (2433KB) (111)
 EXPERIMENTAL INVESTIGATION ON THE EFFECT OF PARTICLES ON THE WAKE OF HEMISPHERIC DISTURBANCE1) Zhang Yantao, Sun Jiao, Gao Tianda, Fan Ying, Chen Wenyi Chinese Journal of Theoretical and Applied Mechanics. 2020 Vol. 52 (3): 728-739 doi: 10.6052/0459-1879-19-353 Abstract (210)   HTML (2)    PDF (20860KB) (72)
 RESEARCH PROGRESS IN HYDRODYNAMICS OF WIND-WAVE COMBINED POWER GENERATION SYSTEM1) Zhou Binzhen,Hu Jianjian,Xie Bin,Ding Boyin,Xia Yingkai,Zheng Xiaobo,Lin Zhiliang,Li Ye Chinese Journal of Theoretical and Applied Mechanics. 2019 Vol. 51 (6): 1641-1649 doi: 10.6052/0459-1879-19-202 Abstract (204)   HTML (8)    PDF (11972KB) (130)
 EXPERIMENTAL STUDY ON SEA ICE TENSILE STRENGTH BASED ON BRAZILIAN TESTS1) Chen Xiaodong, Cui Haixin, Wang Anliang, Ji Shunying Chinese Journal of Theoretical and Applied Mechanics. 2020 Vol. 52 (3): 625-634 doi: 10.6052/0459-1879-20-036 Abstract (201)   HTML (13)    PDF (7384KB) (106)
•  » Semi-analytical solution of seismic response of serially isolated equipment and supporter Yongfeng Du Yanhui Liu Hui Li Chinese Journal of Theoretical and Applied Mechani. 2009 Vol. 41 (3): 440-448 Cited by: Baidu(28)
 » Sub-grid scale combustion models for two-phase combustion flow in an annular combustor Yingwen Yan Jianxing Zhao Jingyong Zhang Yong Liu Chinese Journal of Theoretical and Applied Mechani. 2006 Vol. 38 (6): 776-784 Cited by: Baidu(11)
 » Convergence of integral functional and variational solutions of displacement and stress of soil due to pile-driving Gao Zikun Jianyong Shi Chinese Journal of Theoretical and Applied Mechani. 2009 Vol. 41 (4): 555-562 Cited by: Baidu(10)
 » Topology optimization of continuum structures with materials exhibiting different tensile and compressive properties Hongwu Zhang Chinese Journal of Theoretical and Applied Mechani. 2008 Vol. 40 (5): 646-653 Cited by: Baidu(9)
Former Editors-in-Chief
11th' Editorial Committee

Editor-in-Chief
Wei Yueguang
Editor-in-Chief
Wei Yueguang

Associate Editor-in-Chief
Chen Liqun　　Feng Xiating 　　　 Lu Xiyun
Wang Yuesheng　　　Wang Guangqian