A hybrid contact approach for modeling soil-structure interaction using the material point method

The grid-based multi-velocity field technique has become increasingly popular for simulating the Material Point Method(MPM)in contact problems.However,this traditional technique has some shortcomings,such as(1)early contact and contact penetration can occur when the contact conditions are unsuitable,and(2)the method is not available for contact problems involving rigid-nonrigid materials,which can cause numerical instability.This study presents a new hybrid contact approach for the MPM to address these limitations to simulate the soil and structure interactions.The approach combines the advantages of point-point and point-segment contacts to implement contact detection,satisfying the impenetrability condition and smoothing the corner contact problem.The proposed approach is first validated through a disk test on an inclined slope.Then,several typical cases,such as granular collapse,bearing capacity,and deformation of a flexible retaining wall,are simulated to demonstrate the robustness of the proposed approach compared with FEM or analytical solutions.Finally,the proposed method is used to simulate the impact of sand flow on a deformable structure.The results show that the proposed contact approach can well describe the phenomenon of soil-structure interaction problems.

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.

Gaps,challenges and possible solution for prediction of wheel-rail rolling contact fatigue crack initiation

The prediction of wheel/rail rolling contact fatigue(RCF)crack initiation during railway operations is an important task.Since RCF crack evolution is influenced by many factors,its prediction process is complex.This paper reviews the existing approaches to predict RCF crack initiation.The crack initiation region is predicted by the shakedown map.By combining the shakedown map with various initiation criteria and the critical plane method,the crack initiation life is calculated.The classification,methodologies,theories and applications of these approaches are included in this paper.The advantages and limitations of these methods are analyzed to provide recommendation for RCF crack initiation prediction.This review highlights that wheel/rail dynamic characteristic,complex working conditions,surface defects and wear all affect the RCF crack initiation.The optimal selection of criteria is essential in the crack initiation prediction.Based on the research gap regarding the challenging process of crack initiation prediction detailed in this review,a proposed prediction process of RCF crack initiation is proposed to achieve a more accurate result.

Research on the influence of flexible wheelset rotation effect on wheel rail contact force

Purpose-Under the high-speed operating conditions,the effects of wheelset elastic deformation on the wheel rail dynamic forces will become more notable compared to the low-speed condition.In order to meet different analysis requirements and selecting appropriate models to analyzing the wheel rail interaction,it is crucial to understand the influence of wheelset flexibility on the wheel-rail dynamics under different speeds and track excitations condition.Design/methodology/approach-The wheel rail contact points solving method and vehicle dynamics equations considering wheelset flexibility in the trajectory body coordinate system were investigated in this paper.As for the wheel-rail contact forces,which is a particular force element in vehicle multibody system,a method for calculating the Jacobian matrix of the wheel-rail contact force is proposed to better couple the wheel-rail contact force calculation with the vehicle dynamics response calculation.Based on the flexible wheelset modeling approach in this paper,two vehicle dynamic models considering the wheelset as both elastic and rigid bodies are established,two kinds of track excitations,namely normal measured track irregularities and short-wave irregularities are used,wheel-rail geometric contact characteristic and wheel-rail contact forces in both time and frequency domains are compared with the two models in order to study the influence of flexible wheelset rotation effect on wheel rail contact force.Findings-Under normal track irregularity excitations,the amplitudes of vertical,longitudinal and lateral forces computed by the flexible wheelset model are smaller than those of the rigid wheelset model,and the virtual penetration and equivalent contact patch are also slightly smaller.For the flexible wheelset model,the wheel rail longitudinal and lateral creepages will also decrease.The higher the vehicle speed,the larger the differences in wheel-rail forces computed by the flexible and rigid wheelset model.Under track short-wave irregularity excitations,the vertical force amplitude computed by the flexible wheelset is also smaller than that of the rigid wheelset.However,unlike the excitation case of measured track irregularity,under short-wave excitations,for the speed within the range of 200 to 350 km/h,the difference in the amplitude of the vertical force between the flexible and rigid wheelset models gradually decreases as the speed increase.This is partly due to the contribution of wheelset's elastic vibration under short-wave excitations.For low-frequency wheel-rail force analysis problems at speeds of 350 km/h and above,as well as high-frequency wheel-rail interaction analysis problems under various speed conditions,the flexible wheelset model will give results agrees better with the reality.Originality/value-This study provides reference for the modeling method of the flexible wheelset and the coupling method of wheel-rail contact force to the vehicle multibody dynamics system.Furthermore,by comparative research,the influence of wheelset flexibility and rotation on wheel-rail dynamic behavior are obtained,which is useful to the application scope of rigid and flexible wheelset models.

Influence of railway wheel tread damage on wheel-rail impact loads and the durability of wheelsets

