Dedication on the Occasion of the Thirtieth Anniversary of the Founding of Institute of Applied Mathematics

In March of 1979, Chinese Academy of Sciences （CAS） established, with the approval of the State Council of China, an office for promoting the application of mathematics and Interdisciplinary studies in practice. Later in October of 1979, based on this office CAS established the Institute of Applied Mathematics （IAM）. The first director of IAM was the world-wide famous mathematician, Professor HUA Loo-Keng, and most faculty members of IAM came from Institute of Mathematics within CAS, which was founded in July of 1952 and was also directed by Prof. HUA.

Dynamics and vibration reduction performance of asymmetric tristable nonlinear energy sink

With its complex nonlinear dynamic behavior,the tristable system has shown excellent performance in areas such as energy harvesting and vibration suppression,and has attracted a lot of attention.In this paper,an asymmetric tristable design is proposed to improve the vibration suppression efficiency of nonlinear energy sinks(NESs)for the first time.The proposed asymmetric tristable NES(ATNES)is composed of a pair of oblique springs and a vertical spring.Then,the three stable states,symmetric and asymmetric,can be achieved by the adjustment of the distance and stiffness asymmetry of the oblique springs.The governing equations of a linear oscillator(LO)coupled with the ATNES are derived.The approximate analytical solution to the coupled system is obtained by the harmonic balance method(HBM)and verified numerically.The vibration suppression efficiency of three types of ATNES is compared.The results show that the asymmetric design can improve the efficiency of vibration reduction through comparing the chaotic motion of the NES oscillator between asymmetric steady states.In addition,compared with the symmetrical tristable NES(TNES),the ATNES can effectively control smaller structural vibrations.In other words,the ATNES can effectively solve the threshold problem of TNES failure to weak excitation.Therefore,this paper reveals the vibration reduction mechanism of the ATNES,and provides a pathway to expand the effective excitation amplitude range of the NES.

Theoretical and experimental study of a bi-stable piezoelectric energy harvester under hybrid galloping and band-limited random excitations

In the practical environment,it is very common for the simultaneous occurrence of base excitation and crosswind.Scavenging the combined energy of vibration and wind with a single energy harvesting structure is fascinating.For this purpose,the effects of the wind speed and random excitation level are investigated with the stochastic averaging method(SAM)based on the energy envelope.The results of the analytical prediction are verified with the Monte-Carlo method(MCM).The numerical simulation shows that the introduction of wind can reduce the critical excitation level for triggering an inter-well jump and make a bi-stable energy harvester(BEH)realize the performance enhancement for a weak base excitation.However,as the strength of the wind increases to a particular level,the influence of the random base excitation on the dynamic responses is weakened,and the system exhibits a periodic galloping response.A comparison between a BEH and a linear energy harvester(LEH)indicates that the BEH demonstrates inferior performance for high-speed wind.Relevant experiments are conducted to investigate the validity of the theoretical prediction and numerical simulation.The experimental findings also show that strong random excitation is favorable for the BEH in the range of low wind speeds.However,as the speed of the incoming wind is up to a particular level,the disadvantage of the BEH becomes clear and evident.

Analysis of secondary makeup characteristics of drill collar joint and prediction of downhole equivalent impact torque of Well SDTK1

Based on the three-dimensional elastic-plastic finite element analysis of the 8"(203.2 mm)drill collar joint,this paper studies the mechanical characteristics of the pin and box of NC56 drill collar joints under complex load conditions,as well as the downhole secondary makeup features,and calculates the downhole equivalent impact torque with the relative offset at the shoulder of internal and external threads.On the basis of verifying the correctness of the calculation results by using measured results in Well GT1,the prediction model of the downhole equivalent impact torque is formed and applied in the first extra-deep well with a depth over 10000 m in China(Well SDTK1).The results indicate that under complex loads,the stress distribution in drill collar joints is uneven,with relatively higher von Mises stress at the shoulder and the threads close to the shoulder.For 203.2 mm drill collar joints pre-tightened according to the make-up torque recommended by American Petroleum Institute standards,when the downhole equivalent impact torque exceeds 65 kN·m,the preload balance of the joint is disrupted,leading to secondary make-up of the joint.As the downhole equivalent impact torque increases,the relative offset at the shoulder of internal and external threads increases.The calculation results reveal that there exists significant downhole impact torque in Well SDTK1 with complex loading environment.It is necessary to use double shoulder collar joints to improve the impact torque resistance of the joint or optimize the operating parameters to reduce the downhole impact torque,and effectively prevent drilling tool failure.

