Publications

Abstract not provided.

Abstract not provided.

Kevlar materials make excellent body armor due to their fabric-like flexibility and ultra-high tensile strength. Carbon composites are made up from many layers of carbon AS-4 material impregnated with epoxy. Fiber orientation is bidirectional, orientated at 0° and 90°. They also have ultra-high tensile strength but can be made into relatively hard armor pieces. Once many layers are cut and assembled they can be ergonomicically shaped in a mold during the heated curing process. Kevlar and carbon composites can be used together to produce light and effective body armor. This paper will focus on computer analysis and laboratory testing of a Kevlar/carbon composite cross-section proposed for body armor development. The carbon composite is inserted between layers of Kevlar. The computer analysis was performed with a Lagrangian transversely Isotropic material model for both the Kevlar and Carbon Composite. The computer code employed is AUTODYN. Both the computer analysis and laboratory testing utilized different fragments sizes of hardened steel impacting on the armor cross-section. The steel fragments are right-circular cylinders. Laboratory testing was undertaken by firing various sizes of hardened steel fragments at square test coupons of Kevlar layers and heat cured carbon composites. The V50 velocity for the various fragment sizes was determined from the testing. This V50 data can be used to compare the body armor design with other previously designed armor systems. AUTODYN [1] computer simulations of the fragment impacts were compared to the experimental results and used to evaluate and guide the overall design process. This paper will include the detailed transversely isotropic computer simulations of the Kevlar/carbon composite cross-section as well as the experimental results and a comparison between the two. Conclusions will be drawn about the design process and the validity of current computer modeling methods for Kevlar and carbon composites. Copyright © 2005 by ASME.

Abstract not provided.

The purpose of the present work is to increase our understanding of which properties of geomaterials most influence the penetration process with a goal of improving our predictive ability. Two primary approaches were followed: development of a realistic, constitutive model for geomaterials and designing an experimental approach to study penetration from the target's point of view. A realistic constitutive model, with parameters based on measurable properties, can be used for sensitivity analysis to determine the properties that are most important in influencing the penetration process. An immense literature exists that is devoted to the problem of predicting penetration into geomaterials or similar man-made materials such as concrete. Various formulations have been developed that use an analytic or more commonly, numerical, solution for the spherical or cylindrical cavity expansion as a sort of Green's function to establish the forces acting on a penetrator. This approach has had considerable success in modeling the behavior of penetrators, both as to path and depth of penetration. However the approach is not well adapted to the problem of understanding what is happening to the material being penetrated. Without a picture of the stress and strain state imposed on the highly deformed target material, it is not easy to determine what properties of the target are important in influencing the penetration process. We developed an experimental arrangement that allows greater control of the deformation than is possible in actual penetrator tests, yet approximates the deformation processes imposed by a penetrator. Using explosive line charges placed in a central borehole, we loaded cylindrical specimens in a manner equivalent to an increment of penetration, allowing the measurement of the associated strains and accelerations and the retrieval of specimens from the more-or-less intact cylinder. Results show clearly that the deformation zone is highly concentrated near the borehole, with almost no damage occurring beyond 1/2 a borehole diameter. This implies penetration is not strongly influenced by anything but the material within a diameter or so of the penetration. For penetrator tests, target size should not matter strongly once target diameters exceed some small multiple of the penetrator diameter. Penetration into jointed rock should not be much affected unless a discontinuity is within a similar range. Accelerations measured at several points along a radius from the borehole are consistent with highly-concentrated damage and energy absorption; At the borehole wall, accelerations were an order of magnitude higher than at 1/2 a diameter, but at the outer surface, 8 diameters away, accelerations were as expected for propagation through an elastic medium. Accelerations measured at the outer surface of the cylinders increased significantly with cure time for the concrete. As strength increased, less damage was observed near the explosively-driven borehole wall consistent with the lower energy absorption expected and observed for stronger concrete. As it is the energy absorbing properties of a target that ultimately stop a penetrator, we believe this may point the way to a more readily determined equivalent of the S number.

