Dynamics of a system comprising a pre-stressed orthotropic layer and pre-stressed orthotropic half-plane under the action of a moving load

This paper investigates the dynamic response to a moving load of a system comprising an initially stressed covering layer and initially stressed half-plane, within the scope of the piecewise-homogeneous body model utilizing three-dimensional linearized wave propagation theory in the initially stressed body. It was assumed that the materials of the layer and half-plane are anisotropic (orthotropic), and that the velocity of the line-located moving load is constant as it acts on the free face of the covering layer. The investigations were made for a two-dimensional problem (plane-strain state) under subsonic velocity of the moving load for complete and incomplete contact conditions. Corresponding numerical results were obtained for the stiffer layer and soft half-plane system in which the modulus of elasticity of the covering layer material (for the moving direction of the load) is greater than that of the half-plane material, which was assumed to be isotropic. Numerical results are presented and discussed for the critical velocity and stress distribution for various values of the problem parameters. In particular, it was established that, the critical velocity of the moving load is controlled mainly with a Rayleigh wave speed of a half-plane material and the initial stretching of the covering layer causes to increase these values.     

Dynamics of a system comprising an orthotropic layer and orthotropic half-plane under the action of an oscillating moving load

This paper investigates the dynamic response to a time-harmonic oscillating moving load of a system comprising a covering layer and half-plane, within the scope of the piecewise-homogeneous body model utilizing of the exact equations of the linear theory of elastodynamics. It is assumed that the materials of the layer and half-plane are anisotropic (orthotropic), and that the velocity of the line-located time-harmonic oscillating moving load is constant as it acts on the free face of the covering layer. Our investigations were carried out for a two-dimensional problem (plane-strain state) under subsonic velocity for a moving load in complete and incomplete contact conditions. The corresponding numerical results were obtained for the stiffer layer and soft half-plane system in which the modulus of elasticity of the covering layer material (for the moving direction of the load) is greater than that of the half-plane material. Numerical results are presented and discussed for the critical velocity, displacement and stress distribution for various values of the problem parameters. In particular, it is established that the critical velocity of the moving load is controlled mainly with a Rayleigh wave speed of a half-plane material and the existence of the oscillation of the moving load causes two types of critical velocity to appear: one of which is less, but the other one is greater than that attained for the case where the mentioned oscillation is absent.     

A friction contact stiffness model of fractal geometry in forced response analysis of a shrouded blade

To investigate the nonlinear vibration behavior of a shrouded blade with friction dynamic contact interface, a friction contact stiffness model is proposed to describe the friction force at different rough interfaces and different normal loads. In the proposed model, the friction contact interface is discretized to a series of friction contact pairs and each of them can experience stick, slip, or separate states. Fractal geometry is used to simulate the topography of contact surfaces. The contact stiffness is calculated using the Hertz contact theory and fractal geometry, which is related to contact interfaces parameters including normal load, roughness, Young??s modulus, and Poisson??s ratio. The trajectory tracking method is used to predict the friction force and it is not necessary to judge the transition condition among stick, slip, and separate states. It is suitable for complicated periodic motion of the contact interfaces. The forced response of a real shrouded blade is predicted using the proposed model and the multi-harmonic balance method. The effect of surface roughness, initial normal load, and contact area on the forced response of a shrouded blade is studied. It is shown that contact stiffness increases with normal load and fractal dimension. The resonant amplitude is sensitive to the initial normal load and contact surface roughness. The response can be influenced by the contact area, which is an important parameter for blade designers.     

Low-velocity impact response of sandwich beams with functionally graded core

The problem of low-speed impact of a one-dimensional sandwich panel by a rigid cylindrical projectile is considered. The core of the sandwich panel is functionally graded such that the density, and hence its stiffness, vary through the thickness. The problem is a combination of static contact problem and dynamic response of the sandwich panel obtained via a simple nonlinear spring-mass model (quasi-static approximation). The variation of core Young’s modulus is represented by a polynomial in the thickness coordinate, but the Poisson’s ratio is kept constant. The two-dimensional elasticity equations for the plane sandwich structure are solved using a combination of Fourier series and Galerkin method. The contact problem is solved using the assumed contact stress distribution method. For the impact problem we used a simple dynamic model based on quasi-static behavior of the panel—the sandwich beam was modeled as a combination of two springs, a linear spring to account for the global deflection and a nonlinear spring to represent the local indentation effects. Results indicate that the contact stiffness of the beam with graded core increases causing the contact stresses and other stress components in the vicinity of contact to increase. However, the values of maximum strains corresponding to the maximum impact load are reduced considerably due to grading of the core properties. For a better comparison, the thickness of the functionally graded cores was chosen such that the flexural stiffness was equal to that of a beam with homogeneous core. The results indicate that functionally graded cores can be used effectively to mitigate or completely prevent impact damage in sandwich composites.     

