Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum mechanics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Correlations Between Metallurgical Characterization Studies, Exploratory Mechanical Tests, and Continuum Mechanics Approaches to Constitutive Equations

Austenitic stainless steels, such as types 316 and 304, are widely used as pressure vessel materials in the temperature range of 425 to 650 C. Stainless steel specimens were tested to rupture at two different stress levels sigma and sigma 2 sigma 1 sigma 2) to establish the normal stain-time behavior. A subsequent test was performed in which the specimen was crept at the higher stress (sigma 1) to the beginning of the secondary stage of creep, presumed to be the strain/time conditions at which a steady state microstructure is developed, and then the stress was reduced to the lower level (sigma 2). The associated microstructure, and significance of this microstructure on the creep strain-hardening model for variable uniaxial loads were assesed and found to be consistent with the use of creep-recovery models at high stresses and temperatures and strain-hardening models at low stresses and tempertures.

Moteff, J.↗

On the Fidelity of Computational Models for the Flow of Milled Loblolly Pine: A Benchmark Study on Continuum-Mechanics Models and Discrete-Particle Models

The upstream of bioenergy industry has suffered from unreliable operations of granular biomass feedstocks in handling equipment. Computational modeling, including continuum-mechanics models and discrete-particle models, offers insightful understandings and predictive capabilities on the flow of milled biomass and can assist equipment design and optimization. This paper presents a benchmark study on the fidelity of the continuum and discrete modeling approaches for predicting granular biomass flow. We first introduce the constitutive law of the continuum-mechanics model and the contact law of the coarse-grained discrete-particle model, with model parameters calibrated against laboratory characterization tests of the milled loblolly pine. Three classical granular material flow systems (i.e., a lab-scale rotating drum, a pilot-scale hopper, and a full-scale inclined plane) are then simulated using the two models with the same initial and boundary conditions as the physical experiments. The close agreement of the numerical predictions with the experimental measurements on the hopper mass flow rate, the hopper critical outlet width, the material stopping thickness on the inclined plane, and the dynamic angle of repose, clearly indicates that the two methods can capture the critical flow behavior of granular biomass. The qualitative comparison shows that the continuum-mechanics model outperforms in parameterization of materials and wall friction, and large-scale systems, while the discrete-particle model is more preferred for discontinuous flow systems at smaller scales. Industry stakeholders can use these findings as guidance for choosing appropriate numerical tools to model biomass material flow in part of the optimization of material handling equipment in biorefineries.

Jin, Wencheng↗

The einstein equivalence principle, intrinsic spin and the invariance of constitutive equations in continuum mechanics

The invariance of constitutive equations in continuum mechanics is examined from a basic theoretical standpoint. It is demonstrated the constitutive equations which are not form invariant under arbitrary translational accelerations of the reference frame are in violation of the Einstein equivalane principle. Furthermore, by making use of an analysis based on statistical mechanics, it is argued that any frame-dependent terms in constitutive equations must arise from the intrinsic spin tensor and are negligible provided that the ratio of microscopic to macroscopic time scales is extremely small. The consistency of these results with existing constitutive theories is discussed in detail along with possible avenues of future research.

Speziale, Charles G.↗

A highly stable explicit integration technique for computational continuum mechanics

A user-oriented subroutine package is built around a highly stable explicit integration algorithm for solution of large order systems of ordinary differential equations, as result from discretization of initial-boundary value problems in continuum mechanics. Fast and accurate solutions, for problems in laminar and turbulent, two and three-dimensional viscous flow fields and multi-dimensional transient heat transfer, are presented using this algorithm, as embodied within a general purpose finite element computer program.

Baker, A. J.↗

Post-instability in continuous systems. I - Failure of differentiability of solutions in continuum mechanics

It is pointed out that mathematical models of continua are based on certain assumptions regarding functions which must be at least piece-wise differentiable. The assumption about smoothness of the functions makes it possible to use the mathematical technique of differentiable equations. However, this artificial mathematical limitation follows neither from the principles of mechanics nor from the definition of a continuum. The price paid for such a mathematical convenience is instability (in the class of smooth functions) of the solutions to the corresponding governing equations in some regions of the parameters. A new mathematical technique should, therefore, be developed to describe the solutions which are not necessarily differentiable. The present investigation is concerned with the criteria of applicability of the classical models of continua from the view point of stability of the corresponding solutions, postinstability models derived by reformulation of the original models, and postinstability models in enlarged classes of functions.

