Recently,exploration breakthroughs have been made in the Lower Cretaceous sandstone reservoirs in the Doseo Basin,but the identification of reservoir fluid property is difficult due to variable reservoir lithology,com...Recently,exploration breakthroughs have been made in the Lower Cretaceous sandstone reservoirs in the Doseo Basin,but the identification of reservoir fluid property is difficult due to variable reservoir lithology,complex oil-water contact within and faint responses of the oil zone,which causes the lower accuracy of reservoir fluid property identification with conventional mudlogging and wirelogging techniques.Applying the geochemical logging,fluorescent logging,mud logging and cutting logging technology,in combination with formation test data,this paper distinguishes the crude oil types,analyzes the logging response characteristics of oil zone after water washing,and establishes the interpretation charts and parameter standards for reservoir fluid properties.The crude oil can be divided into two types,namely viscous-heavy and thin-light,based on total hydrocarbon content and component concentration tested by mud logging,features of pyrolysis gas chromatogram and fluorescence spectroscopy.The general characteristics of oil layers experienced water washing include the decrease of total hydrocarbon content and component concentration from mud logging,the decrease of S1 and PS values from geochemical logging,the decrease of hydrocarbon abundance and absence of some light components in pyrolysis gas chromatogram,and the decrease of fluorescence area and intensity from fluorescence logging.According to crude oil types,the cross plots of S1 versus peak-baseline ratio,and the cross plots of rock wettability versus fluorescence area ratio are drawn and used to interpret reservoir fluid property.Meanwhile,the standards of reservoir fluid parameter are established combining with the parameters of PS and the parameters in above charts,and comprehensive multiparameter correlation in both vertical and horizontal ways is also performed to interpret reservoir fluid property.The application in the Doseo Basin achieved great success,improving interpretation ability of fluid property in the reservoir with complex oil-water contact,and also provided technical reference for the efficient exploration and development of similar reservoirs.展开更多
The dynamic characteristics of a high speed angular-contact ceramic ball bearing are studied and compared with that of the steel ball bearing. According to rolling bearing analysis theory, the bearing dynamic equation...The dynamic characteristics of a high speed angular-contact ceramic ball bearing are studied and compared with that of the steel ball bearing. According to rolling bearing analysis theory, the bearing dynamic equations are established and are solved based on Hook-Jeeves's optimization theory on the computer. The results show that the bearing dynamic characteristics mainly depend on the rotational speed and ball material property at high speed. The bearing stiffness, initially decreases obviously and then increases with the increase of rotational speed. The ceramic ball bearing gains significant advantages over the steel ball bearing in high speed applications, such as lower contact stress, smaller deformation, less altering amount of contact angle, decreasing extent of variation of axial and radial stiffness and higher performance stability.展开更多
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 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.展开更多
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.展开更多
To investigate the cage stability of high-speed oil-lubricated angular contact ball bearings, a dynamic model of cages is developed on the basis of Gupta’s and Meeks’ work. The model can simulate the cage motion und...To investigate the cage stability of high-speed oil-lubricated angular contact ball bearings, a dynamic model of cages is developed on the basis of Gupta’s and Meeks’ work. The model can simulate the cage motion under oil lubrication with all six degrees of freedom. Particularly, the model introduces oil-film damping and hysteresis damping, and deals with the collision contact as imperfect elastic contact. In addition, the effects of inner ring rotational speed, the ratio of pocket clearance to guiding clearance and applied load on the cage stability are investigated by simulating the cage motion with the model. The results can provide a theoretical basis for the design of ball bearing parameters.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
Porous materials can be found in a variety of geophysical and engineering applications.The existence of thermal contact resistance at the interface between bilayered saturated porous strata would result in a significa...Porous materials can be found in a variety of geophysical and engineering applications.The existence of thermal contact resistance at the interface between bilayered saturated porous strata would result in a significant temperature difference at the interface.An attempt is made to study the thermo-hydro-mechanical coupling dynamic response of bilayered saturated porous strata with thermal contact resistance and elastic wave impedance.The corresponding analytical solutions for the dynamic response of bilayered saturated porous strata under a harmonic thermal load are derived by the operator decomposition method,and their rationality is verified by comparing them with existing solutions.The influences of thermal contact resistance,thermal conductivity ratio,and porosity ratio on the dynamic response of bilayered saturated porous strata are systematically investigated.Outcomes disclose that with the increase of thermal contact resistance,the displacement,pore water pressure and stress decrease gradually,and the temperature jump at the interface between two saturated porous strata increases.展开更多
The distribution characteristics of the oil-water contact are the basis for the reservoir exploration and development and reserves evaluation. The reservoir with a tilted oil-water contact has a unique formation mecha...The distribution characteristics of the oil-water contact are the basis for the reservoir exploration and development and reserves evaluation. The reservoir with a tilted oil-water contact has a unique formation mechanism, and the understanding of its distribution and formation mechanism will directly affect the evaluations for the reservoir type, well deployment, selection of well pattern and type, determination of test section, and reserves evaluation. Based on the analysis of reservoir characteristics, petrophysical properties and geological structure in 40 reservoirs worldwide with tilted oil-water contacts, the progress of the research on the formation mechanisms of titled oil-water contacts is summarized in terms of the hydrodynamic conditions, reservoir heterogeneity, neotectonic movement and oil-gas exploitation. According to the formation mechanism of tilted oil-water contacts and the needs of exploration research, different aspects of research methods are summarized and classified, such as the calculation of equipotential surfaces for oil and water in the formation, analysis of formation pressure and analysis of reservoir physical properties and so on. Based upon statistical analysis, it is suggested that the degree of the inclination of the oil-water contact be divided based on the dip of oil-water contact(DipTOWC). The tilted oil-water contact is divided into three categories: large dip(DipTOWC≥55 m/km), medium dip(4 m/km≤DipTOWC55 m/km), and small dip(DipTOWC4 m/km). The classification and evaluation method can be combined with structure amplitude and reservoir property. The formation mechanism of domestic and international reservoirs with tilted oil-water contacts are summarized in this paper, which have important significance in guiding the exploration and development of the oilfield with tilted oil-water contacts, reserves evaluation, and well deployment.展开更多
In this study,interface shapes of horizontal oil–water two-phase flow are predicted by using Young-Laplace equation model and minimum energy model.Meanwhile,the interface shapes of horizontal oil–water twophase flow...In this study,interface shapes of horizontal oil–water two-phase flow are predicted by using Young-Laplace equation model and minimum energy model.Meanwhile,the interface shapes of horizontal oil–water twophase flow in a 20 mm inner diameter pipe are measured by a novel conductance parallel-wire array probe(CPAP).It is found that,for flow conditions with low water holdup,there is a large deviation between the model-predicted interface shape and the experimentally measured one.Since the variation of pipe wetting characteristics in the process of fluid flow can lead to the changes of the contact angle between the fluid and the pipe wall,the models mentioned above are modified by considering dynamic contact angle.The results indicate that the interface shapes predicted by the modified models present a good consistence with the ones measured by CPAP.展开更多
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.展开更多
A dynamics formula was established for the flexible cage of high-speed angular contact ball bearing. A modified Craig-Bampton component mode synthetic method was used to establish the formula with regard to the flexib...A dynamics formula was established for the flexible cage of high-speed angular contact ball bearing. A modified Craig-Bampton component mode synthetic method was used to establish the formula with regard to the flexibility of cage and based on a dynamic analysis of angular contact ball bearing,and a rigid-flexible multi-body dynamic analysis program was developed using ADAMS,which is verified by a computation example of Gupta. The results show that it's not likely to keep the rotation smoothness of cage when the ratio of pocket clearance to guiding clearance and the ratio of radial load to axial load become too large or too small. By comparison,the flexible cage runs more smoothly than the rigid cage.展开更多
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.展开更多
The DRM (dynamic contact resistance measurement) in high voltage circuit breakers is a manner of evaluating the internal ageing condition of the chamber. DRM is similar to static contact resistance measurement testi...The DRM (dynamic contact resistance measurement) in high voltage circuit breakers is a manner of evaluating the internal ageing condition of the chamber. DRM is similar to static contact resistance measurement testing, but instead of measuring a single value when the breaker contacts are closed (static value), the ohmic resistance is measured at various contact positions, from the beginning of the contact opening until a complete separation of the contacts. The relationship between the contact resistances of the new circuit breaker and the ageing circuit breaker in operation provides subsidy for the evaluation of both the main and arcing contact conditions. This research aims to analyze the correlation between the various levels of degradation of the contacts and the configuration of the DRM curve. This work considers curve samples from new acceleration tests. breaker chamber contacts and different levels of degradation by展开更多
基金funded by a project entitled exploration field evaluation and target optimization of key basins in Chad and Niger(No.2019D-4308)initiated by the scientific research and technology development project of china national petroleum corporation.
文摘Recently,exploration breakthroughs have been made in the Lower Cretaceous sandstone reservoirs in the Doseo Basin,but the identification of reservoir fluid property is difficult due to variable reservoir lithology,complex oil-water contact within and faint responses of the oil zone,which causes the lower accuracy of reservoir fluid property identification with conventional mudlogging and wirelogging techniques.Applying the geochemical logging,fluorescent logging,mud logging and cutting logging technology,in combination with formation test data,this paper distinguishes the crude oil types,analyzes the logging response characteristics of oil zone after water washing,and establishes the interpretation charts and parameter standards for reservoir fluid properties.The crude oil can be divided into two types,namely viscous-heavy and thin-light,based on total hydrocarbon content and component concentration tested by mud logging,features of pyrolysis gas chromatogram and fluorescence spectroscopy.The general characteristics of oil layers experienced water washing include the decrease of total hydrocarbon content and component concentration from mud logging,the decrease of S1 and PS values from geochemical logging,the decrease of hydrocarbon abundance and absence of some light components in pyrolysis gas chromatogram,and the decrease of fluorescence area and intensity from fluorescence logging.According to crude oil types,the cross plots of S1 versus peak-baseline ratio,and the cross plots of rock wettability versus fluorescence area ratio are drawn and used to interpret reservoir fluid property.Meanwhile,the standards of reservoir fluid parameter are established combining with the parameters of PS and the parameters in above charts,and comprehensive multiparameter correlation in both vertical and horizontal ways is also performed to interpret reservoir fluid property.The application in the Doseo Basin achieved great success,improving interpretation ability of fluid property in the reservoir with complex oil-water contact,and also provided technical reference for the efficient exploration and development of similar reservoirs.
基金TheNaturalScienceFoundationofJiangsuProvince (No .BK2 0 0 2 0 5 9) theKeyTechnologiesR&DProgramofJiangsuProvinceDuringthe 10thFive YearPlanPeriod (No .BE2 0 0 3 0 71) .
文摘The dynamic characteristics of a high speed angular-contact ceramic ball bearing are studied and compared with that of the steel ball bearing. According to rolling bearing analysis theory, the bearing dynamic equations are established and are solved based on Hook-Jeeves's optimization theory on the computer. The results show that the bearing dynamic characteristics mainly depend on the rotational speed and ball material property at high speed. The bearing stiffness, initially decreases obviously and then increases with the increase of rotational speed. The ceramic ball bearing gains significant advantages over the steel ball bearing in high speed applications, such as lower contact stress, smaller deformation, less altering amount of contact angle, decreasing extent of variation of axial and radial stiffness and higher performance stability.
文摘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.
基金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.
基金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.
基金Supported by National Key Technology Research and Development Program of China during the 11th Five-Year Plan Period (No. JPPT-115-189)National Natural Science Foundation of China (No. 50975033)
文摘To investigate the cage stability of high-speed oil-lubricated angular contact ball bearings, a dynamic model of cages is developed on the basis of Gupta’s and Meeks’ work. The model can simulate the cage motion under oil lubrication with all six degrees of freedom. Particularly, the model introduces oil-film damping and hysteresis damping, and deals with the collision contact as imperfect elastic contact. In addition, the effects of inner ring rotational speed, the ratio of pocket clearance to guiding clearance and applied load on the cage stability are investigated by simulating the cage motion with the model. The results can provide a theoretical basis for the design of ball bearing parameters.
基金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.
基金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.
基金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(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.
文摘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.
基金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.
基金Projects(52108347,52178371)supported by the National Natural Science Foundation of ChinaProject(LQ22E080010)supported by the Exploring Youth Project of Zhejiang Natural Science Foundation,China。
文摘Porous materials can be found in a variety of geophysical and engineering applications.The existence of thermal contact resistance at the interface between bilayered saturated porous strata would result in a significant temperature difference at the interface.An attempt is made to study the thermo-hydro-mechanical coupling dynamic response of bilayered saturated porous strata with thermal contact resistance and elastic wave impedance.The corresponding analytical solutions for the dynamic response of bilayered saturated porous strata under a harmonic thermal load are derived by the operator decomposition method,and their rationality is verified by comparing them with existing solutions.The influences of thermal contact resistance,thermal conductivity ratio,and porosity ratio on the dynamic response of bilayered saturated porous strata are systematically investigated.Outcomes disclose that with the increase of thermal contact resistance,the displacement,pore water pressure and stress decrease gradually,and the temperature jump at the interface between two saturated porous strata increases.
文摘The distribution characteristics of the oil-water contact are the basis for the reservoir exploration and development and reserves evaluation. The reservoir with a tilted oil-water contact has a unique formation mechanism, and the understanding of its distribution and formation mechanism will directly affect the evaluations for the reservoir type, well deployment, selection of well pattern and type, determination of test section, and reserves evaluation. Based on the analysis of reservoir characteristics, petrophysical properties and geological structure in 40 reservoirs worldwide with tilted oil-water contacts, the progress of the research on the formation mechanisms of titled oil-water contacts is summarized in terms of the hydrodynamic conditions, reservoir heterogeneity, neotectonic movement and oil-gas exploitation. According to the formation mechanism of tilted oil-water contacts and the needs of exploration research, different aspects of research methods are summarized and classified, such as the calculation of equipotential surfaces for oil and water in the formation, analysis of formation pressure and analysis of reservoir physical properties and so on. Based upon statistical analysis, it is suggested that the degree of the inclination of the oil-water contact be divided based on the dip of oil-water contact(DipTOWC). The tilted oil-water contact is divided into three categories: large dip(DipTOWC≥55 m/km), medium dip(4 m/km≤DipTOWC55 m/km), and small dip(DipTOWC4 m/km). The classification and evaluation method can be combined with structure amplitude and reservoir property. The formation mechanism of domestic and international reservoirs with tilted oil-water contacts are summarized in this paper, which have important significance in guiding the exploration and development of the oilfield with tilted oil-water contacts, reserves evaluation, and well deployment.
基金supported by the National Natural Science Foundation of China(Grant Nos.41974139,41504104,11572220,51527805)Natural Science Foundation of Tianjin,China(19JCYBJC18400)。
文摘In this study,interface shapes of horizontal oil–water two-phase flow are predicted by using Young-Laplace equation model and minimum energy model.Meanwhile,the interface shapes of horizontal oil–water twophase flow in a 20 mm inner diameter pipe are measured by a novel conductance parallel-wire array probe(CPAP).It is found that,for flow conditions with low water holdup,there is a large deviation between the model-predicted interface shape and the experimentally measured one.Since the variation of pipe wetting characteristics in the process of fluid flow can lead to the changes of the contact angle between the fluid and the pipe wall,the models mentioned above are modified by considering dynamic contact angle.The results indicate that the interface shapes predicted by the modified models present a good consistence with the ones measured by CPAP.
基金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.
文摘A dynamics formula was established for the flexible cage of high-speed angular contact ball bearing. A modified Craig-Bampton component mode synthetic method was used to establish the formula with regard to the flexibility of cage and based on a dynamic analysis of angular contact ball bearing,and a rigid-flexible multi-body dynamic analysis program was developed using ADAMS,which is verified by a computation example of Gupta. The results show that it's not likely to keep the rotation smoothness of cage when the ratio of pocket clearance to guiding clearance and the ratio of radial load to axial load become too large or too small. By comparison,the flexible cage runs more smoothly than the rigid cage.
基金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.
文摘The DRM (dynamic contact resistance measurement) in high voltage circuit breakers is a manner of evaluating the internal ageing condition of the chamber. DRM is similar to static contact resistance measurement testing, but instead of measuring a single value when the breaker contacts are closed (static value), the ohmic resistance is measured at various contact positions, from the beginning of the contact opening until a complete separation of the contacts. The relationship between the contact resistances of the new circuit breaker and the ageing circuit breaker in operation provides subsidy for the evaluation of both the main and arcing contact conditions. This research aims to analyze the correlation between the various levels of degradation of the contacts and the configuration of the DRM curve. This work considers curve samples from new acceleration tests. breaker chamber contacts and different levels of degradation by