Dynamic homogenization of a complex geophysical medium by inversion of its near-field seismic response

The present study originated in the forward problem of the prediction of the effects of seismic waves (generated by impulsive deep-down sources) in urban areas. The traditional, numerically-intensive approaches to this problem have not, until now, given rise to simple theoretical paradigms which might explain how and why the (often-destructive) response of cities is conditioned by factors such as the city density, building average height, average building composition, site geometry and composition, and characteristics of the solicitation such as incident angles, polarization and frequency. We propose to homogenize the city in order to simplify, and make possible the understanding of, the site-city-solicitation interaction. This homogenization is treated as an inverse problem, i.e., by which we: (1) generate near-field response data for a ‘real’ city, (2) replace (initially by thought) the city by a homogeneous (surrogate) layer above, and in firm contact, with the underlying site, (3) compute the response of the surrogate layer/site response for various trial constitutive properties, (4) search for the global minimum of the discrepancy between the response data and the various trial parameter responses (5) attribute the homogenized properties of the city to the surrogate layer for which the minimum of the discrepancy is attained. We carry out this five-step procedure for a host of ‘real’ city and solicitation parameters, notably the frequency. The result is that: (i) for low frequencies and/or large city densities, the effective constitutive properties are their static equivalents, i.e., the effective shear modulus is the product of a factor related to the city density with the shear modulus of a generic substructure of the city and the effective complex velocity is equal to the complex velocity of the said generic substructure, (2) at higher frequencies and/or smaller city densities, the effective constitutive properties are dispersive and do not take on a simple mathematical form, with this dispersion compensating for the discordance between the ways the inhomogeneous city structure and the homogeneous surrogate layer respond to the seismic wave. For typical seismic solicitation frequencies, the city, represented as a layer with static homogenized properties, is quite adequate to account for the principal features of the response (notably those of the time-domain response). The model of the layer with dispersive homogenized properties is more suitable to account for such features as resonances due to the excitation of surface wave modes.     

Two-dimensional sliding frictional contact of functionally graded materials

A multi-layered model for sliding frictional contact analysis of functionally graded materials (FGMs) with arbitrarily varying shear modulus under plane strain-state deformation has been developed. Based on the fact that an arbitrary curve can be approached by a series of continuous but piecewise linear curves, the FGM is divided into several sub-layers and in each sub-layers the shear modulus is assumed to be a linear function while the Poisson's ratio is assumed to be a constant. In the contact area, it is assumed that the friction is one of Coulomb type. With this model the fundamental solutions for concentrated forces acting perpendicular and parallel to the FGMs layer surface are obtained. Then the sliding frictional contact problem of a functionally graded coated half-space is investigated. The transfer matrix method and Fourier integral transform technique are employed to cast the problem to a Cauchy singular integral equation. The contact stresses and contact area are calculated for various moving stamps by solving the equations numerically. The results show that appropriate gradual variation of the shear modulus can significantly alter the stresses in the contact zone.     

A Simplified Treatise of the Scott Bond Testing Method

The Scott bond test method has been used extensively in the paper industry over the years as a means to assess the bond strength of paper. The method has been a subject of some controversy lately since it does not always correlate to the sensitivity of the material to fracture by delamination. To gain some further insight into which parameters govern the fracture process in a Scott bond test a simplified approach has been chosen in order to formulate an analytical mathematical/mechanical model of the test. The model is dynamic in the sense that inertia effects are included. The material model utilised is a simple cohesive theory that assumes a linear behaviour between stress and crack opening when the material has started to degrade. This choice of material model makes the mathematical model very nonlinear. In fact, a system of three coupled nonlinear second order partial differential equations have to be solved and adjusted to the correct initial conditions. The material parameters needed for the model are the elastic modulus in the thickness direction, the transverse shear (elastic) modulus, the tensile strength (in the thickness direction) and the fracture work (per unit area) for a delamination crack. To investigate the ability of the model, a Scott bond testing apparatus have been equipped with a piezoelectric load sensor. The load cell was mounted on the apparatus’ pendulum so that the load acting on the sample holder could be recorded during the whole impact stage. This was done for a number of different initial velocities of the pendulum and it is found that the model gives a fair prediction of the contact load.     

Point loading on hyperbolic paraboloidal shells

A combination of simple theory and experiment is used to investigated a complex structural problem. The problem is the point loading of hyperbolic paraboloidal shells. The solution is based on an analogy between the response of these shells to a point load and the localized behavior of a flat plate that rests on an elastic foundation. The subgrade modulus of this analogous foundation is derived from the measured deflection of a point load on shell models. The value of the modulus was found to be about one-fifth that of a spherical shell with the same magnitude of principal curvature.     

