GENERAL FORM OF MATCHING EQUATION OF ELASTIC-PLASTIC FIELD NEAR CRACK LINE FOR MODEⅠCRACK UNDER PLANE STRESS CONDITION

Crack line field analysis method has become an independent method for crack elastic-plastic analysis, which greatly simplifies the complexity of crack elastic-plastic problems and overcomes the corresponding mathematical difficulty. With this method, the precise elastic-plastic solutions near crack lines for variety of crack problems can be obtained. But up to now all solutions obtained by this method were for different concrete problems, no general steps and no general form of matching equations near crack line are given out. With crack line analysis method, this paper proposes the general steps of elastic-plastic analysis near crack line for mode Ⅰ crack in elastic-perfectly plastic solids under plane stress condition, and in turn given out the solving process and result for a specific problem.     

Optimal Elastic Design for Prescribed Maximum Deflection

ABSTRACT

For the problems of the optimal elastic design with prescribed maximum deflection, a variational formulation is proposed, with reference to the one-or two-dimensional bending structures.

Necessary optimum conditions are found, and the physical features of the optimal solutions are discussed for the “absolute” minimum cost problems, and, when dealing with beams, for the solutions with piece-wise constant design function.

Some examples are solved by using numerical methods that are directly derived from the variational formulation.     

Edge impact of an elastic-plastic semi-infinite plate

The edge impact of an elastic-plastic semiinfinite plate subject to conditions of plane strain is investigated analytically and experimentally. The theoretical analysis is based on the strain-rate-independent theory of plastic-wave propagation. The plate is initially unstressed; the boundary condition for the edge of the plate corresponds to constant-velocity longitudinal impact except that the step in velocity has a finite rise time. Calculations are carried out according to both the approximate one-dimensional theory and the two-dimensional theory for an elastic-plastic isotropic work-hardening material. The rise time for both solutions is chosen so as to optimize the agreement between theoretical and experimental strain-time profiles. A numerical solution of the two-dimensional equations is obtained by using a difference method developed by Clifton;¹ the onedimensional approximation is solved by the well-known method of characteristics. The problem was approximated experimentally by the axial collision of 4-in.-diam annealed aluminum thick-walled cylinders with a diameter-to-wall-thickness ratio of ten. For impact velocities of 90, 130 and 160 ips (corresponding to maximum strains of 0.12, 0.22 and 0.36 percent, respectively), dispersive characteristics and maximum strain amplitudes of the strain wave are found to be in good agreement with the theoretical predictions of both solutions. However, the two-dimensional solution indicates that the stresses, strains and velocities in regions of high strain-rate are highly nonuniform across the plate thickness.     

Numerical simulation of flow fields induced by a supersonic projectile moving in tubes

Numerical study of the shocked flows generated by a supersonic projectile released from a launch tube into a big chamber has been performed in this paper. Based on fixed Cartesian grids, the two-dimensional axisymmetric Euler equations are solved by the fifth-order WENO scheme implemented with moving boundary conditions. Both the level set technique and ghost fluid method are used for capturing the moving interface of the projectile implicitly. The numerical results show that complex shock phenomena exist in the transient shock flow, resulting from shock-wave reflection, shock-wave focusing, shock-wave/projectile interaction and shock-wave/contact surface interactions. The relationships between the acceleration of the projectile and the transient shock flow are also discussed in detail.     

Shock-wave equation-of-state studies at Los Alamos

A history of the shock-wave equation-of-state (EOS) studies at Los Alamos is given. Particular emphasis is placed on the pioneering research in the 1950s where many of the experimental techniques and methods of analysis were developed, which we now take for granted. A brief review of shock-wave physics is given, which illustrates important hydrodynamic and thermodynamic concepts. Recent studies on the EOS of Ti are presented with emphasis on the-to- phase transition. VISAR wave profiles for polycrystalline Ni and singlecrystal Ni are presented to determine the strengths of these materials under pressure. Low-density polystyrene foam Hugoniot experiments are described and results analyzed.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1990.     

Termination of the starting problem of dynamic expansion of a spherical cavity in an infinite elastic-perfectly-plastic medium

Perhaps the simplest non-trivial problem in small deformation dynamic plasticity is expansion of a spherical cavity in an infinite elastic-perfectly-plastic medium. Here, example problems are considered with two boundary conditions at the cavity's surface: constant velocity and constant pressure. Attempts to obtain analytical solutions are complicated by the fact that, in general, the elastic-plastic boundary propagates with variable speed. However, it is known that the elastic-plastic boundary propagates at constant speed for the starting problem when the shocks due to the applied loads are large enough to cause inelastic response at the instant they are applied. When the value of the applied pressure equals the shock pressure due to the applied velocity the solutions of the two boundary value problems are initially identical and can be compared. The objective of this paper is to review the literature and to examine the termination conditions for the starting problem. Specifically, the starting problem terminates when either the jump in radial stress at the elastic-plastic boundary or the loading condition for plasticity vanishes there. These termination conditions depend on the applied load and on Poisson's ratio.     

