Large deflection and stability analysis by the direct stiffness method
Matrix numerical procedure by digital computer equipment, and application to elastic stability and deflection problems
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.
Matrix numerical procedure by digital computer equipment, and application to elastic stability and deflection problems
Reproducing in vitro the complex multiscale physical features of human tissues creates novel biomedical opportunities and fundamental understanding of cell–environment interfaces and interactions. While stiffness has been recognized as a key driver of cell behavior, systematic studies on the role of stiffness have been limited to values in the KPa–MPa range, significantly below the stiffness of bone. Here, in this work, a platform enabling the tuning of the stiffness of a biocompatible polymeric interface up to values characteristic of human bone is reported, which are in the GPa range, by using extremely thin polymer films on glass and cross-linking the films using ultraviolet (UV) light irradiation. It is shown that a higher stiffness is related to better adhesion, proliferation, and osteogenic differentiation, and that it is possible to switch on/off cell attachment and growth by solely tuning the stiffness of the interface, without any surface chemistry or topography modification. Since the stiffness is tuned directly by UV irradiation, this platform is ideal for rapid and simple fabrication of stiffness patterns and gradients, thus representing an innovative tool for combinatorial studies of the synergistic effect of tissue environmental cues on cell behavior, and creates new opportunities for next-generation biosensors, single-cell patterning, and lab-on-a-chip devices.
The goal of the project was to develop the necessary analysis tools for a feasibility study of a cable suspended robot system for examining the space shuttle orbiter payload bay radiators These tools were developed to address design issues such as workspace size, tension requirements on the cable, the necessary accuracy and resolution requirements and the stiffness and movement requirements of the system. This report describes the mathematical models for studying the inverse kinematics, statics, and stiffness of the robot. Each model is described by a matrix. The manipulator Jacobian was also related to the stiffness matrix, which characterized the stiffness of the system. Analysis tools were then developed based on the singular value decomposition (SVD) of the corresponding matrices. It was demonstrated how the SVD can be used to quantify the robot's performance and to provide insight into different design issues.
This paper presents a generalized and efficient method for quasi-static analysis of mooring systems, including complex scenarios such as when shared mooring lines interconnect multiple floating wind or wave energy devices. While quasi-static mooring models are well established, most published formulations are focused on specific applications, and no publicly available implementations provide efficient handling of large mooring system networks. The present formulation addresses these gaps by: (1) formulating solutions for edge cases not typically supported by quasi-static models; (2) creating a fully generalized model structure such that any combination of mooring lines, point masses, and floating bodies can be assembled; and (3) deriving analytic expressions for the system Jacobians (stiffness matrices) so that systems with many degrees of freedom can be solved efficiently. These techniques form the theory basis of MoorPy, an open-source mooring analysis library. The model is demonstrated on nine scenarios of increasing complexity with features of interest for offshore renewable energy applications. When compared with steady-state results from a lumped-mass dynamic model, the results show that the quasi-static formulation accurately calculates profiles and tensions and that its analytic approach provides more efficient and reliable computation of system stiffness matrices than finite-differencing methods. These results verify the accuracy of the MoorPy model.
Machine learning (ML) models are increasingly being used in many engineering fields due to the advancements in ML algorithms and availability of high-speed computing power. One of the most popular ML class of models is artificial neural networks (ANN). ML is increasingly being used in the design and analysis of composite materials and structures, specifically in the constitutive modeling of composite materials with the focus on greatly accelerating multiscale analyses of composite materials and structures through development of surrogate models. Towards that end, both Python and MATLAB-based neural nets have been developed to predict initial stiffness and fatigue life of an eight-ply symmetric polymer matrix composite laminate. Two types of neural networks, a Multilayer Perceptron (MLP) and a Recurrent Neural Network (RNN), have been developed for both platforms. Results show that the both neural net types can provide an excellent estimate of initial stiffness as well as fatigue life of eight-ply symmetric polymer matrix composite laminate. RNNs are better able to capture the shape of the fatigue curve of a laminate. This tool can be very useful for system level studies to obtain an estimate of desired properties and life of PMC composite laminates. The associated surrogate models could also be used in composite multiscale analyses to replace the actual physics-based calculations at lower scales and thereby significantly increase the computational efficiency of such analyses and thus make multiscale analyses a viable industrial tool for large scale structural problems.
A mathematical model utilizing the internal state variable concept is proposed for predicting the upper bound of the reduced axial stiffnesses in cross-ply laminates with matrix cracks. The axial crack opening displacement is explicitly expressed in terms of the observable axial strain and the undamaged material properties. A crack parameter representing the effect of matrix cracks on the observable axial Young's modulus is calculated for glass/epoxy and graphite/epoxy material systems. The results show that the matrix crack opening displacement and the effective Young's modulus depend not on the crack length, but on its ratio to the crack spacing.
Randomness in woven ceramic matrix composite architecture has been found to cause large variability in stiffness and strength. The inherent voids are an aspect of the architecture that may cause a significant portion of the variability. A study is undertaken to investigate the effects of many voids of random sizes and distributions. Response surface approximations were formulated based on void parameters such as area and length fractions to provide an estimate of the effective stiffness. Obtaining quantitative relationships between the properties of the voids and their effects on stiffness of ceramic matrix composites are of ultimate interest, but the exploratory study presented here starts by first modeling the effects of voids on an isotropic material. Several cases with varying void parameters were modeled which resulted in a large amount of variability of the transverse stiffness and out-of-plane shear stiffness. An investigation into a physical explanation for the stiffness degradation led to the observation that the voids need to be treated as an entity that reduces load bearing capabilities in a space larger than what the void directly occupies through a corrected length fraction or area fraction. This provides explanation as to why void volume fraction is not the only important factor to consider when computing loss of stiffness.
Polyimide composites are being evaluated for use in lightweight support structures designed to preserve the ideal flow geometry within thin shell combustion chambers of future space launch propulsion systems. Principles of lightweight design and innovative manufacturing techniques have yielded a sandwich structure with an outer face sheet of carbon fiber polyimide matrix composite. While the continuous carbon fiber enables laminated skin of high specific stiffness; the polyimide matrix materials ensure that the rigidity and durability is maintained at operation temperatures of 316 C. Significant weight savings over all metal support structures are expected. The protypical structure is the result of ongoing collaboration, between Boeing and NASA-GRC seeking to introduce polyimide composites to the harsh environmental and loads familiar to space launch propulsion systems. Design trade analyses were carried out using relevant closed form solutions, approximations for sandwich beams/panels and finite element analysis. Analyses confirm the significant thermal stresses exist when combining materials whose coefficients of thermal expansion (CTEs) differ by a factor of about 10 for materials such as a polymer composite and metallic structures. The ramifications on design and manufacturing alternatives are reviewed and discussed. Due to stringent durability and safety requirements, serious consideration is being given to the synergistic effects of temperature and mechanical loads. The candidate structure operates at 316 C, about 80% of the glass transition temperature T(sub g). Earlier thermomechanical fatigue (TMF) investigations of chopped fiber polyimide composites made this near to T(sub g), showed that cyclic temperature and stress promoted excessive creep damage and strain accumulation. Here it is important to verify that such response is limited in continuous fiber laminates.
Four lay-ups of continuous fiber silicon carbide (SCS2) fiber/aluminum matrix composites were tested to assess fatigue mechanisms including stiffness loss when cycled below their respective fatigue limits. The lay-ups were 0 (sub 8), 0(sub 2)/ + or - 45 (sub 2s), 0/90 (sub 2s),and 0/ + or 45/90 (subs). The data were compared with predictions from the author's previously published shakedown model which predicts fatigue-induced stiffness loss in metal matrix composites. A fifth lay-up, + or - 45 (sub 2s), was tested to compare shakedown and fatigue limits. The particular batch of silicon-carbide fibers tested in this program had a somewhat lower modulus (340 GPa) than expected and displayed poor bonding to the aluminum matrix. Good agreement was obtained between the stiffness loss model and the test data. The fatigue damage below the fatigue limit was primarily in the form of matrix cracking. The fatigue limit corresponded to the laminate shakedown for the + or - 45 (sub 2s) laminate.
Four lay-ups of continuous fiber silicon carbide (SCS2) fiber/aluminum matrix composites were tested to assess fatigue mechanisms including stiffness loss when cycled below their respective fatigue limits. The lay-ups were 0 (sub 8), 0(sub 2)/ + or - 45 (sub 2s), 0/90 (sub 2s), and 0/ + or - 45/90 (subs). The data were compared with predictions from the author's previously published shakedown model which predicts fatigue-induced stiffness loss in metal matrix composites. A fifth lay-up, + or - 45 (sub 2s), was tested to compare shakedown and fatigue limits. The particular batch of silicon-carbide fibers tested in this program had a somewhat lower modulus (340 GPa) than expected and displayed poor bonding to the aluminum matrix. Good agreement was obtained between the stiffness loss model and the test data. The fatigue damage below the fatigue limit was primarily in the form of matrix cracking. The fatigue limit corresponded to the laminate shakedown for the + or - 45 (sub 2s) laminate.
Deployable smart acoustic resonance particles have dense cores, compliant matrixes surrounding the cores and stiff outer shells surrounding the matrixes. The particles have mechanical stress sensitivities that provide unique band gap shifts when compressed. Groups of similar particles with similar stress sensitivities and similar band gap shifts are added at different times to hydraulic fluids, as circulated through wells with the fluid and pushed into fractures. A plural, sonic monopole well logging tool is lowered into the well to determine locations and depth of fractures and local pressures by distinct resonance of individual groups.
Carbon-fiber-reinforced silicon carbide matrix (C/SiC) composites processed by chemical vapor infiltration are candidate materials for aerospace thermal structures. Carbon fibers can retain properties at very high temperatures, but they are known to have poor oxidation resistance in adverse, high-temperature environments. Nevertheless, the combination of CVI-SiC matrix with higher stiffness and oxidation resistance, the interfacial coating, and additional surface-seal coating provides the necessary protection to the carbon fibers, and makes the material viable for high-temperature space applications operating under harsh environments. Furthermore, C/SiC composites, like other ceramic matrix composites (CMCs), exhibit graceful non-catastrophic failure because of various inherent energy dissipating mechanisms. The material exhibits nonlinearity in deformation even at very low stress levels. This is the result of the severe matrix microcracking present in the as processed composite because of large differences between the coefficients of thermal expansion of the fiber and the matrix. Utilization of these advanced composites in next generation space vehicles will require innovative structural configurations, updated materials, and refined analyses. Structural safety issues for these vehicles are in direct competition with performance and cost. One would have to quantify the uncertainties associated with the design using formal probabilistic methods. Specifically four fundamental aspects on which analyses are based-- (1) loading conditions, (2) material behavior, (3) geometrical configurations, and (4) structural connections between the composite components and baseline structure--are stochastic in nature. A direct way to formally account for uncertainties is to develop probabilistic structural analysis methods where all participating variables are described by appropriate probability density functions. The present work, however, focuses on analyzing the stochastic material behavior of these advanced composites using formal probabilistic analysis methods. Often, some of the desirable property characteristics that allow composites to offer advantages over conventional structural materials (like tailoring of composite properties) and the complexity are in fact responsible for their greater statistical variability and the requirements for more characterization tests. Composite properties are anisotropic as well, having different properties in different directions. This means that characterization of a property such as stiffness--which will vary greatly depending on the orientation of the fiber relative to the direction of the testing--must be repeated for several different directions and loading conditions. The fabrication process for composites also introduces statistical variations in properties and geometry. A composite part is produced in a number of steps, each of which introduces statistical variability. The matrix is usually produced from a combination of raw materials; and the fiber, which has its own set of properties, is often coated or surface treated, introducing yet another source of variability.
Matrix algebra solution, using coordinate systems transformations, to boundary value problem for structures composed of shell elements
A NASTRAN bulk dataset preprocessor was developed to facilitate the integration of filamentary composite laminate properties into composite structural resizing for stiffness requirements. The NASCOMP system generates delta stiffness and delta mass matrices for input to the flutter derivative program. The flutter baseline analysis, derivative calculations, and stiffness and mass matrix updates are controlled by engineer defined processes under an operating system called CBUS. A multi-layered design variable grid system permits high fidelity resizing without excessive computer cost. The NASCOMP system uses ply layup drawings for basic input. The aeroelastic resizing for stiffness capability was used during an actual design exercise.
This paper extends recent work done by the authors in modeling length scale-dependent damage behavior of ceramic matrix composites (CMCs) to include effects of local anisotropy introduced by matrix cracking. This model captures scale-dependent damage initiation and propagation behavior of the brittle matrix by employing internal state variable (ISV) theory within a multiscale modeling framework to obtain damaged matrix stress/strain constitutive relationships at each length scale. The damage ISV captures the effects of matrix cracking and growth by using fracture mechanics and the self-consistent scheme to determine the reduced stiffness of the cracked matrix. Matrix cracks, which activate when stress intensity factors near manufacturing induced cavities exceed the fracture toughness of the material, are assumed to be transversely isotropic in the plane of the crack, and matrix anisotropy occurs when the damaged stiffness tensor is rotated from the crack plane to the global axes. The crack progression and temporal evolution of the damage ISV are governed by fracture mechanics and crack growth kinetics. The model effectively captures first matrix cracking, which is the first significant deviation from linear elasticity. The nonlinear predictive capabilities of the material model are demonstrated for monolithic silicon carbide (SiC) and a 2D woven five-harness satin (5HS) carbon fiber SiC matrix (C/SiC) CMC.
The effect of heat treatment on material properties of Sic/Ti-15-3 was measured by vibration tests. Heat treatment changes the microstructure, which stiffens the matrix and reduces its damping capacity. Test results illustrate how the changes in matrix material affect the stiffness and damping properties of the composite. Damping was found to be more sensitive than stiffness to microstructural changes in the matrix. Effects of heat treatment temperature and exposure time are presented.
Viscoplastic deformation of shale is frequently observed in many subsurface applications. Many studies have suggested that this viscoplastic behavior is anisotropic—specifically, transversely isotropic—and closely linked to the layered composite structure at the microscale. In this work, we develop a two-scale constitutive model for shale in which anisotropic viscoplastic behavior naturally emerges from semianalytical homogenization of a bilayer microstructure. The microstructure is modeled as a composite of soft layers, representing a ductile matrix formed by clay and organics, and hard layers, corresponding to a brittle matrix composed of stiff minerals. This layered microstructure renders the macroscopic behavior anisotropic, even when the individual layers are modeled with isotropic constitutive laws. Using a common correlation between clay and organic content and magnitude of creep, we apply a viscoplastic modified Cam-Clay plasticity model to the soft layers, while treating the hard layers as a linear elastic material to minimize the number of calibration parameters. We then describe the implementation of the proposed model in a standard material update subroutine. The model is validated with laboratory creep data on samples from three gas shale formations. We also demonstrate the computational behavior of the proposed model through simulation of time-dependent borehole closure in a shale formation with different bedding plane directions.
Water-extractable β-glucan with high molar mass (HM) determines health benefits of oat food. Oat β-glucan was extracted by a standardized in vitro digestion method and co-existing water-extractable polysaccharide (WEP) fraction and its HM-arabinoxylan (HM-AX) subfraction were isolated to identify their highly acid-resistant subunit and investigate molecular interactions between constituent polymers. The WEP and HM-AX samples consisted of arabinoxylans (AXs) (74 and 76 %, respectively), however, cellulose constituted the secondary component (6.6 and 12.8 %, respectively). Multi-detection HPSEC along with specific enzymatic hydrolysis of AXs revealed the presence of the HM-xylan domain (16 and 34 %, respectively) built of numerous single- and multi-component populations with random coil and rod-like conformations, which were embedded in a xylan matrix with spherical conformation and controlled the macromolecular shape. Unlike single-component populations, the multi-component ones were resistant to hydrolytic action of AX-hydrolyzing enzymes and represented the subunits that anchor matrix polysaccharides onto cellulose surface. These results indicate that water-extractable β-glucan macromolecule comprises as integral element a cellulose core with two linking populations, HM-xylan and low molar mass glucomannan, which are surrounded by a feruloylated AX-arabinan-arabinogalactan composite and next laminated by β-glucan matrix. Here, the stiff cellulose-xylan backbone is the basis of HM β-glucan organization, controlled by its cellulose-like segments.