Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Orthotropic Composite”

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

User-defined Material Model for Thermo-mechanical Progressive Failure Analysis

Previously a user-defined material model for orthotropic bimodulus materials was developed for linear and nonlinear stress analysis of composite structures using either shell or solid finite elements within a nonlinear finite element analysis tool. Extensions of this user-defined material model to thermo-mechanical progressive failure analysis are described, and the required input data are documented. The extensions include providing for temperature-dependent material properties, archival of the elastic strains, and a thermal strain calculation for materials exhibiting a stress-free temperature.

Knight, Norman F., Jr.↗

Modeling creep behavior in ceramic matrix composites

Here, in this work, a three-dimensional viscoplasticity formulation with progressive damage is developed and used to investigate the complex time-dependent constituent load transfer and progressive damage behavior in ceramic matrix composites (CMCs) subjected to creep. The viscoplasticity formulation is based on Hill's orthotropic plastic potential, an associative flow rule, and the Norton-Bailey creep power law with Arrhenius temperature dependence. A fracture mechanics-informed isotropic matrix damage model is used to account for CMC brittle matrix damage initiation and propagation, in which two scalar damage variables capture the effects of matrix porosity as well as matrix property degradation due to matrix crack initiation and propagation. The Curtin progressive fiber damage model is utilized to simulate progressive fiber failure. The creep-damage formulation is subsequently implemented as a constitutive model in the generalized method of cells (GMC) micromechanics formulation to simulate time-dependent deformation and material damage under creep loading conditions. The developed framework is used to simulate creep of single fiber SiC/SiC microcomposites. Simulation results are in excellent agreement with experimental and numerical data available in the literature.

36 MATERIALS SCIENCE↗

Thermal stresses of a wind turbine blade made of orthotropic material

This study is to investigate the thermal stress of a wind turbine blade made of wood composite material. First, the governing partial differential equation on heat conduction is stated, then, a finite element procedure using a variational approach is employed for the solution of the governing equation. Thus, the temperature field throughout the blade is determined. Next, based on the temperature field, a finite element procedure using potential energy approach is applied to determine the thermal stress field. A set of results is obtained through the use of a computer, which is considered to be satisfactory.

Fu, Kuan-Chen↗

Minimum-weight design of laminated composite plates for postbuckling performance

Minimum-weight design of simply-supported, symmetrically laminated, thin, rectangular, specially orthotropic laminated plates for buckling and postbuckling strength is investigated. The postbuckling analysis is based on an Marguerre-type energy method. The failure load of laminates is calculated by the maximum strain failure criterion based on the in-plane strains. Design variables are individual layer thicknesses with specified fiber orientations. Optimization with discrete valued design is achieved by introducing additional penalty terms to the regular pseudoobjective function of sequential unconstrained minimization technique. The proposed optimization technique is applied to the design of rectangular laminates with various aspect ratios loaded by axial compressive loads.

Shin, Dong KU↗

Triangular Element For Analyzing Elasticity Of Laminates

Flat triangular element developed for use in finite-element analyses of stress and strain in laminated plates made of such materials as plywood or advanced fiber/epoxy composite materials. Has multiple layers, each of which can have different isotropic or orthotropic elastic properties. Many such elements used in finite-element mesh to calculate stiffness of plate. Formulation of element straight-forward, and calculation of its stiffness matrix simple and fast.

Martin, C. Wayne↗

Modeling and life prediction methodology for Titanium Matrix Composites subjected to mission profiles