Analysis of thermal management and anti-mechanical abuse of multi-functional battery modules based on magneto-sensitive shear thickening fluid

Electric vehicles(EVs)have garnered significant attention as a vital driver of economic growth and environmental sustainability.Nevertheless,ensuring the safety of high-energy batteries is now a top priority that cannot be overlooked during large-scale applications.This paper proposes an innovative active protection and cooling integrated battery module using smart materials,magneto-sensitive shear thickening fluid(MSTF),which is specifically designed to address safety threats posed by lithium-ion batteries(LIBs)exposed to harsh mechanical and environmental conditions.The theoretical framework introduces a novel approach for harnessing the smoothed-particle hydrodynamics(SPH)methodology that incorporates the intricate interplay of non-Newtonian fluid behavior,capturing the fluid-structure coupling inherent to the MSTF.This approach is further advanced by adopting an enhanced Herschel-Bulkley(H-B)model to encapsulate the intricate rheology of the MSTF under the influence of the magnetorheological effect(MRE)and shear thickening(ST)behavior.Numerical simulation results show that in the case of cooling,the MSTF is an effective cooling medium for rapidly reducing the temperature.In terms of mechanical abuse,the MSTF solidifies through actively applying the magnetic field during mechanical compression and impact within the battery module,resulting in 66%and 61.7%reductions in the maximum stress within the battery jellyroll,and 31.1%and 23%reductions in the reaction force,respectively.This mechanism effectively lowers the risk of short-circuit failure.The groundbreaking concepts unveiled in this paper for active protection battery modules are anticipated to be a valuable technological breakthrough in the areas of EV safety and lightweight/integrated design.

Dirac method for nonlinear and non-homogenous boundary value problems of plates

The boundary value problem plays a crucial role in the analytical investigation of continuum dynamics. In this paper, an analytical method based on the Dirac operator to solve the nonlinear and non-homogeneous boundary value problem of rectangular plates is proposed. The key concept behind this method is to transform the nonlinear or non-homogeneous part on the boundary into a lateral force within the governing function by the Dirac operator, which linearizes and homogenizes the original boundary, allowing one to employ the modal superposition method for obtaining solutions to reconstructive governing equations. Once projected into the modal space, the harmonic balance method(HBM) is utilized to solve coupled ordinary differential equations(ODEs)of truncated systems with nonlinearity. To validate the convergence and accuracy of the proposed Dirac method, the results of typical examples, involving nonlinearly restricted boundaries, moment excitation, and displacement excitation, are compared with those of the differential quadrature element method(DQEM). The results demonstrate that when dealing with nonlinear boundaries, the Dirac method exhibits more excellent accuracy and convergence compared with the DQEM. However, when facing displacement excitation, there exist some discrepancies between the proposed approach and simulations;nevertheless, the proposed method still accurately predicts resonant frequencies while being uniquely capable of handling nonuniform displacement excitations. Overall, this methodology offers a convenient way for addressing nonlinear and non-homogenous plate boundaries.

Parametric resonance of axially functionally graded pipes conveying pulsating fluid

Based on the generalized Hamilton's principle,the nonlinear governing equation of an axially functionally graded(AFG)pipe is established.The non-trivial equilibrium configuration is superposed by the modal functions of a simply supported beam.Via the direct multi-scale method,the response and stability boundary to the pulsating fluid velocity are solved analytically and verified by the differential quadrature element method(DQEM).The influence of Young's modulus gradient on the parametric resonance is investigated in the subcritical and supercritical regions.In general,the pipe in the supercritical region is more sensitive to the pulsating excitation.The nonlinearity changes from hard to soft,and the non-trivial equilibrium configuration introduces more frequency components to the vibration.Besides,the increasing Young's modulus gradient improves the critical pulsating flow velocity of the parametric resonance,and further enhances the stability of the system.In addition,when the temperature increases along the axial direction,reducing the gradient parameter can enhance the response asymmetry.This work further complements the theoretical analysis of pipes conveying pulsating fluid.

Natural vibration and critical velocity of translating Timoshenko beam with non-homogeneous boundaries

