考虑几何与材料及静风荷载的非线性因素的大跨径桥梁静风稳定分析法

针对现有的桥梁静风稳定分析方法中存在的问题，提出了增量与内外两重迭代相结合的新方法，并且考虑了结构几何、材料和静风荷载非线性。在上述方法的基础上，编制了桥梁非线性空气静力稳定分析程序BNAP，并进行了相应的算例分析，所得结果表明该方法具有计算稳定和速度快的优点。最后，以一座主跨1000米的斜拉桥为例，分析了结构几何非线性、材料非线性和静风荷载非线性对大跨径桥梁空气静力稳定性的影响。     

构造杂交应力单元的柔度矩阵H对角化方法

证明了杂交元柔度矩阵 H非奇异的充分必要条件是假设应力模式线性无关 ;以及等价应力模式形成相同的杂交元。在此基础上建立了假设应力模式的 Hilbert子空间 ,从而可以利用 Schmidt方法简单地得到等价的正交应力模式 ,实现了柔度矩阵 H对角化 ,使得杂交元形成过程中完全避免了繁杂的矩阵求逆运算 ,提高了杂交元分析的计算效率 ,特别在柔度矩阵不容易显式求逆的材料非线性分析中更具有实际意义     

Numerical study of sliding wear caused by a loaded pin on a rotating disc

A computational approach is proposed to predict the sliding wear caused by a loaded spherical pin contacting a rotating disc, a condition typical of the so-called pin-on-disc test widely used in tribological studies. The proposed framework relies on the understanding that, when the pin contacts and slides on the disc, a predominantly plane strain region exists at the centre of the disc wear track. The wear rate in this plane strain region can therefore be determined from a two dimensional idealisation of the contact problem, reducing the need for computationally expensive three dimensional contact analyses. Periodic unit cell techniques are used in conjunction with a ratchetting-based failure criterion to predict the wear rate in the central plane strain region. The overall three dimensional wear rate of the disc is then determined by scaling the plane strain wear rate with a conversion factor related to the predicted shape of the wear track. The approach is used to predict pin-on-disc test data from an Al-Si coating using a tungsten carbide pin. The predicted results are found to be consistent with measured data.     

Neural network optimization of material composition of a functionally graded material plate at arbitrary temperature range and temperature rise

Summary A neural network model is applied to optimization problems of material compositions for a functionally graded material plate with arbitrarily distributed and continuously varied material properties in the thickness direction. Unsteady temperature distribution is evaluated by taking into account the bounds of the number of the layers. Thermal stress components for an infinite functionally graded material plate are formulated under traction-free mechanical conditions. As a numerical example, a plate composed of zirconium oxide and titanium alloy is considered. In the optimization problem of minimizing the thermal stress distribution, the numerical calculations are carried out making use of the neural network. The optimum material composition is determined by taking into account the effect of temperature-dependence of material properties. The results obtained by neural network and ordinary nonlinear programming method are compared. Received 3 March 1998; accepted for publication 22 May 1998     

Dynamical analysis of cavitation for a transversely isotropic incompressible hyper-elastic medium: Periodic motion of a pre-existing micro-void

In this paper, the dynamical cavitation behavior is analyzed for a sphere composed of a class of transversely isotropic incompressible hyper-elastic materials, where there is a pre-existing micro-void in the interior of the sphere. A second-order non-linear ordinary differential equation that governs the motion of the initial micro-void is obtained by using the boundary conditions. On analyzing the qualitative properties of the solutions of the differential equation, some interesting conclusions are proposed. It is proved that the number of equilibrium points of the differential equation depends on the values of the material parameters, and that the phase diagrams of the equation are closed, smooth and convex trajectories. For any prescribed surface tensile dead-loads, the motion of the initial micro-void undergoes a non-linear periodic oscillation. The dependence of the periodic motion of the initial micro-void on material parameters and the radius of the initial micro-void is examined, and numerical results are also provided. It is worth pointing out that the conclusions in this paper can be used to describe approximately the physical implications of the dynamical formation of a cavity in the sphere.     

Two-dimensional transient dynamic response of orthotropic layered media

In this study, two-dimensional transient dynamic response of orthotropic plane layered media is investigated. The plane multilayered media consist of N different generally orthotropic, homogeneous and linearly elastic layers with different ply angles. In the generally orthotropic layer, representing a ply reinforced by unidirectional fibers with an arbitrary orientation angle, the principal material directions do not coincide with body coordinate axes. The solution is obtained by employing a numerical technique which combines the use of Fourier transform with the method of characteristics. The numerical results are displayed in curves denoting the variations of stress and displacement components with time at different locations. These curves clearly reveal, in wave profiles, the scattering effects caused by the reflections and refractions of waves at the boundaries and at the interfaces of the layers, and also the effects of anisotropy caused by fiber orientation angle. The curves properly predict the sharp variations in the response at the neighborhood of the wave fronts, which shows the power of the numerical technique employed in the study. By suitably adjusting the elastic constants, the results for multilayered media with transversely isotropic layers, or layers with cubic symmetry, or isotropic layers can easily be obtained from the general formulation. Furthermore, solutions for some special cases, including Lamb’s problem for an elastic half-space, are obtained and compared with the available solutions in the literature and very good agreement is found. Preliminary version presented at the Second International Congress on Mechatronics (MECH2K3), Graz, Austria, July 14-17, 2003.     

Hertzian fracture at unloading

Hertzian fracture through indentation of flat float glass specimens by steel balls has been examined experimentally. Initiation of cone cracks has been observed and failure loads together with contact and fracture radii determined at monotonically increasing load but also during unloading phases. Contact of dissimilar elastic solids under decreasing load may cause crack inception triggered by finite interface friction and accordingly the coefficient of friction was determined by two different methods. In order to make relevant predictions of experimental findings, a robust computational procedure has been developed to determine global and local field values in particular at unloading at finite friction. It was found that at continued loading it is possible to specify in advance how the contact domain divides into invariant regions of stick and slip. The maximum tensile stress was found to occur at the free surface just outside the contact contour, the relative distance depending on the different elastic compliance properties and the coefficient of friction. In contrast, at unloading invariance properties are lost and stick/slip regions proved to be severely history dependant and in particular with an opposed frictional shear stress at the contact boundary region. This causes an increase of the maximum tensile stress at the contour under progressive unloading. Predictions of loads to cause crack initiation during full cycles were made based on a critical stress fracture criterion and proved to be favourable as compared to the experimental results.     

INCREMENTAL MICRO-MECHANICAL MODEL OF PLAIN WOVEN FABRIC

I. INTRODUCTION Previous research on woven fabric ignored the micro-weaving structures in fabric and modeled thewoven fabric as a ?exible orthotropic plate[1??3]. However, such a model is unable to predict certainbuckling phenomena of fabric commonly ob…     

Limiting state of an elastoplastic orthotropic plate with a periodic system of collinear cracks

A modified Dugdale model is used to study the fracture of an orthotropic elastoplastic plate with a periodic system of rectilinear cracks. The material of the plate obeys a general yield criterion. The general form of solution is obtained in terms of Kolosov-Muskhelishvili potentials. The size of the plastic zone is expressed in terms of the external load and geometrical parameters. The equations for the determination of the stresses in the plastic zone and the crack opening displacement are derived. The effect of anisotropy on the formation of the plastic zones at the crack tip is examined __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 5, pp. 80–88, May 2007.