Search NASA⌕ Search

SEARCH · Search NASA

Results for “matched load”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Failure Simulation Testing of the Z-1 Spacesuit Titanium Bearing Assemblies

The Z-2 is a candidate for NASA's next generation spacesuit, designed for a range of possible missions with enhanced mobility for spacewalks both on planetary surfaces and in microgravity. Increased mobility was accomplished through innovations in shoulder and hip joints, using a number of new bearings to allow spacesuit wearers to dip, walk, and bend with ease; all important tasks for a planetary explorer collecting samples or traveling over rough terrain. The Advanced Spacesuit Development Team of NASA Johnson Space Center requested that the NASA White Sands Test Facility (WSTF) perform a series failure simulation tests on three titanium bearing assemblies, an elemental part of the joint construction used in new spacesuit designs. This testing simulated two undetected failures within the bearings and as a result the objective of this test program was to evaluate whether a failed or failing bearing could result in ignition of the titanium race material due to friction. The first failure was an inner seal leak sufficient to pressurize the race with +99 percent oxygen. The second failure was an improperly installed or mismatched ball port that created a protrusion in the ball bearing race, partially obstructing the nominal rolling path of each ball bearing. When the spacesuit bearings are assembled, bearing balls are loaded into the assembly via a ball port. The ball port is specific and unique to each bearing assembly (matched pair). The simulated mismatched ball port is a significant source of friction, which would be caused by an assembly error. To evaluate this risk, the bearings were cycled in a simulated worst-case scenario environment, with operational loads, and potential flaw conditions. During test the amount of actuation torque required and heat generated through continuous operation were measured and the bearings were observed for sparks or burning events. This paper provides detailed descriptions of the test hardware, methodology, and results.

de Baca, Richard C.↗
NASA NTRS1991

Mantle viscosity - A comparison of models from postglacial rebound and from the geoid, plate driving forces, and advected heat flux

Models of the radial variation of effective viscosity inferred from the earth's response to surface loads associated with Pleistocene deglaciation are compared to structures inferred from models of geodynamic phenomena associated with convection: the geoid, plate-driving forces, and advected heat flux. While observations of the earth's response to surface loads do not have sufficient resolution to justify more than two viscous layers, adequately matching the observed long-wavelength geoid anomalies associated with density contrasts in the lower mantle (inferred from seismic tomography) and in the upper mantle (inferred from a model of subducted slabs) requires more structure. It is possible to explain the geoid, observed plate velocities, the advected heat flux in the lower mantle, and relative sea-level variations in oceanic regions, all with a mantle with a high-viscosity/elastic lid, an asthenospheric channel of 2 x 10 exp 19 Pa s from 100 to 400-km depth, a 6 x 10 exp 20 Pa s transition zone, and a lower mantle of 6 x 10 exp 21 Pa s. The uplift history of Australia, Fennoscandia, and Laurentia can be explained with an asthenospheric viscosity less than a factor of 10 higher. Lateral variations in lower mantle viscosity are not required. Transient creep appears to be unimportant for the recent response-to-surface loads from Pleistocene deglaciation.

Hager, Bradford H.↗
NASA NTRS2003

Contributions of the NASA Langley Research Center to the DARPA/AFRL/NASA/ Northrop Grumman Smart Wing Program