In most practical engineering applications,the translating belt wraps around two fixed wheels.The boundary conditions of the dynamic model are typically specified as simply supported or fixed boundaries.In this paper,non-homogeneous boundaries are introduced by the support wheels.Utilizing the translating belt as the mechanical prototype,the vibration characteristics of translating Timoshenko beam models with nonhomogeneous boundaries are investigated for the first time.The governing equations of Timoshenko beam are deduced by employing the generalized Hamilton's principle.The effects of parameters such as the radius of wheel and the length of belt on vibration characteristics including the equilibrium deformations,critical velocities,natural frequencies,and modes,are numerically calculated and analyzed.The numerical results indicate that the beam experiences deformation characterized by varying curvatures near the wheels.The radii of the wheels play a pivotal role in determining the change in trend of the relative difference between two beam models.Comparing the results unearths that the relative difference in equilibrium deformations between the two beam models is more pronounced with smaller-sized wheels.When the two wheels are of equal size,the critical velocities of both beam models reach their respective minima.In addition,the relative difference in natural frequencies between the two beam models exhibits nonlinear variation and can easily exceed 50%.Furthermore,as the axial velocities increase,the impact of non-homogeneous boundaries on modal shape of translating beam becomes more significant.Although dealing with non-homogeneous boundaries is challenging,beam models with non-homogeneous boundaries are more sensitive to parameters,and the differences between the two types of beams undergo some interesting variations under the influence of non-homogeneous boundaries.

Dynamical distribution of continuous service time model involving non-Maxwellian collision kernel and value functions

The distribution of continuous service time in call centers is investigated.A non-Maxwellian collision kernel combining two different value functions in the interaction rule are used to describe the evolution of continuous service time,respectively.Using the statistical mechanical and asymptotic limit methods,Fokker–Planck equations are derived from the corresponding Boltzmann-type equations with non-Maxwellian collision kernels.The steady-state solutions of the Fokker–Planck equation are obtained in exact form.Numerical experiments are provided to support our results under different parameters.

Prediction of single cell mechanical properties in microchannels based on deep learning

Traditional methods for measuring single-cell mechanical characteristics face several challenges,including lengthy measurement times,low throughput,and a requirement for advanced technical skills.To overcome these challenges,a novel machine learning(ML)approach is implemented based on the convolutional neural networks(CNNs),aiming at predicting cells'elastic modulus and constitutive equations from their deformations while passing through micro-constriction channels.In the present study,the computational fluid dynamics technology is used to generate a dataset within the range of the cell elastic modulus,incorporating three widely-used constitutive models that characterize the cellular mechanical behavior,i.e.,the Mooney-Rivlin(M-R),Neo-Hookean(N-H),and Kelvin-Voigt(K-V)models.Utilizing this dataset,a multi-input convolutional neural network(MI-CNN)algorithm is developed by incorporating cellular deformation data as well as the time and positional information.This approach accurately predicts the cell elastic modulus,with a coefficient of determination R^(2)of 0.999,a root mean square error of 0.218,and a mean absolute percentage error of 1.089%.The model consistently achieves high-precision predictions of the cellular elastic modulus with a maximum R^(2)of 0.99,even when the stochastic noise is added to the simulated data.One significant feature of the present model is that it has the ability to effectively classify the three types of constitutive equations we applied.The model accurately and reliably predicts single-cell mechanical properties,showcasing a robust ability to generalize.We demonstrate that incorporating deformation features at multiple time points can enhance the algorithm's accuracy and generalization.This algorithm presents a possibility for high-throughput,highly automated,real-time,and precise characterization of single-cell mechanical properties.

Dynamic performance and parameter optimization of a half-vehicle system coupled with an inerter-based X-structure nonlinear energy sink

Inspired by the demand of improving the riding comfort and meeting the lightweight design of the vehicle, an inerter-based X-structure nonlinear energy sink(IXNES) is proposed and applied in the half-vehicle system to enhance the dynamic performance. The X-structure is used as a mechanism to realize the nonlinear stiffness characteristic of the NES, which can realize the flexibility, adjustability, high efficiency, and easy operation of nonlinear stiffness, and is convenient to apply in the vehicle suspension, and the inerter is applied to replacing the mass of the NES based on the mass amplification characteristic. The dynamic model of the half-vehicle system coupled with the IX-NES is established with the Lagrange theory, and the harmonic balance method(HBM) and the pseudo-arc-length method(PALM) are used to obtain the dynamic response under road harmonic excitation. The corresponding dynamic performance under road harmonic and random excitation is evaluated by six performance indices, and compared with that of the original half-vehicle system to show the benefits of the IX-NES. Furthermore, the structural parameters of the IX-NES are optimized with the genetic algorithm. The results show that for road harmonic and random excitation, using the IX-NES can greatly reduce the resonance peaks and root mean square(RMS) values of the front and rear suspension deflections and the front and rear dynamic tire loads, while the resonance peaks and RMS values of the vehicle body vertical and pitching accelerations are slightly larger.When the structural parameters of the IX-NES are optimized, the vehicle body vertical and pitching accelerations of the half-vehicle system could reduce by 2.41% and 1.16%,respectively, and the other dynamic performance indices are within the reasonable ranges.Thus, the IX-NES combines the advantages of the inerter, X-structure, and NES, which improves the dynamic performance of the half-vehicle system and provides an effective option for vibration attenuation in the vehicle engineering.