Titanium matrix composites (TMC) are being evaluated as structural materials for elevated temperature applications in future generation hypersonic vehicles. In such applications, TMC components are subjected to complex thermomechanical loading profiles at various elevated temperatures. Therefore, thermomechanical fatigue (TMF) testing, using a simulated mission profile, is essential for evaluation and development of life prediction methodologies. The objective of the research presented in this paper was to evaluate the TMF response of the (0/90)2s SCS-6/Timetal-21S subjected to a generic hypersonic flight profile and its portions with a temperature ranging from -130 C to 816 C. It was found that the composite modulus, prior to rapid degradation, had consistent values for all the profiles tested. A micromechanics based analysis was used to predict the stress-strain response of the laminate and of the constituents in each ply during thermomechanical loading conditions by using only constituent properties as input. The fiber was modeled as elastic with transverse orthotropic and temperature dependent properties. The matrix was modeled using a thermoviscoplastic constitutive relation. In the analysis, the composite modulus degradation was assumed to result from matrix cracking and was modeled by reducing the matrix modulus. Fatigue lives of the composite subjected to the complex generic hypersonic flight profile were well correlated using the predicted stress in 0 degree fibers.

Mirdamadi, M.↗

Multi-Scale Analysis of Effective Mechanical Properties of Porous 3D Woven Composite Materials

This work aims to compute and analyze the mechanical properties of 3D woven composite materials that share the same yarn structure but differ in their matrix porosity and yarn's fiber volume fraction. Assuming linear elastic behavior, NASA's PuMA software was used to characterize the different composites from their constituents' data, averaging their elastic properties in three main directions to obtain their effective orthotropic mechanical properties. This was performed in three steps: firstly, the porous matrix phase was analyzed at the micro-scale; then, these results were used, along with the mechanical properties of the fibers, to model the yarns at the meso-scale; finally, results from the constituents were homogenized and used to model the mechanical behavior of the unit cell at the macro-scale. The 3D woven composite materials selected for this study are made of a mixture of phenolic resin with porosity values between 0 and 0.6, and infiltrated carbon fiber yarns with fiber volume fractions between 0.5 and 0.8. Results from this study show the effect of different matrix porosities and fiber volume fractions on the macro-scale mechanical properties, and were validated through theoretical and semi-empirical expressions found in the literature.

3D woven materials↗

Photoelastic techniques for the complete determination of stresses in composite structures

Three new methods are proposed for determining the individual values of the principal stresses or strains in an orthotropic birefringent model using transmission photoelastic analysis. One of the methods is an extension of the oblique incidence technique in which the model (or the light beam) is rotated about one of the material symmetry axes. In another method, transmission and reflection photoelastic responses are combined, while in the third method small circular holes are drilled at points of interest and the isochromatic fringe order is measured at selected points on the hole boundary. The three methods are applied to an orthotropic disc under diametral compression and results are compared with strain gage data.

Prabhakaran, R.↗

Fracture and crack growth in orthotropic laminates. Part 1: Analysis of a hybrid, unidirectional laminate with damage

The fracture behavior of unifirectional hybrid (buffer strip) composite laminates is studied. Three particular solutions are discussed: (1) broken fibers in a unidirectional half plane; (2) adjoined half planes of different fiber and matrix properties and (3) the solution of two half planes bounding a third distinct region of finite width. This finite width region represents a buffer strip and the potential of this strip to arrest a crack that originates in one of the half planes is investigated. The analysis is based on a materials modeling approach using the classical shear lag assumption to described the stress transfer between fibers. Explicit fiber and matrix properties of the three regions are retained and changes in the laminate behavior as a function of the relative material properties, buffer strip width and initial crack length are discussed.

Goree, J. G.↗

Tension fatigue analysis and life prediction for composite laminates

A tension fatigue life prediction methodology for composite laminates is presented. Tension fatigue tests were conducted on quasi-isotropic and orthotropic glass epoxy, graphite epoxy, and glass/graphite epoxy hybrid laminates. Edge delamination onset data were used to generate plots of strain energy release rate as a function of cycles to delamination onset. These plots were then used along with strain energy release rate analyses of delaminations initiating at matrix cracks to predict local delamination onset. Stiffness loss was measured experimentally to account for the accumulation of matrix cracks and for delamination growth. Fatigue failure was predicted by comparing the increase in global strain resulting from stiffness loss to the decrease in laminate failure strain resulting from delaminations forming at matrix cracks through the laminate thickness. Good agreement between measured and predicted lives indicated that the through-thickness damage accumulation model can accurately describe fatigue failure for laminates where the delamination onset behavior in fatigue is well characterized, and stiffness loss can be monitored in real time to account for damage growth.

