Tensionless Contact Of A Finite Beam Resting On Reissner Foundation

A more generalized model of a beam resting on a tensionless Reissner foundation is presented. Compared with the Winkler foundation model, the Reissner foundation model is a much improved one. In the Winkler foundation model, there is no shear stress inside the foundation layer and the foundation is assumed to consist of closely spaced, independent springs. The presence of shear stress inside Reissner foundation makes the springs no longer independent and the foundation to deform as a whole. Mathematically, the governing equation of a beam on Reissner foundation is sixth order differential equation compared with fourth order of Winkler one. Because of this order change of the governing equation, new boundary conditions are needed and related discussion is presented. The presence of the shear stress inside the tensionless Reissner foundation together with the unknown feature of contact area/length makes the problem much more difficult than that of Winkler foundation. In the model presented here, the effects of beam dimension, gap distance, loading asymmetry and foundation shear stress on the contact length are all incorporated and studied. As the beam length increases, the results of a finite beam with zero gap distance converge asymptotically to those obtained by the previous model for an infinitely long beam. (C) 2008 Elsevier Ltd. All rights reserved.

Thermal fatigue on pistons induced by shaped high power laser. Part II: Design of spatial intensity distribution via numerical simulation

In the laser induced thermal fatigue simulation test on pistons, the high power laser was transformed from the incident Gaussian beam into a concentric multi-circular pattern with specific intensity ratio. The spatial intensity distribution of the shaped beam, which determines the temperature field in the piston, must be designed before a diffractive optical element (DOE) can be manufactured. In this paper, a reverse method based on finite element model (FEM) was proposed to design the intensity distribution in order to simulate the thermal loadings on pistons. Temperature fields were obtained by solving a transient three-dimensional heat conduction equation with convective boundary conditions at the surfaces of the piston workpiece. The numerical model then was validated by approaching the computational results to the experimental data. During the process, some important parameters including laser absorptivity, convective heat transfer coefficient, thermal conductivity and Biot number were also validated. Then, optimization procedure was processed to find favorable spatial intensity distribution for the shaped beam, with the aid of the validated FEM. The analysis shows that the reverse method incorporated with numerical simulation can reduce design cycle and design expense efficiently. This method can serve as a kind of virtual experimental vehicle as well, which makes the thermal fatigue simulation test more controllable and predictable. (C) 2007 Elsevier Ltd. All rights reserved.

First-passage time of Duffing oscillator under combined harmonic and white-noise excitations

The first-passage time of Duffing oscillator under combined harmonic and white-noise excitations is studied. The equation of motion of the system is first reduced to a set of averaged Ito stochastic differential equations by using the stochastic averaging method. Then, a backward Kolmogorov equation governing the conditional reliability function and a set of generalized Pontryagin equations governing the conditional moments of first-passage time are established. Finally, the conditional reliability function, and the conditional probability density and moments of first-passage time are obtained by solving the backward Kolmogorov equation and generalized Pontryagin equations with suitable initial and boundary conditions. Numerical results for two resonant cases with several sets of parameter values are obtained and the analytical results are verified by using those from digital simulation.

Size effects on the melting of nickel nanowires: a molecular dynamics study

The melting process of nickel nanowires are simulated by using molecular dynamics with the quantum Sutten-Chen many-body force field. The wires studied were approximately cylindrical in cross-section and periodic boundary conditions were applied along their length; the atoms were arranged initially in a face-centred cubic structure with the [0 0 1] direction parallel to the long axis of the wire. The size effects of the nanowires on the melting temperatures are investigated. We find that for the nanoscale regime, the melting temperatures of Ni nanowires are much lower than that of the bulk and are linear with the reciprocal of the diameter of the nanowire. When a nanowire is heated up above the melting temperature, the neck of the nanowire begins to arise and the diameter of neck decreases rapidly with the equilibrated running time. Finally, the breaking of nanowire arises, which leads to the formation of the spherical clusters. (C) 2004 Elsevier B.V. All rights reserved.

Modeling, numerical methods, and simulation for particle-fluid two-phase flow problems