具有梯度界面层的双材料悬臂梁解析解

采用应力函数法,求得了具有弹性模量沿高度线性变化的梯度界面层的双材料悬臂梁在均布载荷作用下的应力和位移解析解。该解可退化为双材料梁、弹性模量沿整个梁高线性变化的梯度梁以及均质材料梁的情况,退化为均质材料梁时与已有结果一致。通过一具体算例将得到的解析解与有限元解进行了比较,两者吻合较好。并讨论了梯度界面层的高度变化对梁中的应力和梁端挠度的影响。结果表明,在梁的总高度不变的情况下,增加梯度界面层的高度可减小弯曲应力和梁端挠度,而对挤压应力和切应力的影响很小。     

用拉格朗日函数分析不同模量静不定桁架内力

在外载荷作用下的不同模量静不定桁架平衡问题,是任意有限多个自变量的多元函数在任意有限多个约束条件下的极值问题,对采用拉格朗日乘数法求解此类极值问题进行了数学证明．通过求解不同模量静不定桁架极限载荷的几个算例,阐述拉格朗日乘数法在计算不同模量静不定桁架极限载荷中的应用．研究结果表明：采用拉格朗日乘数法求解不同模量静不定桁架极限载荷的通用性较强,用拉格朗日乘数法求解不同模量静不定桁架极限载荷的方法不但克服了常规方法需利用几何关系建立协调方程的缺陷,且具有力学概念清晰直观、计算过程简便、便于工程设计人员在实际中掌握和应用．     

三体介质温成形界面接触换热特性的试验研究

温成形摩擦界面模具与工件之间的传热特性对工件质量和模具寿命有重要影响,固体粉末介质导入该摩擦副可实现高温润滑,但其传热特性与传统加工方式的有很大不同.采用稳态法自行设计了三体界面的传热特性试验,研究和分析了界面温度、接触载荷、层厚对带有石墨粉和氧化铝粉润滑层的H62铜合金和45钢之间的三体界面接触换热系数的影响.结果表明:带有石墨粉润滑层的三体界面接触换热系数随温度的增加先升高后降低,随载荷和层厚的增加先缓慢增加后迅速增加;带有氧化铝粉润滑层的三体界面接触换热系数随温度的变化缓慢升高,与载荷基本成线性关系,随层厚的增加而降低.温度改变了固体润滑剂的材料热阻和上下试样表面硬度及氧化层厚度,载荷改变了三体界面实际接触面积和接触属性,层厚决定能否完全隔开上下试样,不同物性固体润滑剂决定了其材料热阻在三体界面接触热阻中的主次关系.     

A Simple Model for a Fluid-Filled Open-Cell Foam

A simple microstructure model is used to describe a fluid-filled open-cell foam. In the simplest case it consists of parallel elastic plates with gaps between them, which are filled with a Newtonian fluid. We assume that the load applied to this model material is uniaxial. The constitutive equation is formulated with the pressure of the fluid as an inner variable. The model yields an evolutional equation for the fluid pressure which itself is a field equation, that is a partial differential equation in time and space coordinates. This differential equation is solved for an instantaneously applied constant load and for a harmonically oscillating load. The solution of the differential equation, in combination with the constitutive equation leads to a relation between mean applied load and global strain of the test specimen. Finally, we obtain the creep compliance and the complex modulus of the foam material, respectively. The influence of different geometries of the foam and of different material behaviour of the matrix and fluid on the creep compliance and the complex modulus is discussed.     

Free vibration of circular cylindrical shell with constrained layer damping

Free vibration characteristics of circular cylindrical shell with passive constrained layer damping (PCLD) are presented. Wave propagation approach rather than finite element method, transfer matrix method, and Rayleigh-Ritz method is used to solve the problem of vibration of PCLD circular cylindrical shell under a simply supported boundary condition at two ends. The governing equations of motion for the orthotropic cylindrical shell with PCLD are derived on the base of Sanders’ thin shell theory. Numerical results show that the present method is more effective in comparison with other methods. The effects of the thickness of viscoelastic core and constrained layer, the elastic modulus ratio of orthotropic constrained layer, the complex shear modulus of viscoelastic core on frequency parameter, and the loss factor are discussed.     

不锈钢表面等离子合金化渗铜层的纳米力学性能研究

为了研究不锈钢渗铜层的纳米力学性能,本文采用等离子表面合金化技术在304不锈钢表面制备了渗铜层。利用纳米压入硬度仪,采用连续刚度测试法,对渗铜层以及不锈钢基体的表面和横截面的纳米力学性能进行了测试,得到了纳米压入过程的载荷-压入深度曲线,发现渗铜层抵抗外载荷的能力低于不锈钢基体。并得到了渗铜层和不锈钢基体的表面、横截面方向的硬度以及杨氏模量,经对比得到渗铜层的杨氏模量和硬度都要比不锈钢基体的低,并且渗铜层的力学性能表现出各向异性。对纳米压入的数据进行分析,发现渗铜层在小尺度压入时硬度和杨氏模量表现出明显的"尺度效应"。     