Blastwalls are often assumed to be the answer for facility protection from malevolent explosive assault, particularly from large vehicle bombs (LVB's). The assumption is made that the blastwall, if it is built strong enough to survive, will provide substantial protection to facilities and people on the side opposite the LVB. This paper will demonstrate through computer simulations and experimental data the behavior of explosively induced air blasts during interaction with blastwalls. It will be shown that air blasts can effectively wrap around and over blastwalls. Significant pressure reduction can be expected on the downstream side of the blastwall but substantial pressure will continue to propagate. The effectiveness of the blastwall to reduce blast overpressure depends on the geometry of the blastwall and the location of the explosive relative to the blastwall.

Abstract not provided.

The authors understanding of multiphase physics and the associated predictive capability for multi-phase systems are severely limited by current continuum modeling methods and experimental approaches. This research will deliver an unprecedented modeling capability to directly simulate three-dimensional multi-phase systems at the particle-scale. The model solves the fully coupled equations of motion governing the fluid phase and the individual particles comprising the solid phase using a newly discovered, highly efficient coupled numerical method based on the discrete-element method and the Lattice-Boltzmann method. A massively parallel implementation will enable the solution of large, physically realistic systems.

Abstract not provided.

Abstract not provided.

Processes that involve particle-laden fluids are common in geomechanics and especially in the petroleum industry. Understanding the physics of these processes and the ability to predict their behavior requires the development of coupled fluid-flow and particle-motion computational methods. This paper outlines an accurate and robust coupled computational scheme using the lattice-Boltzmann method for fluid flow and the discrete-element method for solid particle motion. Results from several two-dimensional validation simulations are presented. Simulations reported include the sedimentation of an ellipse, a disc and two interacting discs in a closed column of fluid. The recently discovered phenomenon of drafting, kissing, and tumbling is fully reproduced in the two-disc simulation.

Techniques for modeling oil well sand production have been developed using the formulations for superquadric discrete elements and Darcy fluid flow. Discrete element models are generated using the new technique of particle cloning. Discrete element sources and sinks allow simulation of sand production from the initial state through the transition to an equilibrium state where particles are created and removed at the same rate.

A discrete element computer program named DMC_BLAST (Distinct Motion Code) has been under development since 1987 for modeling rock blasting (Preece & Taylor, 1989). This program employs explicit time integration and uses spherical or cylindrical elements that are represented as circles in two dimensions. DMC_BLAST calculations compare favorably with data from actual bench blasts (Preece et al, 1993). Coal seam chilling refers to the shattering of a significant portion of the coal leaving unusable fines. It is also refereed to as coal damage. Chilling is caused during a blast by a combination of explosive shock energy and movement of the adjacent rock. Chilling can be minimized by leaving a buffer zone between the bottom of the blastholes and the coal seam or by changing the blast design to decrease the powder factor or by a combination of both. Blast design in coal mine cast blasting is usually a compromise between coal damage and rock fragmentation and movement (heave). In this paper the damage to coal seams from rock movement is examined using the discrete element computer code DMC_BLAST. A rock material strength option has been incorporated into DMC_BLAST by placing bonds/links between the spherical particles used to model the rock. These bonds tie the particles together but can be broken when the tensile, compressive or shear stress in the bond exceeds the defined strength. This capability has been applied to predict coal seam damage, particularly at the toe of a cast blast where drag forces exerted by movement of the overlying rock can adversely effect the top of the coal at the bench face. A simulation of coal mine cast blasting has been performed with special attention being paid to the strength of the coal and its behavior at t he bench face during movement of the overlying material.

Coupled blast-structural computational simulations using supercomputer capabilities will significantly advance the understanding of how complex structures respond under dynamic loads caused by explosives and earthquakes, an understanding with application to the surety of both federal and nonfederal buildings. Simulation of the effects of explosives on structures is a challenge because the explosive response can best be simulated using Eulerian computational techniques and structural behavior is best modeled using Lagrangian methods. Due to the different methodologies of the two computational techniques and code architecture requirements, they are usually implemented in different computer programs. Explosive and structure modeling in two different codes make it difficult or next to impossible to do coupled explosive/structure interaction simulations. Sandia National Laboratories has developed two techniques for solving this problem. The first is called Smoothed Particle Hydrodynamics (SPH), a relatively new gridless method comparable to Eulerian, that is especially suited for treating liquids and gases such as those produced by an explosive. The SPH capability has been fully implemented into the transient dynamics finite element (Lagrangian) codes PRONTO-2D and -3D. A PRONTO-3D/SPH simulation of the effect of a blast on a protective-wall barrier is presented in this paper. The second technique employed at Sandia National Laboratories uses a relatively new code called ALEGRA which is an ALE (Arbitrary Lagrangian-Eulerian) wave code with specific emphasis on large deformation and shock propagation. ALEGRA is capable of solving many shock-wave physics problems but it is especially suited for modeling problems involving the interaction of decoupled explosives with structures.

A discrete element computer program named DMC{underscore}BLAST (Distinct Motion Code) has been under development since 1987 for modeling rock blasting (Preece {ampersand} Taylor, 1989). This program employs explicit time integration and uses spherical or cylindrical elements that are represented as circles in two dimensions (2-D). DMC{underscore}BLAST calculations compare favorably with data from actual bench blasts (Preece et al, 1993). Buffer Choke blasting is commonly used in surface gold mines to break the rock and dilate it sufficiently for ease of digging, with the assumption of insignificant horizontal movement. The blast designs usually call for relatively shallow holes benches ({lt} 11 m) with small blastholes (approx. 165 mm), small burdens and spacings ({lt}5 m), often with 50% or more of the hole stemmed. Control of blast-induced horizontal movement is desired because the ore is assayed in place from the blasthole drill cuttings and digging polygons of ore and waste are laid out before the blast. Horizontal movement at the ore waste boundary can result in dilution of the ore or loss of ore with the waste. The discrete element computer program DMC{underscore}BLAST has been employed to study spatial variation of horizontal rock motion during buffer choke blasting. Patterns of rock motion can be recognized from the discrete element simulations that would be difficult or impossible to recognize in the field (Preece, Tidman and Chung, 1997). Techniques have been developed to calculate ore waste and dilution from the horizontal movement predicted by DMC{underscore}BLAST. Four DMC{underscore}BLAST simulations of buffer blasting have been performed. The blasts are identical except that the burden and spacing are systematically varied which also changes the powder factor. Predictions of ore waste or dilution are made for each burden in the blast, assuming no horizontal movement, to illustrate the spatial variation observed.

Simulation of the effects of explosives on structures is a challenge because the explosive response can best be simulated using Eulerian computational techniques and structural behavior is best modeled using Lagrangian methods. Due to the different methodology of the two computational techniques and code architecture requirements, they are usually implemented in different computer programs. Explosive and structure modeling in two different codes make it difficult or next to impossible to do coupled explosive/structure interaction simulations. Sandia National Laboratories has developed two techniques for solving this problem. The first is called Smoothed Particle Hydrodynamics (SPH), a relatively new gridless method comparable to Eulerian, that is especially suited for treating liquids and gases such as those produced by an explosive. The SPH capability has been fully implemented into the transient dynamics finite element (Lagrangian) codes PRONTO-2D and -3D. A PRONTO-3D/SPH simulation of the effect of a blast on a protective-wall barrier is presented in this paper. The second technique employed at Sandia uses a new code called Zapotec that combines the 3-D Eulerian code CTH and the Lagrangian code PRONTO-3D with minimal changes to either code. CTH and PRONTO-3D are currently executing on the Sandia Terraflops machine (9000 Pentium Pro processors). Eulerian simulations with 100 million cells have been completed on the current configuration of the machine (4500 Pentium Pro processors). The CTH and PRONTO-3D combination will soon be executing in a coupled fashion on this machine.

A discrete element computer program named DMC{_}BLAST (Distinct Motion Code) has been under development since 1987 for modeling rock blasting. This program employs explicit time integration and uses spherical or cylindrical elements that are represented as circles in 2-D. DMC{_}BLAST calculations compare favorably with data from actual bench blasts. The blast modeling capabilities of DMC{_}BLAST have been expanded to include independently dipping geologic layers, top surface, bottom surface and pit floor. The pit can also now be defined using coordinates based on the toe of the bench. A method for modeling decked explosives has been developed which allows accurate treatment of the inert materials (stemming) in the explosive column and approximate treatment of different explosives in the same blasthole. A DMC{_}BLAST user can specify decking through a specific geologic layer with either inert material or a different explosive. Another new feature of DMC{_}BLAST is specification of an uplift angle which is the angle between the normal to the blasthole and a vector defining the direction of explosive loading on particles adjacent to the blasthole. A buffer (choke) blast capability has been added for situations where previously blasted material is adjacent to the free face of the bench preventing any significant lateral motion during the blast.

A Discrete element computer program named DMC (Distinct Motion Code) has been developed for modeling rock blasting. This program employs explicit time integration and uses spherical or cylindrical elements which are represented as circles in 2-D. DMC calculations have been compared with measurements on bench blasts in the field with relatively good comparison. Structural rock mass characteristics have a significant impact on any blast and DMC has not, until now, included these effects. This paper discusses a recently added DMC capability for treating joints and bedding planes in bench blast simulations. Material strength is treated in DMC by creating links between spheres to hold them together. The links can be broken based on any criterion; simple tension, compression and shear are currently employed. Joint sets are treated in DMC by defining the dip of each set toward or away from the bench face along with the joint spacing. Strength links that cross joint planes can have their strength properties modified or they can be deleted. Modification of the link patterns based on joint sets creates distinct blocks of spheres outlined by the intersecting joints. These blocks of spheres move together as a solid unit unless stress and strain conditions within the block indicate that links should be broken. Simulations using this capability show some blocks remaining intact throughout the blast and some being partially or completely broken. When this occurs, the joint pattern is shown to influence the characteristics of the blast. Upon completion of this capability both rock breakage and motion will be modeled during the same simulation. Much work remains to be done on this concept making this paper a progress report on the development of this new capability.

Cast blasting can be designed to utilize explosive energy effectively and economically for coal mining operations to remove overburden material. The more overburden removed by explosives, the less blasted material there is left to be transported with mechanical equipment, such as draglines and trucks. In order to optimize the percentage of rock that is cast, a higher powder factor than normal is required plus an initiation technique designed to produce a much greater degree of horizontal muck movement. This paper compares two blast models known as DMC (Distinct Motion Code) and SABREX (Scientific Approach to Breaking Rock with Explosives). DMC, applies discrete spherical elements interacted with the flow of explosive gases and the explicit time integration to track particle motion resulting from a blast. The input to this model includes multi-layer rock properties, and both loading geometry and explosives equation-of-state parameters. It enables the user to have a wide range of control over drill pattern and explosive loading design parameters. SABREX assumes that heave process is controlled by the explosive gases which determines the velocity and time of initial movement of blocks within the burden, and then tracks the motion of the blocks until they come to a rest. In order to reduce computing time, the in-flight collisions of blocks are not considered and the motion of the first row is made to limit the motion of subsequent rows. Although modelling a blast is a complex task, the DMC can perform a blast simulation in 0.5 hours on the SUN SPARCstation 10--41 while the new SABREX 3.5 produces results of a cast blast in ten seconds on a 486-PC computer. Predicted percentage of cast and face velocities from both computer codes compare well with the measured results from a full scale cast blast.

A 3-D finite element analysis was performed to evaluate the stability of the SPR upper and lower oil storage levels at Weeks Island. The mechanical analysis predicted stresses and strains from which pillar stability was inferred using a fracture criterion developed from previous testing of Weeks Island salt. This analysis simulated the sequential mining of the two levels and subsequent oil fill of the mine. The predicted subsidence rates compare well to those measured over the past few years. Predicted failure mechanisms agree with observations made at the time the mine was being modified for oil storage. The modeling technique employed here treats an infinite array of pillars and is a reasonable representation of the behavior at the center of the mine. This analysis predicts that the lower level pillars, at the center of the mine, have fractured and their stability at this time is questionable. Localized pillar fracturing is predicted and implies that the mine is entering a phase of continual time dependent deterioration. Continued and expanded monitoring of the facility and development of methods to assess and predict its behavior are more important now than ever.

The blast-induced movement and final location of geologic layers that may cause environmental problems can be predicted using discrete element methods. This prediction capability can be used by mine operators to locate the material in the muck pile during excavation which would allow encapsulation to prevent groundwater infiltration.

The computer program, DMC (Distinct Motion Code), which was developed for simulating the rock motion associated with blasting, has been used to study the influence of row delay timing on rock motion. The numerical simulations correspond with field observations in that very short delays (< 50ms) and very long delays (> 300ms) produce a lower percent-cast than a medium delay (100 to 200 ms). The DMC predicted relationship between row delay timing and percent-cast is more complex than expected with a dip in the curve where the optimum timing might be expected. More study is required to gain a full understanding of this phenomenon.

A discrete element computer program, DMC (Distinct Motion Code) was developed to simulate blast-induced rock motion. To simplify the complex task of entering material and explosive design parameters as well as bench configuration, a full-featured graphical interface has been developed. DMC is currently executed on both Sun SPARCstation 2 and Sun SPARCstation 10 platforms and routinely used to model bench and crater blasting problems. This paper will document the design and development of the full-featured interface to DMC. The development of the interface will be tracked through the various stages, highlighting the adjustments made to allow the necessary parameters to be entered in terms and units that field blasters understand. The paper also discusses a novel way of entering non-integer numbers and the techniques necessary to display blasting parameters in an understandable visual manner. A video presentation will demonstrate the graphics interface and explains its use.

A discrete element computer program, DMC (Distinct-Motion Code), has been used for several years to simulate blasting-induced rock motion. Recent enhancements of DMC`s capabilities have included addition of an algorithm that couples together rock motion and gas flow. This allows the user to specify a particular explosive which also specifies equation-of-state and other parameters necessary to model explosive gas flow from the blastwell. Rock loading by the flowing gas is calculated automatically. The mechanism for calculating the rock loading is the subject of this paper. The rock motion effects the gas flow calculation by changing the porosity. DMC is currently being used on a SLTN SPARCstation 2 computer workstation.

Geologic materials are often modeled with discrete spheres because the material is not continuous and discrete spherical models simplify the mathematics. Spherical element models have been created using assemblages of spheres with a specified particle size distribution or by assuming the particles are all the same size and making the assemblage a close-packed array of spheres. Both of these approaches yield a considerable amount of material dilatation upon movement. This has proven to be unsatisfactory for sedimentary rock formations that contain bedding planes where shear movement can occur with minimal dilatation of the interface. A new concept referred to as packing angle has been developed to allow the modeler to build arrays of spheres that are the same size but have the rows of spheres offset from each other. ne row offset is a function of the packing angle and allows the modeler to control the dilatation as rows of spheres experience relative horizontal motion.