The behavior of two immiscible low-viscosity liquids differing in density and viscosity in a vertical flat layer undergoing modulated rotation is experimentally studied.The layer has a circular axisymmetric boundary.I...The behavior of two immiscible low-viscosity liquids differing in density and viscosity in a vertical flat layer undergoing modulated rotation is experimentally studied.The layer has a circular axisymmetric boundary.In the absence of modulation of the rotation speed,the interphase boundary has the shape of a short axisymmetric cylinder.A new effect has been discovered,under the influence of rotation speed modulation,the interface takes on a new dynamic equilibrium state.A more viscous liquid covers the end boundaries of the layer in the form of thin films,which have the shape of round spots of almost constant radius;with increasing amplitude of the velocity modulation,the wetting boundary expands.It is found that upon reaching the critical amplitude of oscillations,the film of a viscous liquid loses stability,and the outer edge of the wetting spot collapses and takes on a feathery structure.It is shown that this threshold is caused by the development of the Kelvin-Helmholtz oscillatory instability of the film.The spreading radius of a spot of light viscous liquid and its stability are studied depending on the rotation rate,amplitude,and frequency of rotation speed modulation.The discovered averaged effects are determined by different oscillatory interaction of fluids with the end-walls of the cell,due to different viscosities.The effect of films forming can find application in technological processes to intensify mass transfer at interphase boundaries.展开更多
A novel single-step method is proposed for the analysis of dynamic response of visco-elastic structures containing non-smooth contactable interfaces. In the method, a two-level algorithm is employed for dealing with a...A novel single-step method is proposed for the analysis of dynamic response of visco-elastic structures containing non-smooth contactable interfaces. In the method, a two-level algorithm is employed for dealing with a nonlinear boundary condition caused by the dynamic contact of interfaces. At the first level, an explicit method is adopted to calculate nodal displacements of global viscoelastic system without considering the effect of dynamic contact of interfaces and at the second level, by introducing contact conditions of interfaces, a group of equations of lower order is derived to calculate dynamic contact normal and shear forces on the interfaces. The method is convenient and efficient for the analysis of problems of dynamic contact. The accuracy of the method is of the second order and the numerical stability condition is wider than that of other explicit methods.展开更多
The impact dynamics of a flexible multibody system is investigated. By using a partition method, the system is divided into two parts, the local impact region and the region away from the impact. The two parts are con...The impact dynamics of a flexible multibody system is investigated. By using a partition method, the system is divided into two parts, the local impact region and the region away from the impact. The two parts are connected by specific boundary conditions, and the system after partition is equivalent to the original system. According to the rigid-flexible coupling dynamic theory of multibody system, system's rigid-flexible coupling dynamic equations without impact are derived. A local impulse method for establishing the initial impact conditions is proposed. It satisfies the compatibility con- ditions for contact constraints and the actual physical situation of the impact process of flexible bodies. Based on the contact constraint method, system's impact dynamic equa- tions are derived in a differential-algebraic form. The contact/separation criterion and the algorithm are given. An impact dynamic simulation is given. The results show that system's dynamic behaviors including the energy, the deformations, the displacements, and the impact force during the impact process change dramatically. The impact makes great effects on the global dynamics of the system during and after impact.展开更多
In order to investigate the mechanical behavior of systems with complex architecture and a large number of contacting bodies,a finite element software,named LiToTac,has been developed by using the object-oriented prog...In order to investigate the mechanical behavior of systems with complex architecture and a large number of contacting bodies,a finite element software,named LiToTac,has been developed by using the object-oriented programming technique.This software,with an interactive graphical user interface,is able to handle highly non-linear problems including multiple contacts and large deformation.More importantly,the contact detection based on a hybrid three-stages methodology can be performed automatically,which is more efficient than the common strategies of pre-defining contact zones in commercial FEM software like ANSYS,ABAQUS,etc.In addition,the contact solver in LiToTac is portable between dynamic and quasi-static codes and can accurately solve contact coupled with friction in a reduced system.Several numerical examples are carried out to illustrate the functionality and capacity of the software package.展开更多
There are five most widely used contact angle schemes in the pseudopotential lattice Boltzmann(LB)model for simulating the wetting phenomenon:The pseudopotential-based scheme(PB scheme),the improved virtualdensity sch...There are five most widely used contact angle schemes in the pseudopotential lattice Boltzmann(LB)model for simulating the wetting phenomenon:The pseudopotential-based scheme(PB scheme),the improved virtualdensity scheme(IVD scheme),the modified pseudopotential-based scheme with a ghost fluid layer constructed by using the fluid layer density above the wall(MPB-C scheme),the modified pseudopotential-based scheme with a ghost fluid layer constructed by using the weighted average density of surrounding fluid nodes(MPB-W scheme)and the geometric formulation scheme(GF scheme).But the numerical stability and accuracy of the schemes for wetting simulation remain unclear in the past.In this paper,the numerical stability and accuracy of these schemes are clarified for the first time,by applying the five widely used contact angle schemes to simulate a two-dimensional(2D)sessile droplet on wall and capillary imbibition in a 2D channel as the examples of static wetting and dynamic wetting simulations respectively.(i)It is shown that the simulated contact angles by the GF scheme are consistent at different density ratios for the same prescribed contact angle,but the simulated contact angles by the PB scheme,IVD scheme,MPB-C scheme and MPB-W scheme change with density ratios for the same fluid-solid interaction strength.The PB scheme is found to be the most unstable scheme for simulating static wetting at increased density ratios.(ii)Although the spurious velocity increases with the increased liquid/vapor density ratio for all the contact angle schemes,the magnitude of the spurious velocity in the PB scheme,IVD scheme and GF scheme are smaller than that in the MPB-C scheme and MPB-W scheme.(iii)The fluid density variation near the wall in the PB scheme is the most significant,and the variation can be diminished in the IVD scheme,MPB-C scheme andMPBWscheme.The variation totally disappeared in the GF scheme.(iv)For the simulation of capillary imbibition,the MPB-C scheme,MPB-Wscheme and GF scheme simulate the dynamics of the liquid-vapor interface well,with the GF scheme being the most accurate.The accuracy of the IVD scheme is low at a small contact angle(44 degrees)but gets high at a large contact angle(60 degrees).However,the PB scheme is the most inaccurate in simulating the dynamics of the liquid-vapor interface.As a whole,it is most suggested to apply the GF scheme to simulate static wetting or dynamic wetting,while it is the least suggested to use the PB scheme to simulate static wetting or dynamic wetting.展开更多
Dynamic wetting plays an important role in the physics of multiphase flow, and has a significant influence on many industrial and geotechnical applications. In this work, a modified smoothed particle hydrodynamics (SP...Dynamic wetting plays an important role in the physics of multiphase flow, and has a significant influence on many industrial and geotechnical applications. In this work, a modified smoothed particle hydrodynamics (SPH) model is employed to simulate surface tension, contact angle and dynamic wetting effects at meso-scale. The wetting and dewetting phenomena are simulated in a capillary tube, where the liquid particles are raised or withdrawn by a shifting substrate. The SPH model is modified by introducing a newly developed viscous force formulation at the liquid-solid interface to reproduce the rate-dependent behaviour of the moving contact line. Dynamic contact angle simulations with the interfacial viscous force are conducted to verify the effectiveness and accuracy of this new formulation. In addition, the influence of interfacial viscous forces with different magnitude on the contact angle dynamics is examined by empirical power-law correlations;the derived constants suggest that the dynamic contact angle changes monotonically with the interfacial viscous force. The simulation results are consistent with experimental observations and theoretical predictions, implying that the interfacial viscous force can be associated with the slip length of flow and the microscopic surface roughness. This work demonstrates that the modified SPH model can successfully account for the rate-dependent effects of a moving contact line, and can be used for realistic multiphase flow simulations under dynamic conditions.展开更多
Under the frame of multibody dynamics, the contact dynamics of elasto-plastic spatial thin beams is numerically studied by using the spatial thin beam elements of absolute nodal coordinate formulation(ANCF). The int...Under the frame of multibody dynamics, the contact dynamics of elasto-plastic spatial thin beams is numerically studied by using the spatial thin beam elements of absolute nodal coordinate formulation(ANCF). The internal force of the elasto-plastic spatial thin beam element is derived under the assumption that the plastic strain of the beam element depends only on its longitudinal deformation.A new body-fixed local coordinate system is introduced into the spatial thin beam element of ANCF for efficient contact detection in the contact dynamics simulation. The linear isotropic hardening constitutive law is used to describe the elasto-plastic deformation of beam material, and the classical return mapping algorithm is adopted to evaluate the plastic strains. A multi-zone contact approach of thin beams previously proposed by the authors is also introduced to detect the multiple contact zones of beams accurately, and the penalty method is used to compute the normal contact force of thin beams in contact. Four numerical examples are given to demonstrate the applicability and effectiveness of the proposed elasto-plastic spatial thin beam element of ANCF for flexible multibody system dynamics.展开更多
Computational e ciency and accuracy always conflict with each other in molecular dynamics(MD) simulations. How to enhance the computational e ciency and keep accuracy at the same time is concerned by each correspondin...Computational e ciency and accuracy always conflict with each other in molecular dynamics(MD) simulations. How to enhance the computational e ciency and keep accuracy at the same time is concerned by each corresponding researcher. However, most of the current studies focus on MD algorithms, and if the scale of MD model could be reduced, the algorithms would be more meaningful. A local region molecular dynamics(LRMD) simulation method which can meet these two factors concurrently in nanoscale sliding contacts is developed in this paper. Full MD simulation is used to simulate indentation process before sliding. A criterion called contribution of displacement is presented, which is used to determine the e ective local region in the MD model after indentation. By using the local region, nanoscale sliding contact between a rigid cylindrical tip and an elastic substrate is investigated. Two two?dimensional MD models are presented, and the friction forces from LRMD simulations agree well with that from full MD simulations, which testifies the e ectiveness of the LRMD simulation method for two?dimensional cases. A three?dimensional MD model for sliding contacts is developed then to show the validity of the LRMD simulation method further. Finally, a discussion is carried out by the principles of tribology. In the discussion, two two?dimensional full MD models are used to simulate the nanoscale sliding contact problems. The results indicate that original smaller model will induce higher equivalent scratching depth, and then results in higher friction forces, which will help to explain the mechanism how the LRMD simulation method works. This method can be used to reduce the scale of MD model in large scale simulations, and it will enhance the computational e ciency without losing accuracy during the simula?tion of nanoscale sliding contacts.展开更多
In nanoscale sliding contact,adhesion effects and adhesive force are predominant,and high friction force will be produced.Friction energy is mainly converted into heat,and the heat will make nanomaterials become soft ...In nanoscale sliding contact,adhesion effects and adhesive force are predominant,and high friction force will be produced.Friction energy is mainly converted into heat,and the heat will make nanomaterials become soft to affect friction behaviors,so it is important to investigate the friction and thermal properties of the nanoscale sliding contacts.A model of a nanoscale sliding contact between a rigid cylindrical tip and an FCC copper substrate is developed by molecular dynamics simulation.The thermal properties of the substrate and the friction behaviors are studied at different sliding velocities and different tip radii.The results show that at a low sliding velocity,the friction force fluctuation is mainly caused by material melting⁃solidification,while at a high sliding velocity the material melting is a main factor for the friction reduction.The average friction forces increase at initial phase and then decrease with increasing sliding velocity,and the average temperature of the substrate increases as sliding velocity increases.Increasing tip radius significantly increases the temperature,while the coupled effects of tip radius and temperature rise make friction force increase slightly.展开更多
Usage of methods for determination of long-term trends of the dynamics of vegetation formation on the background of climate changes becomes more and more important at modern stage of the assessment of natural systems ...Usage of methods for determination of long-term trends of the dynamics of vegetation formation on the background of climate changes becomes more and more important at modern stage of the assessment of natural systems development. This causes to researchers a series of problems from choice of conceptual base to notions and terms of the processes identified and of state of vegetation of different environments. Solution of such a task results inevitably in necessity to correct the understanding of existing processes occurring in the vegetation cover. It allows establishing a direction of their development in the system of natural factors of any territory. As a result, we have base for determination of age, site and role of current state of phytocenoses in successional systems. It is necessary for this to reveal the peculiarities of phytocenoses composition and formation due to climate dynamics and to determine a period of ecosystems homeostasis, especially for phytocenoses of contact natural conditions. Optimal values of phytocenoses diversity, like one of whole biosystems, depend on the amount of resource in the environment, on stability degree and on evolutional development of cenoses due to ecological, micro-evolutional and evolutional processes.展开更多
In order to implement the dynamic characteristic of a dual power-split transmission, a dynamic me- chanics model is built. Firstly, according to the method of theoretical analysis of the tooth contact analysis (TCA)...In order to implement the dynamic characteristic of a dual power-split transmission, a dynamic me- chanics model is built. Firstly, according to the method of theoretical analysis of the tooth contact analysis (TCA) and loaded tooth contact analysis (LTCA), the actual meshing process of each gear pairs is simulated, and the time-varying mesh stiffness excitations are obtained, which can improve the numerical precision. Second- ly, by using the lumped mass method, the bending-torsional coupling three dimensional dynamic model of the dual power-split transmission is established, and the identical dimensionless equations are deduced by elimina- ting the effect of rigid displacement and the method of dimensional normalization. Finally, by the method of the fourth order Runge-Kutta algorithm with variable step lengths, the responses of this system in a frequency domain and time domain are obtained, and the dynamic load change characteristics of each gear pairs are analyzed. The results show that the establishment, solution and analysis of the system dynamics model could provide a basis for the dynamic design, and have an important significance for the dynamic efficiency analysis and dynamic perform- ance optimization design of the dual power-split transmission.展开更多
The effect of involute contact ratio on the torsional vibration behavior ofspur gear-pair is studied analytically through a mass-spring model. The tooth stiffness in model notonly has a relation with time, as many pri...The effect of involute contact ratio on the torsional vibration behavior ofspur gear-pair is studied analytically through a mass-spring model. The tooth stiffness in model notonly has a relation with time, as many prior studies presented, but, more important, with involutecontact ratio (ICR) as well. The ICR embodies its impact on the spur gear's dynamic performancethrough changing the proportion of tooth stiffness when there are n+1 teeth in contact to toothstiffness when there are n teeth in contact. A couple of curves about impact of ICR on the gear'sdynamic performance are presented, and they clearly demonstrate that the model can accuratelydescribe the effects of ICR on dynamic transmission error. Finally, some conclusions useful toreduce vibration and noise of gear-pair are proposed.展开更多
The effect of rail corrugation on the vertical dynamics of railway vehicle coupled with a curved track is investigated in detail with a numerical method when a wheelset is steadily curving. In the calculation of rail ...The effect of rail corrugation on the vertical dynamics of railway vehicle coupled with a curved track is investigated in detail with a numerical method when a wheelset is steadily curving. In the calculation of rail corrugation we consider the combination of Kalkers rolling contact theory modified, a model of material loss on rail running surface, and a dynamics model of railway vehicle coupled with a curved track. In the establishment of the dynamic model, for simplicity, one fourth of the freight car without lateral motions, namely a wheelset and the equivalent one fourth freight car body above it, is considered. The Euler beam is used to model the rails and the track structure under the rails is replaced with equivalent springs, dampers and mass bodies. The numerical results show the great influence of the rail corrugation on the vibration of the parts of the vehicle and the track, and the some characters of rail corrugation in development.展开更多
Soft machine refers to a kind of mechanical system made of soft materials to complete sophisticated missions, such as handling a fragile object and crawling along a narrow tunnel corner, under low cost control and act...Soft machine refers to a kind of mechanical system made of soft materials to complete sophisticated missions, such as handling a fragile object and crawling along a narrow tunnel corner, under low cost control and actuation. Hence, soft machines have raised great challenges to computational dynamics. In this review article, recent studies of the authors on the dynamic modeling, numerical simulation, and experimental validation of soft machines are summarized in the framework of multibody system dynamics. The dynamic modeling approaches are presented first for the geometric nonlinearities of coupled overall motions and large deformations of a soft component, the physical nonlinearities of a soft component made of hyperelastic or elastoplastic materials, and the frictional contacts/impacts of soft components, respectively. Then the computation approach is outlined for the dynamic simulation of soft machines governed by a set of differential-algebraic equations of very high dimensions, with an emphasis on the efficient computations of the nonlinear elastic force vector of finite elements. The validations of the proposed approaches are given via three case studies, including the locomotion of a soft quadrupedal robot, the spinning deployment of a solar sail of a spacecraft, and the deployment of a mesh reflector of a satellite antenna, as well as the corresponding experimental studies. Finally, some remarks are made for future studies.展开更多
When characterizing flows in miniaturized channels, the determination of the dynamic contact angle is important. By measuring the dynamic contact angle, the flow properties of the flowing liquid and the effect of mate...When characterizing flows in miniaturized channels, the determination of the dynamic contact angle is important. By measuring the dynamic contact angle, the flow properties of the flowing liquid and the effect of material properties on the flow can be characterized. A machine vision based system to measure the contact angle of front or rear menisci of a moving liquid plug is described in this article. In this research, transparent flow channels fabricated on thermoplastic polymer and sealed with an adhesive tape are used. The transparency of the channels enables image based monitoring and measurement of flow variables, including the dynamic contact angle. It is shown that the dynamic angle can be measured from a liquid flow in a channel using the image based measurement system. An image processing algorithm has been developed in a MATLAB environment. Images are taken using a CCD camera and the channels are illuminated using a custom made ring light. Two fitting methods, a circle and two parabolas, are experimented and the results are compared in the measurement of the dynamic contact angles.展开更多
Considering the discontinuous characteristics of sea ice on various scales,a modified discrete element model(DEM) for sea ice dynamics is developed based on the granular material rheology.In this modified DEM,a soft...Considering the discontinuous characteristics of sea ice on various scales,a modified discrete element model(DEM) for sea ice dynamics is developed based on the granular material rheology.In this modified DEM,a soft sea ice particle element is introduced as a self-adjustive particle size function.Each ice particle can be treated as an assembly of ice floes,with its concentration and thickness changing to variable sizes under the conservation of mass.In this model,the contact forces among ice particles are calculated using a viscous-elastic-plastic model,while the maximum shear forces are described with the Mohr-Coulomb friction law.With this modified DEM,the ice flow dynamics is simulated under the drags of wind and current in a channel of various widths.The thicknesses,concentrations and velocities of ice particles are obtained,and then reasonable dynamic process is analyzed.The sea ice dynamic process is also simulated in a vortex wind field.Taking the influence of thermodynamics into account,this modified DEM will be improved in the future work.展开更多
Fractal theory provides scale?independent asperity contact loads and assumes variable curvature radii in the contact analyses of rough surfaces, the current research for which mainly focuses on the mechanism study. Th...Fractal theory provides scale?independent asperity contact loads and assumes variable curvature radii in the contact analyses of rough surfaces, the current research for which mainly focuses on the mechanism study. The present study introduces the fractal theory into the dynamic research of gas face seals under face?contacting conditions. Structure?Function method is adopted to handle the surface profiles of typical carbon?graphite rings, proving the fractal con?tact model can be used in the field of gas face seals. Using a numerical model established for the dynamic analyses of a spiral groove gas face seal with a flexibly mounted stator, a comparison of dynamic performance between the Majumdar?Bhushan(MB) fractal model and the Chang?Etsion?Bogy(CEB) statistical model is performed. The result shows that the two approaches induce differences in terms of the occurrence and the level of face contact. Although the approach distinctions in film thickness and leakage rate can be tiny, the distinctions in contact mechanism and end face damage are obvious. An investigation of fractal parameters D and G shows that a proper D(nearly 1.5) and a small G are helpful in raising the proportion of elastic deformation to weaken the adhesive wear in the sealing dynamic performance. The proposed research provides a fractal approach to design gas face seals.展开更多
The existence of clearance in the joints of mechanisms system is inevitable.The movements of the real mechanism are deflection from the ideal mechanism due to the clearances and the motion accuracy is decreased.The ef...The existence of clearance in the joints of mechanisms system is inevitable.The movements of the real mechanism are deflection from the ideal mechanism due to the clearances and the motion accuracy is decreased.The effects of the hinge clearance on the crank and rocker mechanism system are studied.The system dynamics equation with clearance is presented.The contact dynamics model is established using the nonlinear equivalent spring-damp model and the friction effect is considered by using Coulomb friction model.Then the models are incorporated into ADAMS,and based on the model,large numbers numeric simulations are made.The regularity of contact forces in clearance are studied in detail.And the effects of clearance size,clearance friction on the mechanism dynamics characteristic are analyzed.The simulation results can predict the effects of clearance on the mechanism dynamics characteristic preferably.展开更多
Dynamic contact stiffness at the interface between a vibrating rigid sphere and a semi-infinite transversely isotropic viscoelastic solid is investigated. An oscillating force superimposed onto a static compressive fo...Dynamic contact stiffness at the interface between a vibrating rigid sphere and a semi-infinite transversely isotropic viscoelastic solid is investigated. An oscillating force superimposed onto a static compressive force in the vertical direction excites the vibration of a rigid sphere, which causes variable contact radius and contact pressure distribution in the contact region. The assumption of a sufficiently small oscillating force yields a dynamic contact-pressure distribution of a constant contact radius, which gives dynamic contact stiffness at the interface between the rigid sphere and the semi-infinite solid. Numerical calculations show the influence of vibration frequency of the sphere, and elastic constants of the transversely isotropic solid on dynamic contact stiffness, which benefits quantitative evaluation of elastic constants and orientation of single hexagonal grains by resonance-frequency shifts of the oscillator in resonance ultrasound microscopy.展开更多
In order to describe the performance of thin wall bearing on rotor system more accurate,the simplified model of bearing local stiffness was proposed. The load distribution and local contact deformation in angular cont...In order to describe the performance of thin wall bearing on rotor system more accurate,the simplified model of bearing local stiffness was proposed. The load distribution and local contact deformation in angular contact ball bearings were calculated using quasi dynamic calculation method. Based on the relationship of local load to contact deformation,the calculation model of local bearing stiffness was subsequently built to get radial and axial components of local stiffness. Effects of external loads on the local bearing stiffness were analyzed. The results showed that local stiffness in bearings is symmetric to the axis of radial load,and its value has a maximum on the symmetry axis along the radial load direction. External radial and axial load have different effects on local bearing stiffness.展开更多
基金financially supported by the Russian Science Foundation(Project 23-11-00242).
文摘The behavior of two immiscible low-viscosity liquids differing in density and viscosity in a vertical flat layer undergoing modulated rotation is experimentally studied.The layer has a circular axisymmetric boundary.In the absence of modulation of the rotation speed,the interphase boundary has the shape of a short axisymmetric cylinder.A new effect has been discovered,under the influence of rotation speed modulation,the interface takes on a new dynamic equilibrium state.A more viscous liquid covers the end boundaries of the layer in the form of thin films,which have the shape of round spots of almost constant radius;with increasing amplitude of the velocity modulation,the wetting boundary expands.It is found that upon reaching the critical amplitude of oscillations,the film of a viscous liquid loses stability,and the outer edge of the wetting spot collapses and takes on a feathery structure.It is shown that this threshold is caused by the development of the Kelvin-Helmholtz oscillatory instability of the film.The spreading radius of a spot of light viscous liquid and its stability are studied depending on the rotation rate,amplitude,and frequency of rotation speed modulation.The discovered averaged effects are determined by different oscillatory interaction of fluids with the end-walls of the cell,due to different viscosities.The effect of films forming can find application in technological processes to intensify mass transfer at interphase boundaries.
基金The project supported by the National Natural Science Foundation of China(59578032)the Key Project of the Ninth Five-Year Plan(96221030202)
文摘A novel single-step method is proposed for the analysis of dynamic response of visco-elastic structures containing non-smooth contactable interfaces. In the method, a two-level algorithm is employed for dealing with a nonlinear boundary condition caused by the dynamic contact of interfaces. At the first level, an explicit method is adopted to calculate nodal displacements of global viscoelastic system without considering the effect of dynamic contact of interfaces and at the second level, by introducing contact conditions of interfaces, a group of equations of lower order is derived to calculate dynamic contact normal and shear forces on the interfaces. The method is convenient and efficient for the analysis of problems of dynamic contact. The accuracy of the method is of the second order and the numerical stability condition is wider than that of other explicit methods.
基金supported by the National Natural Science Foundation of China(Nos.11132007,11272155,and 10772085)the Fundamental Research Funds for the Central Universities(No.30920130112009)the 333 Project of Jiangsu Province of China(No.BRA2011172)
文摘The impact dynamics of a flexible multibody system is investigated. By using a partition method, the system is divided into two parts, the local impact region and the region away from the impact. The two parts are connected by specific boundary conditions, and the system after partition is equivalent to the original system. According to the rigid-flexible coupling dynamic theory of multibody system, system's rigid-flexible coupling dynamic equations without impact are derived. A local impulse method for establishing the initial impact conditions is proposed. It satisfies the compatibility con- ditions for contact constraints and the actual physical situation of the impact process of flexible bodies. Based on the contact constraint method, system's impact dynamic equa- tions are derived in a differential-algebraic form. The contact/separation criterion and the algorithm are given. An impact dynamic simulation is given. The results show that system's dynamic behaviors including the energy, the deformations, the displacements, and the impact force during the impact process change dramatically. The impact makes great effects on the global dynamics of the system during and after impact.
基金We gratefully acknowledge the financial support of the National Key R&D Program of China(Grant No.2017YFB0703200)the National Natural Science Foundation of China(Grant No.11772274).
文摘In order to investigate the mechanical behavior of systems with complex architecture and a large number of contacting bodies,a finite element software,named LiToTac,has been developed by using the object-oriented programming technique.This software,with an interactive graphical user interface,is able to handle highly non-linear problems including multiple contacts and large deformation.More importantly,the contact detection based on a hybrid three-stages methodology can be performed automatically,which is more efficient than the common strategies of pre-defining contact zones in commercial FEM software like ANSYS,ABAQUS,etc.In addition,the contact solver in LiToTac is portable between dynamic and quasi-static codes and can accurately solve contact coupled with friction in a reduced system.Several numerical examples are carried out to illustrate the functionality and capacity of the software package.
基金sponsored by the National Natural Science Foundation of China under Grant No.52206101Shanghai Sailing Program under Grant No.20YF1431200the Experiments for Space Exploration Program and the Qian Xuesen Laboratory,China Academy of Space Technology under Grant No.TKTSPY-2020-01-01.
文摘There are five most widely used contact angle schemes in the pseudopotential lattice Boltzmann(LB)model for simulating the wetting phenomenon:The pseudopotential-based scheme(PB scheme),the improved virtualdensity scheme(IVD scheme),the modified pseudopotential-based scheme with a ghost fluid layer constructed by using the fluid layer density above the wall(MPB-C scheme),the modified pseudopotential-based scheme with a ghost fluid layer constructed by using the weighted average density of surrounding fluid nodes(MPB-W scheme)and the geometric formulation scheme(GF scheme).But the numerical stability and accuracy of the schemes for wetting simulation remain unclear in the past.In this paper,the numerical stability and accuracy of these schemes are clarified for the first time,by applying the five widely used contact angle schemes to simulate a two-dimensional(2D)sessile droplet on wall and capillary imbibition in a 2D channel as the examples of static wetting and dynamic wetting simulations respectively.(i)It is shown that the simulated contact angles by the GF scheme are consistent at different density ratios for the same prescribed contact angle,but the simulated contact angles by the PB scheme,IVD scheme,MPB-C scheme and MPB-W scheme change with density ratios for the same fluid-solid interaction strength.The PB scheme is found to be the most unstable scheme for simulating static wetting at increased density ratios.(ii)Although the spurious velocity increases with the increased liquid/vapor density ratio for all the contact angle schemes,the magnitude of the spurious velocity in the PB scheme,IVD scheme and GF scheme are smaller than that in the MPB-C scheme and MPB-W scheme.(iii)The fluid density variation near the wall in the PB scheme is the most significant,and the variation can be diminished in the IVD scheme,MPB-C scheme andMPBWscheme.The variation totally disappeared in the GF scheme.(iv)For the simulation of capillary imbibition,the MPB-C scheme,MPB-Wscheme and GF scheme simulate the dynamics of the liquid-vapor interface well,with the GF scheme being the most accurate.The accuracy of the IVD scheme is low at a small contact angle(44 degrees)but gets high at a large contact angle(60 degrees).However,the PB scheme is the most inaccurate in simulating the dynamics of the liquid-vapor interface.As a whole,it is most suggested to apply the GF scheme to simulate static wetting or dynamic wetting,while it is the least suggested to use the PB scheme to simulate static wetting or dynamic wetting.
文摘Dynamic wetting plays an important role in the physics of multiphase flow, and has a significant influence on many industrial and geotechnical applications. In this work, a modified smoothed particle hydrodynamics (SPH) model is employed to simulate surface tension, contact angle and dynamic wetting effects at meso-scale. The wetting and dewetting phenomena are simulated in a capillary tube, where the liquid particles are raised or withdrawn by a shifting substrate. The SPH model is modified by introducing a newly developed viscous force formulation at the liquid-solid interface to reproduce the rate-dependent behaviour of the moving contact line. Dynamic contact angle simulations with the interfacial viscous force are conducted to verify the effectiveness and accuracy of this new formulation. In addition, the influence of interfacial viscous forces with different magnitude on the contact angle dynamics is examined by empirical power-law correlations;the derived constants suggest that the dynamic contact angle changes monotonically with the interfacial viscous force. The simulation results are consistent with experimental observations and theoretical predictions, implying that the interfacial viscous force can be associated with the slip length of flow and the microscopic surface roughness. This work demonstrates that the modified SPH model can successfully account for the rate-dependent effects of a moving contact line, and can be used for realistic multiphase flow simulations under dynamic conditions.
基金supported in part by the National Natural Science Foundation of China (Grants 11290151 and 11221202)supported in part by the Beijing Higher Education Young Elite Teacher Project (Grant YETP1201)
文摘Under the frame of multibody dynamics, the contact dynamics of elasto-plastic spatial thin beams is numerically studied by using the spatial thin beam elements of absolute nodal coordinate formulation(ANCF). The internal force of the elasto-plastic spatial thin beam element is derived under the assumption that the plastic strain of the beam element depends only on its longitudinal deformation.A new body-fixed local coordinate system is introduced into the spatial thin beam element of ANCF for efficient contact detection in the contact dynamics simulation. The linear isotropic hardening constitutive law is used to describe the elasto-plastic deformation of beam material, and the classical return mapping algorithm is adopted to evaluate the plastic strains. A multi-zone contact approach of thin beams previously proposed by the authors is also introduced to detect the multiple contact zones of beams accurately, and the penalty method is used to compute the normal contact force of thin beams in contact. Four numerical examples are given to demonstrate the applicability and effectiveness of the proposed elasto-plastic spatial thin beam element of ANCF for flexible multibody system dynamics.
基金National Natural Science Foundation of China(Grant Nos.51675429,51205313)Fundamental Research Funds for the Central Universities of China(Grant No.3102014JCS05009)111 Project of China(Grant No.B13044)
文摘Computational e ciency and accuracy always conflict with each other in molecular dynamics(MD) simulations. How to enhance the computational e ciency and keep accuracy at the same time is concerned by each corresponding researcher. However, most of the current studies focus on MD algorithms, and if the scale of MD model could be reduced, the algorithms would be more meaningful. A local region molecular dynamics(LRMD) simulation method which can meet these two factors concurrently in nanoscale sliding contacts is developed in this paper. Full MD simulation is used to simulate indentation process before sliding. A criterion called contribution of displacement is presented, which is used to determine the e ective local region in the MD model after indentation. By using the local region, nanoscale sliding contact between a rigid cylindrical tip and an elastic substrate is investigated. Two two?dimensional MD models are presented, and the friction forces from LRMD simulations agree well with that from full MD simulations, which testifies the e ectiveness of the LRMD simulation method for two?dimensional cases. A three?dimensional MD model for sliding contacts is developed then to show the validity of the LRMD simulation method further. Finally, a discussion is carried out by the principles of tribology. In the discussion, two two?dimensional full MD models are used to simulate the nanoscale sliding contact problems. The results indicate that original smaller model will induce higher equivalent scratching depth, and then results in higher friction forces, which will help to explain the mechanism how the LRMD simulation method works. This method can be used to reduce the scale of MD model in large scale simulations, and it will enhance the computational e ciency without losing accuracy during the simula?tion of nanoscale sliding contacts.
基金Sponsored by the National Natural Science Foundation of China(Grant Nos.52075444,51675429)the Key Project of National Natural Science Foundation of China(Grant No.51535009)the Fundamental Research Funds for the Central Universities(Grant No.31020190503004).
文摘In nanoscale sliding contact,adhesion effects and adhesive force are predominant,and high friction force will be produced.Friction energy is mainly converted into heat,and the heat will make nanomaterials become soft to affect friction behaviors,so it is important to investigate the friction and thermal properties of the nanoscale sliding contacts.A model of a nanoscale sliding contact between a rigid cylindrical tip and an FCC copper substrate is developed by molecular dynamics simulation.The thermal properties of the substrate and the friction behaviors are studied at different sliding velocities and different tip radii.The results show that at a low sliding velocity,the friction force fluctuation is mainly caused by material melting⁃solidification,while at a high sliding velocity the material melting is a main factor for the friction reduction.The average friction forces increase at initial phase and then decrease with increasing sliding velocity,and the average temperature of the substrate increases as sliding velocity increases.Increasing tip radius significantly increases the temperature,while the coupled effects of tip radius and temperature rise make friction force increase slightly.
文摘Usage of methods for determination of long-term trends of the dynamics of vegetation formation on the background of climate changes becomes more and more important at modern stage of the assessment of natural systems development. This causes to researchers a series of problems from choice of conceptual base to notions and terms of the processes identified and of state of vegetation of different environments. Solution of such a task results inevitably in necessity to correct the understanding of existing processes occurring in the vegetation cover. It allows establishing a direction of their development in the system of natural factors of any territory. As a result, we have base for determination of age, site and role of current state of phytocenoses in successional systems. It is necessary for this to reveal the peculiarities of phytocenoses composition and formation due to climate dynamics and to determine a period of ecosystems homeostasis, especially for phytocenoses of contact natural conditions. Optimal values of phytocenoses diversity, like one of whole biosystems, depend on the amount of resource in the environment, on stability degree and on evolutional development of cenoses due to ecological, micro-evolutional and evolutional processes.
基金supported by the Natural Science Foundation of China under Grant No.51175423
文摘In order to implement the dynamic characteristic of a dual power-split transmission, a dynamic me- chanics model is built. Firstly, according to the method of theoretical analysis of the tooth contact analysis (TCA) and loaded tooth contact analysis (LTCA), the actual meshing process of each gear pairs is simulated, and the time-varying mesh stiffness excitations are obtained, which can improve the numerical precision. Second- ly, by using the lumped mass method, the bending-torsional coupling three dimensional dynamic model of the dual power-split transmission is established, and the identical dimensionless equations are deduced by elimina- ting the effect of rigid displacement and the method of dimensional normalization. Finally, by the method of the fourth order Runge-Kutta algorithm with variable step lengths, the responses of this system in a frequency domain and time domain are obtained, and the dynamic load change characteristics of each gear pairs are analyzed. The results show that the establishment, solution and analysis of the system dynamics model could provide a basis for the dynamic design, and have an important significance for the dynamic efficiency analysis and dynamic perform- ance optimization design of the dual power-split transmission.
文摘The effect of involute contact ratio on the torsional vibration behavior ofspur gear-pair is studied analytically through a mass-spring model. The tooth stiffness in model notonly has a relation with time, as many prior studies presented, but, more important, with involutecontact ratio (ICR) as well. The ICR embodies its impact on the spur gear's dynamic performancethrough changing the proportion of tooth stiffness when there are n+1 teeth in contact to toothstiffness when there are n teeth in contact. A couple of curves about impact of ICR on the gear'sdynamic performance are presented, and they clearly demonstrate that the model can accuratelydescribe the effects of ICR on dynamic transmission error. Finally, some conclusions useful toreduce vibration and noise of gear-pair are proposed.
基金The project supported by the National Natural Science Foundation of China(59935100)the Foundation of Ph.D Student Education of China(20020613001)the Foundation of the Author of National Excellent Doctoral Dissertation of China(2000048,2002048)
文摘The effect of rail corrugation on the vertical dynamics of railway vehicle coupled with a curved track is investigated in detail with a numerical method when a wheelset is steadily curving. In the calculation of rail corrugation we consider the combination of Kalkers rolling contact theory modified, a model of material loss on rail running surface, and a dynamics model of railway vehicle coupled with a curved track. In the establishment of the dynamic model, for simplicity, one fourth of the freight car without lateral motions, namely a wheelset and the equivalent one fourth freight car body above it, is considered. The Euler beam is used to model the rails and the track structure under the rails is replaced with equivalent springs, dampers and mass bodies. The numerical results show the great influence of the rail corrugation on the vibration of the parts of the vehicle and the track, and the some characters of rail corrugation in development.
基金supported in part by the National Natural Science Foundation of China (Grants 11290150 and 11290151)
文摘Soft machine refers to a kind of mechanical system made of soft materials to complete sophisticated missions, such as handling a fragile object and crawling along a narrow tunnel corner, under low cost control and actuation. Hence, soft machines have raised great challenges to computational dynamics. In this review article, recent studies of the authors on the dynamic modeling, numerical simulation, and experimental validation of soft machines are summarized in the framework of multibody system dynamics. The dynamic modeling approaches are presented first for the geometric nonlinearities of coupled overall motions and large deformations of a soft component, the physical nonlinearities of a soft component made of hyperelastic or elastoplastic materials, and the frictional contacts/impacts of soft components, respectively. Then the computation approach is outlined for the dynamic simulation of soft machines governed by a set of differential-algebraic equations of very high dimensions, with an emphasis on the efficient computations of the nonlinear elastic force vector of finite elements. The validations of the proposed approaches are given via three case studies, including the locomotion of a soft quadrupedal robot, the spinning deployment of a solar sail of a spacecraft, and the deployment of a mesh reflector of a satellite antenna, as well as the corresponding experimental studies. Finally, some remarks are made for future studies.
基金This research was done as part of TEKES-funded PanFlow project and as part of a project OPTIMI funded by the Academy of Finland (grant number 117587) in Micro- and Nanosystems Research Group, Tampere University of Technology, Finland.
文摘When characterizing flows in miniaturized channels, the determination of the dynamic contact angle is important. By measuring the dynamic contact angle, the flow properties of the flowing liquid and the effect of material properties on the flow can be characterized. A machine vision based system to measure the contact angle of front or rear menisci of a moving liquid plug is described in this article. In this research, transparent flow channels fabricated on thermoplastic polymer and sealed with an adhesive tape are used. The transparency of the channels enables image based monitoring and measurement of flow variables, including the dynamic contact angle. It is shown that the dynamic angle can be measured from a liquid flow in a channel using the image based measurement system. An image processing algorithm has been developed in a MATLAB environment. Images are taken using a CCD camera and the channels are illuminated using a custom made ring light. Two fitting methods, a circle and two parabolas, are experimented and the results are compared in the measurement of the dynamic contact angles.
基金Special Fund of Marine Commonweal Industry under contact Nos 201105016 and 201205007supported by National Marine Environment Forecasting Centrethe National Natural Science Foundation of China under contact No.41176012
文摘Considering the discontinuous characteristics of sea ice on various scales,a modified discrete element model(DEM) for sea ice dynamics is developed based on the granular material rheology.In this modified DEM,a soft sea ice particle element is introduced as a self-adjustive particle size function.Each ice particle can be treated as an assembly of ice floes,with its concentration and thickness changing to variable sizes under the conservation of mass.In this model,the contact forces among ice particles are calculated using a viscous-elastic-plastic model,while the maximum shear forces are described with the Mohr-Coulomb friction law.With this modified DEM,the ice flow dynamics is simulated under the drags of wind and current in a channel of various widths.The thicknesses,concentrations and velocities of ice particles are obtained,and then reasonable dynamic process is analyzed.The sea ice dynamic process is also simulated in a vortex wind field.Taking the influence of thermodynamics into account,this modified DEM will be improved in the future work.
基金Supported by China Postdoctoral Science Foundation(Grant No.2017M621458)National Science and Technology Support Plan Projects(Grant No.2015BAA08B02)National Natural Science Foundation of China(Grant No.11632011),National Natural Science Foundation of China(Grant No.11372183)
文摘Fractal theory provides scale?independent asperity contact loads and assumes variable curvature radii in the contact analyses of rough surfaces, the current research for which mainly focuses on the mechanism study. The present study introduces the fractal theory into the dynamic research of gas face seals under face?contacting conditions. Structure?Function method is adopted to handle the surface profiles of typical carbon?graphite rings, proving the fractal con?tact model can be used in the field of gas face seals. Using a numerical model established for the dynamic analyses of a spiral groove gas face seal with a flexibly mounted stator, a comparison of dynamic performance between the Majumdar?Bhushan(MB) fractal model and the Chang?Etsion?Bogy(CEB) statistical model is performed. The result shows that the two approaches induce differences in terms of the occurrence and the level of face contact. Although the approach distinctions in film thickness and leakage rate can be tiny, the distinctions in contact mechanism and end face damage are obvious. An investigation of fractal parameters D and G shows that a proper D(nearly 1.5) and a small G are helpful in raising the proportion of elastic deformation to weaken the adhesive wear in the sealing dynamic performance. The proposed research provides a fractal approach to design gas face seals.
基金Sponsored by Program for Changjiang Scholars and Innovative Research Team in University(Grant No.IRT0520)
文摘The existence of clearance in the joints of mechanisms system is inevitable.The movements of the real mechanism are deflection from the ideal mechanism due to the clearances and the motion accuracy is decreased.The effects of the hinge clearance on the crank and rocker mechanism system are studied.The system dynamics equation with clearance is presented.The contact dynamics model is established using the nonlinear equivalent spring-damp model and the friction effect is considered by using Coulomb friction model.Then the models are incorporated into ADAMS,and based on the model,large numbers numeric simulations are made.The regularity of contact forces in clearance are studied in detail.And the effects of clearance size,clearance friction on the mechanism dynamics characteristic are analyzed.The simulation results can predict the effects of clearance on the mechanism dynamics characteristic preferably.
基金supported by the National Natural Science Foundation of China(No.10602053)the Research Grants from Institute of Crustal Dynamics(Nos.ZDJ2007-2 and ZDJ2007-28).
文摘Dynamic contact stiffness at the interface between a vibrating rigid sphere and a semi-infinite transversely isotropic viscoelastic solid is investigated. An oscillating force superimposed onto a static compressive force in the vertical direction excites the vibration of a rigid sphere, which causes variable contact radius and contact pressure distribution in the contact region. The assumption of a sufficiently small oscillating force yields a dynamic contact-pressure distribution of a constant contact radius, which gives dynamic contact stiffness at the interface between the rigid sphere and the semi-infinite solid. Numerical calculations show the influence of vibration frequency of the sphere, and elastic constants of the transversely isotropic solid on dynamic contact stiffness, which benefits quantitative evaluation of elastic constants and orientation of single hexagonal grains by resonance-frequency shifts of the oscillator in resonance ultrasound microscopy.
基金Sponsored by the National Basic Research Program of China (Grant No. 2013CB632305)the National Natural Science Foundation of China (GrantNo. 51275125)
文摘In order to describe the performance of thin wall bearing on rotor system more accurate,the simplified model of bearing local stiffness was proposed. The load distribution and local contact deformation in angular contact ball bearings were calculated using quasi dynamic calculation method. Based on the relationship of local load to contact deformation,the calculation model of local bearing stiffness was subsequently built to get radial and axial components of local stiffness. Effects of external loads on the local bearing stiffness were analyzed. The results showed that local stiffness in bearings is symmetric to the axis of radial load,and its value has a maximum on the symmetry axis along the radial load direction. External radial and axial load have different effects on local bearing stiffness.