Fractal model of thermal elasto-plastic contact of rough surfaces

Without considering the influence of heat,existing fractal contact models are not applicable to analyze the contacts when the temperature changes.For this problem,the normal load model and the normal stiffness model of thermal elasto-plastic contact of rough surfaces are developed respectively in this paper.The proposed model is based on the normal contact mechanics model of fractal theory of anisotropic and thermal elasto-plastic contact theory which can be used to characterize the rough surface thermodynamic properties.Then the validity of the model is verified.Finally,the influence of main parameters on the total normal load and the whole normal stiffness of thermal elasto-plastic contact at the interface is analyzed by contact simulation.The results show that the total normal load of thermal elasto-plastic contact increases with the increases of temperature.The whole normal stiffness of thermal elasto-plastic contact increases with increasing coefficient of linear expansion,scale factor,temperature difference or fractal dimension,but decreases with increasing fractal roughness.This model expands basic theory and applications of traditional models,and can be used to calculate and analyze the contacts when the temperature changes.

Research on synthetic aperture radar imaging technology of one-dimensional layered rough surfaces

A quick and exact imaging method for one-dimensional layered rough surfaces is proposed in this paper to study the scattering characteristics of a layered medium that exists widely in nature.The boundary integral equations of layered rough surfaces are solved by using the propagation-inside-layer expansion combined with the forward and backward spectral acceleration method（PILE＋FB-SA）,and the back scattering data are obtained.Then,a conventional synthetic aperture radar（SAR） imaging procedure called back projection method is used to generate a two-dimensional（2D） image of the layered rough surfaces.Combined with the relative dielectric permittivity of realistic soil,the random rough surfaces with Gauss spectrum are used to simulate the layered medium with rough interfaces.Since the back scattering data are computed by using the fast numerical method,this method can be used to study layered rough surfaces with any parameter,which has a great application value in the detection and remote sensing areas.

Modeling bidirectional reflection distribution function of microscale random rough surfaces

The radiative properties of three different materials surfaces with one-dimensional microscale random roughness were obtained with the finite difference time domain method(FDTD) and near-to-far-field transformation.The surface height conforms to the Gaussian probability density function distribution.Various computational modeling issues that affect the accuracy of the predicted properties were discussed.The results show that,for perfect electric conductor(PEC) surfaces,as the surface roughness increases,the magnitude of the spike reduces and eventually the spike disappears,and also as the ratio of root mean square roughness to the surface correlation distance increases,the retroreflection becomes evident.The predicted values of FDTD solutions are in good agreement with the ray tracing and integral equation solutions.The overall trend of bidirectional reflection distribution function(BRDF) of PEC surfaces and silicon surfaces is the same,but the silicon's is much less than the former's.The BRDF difference from two polarization modes for the gold surfaces is little for smaller wavelength,but it is much larger for the longer wavelength and the FDTD simulation results agree well with the measured data.In terms of PEC surfaces,as the incident angle increases,the reflectivity becomes more specular.

Scattering from Rough Surfaces Using the"Shooting and Bouncing Rays" Technique

The “shooting and bouncing rays” (SBR) technique is used to analyze the electromagnetic scattering characters of ocean rough surfaces varying with time. Some numerical results are presented and compared with the method of moments, and some factors, such as the incident angle, polarization and frequency are investigated which influence on electromagnetic scattering characters of ocean rough surfaces.

Carpet Cloak Design for Rough Surfaces

Conventional carpet cloak structures have been utilized to conceal the objects located on a planar perfect electric conductor surface. We systematically investigate hiding arbitrarily shaped objects on a rough surface, as a more general and practical scenario. In addition, the required cloak is designed considering different boundary conditions for the surface beneath the object, despite the previous studies. To achieve an invisibility cloak, taking advantage of linear coordinate transformation, a simple homogeneous material is obtained to realize the cloak structure, facilitating the fabrication processes. Numerical simulations validate the performance of the proposed cloaking method. Therefore, the proposed structure is capable of cloaking in more general and complicated scenarios.

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.

Modeling of rough surfaces with given roughness parameters

