Publications

40 Results
Skip to search filters

Stochastic Modeling in a Multimaterial Continuum Mixture Shock Physics Code

Schumacher, Shane C.

Stochastic modelling approaches are presented to capture random effects at multiple time and length scales. Random processes that occur at the microscale produce nondeterministic effects at the macroscale. Here we present three stochastic modeling approaches that describe random processes at microscopic length scales and map these processes to the macroscopic length scale. The first stochastic modeling approach is based upon a particle based numerical technique to solve a Stochastic Differential Equation (SDE) using an arbitrary diffusion process to capture random processes at the microstructural level. The second approach prescribes a Probability Density Function (PDF) for the drift and diffusion of the random variable derived using the forward and backward Kolmogorov equations. This method requires mean and drift evolution PDF transport equations. The third approach is the coupling of multiple random variables which are dependent on each other. The relationship of the PDFs and a coupling function, known as a copula, produces a Joint Probability Density Function (JPDF). These stochastic modeling approaches are implemented into a Multiple Component (MC) shock physics computational code and used to model statistical fracture and reactive flow applications.

More Details

Dynamic shock response of an S2 glass/SC15 epoxy woven fabric composite material system

AIP Conference Proceedings

Key, Christopher T.; Alexander, Charles S.; Harstad, Eric N.; Schumacher, Shane C.

The use of S2 glass/SC15 epoxy woven fabric composite materials for blast and ballistic protection has been an area of on-going research over the past decade. In order to accurately model this material system within potential applications under extreme loading conditions, a well characterized and understood anisotropic equation of state (EOS) is needed. This work details both an experimental program and associated analytical modelling efforts which aim to provide better physical understanding of the anisotropic EOS behavior of this material. Experimental testing focused on planar shock impact tests loading the composite to peak pressures of 15 GPa in both the transverse and longitudinal orientations. Test results highlighted the anisotropic response of the material and provided a basis by which the associated numeric micromechanical investigation was compared. Results of the combined experimental and numerical modeling investigation provided insights into not only the constituent material influence on the composite response but also the importance of the plain weave microstructure geometry and the significance of the microstructural configuration.

More Details

Evaluation of a strain based failure criterion for the multi-constituent composite model under shock loading

EPJ Web of Conferences

Key, Christopher T.; Schumacher, Shane C.; Alexander, Charles S.

This study details and demonstrates a strain-based criterion for the prediction of polymer matrix composite material damage and failure under shock loading conditions. Shock loading conditions are characterized by high-speed impacts or explosive events that result in very high pressures in the materials involved. These material pressures can reach hundreds of kbar and often exceed the material strengths by several orders of magnitude. Researchers have shown that under these high pressures, composites exhibit significant increases in stiffness and strength. In this work we summarize modifications to a previous stress based interactive failure criterion based on the model initially proposed by Hashin, to include strain dependence. The failure criterion is combined with the multi-constituent composite constitutive model (MCM) within a shock physics hydrocode. The constitutive model allows for decomposition of the composite stress and strain fields into the individual phase averaged constituent level stress and strain fields, which are then applied to the failure criterion. Numerical simulations of a metallic sphere impacting carbon/epoxy composite plates at velocities up to 1000 m/s are performed using both the stress and strain based criterion. These simulation results are compared to experimental tests to illustrate the advantages of a strain-based criterion in the shock environment.

More Details

Improved understanding of the dynamic response in anisotropic directional composite materials through the combination of experiments and modeling

Journal of Physics: Conference Series

Alexander, Charles S.; Key, C.T.; Schumacher, Shane C.

Recently there has been renewed interest in the dynamic response of composite materials; specifically low density epoxy matrix binders strengthened with continuous reinforcing fibers. This is in part due to the widespread use of carbon fiber composites in military, commercial, industrial, and aerospace applications. The design community requires better understanding of these materials in order to make full use of their unique properties. Planar impact testing was performed resulting in pressures up to 15 GPa on a unidirectional carbon fiber - epoxy composite, engineered to have high uniformity and low porosity. Results illustrate the anisotropic nature of the response under shock loading. Along the fiber direction, a two-wave structure similar to typical elastic-plastic response is observed, however, when shocked transverse to the fibers, only a single bulk shock wave is detected. At higher pressures, the epoxy matrix dissociates resulting in a loss of anisotropy. Greater understanding of the mechanisms responsible for the observed response has been achieved through numerical modeling of the system at the micromechanical level using the CTH hydrocode. From the simulation results it is evident that the observed two-wave structure in the longitudinal fiber direction is the result of a fast moving elastic precursor wave traveling in the carbon fibers ahead of the bulk response in the epoxy resin. Similarly, in the transverse direction, results show a collapse of the resin component consistent with the experimental observation of a single shock wave traveling at speeds associated with bulk carbon. Experimental and simulation results will be discussed and used to show where additional mechanisms, not fully described by the currently used models, are present. © Published under licence by IOP Publishing Ltd.

More Details

Composite layering technique for use in a Eulerian shock physics code

AIP Conference Proceedings

Schumacher, Shane C.

The high strength and low density characteristics of fiber reinforced composite materials have made them applicable to a large variety of applications. As these applications grow, their performance in high strain rate shock environments has increased. The modeling and simulation of such materials is difficult due to their anisotropic behavior and complex internal geometries. Fiber reinforced composite materials consist of a collection of layers that create a laminate. Each layer is typically transverse isotropic or orthotropic consisting of a fiber and matrix material. One approach is to explicitly model each layer, while accurate, this is often not feasible for full system calculations as the laminate layer count increases in size. Additionally, modeling each layer given the finite thickness proves to be a challenging process and typically a smearing approach is used to represent the laminate response removing the identity and material response of each layer. The creation of a layering capability is a good compromise between the inaccuracy of smearing and the computational cost of explicitly modeling each layer. The layering is done using a sub-grid technique in an individual grid cell. The grid cell is partitioned based on layer location in the laminate and the material deformation. The volume occupied by the given layer is computed and the layer calculates a material response based on the cell strain field. The resulting material stress and state variables are volume weighted with the remaining layers in the given grid cell yielding a cell response. The result is a technique that requires less computation time than modeling each layer while increasing the accuracy over smeared approximations. © 2012 American Institute of Physics.

More Details
40 Results
40 Results