Hydrodynamics of the interaction between bodies in a viscous fluid in the vicinity of their contact at low and high Reynolds numbers

An interaction between bodies in the vicinity of their contact is investigated. Stresses determining the main part of the forces acting on the bodies in motion are produced in a narrow gap between the bodies in the vicinity of their contact. In many cases the velocity and pressure fields in the vicinity of the contact can be determined and the main asymptotics for the hydrodynamic interaction force can be constructed in the small distance between the surfaces. An overview of the problems solved using this approach is presented. For certain problems new formulations are given. The plausibility of the results is confirmed by comparing with the available exact particular solutions and the experimental data. Owing to the restrictions on the size of the paper, only two-dimensional problems are considered, although the approach developed can be applied to the solution of three-dimensional problems as well.     

Bounds to bifurcation stresses in solids with non-associated plastic flow law at finite strain

In the present paper, Hill's theory of bifurcation and stability in solids obeying normality is generalized to include a non-associated flow law. A one-parameter family of linear comparison solids has been found that admits a potential and has the property that if uniqueness is certain for the comparison solid then bifurcation and instability are precluded for the underlying elastic-plastic solid. The uniqueness criterion derived may be used as a device to determine lower bounds to the magnitudes of primary bifurcation and instability stresses which are ordinarily unknown. A second linear solid is introduced whose constitutive relations have the same form as the elastic-plastic solid “in loading”. The first eigenstate of this solid gives an upper bound to the primary bifurcation state of the underlying elastic-plastic solid. The search for the genuine primary bifurcation state is therefore replaced by a search for upper and lower bounds in the situation when normality fails to hold. The theory is applied to problems of homogeneous stress states.     

A marker-and-cell approach to viscoelastic free surface flows using the PTT model

This work is concerned with the development of a numerical method capable of simulating two-dimensional viscoelastic free surface flows governed by the non-linear constitutive equation PTT (Phan-Thien–Tanner). In particular, we are interested in flows possessing moving free surfaces. The fluid is modelled by a marker-and-cell type method and employs an accurate representation of the fluid surface. Boundary conditions are described in detail and the full free surface stress conditions are considered. The PTT equation is solved by a high order method which requires the calculation of the extra-stress tensor on the mesh contour. The equations describing the numerical technique are solved by the finite difference method on a staggered grid. In order to validate the numerical method fully developed flow in a two-dimensional channel was simulated and the numerical solutions were compared with known analytic solutions. Convergence results were obtained throughout by using mesh refinement. To demonstrate that complex free surface flows using the PTT model can be computed, extrudate swell and a jet flowing onto a rigid plate were simulated.     

General forms of elastic-plastic matching equations for mode-Ⅲ cracks near crack line

Crack line analysis is an effective way to solve elastic-plastic crack problems.Application of the method does not need the traditional small-scale yielding conditions and can obtain sufficiently accurate solutions near the crack line. To address mode-Ⅲ crack problems under the perfect elastic-plastic condition,matching procedures of the crack line analysis method are summarized and refined to give general forms and formulation steps of plastic field,elastic-plastic boundary,and elastic-plastic matching equations near the crack line. The research unifies mode-Ⅲ crack problems under different conditions into a problem of determining four integral constants with four matching equations.An example is given to verify correctness,conciseness,and generality of the procedure.     

General forms of elastic-plastic matching equations for mode-Ⅲ cracks near crack line

Crack line analysis is an effective way to solve elastic-plastic crack problems. Application of the method does not need the traditional small-scale yielding conditions and can obtain sufficiently accurate solutions near the crack line. To address mode- Ⅲ crack problems under the perfect elastic-plastic condition, matching procedures of the crack line analysis method axe summarized and refined to give general forms and formulation steps of plastic field, elastic-plastic boundary, and elastic-plastic matching equations near the crack line. The research unifies mode-III crack problems under different conditions into a problem of determining four integral constants with four matching equations. An example is given to verify correctness, conciseness, and generality of the procedure.     

饱和多孔介质分析解的唯一性与应变局部化分岔

基于不连续性分岔基本理论导出了静态非渗流状态下弹塑性饱和多孔介质应变局部化发生时的临界硬/软化模量,利用二阶功正定性原理研究了两相问题分析解的唯一性问题,并给出了基于主轴空间下解的显式表达式。研究工作表明,在静态非渗流状态下,弹塑性饱和多孔介质分析的唯一性与应变局部化发生的临界条件除了在量值上与单相介质有着明显的不同外两者之间还有许多一致的特性,这些一致的特性对问题的分析是十分重要的。     

Investigation of hydrodynamic instability of CO<Subscript>2</Subscript> injection into an aquifer

The two-dimensional problem of supercritical carbon dioxide injection into an aquifer is solved. Shocks and rarefaction waves propagating in a sequence from an injection well into the formation are described within the framework of a complete nonisothermal model of flows in a porous medium. In the approximation of isothermal immiscible water and carbon dioxide flow the hydrodynamic stability of the leading displacement front is investigated for various reservoir pressures and temperatures. The parameters of unstable fronts are determined using a sufficient instability condition formulated in analytic form. The approximate analytic results are supported by the direct numerical simulation of CO₂ injection using the complete model in which thermal effects and phase transitions are taken into account.     