A family of quasi-cubic blended splines and applications

A class of quasi-cubic B-spline base functions by trigonometric polynomials are established which inherit properties similar to those of cubic B-spline bases. The corresponding curves with a shape parameter a, defined by the introduced base functions, include the B-spline curves and can approximate the B-spline curves from both sides. The curves can be adjusted easily by using the shape parameter a, where dpi（a,t） is linear with respect to da for the fixed t. With the shape parameter chosen properly, the defined curves can be used to precisely represent straight line segments, parabola segments, circular arcs and some transcendental curves, and the corresponding tensor product surfaces can also represent spherical surfaces, cylindrical surfaces and some transcendental surfaces exactly. By abandoning positive property, this paper proposes a new C^2 continuous blended interpolation spline based on piecewise trigonometric polynomials associated with a sequence of local parameters. Illustration showed that the curves and surfaces constructed by the blended spline can be adjusted easily and freely. The blended interpolation spline curves can be shape-preserving with proper local parameters since these local parameters can be considered to be the magnification ratio to the length of tangent vectors at the interpolating points. The idea is extended to produce blended spline surfaces.

Nonlinear dynamics of axially moving viscoelastic Timoshenko beam under parametric and external excitations

This investigation focuses on the nonlinear dynamic behaviors in the trans- verse vibration of an axiMly accelerating viscoelastic Timoshenko beam with the external harmonic excitation. The parametric excitation is caused by the harmonic fluctuations of the axial moving speed. An integro-partial-differential equation governing the transverse vibration of the Timoshenko beam is established. Many factors are considered, such as viscoelasticity, the finite axial support rigidity, and the longitudinally varying tension due to the axial acceleration. With the Galerkin truncation method, a set of nonlinear ordinary differential equations are derived by discretizing the governing equation. Based on the numerical solutions, the bifurcation diagrams are presented to study the effect of the external transverse excitation. Moreover, the frequencies of the two excitations are assumed to be multiple. Further, five different tools, including the time history, the Poincaré map, and the sensitivity to initial conditions, are used to identify the motion form of the nonlinear vibration. Numerical results also show the characteristics of the quasiperiodic motion of the translating Timoshenko beam under an incommensurable re- lationship between the dual-frequency excitations.

NONLINEAR DYNAMICS OF AXIALLY ACCELERATING VISCOELASTIC BEAMS BASED ON DIFFERENTIAL QUADRATURE

This paper investigates nonlinear dynamical behaviors in transverse motion of an axially accelerating viscoelastic beam via the differential quadrature method. The governing equation, a nonlinear partial-differential equation, is derived from the viscoelastic constitution relation using the material derivative. The differential quadrature scheme is developed to solve numerically the governing equation. Based on the numerical solutions, the nonlinear dynamical behaviors are identified by use of the Poincare map and the phase portrait. The bifurcation diagrams are presented in the case that the mean axial speed and the amplitude of the speed fluctuation are respectively varied while other parameters are fixed. The Lyapunov exponent and the initial value sensitivity of the different points of the beam, calculated from the time series based on the numerical solutions, are used to indicate periodic motions or chaotic motions occurring in the transverse motion of the axially accelerating viscoelastic beam.

Suppression of multiple modal resonances of a cantilever beam by an impact damper