Modeling of rough surfaces with given roughness parameters is studied,where surfaces with Gaussian height distribution and exponential auto-correlation function(ACF) are concerned.A large number of micro topography samples are randomly generated first using the rough surface simulation method with FFT.Then roughness parameters of the simulated roughness profiles are calculated according to parameter definition,and the relationship between roughness parameters and statistical distribution parameters is investigated.The effects of high-pass filtering with different cut-off lengths on the relationship are analyzed.Subsequently,computing formulae of roughness parameters based on standard deviation and correlation length are constructed with mathematical regression method.The constructed formulae are tested with measured data of actual topographies,and the influences of auto-correlation variations at different lag lengths on the change of roughness parameter are discussed.The constructed computing formulae provide an approach to active modeling of rough surfaces with given roughness parameters.

Wear Resistance of 3Cr2W8V Rough Surfaces

Three types of rough surface were processed by laser irradiation on the 3Cr2W8V material hot-work die steel surface. The wear experiments with smooth surface and rough surface samples were repeated on the pin-tray wear machine. According to the wear results, we studied the regularity of wear resistance of different rough surface samples. The results indicated that bionic rough surface can improve the wear resistance of the material and the wear resistance can be increased 1-2 times, compared with the smooth surface. Also, the wear resistance of the rough surface was affected by laser current and duration of impulse. The bigger the laser current or the impulse duration, the better is the wear resistance. When the distance between the same kind of units which are distributed on the surfaces is changed, the wear resistance changes. The wear resistance of a bionic rough surface on which the grid units were distributed at spacing of 1 mm was the best. And we designed the wear models.

Preload-responsive adhesion of microfibre arrays to rough surfaces

Adhesion of bio-inspired microfibre arrays to a rough surface is studied theoretically. The array consists of vertical elastic rods fixed on a rigid backing layer, and the surface is modeled by rigid steps with a normally distributed height. Analytical expressions are obtained for the adhesion forces in both the approach and retraction processes. It is shown that, with the increasing preload, the pull-off force increases at first and then attains a plateau value. The results agree with the previous experiments and are expected helpful in adhesion control of the array in practical applications.

Study of MPI based on parallel MOM on PC clusters for EM-beam scattering by 2-D PEC rough surfaces

This paper firstly applies the finite impulse response filter （FIR） theory combined with the fast Fourier transform （FFT） method to generate two-dimensional Gaussian rough surface. Using the electric field integral equation （EFIE）, it introduces the method of moment （MOM） with RWG vector basis function and Galerkin＇s method to investigate the electromagnetic beam scattering by a two-dimensional PEC Gaussian rough surface on personal computer （PC） clusters. The details of the parallel conjugate gradient method （CGM） for solving the matrix equation are also presented and the numerical simulations are obtained through the message passing interface （MPI） platform on the PC clusters. It finds significantly that the parallel MOM supplies a novel technique for solving a two-dimensional rough surface electromagnetic-scattering problem. The influences of the root-mean-square height, the correlation length and the polarization on the beam scattering characteristics by two-dimensional PEC Gaussian rough surfaces are finally discussed.

Variational Approach to Scattering by Inhomogeneous Layers Above Rough Surfaces

In this paper, we study, via variational methods, the problem of scattering of time harmonic acoustic waves by unbounded inhomogeneous layers above a sound soft rough surface. We first propose a variational formulation and exploit it as a theoretical tool to prove the well-posedness of this problem when the media is non-absorbing for arbitrary wave number and obtain an estimate about the solution, which exhibit explicitly dependence of bound on the wave number and on the geometry of the domain. Then, based on the non-absorbing results, we show that the variational problem remains uniquely solvable when the layer is absorbing by means of a priori estimate of the solution. Finally, we consider the finite element approximation of the problem and give an error estimate.

VIBRATION OF CANTILEVERS WITH ROUGH SURFACES

We study the effect of surface roughness on the resonance frequency of micro-cantilever sensors. The analysis demonstrates that surface roughness can enhance, decrease or even annul the effect of surface stress on the resonance frequency, depending on the surface inclination angle and the Poisson ratio of the coating film on the cantilever.

Adhesive Contact of Elastic Solids with Self-Affine Fractal Rough Surfaces

The elastic adhesive contact of self-affine fractal rough surfaces against a rigid flat is simulated using the finite element method. An array of nonlinear springs, of which the force-separation law obeys the Lennard–Jones potential, is introduced to account for the interfacial adhesion. For fractal rough surfaces, the interfacial interaction is generally attractive for large mean gaps, but turns repulsive as the gap continuously shrinks. The interfacial interactions at the turning point corresponding to the spontaneous contact are shown for various surfaces. For relatively smooth surfaces, the probability density distributions of repulsion and attraction are nearly symmetric. However, for rougher surfaces, the simulation results suggest a uniform distribution for attraction but a monotonously decreasing distribution with a long tail for repulsion. The pull-off force rises with increasing ratio of the work of adhesion to the equilibrium distance, whereas decreases for solids with a higher elastic modulus and a larger surface roughness. The current study will be helpful for understanding the adhesion of various types of rough solids.