An overview of the contributions of the NASA Langley Research Center (LaRC) to the DARPA/AFRL/NASA/ Northrop Grumman Corporation (NGC) Smart Wing program is presented. The overall objective of the Smart Wing program was to develop smart** technologies and demonstrate near-flight-scale actuation systems to improve the aerodynamic performance of military aircraft. NASA LaRC s roles were to provide technical guidance, wind-tunnel testing time and support, and Computational Fluid Dynamics (CFD) analyses. The program was divided into two phases, with each phase having two wind-tunnel entries in the Langley Transonic Dynamics Tunnel (TDT). This paper focuses on the fourth and final wind-tunnel test: Phase 2, Test 2. During this test, a model based on the NGC Unmanned Combat Air Vehicle (UCAV) concept was tested at Mach numbers up to 0.8 and dynamic pressures up to 150 psf to determine the aerodynamic performance benefits that could be achieved using hingeless, smoothly-contoured control surfaces actuated with smart materials technologies. The UCAV-based model was a 30% geometric scale, full-span, sting-mounted model with the smart control surfaces on the starboard wing and conventional, hinged control surfaces on the port wing. Two LaRC-developed instrumentation systems were used during the test to externally measure the shapes of the smart control surface and quantify the effects of aerodynamic loading on the deflections: Videogrammetric Model Deformation (VMD) and Projection Moire Interferometry (PMI). VMD is an optical technique that uses single-camera photogrammetric tracking of discrete targets to determine deflections at specific points. PMI provides spatially continuous measurements of model deformation by computationally analyzing images of a grid projected onto the model surface. Both the VMD and PMI measurements served well to validate the use of on-board (internal) rotary potentiometers to measure the smart control surface deflection angles. Prior to the final entry, NASA LaRC also performed three-dimensional unstructured Navier Stokes CFD analyses in an attempt to predict the potential aerodynamic impact of the smart control surface on overall model forces and moments. Eight different control surface shapes were selected for study at Mach = 0.6, Reynolds number = 3.25 x 10(exp 6), and + 2 deg., 3 deg., 8 deg., and 10 deg.model angles-of-attack. For the baseline, undeflected control surface geometry, the CFD predictions and wind-tunnel results matched well. The agreement was not as good for the more complex aero-loaded control surface shapes, though, because of the inability to accurately predict those shapes. Despite these results, the NASA CFD study served as an important step in studying advanced control effectors.

Florance, Jennifer P.↗
NASA NTRS2021

Less Atmospheric Radiative Heating By Dust Due to the Synergy of Coarser Size and Aspherical Shape

Mineral dust aerosols cool and warm the atmosphere by scattering and absorbing solar (shortwave: SW) and thermal (longwave: LW) radiation. However, significant uncertainties remain in dust radiative effects, largely due to differences in the dust size distribution and spectral optical properties simulated in Earth system models. Dust models typically underestimate the coarse dust load (more than 2.5 µm in diameter) and assume a spherical shape, which leads to an overestimate of the fine dust load (less than 2.5 µm) after the dust emissions in the models are scaled to match observed dust aerosol optical depth at 550 nm (DAOD550). Here, we improve the simulated dust properties with data sets that leverage measurements of size-resolved dust concentration, asphericity factor, and refractive index in a coupled global chemical transport model with a radiative transfer module. After the adjustment of size-resolved dust concentration and spectral optical properties, the global and annual average of DAOD550 from the simulation increases from 0.023 to 0.029 and falls within the range of a semi-observationally based estimate (0.030 ± 0.005). The reduction of fine dust load after the adjustment leads to a reduction of the SW cooling at the top of the atmosphere (TOA). To improve agreement against a semi-observationally based estimate of the radiative effect efficiency at TOA, we find that a less absorptive SW dust refractive index is required for coarser aspherical dust. Thus, only a minor difference is estimated for the net global dust radiative effect at TOA (−0.08 vs. −0.00 W m−2 on a global scale). Conversely, our sensitivity simulations reveal that the surface warming is substantially enhanced near the strong dust source regions (less cooling to −0.23 from −0.60 W m−2 on a global scale). Thus, less atmospheric radiative heating is estimated near the major source regions (less heating to 0.15 from 0.59 W m−2 on a global scale), because of enhanced LW warming at the surface by the synergy of coarser size and aspherical shape.

simulated dust properties↗
NASA NTRS1972

High sensitivity infrared 10.6 micrometer heterodyne receiver development

The results are presented for a study on the design of an infrared 10.6-micrometer quantum-noise-limited optical receiver subsystem. Performance measurements of the HgCdTe photomixer preamplifier combination were carried out for photomixer temperatures up to 152 K and a photomixer frequency response of up to 420 MHz was obtained. Results of temperature and bias cycling of HgCdTe photomixers are reported. Design considerations for an operational 10.6 micrometer heterodyne receiver are presented. These consist of design data on required laser LO illumination, heat load levels for photomixer cooler, photomixer uniformity and the effects of photomixer impedance match on receiver sensitivity. Analysis and measurements of 10.6 micrometer heterodyne detection in an extrinsic photoconductive (p-type) HgCdTe photomixer are also presented.

Source record↗
NASA NTRS1988

Flexural behavior of the continental lithosphere in Italy - Constraints imposed by gravity and deflection data