In this paper, we study the issues of modeling, numerical methods, and simulation with comparison to experimental data for the particle-fluid two-phase flow problem involving a solid-liquid mixed medium. The physical situation being considered is a pulsed liquid fluidized bed. The mathematical model is based on the assumption of one-dimensional flows, incompressible in both particle and fluid phases, equal particle diameters, and the wall friction force on both phases being ignored. The model consists of a set of coupled differential equations describing the conservation of mass and momentum in both phases with coupling and interaction between the two phases. We demonstrate conditions under which the system is either mathematically well posed or ill posed. We consider the general model with additional physical viscosities and/or additional virtual mass forces, both of which stabilize the system. Two numerical methods, one of them is first-order accurate and the other fifth-order accurate, are used to solve the models. A change of variable technique effectively handles the changing domain and boundary conditions. The numerical methods are demonstrated to be stable and convergent through careful numerical experiments. Simulation results for realistic pulsed liquid fluidized bed are provided and compared with experimental data. (C) 2004 Elsevier Ltd. All rights reserved.

Diffusion dominated process for the crystal growth of a binary alloy

The pure diffusion process has been often used to study the crystal growth of a binary alloy in the microgravity environment. In the present paper, a geometric parameter, the ratio of the maximum deviation distance of curved solidification and melting interfaces from the plane to the radius of the crystal rod, was adopted as a small parameter, and the analytical solution was obtained based on the perturbation theory. The radial segregation of a diffusion dominated process was obtained for cases of arbitrary Peclet number in a region of finite extension with both a curved solidification interface and a curved melting interface. Two types of boundary conditions at the melting interface were analyzed. Some special cases such as infinite extension in the longitudinal direction and special range of Peclet number were reduced from the general solution and discussed in detail.

Analysis of flexural vibration of viscoelastically damped sandwich plates

The simplified governing equations and corresponding boundary conditions of flexural vibration of viscoelastically damped unsymmetrical sandwich plates are given. The asymptotic solution of the equations is then discussed. If only the first terms of the asymptotic solution of all variables are taken as an approximate solution, the result is identical with that obtained from the Modal Strain Energy (MSE) Method. As more terms of the asymptotic solution are taken, the successive calculations show improved accuracy. With the natural frequencies and the modal loss factors of a damped sandwich plate known, one can calculate the response of the plate to various loads providing a reliable basis for engineering design.

Study on cracked spherical shells

The local-global anatysis method is systematically extended to the fracture analysis of spherical shells. On the basis of the shallow shell theory, which takes into account transverse shear deformations, governing equations for cracked spherical shells expressed in displacement and stress functions f, F and φ are proposed, and then a general solution including Modes, Ⅰ, Ⅱ, Ⅲ for stress-strain fields at crack tip in a spherical shell is obtained, which plays the same role as Williams's expansion in plane elasticity. The numerical results for finite-size spherical shells under different boundary conditions have been obtained. Furthermore, the bulging factors are analyzed with regard to shearing stiffness and an approximate formula is given.

考虑刚度及边界条件的索力精确求解

目前索力的测试大多不考虑结构刚度的影响，这导致测试的结果偏大。本文中不仅考虑了索的刚度，而且将索的边界条件假定为简支以及简支和固支耦合两种情况，得到不同假设下的理论公式，分析了不同方法下索力的区别，并利用实测的数据进行分析，给出了末考虑刚度情况下和仅考虑为简支时的误差，证明在索力测试中，特别是短索的测试中，特别是短索的测试中刚度以及边界条件考虑的必要性。

Statistical Simulation of Thin-Film Bearings

The information preservation (IP) method and the direct simulation Monte Carlo (DSMC) method are used to simulate the gas flows between the write/read head and the platter of the disk drive (the slider bearing problem). The results of both methods are in good agreement with numerical solution of the Reynolds equation in the cases studied. However, the DSMC method owing to the problem of large sample size demand and the difficulty in regulating boundary conditions at the inlet and outlet was able to simulate only short bearings, while IP simulates the bearing of authentic length ~1000 m ? and can provide more detailed flow information.

Transient Flows in Low Permeability Porous Media with Composite Start-Up Pressure Gradient

A new mathematical model for the transient flow in the composite low permeability is established. It is solved by FEM with different boundary conditions such as infinite, circular closed and constant pressure boundary conditions. The typical curves for transient wellbore pressure have been presented. It is shown that the pressure and pressure derivative curves with composite start-up pressure gradients have different slopes which are depended on the start-up pressure gradients and the mobility radios in different regions. The boundary effects are the same as the normal reservoirs without start-up pressure gradients. The study provides a new tool to analyze the transient pressure test data in the low permeability reservoir.

Edge and Size Effects in Micro Composite Structure

