Three-dimensional analytical elasticity solution for loaded functionally graded coated half-space

In the present paper, a three-dimensional problem of elasticity of normal and tangential loading of surface of the functionally graded coated half-space is considered. In case when Poisson's ratio is constant and the Young's modulus is a power or exponential function of the distance from the surface of the half-space, analytical solution using Fourier transform is obtained. Stress field due to Hertz contact pressure in an elliptical region are studied as a function of the parameter b/a (where a and b are axes of the contact ellipse) and coating thickness.     

Boundary film shear elastic modulus effect in hydrodynamic contact. Part II: influence of operational parameters

The present paper is the subsequent research of the first part (Theor Comput Fluid Dyn, 2009). It investigates the boundary film shear elastic modulus effect in a hydrodynamic contact in different operating conditions. The hydrodynamic contact is one-dimensional, composed of two parallel plane surfaces, which are respectively rough rigid with rectangular micro projections in profile periodically distributed on the surface and ideally smooth rigid. The whole contact consists of cavitated area and hydrodynamic area. The hydrodynamic area consists of many micro Raleigh bearings which are discontinuously and periodically distributed in the contact. The hydrodynamic contact in a micro Raleigh bearing consists of boundary film area and fluid film area which, respectively, occur in the outlet and inlet zones. In boundary film area, the film slips at the upper contact surface due to the limited shear stress capacity of the film–contact interface, while the film does not slip at the lower contact surface due to the shear stress capacity of the film–contact interface large enough. In boundary film area, the viscosity, density, and shear elastic modulus of the film are varied across the film thickness due to the film–contact interactions, and their effective values are used in modeling which depends on the film thickness. In fluid film area, the film does not slip at either of the contact surfaces, and the shear elastic modulus of the film is neglected. It is found from the simulation results that the boundary film shear elastic modulus influences are normally negligible on the mass flow through the contact, the carried load of the contact and the overall film thickness of the contact, and the boundary film shear elastic modulus would normally influence the local film thickness in an elastic contact when the local film thickness is on the film molecule diameter scale. It is also found that the boundary film shear elastic modulus effect has the tendency of being increased with the reduction of the width of a micro contact. It is increased with the reduction of the boundary film–contact interfacial shear strength or with the increase of the critical boundary film thickness, while it is strongest at certain values of the contact surface roughness, the width ratio of fluid film area to boundary film area, and the lubricant film shear elastic modulus.

---

     

Micromechanical study of elastic moduli of loose granular materials

In micromechanics of the elastic behaviour of granular materials, the macro-scale continuum elastic moduli are expressed in terms of micro-scale parameters, such as coordination number (the average number of contacts per particle) and interparticle contact stiffnesses in normal and tangential directions. It is well-known that mean-field theory gives inaccurate micromechanical predictions of the elastic moduli, especially for loose systems with low coordination number. Improved predictions of the moduli are obtained here for loose two-dimensional, isotropic assemblies. This is achieved by determining approximate displacement and rotation fields from the force and moment equilibrium conditions for small sub-assemblies of various sizes. It is assumed that the outer particles of these sub-assemblies move according to the mean field. From the particle displacement and rotation fields thus obtained, approximate elastic moduli are determined. The resulting predictions are compared with the true moduli, as determined from the discrete element method simulations for low coordination numbers and for various values of the tangential stiffness (at fixed value of the normal stiffness). Using this approach, accurate predictions of the moduli are obtained, especially when larger sub-assemblies are considered. As a step towards an analytical formulation of the present approach, it is investigated whether it is possible to replace the local contact stiffness matrices by a suitable average stiffness matrix. It is found that this generally leads to a deterioration of the accuracy of the predictions. Many micromechanical studies predict that the macroscopic bulk modulus is hardly influenced by the value of the tangential stiffness. It is shown here from the discrete element method simulations of hydrostatic compression that for loose systems, the bulk modulus strongly depends on the stiffness ratio for small stiffness ratios.     

