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.
The behavior of thin square perforated plates under the action of uniform shear deformation is studied experimentally and analytically using finite element analysis. Elastic Shear buckling strength is established as a function of the diameter of a round, centrally located hole in the plate. Post buckling behavior and the behavior of perforated plates with various ring stiffeners are also studied experimentally.
This paper considers the stress concentration problems of two types of cracked composite systems: (1) a composite system with a broken fiber (a penny-shaped crack problem), and (2) a composite system with a cracked matrix (an annular crack problem). The cracked composite systems are modeled with triangular and trapezoidal ring finite elements. Using NASTRAN (NASA Structural Analysis) finite element computer program, the stress and deformation fields in the cracked composite systems are calculated. The effect of fiber-matrix material combination on the stress concentrations and on the crack opening displacements is studied.
System of three computer programs developed for use with BOPACE finite-element program for analyzing cumulative plastic deformation of structures subject to high thermomechanical load cycles.
Based on theories of laminate anisotropic elasticity and interlaminar fracture, the complete solution structure associated with a composite delamination is determined. Fracture mechanics parameters characterizing the interlaminar crack behavior are defined from asymptotic stress solutions for delaminations with different crack-tip deformation configurations. A numerical method employing singular finite elements is developed to study delaminations in fiber composites with any arbitrary combinations of lamination, material, geometric, and crack variables. The special finite elements include the exact delamination stress singularity in its formulation. The method is shown to be computationally accurate and efficient, and operationally simple. To illustrate the basic nature of composite delamination, solutions are shown for edge-delaminated (0/-0/-0/0) and (+ or - 0/+ or - 0/90/90 deg) graphite-epoxy systems under uniform axial extenstion. Three-dimensional crack-tip stress intensity factors, associated energy release rates, and delamination crack-closure are determined for each individual case. The basic mechanics and mechanisms of composite delamination are studied, and fundamental characteristics unique to recently proposed tests for interlaminar fracture toughness of fiber composite laminates are examined.
Based on theories of laminate anisotropic elasticity and interlaminar fracture, the complete solution structure associated with a composite delamination is determined. Fracture mechanics parameters characterizing the interlaminar crack behavior are defined from asymptotic stress solutions for delaminations with different crack-tip deformation configurations. A numerical method employing singular finite elements is developed to study delaminations in fiber composites with any arbitrary combinations of lamination, material, geometric, and crack variables. The special finite elements include the exact delamination stress singularity in its formulation. The method is shown to be computationally accurate and efficient, and operationally simple. To illustrate the basic nature of composite delamination, solutions are shown for edge-delaminated (0/-0/-0/0) and (+ or - 0/+ or - 0/90/90 deg) graphite-epoxy systems under uniform axial extension. Three-dimensional crack-tip stress intensity factors, associated energy release rates, and delamination crack-closure are determined for each individual case. The basic mechanics and mechanisms of composite delamination are studied, and fundamental characteristics unique to recently proposed tests for interlaminar fracture toughness of fiber composite laminates are examined. Previously announced in STAR as N84-13222
A numerical simulation was performed, by means of a large-scale finite element code capable of handling large deformations and/or nonlinear behavior, to investigate the suitability of the nonlinear pulse-control algorithm to suppress the vibrations induced in the Spacecraft Control Laboratory Experiment (SCOLE) components under realistic maneuvers. Among the topics investigated were the effects of various control parameters on the efficiency and robustness of the vibration control algorithm. Advanced nonlinear control techniques were applied to an idealized model of some of the SCOLE components to develop an efficient algorithm to determine the optimal locations of point actuators, considering the hardware on the SCOLE project as distributed in nature. The control was obtained from a quadratic optimization criterion, given in terms of the state variables of the distributed system. An experimental investigation was performed on a model flexible structure resembling the essential features of the SCOLE components, and electrodynamic and electrohydraulic actuators were used to investigate the applicability of the control algorithm with such devices in addition to mass-ejection pulse generators using compressed air.
The aeromechanical stability of a 1/8th Froude scaled bearingless rotor model was investigated experimentally in a wind tunnel. Both shaft-fixed and shaft-free conditions were examined to study the aeroelastic stability of a bearingless rotor without the incorporation of auxiliary dampers. This wind tunnel investigation generated a set of stability data for four different advance ratios, and a wide range of collective pitch settings. Theoretical analysis was performed using the newly developed University of Maryland Advanced Rotorcraft Code (UMARC). For analysis, the blade is modeled as an elastic beam undergoing flap bending, lag bending, elastic twist, and axial deformation. Blade response is calculated using a finite element method in time. Nonlinear aerodynamic effects are included by using a semiempirical stall modeling. The linearized periodic rotor perturbation equations in the nonrotating frame are solved for stability roots using Floquet transition matrix theory, as well as constant coefficient approximation. The predicted results are compared with the experimental data.
A simple method is developed for solving plane-plastic-stress problems with axial symmetry in the strain-hardening range which is based on the deformation theory of plasticity employing the finite-strain concept. The equations defining the problems are first reduced to two simultaneous nonlinear differential equations involving two dependent variables: (a) the octahedral shear strain, and (b) a parameter indicating the ratio of principal stresses. By multiplying the load and dividing the radius by an arbitrary constant, it is possible to solve these problems without iteration for any value of the modified load. The constant is determined by the boundary condition. This method is applied to a circular membrane under pressure, a rotating disk without and with a central hole, and an infinite plate with a circular hole. Two materials, inconel x and 16-25-6, the octahedral shear stress-strain relations of which do not follow the power law, are used. Distributions of octahedral shear strain, as well as of principal stresses and strains, are obtained. These results are compared with the results of the same problems in the elastic range.
A three-dimensional finite element program called VISCOPAC was developed and used to conduct a micromechanics analysis of titanium metal matrix composites. The VISCOPAC program uses a modified Eisenberg-Yen thermo-viscoplastic constitutive model to predict matrix behavior under thermomechanical fatigue loading. The analysis incorporated temperature-dependent elastic properties in the fiber and temperature-dependent viscoplastic properties in the matrix. The material model was described and the necessary material constants were determined experimentally. Fiber-matrix interfacial behavior was analyzed using a discrete fiber-matrix model. The thermal residual stresses due to the fabrication cycle were predicted with a failed interface, The failed interface resulted in lower thermal residual stresses in the matrix and fiber. Stresses due to a uniform transverse load were calculated at two temperatures, room temperature and an elevated temperature of 650 C. At both temperatures, a large stress concentration was calculated when the interface had failed. The results indicate the importance of accuracy accounting for fiber-matrix interface failure and the need for a micromechanics-based analytical technique to understand and predict the behavior of titanium metal matrix composites.
Time-accurate unsteady flow simulations are being performed supporting the SRM T+68sec pressure 'spike' anomaly investigation. The anomaly occurred in the RH SRM during the STS-54 flight (STS-54B) but not in the LH SRM (STS-54A) causing a momentary thrust mismatch approaching the allowable limit at that time into the flight. Full-motor internal flow simulations using the USA-2D axisymmetric code are in progress for the nominal propellant burn-back geometry and flow conditions at T+68-sec--Pc = 630 psi, gamma = 1.1381, T(sub c) = 6200 R, perfect gas without aluminum particulate. In a cooperative effort with other investigation team members, CFD-derived pressure loading on the NBR and castable inhibitors was used iteratively to obtain nominal deformed geometry of each inhibitor, and the deformed (bent back) inhibitor geometry was entered into this model. Deformed geometry was computed using structural finite-element models. A solution for the unsteady flow has been obtained for the nominal flow conditions (existing prior to the occurrence of the anomaly) showing sustained standing pressure oscillations at nominally 14.5 Hz in the motor IL acoustic mode that flight and static test data confirm to be normally present at this time. Average mass flow discharged from the nozzle was confirmed to be the nominal expected (9550 lbm/sec). The local inlet boundary condition is being perturbed at the location of the presumed reconstructed anomaly as identified by interior ballistics performance specialist team members. A time variation in local mass flow is used to simulate sudden increase in burning area due to localized propellant grain cracks. The solution will proceed to develop a pressure rise (proportional to total mass flow rate change squared). The volume-filling time constant (equivalent to 0.5 Hz) comes into play in shaping the rise rate of the developing pressure 'spike' as it propagates at the speed of sound in both directions to the motor head end and nozzle. The objectives of the present analysis are to: (1) capture the dynamic responses of the motor combustion gas flow to correlate with available low-frequency (less than 12.5 sample/sec) data and (2) observe the high-frequency (up to 50 Hz) characteristics of the response to determine any potentials for dynamic coupling.
A computer program implements a finite-element-based method of predicting the deformations of thin aerospace structures made of isotropic materials or anisotropic fiber-reinforced composite materials. The technique and corresponding software are applicable to thin shell structures in general and are particularly useful for analysis of thin beamlike members having open cross-sections (e.g. I-beams and C-channels) in which significant warping can occur.
Ten commercial laboratories participated in an interlaboratory study to establish the repeatability and reproducibility of compression strength tests conducted according to ASTM International Standard Test Method E9. The test employed a cylindrical aluminum AA2024-T351 test specimen. Participants measured elastic modulus and 0.2 % offset yield strength, YS(0.2 % offset), using an extensometer attached to the specimen. The repeatability and reproducibility of the yield strength measurement, expressed as coefficient of variations were cv(sub r)= 0.011 and cv(sub R)= 0.020 The reproducibility of the test across the laboratories was among the best that has been reported for uniaxial tests. The reported data indicated that using diametrically opposed extensometers, instead of a single extensometer doubled the precision of the test method. Laboratories that did not lubricate the ends of the specimen measured yield stresses and elastic moduli that were smaller than those measured in laboratories that lubricated the specimen ends. A finite element analysis of the test specimen deformation for frictionless and perfect friction could not explain the discrepancy, however. The modulus measured from stress-strain data were reanalyzed using a technique that finds the optimal fit range, and applies several quality checks to the data. The error in modulus measurements from stress-strain curves generally increased as the fit range decreased to less than 40 % of the stress range.
The performance of a three-dimensional continuum damage mechanics (CDM) model for intralaminar damage, coupled with a cohesive zone model for delamination, is examined for damage prediction in cross-ply and quasi-isotropic laminated plates subjected to quasi-static transverse indentation. Details of the CDM model based on the deformation gradient decomposition (DGD) method and the finite element (FE) modeling strategy are presented and discussed. Critical features of the damage model and FE modeling approach are assessed by comparing analysis predictions to experimental data for cross-ply and quasi-isotropic panels in which damage is dominated by complex interactions between matrix cracks and delaminations.
Efficient, semi-analytical models of bonded joints are utilized for initial joint design and sizing due to their convenience and speed. Many of these models are based on beam and plate theories, which are most applicable when in-plane stresses/strains dominate. While this can be a sensible assumption for obtaining the deformations of the joint, out of plane (transverse) stresses are often of critical importance for failure in bonded joints composed of fiber reinforced composite materials. Since through-thickness strains are neglected in these semi-analytical theories, the present study shows a method to back-calculate the out of plane stresses based on stress balance equations, which use the in-plane stresses and displacements that have already been obtained. The method was implemented within the Joint Element Designer (JED) software, and verification and validation examples have been shown. As is typical for these types of models, corner stress singularities are not captured by the JED model, but the general trends of the stress and deformation fields agreed with a higher fidelity finite element model. As such, this method for obtaining out of plane stresses, when coupled with the JED code, is very useful for the complete stress analysis of bonded composite joints.
A weak formulation in structural analysis that provides well conditioned matrices suitable for iterative solutions is presented. A mixed formulation ensures the proper representation of the problem and the constitutive relations are added in a penalized form. The problem is solved by a double conjugate gradient algorithm combined with an element by element approximate factorization procedure. The double conjugate gradient strategy resembles Uzawa's variable-length type algorithms the main difference is the presence of quadratic terms in the mixed variables. In the case of shear deformable beams these terms ensure that the proper finite thickness solution is obtained.
The in-flight elastic wing twist of a fighter-type aircraft was studied to provide for an improved on-board real-time computed prediction of pointing variations of three wing store stations. This is an important capability to correct sensor pod alignment variation or to establish initial conditions of iron bombs or smart weapons prior to release. The original algorithm was based upon coarse measurements. The electro-optical Flight Deflection Measurement System measured the deformed wing shape in flight under maneuver loads to provide a higher resolution database from which an improved twist prediction algorithm could be developed. The FDMS produced excellent repeatable data. In addition, a NASTRAN finite-element analysis was performed to provide additional elastic deformation data. The FDMS data combined with the NASTRAN analysis indicated that an improved prediction algorithm could be derived by using a different set of aircraft parameters, namely normal acceleration, stores configuration, Mach number, and gross weight.
A static aeroelastic analysis capability that can predict aerodynamic loads for the deformed shape of the composite wing has been developed. The finite element method (FEM) was used for composite plate structural analysis, and the linear vortex lattice method (VLM) was used for steady aerodynamic analysis. The final deformed shape of the wing due to applied forces is determined by iterative manner using FEM and VLM. FEM and VLM analysis are related by a surface spline interpolation procedure. The wing with Gr/Ep composite material has been investigated to see the wing deformation effect. Aerodynamic load change due to wing flexibility has been investigated. Also, the effect of fiber orientation and sweep angle on the deformation pattern and aerodynamic coefficients are examined. For a certain fiber orientation, the deflection and aerodynamic loading of the composite wing is very much reduced. The swept forward wing has more significant effect of wing flexibility on aerodynamic coefficient than the swept back wing does.