Obrien, T. K.↗

Shear buckling of specially orthotropic plates with centrally located cutouts

There is significant industry demand for a method of analyzing the shear buckling of composite plates with cutouts. This method should be able to easily accommodate frequent changes in model design; inflexibility being the major drawback to current finite element methods. The approach taken here is broken into two problems, prebuckling and buckling. To solve for the prebuckling stresses, complex variable equations are used in conjunction with boundary collocation. The least squares approach is utilized to improve the accuracy of the results. The buckling problem is solved using the Ritz method. A product of the Ritz method is a complicated integral equation which is solved using numerical integration. To date, the aforementioned method of determining the prebuckling stresses was verified against infinite plate theory and finite elements. Recent efforts have focused on reducing run time, improving accuracy and adding more boundary condition choices.

Klang, Eric C.↗

Shear buckling of specially orthotropic plates with centrally located cutouts

There is significant industry demand for a method of analyzing the shear buckling of composite plates with cutouts. This method should be able to easily accommodate frequent changes in model design; inflexibility is the major drawback to current finite element methods. The approach taken here is broken into two problems, prebuckling and buckling. To solve for the prebuckling stresses, complex variable equations are used in conjunction with boundary collocation. The least squares approach is utilized to improve the accuracy of the results. The buckling problem is solved using the Ritz method. A product of the Ritz method is a complicated integral equation which is solved using numerical integration. To date, the aforementioned method of determining the prebuckling stresses was verified against infinite plate theory and finite elements.

Owen, Vicki L.↗

Current and projected use of carbon composites in United States aircraft

Carbon composite materials are finding limited use in both civil and military aircraft structures to exploit their weight saving potential for improved performance or fuel efficiency. Since these needs are growing, and a manufacturing cost savings potential is also recognized, wider use of composites in the near future may be expected. Carbon composites generally involve fiber volume fractions in excess of 60 percent in a variety of orthotropic sandwich or solid laminates. Bidirectional woven carbon cloth, common in recently designed structure, may provide higher, more uniform laminate conductivities than tape.

Leonard, R.↗

Application of a trigonometric finite difference procedure to numerical analysis of compressive and shear buckling of orthotropic panels

A numerical analysis developed for the buckling of rectangular orthotropic layered panels under combined shear and compression is described. This analysis uses a central finite difference procedure based on trigonometric functions instead of using the conventional finite differences which are based on polynomial functions. Inasmuch as the buckle mode shape is usually trigonometric in nature, the analysis using trigonometric finite differences can be made to exhibit a much faster convergence rate than that using conventional differences. Also, the trigonometric finite difference procedure leads to difference equations having the same form as conventional finite differences; thereby allowing available conventional finite difference formulations to be converted readily to trigonometric form. For two-dimensional problems, the procedure introduces two numerical parameters into the analysis. Engineering approaches for the selection of these parameters are presented and the analysis procedure is demonstrated by application to several isotropic and orthotropic panel buckling problems. Among these problems is the shear buckling of stiffened isotropic and filamentary composite panels in which the stiffener is broken. Results indicate that a break may degrade the effect of the stiffener to the extent that the panel will not carry much more load than if the stiffener were absent.

Stein, M.↗

Two-dimensional modelling of compressive failure in delaminated laminates

An analytical model is developed to assess the compressive strength criticality of near-surface interlaminar defects in laminated composites. The delaminated region is elliptic in shape, separating a thick isotropic plate from a thin orthotropic layer whose material axes coincide with the ellipse axes. The growth conditions and growth behavior of this defect are studied by breaking the overall problem into an elastic stability problem and a fracture problem. Post-buckling solution for the elliptic section is obtained using the Rayleigh-Ritz method while an energy balance criterion based on a self-similar disbond growth governs the fracture. The parameters controlling the growth or arrest of the delamination damage are identified as the fracture energy, disbond depth and elastic properties of the materials from both sides of the delaminating interface. By varying the degree of material anisotropy relative to the loading axis a range in growth behavior was found including stable or unstable crack growth parallel to or normal to the loading axis.