Boundary film shear elastic modulus effect in hydrodynamic contact. Part I: theoretical analysis and typical solution

Boundary film shear elastic modulus effect is analyzed in a hydrodynamic contact. The contact is one-dimensional composed of two parallel plane surfaces, which are, respectively, rough rigid with rectangular micro projections in profile periodically distributed on the surface and ideally smooth rigid. The whole contact is consisted of cavitated area and hydrodynamic area. The hydrodynamic area consists of many micro Raleigh bearings which are discontinuously and periodically distributed in the contact. Analysis is thus carried out for a micro Raleigh bearing in this contact. The hydrodynamic contact in this micro Raleigh bearing consists of boundary film area and fluid film area which, respectively, occur in the outlet and inlet zones. In boundary film area, the film slips at the upper contact surface due to the limited shear stress capacity of the film–contact interface, while the film does not slip at the lower contact surface due to the shear stress capacity large enough at the film–contact interface. In boundary film area, the viscosity, density and shear elastic modulus of the film are varied across the film thickness due to the film–contact interactions, and their effective values are used in modeling, which depend on the film thickness. The analytical approach proposed by Zhang (J Mol Liq 128:60–64, 2006) and Zhang et al. (Int J Fluid Mech Res 30:542–557, 2003) is used for boundary film area. In fluid film area, the film does not slip at either of the contact surfaces, and the shear elastic modulus of the film is neglected. Conventional hydrodynamic analysis is used for fluid film area. The present paper presents the theoretical analysis and a typical solution. It is found that for the simulated case the boundary film shear elastic modulus effects on the mass flow through the contact, the overall film thickness of the contact and the carried load of the contact are negligible but the boundary film shear elastic modulus effect on the local film thickness of the contact may be significant when the boundary film thickness is on the 1 nm scale and the contact surfaces are elastic. In Part II will be presented detailed results showing boundary film shear elastic modulus effects in different operating conditions.

---

     

Effective elastic properties of porous materials: Homogenization schemes vs experimental data

In this paper, we focus on the prediction of elastic moduli of isotropic porous materials made of a solid matrix having a Poisson's ratio vm of 0.2. We derive simple analytical formulae for these effective moduli based on well-known Mean-Field Eshelby-based Homogenization schemes. For each scheme, we find that the normalized bulk, shear and Young's moduli are given by the same form depending only on the porosity p. The various predictions are then confronted with experimental results for the Young's modulus of expanded polystyrene (EPS) concrete. The latter can be seen as an idealized porous material since it is made of a bulk cement matrix, with Poisson's ratio 0.2, containing spherical mono dispersed EPS beads. The Differential method predictions are found to give a very good agreement with experimental results. Thus, we conclude that when vm=0.2, the normalized effective bulk, shear and Young's modulus of isotropic porous materials can be well predicted by the simple form (1 − p)² for a large range of porosity p ranging between 0 and 0.56.     

齿向修形对滤波减速器润滑性能的影响分析

综合考虑了滤波减速器齿向修形参数、真实齿面粗糙度和瞬态效应等因素,建立了轮齿混合润滑数学模型,数值计算了不同修形参数值对应不同啮合点的最大压力和中心膜厚,分析了齿面粗糙度和转速对润滑性能的影响.结果表明:修形参数r和Ry均存在一个优化范围,使得轮齿表面最大油膜压力显著降低,边缘效应弱化,而中心膜厚则随着r和Ry的增大而逐渐增大;未修形轮齿边缘油膜压力受粗糙度的影响而急剧增大,边缘效应更加显著,修形后轮齿的边缘效应得到了明显改善,因此,轮齿修形也因粗糙表面的存在而显得更加重要;随着转速逐渐降低,轮齿表面的平均油膜厚度逐渐变小,接触比逐渐增大,轮齿表面由弹流润滑逐渐转为混合润滑,最后演变为边界润滑.     