Slow motions of a circular cylinder experiencing slip near a plane wall

A theoretical study is presented for the two-dimensional creeping flow caused by a long circular cylindrical particle translating and rotating in a viscous fluid near a large plane wall parallel to its axis. The fluid is allowed to slip at the surface of the particle. The Stokes equations for the fluid velocity field are solved in the quasi-steady limit using cylindrical bipolar coordinates. Semi-analytical solutions for the drag force and torque acting on the particle by the fluid are obtained for various values of the slip coefficient associated with the particle surface and of the relative separation distance between the particle and the wall. The results indicate that the translation and rotation of the confined cylinder are not coupled with each other. For the motion of a no-slip cylinder near a plane wall, our hydrodynamic drag force and torque results reduce to the closed-form solutions available in the literature. The boundary-corrected drag force and torque acting on the particle decrease with an increase in the slip coefficient for an otherwise specified condition. The plane wall exerts the greatest drag on the particle when its migration occurs normal to it, and the least in the case of motion parallel to it. The enhancement in the hydrodynamic drag force and torque on a translating and rotating particle caused by a nearby plane wall is much more significant for a cylinder than for a sphere.     

On the application of pseudospectral FFT techniques to non-periodic problems

The reduction-to-periodicity method using the pseudospectral fast Fourier transform (FFT) technique is applied to the solution of non-periodic problems, including the two-dimensional incompressible Navier–Stokes equations. The accuracy of the method is explored by calculating the derivatives of given functions, one- and two-dimensional convective-diffusive problems, and by comparing the relative errors due to the FFT method with a second-order finite difference (FD) method. Finally, the two-dimensional Navier–Stokes equations are solved by a fractional step procedure using both the FFT and the FD methods for the driven cavity flow and the backward-facing step problems. Comparisons of these solutions provide a realistic assessment of the FFT method.     

Solution of the plane stress problems of strain-hardening materials described by power-law using the complex pseudo-stress function

In the present paper,the compatibility equation for the plane stress problems of power-law materials is transformed into a biharmonic equation by introducing the so-calledcomplex pseudo-stress function,which makes it possible to solve the elastic-plastic planestress problems of strain hardening materials described by power-law using the complexvariable function method like that in the linear elasticity theory.By using this generalmethod,the close-formed analytical solutions for the stress,strain and displacementcomponents of the plane stress problems’of power-law materials is deduced in the paper,which can also be used to solve the elasto-plastic plane stress problems of strain-hardeningmaterials other than that described by power-law.As an example,the problem of a power-law material infinite plate containing a circular hole under uniaxial tension is solved byusing this method,the results of which are compared with those of a known asymptoticanalytical solution obtained by the perturbation method.     

Path-independentJ-integral and its dual form in elastic-plastic solids with finite displacements

In this paper, based on energy variational principles of elastic-plastic solids, the path-independentJ-integral and its dual form in elastic-plastic solids with finite displacements are presented. Whose testification is given there after.     

Reflection of shock wave from a compression corner in a particle-laden gas region

We investigated in this paper the progression of a shock-wave reflected from a compression corner in a particle-laden gas medium using a TVD class numerical technique and a MacCormack scheme. For a gas-only flow, the numerical results agreed well with the existing experimental data, suggesting that the gas phase is correctively solved. The effect of particle size and mass fraction ratio is investigated for a dilute gas-particle flow. It has been shown that the shock-wave diffraction and the flow configuration after the shock can become remarkably different from the gas-only flow depending on the particle parameters. Relaxation phenomenon due to the momentum drag and the heat exchange between the gas and the particle phases is explained.Graduate Student of Korea Advanced Institute of Science and TechnologyThis article was processed using Springer-Verlag T_EX Shock Waves macro package 1990.     

Constrained moving least‐squares immersed boundary method for fluid‐structure interaction analysis

A numerical method is presented for the analysis of interactions of inviscid and compressible flows with arbitrarily shaped stationary or moving rigid solids. The fluid equations are solved on a fixed rectangular Cartesian grid by using a higher‐order finite difference method based on the fifth‐order WENO scheme. A constrained moving least‐squares sharp interface method is proposed to enforce the Neumann‐type boundary conditions on the fluid‐solid interface by using a penalty term, while the Dirichlet boundary conditions are directly enforced. The solution of the fluid flow and the solid motion equations is advanced in time by staggerly using, respectively, the third‐order Runge‐Kutta and the implicit Newmark integration schemes. The stability and the robustness of the proposed method have been demonstrated by analyzing 5 challenging problems. For these problems, the numerical results have been found to agree well with their analytical and numerical solutions available in the literature. Effects of the support domain size and values assigned to the penalty parameter on the stability and the accuracy of the present method are also discussed.