Measuring the Elastic Properties of Freestanding Thick Films Using a Nanoindenter-Based Bending Test

Films with thickness ranging from 10 to 100 μm are increasingly being used as the structural components of microelectromechanical systems (MEMS). Measuring the mechanical properties of these thick films is essential for enabling the design of MEMS with high performance and sufficient reliability. In this paper, we present a simple and convenient method for measuring the elastic modulus of thick films by loading a clamped circular film using a spherical tip. The test is implemented using a commercial nanoindenter so that the load and displacement can be measured with resolution of micronewtons and nanometers, respectively. Robust protocols have been developed for implementing the test within the constraints imposed by the nanoindenter. A crucial component of these protocols is a method for selecting loads to ensure deformation in the elastic bending regime and to minimize the relative contribution of contact indentation. The accuracy and utility of the nanoindenter-based bending test are discussed using measurements on thick films of aluminum and a standard epoxy.     

Surface instability of an elastic half space with material properties varying with depth

If a body with a stiffer surface layer is loaded in compression, a surface wrinkling instability may be developed. A bifurcation analysis is presented for determining the critical load for the onset of wrinkling and the associated wavelength for materials in which the elastic modulus is an arbitrary function of depth. The analysis leads to an eigenvalue problem involving a pair of linear ordinary differential equations with variable coefficients which are discretized and solved using the finite element method.The method is validated by comparison with classical results for a uniform layer on a dissimilar substrate. Results are then given for materials with exponential and error-function gradation of elastic modulus and for a homogeneous body with thermoelastically induced compressive stresses.     

Viscoelastic composite with a cracked layer under anti-plane shear

A viscoelastic cracked sandwich composite under anti-plane load is presented via the application of the Laplace transform and the complex variable formulation. Using a special technique of analytical continuation associated with the method of image, the triple-layered problem is reduced to a singular function by simple algebraic manipulations. The superposition method and the singular integral equation is applied to deal with the concentrated load and crack problem. Finally, some typical viscoelastic models and its corresponding stresses and stress intensity factors are also discussed. The results show that the time dependent stress intensity factor is affected by the relative strength of creep compliances of each layer.     

Analytical solutions for bending and flexure of helically reinforced cylinders

There has been a great deal of interest in the problems of modelling cables and ropes. A recent review by Cardou and Jolicoeur [Appl. Mech. Rev. 50 (1997) 1] considers the modelling of a cable which consists of a central core surrounded by one or several helically wound wire layers. One approach has been to consider the deformations of an individual helical wire and to synthesise the model of a cable by using contact conditions between the various wires. Other authors have adopted a continuum approach regarding each layer as a transversely isotropic material whose principal direction is along a helix surrounding the central axis of the cable. In each layer the helix angle is constant so that, when referred to cylindrical polar co-ordinates, the cylinder has a constant stiffness matrix in each layer. The intention in this paper is to use the continuum approach and describe the analytical solutions that govern the simple bending, flexure, or bending under a uniform load, of an anisotropic elastic cylinder consisting of a single material of this type. The extension of this work to a composite cylinder consisting of several concentric layers, surrounding a central core, which are either bonded together or make a frictionless contact, is briefly described.     

基于SWT方法的钢绞线索微动疲劳特性分析

为得到钢绞线索丝间接触区的应力场分布并预测微动疲劳裂纹萌生位置和微动疲劳寿命,本文利用参数化方法建立了精细化的钢绞线拉索有限元模型,包括整索模型和不同层丝间接触区域的局部精细化子模型．分析了钢绞线索在两种交变荷载工况下的应力场变化情况,并基于多轴疲劳SWT(Smith-Watson-Topper)临界平面法进行了疲劳特性分析和疲劳寿命预测．主要结论如下：钢绞线索内接触区边缘处的微动幅值较大,中心处几乎没有相对滑动,微动疲劳的初始裂纹萌生点位于接触区域边缘;经不同区域子模型分析比较,在轴向循环荷载作用下,外层钢丝的接触区域比内层钢丝更易发生微动疲劳损伤;在横向位移循环荷载作用下,同层钢丝因位置角度不同而产生了较大的疲劳特性差异,且相比轴向循环拉伸,该工况下最不利单丝的微动疲劳寿命更低;与非接触区域相比,接触区的疲劳寿命大幅降低,微动现象对钢绞线索的抗疲劳性能有明显降低作用．     

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.

---