Thermal-contact capacity of one-dimensional attractive Gaudin-Yang model

Tan's contact C is an important quantity measuring the two-body correlations at short distances in a dilute system.Here we make use of the technique of exactly solved models to study the thermal-contact capacity K_(T),i.e.,the derivative of C with respect to temperature in the attractive Gaudin-Yang model.It is found that K_(T) is useful in identifying the low temperature phase diagram,and using the obtained analytical expression of K_(T),we study its critical behavior and the scaling law.Especially,we show K_(T) versus temperature and thus the non-monotonic tendency of C in a tiny interval,for both spin-balanced and imbalanced phases.Such a phenomenon is merely observed in multi-component systems such as SU(2)Fermi gases and spinor bosons,indicating the crossover from the Tomonaga-Luttinger liquid to the spin-coherent liquid.

Calculation of high-speed abrasion forces for ceramic coatings based on the Hertzian contact model

Aero engine seal coatings can effectively improve the air tightness of aircraft engines and increase fuel efficiency.However,due to the frictional forces between the blades and the coating,the coating often flakes off,resulting in damage to the blades and causing eco-nomic losses.Therefore,it is necessary to analyze the friction between the blades and the coating.In this paper,three ceramic-based high-temperature seal coatings with different polyphenylene ester contents and a pure Yttria-stabilised zirconia coating were prepared by atmo-spheric plasma spraying(APS).The hardness and modulus of elasticity of the coated surfaces were obtained by hardness and modulus of elasticity tests,and the coatings were subjected to high-speed touch abrasion tests.The Hertzian contact model was used to calculate the maximum normal contact load on the coating during the process.The test values were compared with the theoretical values and it was found that the calculated values were always greater than the test values.In order to make the Hertzian contact model more accurate in calculating the touching and abrasion forces,the contact coefficients in the Hertzian contact model were optimized.Substituting the optimized coeffi-cients into the Hertzian contact model,the results show that the calculated results after optimizing the coefficients are much closer to the test values,with deviations from the test values ranging from 1%to 38%.

Wheel-rail contact model for railway vehicle-structure interaction applications:development and validation

An enhancement in the wheel-rail contact model used in a nonlinear vehicle-structure interaction(VSI)methodology for railway applications is presented,in which the detection of the contact points between wheel and rail in the concave region of the thread-flange transition is implemented in a simplified way.After presenting the enhanced formulation,the model is validated with two numerical applications(namely,the Manchester Benchmarks and a hunting stability problem of a sus-pended wheelset),and one experimental test performed in a test rig from the Railway Technical Research Institute(RTRI)in Japan.Given its finite element(FE)nature,and contrary to most of the vehicle multibody dynamic commercial software that cannot account for the infrastructure flexibility,the proposed VSI model can be easily used in the study of train-bridge systems with any degree of complexity.The validation presented in this work proves the accuracy of the proposed model,making it a suitable tool for dealing with different railway dynamic applications,such as the study of bridge dynamics,train running safety under different scenarios(namely,earthquakes and crosswinds,among others),and passenger riding comfort.

Contact Angle Prediction Model for Underwater Oleophobic Surfaces Based on Multifractal Theory

Traditional microstructure scale parameters have difficulty describing the structure and distribution of a roughmaterial’s surface morphology comprehensively and quantitatively. This study constructs hydrophilic and underwateroleophobic surfaces based on polyvinylidene fluoride (PVDF) using a chemical modification method, and the fractaldimension and multifractal spectrum are used to quantitatively characterize the microscopic morphology. A new contactangle prediction model for underwater oleophobic surfaces is established. The results show that the fractal dimension ofthe PVDF surface first increases and then decreases with the reaction time. The uniformity characterized by the multifractalspectrum was generally consistent with scanning electron microscope observations. The contact angle of water droplets onthe PVDF surface is negatively correlated with the fractal dimension, and oil droplets in water are positively correlated.When the fractal dimension is 2.0975, the new contact angle prediction model has higher prediction accuracy. Themaximum and minimum relative deviations of the contact angle between the theoretical and measured data are 18.20%and 0.72%, respectively. For water ring transportation, the larger the fractal dimension and spectral width of the materialsurface, the smaller the absolute value of the spectral difference, the stronger the hydrophilic and oleophobic properties, andthe better the water ring transportation stability.