Zak, M.↗

An MPMD approach coupling electromagnetic continuum mechanics approximations in ALEGRA

In this work, two complementary approximations for describing aspects of continuum electromagnetics in moving media are discussed: electroquasistatic and magnetoquasistatic. Each has been implemented in the finite element shock code ALEGRA for modeling dynamic electromechanical phenomena on typical engineering time scales, with fully integrated circuit coupling. The approximations can be obtained by consistent asymptotic balancing of Maxwell’s equations relative to timescales associated with magnetic diffusion, charge relaxation, and electromagnetic wave propagation. In ALEGRA, the electroquasistatic approximation is used for ferroelectric (FE) modeling, while the magnetoquasistatic approximation is used for magnetohydrodynamic (MHD) modeling. In this paper we introduce for the first time a detailed derivation of a useful quasi-steady “low-R m ” variant of the MHD approximation applicable for cases, such as with detonators, where the thermodynamic pressure arising from Joule heating dominates over magnetic forces. An additional purpose of this paper is to present a coupling mode using Multiple Program-Multiple Data (MPMD) message passing communication that allows the user to run 3D FE problems together with 2D and/or 3D MHD problems with the respective simulation domains coupled through a common circuit equation. The MPMD coupling capability is used here to model the dynamic coupling of a notional ferroelectric generator with an RP-87 exploding bridgewire detonator. The simulated bridgewire heats up and bursts under current generated by simulated depoling of the ferroelectric generator, as a demonstration of the MPMD capability.

42 ENGINEERING↗

A space-time tensor formulation for continuum mechanics in general curvilinear, moving, and deforming coordinate systems

Tensor methods are used to express the continuum equations of motion in general curvilinear, moving, and deforming coordinate systems. The space-time tensor formulation is applicable to situations in which, for example, the boundaries move and deform. Placing a coordinate surface on such a boundary simplifies the boundary condition treatment. The space-time tensor formulation is also applicable to coordinate systems with coordinate surfaces defined as surfaces of constant pressure, density, temperature, or any other scalar continuum field function. The vanishing of the function gradient components along the coordinate surfaces may simplify the set of governing equations. In numerical integration of the equations of motion, the freedom of motion of the coordinate surfaces provides a potential for enhanced resolution of the continuum field function. An example problem of an incompressible, inviscid fluid with a top free surface is considered, where the surfaces of constant pressure (including the top free surface) are coordinate surfaces.

Avis, L. M.↗

A mathematical model of post-instability in continuum mechanics

The post-instability can appear in the form of failure of hyperbolicity or in the form of cascade instability (stretching of vorticity in an inviscid fluid or instability of the Navier-Stokes equations). All the cases are characterized by an unlimited decrease of the scale of the motions, in the course of which the derivatives of the corresponding functions tend to infinity, though the functions themselves remain finite. Such an instability shows that the corresponding state of a continuum cannot be properly simulated by smooth functions. Hence, the original model must be corrected by giving up the requirement about differentiability and enlarging the class of functions. Applications to a model of an inviscid fluid and to a model of turbulence are discussed.

Zak, M.↗

A large deformation multiphase continuum mechanics model for shock loading of soft porous materials

A large deformation, coupled finite-element (FE) model is developed to simulate the multiphase response of soft porous materials subjected to high strain-rate loading. The approach is based on the theory of porous media (TPM) at large deformations. Simplifications to the one-dimensional regime studied in the numerical simulations follow. An overview of several different time integration schemes is presented for the purpose of solving the nonlinear dynamic coupled balance of momenta (mixture and fluid) and balance of mass of the mixture equations. Numerical examples are presented for (i) verification against closed-form analytical solutions assuming small loads, (ii) demonstrating large deformation effects at high strain-rate, and (iii) showing differences in deformations between a single-phase elastodynamics model with occluded compressible pore fluid and a multiphase poroelastodynamics model at high strain-rate. The multiphase model shows that the relative motion of the pore fluid significantly dampens the deformation response of the solid skeleton as compared to the single-phase model, and makes it possible to extract quantitative values for the stresses of the different constituents, thereby allowing one to form preliminary conclusions about the onset of damage in the solid skeleton. The novelty of the current work is developing a multiphase, large deformation, mixture theory numerical model for high strain-rate loading of soft porous materials. It was discovered that explicit, adaptive time-stepping Runge–Kutta schemes offer high accuracy at relatively low cost when compared to traditional implicit or explicit central difference time-stepping schemes for shock-like loadings. Here, shock viscosity is added to the mixture momentum balance equation to regularize the shock front, and a stabilization term is added to the mixture mass balance equation to stabilize equal order interpolation finite elements for the coupled finite element solution of multiphase materials.