Dynamic wheel-rail contact forces induced by a severe form of wheel tread damage have been measured by a wheel impact load detector during full-scale field tests at different vehicle speeds.Based on laser scanning,the measured three-dimensional damage geometry is employed in simulations of dynamic vehicle-track interaction to calibrate and verify a simulation model.The relation between the magnitude of the impact load and various operational parameters,such as vehicle speed,lateral position of wheel-rail contact,track stiffness and position of impact within a sleeper bay,is investigated.The calibrated model is later employed in simulations featuring other forms of tread damage;their effects on impact load and subsequent fatigue impact on bearings,wheel webs and subsurface initiated rolling contact fatigue of the wheel tread are assessed.The results quantify the effects of wheel tread defects and are valuable in a shift towards condition-based maintenance of running gear,and for general assessment of the severity of different types of railway wheel tread damage.

Existence and numerical approximation of a solution to frictional contact problem for electro-elastic materials

In this paper,a frictional contact problem between an electro-elastic body and an electrically conductive foundation is studied.The contact is modeled by normal compliance with finite penetration and a version of Coulomb’s law of dry friction in which the coefficient of friction depends on the slip.In addition,the effects of the electrical conductivity of the foundation are taken into account.This model leads to a coupled system of the quasi-variational inequality of the elliptic type for the displacement and the nonlinear variational equation for the electric potential.The existence of a weak solution is proved by using an abstract result for elliptic variational inequalities and a fixed point argument.Then,a finite element approximation of the problem is presented.Under some regularity conditions,an optimal order error estimate of the approximate solution is derived.Finally,a successive iteration technique is used to solve the problem numerically and a convergence result is established.

Multi-point Contact of the High-speed Vehicle-turnout System Dynamics

The wheel-rail contact problems, such as the number, location and the track of contact patches, are very important for optimizing the spatial structure of the rails and lowering the vehicle-turnout system dynamics. However, the above problems are not well solved currently because of having the difficulties in how to determine the multi-contact, to preciously present the changeable profiles of the rails and to establish an accurate spatial turnout system dynamics model. Based on a high-speed vehicle-turnout coupled model in which the track is modeled as flexible with rails and sleepers represented by beams, the line tracing extreme point method is introduced to investigate the wheel-rail multiple contact conditions and the key sections of the blade rail, longer nose rail, shorter rail in the switch and nose rail area are discretized to represent the varying profiles of rails in the turnout. The dynamic interaction between the vehicle and turnout is simulated for cases of the vehicle divergently passing the turnout and the multi-point contact is obtained. The tracks of the contact patches on the top of the rails are presented and the wheel-rail impact forces are offered in comparison with the contact patches transference on the rails. The numerical simulation results indicate that the length of two-point contact occurrence of a worn wheel profile and rails is longer than that of the new wheel profile and rails; The two-point contact definitely occurs in the switch and crossing area. Generally, three-point contact doesn’t occur for the new rail profile, which is testified by the wheel-rails interpolation distance and the first order derivative function of the tracing line extreme points. The presented research is not only helpful to optimize the structure of the turnout, but also useful to lower the dynamics of the high speed vehicle-turnout system.

NUMERICAL SIMULATION OF TWO-POINT CONTACT BETWEEN WHEEL AND RAIL

The elastic-plastic contact problem with rolling friction of wheel-rail is solved using the FE parametric quadratic programming method. Thus, the complex elastic-plastic contact problem can be calculated with high accuracy and efficiency, while the Hertz＇s hypothesis and the elastic semi-space assumption are avoided. Based on the ‘one-point＇ contact calculation of wheel-rail, the computational model of ‘two-point＇ contact are established and calculated when the wheel flange is close to the rail. In the case of ‘two-point＇ contact, the changing laws of wheelrail contact are introduced and contact forces in various load cases are carefully analyzed. The main reason of wheel flange wear and rail side wear is found. Lubrication computational model of the wheel flange is constructed. Comparing with the result without lubrication, the contact force between wheel flange and rail decreases, which is beneficial for reducing the wear of wheel-rail.

Transition from tunneling regime to local point contact realized on Ba0.6K0.4Fe2As2 surface

Using scanning tunneling spectroscopy, we studied the transition from tunneling regime to local point contact on the iron-based superconductor Ba0.6K0.4Fe2As2. By gradually reducing the junction resistance, a series of spectra were obtained with the characteristics evolving from single-particle tunneling into Andreev reflection. The spectra can be well fitted to the modified Blonder–Tinkham–Klapwijk（BTK） model and exhibit significant changes of both spectral broadening and orbital selection due to the formation of point contact. The spatial resolution of the point contact was estimated to be several nanometers, providing a unique way to study the inhomogeneity of unconventional superconductors on such a scale.