A new 3D plastoelastohydrodynamic lubrication model for rough surfaces

Plastoelastohydrodynamic lubrication of rough surfaces(R-PEHL)is a cutting-edge area of research in interface fluid-structure coupling analysis.The existing R-PEHL model calculates the elastic-plastic deformation of rough surface by the Love equation in a semi-infinite space smooth surface,which deviates from the actual surface.Therefore,it is an innovative work to study the exact solution of elastic-plastic deformation of rough surface and its influence on the solution results of R-PEHL model.In this paper,a new contact calculation model of plastoelastohydrodynamic lubrication(PEHL)with three-dimensional(3D)rough surface is proposed by integrating numerical method of EHL and finite element method.The new model eliminates an original error introduced by the assumption of semi-infinite space in contact calculation,providing wide applicability and high accuracy.Under the given rough surfaces and working conditions,the study reveals that:(1)the oil film pressure calculated by the new model is lower than that of the smooth surface in semi-infinite space by 200–800 MPa;(2)the Mises stress of the new model is 2.5%–26.6%higher than that of the smooth surface in semi-infinite space;(3)compared with the semi-infinite space assumption,the rough surface plastic deformation of the new model is increased by 71%–173%,and the local plastic deformation singularity may appear under the semi-infinite space assumption;(4)the plastic deformation caused by the first contact cycle on the rough surface of the new model accounts for 66.7%–92.9%of the total plastic deformation,and the plastic deformation of the semi-infinite space accounts for 50%–83.3%.This study resolves the contradiction between the smooth surface assumption and the rough surface in the existing R-PEHL model,establishing a solid logic foundation for the accurate solution of R-PEHL model.

THE NYSTROM METHOD FOR ELASTIC WAVE SCATTERINGBY UNBOUNDED ROUGH SURFACES

We consider a numerical algorithm for the two-dimensional time-harmonic elastic wave scattering by unbounded rough surfaces with Dirichlet boundary condition.A Nystr¨om method is proposed for the scattering problem based on the integral equation method.Convergence of the Nystr¨om method is established with convergence rate depending on the smoothness of the rough surfaces.In doing so,a crucial role is played by analyzing the singularities of the kernels of the relevant boundary integral operators.Numerical experiments are presented to demonstrate the effectiveness of the method.Mathematics subject classification:35P25,45P05.

Theoretical Modeling and Surface Roughness Prediction of Microtextured Surfaces in Ultrasonic Vibration-Assisted Milling

Textured surfaces with certain micro/nano structures have been proven to possess some advanced functions,such as reducing friction,improving wear and increasing wettability.Accurate prediction of micro/nano surface textures is of great significance for the design,fabrication and application of functional textured surfaces.In this paper,based on the kinematic analysis of cutter teeth,the discretization of ultrasonic machining process,transformation method of coordinate systems and the cubic spline data interpolation,an integrated theoretical model was established to characterize the distribution and geometric features of micro textures on the surfaces machined by different types of ultrasonic vibration-assisted milling(UVAM).Based on the theoretical model,the effect of key process parameters(vibration directions,vibration dimensions,cutting parameters and vibration parameters)on tool trajectories and microtextured surface morphology in UVAM is investigated.Besides,the effect of phase difference on the elliptical shape in 2D/3D ultrasonic elliptical vibration-assisted milling(UEVAM)was analyzed.Compared to conventional numerical models,the method of the cubic spline data interpolation is applied to the simulation of microtextured surface morphology in UVAM,which is more suitable for characterizing the morphological features of microtextured surfaces than traditional methods due to the presence of numerous micro textures.The prediction of surface roughness indicates that the magnitude of ultrasonic amplitude in z-direction should be strictly limited in 1D rotary UVAM,2D and 3D UEVAM due to the unfavorable effect of axial ultrasonic vibration on the surface quality.This study can provide theoretical guidance for the design and fabrication of microtextured surfaces in UVAM.

Shear sliding of rough-walled fracture surfaces under unloading normal stress