Engineering↗

Green's function fast multipole method for continuum mechanics (SM-FMM)

Solid Mechanics Fast Multipole Method based on elastic Green's function and accelerated using FFTs The code calculates the mechanical fields (stress and strain) for a heterogeneous elasto-plastic material unit cell under quasi-static conditions (no dynamic effects). The fields are calculated using a Green's function method, where the strain field is given by a discrete convolution of Green's operator with an auxiliary stress field. The convolution is calculated using the fast multipole method.

Zecevic, Miroslav↗

Time-dependent density functional theory of charge, energy, and spin transport, and quantum continuum mechanics (Final Report)

The overarching goal of this project, which unfolded over a period of seventeen years, including two no-cost extensions in the last two years, was to extend and deepen our ability to use time-dependent density functional theory and related theoretical methods for the description of the transport of charge, spin, or energy in quantum electronic systems at the nanoscale. The output of our efforts is contained in 52 papers, listed below, which we have published since 2006. There have also been many presentations at international conferences, workshops, colloquia and seminars. This body of work has been carried out with the essential input of 8 postdoctoral students (who were supported by the grant) and several external senior collaborators and their postdocs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

3D Progressive Damage Modeling for Laminated Composite Based on Crack Band Theory and Continuum Damage Mechanics

A simple continuum damage mechanics (CDM) based 3D progressive damage analysis (PDA) tool for laminated composites was developed and implemented as a user defined material subroutine to link with a commercially available explicit finite element code. This PDA tool uses linear lamina properties from standard tests, predicts damage initiation with an easy-to-implement Hashin-Rotem failure criteria, and in the damage evolution phase, evaluates the degradation of material properties based on the crack band theory and traction-separation cohesive laws. It follows Matzenmiller et al.'s formulation to incorporate the degrading material properties into the damaged stiffness matrix. Since nonlinear shear and matrix stress-strain relations are not implemented, correction factors are used for slowing the reduction of the damaged shear stiffness terms to reflect the effect of these nonlinearities on the laminate strength predictions. This CDM based PDA tool is implemented as a user defined material (VUMAT) to link with the Abaqus/Explicit code. Strength predictions obtained, using this VUMAT, are correlated with test data for a set of notched specimens under tension and compression loads.

Wang, John T.↗

Implementation of a Matrix Crack Spacing Parameter in a Continuum Damage Mechanics Finite Element Model

Continuum Damage Mechanics (CDM) based progressive damage and failure analysis (PDFA) methods have demonstrated success in a variety of finite element analysis (FEA) implementations. However, the technical maturity of CDM codes has not yet been proven for the full design space of composite materials in aerospace applications. CDM-based approaches represent the presence of damage by changing the local material stiffness definitions and without updating the original mesh or element integration schemes. Without discretely representing cracks and their paths through the mesh, damage in models with CDM-based materials is often distributed in a region of partially damaged elements ahead of stress concentrations. Having a series of discrete matrix cracks represented by a softened region may affect predictions of damage propagation and, thus, structural failure. This issue can be mitigated by restricting matrix damage development to discrete, fiber-aligned rows of elements; hence CDM-based matrix cracks can be implemented to be more representative of discrete matrix cracks. This paper evaluates the effect of restricting CDM matrix crack development to discrete, fiber-aligned rows where the spacing of these rows is controlled by a user-defined crack spacing parameter. Initially, the effect of incrementally increasing matrix crack spacing in a unidirectional center notch coupon is evaluated. Then, the lessons learned from the center notch specimen are applied to open-hole compression finite element models. Results are compared to test data, and the limitations, successes, and potential of the matrix crack spacing approach are discussed.

Hyder, Imran↗

