Parameterization of Battery Electrothermal Models Coupled With Finite Element Flow Models for Cooling
Not provided.
SEARCH · Search NASA
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.
Not provided.
Solar sails employ a unique form of propulsion, gaining momentum from incident and reflected photons. However, the momentum transferred by an individual photon is extremely small. Consequently, a solar sail must have an extremely large surface area and also be extremely light. The flexibility of the sail then must be considered when designing or evaluating control laws. In this paper, solar sail flexibility and its influence on control effectiveness is considered using idealized two-dimensional models to represent physical phenomena rather than a specific design. Differential equations of motion are derived for a distributed parameter model of a flexible solar sail idealized as a rotating central hub with two opposing flexible booms. This idealization is appropriate for solar sail designs in which the vibrational modes of the sail and supporting booms move together allowing the sail mass to be distributed along the booms in the idealized model. A reduced analytical model of the flexible response is considered. Linear feedback torque control is applied at the central hub. Two translational disturbances and a torque disturbance also act at the central hub representing the equivalent effect of deflecting sail shape about a reference line. Transient simulations explore different control designs and their effectiveness for controlling orientation, for reducing flexible motion and for disturbance rejection. A second model also is developed as a two-dimensional "pathfinder" model to calculate the effect of solar sail shape on the resultant thrust, in-plane force and torque at the hub. The analysis is then extended to larger models using the finite element method. The finite element modeling approach is verified by comparing results from a two-dimensional finite element model with those from the analytical model. The utility of the finite element modeling approach for this application is then illustrated through examples based on a full finite element model.
Understanding stress distributions in solid materials is complicated by the fact that most materials are polycrystalline in nature, with each crystal having an elastically anisotropic reaction to force. This study is to gain a better understanding on how external forces placed on a polycrystal are related to internal reactions within and between the grains. The hypothesis is that the stress distribution in porous to fully dense materials are self-organized based on strong contacts between and within the individual grains created by force chains. Force chains, commonly known in loaded granular materials, and could be the phenomenon that connect micro to macro deformation. Scale bridging measurements conducted through Raman spectroscopy, Atomic Force Microscopy, and Digital Image Correlation will be used to create stress maps, modulus maps, and elastic strain maps across a variety of geological and pharmaceutical polycrystals. When possible, the resulting maps will be compared to current homogenization schemes and a full field models. Finite element modeling will be used to assess whether the patterning seen in the experimental map is a reasonable approximation based on the orientation data of the samples used. A minimum of three publications is projected to be accomplished focusing each on a different method to experimentally test and analyze stress distributions.
This report describes the methodology used for the safety basis analysis of the Building (Bldg.) 3525 casks and presents preliminary results for a selected geometry. The specific problem analyzed is a free fall load case onto an unyielding target at various orientations. The various orientations were analyzed to ensure that the packages were subjected to the orientation that causes the most damage. The analysis approach, material modeling, finite element modeling, and results for the selected case are presented, along with a brief discussion of conservatism and uncertainties. The results for the analyzed load cases show that through selected administrative controls, no breach of the container or loss of shielding will occur in the package.
Design for prevention of aeroelastic instability (that is, the critical speeds leading to aeroelastic instability lie outside the operating range) is an integral part of the wing design process. Availability of the sensitivity derivatives of the various critical speeds with respect to shape parameters of the wing could be very useful to a designer in the initial design phase, when several design changes are made and the shape of the final configuration is not yet frozen. These derivatives are also indispensable for a gradient-based optimization with aeroelastic constraints. In this study, flutter characteristic of a typical section in subsonic compressible flow is examined using a state-space unsteady aerodynamic representation. The sensitivity of the flutter speed of the typical section with respect to its mass and stiffness parameters, namely, mass ratio, static unbalance, radius of gyration, bending frequency, and torsional frequency is calculated analytically. A strip theory formulation is newly developed to represent the unsteady aerodynamic forces on a wing. This is coupled with an equivalent plate structural model and solved as an eigenvalue problem to determine the critical speed of the wing. Flutter analysis of the wing is also carried out using a lifting-surface subsonic kernel function aerodynamic theory (FAST) and an equivalent plate structural model. Finite element modeling of the wing is done using NASTRAN so that wing structures made of spars and ribs and top and bottom wing skins could be analyzed. The free vibration modes of the wing obtained from NASTRAN are input into FAST to compute the flutter speed. An equivalent plate model which incorporates first-order shear deformation theory is then examined so it can be used to model thick wings, where shear deformations are important. The sensitivity of natural frequencies to changes in shape parameters is obtained using ADIFOR. A simple optimization effort is made towards obtaining a minimum weight design of the wing, subject to flutter constraints, lift requirement constraints for level flight and side constraints on the planform parameters of the wing using the IMSL subroutine NCONG, which uses successive quadratic programming.
Triply Periodic Minimal Surface (TPMS) lattice structures have been of increasing interest due to their light weighting, enhanced mechanical properties, and energy absorption characteristics for automotive and biomedical applications. With the advent of additive manufacturing and geometric modeling software, TPMS lattices with complex geometries can be realized. In this work, TPMS lattice structures were fabricated with PLA using fused filament fabrication (FFF) and their dynamic properties are characterized through drop tower experiments. Although lightweight TPMS lattices are beneficial for their impact absorption capability, most of the existing works are limited to quasi-static compression and dynamic impact tests are rarely performed. The current study investigates the stress-strain and energy absorption characteristics of TPMS lattices through drop tower testing and numerical modeling. Finite element modeling for TPMS lattices is carried out to validate the experimental responses. The mechanical properties, deformation, and failure mechanisms of TPMS lattices under dynamic impact are summarized for potential future applications.
A fully coupled thermal–microstructural–mechanical finite element modelling framework is developed to investigate the distortion and residual stresses during directed energy deposition (DED) of multi-phase Ti–6Al–4V alloy. The Johnson–Cook constitutive model is used to predict the yield strength of each phase as a function of strain, strain rate and temperature where the flow stress is calculated by a linear mixing rule based on the volumetric phase fractions. A thin-walled rectangular sample is chosen as the reference geometry and the results are compared with experimentally measured in situ thermal history and distortion data, where a reasonable agreement is achieved. The proposed modelling framework with physics-based material constitutive model provides useful information for a better understanding of process–microstructure–property relations in additive manufacturing by DED.
The AH-1G NASTRAN finite element model (FEM) is described and the correlations with measured data that were conducted to verify the model are summarized. Comparisons of the AH-1G NASTRAN FEM calculations with measured data include the following: (1) fuselage and tailboom static load deflection (stiffness) testing, (2) airframe ground vibration testing (0-30 H<), (3) airframe flight vibration testing (main rotor, 2,4, and 6/rev), and (4) tailboom effective skin static testing. A description of the modeling rationale and techniques used to develop the NASTRAN FEM is presented in conjunction with all previous correlation work. In general, the correlations show good agreement between analysis and test in stiffness and vibration response through 15 to 20 Hz. For higher frequencies (equal to or greater than 4/rev (21.6 Hz)), the vibration responses generally did not agree well. Also, the lateral (2/rev (10.8 Hz)) flight vibration responses were much lower in the FEM than test, indicating that there is a significant excitation source other than at the main rotor hub that is affecting the lateral vibrations, such as downwash impingement on the vertical tail.
Bayesian finite element (FE) model updating using direct model evaluations of large-scale high-fidelity FE models is extremely computationally expensive. Surrogate models can be used as fast emulators of FE models to accelerate the model calibration process. The physics/mechanics-based FE models are still the underpinning behind the surrogate models. Here, this paper evaluates the loss in accuracy and the gain in computational time while performing Bayesian model updating by using surrogate model evaluations compared to using direct FE model evaluations. This evaluation is crucial before entirely relying on surrogate models in model updating for structural health monitoring (SHM) and damage prognosis (DP) purposes. This paper also demonstrates Bayesian updating and surrogate model construction of large-scale high-fidelity FE models of infrastructure systems. In this regard, the miter gate structural system is considered as the testbed structure. Three predominant damage modes (loss of contact between gate and wall, loss of thickness due to corrosion, and loss of tension in the diagonal rods) are considered for model updating purposes. Bayesian model updating is performed using direct FE evaluations by leveraging parallel computing. Two types of surrogates, namely polynomial chaos expansion (PCE) and Gaussian process regression (GPR), are developed for the miter gate. Model updating is performed again using the trained surrogate models, and the updating results are compared with their counterparts obtained using the direct FE evaluation results. The posterior distribution of the FE model parameters obtained using the trained surrogates are sufficiently accurate with respect to the posterior obtained utilizing the direct FE evaluations. In addition, an approximate 4-fold decrease in the computational time was observed when using surrogate model evaluations instead of direct FE evaluations for model updating.
Creep strength-enhanced ferritic (CSEF) steels such as Grade 91 (9Cr-1Mo-V) and Grade 92 (Fe-9Cr-2W-0.5Mo) steels are widely used in the fossil-fuel-fired and nuclear power plants. The weld integrity of these steels is crucial for power plants' safe and reliable operations. Due to harsh service conditions, the steel weld can become susceptible to environmental degradation. Field welding repair is used to restore the degraded weld’s performance where a controlled temper-bead welding technique is commonly used to temper the freshly formed martensite during welding. However, knowledge of weld repairability is limited and experimental trial and error optimization to achieve desired microstructure and joint properties is expensive and time-consuming. Many existing computational models, e.g., finite element models, are limited to solving heat conduction equation and ignoring convective heat transfer due to molten metal flow. These models can result in over-prediction of peak temperatures of weld pool and heat-affected zone (HAZ), which in turn can affect the accuracy of tempering prediction. Moreover, these finite element models require an input of the deposit profiles in advance and thus limits the usability of these models. Here, a molten pool-based, multi-pass multi-layer model has been developed based on computational fluid dynamics (CFD) approach with the Volume of Fluid (VOF) method. The model calculates the bead formation, thereby eliminating the need for pre-determined bead profiles required by finite element models. For computational efficiency, a coordinate system attached to the moving heat source is utilized. A subroutine is developed to convert the temperature profiles in the reference frame stationary to the heat source to that stationary to the workpiece. The converted thermal cycles are then imported into a microstructure model to compute the tempering kinetics and resultant hardness using a Johnson-Mehl-Avrami-Kolmogorov (JMAK), and modified Grange-Baughman parameter. The modeling approach is first developed and validated on single- and multi-pass deposition of stainless steel filler metal onto a SA-533 high strength steel substrate. The models are then applied to a multi-pass V-groove repair weld of Grade 91 steel plate as well as directed energy deposition of Grade 92 steel. Non-destructive characterization of microstructures was performed on Grade 91 and 92 steel welds. Two welding processes, cold metal transfer (CMT) and flux-cored arc welding (FCAW), were investigated for the Grade 91 steel weld samples. For the Grade 92 weld samples, three different heat inputs (low, medium, and high) of gas tungsten arc welding (GTAW) were utilized to replicate traditional field welding processes. The non-destructive evaluation (NDE) method used for this research was immersion ultrasonic testing (UT) using a micro-resolution ultrasonic imaging methodology specifically designed to operate in the through-transmission configuration operating at 20 MHz of frequency. The system used a focused ultrasonic beam spot size diameter between 250-300 μm, and a 6 μm laser vibrometer spot size for detection, to produce highly defined images with longitudinal and mode-converted shear waves. From the micro-resolution ultrasonic C-scan images, three microstructural regions, i.e., weld metal (WM), HAZ, and base metal (BM), were clearly identifiable. Various levels of ultrasonic amplitudes distributed over the three regions were correlated with electron beam backscattered diffraction (EBSD) images using grain size, grain boundaries, and dislocation densities. The results showed that areas with relatively higher ultrasonic amplitude levels were associated with smaller grains and higher dislocation densities, while areas with lower amplitude levels were associated with larger grains and lower dislocation densities. In addition, ultrasonic velocity data obtained across the three different weld microstructural regions of Grade 91 test samples were correlated with optical metallographic images and hardness measurements. The results showed distinctive decreases in ultrasonic velocity and hardness over the HAZ region, where weld failures often occur during service.
Finite element models are often developed for load validation, structural certification, response predictions, and to study alternate design concepts. In rare occasions, models developed with a nominal set of parameters agree with experimental data without the need to update parameter values. Today, model updating is generally heuristic and often performed by a skilled analyst with in-depth understanding of the model assumptions. Parameter uncertainties play a key role in understanding the model update problem and therefore probabilistic analysis tools, developed for reliability and risk analysis, may be used to incorporate uncertainty in the analysis. In this work, probability analysis (PA) tools are used to aid the parameter update task using experimental data and some basic knowledge of potential error sources. Discussed here is the first application of PA tools to update parameters of a finite element model for a composite wing structure. Static deflection data at six locations are used to update five parameters. It is shown that while prediction of individual response values may not be matched identically, the system response is significantly improved with moderate changes in parameter values.
AutoCAD-to-GIFTS translator program, ACTOG, developed to facilitate quick generation of small finite-element models using CASA/GIFTS finite-element modeling program. Reads geometric data of drawing from Data Exchange File (DXF) used in AutoCAD and other PC-based drafting programs. Geometric entities recognized by ACTOG include points, lines, arcs, solids, three-dimensional lines, and three-dimensional faces. From this information, ACTOG creates GIFTS SRC file, which then reads into GIFTS preprocessor BULKM or modified and reads into EDITM to create finite-element model. SRC file used as is or edited for any number of uses. Written in Microsoft Quick-Basic (Version 2.0).
Computational finite element (FE) models are used in suited astronaut injury risk assessments; however, these models’ verification, validation, and credibility (VV&C) procedures for simulating injuries in altered gravity environments are limited. Our study conducts VV&C assessments of THUMS and Elemance whole-body FE models for predicting suited astronaut injury biomechanics using eight credibility factors, as per NASA-STD-7009A. Credibility factor ordinal scores are assigned by reviewing existing documentation describing VV&C practices, and credibility sufficiency thresholds are assigned based on input from subject matter experts. Our results show the FE models are credible for suited astronaut injury investigation in specific ranges of kinematic and kinetic conditions correlating to highway and contact sports events. Nevertheless, these models are deficient when applied outside these ranges. Several credibility elevation strategies are prescribed to improve models’ credibility for the NASA-centric application domain.
Computational finite element (FE) models are used in suited astronaut injury risk assessments; however, these models’ verification, validation, and credibility (VV&C) procedures for simulating injuries in altered gravity environments are limited. Our study conducts VV&C assessments of THUMS and Elemance whole-body FE models for predicting suited astronaut injury biomechanics using eight credibility factors, as per NASA-STD-7009A. Credibility factor ordinal scores are assigned by reviewing existing documentation describing VV&C practices, and credibility sufficiency thresholds are assigned based on input from subject matter experts. Our results show the FE models are credible for suited astronaut injury investigation in specific ranges of kinematic and kinetic conditions correlating to highway and contact sports events. Nevertheless, these models are deficient when applied outside these ranges. Several credibility elevation strategies are prescribed to improve models’ credibility for the NASA-centric application domain.
A mesoscale finite element material model is proposed to analyze structures that fail by the fiber kinking damage mode. To evaluate the assumptions of the mesoscale model, the results were compared with those of a high-fidelity micromechanical model. A direct comparison between the two models shows remarkable correlation, indicating that the key features of the fiber kinking phenomenon are appropriately accounted for in the mesoscale model. The mesoscale model is applied to structural analysis cases to demonstrate the capabilities of the model. A verification study is conducted with an unnotched compression specimen and preliminary validation is demonstrated with a notched compression specimen. The results show that the model is successful at representing the kinematics of fiber kinking while at the same time highlighting the need for further verification and validation.
The automotive industry has used computational methods for crashworthiness since the early 1970's. These methods have ranged from simple lumped parameter models to full finite element models. The emergence of the full finite element models in the mid 1980's has significantly altered the research direction. However, there remains a need for both simple, rapid modeling methods and complex detailed methods. Some directions for continuing research are discussed.
A sensitivity-based methodology for improving the finite element model of a given structure using test modal data and a few sensors is presented. The proposed method searches for both the location and sources of the mass and stiffness errors and does not interfere with the theory behind the finite element model while correcting these errors. The updating algorithm is derived from the unconstrained minimization of the squared L sub 2 norms of the modal dynamic residuals via an iterative two-step staggered procedure. At each iteration, the measured mode shapes are first expanded assuming that the model is error free, then the model parameters are corrected assuming that the expanded mode shapes are exact. The numerical algorithm is implemented in an element-by-element fashion and is capable of 'zooming' on the detected error locations. Several simulation examples which demonstate the potential of the proposed methodology are discussed.
In this study, laser powder bed fusion (L-PBF) additive manufacturing (AM) is a key enabling technology to manufacture highly complex and integrated metallic structures. In L-PBF AM process, the melting of the metal powders and the layers underneath can be governed by either “conduction mode” or “keyhole mode”, with the keyhole mode reportedly leading to porosity and decreased strength and ductility by many studies. In part scale simulations, finite element (FE) model is often used to study the temperature distribution during printing and to predict the residual stress, where a volumetric heat flux with a Gaussian or a double ellipsoidal (Goldak) distribution is often applied as the laser heat source. However, the above heat source models can only capture the melt pool shape in the conduction mode, and fail to capture the transition to keyhole melting mode when the process parameters change. To overcome this inaccuracy, an extended Goldak heat source model is proposed by introducing a laser penetration term as a function of laser parameters obtained from a Gaussian-Process (GP) model. The model is validated by “2D pad” AlSi10Mg L-PBF experiments under a wide range of laser power, scan speed, and laser focus offset, and the results show the model successfully captures the measured melt pool shape in all conditions.