Numerical Stability and Accuracy of Contact Angle Schemes in Pseudopotential Lattice Boltzmann Model for Simulating Static Wetting and Dynamic Wetting

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.

Heterogeneities of grain boundary contact for simulation of laboratoryscale mechanical behavior of granitic rocks

From a practical point of view,grain structure heterogeneities are key parameters that control the rock response and still remains a challenge to incorporate in a quantitative manner.One of the less discussed topics in the context of the grain-based model(GBM)in the particle flow code(PFC)is the contact heterogeneities and the appropriate contact model to mimic the grain boundary behavior.Generally,the smooth joint(SJ)model and linear parallel bond(LPB)model are used to simulate the grain boundary behavior.However,the literature does not document the suitability of different models for specific problems.Another challenge in implementing GBM in PFC is that only a single bonding parameter is used at the grain boundaries.The aim of this study is to investigate the responses of a laboratory-scale specimen with SJ and LPB models,considering grain boundary heterogeneous and homogeneous contact parameters.Uniaxial and biaxial compression tests are performed to calibrate the response of Creighton granite.The stressestrain curves,volumetric dilation,inter-crack(crack in the grain boundary),and intra-crack(crack within the grain)development,and failure patterns associated with different contact models are examined.It was found that both the SJ and LPB models can reproduce the pre-peak behavior observed for a granitic rock type.However,the LPB model is unable to reproduce the post-peak behavior.Due to the large interlocking effect originating from the balls in contact and the ball size in the LPB model,local dilation is induced at the grain boundaries.This overestimates the volumetric dilation and residual shear strength.The LPB model tends to result in discontinuous inter-cracks and stress localization in the rock specimen,resulting in fine fragments at the rock surface during failure.

Experiment and constitutive modelling of creep deformation in the frozen silt-concrete interface

To ensure the long-term safety and stability of bridge pile foundations in permafrost regions,it is necessary to investigate the rheological effects on the pile tip and pile side bearing capacities.The creep characteristics of the pile-frozen soil interface are critical for determining the long-term stability of permafrost pile foundations.This study utilized a self-developed large stress-controlled shear apparatus to investigate the shear creep characteristics of the frozen silt-concrete interface,and examined the influence of freezing temperatures(−1,−2,and−5°C),contact surface roughness(0,0.60,0.75,and 1.15 mm),normal stress(50,100,and 150 kPa),and shear stress on the creep characteristics of the contact surface.By incorporating the contact surface’s creep behavior and development trends,we established a creep constitutive model for the frozen silt-concrete interface based on the Nishihara model,introducing nonlinear elements and a damage factor.The results revealed significant creep effects on the frozen silt-concrete interface under constant load,with creep displacement at approximately 2-15 times the instantaneous displacement and a failure creep displacement ranging from 6 to 8 mm.Under different experimental conditions,the creep characteristics of the frozen silt-concrete interface varied.A larger roughness,lower freezing temperatures,and higher normal stresses resulted in a longer sample attenuation creep time,a lower steady-state creep rate,higher long-term creep strength,and stronger creep stability.Building upon the Nishihara model,we considered the influence of shear stress and time on the viscoelastic viscosity coefficient and introduced a damage factor to the viscoplasticity.The improved model effectively described the entire creep process of the frozen silt-concrete interface.The results provide theoretical support for the interaction between pile and soil in permafrost regions.

Fractal Prediction Model of Thermal Contact Conductance of Rough Surfaces

The thermal contact conductance problem is an important issue in studying the heat transfer of engineering surfaces, which has been widely studied since last few decades, and for predicting which many theoretical models have been established. However, the models which have been existed are lack of objectivity due to that they are mostly studied based on the statistical methodology characterization for rough surfaces and simple partition for the deformation formats of contact asperity. In this paper, a fractal prediction model is developed for the thermal contact conductance between two rough surfaces based on the rough surface being described by three-dimensional Weierstrass and Mandelbrot fractal function and assuming that there are three kinds of asperity deformation modes： elastic, elastoplastic and fully plastic. Influences of contact load and contact area as well as fractal parameters and material properties on the thermal contact conductance are investigated by using the presented model. The investigation results show that the thermal contact conductance increases with the increasing of the contact load and contact area. The larger the fractal dimension, or the smaller the fractal roughness, the larger the thermal contact conductance is. The thermal contact conductance increases with decreasing the ratio of Young＇s elastic modulus to the microhardness. The results obtained indicate that the proposed model can effectively predict the thermal contact conductance at the interface, which provide certain reference to the further study on the issue of heat transfer between contact surfaces.

