Elasto-viscoplastic-viscous theories of granular and porous media

The mechanics of granular and porous media is considered in the light of the modern theories of structured continuum. The basic laws of motion are presented and several constitutive relations are derived. The special case of elastic porous media is considered in detail and the basic field equations are derived and the possible application of the results to soil dynamics is pointed out. The theory of the flow of granular media is also considered and basic equations of motion are derived where the results of Goodman and Cowin are recovered. The viscoplastic flow of porous media is studied and the possible application to soil mechanics is also discussed.     

Dynamic stress intensity factors around two rectangular cracks in an infinite elastic plate under impact load

A dynamic problem for two equal rectangular cracks in an infinite elastic plate is considered. The two cracks are placed perpendicular to the plane surfaces of the plate. An incoming shock tensile stress is returned by the cracks. In the Laplace transform domain, the boundary conditions at the two sides of the plate are satisfied using the Fourier transform technique. The mixed boundary conditions are reduced to dual integral equations. Crack displacement is expanded in a series of functions which are zero outside of the cracks. The unknown coefficients in the series are determined by the Schmidt method. The stress intensity factors are defined in the Laplace transform domain and these are inverted using a numerical method.     

Experimental and numerical investigation of a viscoplastic Carbopol gel injected into a prototype 3D mold cavity

A numerical model for free-surface flow of a viscoplastic liquid into a cavity is presented. This flow is regarded as a basic model of injection molding, which is a widely used processing technology. Model experiments of the injection process are performed with a water-based gel with shear-thinning behavior. The filling process is visualized by tracing the free surface of the gel within the cavity. Filling times of the cavity are deduced from the experimental observations. The filling process is also analyzed by means of numerical simulation.The flow equations are integrated according to the finite-volume method. The volume-of-fluid method is employed in order to describe the flow of two incompressible, immiscible phases, the phase interface is resolved by the method of geometric reconstruction or alternatively by the method of surface compression. The Herschel–Bulkley model is used in order to describe the shear-thinning behavior of the gel and the effects of a yielding point. The governing equations of the flow are solved by means of the commercial code Fluent as well as the Open Source code OpenFOAM.The results of the numerical simulations are analyzed in detail. They are compared with the experimental findings. Cavity filling times in the experiments and the simulations are in good agreement. Different patterns of the filling flow depending on the injection parameters are evident in the experiments and the simulations. They are characterized and arranged with respect to the similarity parameters of the flow configuration. Again, the results of the simulation are found to agree well with the experimental observations.     

渭河盆地地裂缝发育基本特征

渭河盆地由于其特殊的地质环境,是地裂缝发育较多的地区。本文在渭河盆地地裂缝详细调查的基础上,简要论述了地裂缝的分布概况,详细研究了渭河盆地地裂缝的发育基本特征,并初步探讨了地裂缝的成因。渭河盆地196条地裂缝大都密集分布在断层近侧,与断层走向有明显的一致性和相关性。地裂缝基本特征主要是形态特征和空间分布。地裂缝平面形态有直线、S型和锯齿状;剖面形态为上宽下窄,并逐渐消失;倾角较陡。地裂缝空间分布上有方向性、开启性和网络性。渭河盆地地裂缝主要是构造地裂缝,受构造断裂控制;地下水、地表水、采空和黄土湿陷性只是诱导因素。     

High accuracy non-equidistant method for singular perturbation reaction-diffusion problem

Singular perturbation reaction-diffusion problem with Dirichlet boundary condition is considered. It is a multi-scale problem. Presence of small parameter leads to boundary layer phenomena in both sides of the region. A non-equidistant finite difference method is presented according to the property of boundary layer. The region is divided into an inner boundary layer region and an outer boundary layer region according to transition point of Shishkin. The steps sizes are equidistant in the outer boundary layer region. The step sizes are gradually increased in the inner boundary layer region such that half of the step sizes are different from each other. Truncation error is estimated. The proposed method is stable and uniformly convergent with the order higher than 2. Numerical results are given, which are in agreement with the theoretical result.     

Comparison of waveguide properties of curved versus straight planar elastic layers

Dispersion equations are solved for the in-plane and anti-plane wave propagation in planar elastic layer with constant curvature. The classical Lamé formulation of displacements via elastic potentials is applied and appropriate simplifications are employed. The dispersion diagrams in each case are compared with their counterparts for a straight layer, e.g., the classical Rayleigh–Lamb solution. The curvature-induced symmetry-breaking effects are investigated for layers with symmetric boundary conditions. The role of curvature is also investigated in the cases, when the boundary conditions are not symmetrical. The elementary Bernoulli–Euler theory is employed to analyze the wave guide properties of a curved planar elastic beam in its in-plane deformation. The validity range of the Bernoulli–Euler theory is assessed via comparison of dispersion diagrams.     

Stability of the in-plane bending mode of thin-walled framed structures

