Top Banner
FEMCI Workshop - May 18th, 2000 - S. Lienard 1 A Finite Element Modeling Technique for Dynamic Analyses of Preloaded Large Thin Film Membrane Structures Sebastien Lienard - John Johnston NASA - Goddard Space Flight Center May 18th, 2000
16

A Finite Element Modeling Technique for Dynamic … Finite Element Modeling Technique for Dynamic Analyses of Preloaded ... NL Stiffness matrix NASTRAN SOL3 ... Mathcad Frequency domain

May 03, 2018

Download

Documents

votuyen
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: A Finite Element Modeling Technique for Dynamic … Finite Element Modeling Technique for Dynamic Analyses of Preloaded ... NL Stiffness matrix NASTRAN SOL3 ... Mathcad Frequency domain

FEMCI Workshop - May 18th, 2000 - S. Lienard 1

A Finite Element Modeling Techniquefor Dynamic Analyses of Preloaded

Large Thin Film Membrane Structures

Sebastien Lienard - John JohnstonNASA - Goddard Space Flight Center

May 18th, 2000

Page 2: A Finite Element Modeling Technique for Dynamic … Finite Element Modeling Technique for Dynamic Analyses of Preloaded ... NL Stiffness matrix NASTRAN SOL3 ... Mathcad Frequency domain

FEMCI Workshop - May 18th, 2000 - S. Lienard 2

Overview

• Introduction• Sunshield Mechanical Design• Modeling Challenges• Stress Analysis - Wrinkle Formation• Wrinkle Pattern• Cable Network Method• Preloading• Solving Process• ISIS Modeling Environment• Dynamic Results• Modeling Summary• Closing remarks• References

Page 3: A Finite Element Modeling Technique for Dynamic … Finite Element Modeling Technique for Dynamic Analyses of Preloaded ... NL Stiffness matrix NASTRAN SOL3 ... Mathcad Frequency domain

FEMCI Workshop - May 18th, 2000 - S. Lienard 3

Introduction

Inflatable Sunshield In Space Flight Experiment - Overview

Page 4: A Finite Element Modeling Technique for Dynamic … Finite Element Modeling Technique for Dynamic Analyses of Preloaded ... NL Stiffness matrix NASTRAN SOL3 ... Mathcad Frequency domain

FEMCI Workshop - May 18th, 2000 - S. Lienard 4

Sunshield Mechanical Design

• Container– Stores and restrains shield components

during launch phase– Interfaces with deployable mast– Includes inflation system and electronics

• Thin film membranes (4 layers)– Thermal shield– 13microns thick Kapton

• Inflatable booms (4)– Support the membranes

• Ladder structures– Maintain membrane spacing

• Constant force springs– Apply tension to the membranes

Page 5: A Finite Element Modeling Technique for Dynamic … Finite Element Modeling Technique for Dynamic Analyses of Preloaded ... NL Stiffness matrix NASTRAN SOL3 ... Mathcad Frequency domain

FEMCI Workshop - May 18th, 2000 - S. Lienard 5

Modeling Challenges

• Modeling of multi-layer, thin-film sunshield structuresis challenging due to the negligible bending stiffnessexhibited by thin-film membranes.

• Preloading is required to develop out-of-planestiffness in the membranes, and must be accountedfor in dynamic analysis.

• Additionally, the presence of wrinkles alters thestructural behavior of the membranes.

Page 6: A Finite Element Modeling Technique for Dynamic … Finite Element Modeling Technique for Dynamic Analyses of Preloaded ... NL Stiffness matrix NASTRAN SOL3 ... Mathcad Frequency domain

FEMCI Workshop - May 18th, 2000 - S. Lienard 6

Stress Analysis - Wrinkle formation

• Flat half membrane layer with boundary condition ofsymmetry.

Major principal stress Minor principal stress

Positive stresses

Negative stresses:wrinkled regions

Singularities occurduring dynamic

solving process dueto negative stressesin plate/membrane

elements!

Page 7: A Finite Element Modeling Technique for Dynamic … Finite Element Modeling Technique for Dynamic Analyses of Preloaded ... NL Stiffness matrix NASTRAN SOL3 ... Mathcad Frequency domain

FEMCI Workshop - May 18th, 2000 - S. Lienard 7

Wrinkle Pattern

• Visual assumptions:

Thin film membranes subjectto discrete tensile loadsexhibit global wrinklingpatterns along straight linesemanating from theload points.