Elastoplastic Contact Mechanics Model of Rough Surface Based on Fractal Theory

Because the result of the MB fractal model contradicts with the classical contact mechanics, a revised elastoplastic contact model of a single asperity is developed based on fractal theory. The critical areas of a single asperity are scale dependent, with an increase in the contact load and contact area, a transition from elastic, elastoplastic to full plastic deformation takes place in this order. In considering the size distribution function, analytic expression between the total contact load and the real contact area on the contact surface is obtained. The elastic, elastoplastic and full plastic contact load are obtained by the critical elastic contact area of the biggest asperity and maximun contact area of a single asperity. The results show that a rough surface is firstly in elastic deformation. As the load increases, elastoplastic or full plastic deformation takes place. For constant characteristic length scale G, the slope of load-area relation is proportional to fractal dimension D. For constant fractal dimension D, the slope of load-area relation is inversely proportional to G. For constant D and G, the slope of load-area relation is inversely proportional to property of the material ~b, namely with the same load, the material of rough surface is softer, and the total contact area is larger. The contact mechanics model provides a foundation for study of the friction, wear and seal performance of rough surfaces.

Modeling and Calculation of Impact Friction Caused by Corner Contact in Gear Transmission

Corner contact in gear pair causes vibration and noise,which has attracted many attentions.However,teeth errors and deformation make it difficulty to determine the point situated at corner contact and study the mechanism of teeth impact friction in the current researches.Based on the mechanism of corner contact,the process of corner contact is divided into two stages of impact and scratch,and the calculation model including gear equivalent error-combined deformation is established along the line of action.According to the distributive law,gear equivalent error is synthesized by base pitch error,normal backlash and tooth profile modification on the line of action.The combined tooth compliance of the first point lying in corner contact before the normal path is inversed along the line of action,on basis of the theory of engagement and the curve of tooth synthetic complianceload-history.Combined secondarily the equivalent error with the combined deflection,the position standard of the point situated at corner contact is probed.Then the impact positions and forces,from the beginning to the end during corner contact before the normal path,are calculated accurately.Due to the above results,the lash model during corner contact is founded,and the impact force and frictional coefficient are quantified.A numerical example is performed and the averaged impact friction coefficient based on the presented calculation method is validated.This research obtains the results which could be referenced to understand the complex mechanism of teeth impact friction and quantitative calculation of the friction force and coefficient,and to gear exact design for tribology.

Contact Analysis and Modeling of Standing Wave Linear Ultrasonic Motor

A contact model for describing the contact mechanics between the stator and slider of the standing wave linear ultrasonic motor was presented. The proposed model starts from the assumption that the vibration characteristics of the stator is not affected by the contact process. A modified friction models was used to analyze the contact problems. Firstly, the dynamic normal contact force, interface friction force, and steady-state characteristics were analyzed. Secondly, the influences of the contact layer material, the dynamic characteristics of the stator, and the pre-load on motor performance were simulated. Finally, to validate the contact model, a linear ultrasonic motor based on in-plane modes was used as an example. The corresponding results show that a set of simulation of motor performances based on the proposed contact mechanism is in good agreement with experimental results. This model is helpful to understanding the operation principle of the standing wave linear motor and thus contributes to the design of these types of motor.

Influence of particle contact models on soil response of poorly graded sand during cavity expansion in discrete element simulation