Effect of Surface Texturing on Lubricating Condition under Point Contact Using Numerical Analysis

The contact fatigue life of machine elements is affected by pitting, wear and so on, under heavier loading conditions. Increasing the fatigue life requires mainly the improvements of lubricating condition, operating condition and materials. In order to improve the lubricating condition, it is necessary to investigate the relation of the microscopic surface texturing and the contact modes of machine elements. In this paper, thus, the pressure and oil film thickness of the contact between sphere and the plate with 5 kinds surface texturing were calculated using a commercial software based on Reynolds equation. There was sufficient evidence to suggest that the dimple shape was the optimum texturing to increase the lubricating condition.

An integrated approach for the optimization of wheel-rail contact force measurement systems

A comprehension of railway dynamic behavior implies the measure of wheel-rail contact forces which are affected by disturbances and errors that are often difficult to be quantified. In this study, a benchmark test case is proposed, and a bogie with a layout used on some European locomotives such as SIEMENS El90 is studied. In this layout, an additional shaft on which brake disks are installed is used to transmit the braking torque to the wheelset through a single-stage gearbox. Using a mixed approach based on finite element techniques and statistical considerations, it is possible to evaluate an optimal layout for strain gauge positioning and to optimize the measurement system to diminish the effects of noise and disturbance. We also conducted preliminary evaluations on the precision and frequency response of the proposed system.

Contact Force Distribution and Static Load-carrying Capacity of Large Size Double Row Four-point Contact Ball Bearing

Clearance not only affects the startup torque,rotation precision and stiffness of bearing,but also affects the load distribution,load-carrying capacity and life of bearing.A computational model in which the clearance of bearing is first included is presented for determining the contact force distribution and static load-carrying capacity of a double row four-point contact ball bearing which is subjected to the combined radial,axial and overturning moment loadings.The relation between the negative axial clearance and the contact force distribution is analyzed.The static load-carrying capacity curves are established,and the effects of the changes in negative axial clearance,curvature radius coefficient of raceway groove and initial contact angle on the static load-carrying capacity are analyzed.The results show that,with the increase in the absolute value of negative clearance,the maximum contact load decreases first and then increases.The clearance values in the range of 0.2 mme0 mm have little effect on the static load-carrying capacity of bearing.With the increase in the curvature radius coefficient of raceway groove and the decrease in the initial contact angle,the static load capacity of bearing decreases.

Magnetic vortex gyration mediated by point-contact position

Micromagnetic simulation is employed to study the gyration motion of magnetic vortices in distinct permalloy nanodisks driven by a spin-polarized current. The critical current density for magnetic vortex gyration, eigenfrequency, trajectory, velocity and the time for a magnetic vortex to obtain the steady gyration are analyzed. Simulation results reveal that the magnetic vortices in larger and thinner nanodisks can achieve a lower-frequency gyration at a lower current density in a shorter time. However, the magnetic vortices in thicker nanodisks need a higher current density and longer time to attain steady gyration but with a higher eigenfrequency. We also find that the point-contact position exerts different influences on these parameters in different nanodisks, which contributes to the control of the magnetic vortex gyration. The conclusions of this paper can serve as a theoretical basis for designing nano-oscillators and microwave frequency modulators.

Distinction between critical current effects and intrinsic anomalies in the point-contact Andreev reflection spectra of unconventional superconductors

In this work,we discuss the origin of several anomalies present in the point-contact Andreev reflection spectra of（Li1-xFex）OHFeSe,LiTi2O4,and La2-xCexCuO4.While these features are similar to those stemming from intrinsic superconducting properties,such as Andreev reflection,electron-boson coupling,multigap superconductivity,d-wave and p-wave pairing symmetry,they cannot be accounted for by the modified Blonder–Tinkham–Klapwijk（BTK） model,but require to consider critical current effects arising from the junction geometry.Our results point to the importance of tracking the evolution of the dips and peaks in the differential conductance as a function of the bias voltage,in order to correctly deduce the properties of the superconducting state.

A solution to the hydrodynamic lubrication of a circular point contact sliding over a flat surface with cavitation

This letter presents an analytical solution to the hydrodynamic lubrication of a circular point contact sliding over a flat surface with cavitation.The solution is found by solving the Reynolds equation with Reynolds boundary condition for cavitation.The cavitation boundary is shown to be straight lines directed 108.4°against the sliding direction.The result is experimentally verified in the limit of large values of viscosity,sliding velocity and radius of a spherical ball.The solution raises questions about the coupling between cavitation and film rupture and can be used as an independent check on the validity of numerical solutions.

Simulation on derailment base on a new wheel-rail contact method