One tenth scale test article

Page 8: A Finite Element Modeling Technique for Dynamic … Finite Element Modeling Technique for Dynamic Analyses of Preloaded ... NL Stiffness matrix NASTRAN SOL3 ... Mathcad Frequency domain

FEMCI Workshop - May 18th, 2000 - S. Lienard 8

Cable Network Method

• The cable network method was developed specifically to modelpretensioned, wrinkled membrane structures.

– Technique originated by M. Mikulas/U. Colorado-Boulder– Further development by S.Lienard/NASA GSFC

• Based on the established principal that load transfer in wrinkled regionstakes place along wrinkle lines.

• The membrane is meshed with a network of preloaded ‘cables’ mappedto the wrinkle pattern of the structure.

– Longitudinal cables are oriented along the wrinkle pattern (load path).– Transverse cables act as a connection between cables and represent the

mass distribution in the structure.– This approach provides an approximate representation of the load paths

and mass distribution in the structure.

• Method is limited in that it requires prior knowledge of the wrinklepattern to generate the cable network and does not account for in-planeshear effects.

Page 9: A Finite Element Modeling Technique for Dynamic … Finite Element Modeling Technique for Dynamic Analyses of Preloaded ... NL Stiffness matrix NASTRAN SOL3 ... Mathcad Frequency domain

FEMCI Workshop - May 18th, 2000 - S. Lienard 9

Cable Network Method - cont.

• The cable network method has been utilized to model the ISISsunshield and the one-tenth scale NGST sunshield model.

• Validation efforts use the one-tenth scale model ground tests toprovide data for model correlation.– Comparison of cable network model predictions and preliminary

test results shows good agreement.– Further testing is currently underway.

Finite Element Mesh of a membrane layerfor 1/10th Scale Model

Developed using Cable Network Method

Connector (transverse cables)connect longitudinal cables and

provide uniform mass distribution

Longitudinal (inner andouter) cables carry loads

Page 10: A Finite Element Modeling Technique for Dynamic … Finite Element Modeling Technique for Dynamic Analyses of Preloaded ... NL Stiffness matrix NASTRAN SOL3 ... Mathcad Frequency domain

FEMCI Workshop - May 18th, 2000 - S. Lienard 10

Preloading

• Static: Forces– Tension per membrane layer: 1.425N– Compression per boom: 5.7N

• Note: Degree of freedom Tx of the rigid element isconstrained for dynamic analysis

Page 11: A Finite Element Modeling Technique for Dynamic … Finite Element Modeling Technique for Dynamic Analyses of Preloaded ... NL Stiffness matrix NASTRAN SOL3 ... Mathcad Frequency domain

FEMCI Workshop - May 18th, 2000 - S. Lienard 11

Solving Process

• Static: Geometric Nonlinear– The preload produces large stiffness change that has to be

applied using iterative process to generate accuratestrain/stress distribution.

– Export of the updated stiffness matrix representing state ofstrain energy present in the structure.

• Dynamic: Modal, Frequency Response, Transient Response– Dynamic response must be calculated using an accurate

representation of the state of strain energy in the membrane.– Import the updated stiffness matrix from static analysis

Page 12: A Finite Element Modeling Technique for Dynamic … Finite Element Modeling Technique for Dynamic Analyses of Preloaded ... NL Stiffness matrix NASTRAN SOL3 ... Mathcad Frequency domain

FEMCI Workshop - May 18th, 2000 - S. Lienard 12

ISIS Modeling Environment

Preload(Static)

NASTRANSOL 1 - NL

Stiffnessmatrix

NASTRANSOL3

NASTRANSOL11

NASTRANSOL12

Modal results� Natural freq.� Mode shapes� Modal mass participation

DampingDynamic inputs(Time domain)

Frequency domain plots (FRF,PSD)

Time domain results(e.g. acc=f(t) )

Important Modes Amplification Expected measures(Flight simulations)

Matlab orMathcad

Frequency domain plots (FRF,PSD)

Shuttle NoiseDynamic inputs(Freq. domain)

Cable mesh

Shuttlemaneuvers

Page 13: A Finite Element Modeling Technique for Dynamic … Finite Element Modeling Technique for Dynamic Analyses of Preloaded ... NL Stiffness matrix NASTRAN SOL3 ... Mathcad Frequency domain

FEMCI Workshop - May 18th, 2000 - S. Lienard 13