The discrete element method(DEM) has been extensively adopted to investigate many complex geotechnical related problems due to its capability to incorporate the discontinuous nature of granular materials. In particular, when simulating large deformations or distortion of soil(e.g. cavity expansion),DEM can be very effective as other numerical solutions may experience convergence problems. Cavity expansion theory has widespread applications in geotechnical engineering, particularly to the problems concerning in situ testing, pile installation and so forth. In addition, the behaviour of geomaterials in a macro-level is utterly determined by microscopic properties, highlighting the importance of contact models. Despite the fact that there are numerous contact models proposed to mimic the realistic behaviour of granular materials, there are lack of studies on the effects of these contact models on the soil response.Hence, in this study, a series of three-dimensional numerical simulations with different contact constitutive models was conducted to simulate the response of sandy soils during cylindrical cavity expansion. In this numerical investigation, three contact models, i.e. linear contact model, rolling resistance contact model,and Hertz contact model, are considered. It should be noted that the former two models are linear based models, providing linearly elastic and frictional plasticity behaviours, whereas the latter one consists of nonlinear formulation based on an approximation of the theory of Mindlin and Deresiewicz. To examine the effects of these contact models, several cylindrical cavities were created and expanded gradually from an initial radius of 0.055 m to a final radius of 0.1 m. The numerical predictions confirm that the calibrated contact models produced similar results regarding the variations of cavity pressure, radial stress, deviatoric stress, volumetric strain, as well as the soil radial displacement. However, the linear contact model may result in inaccurate predictions when highly angular soil particles are involved. In addition, considering the excessive soil displacement induced by the pile installation(i.e. cavity expansion), a minimum distance of11 a(a is the cavity radius) is recommend for practicing engineers to avoid the potential damages to the existing piles and adjacent structures.

CONVERGENCE RATE OF SOLUTIONS TO STRONG CONTACT DISCONTINUITY FOR THE ONE-DIMENSIONAL COMPRESSIBLE RADIATION HYDRODYNAMICS MODEL

This paper is concerned with a singular limit for the one-dimensional compress- ible radiation hydrodynamics model. The singular limit we consider corresponds to the physical problem of letting the Bouguer number infinite while keeping the Boltzmann number constant. In the case when the corresponding Euler system admits a contact discontinuity wave, Wang and Xie （2011） [12] recently verified this singular limit and proved that the solution of the compressible radiation hydrodynamics model converges to the strong contact 1 discontinuity wave in the L∞-norm away from the discontinuity line at a rate of ε1/4, as the reciprocal of the Bouguer number tends to zero. In this paper, Wang and Xie＇s convergence rate is improved to ε7/8 by introducing a new a priori assumption and some refined energy estimates. Moreover, it is shown that the radiation flux q tends to zero in the L∞-norm away from the discontinuity line, at a convergence rate as the reciprocal of the Bouguer number tends to zero.

Determination of the Normal Contact Stiffness and Integration Time Step for the Finite Element Modeling of Bristle-Surface Interaction

In finite element modeling of impact,it is necessary to define appropriate values of the normal contact stiffness,Kn,and the Integration Time Step(ITS).Because impacts are usually of very short duration,very small ITSs are required.Moreover,the selection of a suitable value of Kn is a critical issue,as the impact behavior depends dramatically on this parameter.In this work,a number of experimental tests and finite element analyses have been performed in order to obtain an appropriate value of Kn for the interaction between a bristle of a gutter brush for road sweeping and a concrete surface.Furthermore,a suitable ITS is determined.The experiments consist of releasing a steel bristle that is placed vertically at a certain distance from a concrete surface and tracking the impact.Similarly,in the finite element analyses,a beam is modeled in free fall and impacting a surface;contact and target elements are attached to the beam and the surface,respectively.The results of the experiments and the modeling are integrated through the principle of conservation of energy,the principle of linear impulse and momentum,and Newton’s second law.The results demonstrate that,for the case studied,Kn and the impact time tend to be independent of the velocity just before impact and that Kn has a very large variation,as concrete is a composite material with a rough surface.Also,the ratio between the largest height of the bristle after impact and the initial height tends to be constant.

Contact Model of Revolute Joint with Clearance Based on Fractal Theory

The contact sti ness of a mechanical bonding surface is an important parameter in determining the normal and radial contact force. To improve the calculation accuracy of the contact force model, the surface roughness of the bonding surface and the energy loss that necessarily occurs during the impact process should be considered com?prehensively. To study the normal contact force of a revolute joint with clearance more accurately in the case of dry friction, a nonlinear sti ness coe cient model considering the surface roughness was established based on fractal theory, which considers the elastic, elastic?plastic, and plastic deformations of the asperities of the contact surface during the contact process. On this basis, a modified nonlinear spring damping model was established based on the Lankarani–Nikravesh contact force model. The laws influencing the surface roughness, recovery coe cient, initial velocity, and clearance size on the impact force were revealed, and were compared with the Lankarani–Nikravesh model and a hybrid model using MATLAB. The maximum impact force was obtained using a modified contact force model under di erent initial velocities, di erent clearances, and di erent degrees of surface roughness, and the calculated results were then compared with the experiment results. This study indicates that the modified model can be used more widely than other models, and is suitable for both large and small clearances. In particular, the modified model is more accurate when calculating the contact force of a revolute joint with a small clearance.