Fabric-enriched continuum breakage mechanics (F-CBM )

A fabric-enriched continuum breakage mechanics (F-CBM) framework is developed to examine the relation between elastic and inelastic anisotropy in granular materials subjected to high pressure. For this purpose, emphasis is given to the influence of the anisotropic fabric on the plastic cap regime – that is, a macroscopic signature of the inelasticity of granular materials widely dependent on the particle crushability. Anisotropic behaviour prior to yielding is introduced through a symmetric second-order fabric tensor embedded in the expression of the elastic energy potential. It is shown that, thanks to its energy-based formulation, the F-CBM framework provides a platform to predict the rotation and distortion of the yield cap of granular materials as an outcome of anisotropic strain energy storage prior to grain crushing. Parametric analyses indicate that both the degree of anisotropy and the orientation of the fabric have an impact on the energy release resulting from grain crushing and contribute to the resulting stress–strain behaviour. It is shown that, despite the use of very few parameters, the proposed framework can accurately predict the anisotropic yielding and inelastic stress–strain response solely on the basis of the directional properties inferred from measurements of elastic anisotropy.

Engineering↗

Application of the Continuum Damage Mechanics Wilshire-Cano-Stewart (WCS) Model

In this study, the applications of the continuum damage mechanics-based Wilshire-Cano-Stewart (WCS) model are explored to predict rupture time, minimum-creep-strain-rate (MCSR), damage, damage evolution, and creep deformation. Increase knowledge in manufacturing methods has pushed the limit of material science and the development of new materials. Conventional testing is required to qualify materials against creep which according to the ASME B&PV III code, 10,000+ hours of experiments are necessary for each heat before materials are put into service. This process is costly and not feasible for new materials. As an alternative, models have been employed to predict creep behaviors and reduce the amount of time necessary for material qualification. Many models have been developed to predict distinct creep behaviors and the question of which model is best remains. Amongst current models, the WCS model has emerge with the ability to predict multiple behaviors using an explicit analytical approach with the ability to predict long-term creep. In this study the novel continuum damage mechanics WCS is employed in multiple applications. The goals of the study are (a) to discuss and determine the framework of the WCS model and validated it mathematically and using parametric simulations, (b) applied the model to accelerated creep data to show the capabilities of the model with non-conventional data, and (c) applied a novel numerical method, the datum temperature method (DTM) to show the model extrapolations and interpolations capabilities with limited and reduce data sets. To accomplish these goals, data is gathered for alloy P91 and Inconel 718 to develop and post-audit validate the model. The benefits of using the WCS model is that it provides an explicit stress and temperature dependency ideal for extrapolations, the coupled equations are suitable for finite element analysis (FEA) implementation, and it follows an explicit calibration approach. The model also proves that it can be applied to accelerated testing data and using the DTM.

Cano, Jaime A↗

Differential continuum damage mechanics models for creep and fatigue of unidirectional metal matrix composites

Three multiaxial isothermal continuum damage mechanics models for creep, fatigue, and creep/fatigue interaction of a unidirectional metal matrix composite volume element are presented, only one of which will be discussed in depth. Each model is phenomenological and stress based, with varying degrees of complexity to accurately predict the initiation and propagation of intergranular and transgranular defects over a wide range of loading conditions. The development of these models is founded on the definition of an initially transversely isotropic fatigue limit surface, static fracture surface, normalized stress amplitude function and isochronous creep damage failure surface, from which both fatigue and creep damage evolutionary laws can be obtained. The anisotropy of each model is defined through physically meaningful invariants reflecting the local stress and material orientation. All three transversely isotropic models have been shown, when taken to their isotropic limit, to directly simplify to previously developed and validated creep and fatigue continuum damage theories. Results of a nondimensional parametric study illustrate (1) the flexibility of the present formulation when attempting to characterize a large class of composite materials, and (2) its ability to predict anticipated qualitative trends in the fatigue behavior of unidirectional metal matrix composites. Additionally, the potential for the inclusion of various micromechanical effects (e.g., fiber/matrix bond strength, fiber volume fraction, etc.), into the phenomenological anisotropic parameters is noted, as well as a detailed discussion regarding the necessary exploratory and characterization experiments needed to utilize the featured damage theories.

Arnold, S. M.↗