Frequency Response FunctionOut-of-plane direction

-20

-15

-10

-5

0

5

10

15

20

0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0

Frequency (Hz)

Acc

eler

atio

n (m

/s2)

Long boom tip Medium boom tip

Short boom tips

Constant 50mg input along out-of-plane axis from 0 to 10Hz

Dynamic Results

2.17Hz: Long side membrane mode

3.09Hz: Short side membrane mode

3.26Hz: Long boom mode

4.61Hz: Medium boom mode

Page 14: A Finite Element Modeling Technique for Dynamic … Finite Element Modeling Technique for Dynamic Analyses of Preloaded ... NL Stiffness matrix NASTRAN SOL3 ... Mathcad Frequency domain

FEMCI Workshop - May 18th, 2000 - S. Lienard 14

Modeling Summary

• Several techniques for modeling the structural behavior ofpretensioned, wrinkled membrane structures exist.

• The ISIS experiment is modeled using the cable network method.– Most mature technique for modeling wrinkled, pretensioned membranes– Model shows good correlation with preliminary test results. Additional

tests underway to fully validate the technique.

Technique Pros Cons MaturityStandard Elements(membrane or plate,neglecting wrinkles)

§ Easy implementation§ In-plane shear and

thermal effects modeled

§ Wrinkling effects not modeled§ Potential for numerical

problems for dynamic analysis§ Dynamic results due not

convergene as mesh size isrefined

§ Fully developed§ Implemented in NASTRAN§ Model validation needed

Cable NetworkMethod

§ Easy implementation§ Quick solving time

§ Requires knowledge of wrinklegeometry (test required)

§ No in-plane effects modeled

§ Fully developed§ Implemented in NASTRAN§ Model validation underway

Iterative MembraneProperties Method

§ Predicts wrinkle regiongeometry

§ In-plane shear andthermal effects modeled

§ Requires relatively fine meshin wrinkle regions

§ Iterative solution required§ Long solving time

§ Under development§ Dynamics analysis not

developed yet§ Implemented in NASTRAN

(Requires external code)§ Model validation underway

Page 15: A Finite Element Modeling Technique for Dynamic … Finite Element Modeling Technique for Dynamic Analyses of Preloaded ... NL Stiffness matrix NASTRAN SOL3 ... Mathcad Frequency domain

FEMCI Workshop - May 18th, 2000 - S. Lienard 15

Closing Remarks

• Modeling technique– Fully developed in NASTRAN

• Validation underway– Ground test article for T/V testing (3.2x1.4x0.1m)

• Test in air (Early June 2000)• Test in vacuum (Late June 2000)• Model validation/correlation (Summer 2000)

– Flight experiment for testing in Space (11.2x4.9x0.3m)• Flight planned for May 2001

Page 16: A Finite Element Modeling Technique for Dynamic … Finite Element Modeling Technique for Dynamic Analyses of Preloaded ... NL Stiffness matrix NASTRAN SOL3 ... Mathcad Frequency domain

FEMCI Workshop - May 18th, 2000 - S. Lienard 16

References

• ISIS Flight Experiment– Linda Pacini and Michael C. Lou, “Next Generation Space Telescope (NGST)

Pathfinder Experiment: Inflatable Sunshield In Space (ISIS),” October 1999, SAE 1999-01-5517.

– Michael L. Adams, et. al, "Design and Flight Testing of an Inflatable Sunshield for theNext Generation Space Telescope (NGST)," April 3-6, 2000, AIAA-2000-1797.

– Sebastien Lienard, John Johnston, et. al, “Analysis and Ground Testing for Validationof the Inflatable Sunshield In Space (ISIS) Experiment,” 41st AIAA Structures,Structural Dynamics, and Materials Conference, Atlanta, GA, Paper No. AIAA-2000-1638, April 2000.

• Modeling and Analyses of Wrinkled Membrane Structures– Adler, A.L., Mikulas, M.M., and Hedgepeth, J.M., “Static and Dynamic Analysis of

Partially Wrinkled Membrane Structures,” 41st AIAA Structures, Structural Dynamics,and Materials Conference, Atlanta, GA, Paper No. AIAA-2000-1810, April 2000.

– Lienard, S.L., “Characterization of Large Thin Film Membrane Dynamic Behavior withUAI-NASTRAN Finite Element Solver,” SAE Paper 199-01-5518, October 1999.