Through high-precision engraving,self-affine sandstone joint surfaces with various joint roughness coefficients(JRC=3.21e12.16)were replicated and the shear sliding tests under unloading normal stress were conducted regarding various initial normal stresses(1e7 MPa)and numbers of shearing cycles(1 e5).The peak shear stress of fractures decreased with shear cycles due to progressively smooth surface morphologies,while increased with both JRC and initial normal stress and could be verified using the nonlinear Barton-Bandis failure criterion.The joint friction angle of fractures exponentially increased by 62.22%e64.87%with JRC while decreased by 22.1%e24.85%with shearing cycles.After unloading normal stress,the sliding initiation time of fractures increased with both JRC and initial normal stress due to more tortuous fracture morphologies and enhanced shearing resistance capacity.The surface resistance index(SRI)of fractures decreased by 4.35%e32.02%with increasing shearing cycles due to a more significant reduction of sliding initiation shear stress than that for sliding initiation normal stress,but increased by a factor of 0.41e1.64 with JRC.After sliding initiation,the shear displacement of fractures showed an increase in power function.By defining a sliding rate threshold of 5105 m/s,transition from“quasi-static”to“dynamic”sliding of fractures was identified,and the increase of sliding acceleration steepened with JRC while slowed down with shearing cycles.The normal displacement experienced a slight increase before shear sliding due to deformation recovery as the unloading stress was unloaded,and then enhanced shear dilation after sliding initiation due to climbing effects of surface asperities.Dilation was positively related to the shear sliding velocity of fractures.Wear characteristics of the fracture surfaces after shearing failure were evaluated using binary calculation,indicating an increasing shear area ratio by 45.24%e91.02%with normal stress.

Slip boundary conditions for rough surfaces

The velocity slip and temperature jump for a two-dimensional rough plate under hypersonic conditions were analyzed using the Direct Simulation Monte Carlo(DSMC)method.Surface roughness was explicitly modeled by introducing various structures on the flat plate.The influences of relative roughness height,which involves the roughness height,roughness spacing,incoming velocity,and the degree of rarefaction,were analyzed and discussed.It is found that with the increase of the relative roughness height,the jump temperature increases,while the slip velocity decreases gradually.The effects of surface roughness on the slip coefficients can be attributed to the change of accommodation coefficients.A new slip model for rough surfaces was established in this paper,which accounts for the coupling effects of gas rarefaction and surface roughness,without the effort to model the surface roughness explicitly.The nitrogen flows in the microchannel,and flows over a blunt cone and an axisymmetric bi-conic body,were simulated under the modified and conventional slip boundary conditions,respectively.The numerical solutions were validated with experimental data.It can be safely concluded that compared with the traditional first-order slip boundary conditions,the modified slip model improves the accuracy of macroscopic properties,especially the heat transfer coefficient.

Numerical Reconstruction of Locally Rough Surfaces with a Newton Iterative Algorithm

In this paper,we propose a Newton iterative algorithm to numerically reconstruct a locally rough surface with Dirichlet and impedance boundary conditions by near-field measurements of acoustic waves.The algorithm relies on the Frechet differentiability analysis of the locally rough surface scattering problem,which is established by reducing the original model into an equivalent boundary value problem with compactly supported boundary data.With a slight modification,the algorithm can be also extended to reconstruct the local perturbation of a non-local rough surface.Finally,numerical results are presented to illustrate the effectiveness of the inversion algorithm with the multi-frequency data.

A Stochastic Gradient Descent Method for Computational Design of Random Rough Surfaces in Solar Cells

In this work,we develop a stochastic gradient descent method for the computational optimal design of random rough surfaces in thin-film solar cells.We formulate the design problems as random PDE-constrained optimization problems and seek the optimal statistical parameters for the random surfaces.The optimizations at fixed frequency as well as at multiple frequencies and multiple incident angles are investigated.To evaluate the gradient of the objective function,we derive the shape derivatives for the interfaces and apply the adjoint state method to perform the computation.The stochastic gradient descent method evaluates the gradient of the objective function only at a few samples for each iteration,which reduces the computational cost significantly.Various numerical experiments are conducted to illustrate the efficiency of the method and significant increases of the absorptance for the optimal random structures.We also examine the convergence of the stochastic gradient descent algorithm theoretically and prove that the numerical method is convergent under certain assumptions for the random interfaces.