Impact dampers are usually used to suppress single mode resonance. The goal of this paper is to clarify the difference when the impact damper suppresses the resonances of different modes. A cantilever beam equipped with the impact damper is modeled. The elastic contact of the ball and the cantilever beam is described by using the Hertz contact model. The viscous damper between the ball and the cantilever beam is modeled to consume the vibrational energy of the cantilever beam. A piecewise ordinary differential-partial differential equation of the cantilever beam is established, including equations with and without the impact damper. The vibration responses of the cantilever beam with and without the impact damper are numerically calculated. The effects of the impact absorber parameters on the vibration reduction are examined. The results show that multiple resonance peaks of the cantilever beam can be effectively suppressed by the impact damper. Specifically, all resonance amplitudes can be reduced by a larger weight ball. Moreover, the impacting gap is very effective in suppressing the vibration of the cantilever beam. More importantly, there is an optimal impacting gap for each resonance mode of the cantilever beam, but the optimal gap for each mode is different.

New Family of Exact Solutions and Chaotic Soltions of Generalized Breor-Kaup System in （2＋1）-Dimensions via an Extended Mapping Approach

Starting from an extended mapping approach, a new type of variable separation solution with arbitrary functions of generalized (2+1)-dimensional Broer-Kaup system (GBK) system is derived. Then based on the derived solitary wave solution, we obtain some specific chaotic solitons to the (2+1)-dimensional GBK system.

COMPUTATION OF VECTOR VALUED BLENDING RATIONAL INTERPOLANTS

As we know, Newton's interpolation polynomial is based on divided differ-ences which can be calculated recursively by the divided-difference scheme while Thiele'sinterpolating continued fractions are geared towards determining a rational functionwhich can also be calculated recursively by so-called inverse differences. In this paper,both Newton's interpolation polynomial and Thiele's interpolating continued fractionsare incorporated to yield a kind of bivariate vector valued blending rational interpolantsby means of the Samelson inverse. Blending differences are introduced to calculate theblending rational interpolants recursively, algorithm and matrix-valued case are dis-cussed and a numerical example is given to illustrate the efficiency of the algorithm.

ZERO DISSIPATION LIMIT OF THE COMPRESSIBLE HEAT-CONDUCTING NAVIER-STOKES EQUATIONS IN THE PRESENCE OF THE SHOCK

The zero dissipation limit of the compressible heat-conducting Navier–Stokes equations in the presence of the shock is investigated. It is shown that when the heat conduction coefficient κ and the viscosity coefficient ε satisfy κ = O（ε）, κ/ε≥ c 〉 0, as ε→ 0 （see （1.3））, if the solution of the corresponding Euler equations is piecewise smooth with shock wave satisfying the Lax entropy condition, then there exists a smooth solution to the Navier–Stokes equations, which converges to the piecewise smooth shock solution of the Euler equations away from the shock discontinuity at a rate of ε. The proof is given by a combination of the energy estimates and the matched asymptotic analysis introduced in [3].

Vibration reduction evaluation of a linear system with a nonlinear energy sink under a harmonic and random excitation

The nonlinear behaviors and vibration reduction of a linear system with a nonlinear energy sink(NES)are investigated.The linear system is excited by a harmonic and random base excitation,consisting of a mass block,a linear spring,and a linear viscous damper.The NES is composed of a mass block,a linear viscous damper,and a spring with ideal cubic nonlinear stiffness.Based on the generalized harmonic function method,the steady-state Fokker-Planck-Kolmogorov equation is presented to reveal the response of the system.The path integral method based on the Gauss-Legendre polynomial is used to achieve the numerical solutions.The performance of vibration reduction is evaluated by the displacement and velocity transition probability densities,the transmissibility transition probability density,and the percentage of the energy absorption transition probability density of the linear oscillator.The sensitivity of the parameters is analyzed for varying the nonlinear stiffness coefficient and the damper ratio.The investigation illustrates that a linear system with NES can also realize great vibration reduction under harmonic and random base excitations and random bifurcation may appear under different parameters,which will affect the stability of the system.

Primary resonance of traveling viscoelastic beam under internal resonance

Under the 3：1 internal resonance condition, the steady-state periodic response of the forced vibration of a traveling viscoelastic beam is studied. The viscoelastic behaviors of the traveling beam are described by the standard linear solid model, and the material time derivative is adopted in the viscoelastic constitutive relation. The direct multi-scale method is used to derive the relationships between the excitation frequency and the response amplitudes. For the first time, the real modal functions are employed to analytically investigate the periodic response of the axially traveling beam. The unde- termined coefficient method is used to approximately establish the real modal functions. The approximate analytical results are confirmed by the Galerkin truncation. Numerical examples are presented to highlight the effects of the viscoelastic behaviors on the steady-state periodic responses. To illustrate the effect of the internal resonance, the energy transfer between the internal resonance modes and the saturation-like phenomena in the steady-state responses is presented.