Because of the load transfer effect of interface layer, the stress distribution inside the composite structure of film/substrate can be very different from the Timoshenko's model. In this paper, we give the derivation and analysis of such load transfer effect of shear-lag (S-L) model. The micro-structure size (boundary conditions) effect together with interface load transfer effect becomes more and more important as the microstructure size including the three dimensions of thickness, width and length shrinks. The microstructure size is also responsible for the so-called edge-induced stress. The edge effect and difference of S-L model and Timoshenko model are also demonstrated.

Information Preservation (IP) Method in Simulation

This paper reviews firstly methods for treating low speed rarefied gas flows: the linearised Boltzmann equation, the Lattice Boltzmann method (LBM), the Navier-Stokes equation plus slip boundary conditions and the DSMC method, and discusses the difficulties in simulating low speed transitional MEMS flows, especially the internal flows. In particular, the present version of the LBM is shown unfeasible for simulation of MEMS flow in transitional regime. The information preservation (IP) method overcomes the difficulty of the statistical simulation caused by the small information to noise ratio for low speed flows by preserving the average information of the enormous number of molecules a simulated molecule represents. A kind of validation of the method is given in this paper. The specificities of the internal flows in MEMS, i.e. the low speed and the large length to width ratio, result in the problem of elliptic nature of the necessity to regulate the inlet and outlet boundary conditions that influence each other. Through the example of the IP calculation of the microchannel (thousands long) flow it is shown that the adoption of the conservative scheme of the mass conservation equation and the super relaxation method resolves this problem successfully. With employment of the same measures the IP method solves the thin film air bearing problem in transitional regime for authentic hard disc write/read head length ( ) and provides pressure distribution in full agreement with the generalized Reynolds equation, while before this the DSMC check of the validity of the Reynolds equation was done only for short ( ) drive head. The author suggests degenerate the Reynolds equation to solve the microchannel flow problem in transitional regime, thus provides a means with merit of strict kinetic theory for testing various methods intending to treat the internal MEMS flows.

Information preservation (IP) method in simulation of internal rarefied gas flows in MEMS

This paper reviews firstly methods for treating low speed rarefied gas flows: the linearised Boltzmann equation, the Lattice Boltzmann method (LBM), the Navier-Stokes equation plus slip boundary conditions and the DSMC method, and discusses the difficulties in simulating low speed transitional MEMS flows, especially the internal flows. In particular, the present version of the LBM is shown unfeasible for simulation of MEMS flow in transitional regime. The information preservation (IP) method overcomes the difficulty of the statistical simulation caused by the small information to noise ratio for low speed flows by preserving the average information of the enormous number of molecules a simulated molecule represents. A kind of validation of the method is given in this paper. The specificities of the internal flows in MEMS, i.e. the low speed and the large length to width ratio, result in the problem of elliptic nature of the necessity to regulate the inlet and outlet boundary conditions that influence each other. Through the example of the IP calculation of the microchannel (thousands m ? long) flow it is shown that the adoption of the conservative scheme of the mass conservation equation and the super relaxation method resolves this problem successfully. With employment of the same measures the IP method solves the thin film air bearing problem in transitional regime for authentic hard disc write/read head length ( 1000 L m ? = ) and provides pressure distribution in full agreement with the generalized Reynolds equation, while before this the DSMC check of the validity of the Reynolds equation was done only for short ( 5 L m ? = ) drive head. The author suggests degenerate the Reynolds equation to solve the microchannel flow problem in transitional regime, thus provides a means with merit of strict kinetic theory for testing various methods intending to treat the internal MEMS flows.

Computation of turbulence-generated noise by large-eddy simulation

The hybrid method of large eddy simulation (LES) and the Lighthill analogy is being developed to compute the sound radiated from turbulent flows. The results obtained from the hybrid method are often contaminated by the absence of small scales in LES, since the energy level of sound is much smaller than that of turbulent flows. Previous researches investigate the effects of subgrid sacle (SGS) eddies on the frequency spectra of sound radiated by isotropic turbulence and suggest a SGS noise model to represent the SGS contributions to the frequency spectra. Their investigations are conducted in physical space and are unavoidably influenced by boundary conditions. In this paper, we propose to perform such calculations in Fourier space so that the effects of boundary conditions can be correctly treated. Posteriori tests are carried out to investigate the SGS contribution to the sound. The results obtained recover the -7/2 law within certain wave-number ranges, but under-estimate the amplitudes of the frequency spectra. The reason for the underestimation is also discussed.