The paper outlines an approach to solving the stability problem for framed structures under arbitrary transverse loading. The available methods are limited by one law of variation in the bending moment responsible for loss of stability. The equilibrium equations for a thin-walled bar are integrated assuming that the bending moment is constant. The solution of the Cauchy problem is given in normal form. The arbitrary varying bending moment is approximated by a piecewise-constant function, which will be a little different from the original if the bar is partitioned into a great number of segments. The equations of the boundary-value problem for a discretized framed structure are derived using the boundary-element method. The critical forces and moments are determined from a transcendent equation. Numerical solutions are presented to demonstrate the high accuracy and efficiency of the approach. The solutions of test problems are in agreement with those obtained by Timoshenko     

Development and application of radiowave photoelasticity for the fracture analysis of dielectrics

The radiowave photoelasticity method is developed and is justified for use in the fracture mechanics of dielectrics. The damage tensors of a material with microcracks are considered. The basic equations of the method are given. The time-dependence of microcrack density is analyzed. Numerical results are presented and analyzed     

UNCONDITIONAL STABLE SOLUTIONS OF THE EULER EQUATIONS FOR TWO－AND THREE－D WINGS IN ARBITRARY MOTION

ＵＮＣＯＮＤＩＴＩＯＮＡＬＳＴＡＢＬＥＳＯＬＵＴＩＯＮＳＯＦＴＨＥＥＵＬＥＲＥＱＵＡＴＩＯＮＳＦＯＲＴＷＯ－ＡＮＤＴＨＲＥＥ－ＤＷＩＮＧＳＩＮＡＲＢＩＴＲＡＲＹＭＯＴＩＯＮＧａｏＺｈｅｎｇｈｏｎｇ（高正红）（ＲｅｃｅｉｖｅｄＪａｎ．１２,１９９５,Ｃ...     

Two degree-of-freedom oscillator coupled to a non-ideal source

In the paper the one-mass two degree-of-freedom system with non-ideal excitation is considered. The resonance motion of the system is investigated. The mathematical model of the system contains three coupled second order differential equations. In the paper an analytical solving procedure is developed. The steady-state motion and the criteria for stability of solutions are developed. Two special cases of motion depending on the frequency properties of the system are studied. When the frequency properties in both orthogonal direction are equal there is only one resonance. If the frequency in one direction is two times higher than in other two different resonances occur: one in x and the other in y direction. The conditions for jump phenomena and for Sommerfeld effect are presented. The analytically obtained solutions are compared with numerical ones. They show good agreement.     

Long-wave instability of a plane multicomponent mixture layer with the Soret effect

Long-wave instability that arises in a plane multicomponent mixture layer with the Soret effect is investigated. The layer is heated from above or below in the gravity field. The solution of the linearized equations of motion is sought in the form of a series in powers of the wave number assumed to be small. For the mixture with an arbitrary number of components the critical parameters of instability are determined. For binary and ternary mixtures the stability maps are plotted. The regions where the long- and short wave perturbations are the most dangerous are distinguished in the space of physical parameters. The analytical results are compared with the numerical solution of the linearized equations over a wide wave number range. As an example, instability in a ternary hydrocarbon mixture is considered.     

Numerical analysis of breaking waves using the moving particle semi-implicit method

The numerical method used in this study is the moving particle semi-implicit (MPS) method, which is based on particles and their interactions. The particle number density is implicitly required to be constant to satisfy incompressibility. A semi-implicit algorithm is used for two-dimensional incompressible non-viscous flow analysis. The particles whose particle number densities are below a set point are considered as on the free surface. Grids are not necessary in any calculation steps. It is estimated that most of computation time is used in generation of the list of neighboring particles in a large problem. An algorithm to enhance the computation speed is proposed. The MPS method is applied to numerical simulation of breaking waves on slopes. Two types of breaking waves, plunging and spilling breakers, are observed in the calculation results. The breaker types are classified by using the minimum angular momentum at the wave front. The surf similarity parameter which separates the types agrees well with references. Breaking waves are also calculated with a passively moving float which is modelled by particles. Artificial friction due to the disturbed motion of particles causes errors in the flow velocity distribution which is shown in comparison with the theoretical solution of a cnoidal wave. © 1998 John Wiley & Sons, Ltd.     

Effect of vertical stacking dies on flow behavior of epoxy molding compound during encapsulation of stacked-chip scale packages

This paper presents three dimensional (3D) simulation of flow visualization in the encapsulation of stacked-chip scales packages (S-CSP), using finite volume method. The S-CSP model is constructed using GAMBIT and simulated using FLUENT CFD software. The epoxy molding compound is Hitachi CEL-9200 XU (LF) and its flow is assumed laminar and incompressible. Cross viscosity model and volume of fluid technique are applied for flow front tracking of the encapsulant. The meshing is performed using tetrahedral elements and the discretization is done by first order upwind scheme. SIMPLE algorithm is selected for solving the governing equations. The top view and 3D view of simulation flow front profiles in the encapsulation process are presented. The percentage of filled volume versus filling time, viscosity versus shear rate and number of voids versus rows of stacked die are plotted. The temperature and pressure distributions within the mold cavity during the encapsulation process are also observed. Further, the possibility and cause of void formation during the encapsulation process are analyzed and discussed in detail. The number of vertical stacking dies and horizontal rows of packages are found to be crucial in the void formation. The numerical results are compared with previous experimental results and found in good conformity.     