Determining the Elastic Modulus of Compliant Thin Films Supported on Substrates from Flat Punch Indentation Measurements

Instrumented indentation is a technique that can be used to measure the elastic properties of soft thin films supported on stiffer substrates, including polymer films, cellulosic sheets, and thin layers of biological materials. When measuring thin film properties using indentation, the effect of the substrate must be considered. Most existing models for determining the properties of thin films from indentation measurements were developed for metal and dielectric films bonded to semiconductor substrates and have been applied to systems with film-substrate modulus ratios between 0.1 and 10. In the present work, flat punch indentation of a thin film either bonded to or in contact with a substrate is examined using finite element modeling. A broad range of film-substrate modulus ratios from 0.0001 to 1 are investigated. As the substrate is effectively rigid compared to the film when the film-substrate modulus ratio is less than 0.0001, the results are also useful for understanding systems with lower film-substrate modulus ratios. The effects of the contact radius, film thickness, elastic properties, and friction between the film and the substrate on the measured stiffness were quantified using finite element modeling in order to understand how the elastic properties of the film can be extracted from indentation measurements. A semi-analytical model was developed to describe the finite element modeling results and facilitate the use of the results to analyze experimental measurements. The model was validated through analysis of indentation measurements of thin polyethylene sheets that were supported on substrates of various stiffness.     

The electro-mechanical response of highly compliant substrates and thin stiff films with periodic cracks

We present results that describe the mechanical response of highly compliant substrates coated with ultra-thin stiff films, with thickness and elastic moduli differences spanning four orders of magnitude. Dimensional analysis based on shear-lag models of cracked films is used to identify key parameters that control the effective elastic properties of the cracked multi-layer, crack opening displacements, and the steady-state energy release rate for channeling crack formation. Analytical forms that describe multi-layer response in terms of film properties and crack spacing are presented and corroborated with numerical models for linear elastic materials. A key result is that the energy release rate scales with 1/(1 − α), where α is one of the Dundurs’ parameters describing elastic mismatch. The results can also be used to evaluate the performance of electrostrictive actuators comprised of cracked blanket electrodes and elastomer dielectrics. In this scenario, an interesting result is that ultra-thin cracked films can continue to distribute charge, since crack openings may be small enough to allow breakdown in air at typical operating voltages.     

Useful zero friction simulations for assessing MBS codes Pascal's formula giving wheelsets frequency for zero wheel-rail friction

Present study provides a simple analytical formula, the “Klingel-like formula” or “Pascal's Formula” that can be used as a reference to test some results of existing railway codes and specifically those using rigid contact. It develops properly the 3D Newton-Euler equations governing the 6 degrees of freedom (DoF) of unsuspended loaded wheelsets in case of zero wheel-rail friction and constant conicity. Thus, by solving numerically these equations, we got pendulum like harmonic oscillations of which the calculated angular frequency is used for assessing the accuracy of the proposed formula so that it can in turn be used as a fast practical target for testing multi-body system (MBS) railway codes. Due to the harmonic property of these pendulum-like oscillations, the square ω² of their angular frequency can be made in the form of a ratio K/M where K depends on the wheelset geometry and load and M on its inertia. Information on K and M are useful to understand wheelsets behavior. The analytical formula is derived from the first order writing of full trigonometric Newton-Euler equations by setting zero elastic wheel-rail penetration and by assuming small displacements. Full trigonometric equations are numerically solved to assess that the formula provides ω² inside a 1% accuracy for usual wheelsets dimensions. By decreasing the conicity down to 1 × 10^?4 rad, the relative formula accuracy is under 3 × 10^?5. In order to test the formula reliability for rigid contact formulations, the stiffness of elastic contacts can be increased up to practical rigidity (Hertz stiffness × 1000).     

Initiation of fracture at the interface corner of bi-material joints