Adhesive contact:from atomistic model to continuum model

Two types of Lennard-Jones potential are widely used in modeling adhesive contacts. However, the relationships between the parameters of the two types of Lennard-Jones potential are not well defined. This paper employs a self- consistent method to derive the Lennard-Jones surface force law from the interatomic Lennard-Jones potential with emphasis on the relationships between the parameters. The ei^ect of using correct parameters in the adhesion models is demonstrated in single sphere-flat contact via continuum models and an atomistic model. Furthermore, the adhesion hysteresis behaviour is investigated, and the S-shaped force-distance relation is revealed by the atomistic model. It shows that the adhesion hysteresis loop is generated by the jump-to-contact and jump-off-contact, which are illustrated by the S-shaped force-distance curve.

Modeling Methods and Test Verification of Root Insert Contact Interface for Wind Turbine Blade

Two modeling methods of the root insert for wind turbine blade are presented,i.e.,the local mesh optimization method(LMOM)and the global modeling method(GMM).Based on the optimized mesh of the local model for the metal contact interface,LMOM is proposed to analyze the load path and stress distribution characteristics,while GMM is used to calculate and analyze the stress distribution characteristics of the resin layer established between the bushing and composite layers of root insert.To validate the GMM,a tension test is carried out.The result successfully shows that the shear strain expresses a similar strain distribution tendency with the GMM′s results.

Mathematic model and tooth contact analysis of a new spiral non-circular bevel gear

A novel spiral non-circular bevel gear that could be applied to variable-speed driving in intersecting axes was proposed by combining the design principles of non-circular bevel gears and the manufacturing principles of face-milling spiral bevel gears.Unlike straight non-circular bevel gears,spiral non-circular bevel gears have numerous advantages,such as a high contact ratio,high intensity,good dynamic performance,and an adjustable contact region.In addition,while manufacturing straight non-circular bevel gears is difficult,spiral non-circular bevel gears can be efficiently and precisely fabricated with a 6-axis bevel gear cutting machine.First,the generating principles of spiral non-circular bevel gears were introduced.Next,a mathematical model,including a generating tooth profile,tooth spiral,pressure angle,and generated tooth profile for this gear type was established.Then the precision of the model was verified by a tooth contact analysis using FEA,and the contact patterns and stress distributions of the spiral non-circular bevel gears were investigated.

Aggregate gradation design of asphalt mixture with stone-to-stone contact based on fuller's model

Based on Fuller＇s model, a gradation design is presented for the load-carrying capacity of an asphalt mixture with stone-to-stone contact. The coarse aggregate gradation, fine aggregate and filler gradation are respectively designed for their different effects on an asphalt mixture with stone-to-stone contact. The equation of coarse aggregate gradation is deduced from the Fuller＇s model. Fine aggregate gradation is calculated by using Fuller＇s equation. The asphalt mixture with the 16 mm maximum size of aggregate is verified. When the coefficient of coarse aggregate gradation is 0.35 and fine aggregate gradation is 0.25, the volume parameters of asphalt mixture meet the technical specifications. Moreover, the high temperature stability and crack resistance at low temperature are all well.

Analytical modeling of complex contact behavior between rock mass and lining structure

Based on the nondestructive test data of operating railway tunnels in China, this paper summarizes the basic characteristics of the complex contact behavior between the rock mass and lining structure. The contact modes are classified into dense contact, local non-contact, and loose contact. Subsequently, the corresponding mechanical model for each contact mode is developed according to its mechanical characteristics using the complex variable method. In the proposed mechanical model, a special algorithm is introduced to detect whether the local non-contact zone is re-contacted. Besides, a novel conformal mapping method is also proposed to accurately calculate the mechanical response of the concrete lining. Finally, the accuracy of the proposed method is verified by comparing it with the finite element method(FEM). Several parameter investigations are conducted to analyze the effects of different contact modes on the rock-lining interaction. The results show that:(i) the height of the local noncontact area does not have a significant effect on the contact stress distribution if no re-contact occurs;(ii) backfill grouting can reduce the local stress concentration caused by poor contact modes;and(iii) reducing the friction coefficient of the interface can lead to a more uniform distribution of internal forces in the concrete lining.