Slip and induced magnetic field effects on peristaltic transport of Johnson-Segalman fluid

The peristaltic flow of a Johnson-Segalman fluid in a planar channel is investigated in an induced magnetic field with the slip condition.The symmetric nature of the flow in a channel is utilized.The velocity slip condition in terms of shear stresses is considered.The mathematical formulation is presented,and the equations are solved under long wavelength and low Reynolds number approximations.The perturbation solutions are established for the pressure,the axial velocity,the micro-rotation component,the stream function,the magnetic-force function,the axial induced magnetic field,and the current distribution across the channel.The solution expressions for small Weissenberg numbers are derived.The flow quantities of interest are sketched and analyzed.     

Interphase momentum interaction effects in the averaged multifield model: Part I: Void propagation in bubbly flows

The volume-averaged form of the linear momentum conservation equations for twophase flow is examined to clarify momentum interaction effects between phases. The case of an accelerating sphere of varying radius in an accelerating fluid is used to derive the form of the interphase force terms. The analysis is extended to assemblies of noninteracting spheres and an interphase force term related to the spatial gradient of phase volume fraction is seen to arise. For the case of bubbly flow, two real characteristics are obtained for dispersed phase volume fractions less than about 0.25. If the term involving the spatial gradient of the phase volume fraction is neglected, then the characteristics are always complex for velocity differences between the phases. The interphase force model is applied to predict experiments on void propagation in bubbly flows. There are no adjustable constants in the model. The experimental data were obtained in our laboratory using cross correlation of signals from a pair of gamma densitometers. The predictions are in excellent agreement with the data. In addition, the predictions are compared with data from several other laboratories, taken over different sets of flow conditions. The predictions are again in close agreement with the data.     

Transient response of a moving spherical shell in an acoustic medium

A ring-stiffened spherical shell is submerged in an acoustic medium. The shell is thin and elastic. The acoustic medium is inviscid, irrotational and compressible. The center of mass of the shell is subjected to a translational acceleration which is an arbitrary function of time. The absolute displacements of the shell are expressed in terms of the relative displacements and the displacement of the base of the shell, base being defined as the rigid ring placed at the equator. The motion of the acoustic medium is governed by the wave equation. The transient response of the shell is investigated numerically. The results are compared with the results of the in-vacuo response. The effects of the plane wave approximation and the base velocity on the transient response of the shell are studied. The numerical results show that the plane wave approximation accurately predicts the response of the shell in the acoustic medium for short times after excitation. The displacements of the shell in fluid are larger than those in vacuo. But when the base of the shell is restrained from translating, the displacements in fluid are smaller than those in vacuo. Therefore, base translation has a very significant effect on the transient response of the shells submerged in an acoustic medium.     

Numerical analysis of nonlinear modes of piecewise linear systems torsional vibrations

The nonlinear modes of essentially nonlinear piecewise-linear finite degrees of freedom systems are calculated by the numerical methods, which are suggested in this paper. The basis of these methods is numerical solutions of the equations of the systems motions in configuration space. The numerical method for the nonlinear modes of essentially nonlinear piecewise-linear systems forced vibrations is suggested. The basis of this approach is the combination of the Rauscher method and the calculations of the autonomous system nonlinear modes. The nonlinear modes of the diesel engine transmission torsional vibrations are analyzed numerically. The vibrations are described by essentially nonlinear piecewise-linear system with fifteen degrees of freedom. The NNMs of this system forced vibrations are observed in the resonance regions. Both NNMs and the motions, which are essentially differ from NNMs, are observed in the distance from the resonances. NNMs of the forced vibrations of the systems with dissipation are close to NNMs of the system without dissipation.     

Holographic moiré in real time

The body of knowledge necessary to observe holographic-moiré patterns in real time is introduced. The basic factors influencing fringe visibility in holographic moiré are analyzed and expressions to evaluate fringe visibility for any given displacement and deformation are given. The application of the introduced theory in the case of real-time observation is discussed. It is shown that the maximum benefits of this technique are achieved by combining it with closed-circuit TV. Several examples of application are given.     

Onset of Thermal Convection in a Vertical Porous Cylinder with a Partly Conducting and Partly Penetrative Cylinder Wall

The onset of thermal convection in a vertical porous cylinder in three dimensions is investigated analytically. Top and bottom of the cylinder are set to be perfectly heat conducting and impermeable, and is uniformly heated from below. The convection problem is solved for a cylinder wall that is partly conducting and partly penetrative. The expressions for semi-conduction and semi-penetration are based on a porous medium separated from its surroundings by a thin wall. The eigenvalue problem is split into two Helmholtz equations, and the results are expressed by Bessel functions in the radial direction. Comparisons are made with existing solutions for the limit cases of a closed cylinder wall that is either conducting or insulating. Two different models are compared for the kinematic limit condition of an open boundary.