The simulation package for special research on derailment of high speed vehicle is established.The process of derailment is different from other behaviors of vehicle dynamics because of large lateral displacement of wheelsets.To get correct results,a new fast algorithm to computing contact force is adopted and the exact geometry analysis is necessary to judge derailment happened.Variation of contact condition and coefficient of friction with speeds are also considered into vehicle-track coupled model.The structure of the package is presented in detail.The results are particular emphasis on investigation influence of maximum track defect,critical vehicle speed and various contact condition on derailment.The simulation can also be used to define the most risk factor leading to derailment.

Modeling of Energy Processes in Wheel-Rail Contacts Operating under Influence of Periodic Discontinuous Forces

In this paper we present new numerical simulation approaches for determining the energy processes under periodic conditions caused by time-discontinuous forces in the wheel-rail contacts. The main advantage of the presented method is the total elimination of frequency analysis, which in effect introduces important simplifications in the identification of the effects in the contact. The second important feature is the fact that the method is based on the analysis of appropriate loops on the energy phase plane leading to an easy estimation of the rail strength through the evaluation of the loop’s area. That model based simulation in the applied dynamics relies on advanced methods for model setup, robust and efficient numerical solution techniques and powerful simulation tools for practical applications. Fundamental properties of contact displacements of the rail surface have been considered on the basis of the newly established method. The contact zone between railway wheels and the rail surfaces made of bulk materials is perceived as strong enough to resist the normal (vertical) forces introduced by heavy loads and the dynamic response induced by track and wheel irregularities. The analysis is carried out for a wheel running on an elastic rail rested on sleepers arranged on completely rigid foundation. The equations of displacement motion are established through the application of the Lagrange equations approach. The established model of the wheel-rail contact dynamics has been applied to that same roll plane but with taking into account a nonlinear characteristic of the sleeper with respect to the ground. Attention then is focused completely on the modeling of the energy absorbed by the rail. The applied method employs the energy state variables as time functions leading to determine the susceptibility of a given contact on the strength induced by the rail roll.

Shot Noise Suppression in a Quantum Point Contact with Short Channel Length

An experimental study on the current shot noise of a quantum point contact with short channel length is reported. The experimentally measured maximum energy level spacing between the ground and the first excited state of the device reached up to 7.5meV, probably due to the hard wall confinement by using shallow electron gas and sharp point contact geometry. The two-dimensionM non-equilibrium shot noise contour map shows noise suppression characteristics in a wide range of bias voltage. Fano factor analysis indicates spin-polarized transport through a short quantum point contact.

The 0.7-Anomaly in Quantum Point Contact;Many-Body or Single-Electron Effect?

Apart from usual quantization steps on the ballistic conductance of quasi-one-dimensional conductor, an additional plateau-like feature appears at a fraction of about 0.7 below the first conductance step in GaAs-based quantum point contacts (QPCs). Despite a tremendous amount of research on this anomalous feature, its origin remains still unclear. Here, a unique model of this anomaly is proposed relying on fundamental principles of quantum mechanics. It is noticed that just after opening a quasi-1D conducting channel in the QPC a single electron travels the channel at a time, and such electron can be—in principle—observed. The act of observation destroys superposition of spin states, in which the electron otherwise exists, and this suppresses their quantum interference. It is shown that then the QPC-conductance is reduced by a factor of 0.74. “Visibility” of electron is enhanced if the electron spends some time in the channel due to resonant transmission. Electron’s resonance can also explain an unusual temperature behavior of the anomaly as well as its recently discovered feature: oscillatory modulation as a function of the channel length and electrostatic potential. A recipe for experimental verification of the model is given.

Strength and failure characteristics of marble spheres subjected to paired point loads

Failure of irregular rock samples may provide implications in the rapid estimation of rock strength,which is imperative in rock engineering practice.In this work,analytical,experimental and numerical investigations were carried out to study the mechanical properties and failure characteristics of rock spheres under paired point loads.Analytical solutions indicted that with the increase in sample size(contact angle)and decrease in Poisson’s ratio,the uneven tensile stress in theta direction decreased.Then laboratory experiments were carried out to investigate the load characteristics and failure mode of spherical marble samples with different sizes subjected to a pair of diametral point loads.The discrete element method(DEM)was adopted to study the failure process of rock spheres.The effect of the sphere diameter on the point load contact angle was examined in terms of peak load,crushed zone distribution and energy dissipation.Experimental and numerical results showed that the samples primarily fail in tension,with crushed zones formed at both loading points.With increase in the sample size,the contact angle,crushed area and total work increase.As the specimen diameter increases from 30 mm to 50 mm,the peak load on the specimen increases from 3.6 kN to 8.8 kN,and the percentage of crushed zone(ratio of crushing zone to sample radius,d/r)increased from 0.191 to 0.262.The results of the study have implications for understanding the failure of irregular rock specimens under point loading conditions and their size effects.