SEARCH · Engineering Papers

Results for “Finite amplitude method”

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.

SourceResearch connectionFrom yearThrough year

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

184 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS2018

Global High Resolution Crustal Magnetic Field at the Surface of the Moon from Low-Altitude Lunar Prospector Magnetic Gradient Data

We derived new vector gradients based models of crustal magnetic field at the lunar surface with data from the Lunar Prospector (LP) satellite using two model parameterization approaches: a global set of 35820 1° spaced (~30 km) equal area monopoles at 20 km below the surface (O’Brien and Parker, 1994; Olsen et al., 2017) and combined results of subsets of 100000 0.66° spaced monopoles at the same depth. We use the scheme of iteratively reweighted least-squares inversion to compute the initial model. Then the amplitudes of these monopoles are determined by minimizing the misfit to the components together with the global average of |Br| at the ellipsoid surface (i.e. applying a L1 model regularization of Br). In previous approaches using vector fields for modeling, we found that external field contamination leads to spurious anomalies in the downward continued field models even with stringent data selection criteria and ad-hoc noise removal techniques (e.g., satellite’s position in the Moon’s wake w.r.t. the solar wind and in the Earth’s magnetotail, internal/external dipoles fields removal, low-order polynomial removal, joint equivalent source cross-validation technique and visually removing remaining anomalous segments). On the other hand, with the use of gradients-only data (along-track first finite differences), we were able to completely bypass the ad-hoc techniques. Similar processing of Kaguya magnetic data, which have only higher altitude coverage over the most of the Moon except in the region of the South Pole–Aitken (SPA) basin, completely misses some of the anomalies seen in the Lunar Prospector data. The combined Lunar Prospector and Kaguya gradient-based models also severely degrade the derivation of the anomaly fields in many regions. Euler analysis of isolated anomaly features from the Reiner Gamma swirl suggests top depths of about 0.3 to 1.5 km and center depths of 10-14 km; in the region between Stein and Vallier craters north of the SPA basin our analysis suggests top depths of around 1.5 km and center depths of 13-15 km. With the spectral depth determination techniques, the SPA basin region yields depths to the base of magnetization ranging between 15 and 40 km. Three-dimensional modeling and the bulk magnetization determinations of the sources constrained by the Euler and spectral methods is underway.

Gradient↗

NASA NTRS2004

Reduced-Order Blade Mistuning Analysis Techniques Developed for the Robust Design of Engine Rotors

The primary objective of this research program is to develop vibration analysis tools, design tools, and design strategies to significantly improve the safety and robustness of turbine engine rotors. Bladed disks in turbine engines always feature small, random blade-to-blade differences, or mistuning. Mistuning can lead to a dramatic increase in blade forced-response amplitudes and stresses. Ultimately, this results in high-cycle fatigue, which is a major safety and cost concern. In this research program, the necessary steps will be taken to transform a state-of-the-art vibration analysis tool, the Turbo-Reduce forced-response prediction code, into an effective design tool by enhancing and extending the underlying modeling and analysis methods. Furthermore, novel techniques will be developed to assess the safety of a given design. In particular, a procedure will be established for using eigenfrequency curve veerings to identify "danger zones" in the operating conditions--ranges of rotational speeds and engine orders in which there is a great risk that the rotor blades will suffer high stresses. This work also will aid statistical studies of the forced response by reducing the necessary number of simulations. Finally, new strategies for improving the design of rotors will be pursued. Several methods will be investigated, including the use of intentional mistuning patterns to mitigate the harmful effects of random mistuning, and the modification of disk stiffness to avoid reaching critical values of interblade coupling in the desired operating range. Recent research progress is summarized in the following paragraphs. First, significant progress was made in the development of the component mode mistuning (CMM) and static mode compensation (SMC) methods for reduced-order modeling of mistuned bladed disks (see the following figure). The CMM method has been formalized and extended to allow a general treatment of mistuning. In addition, CMM allows individual mode mistuning, which accounts for the realistic effects of local variations in blade properties that lead to different mistuning values for different mode types (e.g., mistuning of the first torsion mode versus the second flexural mode). The accuracy and efficiency of the CMM method and the corresponding Turbo-Reduce code were validated for an example finite element model of a bladed disk.

Min, James B.↗

NASA NTRS2013

Engineering Evaluation of International Low Impact Docking System Latch Hooks