Within the context of linear elasticity, a stress singularity of the form Hr^λ−1 may exist at the interface corner of a bi-material joint, where r is the radial distance from the corner, H is the stress intensity factor and λ−1 is the order of the singularity. Recent experimental results in the literature support the use of a critical value of the intensity factor H=H_c as a fracture initiation criterion at the interface corner. In this paper, we examine the validity and limitations of this criterion for predicting the onset of fracture in a butt joint consisting of a thin layer of an elastic-plastic adhesive layer sandwiched between two elastic adherends. The evolution of plastic deformation at the corner is determined theoretically and by the finite element method, and the solution is compared with the extent of the elastic singular field. It is shown that H_c is a valid fracture parameter if h>B(H_c/σ_Y)^1/(1−λ) where the non-dimensional constant B=100 for β=0 and B=13 for β=α/4. Here, h is the thickness of the adhesive layer, σ_Y is the uniaxial yield stress of the bulk adhesive and (α,β) are Dundurs’ parameters (Dundurs, J., J. Appl. Mech. 36 (1969) 650). Experimental results for aluminium/epoxy/aluminium and brass/solder/brass sandwiched joints are used to assess the role of plastic deformation on the validity of the failure criterion.     

边界膜剪切弹性模量对于微接触性能的影响

建立接触模型,理论分析了微接触中边界膜剪切弹性模量对于接触性能的影响. 接触区由两平行平面形成,属一维接触. 上接触表面为粗糙表面,具有矩形微凸体. 下接触表面为光滑平面. 两接触表面均处理成刚性表面. 微接触区中充满流体. 它分成两个子区,在微接触的出口区由于极小的接触间隙充满边界膜,在微接触的入口区由于接触间隙较大充满流体膜. 边界膜和流体膜行为决定整个微接触性能. 当膜厚较大时,这里边界膜可看成纳米级薄膜. 由于上接触表面处有限的剪应力承受能力,边界膜可于上接触表面滑移. 设下接触表面处剪应力承受能力很大而边界膜在下接触表面不滑移. 由于边界膜-接触表面间相互作用,边界膜黏度、密度和剪切弹性模量均沿膜厚变化,在理论分析中使用它们的等效值,这些值与边界膜厚度有关. 流体膜在两个接触表面均不发生滑移,分析中不考虑流体膜剪切弹性模量. 流体膜采用传统分析法. 给出了理论分析和若干变工况参数下的计算结果.      

Phenomenological study of parabolic and spherical indentation of elastic-ideally plastic material

A phenomenological study of parabolic and spherical indentation of elastic ideally plastic materials was carried out by using precise results of finite elements calculations. The study shows that no “pseudo-Hertzian” regime occurs during spherical indentation. As soon as the yield stress of the indented material is exceeded, a deviation from the, purely elastic Hertzian contact behaviour is found. Two elastic–plastic regimes and two plastic regimes are observed for materials of very large Young modulus to Yield stress ratio, E/σ_y. The first elastic–plastic regime corresponds to a strong evolution of the indented plastic zone. The first plastic regime corresponds to the commonly called “fully plastic regime”, in which the average indentation pressure is constant and equal to about three times the yield stress of the indented material. In this regime, the contact depth to penetration depth ratio tends toward a constant value, i.e. h_c/h = 1.47. h_c/h is only constant for very low values of yield strain (σ_y/E lower than 5 × 10^?6) when aE¹/Rσ_y is higher than 10,000. The second plastic regime corresponds to a decrease in the average indentation pressure and to a steeper increase in the pile-up. For materials with very large E/σ_y ratio, the second plastic regime appears when the value of the non-dimensional contact radius a/R is lower than 0.01. In the case of spherical and parabolic indentation, results show that the first plastic regime exists only for elastic-ideally plastic materials having an E/σ_y ratio higher than approximately 2.000.     

A study on the bending stiffness of single-walled carbon nanotubes and related issues