Chai, H.↗

The initiation and growth of delaminations induced by matrix microcracks in laminated composites

A recent variational mechanics analysis of microcracking damage in cross-ply laminates of the form /(S)/90n/s, where (S) is any orthotropic sublaminate much stiffer than /90n/, has been extended to account for the presence of delaminations emanating from the tips of microcracks in the /90 2n/T sublaminate. The new two-dimensional stress analysis is used to calculate the total strain energy, effective modulus, and longitudinal thermal expansion coefficient for a laminate having microcracks and delaminations. These results are used to calculate the energy release rate for the initiation and growth of a delamination induced by a matrix microcrack. At low crack densities, /(S)/90n/s laminates are expected to fail by microcracking and to show little or no delamination. At some critical crack density, which is a function of laminate structure and material properties, the energy release rate for delamination exceeds that for microcracking and delamination is predicted to dominate over microcracking. A quasi-three-dimensional model is used to predict the propagation of arbitrarily shaped delamination fronts. All predictions agree with experimental observations.

Nairn, J. A.↗

Separation of crack extension modes in orthotropic delamination models

In the analysis of an interface crack between dissimilar elastic materials, the mode of crack extension is typically not unique, due to oscillatory behavior of near-tip stresses and displacements. This behavior currently limits the applicability of interfacial fracture mechanics as a means to predict composite delamination. The Virtual Crack Closure Technique (VCCT) is a method used to extract mode 1 and mode 2 energy release rates from numerical fracture solutions. The mode of crack extension extracted from an oscillatory solution using the VCCT is not unique due to the dependence of mode on the virtual crack extension length, Delta. In this work, a method is presented for using the VCCT to extract Delta-independent crack extension modes for the case of an interface crack between two in-plane orthotropic materials. The method does not involve altering the analysis to eliminate its oscillatory behavior. Instead, it is argued that physically reasonable, Delta-independent modes of crack extension can be extracted from oscillatory solutions. Knowledge of near-tip fields is used to determine the explicit Delta dependence of energy release rate parameters. Energy release rates are then defined that are separated from the oscillatory dependence on Delta. A modified VCCT using these energy release rate definitions is applied to results from finite element analyses, showing that Delta-independent modes of crack extension result. The modified technique has potential as a consistent method for extracting crack extension modes from numerical solutions. The Delta-independent modes extracted using this technique can also serve as guides for testing the convergence of finite element models. Direct applications of this work include the analysis of planar composite delamination problems, where plies or debonded laminates are modeled as in-plane orthotropic materials.

Beuth, Jack L.↗

Fracture and crack growth in orthotropic laminates

A mathematical model based on the classical shear-lag assumptions is used to study the residual strength and fracture behavior of composite laminates with symmetrically placed buffer strips. The laminate is loaded by a uniform remote longitudinal tensile strain and has initial damage in the form of a transverse crack in the parent laminate between buffer strips. The crack growth behavior as a function of material properties, number of buffer-strip plies, spacing, width of buffer strips, longitudinal matrix splitting, and debonding at the interface is studied. Buffer-strip laminates are shown to arrest fracture and increase the residual strengths significantly over those of one material laminates, with S-glass being a more effective buffer strip material than Kevlar in increasing the damage tolerance of graphite/epoxy panels. For a typical graphite/epoxy laminate with S-glass buffer-strips, the residual strength is about 2.4 times the residual strength of an all graphite/epoxy panel with the same crack length. Approximately 50% of this increase is due to the S-glass/epoxy buffer-strips, 40% due to longitudinal splitting of the buffer strip interface and 10% due to bonding.

Goree, James G.↗