The international Low Impact Docking System (iLIDS) provides a structural arrangement that allows for visiting vehicles to dock with the International Space Station (ISS) (Fig 1). The iLIDS docking units are mechanically joined together by a series of active and passive latch hooks. In order to preserve docking capability at the existing Russian docking interfaces, the iLIDS latch hooks are required to conform to the existing Russian design. The latch hooks are classified as being fail-safe. Since the latch hooks are fail-safe, the hooks are not fracture critical and a fatigue based service life assessment will satisfy the structural integrity requirements. Constant amplitude fatigue testing to failure on four sets of active/passive iLIDS latch hooks was performed at load magnitudes of 10, 11, and 12 kips. Failure analysis of the hook fatigue failures identified multi-site fatigue initiation that was effectively centered about the hook mid-plane (consistent with the 3D model results). The fatigue crack initiation distribution implies that the fatigue damage accumulation effectively results in a very low aspect ratio surface crack (which can be simulated as thru-thickness crack). Fatigue damage progression resulted in numerous close proximity fatigue crack initiation sites. It was not possible to determine if fatigue crack coalescence occurs during cyclic loading or as result of the fast fracture response. The presence of multiple fatigue crack initiation sites on different planes will result in the formation of ratchet marks as the cracks coalesce. Once the stable fatigue crack becomes unstable and the fast fracture advances across the remaining ligament and the plane stress condition at a free-surface will result in failure along a 45 deg. shear plane (slant fracture) and the resulting inclined edge is called a shear lip. The hook thickness on the plane of fatigue crack initiation is 0.787". The distance between the shear lips on this plane was on the order of 0.48" and it was effectively centered about the mid-plane of the section. The numerous ratchet marks between the shear lips on the fracture initiation plane are indicative of multiple fatigue initiation sites within this region. The distribution of the fatigue damage about the centerline of the hook is consistent with the analytical results that demonstrate peak stress/strain response at the mid-plane that decreases in the direction of the hook outer surfaces. Scanning electron microscope images of the failed sections detected fatigue crack striations in close proximity to the free surface of the hook radius. These findings were documented at three locations on the fracture surface : 1) adjacent to the left shear lip, 2) adjacent to the right shear lip, and 3) near the centerline of the section. The features of the titanium fracture surface did not allow for a determination of a critical crack size via identification of the region where the fatigue crack propagation became unstable. The fracture based service life projections where benchmarked with strain-life analyses. The strainrange response in the hook radius was defined via the correlated finite element models and the modified method of universal slopes was incorporated to define the strain-life equation for the titanium alloy. The strain-life assessment confirmed that the fracture based projections were reasonable for the loading range of interest. Based upon the analysis and component level fatigue test data a preliminary service life capability for the iLIDS active and passive hooks of 2 lifetimes is projected (includes a scatter factor of 4).

Martinez, J.↗

NASA NTRS2018

Implementation of J-A Methodology Elastic-Plastic Crack Instability Analysis Capability into the WARP-3D Code

Characterization of the near crack-tip stress/strain fields is the foundation of fracture mechanics. The description of the near tip stress field and the prediction of when fracture occurs is well established for brittle materials that exhibit linear elastic behavior. However, in ductile materials or conditions that violate linear elastic assumptions (Aluminum alloys, Al 2024-T3, Al 2024- T351 etc.), the elastic-plastic crack-tip stress fields are characterized by the Hutchison-Rice-Rosengren (HRR) field. The J-integral is commonly used to characterize amplitude of the HRR field under elastic-plastic conditions. The J-integral has been demonstrated for crack-tip fields that are under high constraint conditions (i.e., small-scale plasticity where the J-dominance is maintained). However, as the external load increases, yielding changes from small- to largescale plasticity and usually a loss of constraint (i.e., reduction in the triaxial stress field along the crack front). The loss of constraint leads to the deviation of the crack-tip stress fields from that given by the HRR field. Hence, the J-dominance will be gradually lost and additional parameter(s) are required to quantify the crack-tip stress fields and predict fracture behavior. The assessment objectives were to: 1) implement a two-parameter (i.e., J-A) fracture criterion into an elastic-plastic three-dimensional (3D) finite element analysis (FEA), 2) validate the implementation by comparison with the A parameter from literature data, 3) conduct material characterization tests to quantify the material behavior and provide fracture data for validation of the J-A fracture criteria, and (4) perform evaluations to establish if the J-A criteria can be used to predict fracture in a ductile metallic material (e.g., aluminum alloys). The A parameter in these criteria is the second parameter in a three-term elastic-plastic asymptotic expansion of the neartip stress behavior. A series of extensive FEAs were performed using WARP3D software package to obtain solutions for the A parameter for different specimen configurations. The methodology needed for the estimation of the A parameter in the asymptotic expansion was developed and implemented using Matlab®. A user material (UMAT) routine was used to model the material stress-strain response using a Ramberg-Osgood power law with a hardening exponent (n) and a material coefficient (alpha). This UMAT routine was successfully implemented in WARP3D software and validated through comparison with the experimental data. Three configurations were extracted from published results: 1) center cracked plate (CCP), 2) single edge-cracked plate (SECP), and 3) double edge-cracked plate (DECP). These configurations and four other configurations (three-hole tension (THT)), three-point bend (3PTB), three-hole compact tension (3PCT), and compact tension (CT)) were analyzed to verify the methodology that was developed and implemented into WARP3D. Solutions of the A parameter were obtained for remote tension loading conditions that started with small-scale yielding and continued into the large-scale plasticity regime. The results indicate that the methodology developed can be used to calculate the elastic-plastic J-A parameters for test specimens with a range of crack geometries, material strain hardening behaviors, and loading conditions. The J-A parameters were implemented as fracture criteria and used to predict the test results. For comparison, other fracture criteria were used to predict the same test results. Major findings include: The A constraint parameter A varies with specimen type and applied load thus accurate determination is crucial in predicting the failure load, and the A parameter is asymptotic as the failure load is approached, making an accurate determination difficult (i.e., small differences in the A parameter can cause large variations in failure load) for materials exhibiting elastic-plastic behavior. The failure predictions from J-A methodology were more accurate than the traditionally used KC and J methods, and have comparable scatter to that observed when using the crack-tip opening angle (CTOA) method. However, the J-A methodology requires considerable effort (expertise level and labor) to implement and to evaluate the A parameter for different specimen types and materials, or to apply this methodology to part-through crack (e.g., 3D problems) structural applications.

Hamm, Kenneth R., Jr.↗