Single-walled carbon nanotubes (SWCNTs) are usually modeled as elastic tubes and their bending stiffness D is often related to their axial stretching modulus E (Young's modulus) as in mechanics of materials (i.e. D=EI where I is the moment of inertia of the tube). However, recent studies show that large discrepancies may exist when this relationship is used to predict Young's modulus of carbon nanotubes (CNTs) through bending dominated deformations. In the present paper, the bending stiffness of SWCNTs and some related issues are investigated by the combined use of the molecular-mechanics (M-M) model and the deformation mapping technique. Based on the analysis results, the contradictions mentioned above can be explained well. Furthermore, an analytical expression for the bending stiffness of SWCNTs is also presented. It shows that the bending stiffness of a SWCNT is approximately proportional to the cube of its radius which agrees well with the existing molecular dynamics simulation and continuum theory based results.     

A critical state sand model with elastic–plastic coupling

Soil elastic moduli are highly pressure-dependent. Experimental findings have indicated that the elastic shear modulus of sands depends on p^χ, where p is mean principal effective stress and χ is a non-dimensional parameter. χ practically remains unchanged for shear strains less than 10⁻⁵ where the mechanical behavior is purely elastic. However, experiments have revealed that the emergence of plasticity for shear strains larger than 10⁻⁵ provokes a gradual increase in χ. Technically, this observation is an elastic–plastic coupling effect in which plasticity causes to change the elastic characteristics. Here, this issue is considered in hyper-elasticity framework in conjunction with a critical state compatible bounding surface plasticity platform for granular soils. To this aim, constitutive equations linking χ to a proper kinematic hardening parameter are presented. Then, using the proposed approach, a hyper-elastic theory is modified to consider the mentioned elastic–plastic coupling effect in the whole domain of the elastoplastic behavior. Adopting the improved hyper-elasticity necessitates the modification of a number of basic plasticity platform elements. In this regard, dilatancy and plastic hardening modulus of the bounding surface platform are modified. Successful performance of the modified constitutive model is presented against experimental data of loading/unloading triaxial tests.     

Macroscopic elastic properties of regular lattices

In this paper we study analytically the elastic properties of the 2-D and 3-D regular lattices consisting of bonded particles. The particle-scale stiffnesses are derived from the given macroscopic elastic constants (i.e. Young's modulus and Poisson's ratio). Firstly a bonded lattice model is presented. This model permits six kinds of relative motion and corresponding forces between each bonded particle pair. By comparing the strain energy distributions between the discrete lattices and the continuum, the explicit relationship between the microscopic and macroscopic elastic parameters can be obtained for the 2-D hexagonal lattice and the 3-D hexagonal close-packed and face-centered cubic structures. The results suggest that the normal stiffness is determined by Young's modulus and the particle size (in 3-D), and that the ratio of the shear to normal stiffness is related to Poisson's ratio. Rotational stiffness depends on the normal stiffness, shear stiffness and particle sizes. Numerical tests are carried out to validate the analytical results. The results in this paper have theoretical implications for the calibration of the spring stiffnesses in the Discrete Element Method.     

The elastic modulus of particulate composites using the concept of a mesophase

A theoretical model for the evaluation of the elastic modulus in particulate composites has been developed. The method takes into account the existence of a mesophase between main phases, which constitutes an important parameter influencing the behaviour of a composite material. This layer between the matrix and filler develops different physico-chemical properties from those of the constituent phases and variable ones along its thickness. The effect of the progressive variation of the elastic modulus of the mesophase on the modulus of the composite was estimated by applying various simple laws of variation. Convenient laws of variation were introduced, varying from a simple one, assuming a linear law, to a more refined one using a parabolic law. Experimental results with particulates, based on iron-filled epoxy composites, compared satisfactorily with other models. However, the model based on a parabolic law was superior to all others on physical grounds.List of symbols E elastic modulus - v Poisson's ratio - volume fraction - r radius - V volume - K bulk modulus - r thickness of mesophase - parameter which depends on the sudden changes of the heat capacity for the filled and unfilled polymer C _p ^f and C _p ⁰      

Reprint of “Size dependence and orientation dependence of elastic properties of ZnO nanofilms” [In. J. Solids Struct. 45 (2008) 1730–1753]

The elastic properties of ZnO nanofilms with different film thickness, surface orientations and loading directions are investigated by using molecular mechanics (MM) method. The size dependence of elastic properties is relevant to both the film surface crystallographic orientation and loading direction. Both atomic structure analysis and energy calculation are employed to identify the mechanisms of the size-dependent elastic properties, under different loading directions and surface orientations. Upon small axial deformation, the relationship between intralayer and interlayer bond length variation and film elastic stiffness is established; it is found that the atomic layers with larger bond length variation have higher elastic stiffness. The strain energies of atomic layers of ZnO nanofilm and bulk are decoupled, from which the stiffness of film surface, intralayers, and interlayers are derived and compared with their bulk counterparts. The surface stiffness is found to be much lower than that of the interior layers and bulk counterpart, and with the decrease of film thickness, the residual tension-stiffened interior atomic layers are the main contributions of the increased elastic modulus of ZnO nanofilms.     

Size dependence and orientation dependence of elastic properties of ZnO nanofilms

The elastic properties of ZnO nanofilms with different film thickness, surface orientations and loading directions are investigated by using molecular mechanics (MM) method. The size dependence of elastic properties is relevant to both the film surface crystallographic orientation and loading direction. Both atomic structure analysis and energy calculation are employed to identify the mechanisms of the size-dependent elastic properties, under different loading directions and surface orientations. Upon small axial deformation, the relationship between intralayer and interlayer bond length variation and film elastic stiffness is established; it is found that the atomic layers with larger bond length variation have higher elastic stiffness. The strain energies of atomic layers of ZnO nanofilm and bulk are decoupled, from which the stiffness of film surface, intralayers, and interlayers are derived and compared with their bulk counterparts. The surface stiffness is found to be much lower than that of the interior layers and bulk counterpart, and with the decrease of film thickness, the residual tension-stiffened interior atomic layers are the main contributions of the increased elastic modulus of ZnO nanofilms.     

用能量法分析层状材料弹性接触问题

利用能量法分析了层状材料（薄膜／基体）弹性接触问题，得到了具有一阶精度的闭合解，给出了求解薄膜弹性模量和泊松比的表达式，并与有限元的数值解进行了比较。二者比较结果表明：在工程材料范围内，理论解与数值解相差在6％以内；同时表明单相材料中剪切模量与弹性模量之间的关系也适用层状材料中的薄膜材料。在数值解的基础上，讨论了薄膜厚度与压头半径的比值对求解精度的影响，发现此比值对精度影响不大。通过对层状材料等效泊松比与等效弹性模量的定义，给出了用压痕实验测定薄膜泊松比与弹性模量的方法。     

Modification of the elastic properties of nanostructures with surface charges in applied electric fields

The influences of applied electric fields and surface charges on elastic modulus of nanostructures such as nanowires and nanofilms are investigated within the framework of classic continuum mechanics. Under an applied electric field, the surfaces of structures are subjected to the electrostatic forces (negative pressure) along the direction of the electric field, and the resulting surface charges also change the surface mechanical properties due to the Hellman–Feynman (H–F) forces. Through incorporating the surface energy from the negative pressure and the H–F forces into surface free energy, the exact and analytical expressions of the effective elastic modulus of nanowires and nanofilms are addressed by considering the surface energy effects on the elastic modulus of nanostructures, which involves the contribution of the applied electric field and surface charges. The numerical results indicate that applied electric fields parallel to the axis of the nanowire and nanofilms enhance the transverse Young's modulus while reducing axial modulus of nanostructures. The effective modulus of nanowires and nanofilms with lateral surface charges depends on the surface charges density and the sign of the charges. In addition, the effect of electric field and surface charges on Young's moduli of nanowires and nanofilms has been found to be sensitive to structural geometric dimensions such as the thickness of the film and the diameter of the wire.