This paper presents a simple model that describes the geometry of subducted continental lithosphere as the result of flexural loading of a thin elastic sheet. The method involves a procedure for systematically matching the plate flexure and gravity data to the best-fitting two-dimensional flexural geometry obtained by a thin elastic sheet of specified elastic strength. It is shown that the model can explain both the deflection of the lithosphere and the gravity anomalies observed along four profiles across the Apennine convergent system in Italy. The results of calculations suggest that in subduction systems, like the Apennines, that are characterized by backarc extension, forces acting on subducted lithosphere at depth can be expressed at the surface in the geometry of the foredeep basin that forms adjacent to the subduction zone.

Royden, Leigh↗
NASA NTRS1994

Probabilistic Simulation of Stress Concentration in Composite Laminates

A computational methodology is described to probabilistically simulate the stress concentration factors (SCF's) in composite laminates. This new approach consists of coupling probabilistic composite mechanics with probabilistic finite element structural analysis. The composite mechanics is used to probabilistically describe all the uncertainties inherent in composite material properties, whereas the finite element is used to probabilistically describe the uncertainties associated with methods to experimentally evaluate SCF's, such as loads, geometry, and supports. The effectiveness of the methodology is demonstrated by using is to simulate the SCF's in three different composite laminates. Simulated results match experimental data for probability density and for cumulative distribution functions. The sensitivity factors indicate that the SCF's are influenced by local stiffness variables, by load eccentricities, and by initial stress fields.

Chamis, C. C.↗
NASA NTRS2007

Global Aerosol Optical Models and Lookup Tables for the New MODIS Aerosol Retrieval over Land

Since 2000, MODIS has been deriving aerosol properties over land from MODIS observed spectral reflectance, by matching the observed reflectance with that simulated for selected aerosol optical models, aerosol loadings, wavelengths and geometrical conditions (that are contained in a lookup table or 'LUT'). Validation exercises have showed that MODIS tends to under-predict aerosol optical depth (tau) in cases of large tau (tau greater than 1.0), signaling errors in the assumed aerosol optical properties. Using the climatology of almucantur retrievals from the hundreds of global AERONET sunphotometer sites, we found that three spherical-derived models (describing fine-sized dominated aerosol), and one spheroid-derived model (describing coarse-sized dominated aerosol, presumably dust) generally described the range of observed global aerosol properties. The fine dominated models were separated mainly by their single scattering albedo (omega(sub 0)), ranging from non-absorbing aerosol (omega(sub 0) approx. 0.95) in developed urban/industrial regions, to neutrally absorbing aerosol (omega(sub 0) approx.90) in forest fire burning and developing industrial regions, to absorbing aerosol (omega(sub 0) approx. 0.85) in regions of savanna/grassland burning. We determined the dominant model type in each region and season, to create a 1 deg. x 1 deg. grid of assumed aerosol type. We used vector radiative transfer code to create a new LUT, simulating the four aerosol models, in four MODIS channels. Independent AERONET observations of spectral tau agree with the new models, indicating that the new models are suitable for use by the MODIS aerosol retrieval.

Levy, Robert C.↗
NASA NTRS1974

Single-stage experimental evaluation of tandem-airfoil rotor and stator blading for compressors, part 8

An experimental investigation was conducted with an 0.8 hub/tip ratio, single-stage, axial flow compressor to determine the potential of tandem-airfoil blading for improving the efficiency and stable operating range of compressor stages. The investigation included testing of a baseline stage with single-airfoil blading and two tandem-blade stages. The overall performance of the baseline stage and the tandem-blade stage with a 20-80% loading split was considerably below the design prediction. The other tandem-blade stage, which had a rotor with a 50-50% loading split, came within 4.5% of the design pressure rise (delta P(bar)/P(bar) sub 1) and matched the design stage efficiency. The baseline stage with single-airfoil blading, which was designed to account for the actual rotor inlet velocity profile and the effects of axial velocity ratio and secondary flow, achieved the design predicted performance. The corresponding tandem-blade stage (50-50% loading split in both blade rows) slightly exceeded the design pressure rise but was 1.5 percentage points low in efficiency. The tandem rotors tested during both phases demonstrated higher pressure rise and efficiency than the corresponding single-airfoil rotor, with identical inlet and exit airfoil angles.

Brent, J. A.↗
NASA NTRS1972

Buckling analysis for structural sections and stiffened plates reinforced with laminated composites.

A classical buckling analysis is developed for stiffened, flat plates composed of a series of linked flat plate and beam elements. Plates are idealized as multilayered orthotropic elements; structural beads and lips are idealized as beams. The loaded edges of the stiffened plate are simply supported and the conditions at the unloaded edges can be prescribed arbitrarily. The plate and beam elements are matched along their common junctions for displacement continuity and force equilibrium in an exact manner. Offsets between elements are considered in the analysis. Buckling under uniaxial compressive load for plates, sections and stiffened plates is investigated. Buckling loads are found as the lowest of all possible general and local failure modes and the mode shape is used to determine whether buckling is a local or general instability. Numerical correlations with existing analysis and test data for plates, sections and stiffened plates including boron-reinforced structures are discussed. In general, correlations are reasonably good.

Viswanathan, A. V.↗
NASA NTRS1971

Buckling analysis for axially compressed flat plates, structural sections, and stiffened plates reinforced with laminated composites

A classical buckling analysis is developed for stiffened, flat plates composed of a series of linked plate and beam elements. Plates are idealized as multilayered orthotropic elements. Structural beads and lips are idealized as beams. The loaded edges of the stiffened plate are simply-supported and the conditions at the unloaded edges can be prescribed arbitrarily. The plate and beam elements are matched along their common junctions for displacement continuity and force equilibrium in an exact manner. Offsets between elements are considered in the analysis. Buckling under uniaxial compressive load for plates, sections, and stiffened plates is investigated. Buckling loads are the lowest of all possible general and local failure modes, and the mode shape is used to determine whether buckling is a local or general instability. Numerical correlations with existing analysis and test data for plates, sections, and stiffened plates including boron-reinforced structures are discussed. In general correlations are reasonably good.

Viswanathan, A. V.↗
NASA NTRSDate not supplied

Challenges Associated with In-Situ Calibration of Load Cells in Force Limited Vibration Testing

The difference in mounting configuration between flight and test can significantly impact the effectiveness of the test in environmental vibration testing. Many tests are performed with large electrodynamic shakers, which utilize interfaces that seek to replicate a fixed base, such as slip tables and head expanders. This fixed base configuration is rarely seen in flight configurations; rather a more realistic configuration would include a flexible mounting structure with its own compliance and dynamics. This causes significant over and under tests in various frequency bands depending on the differences between the test article and fixture dynamics. The traditional way of avoiding these high loads is to limit the acceleration responses at multiple locations on the test article. However, the effectiveness of this approach is highly dependent upon the validity of the test article’s analytical in order to derive accurate acceleration response limit specifications. Also this technique requires limiting the acceleration responses at many locations throughout the test article, which may not be implementable due to such things as access and cleanliness issues An improved environmental vibration testing technique known as force limiting incorporates measurements of the forces between the test article and shaker system interface and limiting them to a specification that more accurately replicates the interface impedance of the structure the test article will be mounted to in flight. In effect this transforms the high mechanical impedance at the test article to shaker interface to more closely match the mechanical impedance of the flight interface, which avoids producing the unrealistically high interface loads. Typically force gauges or load cells are used to measure these interface forces. However, utilizing load cells can present a multitude of challenges depending upon such things as their installation method, geometric layout, and test fixture setup. Regardless, it is extremely important to perform an in-situ calibration of the load cells prior to vibration testing at any significant levels. This paper will discuss the challenges associated with utilizing load cells during the NASA Evolutionary Xenon Thruster – Commercial (NEXT-C) gridded ion thruster proto-flight vibration test performed at the NASA Glenn Research Center’s Structural Dynamics Laboratory.

Kenneth J Pederson↗
NASA NTRSDate not supplied

Challenges Associated with In-Situ Calibration of Load Cells in Force Limited Vibration Testing

The difference in mounting configuration between flight and test can significantly impact the effectiveness of the test in environmental vibration testing. Many tests are performed with large electrodynamic shakers, which utilize interfaces that seek to replicate a fixed base, such as slip tables and head expanders. This fixed base configuration is rarely seen in flight configurations; rather a more realistic configuration would include a flexible mounting structure with its own compliance and dynamics. This causes significant over and under tests in various frequency bands depending on the differences between the test article and fixture dynamics. The traditional way of avoiding these high loads is to limit the acceleration responses at multiple locations on the test article. However, the effectiveness of this approach is highly dependent upon the validity of the test article’s analytical in order to derive accurate acceleration response limit specifications. Also this technique requires limiting the acceleration responses at many locations throughout the test article, which may not be implementable due to such things as access and cleanliness issues An improved environmental vibration testing technique known as force limiting incorporates measurements of the forces between the test article and shaker system interface and limiting them to a specification that more accurately replicates the interface impedance of the structure the test article will be mounted to in flight. In effect this transforms the high mechanical impedance at the test article to shaker interface to more closely match the mechanical impedance of the flight interface, which avoids producing the unrealistically high interface loads. Typically force gauges or load cells are used to measure these interface forces. However, utilizing load cells can present a multitude of challenges depending upon such things as their installation method, geometric layout, and test fixture setup. Regardless, it is extremely important to perform an in-situ calibration of the load cells prior to vibration testing at any significant levels. This paper will discuss the challenges associated with utilizing load cells during the NASA Evolutionary Xenon Thruster – Commercial (NEXT-C) gridded ion thruster proto-flight vibration test performed at the NASA Glenn Research Center’s Structural Dynamics Laboratory.

K. J. Pederson↗
NASA NTRS1975

Constant-load versus heart rate-targeted exercise - Responses of systolic intervals

Various systolic intervals were measured prior to and during heart rate-targeted bicycle ergometer exercise. There were striking similarities within each matched exercise set for Q-Im, isovolumetric contraction time, preejection period (PEP), and PEP/left ventricular ejection time (LVET). LVET was significantly shorter for rate-targeted exercise. It is concluded that either constant-load or rate-targeted bicycle ergometry may be used with the choice of method determined by the purpose of the protocol, and that systolic intervals (except LVET) should not be much altered owing to the method chosen.

Lance, V. Q.↗
NASA NTRS2011

Development of Benchmark Examples for Delamination Onset and Fatigue Growth Prediction

An approach for assessing the delamination propagation and growth capabilities in commercial finite element codes was developed and demonstrated for the Virtual Crack Closure Technique (VCCT) implementations in ABAQUS. The Double Cantilever Beam (DCB) specimen was chosen as an example. First, benchmark results to assess delamination propagation capabilities under static loading were created using models simulating specimens with different delamination lengths. For each delamination length modeled, the load and displacement at the load point were monitored. The mixed-mode strain energy release rate components were calculated along the delamination front across the width of the specimen. A failure index was calculated by correlating the results with the mixed-mode failure criterion of the graphite/epoxy material. The calculated critical loads and critical displacements for delamination onset for each delamination length modeled were used as a benchmark. The load/displacement relationship computed during automatic propagation should closely match the benchmark case. Second, starting from an initially straight front, the delamination was allowed to propagate based on the algorithms implemented in the commercial finite element software. The load-displacement relationship obtained from the propagation analysis results and the benchmark results were compared. Good agreements could be achieved by selecting the appropriate input parameters, which were determined in an iterative procedure.

Krueger, Ronald↗
NASA NTRS2008

Different Levels of Eccentric Resistance during Eight Weeks of Training Affect Muscle Strength and Lean Tissue Mass

Coupling concentric and eccentric muscle contractions appears to be important in the development of muscle strength and hypertrophy. The interim Resistive Exercise Device (iRED) currently used aboard the International Space Station does not seem to be as effective as free weight training in ambulatory subjects and has not completely protected against muscular deconditioning due to space flight. The lack of protection during space flight could be caused by iRED's proportionally lower eccentric resistance (60-70%) compared to concentric resistance. PURPOSE: To determine the effects of 8 wks of lower body resistive exercise training using five levels of eccentric resistance on muscle strength and lean tissue mass. METHODS: Forty untrained males (34.9 +/- 7 yrs, 80.9 +/- 9.8 kg, 178.2 +/- 7.1 cm; mean +/- SD) completed three 1-repetition maximum (1-RM) strength tests for both the supine leg press (LP) and supine heel raise (HR) prior to training; subjects were matched for LP strength and randomly assigned to one of five training groups. Concentric load (% 1-RM) was constant across groups during training, but each group trained with different levels of eccentric load (0%, 33%, 66%, 100%, or 138% of concentric). Subjects trained 3 d / wk for 8 wks using a periodized program for LP and HR based on percentages of the highest pre-training 1-RM. LP and HR 1-RM and leg lean mass (LLM; assessed by DEXA) were measured pre- and post-training. A two-way ANOVA was used to analyze all dependent measures. Tukey's post hoc tests were used to test significant main effects. Within group pre- to post-training changes were compared using paired t-tests with a Bonferroni adjustment. Statistical significance was set a priori at p 0.05. All data are expressed as mean +/- SE. RESULTS: LP 1-RM strength increased significantly in all groups pre- to post-training. The 138% group increase (20.1 +/- 3.7%) was significantly greater than the 0% (7.9 +/- 2.8%), 33% (7.7 +/- 4.6%), and 66% (7.5 +/- 4.3%) groups. All groups significantly increased HR strength pre- to posttraining (33%: 7.5 +/- 6.1%; 66%: 6.6 +/- 3.7%; 100%: 12.2 +/- 1.8%; 138%: 11.0 +/- 6.4%) except for the 0% (4.9 +/- 9.1%) group. There were no differences between groups. LLM increased significantly pre- to post-training in only the 138% group; there were no differences between groups. CONCLUSIONS: Eight wks of lower body resistive exercise training with eccentric overload resulted in greater increases in LP strength than training with eccentric loads of 66% or less. Post-training HR strength was not affected by eccentric training load, perhaps because of the predominance of Type I fibers typical in the gastrocnemius. Only 138% eccentric training significantly increased LLM. PRACTICAL APPLICATIONS: For athletes or others desiring to maximize muscle strength and hypertrophy gains, training with eccentric loads greater than 100% of concentric resistance will provide greater increases in muscle strength and lean tissue mass in some muscle groups. In a rehabilitation or geriatric exercise setting that places primary emphasis on program adherence and moderate strength gains, training with an eccentric underload may provide strength increases comparable to those of traditional 1:1 training but with less muscle soreness and physiologic insult to the patient, but this has yet to be proven.

English, K. L.↗
NASA NTRS2017

Thermal Modeling of Mars Ground for Surface Missions

Thermal analysis packages are capable of including ground temperature effects for orbital and surface thermal analyses. In particular, Thermal Desktop® offers the option of specifying ground temperatures as a function of time for planetary surface modeling. While suitable for many cases, this approach is not sufficient if an object has local interactions with the ground that could significantly affect the ground temperature. Ground modeling is necessary for the Mars 2020 rover thermal design and analysis since shadowing and heat dissipation from the rover’s Multi-Mission Radioisotope Thermoelectric Generator (MMRTG) can result in significant temperature deviations of the local ground. The Mars 2020 thermal team is explicitly modeling the Martian ground so that these local temperature effects can be captured. The upper portion of the ground is modeled, and material and optical properties of the ground are varied in order to match data collected from orbiting science instruments. Atmospheric surface temperature, sky temperatures, and solar loads from a Mars General Circulation Model (GCM) are used as boundary conditions, resulting in a ground surface temperature profile consistent with the GCM predictions. The rover model is then placed on this modeled ground so that the effects of shadowing and MMRTG dissipation on the ground temperature can be captured.

Farias, Edgardo↗
NASA NTRS1993

Probabilistic simulation of stress concentration in composite laminates

A computational methodology is described to probabilistically simulate the stress concentration factors in composite laminates. This new approach consists of coupling probabilistic composite mechanics with probabilistic finite element structural analysis. The probabilistic composite mechanics is used to probabilistically describe all the uncertainties inherent in composite material properties while probabilistic finite element is used to probabilistically describe the uncertainties associated with methods to experimentally evaluate stress concentration factors such as loads, geometry, and supports. The effectiveness of the methodology is demonstrated by using it to simulate the stress concentration factors in composite laminates made from three different composite systems. Simulated results match experimental data for probability density and for cumulative distribution functions. The sensitivity factors indicate that the stress concentration factors are influenced by local stiffness variables, by load eccentricities and by initial stress fields.

Chamis, C. C.↗