Guideline For Evaluation of Finite Elements and Results

NTIS #PB96-1 53077
SSC-387 GUIDELINE FOR EVALUATION OF

FINITE ELEMENTS AND RESULTS
This document has been approved for public release and salq its distribution is unlimited
SHIP
STRUCTURE
COMMITTEE
SHIP STRUCTUR=OMMITTEE The SHIP STRUCTURE COMMITTEE is constituted to prosecute a research program to improve the hull structures of ships and other marine structures by an extension of knowledge pertaining to design, materials, and methods of construction. RADM J. C. Card, USCG (Chairman) Chief, O~ce of Marine Safety, Security and Environmental Protection U.S. Coast Guard Mr. Thomas H. Peirce Marine Research and Development Coordinator Transportation Development Center Transport Canada Mr. Robert McCarthy Director, Survivability and Structural Integrity Group (SEA 03P) Naval Sea Systems Command Mr. Edwin B. Schimler Associate Administrator for Shipbuilding and Technology Development Maritime Administration Dr. Donald Liu Senior Vice President American Bureau of Shipping
Mr. Thomas Connors Acting Director of Engineering (N7) Military Sealift Command
Dr. Ross Grahm Head, Hydronautics Section Defence Research Establishment-Atlantic TFC HNICAL REPRESE NTATIVE
FXEC UTIVE DIRECTOH CDR Stephen E. Sharpe, USCG U, S, Coast Guard
CONTRACTING
OFFICFR
Mr. William J. Siekierka Naval Sea Systems Command
~HIP STRI ICTIJRFSI IRCOMMIT17=F The SHIP STRUCTURE SUBCOMMllTEE acts for the Ship Structure Committee on technical matters by providing technical coordination for determinating the goals and objectives of the program and by evaluating and interpreting the results in terms of structural design, construction, and operation. MARITIME ADMINISTRATION
MILITARY SEALIFT COMMAND Mr. Mr. Mr. Mr. Robert E. Van Jones (Chairman) Rickard A. Anderson Michael W. Touma Jeffrey E, Beach BUREAU OF SHIPPING
u. s,
COAST GUARD
Mr. Frederick Seibold Mr. Richard P. Voelker Mr. Chao H. Lin Dr. Walter M. Maclean NAVAL SEA SYSTEMS Mr. Mr. Mr. Mr. W. Thomas Packard Charles L Null Edward Kadala Allen H. Engle COMMAND
CAPT George Wright Mr. Walter Lincoln Mr. Rubin Sheinberg
AMERICAN Mr. Mr. Mr. Mr.
TRANSPORT Mr. Mr. Mr. Mr.
CANADA
Glenn Ashe John F, ConIon Phillip G, Rynn William Hanzelek
John Grinstead Ian Bayly David L. Stocks Peter llmonin
DEFENCE ~NTIC
R SEARC
BISMNT
Dr. Neil Pegg LCDR Stephen Gibson Dr. Roger Hollingshead Mr. John Porter STRUCTURE SOCIETY OF NAVAL ARCHITECTS MARINE ENGINEERS Dr, William Sandberg CA:::~yC;~E;ORGMl; Dr, William R. Tyson EWLS AND SUBCOMMITl_EE LIAISON MEMBERS OF SCIENCES -
NATIONAL ACADEMY ~D Dr. Robert Sielski NATIONAI ACADEMY
AND
OF SC lNS CE F
URFS
Dr. John Landes WELDING RESEARCH Dr. Martin Prager COUNCIL
u. s. NAVAL ACADEMY Dr. Ramswar Bhattacharyya

U. S. ~~R~NT Dr. C, B. Kim
MARINE ACADEMY
~MFRICAN IRON ANiT STFFI Mr. Alexander D. Wilson C)FFICF ~F NA AL RESEA RCH Dr. Yapa D. S. ;ajapaske
INSTITUT E
U. S, COAST GUAR17 ARA~FMy LCDR Bruce R. Mustain U. S. TECHNI %LAPIVSORY GROU P TO THE INTERNATIONAL STANDARDS ORGANIZATION CAPT Charles Piersall
E OF TE CHNOLOGY CAPT Alan J. Brown
STUDENT MEMBER Mr. Jason Miller Massachusetts Institute of Technology
\..,-
RECENT SHIP STRUCTURE
COMMITTEE
PUBLICATIONS
Ship Structure Committee Publications - A Special Biblioclraphv This bibliography of SSC reports may be downloaded from the internet at: http: //www.starsoftware. com/uscgnmc/nmc/sscl /index.htm SSC-386 SSC-385 SSC-384 SSC-383 SSC-382 Ships Maintenance Project R. Bea, E. Cramer, R. Schulte-Strauthaus, Mayoss, K. Gallion, K. Ma, R. Holzman, L. Demsetz 1995 Hydrodynamic Impact on Displacement Ship Hulls -An the State of the Art J. Daidola, V. Mishkevich 1995 Assessment R.
of
Post-Yield Stren@h of Icebreakirm Ship Structural Members C. DesRochers, J. Crocker, R. Kumar, D. Brennan, B. Dick, S. Lantos @timum Weld-Metal Dexter and M. Ferrell Strength for Hi~h Strength Steel Structures 1995
1995 R.
Reexamination of Desiqn Criteria for Stiffened Plate PaneIs by D. Ghose and N. Nappi 1995 Residual Strenqth of Damaaed Ghose, N. Nappi 1995 Ship Structural B. Bea 1995 lnteq~y Marine Structures by C. Wiernicki, D.
SSC-381 SSC-380 SSC-379 SSC-378
Information
Svstem by R. Schulte-Strathaus,
Improved Ship Hull Structural Details Relative to Fatique by K. Stambaugh, F. Lawrence and S. Dimitriakis 1994 The Role of Human Error in Desire, Marine Structures by R. Bea 1994 Construction and Reliability of
SSC-377
Hull Structural Concepts For Improved Producibility J. Parente, and W. Robinson 1994
by J. Daidola,
SSC-376 SSC-375 SSC-374 SSC-373 SSC-372 SSC-371
Ice Load Impact Study on the NSF R/V Nathanial B. Palmer by J. St. John and P. Minnick 1995 Uncertainty in Strenqth Models for Marine Structures E. Nikolaidis, B. Ayyub, G. White, P. Hess 1994 by O. Hughes, of Desicm and 1994 1994
Effect of Hiqh Strenqth Steels on Strenqth Considerations Construction Details of Shi~ by R. Heyburn and D. Riker Loads and Load Combinations
by A. Mansour and A. Thayamballi A State of the Art Summary by
Maintenance of Marine Structures: S. Hutchinson and R. Bea 1993
Establishment of a Uniform Format for Data Reporting of Structural Material Properties for Reliability Analysis by N. Pussegoda, L. Malik, and A. Dinovitzer 1993 , Underwater Underwater Repair Procedures for Ship Hulls (Fatique and Ductility of Wet Welds) by K. Grubbs and C. Zanis 1993
SSC-370
COMMIITEE
ON MARINE STRUCTURES
Commission
on Engineering
and Technical
Systems
National Academy of Sciences - National Research Council
The COMMllTEE interagency
ON MARINE STRUCTURES
has technical cognizance
over the
Ship Structure Committees of Tennessee,
research program. Knoxville, TN
John Landes, University
Howard M. Bunch, University Dale G. Karr, University Andrew Kendrick, John Niedzwecki,
of Michigan, Ann Arbor, Ml Consultant, Houston, TX
Bruce G. Collipp, Marine Engineering NKF Services,
of Michigan, Ann Arbor, Ml Montreal, Quebec College Station, TX State University, DC DC University Park, PA
Texas A & M University, Pennsylvania
Barbara A. Shaw, Chairman,
Robert Sielski, National Research Council, Washington, Stephen E. Sharpe, Ship Structure Committee,
Washington,
DESIGN WORK GROUP John Niedzwecki, Chairman, Texas A&M University, College Station, TX Bilal Ayyub, University of Maryland, MS CA College Park, MD
Ovide J. Davis, Pascagoula, Maria Celia Ximenes,
Chevron Shipping Co., San Francisco,
MATERIALS Barbara A. Shaw, Chairman, David P. Edmonds, John F. McIntyre, Advanced Harold S. Reemsnyder, Pennsylvania
WORK GROUP
, University Park, PA
State University, OH
Edison Welding Institute, Columbus, Polymer Sciences, Avon, OH Bethlehem
Steel Corp., Bethlehem, Bethlehem, PA
PA
Bruce R. Somers, Lehigh University,
f ,.,.,, ,,, J
-.,,
....+.,. ...J
Member
Agencies:
c
Committee March 1996
Address
Cerreswndence
to:
Arneri&n Bureau of Shipping Defence Research EWblishmentAtiarIttc Maritime Administration Milita Sealifi Command Navti Sea !Jystems Command TransportCanada United States Coast Guard
Ship Structure Committee

An Interagency Advisory 7
Executive Director StructureCommitke U.S. Coast Guard (G-MMS/SSC) 2100 Second Street, S,W. Washin ton, D.C, 20593-0001 Ph:(2027 267-0003 Fex4202) 267-4616
Ship
SSC-387 SR-1364
GUIDELINE
FOR EVALUATION
OF FINITE
ELEMENTS
AND RESULTS
The use of finite element analysis (FEA ) techniques has grown drastically in the last decade. Several structural failures have demonstrated that, if not used properly, the FEA may mislead the with erroneous results. have become so designer The programs user friendly, that engineers little previous with design experience may use them and commit fundamental mistakes, which can result in inadequate strength in the structure. This project intends to reduce the possibility of this human error occurring in design and analysis of ship structures. It in checklists and discussions, provides, a means to review FEA output to ensure the analysis is prepared appropriately for the intended situation. This is no substitute for solid education, enhanced by the experience of the impact of modeling choices on results. The document is to be construed as a guideline to assist the analyst and reviewer in determining deficiencies in an FEA ; it is not a substitute for technical qualifications. This report supports the Coast Guards new program for Prevention Through People which addresses the human error causes of marine casualties.
Rear Admi al, U.S. Coast Chairman, Ship Structure C{ $//y f; 5
w{
.!
(J
-.
,,!? ,
I 5,,
:.:k.
Technical Report Documentation Page

ReportNo. 2. Government AccessionNo. 3. RecipientsCatalogNo.
SSC-387
Title and Subtitle GUIDELINES FOR EVALUATION FINITE ELEMENT ANALYSIS
PB96-153077
5. OF SHIP STRUCTURAL 6. ReportDate December 1995 PerForming OrganizationCode
Author(s) R.1. Basu, K.J. Kirkhope, J. Srinivasan PerForming OrganizationNameand Address MIL Systems Engineering 200-1150 Morrison Drive Ottawa, Ontario, Canada K2H 8S9 2. SponsoringAgencyNameand Address Ship Structure Committee US Coast Guard 2100 Second Street, SW Washington, DC, USA 20593
5.
8.
PerForming OrganizationReportNo. SR-1364
10. Work Unit No. (TRAIS)
11. Contmctor Grant No.
13. Type of Reportand PeriodCovered Final 14. SponsoringAgencyCode G-M
SupplementaryNotes Sponsored by the Ship Structure Committee and its member agencies.
6. Abstract Finite element analysis (FEA) is the most common structural analysis tool in use today. In marine industries, the use of this technique is becoming more widespread in the design, reliability analysis and performance evaluation of ship structures. Users of FEA have considerable freedom in designing the finite element model, exercising it and interpreting the results. Key components of this process include the selection of the computer program, the determination of the loads and boundaty conditions, development of the engineering model, choice of elements and the design of the mesh. A consequence of this freedom is that significant variability in FEA results can be obtained depending on the assumptions and modelling practices adopted by the analyst. A special dificulty is faced by those who have the responsibility for assessing and approving FEAs. Unsatisfactory analysis is not always obvious and the consequences usually will not manifest themselves until the vessel is in service. The individual concerned may not be an expert in FEA, or familiar with the software package used, and will face a dilemma when coming to judge the acceptability, or othetwise, of the results of the FEA. In response to the difficulty faced by those who evaluate FEAs, a systematic and practical methodology has been developed to assess the validity of the FEA results based on the choice of analysis procedure, type of elemenffs, model size, boundary conditions, load application, etc. In support of this methodology, a selection of finite element models that illustrate variations in FEA modelling practices are also presented. Benchmark tests have also been developed which can be used to evaluate the capabilities of FEA software packages to analyze several typical ship structure problems. 7. KeyWords Finite Element Method, Ship Structure, Structural Analysis (Engineering), Quality Assessment
19.
18. DistributionStatementD~s t~ibution
unlimited
Available from: National Technical Information Service Sprirmfield, VA 22161

21. No. of Pages 262 ,, ~,> / L2 22. Price
SecurityClassIf.(of this report) Unclassified
20. SecurityClassification this page) (of Unclassified Reproduction completedpageauthorized of
$36.50Paper $17.50Microfi he
Form DOT F 1700.7 (8-72)
METRIC CONVERSION CARD

Approximate Conversionsto Meaic Measures SX
Swnboi When YouKnow MUMDIVv b
Cofiversions fmm Meiric Measures

Symbol &~
ToFind
: yd
mi
in2 ft2
ydz
mi 2
Oz lb
LENG-ti 2.5 centimeter cm~ centimeters 30 cm ~ d 0.9 meters kilometers 1.6 L= AREA inches S~U~ 6.5 squareentimeters c cmz 0.09 squareetem m mz = square fee& m2 ~ Squareards y 0.8 squareeters m square mfies 2.6 squm kilometers kmz N = a&es 0.4 hictares ha MASS (weight) ounces 28 0.45 E&aIns rounds RQ = hort tons 0.9 meficion [inches feet jug
(2000 lb)
Symbol When YouKnow Multiply by
To Find
Symbol
LENGTH
mm
milhneters
0.04
cm m m km cmz mz Ianz ha
kg t
teaspoons ~;p tabl&poons cubic nches i fl oz fluid ounces

tSQ c pt qt gal ft3 cups pints quarts gallons
VOLUME 5 milWers milliliters 15 milliliters 16 millditers 30 0.24 liters 0.47 liters 0.95 liters
liters
ML
ML mL ML
ML
mL
L L L L
m3 m3
L L L
m3 m3
A:; AREA squme centimeters 0.16 square inches square meters 1.2 square yalds square kilometers 0.4 square miles 2,5 acres hectares (10,000 m2) MASS (weipht) 0.035 ounces 2.2 pounds f&l&ns 1s short tons metric ton (1,000 kg) VOLUME 0.03 fluid ounces milliliters 0.06 CUb~C inches milliliters 2.1 pints liters I.06 quarts liters 0.26 gallons Mm Cubic feet 35 cubic meters 1.3 cubic vards cub~cmeters
TEMPERATURE {exact)
centimeters meters meters kilometem
0.4 3.3
inches inches feet yards miles
in in ft yd mi
in2
ydz
mi2
Oz
lb
fl Oz
in3
pt qt gal
fi3 3 V(j
~dj ,--~, ---F
;:3 cubic eet f cubic ekm m m cubic ards y 0.76 cubic etem TEMPERATURE (exact) @~eS depes subtract 32, Fahrenheit multiply by 513 Cels[us
m~ - ~
degrees Celsius
-40 -20 I 0 32
muldoli ~] 9/5, degtees Fahrenheit

20 I 37 I 60 I 80 I
c F
100
I 40
II
I
160
I
212 water boils
80 98.6
temperature
TABLE OF CONTENTS
PART 1 PROJECT 1,0
OVERVIEW.,....,,,
.,
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ,,, . . . . . . . . . ,,, . . . . . . . . . . . , . . . . . . . . . . . .,, . . . ...8..... . . . . . . . . . . . . . . . . . . . . . . . . .,,. . . . . . . . . . . . . . . . . . . . . . .
1-1 1-1 1-1 1-2 1-2 1-3 1-3 1-4 1-4
INTRODUCTION c,, . . . . . . . . . . . . . . . . Background . . . . . . . . . . . . . . . . . 1.1 Scope . . . . . . . . . . . . . . . . . . . . . 1.2 Overview of Report . . . . . . . . . . . . 1.3 About the Guidelines . . . . . . . . . . . 1.4 1.5 Using the Guidelines . . . . . . . . . . . The Guidelines As Quality Procedures 1.6 Where to Get Further Information . . 1.7
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ,,,,...
PART 2 ASSESSMENT 1,0
METHODOLOGY
FOR FINITE ELEMENT
ANALYSIS
. . . . . . . . . . . .
2-1 2-4 2-4 2-5 2-6 2-7
PRELIMINARY CHECKS.,,,,,, ,,, ,,, ,,, ,s, ,, s,,,,,,,,, Documentation Requirements. . . . . . . . . . . . . . . . . 1.1 Job Specification Requirements . . . . . . . . . . . . . . . . 1.2 Finite Element Analysis Sof-tware Requirements . . . . . 1.3 1.4 Contractor/PersonnelQualification Requirements . . . ENGINEERING MODELCHECKS. . . . 2,1 Analysis Type and Assumptions Geometry Assumptions .iiii 2,2 Material Properties . . . . . . . . 2.3 Stiffness and Mass Properties 2.4 Dynamic DegreesofFreedom 2.5 2.6 Loads and Boundary Conditions . . . . ,,, . . . . . . . . . . . . . . . ,,, . . . . . . . . . . . . . . . . . . . . . . . .,.,,,,.,,..,,.. . . . . . . . . . . . . . . . . . . . . .
. .,.,,.,,. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.0
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. .,,,,,,..,,,,,,,..,,. . .,, . . . . . . ,,, , , ,,, . . . . . . . . . . . . . . . . . . , . . . . ,,, , , ,,, . . . . . . . . . . . . .,,...... . . . . . . . . . . . . , .,.,,,,.. . . . . . . ,,, . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . 2-8 . . 2-8 . . 2-9 . 2-1o . 2-11 . 2-13 . . . . . . . . 2-14 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2-15 2-15 2-16 2-18 2-19 2-20 2-21 2-21 2-22 2-23 2-24 2-25
3,0
FINITE ELEMENT MODELCHECKS ,,, ,,, ,, ., Element Types . . . . . . . . . . . . . . . . . . 3,1 Mesh Design,,,,,.,...,.. . . . . . . . 3,2 Substructures and SubmodeIling ,,, ,,, 3,3 FE Model Loads and Boundary Conditions 3.4 Solution Options and Procedures ,,, ,,, 3,5 FINITE ELEMENT RESULTSCHECKS . General Solution Checks,,,,, 4,1 Postprocessing Methods.,,, 4,2 Displacement Results .,,.iii 4,3 Stress Results . . . . . . . . . . . . 4.4 4.5 Other Results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.0
.- .. ,/
;. ,,
,,
,,
.. ..-
5.0
CONCLUSIONS CHECKS..,,,, . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5.1 FEAResults and Acceptance Criteria . . . . . . . . . . . . . . . . . . . . . . . . . . 5.2 Load Assessment, . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.3
5.4 5.5 PART 3
Strength/ResistanceAssessment
. . . . . . . . . . . . . . . . . . . . . . .,,,..
2-26 2-26 2-27 2-28 2-29 2-30
Accuracy Assessment, . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Overall Assessment, ,,,,... . .,,,,,,,........,,,,,. . . . . . . .
GUIDELINES 1.0
FOR ASSESSING
FINITE ELEMENT
MODELS . . . . . . . . . . . . . . . . . . . .
AND RESULTS . . . . . . . , . . . . . . . , . . . , . . . , . . . . . . . . . . . . , . . .
. . . . . . . .,,.....3-1 . . . . . . . . . . . . . . . . . . . . . . . .
3-1
PRELIMINARY CHECKS.,,,,,. . . . . . . . . . . . 1.1 Documentation Requirements . . . . . . . . . Job Specifica~ion Requirements . . . . . . . . 1,2 Finite Element Software Requirements . . . 1.3 1,4 1.5 Reasons for Using A Particular FEASottware
3-1 3-2 3-3 3-4 3-4 3-5 3-5 3-7 3-7
Package
. . . . , . , , , i , . . . .
Personnel Competence.,,,,, . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1.5.1 Academic and Professional Qualifications . . . . . . . , , , , , . . . . . . . 1,5.2 Training and Experience . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ,, .,,... . . . . . . . . . . . . . . . . . . . Fluid . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . , . . . . . , . . . . . , . . . . . , . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.0
ENGINEERING MODELCHECKS. . . . . . . . Analysis Type and Assumptions . . 281 2.2 Geometry Assumptions..,., . . . 2.3 Material Properties,,,,,,,. . . . . 2,3,1 Composite Materials, ,, . . Stiffness and Mass Properties . . . . 2,4 2.4,1 Mass and Dynamic Problems 2,4.2 Thelnfluence of Surrounding 2.5 Dynamic Degrees of Freedom . . . . 2.6 Loads and Boundary Conditions.. .
...3-8 . . 3-10 . . 3-11 . . 3-12 . . 3-12 . . 3-13 3-15 3-16 3-18 3-19 3-20 3-20 3-21 3-24 3-25 3-26 3-26 3-27 3-28 3-31 3-31 3-32 3-35 3-35
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . , , . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . , . . . . . . . . . . . . . . . . . , . . . . . . . . . . . . . . . . . , . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ...3-13 . . . . ,,. . . . . . . . . . . ...3-22 , , , . . . . . . . . . . . . . . . . . . . . . .. . . . . . . . . . . . . . . . . . .
3.0
FINITE ELEMENT MODEL CHECKS.. . . . . . . . . . Element Types . . . . . . . . . . . . . . . . . . . 3.1 3.1.1 Structural Action to be Modelled .,, Mesh Design, ,,, . . . . . . . . . . . . . . . . . 3.2 3.2.1 Mesh Density, . . . . . . . . . . . . . . 3.2.2 Element Shape Limitations . . . . . . 3.2.3 Mesh Transitions . . . . . . . . . . . . . 3.2.4 Stiffness Ratio of Adjacent Structure 3.3 3,2,5 Miscellaneous Problems,,, . . . . Substructures and Submodelling, . . . . . 3.3,1 Substructuring . . . . . . . . . . . . . 3.3.2 Static Condensation . . . . . . . . . 3.3.3 Two-Stage Analysis, . . . . . . . . . Loads and Boundary Conditions.. . . . . . Minimum Support Conditions . . . 3.4.1 Boundary Conditions for Simulating 3.4.2 Constraints . . . . . . . . . . . . . . . . 3,4.3 Loads - General, ,,, . . . . . . . . . 3.4.4 . . . . . . .
3.4
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Symmetry . . . . . . . . . . . . . . . .
. . . . . . . , . . . .
. . . . . . . . . . . . . . ., . . . .
ii
3.4.5 3.4.6 3.4.7 3.4.8
Loads Loads Loads Loads
Nodal Force and Prescribed Displacement . . Nodal Temperature . . . . . . . . . . . . . . . . . . Face Pressure . . . . . . . . . . . . . . . . . . . . . Edge Loads . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . .
. . . . . . . . .
. . . . . . . . .
. . . . . . . . .
. . . . . . . . .
. . . . . . . . .
. . . . . . . . .
. . . . . . . . .
. . . . . . . . .
. . . . . . . . .
3-35 3-36 3-36 3-39 3-39 3-40 3-40 3-40 3-41 3-41 3-42 3-42 3-42 3-42 3-42 3-42 3-43 3-43 3-43 3-44 3-44 3-45 3-46 3-48 3-48
3.5
3,4.9 Loads -Thermal . . . . .. t.. . . . 3,4i10Gravity and Acceleration.. . . . . Solution Options and Procedures. . . . .. . 3.5.1 Static Analysis . . . . . . . . . . . . . 3.5.2 Dynamic Analysis . . . . . . . . . . . 3i5.3Buckling Analysis . . .. i..
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.0
FINITE ELEMENT RESULTS CHECKS . . . 4.1 General Solution Checks . . . . . . . 4,1.1 Errors & Warnings . . . . . . 4.1.2 Mass and Centre of Gravity 4.1.3 Self-Consistency . . . . . . . 4.1.4 Static Balance . . . . . . . . . 4,1.5 Defaults, , . . . . . . . . . . . 4,2 4.3 4.4 4.1.6 Checklist . . . . . . . . Postprocessing Methods . . . Displacement Results, . . . . Stress Results, . . . . . . . . . . . . . . . . . . . . . . . . .
4.5
484.1 Stress Components . .. iii.... . 4,4.2 Average and Peak Stresses . . . . . Other Results, .,, . . . . . . . . . . . . . . . . 4.5.1 Natural Frequencies and Modes . .
5.0
CONCLUSIONS CHECKS . . . . . . . . . . . i i . . . . 5,1 FEAResults and Acceptance Criteria . . . 5.2 Load Assessment . . . . . . . . . . . . . . . . 5.3 Strength/Resistance Assessment.. . . . . 5,4 Accuracy Assessment . . . . . . . . . . . . . 5.5 Overall Assessment, i i........ . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3-50 3-50 3-51 . . . . . . . . . . . . . . . . . . . . . . . 3-51 . . . . . . . . . . . . . . . . . . . . . . . 3-51 . . . . . . . . . . . . . . . . . . . . . . . 3-52
PART 4 BENCHMARK 1.0 2.0 INTRODUCTION PROBLEMS FOR ASSESSING FEA SOFTWARE .,, , . . . . . . . . . . . 4-1
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ...4-1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..4-4 ..4-4 .. 4-5 . . 4-6
THE BENCHMARK PROBLEMS . . . . . . . . . . . 2.1 BM-l Reinforced Deck Opening,.. . . 2.2 BM-2 Stiffened Panel . . . . . . . . . . . . 2.3 BM-3Vibration isolation System . . . . 2.4 BM-4 Mast Structure . . . . . . . . . . . . 2.5 BM-5Bracket Connection Detail. . . . THE BENCHMARK APPLICATION TEST FEA PROGRAMS
. . . . . . . . . . . . . . . . . . . . . . . . ..4-7 . . . . . . . . . . . . . . . . . . . . . . . ...4-8 4-9 4-9
3.0 4.0
. . . . . . . . . . . . . . . . . . . . . . . . . . . . FEA SOFTWARE . . . . . . . . . .
OF BENCHMARKS
FOR ASSESSING
... Ill
,,/, -.
,,.. j
PART 5 CONCLUSIONS PART 6 REFERENCES Appendix Appendix Appendix Appendix A B C D .,, ,,, ,,, . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . of Finite Element Methodology Practices Models and Results . . . 6-1 A-1 B-1 C-1 D-1 AND RECOMMENDATIONS . . . . . . . . . . . . . . . . . . . . . . . . . . . 5-1
Evaluation
Forms for Assessment of Assessment
Example Application Examples
. . . . . . . . . . . . . . . . . . and Results FEA Software . . . . . . . . . . . . . . . .
of Variations
in FEAModelling
Ship Structure
Benchmark
Problems for Assessing
iv
! .
...3.
J<
ACKNOWLEDGEMENTS
The authors gratefully
acknowledge
the contributions
of Mr. Aaron Dinovitzer
of Fleet
Technologies Limited for his work on the ALGOR benchmarks presented in Appendix D. The authors also wish to thank Canarctic Shipping Limited, and in particular Mr. John McCallum, for permission to use the Arctic tanker example presented in Appendix B.
PART 1 PROJECT OVERVIEW
1.0 1.1
INTRODUCTION
Background Finite element analysis (FEA) isthemost common structural analysis tool in use today. Great strides have been made in theoretical and computational aspects of FEA. This has been accompanied by phenomenal advances in computer technology, both in hardware and software, together with a rapid reduction in the cost of this technology. A consequence of this is a dramatic increase in the affordability of, and accessibility to, finite element technology, In marine industries the use of this technique is becoming more widespread in the design, reliability analysis, and performance evaluation of ship structures, Finite element analysis is a powerful and flexible engineering analysis tool which allows the analyst considerable freedom in designing the finite element model, exercising it and interpreting the results. Key components of this process include the selection of the computer program, the determination of the loads and boundary conditions, development of the mathematical model, choice of elements, and the design of the mesh. Numerous decisions are made by the analyst during this process. Results from FEAs for the same structure performed by different individuals or organizations may differ significantly as a result of differences in the assumptions and modelling procedures employed. Unsatisfactory analysis is not themselves until the vessel is modifications required at this than would be the case if the A special difficulty always obvious and the consequences may not manifest in service, Design changes and any structural stage are generally much more expensive to implement deficiency was discovered earlier. for assessing and
is faced by those who have the responsibility
approving FEAs. The individual concerned may not be an expert in FEA, or familiar with the software package used, and will face a dilemma when coming to judge the acceptability, or otherwise, of the results of the FEA. This may require the evaluator to incur further cost and time in the attempt to assure satisfactory FEA results. In response to the difficulty faced by those who evaluate FEAs a systematic and practical methodology is required to rapidly assess the validity of the FEA results based on the choice of analysis procedure, type of element/s, model sizer boundary conditions, load application etc. In support of this methodology a selection of finite element models that illustrate good modelling practice are also required. In addition benchmark tests are required to allow the validation of new FEA software packages, or packages that have undergone significant modification.
1-1 ,.-. .
1.2
Scope The scope of the guidelines is confined to linear elastic static and dynamic analysis of
surface ship structures using FEA. The treatment of dynamic analysis is limited to natural frequency and mode calculation. The emphasis is on the structural assembly level rather than on local details, or on the total ship, Only FEA of structures composed of isotropic materials is addressed, therefore excluding fibre reinforced plastics and wood, Despite these limitations the guidelines are applicable to the vast majority of ship structure FEAs.
1.3
Overview
of Report in six parts and four appendices as follows:
The report is structured Part 1:
Project Overview This part introduces the document, and provides the background for the methodologies developed for assessing FEAs and FEA software which are described in subsequent Parts. Assessment Methodology for Finite Element Analysis Part 2: This part presents a systematic methodology for assessing FEAs. Appendix A contains forms that can be used for the evaluation process. Appendix B presents an example of a FEA and its evaluation. Part 3: Guidelines for Assessing Finite Element Models and Results
This part provides guidance in support of the methodology presented in Part 2, It is a comprehensive description of good FEA practice. As an aid to the assessment of FEA models and results some FEAs, typical of ship structures, are presented in Appendix C. These examples are designed to illustrate the influence on the results of varying certain model parameters, Benchmark Problems for Assessing FEA Software Part 4: The assessment methodology described in Part 2 includes a requirement that suitable FEA software is used. In support of the assessment new, or significantly modified, FEA should be evaluated in regard to its suitability for ship structure FEA, The benchmark problems and results presented in Part 4 are for this purpose. The benchmark problems are presented in Appendix D.
Conclusions and Recommendations Part 5: This part summarizes observations and insights gained, in the course of this project, into the process of evaluating finite element models and results, and FEA software. Also presented is a summary of where effort should be directed to further improve the methodologies in response to likely future trends in finite element technology, Part 6: Appendix A References Evaluation Forms for Assessment of Finite Element Models and Results
1-2
Appendix Appendix Appendix
B C D
Example Application Examples
of Assessment
Methodology Practices and Results
of Variations
in Fea Modelling for Assessing
Ship Structure
Benchmarks
Fea Software
1.4
About the Guidelines The purpose of the guidelines evaluating finite element presented in this document is to provide a method for
models and results, and also FEA software, to any FEA and it is difficult to assess quality unless the FEA assessment methodology is
There are many attributes
has been comprehensively documented and a systematic applied, This volume presents such a methodology, The methodology 1. is presented in three levels: of attributes
Level 1 comprises
a checklist
of the FEA that need to be evaluated
as part of the assessment 2.
process. breakdown of the checklist provided under as a summary of the Level 2 assessment. finite element modelling practice. The for
Level 2 comprises a more detailed Level 1 can be regarded Level 1. Level 3 contains guidelines
38
on acceptable
guidelines are cross referenced with the Level 2 checklists. During the assessment process the evaluator may, if required, refer to Level 3 guidelines advice. For simple FEAs, an experienced evaluator can probably perform the assessment without referring to Level 2 checklists, The methodology is structured to allow the evaluator to apply the methodology at the appropriate level of detail. The reader is referred to Figure 2-1 i 1 in Part 2 for a graphical overview of the methodology.
In addition to presenting an assessment methodology and suppofiing material, this report presents benchmark problems for assessing the quality of the FEA software and its suitability for ship structural analysis. 1.5 Using the Guidelines The primary audience for these guidelines is evaluators of FEAs, The guidelines assume that the evaluator is trained in ship structural analysis and design, but is not necessarily expert in FEA, Ideally the guidelines would be provided as part of the job of work, statement of requirements, etc.) to the analysts. could then be viewed as acceptance criteria for the work. requirements listed in the guidelines could then be used to required, specifications (or statement The Level 1 and 2 guidelines The documentation stipulate the documentation
1-3
The methodology can be used for conducting reviews which could then be used to provide intermediate and final approvals. For this purpose each of the five areas of a FEA shown in Figure 2-1.1 would be treated as a phase in the project. Reviews could be held at the end of each phase, or less frequently for smaller projects. Depending on the outcome of the review, approval to proceed to the next stage could be given, or, in the case of serious deficiencies rework would be required, Most FEAs will be iterative in character. This applies particularly to analyses performed in support of design tasks. The iterative nature also applies to certain aspects of the analysis itself, Some modelling decisions can only be validated during evaluation of the results. To facilitate this, the methodology is presented as a step-by-step therefore, can accommodate iterations where necessary, 1.6 The Guidelines As Quality Procedures process, and
The guidelines presented in this document incorporate several elements of a quality system as it pertains to FEA and, as such, could be incorporated in an organizations quality system for FEA,
The requirements for such a system have been developed under the direction of the National Agency for Finite Element Methods and Standards (NAFEMS) Quality Assurance Working Group. These requirements for Standardization) are intended 9001. as a supplement to ISO (International 1.7 Where Organization
to Get Further Information provided in the guidelines information is self-contained, is required. The reader is referred to a there may be
While the information circumstances when
more detailed
There are many texts that describe
FEA and theory.
comprehensive bibliography of books and monographs on finite element technology. Besides these texts there are several publications more suited for engineering office use, These include The following reader may wish to consult: guidelines and application-oriented texts that the
13 PIAIJER, J. FL, What Every Engineering Should Kno w About Analysis, Marcel Dekker, Inc., New York, 1988,
Finite Element
MEYER, C. (Ed.), Finite Element Idealization for Linear Elastic Static and D ynamic Analysis of Structures in Engineering Practice, American Society of Civil Engineers, New York, 1987. . NAFEMS, Guidelines to Finite Element Practice, National Agency for Finite
Validation
Quality System Supplement to ISO 9001 Relating to Finite Element Analysis in the Design and of Engineering Products, Ref: ROOI 3, NAFEMS, East Kilbride, Glasgow, UK, 1990.
2 A, K. Noor, Bibliography of books and monographs on finite element technology, Applied mechanics Review, Vol. 44, No. 6, June 1991. 1-4
.,.,., ,
Element Methods and Standards, Glasgowr UK, August 1984. . STEELE, J. E., Applied 1989.
National
Engineering
Laboratory,
East Kilbride,
Finite Element Modelling,
Marcel
Dekker,
Inc., New York,
1-5
ASSESSMENT
PART 2 METHODOLOGY FOR FINITE ELEMENT
ANALYSIS
The methodology developed for evaluating finite element analyses of ship structures is The evaluation is carried out at two levels conducted in parallel. presented in Figure 2-1,1. The highest level (Level 1 ) addresses general aspects of the finite element analysis (FEA) broken down into five main areas: 1, 2, 3. 4. 5. Preliminary Checks, Engineering Model Checks, Finite Element Model Checks, Finite Element Results Checks, Conclusions Checks. Evaluation of (Level 2) aspects
and
These are identified in each of the five main boxes shown in Figure 2-1.1. each of these general aspects in &urn requires that certain related detailed
be checked, The Level 2 aspects to be checked are listed within the main boxes and are presented in detail in separate tables that form the core of the evaluation process. The Level 2 tables contain many detailed questions regarding specific aspects of the FEA. The way the methodology is intended to be used is described as follows. The evaluator will begin by assembling the analysis documentation and perhaps computer files of the finite element (FE) model and results. The evaluation then begins with the Preliminary Checks contained in Box 1 of Figure 2-1.1, The first of the preliminary checks involve assessment of the contents of the analysis documentation (1,1 Documentation). To perform this assessment, the evaluator refers to the table entitled l. 1 Documentation Requirements. This table asks the evaluator to check that the documentation contains information that is essential for the FEA evaluation. The table also refers the evaluator to Part 3 Section 1.1 of the guideline should further explanation or guidance be necessary. If an item is contained in the documentation, the evaluator should place a check mark (d) in the corresponding box under the Resu/t column. If an item is not included with the documentation, the evaluator may enter a cross (X) in the result box, or NA (for Not Applicable), or ? (for further information required). After checking off each item in the table, the evaluator is asked to answer Question 1.1 at the bottom of the page. The answer will be based on the evaluators assessment of each item listed in the table in Section 2-1 i 1, The evaluator should place the answer in the result box to the right of the question, and then transfer it to the corresponding result box in Figure 2-1.1. It is suggested that the same format of answers be used (eg. #, X, A!A, or ?). The table in Section 2-1,1 also includes spaces for the evaluator to enter comments regarding specific At the end of the evaluation process, and overall aspects of the documentation contents. these comments will provide the evaluator with reminders of specific aspects of the FEA that were good, bad, or not explained well. The evaluator may refer to these comments to seek further explanation or clarification from the contractor / analyst (perhaps at a review meeting, or during a telephone conversation) before deciding on the final acceptability of the FEA. Having completed the first of the preliminary checks, the evaluator then proceeds to the second set of checks entitled 1.2 Job Specification Requirements, In a manner similar
2-1
to the previous checks, the evaluator will refer to the table in Section 2-1.2 and perform checks 1 .2.1 to 1 .2.7 which are aimed at verifying that the analysis covers the main requirements and objectives of the job specification (or contract, or statement of work, etc.). Based on the results of these checks, the evaluator should answer Question 1,2 and This procedure is repeated for the other Preliminary Checks enter the result in Figure 2-1.1. (i.e. 1,3 FEA Software, and 1,4 Contractor/ Analyst Qualifications). Having answered all of the Level 2 questions for Part 1 Preliminary Checks and entered the results into the appropriate box in Figure 2-1.1, the evaluator is then asked the question Preliminary checks are acceptable?. The evaluator should check the Yes or No box below this question based on an assessment of the results of the Level 2 preliminary checks. If the answer is NO, then the FEA is very likely not acceptable since it does not meet certain basic requirements. The evaluator may therefore choose to terminate the evaluation at this point. Otherwise, the answer is Yes and the FEA has passed the preliminary checks and the evaluator is instructed to proceed to the next major aspect of the evaluation, entitled 2 - Engineering Model Checks.
The evaluation process continues as described above for each of the five main areas identified in Figure 2-1.1. At the end of this process, the evaluator will check either the oval box entitled FE analysis is Acceptable, or the one entitled FE analysis is Not Acceptable Ideally, depending on the outcome of the assessment checks, methodology
at the start of the job, the contractor
would be given the assessment
as part of the job specification, This will encourage self-checking provided by the contractor to the customer is complete.
and ensure that the data
A set of blank forms is provided in Appendix A. The forms are in a format that can be used in an engineering office environment. The forms are based on the forms in Part 2 with additional space provided for project information,
2-2
1- Prellmlnafy Cheaks
1.1 Documentation Performthese checksto mssure thatthe analyaisdocumentation,ob j speclfmstion, EA sotlware,and F wntmator I analystqualfi=tion requirementshave been addressed. 1.2 Job Specification 1.3 Flnlte Element Analysis Software 1.4 Contractor /Analyst Qualifications
Result Preliminarychecks are acceptable? Yes No A
I
&
~y Result 2- Engineering Model Checks 2.1 Analyaia Type&Assumptions 2.2 Geometry 2.3 Material Pmpmtiaa 2.4 Stiffness & Maaa Properties 2.5 Dynamic Degrees of Freedom 2.6 Loads & Boundary Conditions Performfhse checksto enaurathat the assumptionsused to developthe engineeringmodel of me problemare reasonable. Engineeringmodel is accspfable7 Yes No
No~
yea~
I
n
kNo
3- Finite Element Model Checks
..
Performthese checksto ensurethat the finiteelement model Is an adequate interpretation me engineeringmodel, of
3.2 Mesh Design 3.3 Sub-tructuras and Submodels 3.4 FE Loads& Boundary Condltfons 3.5 FE Solutlon Options & Procedures w I 1Yes No
q 4. FhshsElement Analysis Resulk Cheoke 4.1 General Solution Checks Perfonm these checks to ensurethat the finiteelement resultsare calculated,pmcassed and presentedin a mannerconsistentwrn me analysis requirements. 4.2 Peat Processing Methods 4.3 Displacement Results 4.4 Stres- Raaults 4.6 Dther Results
Reautt
53
Yes No Result \ Conduaionaof me analysiaare acceptable ? Yes # I I No
Finiteelement resultsare acceptable ?
*
5. cnrlcl iion$
Checks
Performthese checks to ensure that adequate considerationofthe Ioada, absngth,awaptanca titetia, FE model,and resultsaccurecyare includedin arrivingat me wndusions fromme finiteelement analysis,
5.1 FE Rosulk & Acceptance Criteria 5.2 Loads Assessment 5.3 Strength I Reslstence Aeeessmsmt 5.4 Accurecy Assessment S.5 Overall Assessment
FIGURE 2-1.1
FE analysis is
FE analysis is
Acceptable
o Methodology Chart
Not Acceptable
Overall Evaluation
2-3
1.0 1.1
PRELIMINARY Documentation
CHECKS Requirements comprehensive assessment submitted, Refer to Guideline Section 3-1.1 of a FEA, cenain essential information must
In order to perform be provided
in the documentation
Finite Element Analysis Assessment
Check
Result
Comments
1.1,1
Has the following information been provided in the FEA documentation? a) b) c) d) e) f) g) h) i) j) k) 1) Objectives and scope of the analysis. and acceptance criteria.
Analysis requirements FEA software used.
II
Description of physical problem. Description of engineering model, Type of analysis, System of units, Coordinate axis systems,
Description of FEA model, Plots of full FEA model and local details. Element types and degrees of freedom per node. Material properties,
m) Element properties (stiffness & mass properties). n) o) p) q) r) FE loads and boundary conditions. Description and presentation Assessment of the FEA results,
of accuracy of the FEA results,
Conclusions of the analysis. List of references. 1 Result I
Based on the above checks answer Question 1.1 and enter result in Figure 1.0. 1.1 Is the level of documentation sufficient to perform an assessment of the FEA?
Comments
2-4
1.2
Job Specification
Requirements
Perform these checks to ensure that the analysis addresses the objectives, scope, requirements and intent of the job specification (eg. contract document, work specification, statement of work, etc.). Refer To Guideline Section 3-1.2
Finite Element Assessment 1.2.1 Is the job specification referenced 1.2.2
Check and
Result
Comments
identified
in the analysis documentation? and scope of the analysis with 3-1,2 3-1.2
Are the objectives
clearly stated and are they consistent those of the job specification? 1.2.3 Are the analysis requirements and are they consistent job specification? 1.2.4 If certain requirements
clearly stated
with those of the of the job (such 3-1.2
specification have not been addressed as certain load cases), has adequate justification been given? 1.2.5 Are the design / acceptance
criteria clearly with those of for using
3-1.2
stated and are they consistent the job specification? 1.2.6 Is there reasonable justification FEA for this problem?
3-1,2
1.2.7
Has advantage been taken of any previous experimental, analytical, or numerical works that are relevant to this problem?
3-1.2
Based on the above checks answer Question 1.2 and enter result in F[qure 1.0.
I I
Result
I 1.2
Does the analysis address the job specification
requirements?
Comments
2-5
,, , J ,.\=.,,.
.
;) ,.,
1.3
Finite Element
Analysis
Software
Requirements for ship
The FEA software should meet certain minimum standards to be considered acceptable structural analysis applications.
Check
Refer To Guideline Section 3-1,3
Result
Comments
1.3.1
Is the FEA software on the list of approved programs for ship structural analysis applications?
If the answer to Check 1.3.1 is Y, you may skip Checks 1.3.2 and 1.3.3. 1.3.2 Are the capabilities and limitations of the FEA software used to perform the required analysis stated in the analysis documentation? 1.3.3 Is evidence of this capability documented and available for review (egi verification manual, results of ship structure FEA benchmark tests, previous approved FEA of similar problems)? Does the vendor of the FEA software have a quality system to ensure that appropriate standards are maintained in software development and maintenance. 3-1,3 3-1.4
1.3.4
Based on the above checks answer Question 1.3 and enter result in Fiaure 1.0. 1.3 Is the FEA software qualified to perform the required analysis?
Comments
NOTE: Part 4 of this report presents benchmark problems for the purpose of assessing the quality and suitability of FEA software for performing ship structural analysis. On its own, successful performance of the candidate FEA software in exercising the benchmark problems is not sufficient evidence of the quality and suitability of the software. The assessor should, in addition, be able to answer the other questions in the table above affirmatively.
2-6
1,4
Contractor
/ Personnel
Qualification
Requirements
The contractor and contractor personnel should possess certain minimum qualifications for performing ship structure FEA, In addition, the contractor should have a Quality Assurance
(QA) system in place to ensure that proper management,

procedures have been applied in the analysis.
administrative
and checking
Finite Element Assessment Check
Refer To Guideline Section 3-1.5
Result
Comments
1.4,1
Do the contractor personnel have adequate academic training and experience qualifications to perform finite element analysis?
1.4.2 Do the contractor personnel have adequate engineering experience qualifications for performing ship structural design or analysis? 1.4.3 Do the contractor and contractor personnel have adequate professional certification qualifications? 1.4.4 Does the contractor have a working system of Quality Assurance (QA) procedures and checks that are satisfactory for the requirement?
3-1.5
3-1.5
3-1.5
1.4.5 Do the contractor personnel have adequate experience with the FEA software used for the analysis?
3-1.5
Based on the above checks answer Question 1.4 and enter result in Fiaure 1.0. I 1.4 Is the contractor adequately qualified for performing ship structure FEA?
m II
Comments
2-7
2.0 2.1
ENGINEERING Analysis Type
MODEL
CHECKS
and Assumptions
Perform these checks to ensure that the assumptions used in developing the engineering model or idealization of the physical problem are adequate. Refer To Guideline Section 3-2,1
Finite Element Analysis Assessment Check 2.1.1 Does the engineering model employ enough dimensions and freedoms to describe the structural behaviour (eg, 1-D, 2-D, or 3-D)? Does the engineering model address the appropriate scale of response for the problem (eg. global, intermediate, 2,1.3 or local response)? Is the type of analysis appropriate for the type of response and loading of interest (eg. linear, static, dynamic, buckling analysis)? Does the engineering model address all the required results parameters (eg: stress, displacement, frequency, buckling load)? Are all assumptions affecting the choice of engineering model and analysis type justified (watch for non-standard assumptions)? Is the level of detail, accuracy or conservatism of the engineering model appropriate for the criticality of the analysis and type of problem? Does the analysis employ a consistent set of units? Does the analysis employ a consistent global coordinate axis system?
Result
Comments
2.1.2
3-2.1
3-2.1
2.1.4
3-2.1
2.1.5
3-2.1
2.1.6
3-2,1
2.1.7
3-2.1
2.1.8
3-2.1
Based on the above checks answer Question 2.1 and enter result in Figure 1.0. Are the assumptions of the type of analysis and engineering model acceptable? Comments
Result
2-8 -------$ .;
2.2
Geometry
Assumptions have been followed for
Perform the following checks to ensure that correct procedures defining the geometric properties of the structure.
Check
Result
Comments
2.2.1
Does the extent of the model geometry cepture the main structural actions, load paths, and response parameters of interest? Are correct assumptions used to reduce the extent of model geometry (eg. symmetry, boundary conditions at changes in stiffness)? Will the unmodelled structure (ie. outside the boundaries of the engineering model) have an acceptably small influence on the results? Are the effects of geometric simplifications (such as omitting local details, cut-outs, etc. ) on the accuracy of the analysis acceptable ?
2.2.2
3-2,2
2.2.3
3-2.2
2.2.4
3-2,2
2.2.5
For local detail models, have the aims of St. Venantts principle been satisfied?
3-2.2
2.2.6
Do the dimensions defining the engineering model geometry adequately correspond to the dimensions of the structure? For buckling analysis, does the geometry adequately account for discontinuities and imperfections affecting buckling capacity?
3-2.2
2.2.7
3-2,2
Based on the above checks answer Question 2.2 and enter result in Figure 7.0. 2.2 Are the geometry assumptions in the engineering model acceptable?
Result
Comments
2-9
...
,>
2.3
Material
Properties for defining
Perform the following checks to ensure that correct procedures have been followed the material properties of the structure.
Check
2.3.1
Are all materials of structural importance to the problem accounted for in the engineering model? Are the assumed behaviors valid for each material (egi linear elastic, isotropic, anisot,ropic, orthotropic) ? Are the required material parameters defined for the type of analysis (eg. E, v, etc.)?
2.3.2
3-2.3
2.3.3
3-2.3
2.3.4
Are orthotropic and / or layered properties defined correctly for non-isotropic materials such as wood and composites? Are orthotropic properties defined correctly where material orthotropy is used to simulate structural orthotropy (eg. stiffened panels)? If strain rate effects are expected to be significant for this problem, are they accounted for in the material properties data? Are the values of the materials properties data traceable to an acceptable source or reference (eg. handbook, mill certificate, coupon tests)? Are the units for the materials properties data consistent with the system of units adopted for other Darts of the analvsis?
3-2.3
2.3.5
3-2.3
2.3.6
3-2.3
2.3.7
3-2.3
2.3.8
3-2.3
Based on the above checks answer Question 2.3 and enter result in Figure 1.0. 2.3 Are the assumptions and data defining the material properties acceptable?
1 Result I
Comments
210
.,.
1 -, ,,,
2.4
Stiffness
and Mass
Properties
have been followed for
Perform the following checks to ensure that correct procedures defining the stiffness and mass properties of the structure. Refer To Guideline Section 3-2,4
Finite Element Analysis Assessment Check 2.4.1 Are all components that have significant effect on the stiffness of the structure accounted for in the engineering model ? Are the assumed stiffness behaviors valid for each structural component (eg. linear, membrane, bending, shear, torsion, etc.)? Are the required stiffness parameters defined for each component, eg. : Truss members - A - A, IW, IZZ,other Beams, bars - t (uniform or varying) Plates, shells - K (axial or rotational) Springs Do the section properties of stiffeners (where modelled with beams) include correct allowances for the effective plate widths? If torsion flexibility is expected to be important, are torsion flexibility parameters correctly defined for beam sections? If shear flexibility is expected to be important, are shear flexibility parameters correctly defined for beam and/or plate elements?
Result
Comments
2.4.2
3-2.4
2.4.3
3-2,4
2.4.4
3-2.4
2.4.5
3-2,4
2.4.6
3-2.4
If mass or inetiial effects are not applicableto this problem. proceed to Check 2.4.13 on the following page. 2.4.8 Are all components that have significant effect on the mass of the structure accounted for in the engineering model? Have material properties data for density been defined (see also Check 2.3.3)? Has the added mass of entrained water been adequately accounted for with structure partially or totally submerged under water? Are lumped mass representations of structural mass and / or equipment correctly consolidated and located? If rotational inertia is expected to be important, are mass moments of inertia properties correctly defined for masses? 2-11 3-2,4
2,4.9
3-2.4
2.4.10
3-2.4
2.4.11
3-2,4
2.4.12
3-2.4
Finite Element Analysis Assessment 2.4.13
Check
Result
Comments
Are the values of the stiffness and mass properties data supported by acceptable calculations and / or references? If relevant, has fluid-structure interaction been accounted for? Has the added mass been included in the model? Are the units for the stiffness and mass properties data consistent with the system of units for other parts of the analysis?
2.4.14
3-2,4
2.4.15
3-2.4
Based on the above checks answer Question 2.4 and enter result in Figure 1.0. 2.4 Are the assumptions and data defining stiffness and mass properties acceptable? Comments
Result
2-12
. ....
2.5
Dynamic In dynamic
Degrees
of Freedom it is often desirable or necessary to reduce the size of the problem by
analyses,
reducing the number of dynamic degrees of freedom (dof). Perform these checks to ensure that the correct procedures have been followed for selecting dynamic degrees of freedom. If the analysis is not a reduced dynamic analysis, you may proceed directly to Part 2.6.
Finite Element Analysis Assessment Check
Result
Comments
2.5.1
Are dynamic dof defined in enough directions to model the anticipated dynamic response behaviour of the structure? Are the number of dynamic dof at least three times the highest mode required (eg. if 30 modes required, need at least 90 dof)? Are the dynamic dof located where the highest modal displacements are anticipated?
2.5.2
3-2,5
2.5.3
3-2.5
2.5.4
Are the dynamic dof located where the highest mass-to-stiffness ratios occur for the structure? Are dynamic dof located at points where forces or seismic inputs are to be applied for dynamic response analyses? Are the number of dynamic dof such that at least 90% of the structural mass is accounted for in the reduced model in each direction?
3-2,5
2.5.5
3-2.5
2.5.6
3-2.5
Based on the above checks answer Question 2.4 and enter result in Figure 7.0. 2.5 Are the assumptions i Comments and data defining dynamic degrees of freedom acceptable?
Result
2-13
,. ..
2.6
Loads and Boundary
Conditions have been followed for defining
Perform the following checks to ensure that correct procedures the loads and boundary conditions of the problem.
Check
Result
Comments
2.6.1
Are all required loadings / load cases accounted for, and has sufficient justification been provided for omitting certain loadings? Are the loading assumptions stated clearly and are they justified?
2.6.2
3-2.6
2.6.3
Has an assessment been made of the accuracy and / or conservatism of the loads?
3-2,6
2,6.4
Are the procedures for combining loads / load cases (eg. superposition) adequately described and are they justified? Have the boundary conditions assumptions been stated clearly and are they justified?
3-2.6
2.6.5
3-2.6
2.6.6
Do the boundary conditions adequately the anticipated structural behaviour?
reflect
3-2.6
2.6,7
Has an assessment been made of the accuracy of the boundary conditions, and if they provide a lower or upper bound solution?
3-2.6
Based on the above checks answer Question 2.6 and enter result in Figure 7.0. I 2.6 Are the assumptions and data defining loads and boundary conditions reasonable? Comments
! I
Result I I
2-14
./---.. k-
3.0 3.1
FINITE Element
ELEMENT Types
MODEL
CHECKS
Perform these checks to ensure that the correct types of elements have been used to model the problem. To assist in this process a checklist is provided in Part 3, Section 3, paragraph 3.1, Refer To Guideline Section 3-3.1
Finite Element Analysis Assessment Check 3.1.1 Are all of the different types of elements used in the FEA model identified and referenced in the analysis documentation? Are the element types available in the FEA software used appropriate to ship structural analysis? Do the element types support the kind of analysis, geometry, materials, and loads that are of importance for this problem? [f required, do the selected beam element types include capabilities to model transverse shear and / or torsional flexibility behaviour? If required, do the selected beam element types include capabilities to model tapered, off-set or unsymmetric section properties? If required, do the selected beam element types include capabilities for nodal dof end releases (eg. to model partial pinned joints)? If required, do the selected plate element types include capabilities to model out-ofplane loads and bending behaviour? [f required, do the selected plate element types include capabilities to model transverse shear behaviour (ie, thick plate behavior)? If the model is 2-D, are the selected element types (or options) correct for plane stress or plane strain (whichever case applies)? If required, can the selected element types model curved surfaces or boundaries to an acceptable level of accuracy?
Result
Comments
3.1.2
3-3,1
3.1.3
3-3.1
3.1.4
3-3,1
3.1.5
3-3,1
3.1.6
3-3,1
3.1.7
3-3.1
3.1.8
3-3.1
3.1.9
3-3.1
3.1.10
3-3.1
Based on the above checks answer Question 3.1 and enter result in FIqure 1.0. I 3.1 Ara the types of elements used in the FEA model acceptable?
I II
Result
Comments
2-15
3.2
Mesh
Design
As the finite element method is essentially a piece-wise approximation technique, the accuracy is very largely dependant on the mesh design, Perform the following checks to ensure that the finite element mesh is acceptable, Refer To Guideline Section 3-3,2
Check
Result
Comments
Does the mesh design adequately reflect the geometry of the problem (eg. overall geometry, stiffener locations, details, etc.)? Does the mesh design adequately reflect the anticipated structural response (eg, stress gradients, deflections, mode shapes)? Are nodes and elements correctly located for applying loads, support and boundary constraints, and connections to other parts? Does the analysis documentation state or show that there are no illegal elements in the model (ie. no element errors or warnings)? Are the element shapes in the areas of interest acceptable for the types element used and degree of accuracy required? Are mesh transitions from coarse regions to areas of refinement acceptably gradual? Are element aspect ratios acceptable, particularly near and at the areas of interest? Are element taper or skew angles acceptable, particularly near and at the areas of interest? If flat shell elements are used to model curved surfaces, are the curve angles < 10 for stresses, or < 15 for displacement results? If flat shell elements are used for double or tapered curve surfaces, is warping avoided (eg. small curve angles, use of triangles)? Is the mesh free of unintentional gaps or cracks, overlapping or missing elements? Is proper node continuity maintained between adjacent elements (also continuity between beam and plate elements in stiffened panels)?
3.2.2
3-3,2
3.2.3
3-3.2
3.2.4
3-3.2
3.2.5
3-3.2
3.2.6
3-3.2
3.2.7
3-3.2
3.2.8
3-3,2
3.2.9
3-3.2
3.2.10
3-3.2
3.2.11
3-3.2
3.2.12
3-3.2
2-16
..
Check
Result
Comments
3.2.13
Are the orientations of the beam element axes correct for the defined section properties?
3.2.14
Are differences in rotational dof / moment continuity for different element types accounted for (eg, beam joining solid)? Are the outward normals for plate / shell elements of a surface in the same direction?
3-3.2
3.2.15
3-3.2
Based on the above checks answer Question 3.2 and enter result in Fiaure 1.0. 3.2 Is the design of the finite element mesh acceptable?
Comments
2-17
3.3
Substructures
and Submodelling
Substructuring or submodelling techniques may be employed to reduce the size of the problem for computing and / or to take advantage of repetitive geometry in the structure. Perform the following checks to ensure that the acceptable procedures have been followed. Refer To Guideline Section 3-3.3
Check
Result
Comments
Is the overall substructure or submodelling scheme or procedure adequately described in the analysis documentation? Are all individual substructure models, global models and refined submodels identified and described in the analysis documentation? Are the master nodes located correctly and are the freedoms compatible for linking the substructures? Are the master nodes located correctly for application of loads and boundary conditions upon assembly of the overall model? Are loads and boundary conditions applied at the substructure level consistent with those of the overall model? Does the boundary of the refined submodel match the boundary of coarse elements / nodes in the global model at the region of interest? Is the boundary for the submodel at a region of relatively low stress gradient or sufficiently far away from the area of primary interest? Does the refined submodel correctly employ forces and / or displacements from the coarse model as boundary conditions? Does the submodel include all other loads applied to the global model (eg. surface pressure, acceleration loads, etc.)? Have stiffness differences between the coarse global mesh and refined submodel mesh been adequately accounted for?
3.3.2
3-3,3
3.3,3
3-3.3
3.3.4
3-3.3
3.3.5
3-3,3
3.3.6
3-3.3
3.3.7
3-3,3
3.3,8
3-3,3
3.3.9
3-3.3
3.3.10
3-3,3
Based on the above checks answer Question 3.3 and enter result in Fiqure 1.0. I 3.3 Are the substructuring or submodelling procedures acceptable~
m
II
Comments
2-1s
3.5
Solution
Options
and Procedures techniques or
Perform the following checks to ensure that correct solution options, procedures have been used for the finite element model.
Check
Result
Comments
3.5.1
Have any special solution options and procedures been used and, if so, have they been documented? If non-standard options been invoked have they been documented and the reasons for their use been explained? If the problem is a dynamic analysis is the method for eigenvalue and mode extraction atmropriate?
3.5.2
3-3.5
3.5.3
3-3.5
Based on the above checks answer Question 3.5 and enter result in Fiaure 7.0. 3.5 Are the solution options and procedures followed for the FEA acceptable? Comments
G I
2-20
3.4
FE Model
Loads and Boundary
Conditions for
Perform the following checks to ensure that correct procedures have been followed defining the loads and boundary conditions of the finite element model.
1
Check
Result
Comments
3.4.1
Are point load forces applied at the correct node locations on the structure and are they the correct units, magnitude, and direction? Are distributed loads applied at the correct locations on the structure and are they the correct units, magnitude and direction? Are surface pressure loads applied at the correct locations on the structure and are they the correct units, magnitude and direction? Are translational accelerations in the correct units, and do they have the correct magnitude and direction? Are rotational accelerations the correct units, magnitude and direction and about the correct centre of rotation? Are prescribed displacements applied at the correct locations on the structure and are they the correct units, magnitude and direction. Are the displacement boundary conditions applied at the correct node locations?
3.4,2
3-3.4
3.4.3
3-3,4
3.4.4
3-3.4
3.4.5
3-3.4
3.4,6
3-3.4
3.4.7
3-3.4
I
Based on the above checks answer Question 3.4 and enter result in Figure 1.0. h 3.4 Are the FE loads and boundary conditions applied correctly? I I Result
Comments
2-19
,,
,..,
\
,,
U/
4.0 4.1
FINITE General
ELEMENT Solution
RESULTS Checks
CHECKS
Perform these checks to expose any gross errors. Most programs output values of gross parameters associated with the solution process, These parameters typically include summed applied loads and reactions, total mass, position of centre of gravity, etc.
Result
Comments
4.1.1
Are all error and warning messages issued by the software reviewed and understood?
4.1.2
Is the magnitude of mass of the finite element model approximately as expected?
3-4,1
4.1.3
Is the location of centre of gravity of the model, as calculated by the program, reasonable? Are the applied forces in equilibrium with the applied reactions?
3-4.1
4.1.4
3-4.1
Result Based on the above checks answer Question 4.1 and enter result in Figure 1.0. 4.1 Are the general solution parameters acceptable? Comments
2-21
... .-,,,
.l~
4.2
Post Processing
Methods used by the program to
Perform these checks to ensure that the methods, and their limitations, post-process the results are understood.
IE!!lResul
3-4.2 3-4.2
Comments
4.2.1
Are the methods for reducing analysis results described (eg. calculation of safety factors and other parameters calculated by manipulating raw output)? Are the methods for correcting FE results described (@g, correction factors, smoothing factors)?
4.2.2
Based on the above checks answer Question 4.2 and enter result in Figure 1.0. 4.2 Is the methodology used for post processing the results satisfactory? I
Result
Comments
2-22
-L,,
b .-
4.3
Displacement
Results results are consistent with expectations.
Perform these checks to ensure that the displacement
Check
Result
Comments
4.3.1
Are the displacement discussed?
results described and
4.3.2
Are plots of the deformed structure (or mode shape) presented?
3-4.3
4.3.3
Are the directions of displacements consistent with the geometry, loading and boundary conditions? Do the magnitudes sense? of displacements make
3-4.3
4.3.4
3-4.3
4.3.5
Is the deformed shape (or mode shape) smooth and continuous in area of interest?
3-4.3
4.3.6
Are unintentional slits or cuts (indicating elements not connected where they should be) absent?
3-4.3
Based on the above checks answer Question 4.3 and enter result in Figure 1.0. 4.3 Are displacement results consistent with expectations?
[ I
Result -
Comments
2-23
,. . .,,,,
4.4
Stress
Results
Perform these checks to ensure that the stress results are consistent with expectations. Refer To Guideline Section 3-4.4
Check
Result
Comments
Are the stress results described and discussed? Are stress contour plots presented? In the stress plots are the stress parameters or components defined (eg. crX,cfY,TXY, We.)? Is the method of smoothing stress results, or averaging stress results described (eg. element stresses vs nodal average stresses)? Are the units of stress parameters consistent? Are the magnitudes with intuition? of stresses consistent
4.4.2
3-4.4
4.4.3
3-4,4
4.4.4
3-4.4
4.4.5
3-4.4
4.4.6
In cases where there are adjacent plate elements with different thicknesses does the method for averaging stresses account for the differences? Are the stress contours smooth and continuous, particularly in region of primary interest ? Are the stress contours at boundaries consistent with the boundary conditions applied (eg, stress contours perpendicular to boundary if symmetry be)? Are stresses local to the applied loads reasonable? Are there areas in which stresses are above yield (which would invalidate linear elastic analvsis)?
3-4.4
4.4.7
3-4.4
4.4.8
3-4.4
4.4.9
3-4.4
4.4.10
3-4,4
Based on the above checks answer Question 4.4 and enter result in Figure 1.0. 4.4 Are stress results consistent with expectations?
[ I
Result
2-24
f,,< .,.,_,. -
4.5
Other
Results
Perform these checks to ensure that other types of results from the FEA are consistent with expectations. Refer To Guideline Section 3-4,5
Check
Result
Comments
4.5.1
Are the frequencies units?
expressed in correct
4,5.2
Are the magnitudes of natural frequencies consistent with the type of structure and mode number? Are the mode shapes smooth?
3-4.5
4.5.3
3-4.5
Based on the above checks answer Question 4.5 and enter result in Figure 1.0. h 4.5 Are dynamics results consistent with expectations?
I I
Result
Comments
2-25
. ..
5.0 5.1
CONCLUSIONS FEA Results
CHECKS and Acceptance Criteria with
Perform these checks to ensure that the results are in a form suitable for comparison specified acceptance criteria,
Result
Comments
5.1.1
Are the results summarised in a manner that allows comparisons with acceptance criteria, or alternative solutions or data? Are satisfactory explanations provided where the results do not meet acceptance criteria, or where they differ significantly from other comparable solutions or data?
5.1.2
35.1
Based on the above checks answer Question 5.1 and enter result in Figure 1.0. i 5.1 Are the results presented in sufficient detail to allow comparison with acceptance criteria?
Result
Comments
2-26
5.2
Load Assessment Perform these checks and evaluations accuracy, are understood. to ensure that the loads applied in the FEA, and their
Result
Comments
5.2.1
Has an assessment been made of the accuracy or degree of conservatism of the loads used in the FE model with respect to the following aspects : a) types of loads / load cases that were included and excluded
b) basis or theory used to derive loads (eg. linear strip theory for sea motion loads, base acceleration vs DRS for shock, drag coefficients for wind loads, etc.) c) d) e) f) g) magnitudes of loads
loading directions included / excluded load combinations load factors boundary conditions
Based on the above checks answer Question 5.2 and enter result in Figure 1.0. 5.2 Are the accuracy and conservatism, understood? Comments or otherwise, of the applied loading modelled
Result
2-27
-.,
5.3
Strength
/ Resistance
Assessment to ensure that an adequate assessment of the
Perform these checks and evaluations capability of the structure
has been made.
Check
Result
Comments
5.3.1
Has an assessment been made of the accuracy or degree of conservatism of the strength or resistance of the modelled structure with respect to the following aspects : a)
failure theory, failure criteria, allowable stresses, safety factors, etc section properties material properties allowances for imperfection, manufacturing tolerances allowances for corrosion misalignment,
b) c) d)
e)
Based on the above checks answer Question 5.3 and enter result in Figure 1.0. 5.3 Has an adequate Comments assessment been made of the capability of the structure?
l==
2-28
,....
5.4
Accuracy
Assessment to ensure that an attempt has been made to assess the
The checks listed below are intended accuracy of the FEA.
Result
Comments
5.4.1
Has an assessment been made of the scale of FE model and its level of detail and complexity? Have the types of behaviour modelled and not modelled (eg. membrane only instead of membrane plus bending) been assessed? Has the influence of mesh refinement accuracy been considered? on
5.4.2
3-5,4
5.4.3
3-5.4
5.4,4
Has a comparison with other results (eg. other solutions, experiment, etc. ) been made?
3-5.4
5.4.5
Based on the above has an overall assessment of the accuracy of the relevant results been made?
3-5.4
Based on the above checks answer h 5.4
Question 5.4 and enter result in Fiaure 1.0.
I I
m I
I
1
Has an adequate assessment of the accuracy of the analysis been made?
Comments
2-29
.. ...
%, , /, j,
?.J
5.5
Overall
Assessment
The checks listed below are to ensure that the overall conclusions and recommendations resulting from the FEA have been presented and are generally satisfactory.
Result
Comments
5.5.1
Are conclusions from the FEA provided, and are they consistent with the material presented? If appropriate has a way ahead or potential solutions been presented?
5.5.2
3-5,5
5.5.3
Based on consideration of all previous checks is the overall assessment that the FEA is acceptable?
3-5.5
Based on the above checks answer Question 5.5 and enter result in Egure 7.0. 5.5 Is the finite element analysis assessed generally satisfactory? Comments
Result
2-30
GUIDELINES FOR ASSESSING
PART 3 FINITE ELEMENT MODELS AND RESULTS
The guidelines recommended below are structured to match the Assessment Methodology described in Part 2, Therefore, the guidelines are grouped under the same five sections: 1. 2. 3. 4. 5. 1.0 Preliminary Checks Engineering Model Checks Finite Element Model Checks Finite Element Results Checks Conclusions Checks CHECKS
PRELIMINARY
This section describes the checks that need to be undertaken to ensure that the finite element analysis (FEA) satisfies certain basic requirements. The first requirement before evaluating an FEA is to ensure that there is sufficient documentation provided with the analysis. This step should ensure the analysis addresses the objectives, scope, It is necessary to establish that the tools and requirements of the work specification. the analyst uses in the FEA are adequate and appropriate to the analysis; this applies particularly to the software used. Finally, the analyst should be appropriately trained and should have sufficient experience. 1.1 Documentation Requirements is an essential part of any FEA. The documentation submitted party
Proper documentation
should be sufficient to allow a through evaluation of the FEA. The complete documentation package, which can be defined as that required by an independent to reproduce the analysis, should be available and submitted if required by the evaluator. The complete documentation would typically include:
q q q q
9
project data scope and objectives of the analysis list of reference documentation drawings and sketches of the subject structure description of the engineering model rationale for using FEA software and hardware used in the analysis description of the finite element model assumptions used in the analysis description of the results assessment of accuracy of the results conclusions and recommendations or textual form depending
q q q q q q q
The input and output data should be presented in graphical on what is the most convenient for evaluation purposes.
3-1
The documentation requirements listed in Part 2, Section 1- Para 1.1, are the minimum required. In general, any additional information considered necessary for a complete evaluation should also be provided. Plots should be properly annotated to show the location of the subject structure in the ship (eg,, frame numbers, deck numbers etc.), axes to orient the model, location of equipment supported by the structure, and the position of major structural features that define boundaries (eg. bulkheads), All symbols used in the plots should be defined either on the plots or in the body of the report. 1.2 Job Specification Requirements according
The purpose of this check is to ensure that the analysis has been undertaken to the requirements of the job specification. This can be done ,only if the
documentation provided addresses every requirement of the job specification. It is not possible to list all such requirements, but at least the following items should be addressed:
q q q
definition of the problem scope and objectives of the analysis all relevant documentation such as drawings, define the subject structure and loading any previous analyses, service experience subject structure acceptance criteria (eg. allowable
sketches
and reports to completely data related to the
and experimental
stress in an analysis in support of a design)
It is expected that the analyst has carefully read the job specifications and followed it as closely as possible. Deviations from the specifications, if any, should be identified and justified. All reference documents should be identified. If the job specification does not specifically call for a FEA, then the analyst should
explain the rationale for using FEA in preference to another method of structural h is also expected that the analyst is aware analysis, or in preference to experiments. of any previous related studies and their outcome. The selection of FEA as the preferred method of structural analysis will depend on many Features of the problem that should be discussed features of the engineering problem, include, but are not limited to, the following:
q q
purpose of analysis; complexity of the structural redundancy of structural assessment of expected accuracy of known input suitability, or otherwise,
form;
9 .
q q
system; accuracy; variables such as loads, material of hand calculation methods.
properties,
etc.; and
3-2
1.3
Finite Element
Software
Requirements software systems on the market, Most are intended for
There are many finite element
general purpose FEAs, while others are specialist in nature. Ship structure FEA is, to a certain extent, specialized in nature and therefore not all FEA software will perform adequately. It is essential to establish that the software chosen for the job has the In addition it is necessary to ensure that the software has been required capabilities. verified and validated, Commercial maintaining finite element analysis systems are large and complex. Developing and such systems require systematic methods to be applied to the design and
development of the code, the testing, the verification and validation of the code, and the configuration management of the software system. Reputable software vendors rely on quality systems to ensure that the relevant processes that comprise the development and maintenance of the software of FEA software should include an assessment There are several ways in which validating FEA software include:
q
are properly controlled. The evaluation of the vendors quality system. can be validated. The methods for
FEA software
b
q
independent analysis experimental results service experience
Many finite element software vendors publish verification examples, Generally the verification examples are based on problems with closed form solutions. The analytical results are compared with those obtained by exercising the finite element code, While a comprehensive set of satisfactory code it does not constitute proof. verification examples is convincing evidence of good Verification examples based on problems based on
closed-form solutions are necessarily simple and the finite elements models are generally not too demanding on the software. It is necessary, therefore, to employ additional methods to validate the software. An additional validation method is to use benchmark problems that, while simple, are more representative of typical structure, In contrast to the type of verification example mentioned above, benchmark problems can be designed to use combinations of element types, element shapes that vary from the ideal, complex boundary conditions, multiple load cases etc. to test the software, These problems more closely relate to the way in which the software will be used in practice. Closed form solutions are generally not available for benchmark problems. However, results from other well-established FEA software could be regarded as an example of an independent analysis. If results from several other FEA software systems are consistent, or where any differences can be rationalized, then these results can be regarded as benchmarks. Any significant differences between benchmark results and those obtained from the candidate FEA software system would be an indication of unsatisfactory performance.
3-3
Depending on the size of the organization and the volume of FEA work, it may be useful to maintain a register of FEA software validated based on satisfactory performance using the methods outlined above. Alternatively this function could be performed by a body representative of the industry such as a professional society. In the absence of such an arrangement at present, benchmark problems typical of ship structures have been formulated and the results documented in Part 4 of this report. These benchmark problems could be used to evaluate candidate FEA software. If the contractor has documented evidence (based on previous applications of the software to ship structural analysis problems) that the software is capable of performing the required analysis, this requirement may be waived at the discretion of the evaluator. Successful necessary, performance of the candidate FEA software on the benchmark problems but not sufficient, condition for approving the software. The software is a
should also satisfy requirements outlined in the opening paragraphs of this section particularly in regard to requirements for the vendors quality system. 1.4 Reasons for Using A Particular FEA Software Package
It is recognized that the contractor will prefer to use FEA software packages that are readily available and that the analyst has experience with, However, the contractor should make an assessment of the suitability of the selected FEA software for the analysis under consideration. The items that should be discussed include the following:
q
availability
of required element
types
. . . . s 1.5
availability of required material types availability of required load types capability of the software to perform required analysis preprocessing and postprocessing capabilities support from vendors Competence and checking the analysis must meet minimum training and The following aspects of personnel background will need
Personnel
The personnel performing experience requirements. assessment: .

q q
formal academic or professional qualifications engineering expertise in design and analysis of ship structures relevant experience in the modelling and analysis of design problems using the finite element method familiarity with, and appreciation of, the limitations of the particular software employed
Personnel are grouped in two categories: analyst and checker, The analyst is a person who undertakes the FEA, The checker performs independent checks of the analysts work, and certifies the quality of the work.
3-4
~,. ... . .,
The contractor should satisfy the client that the analyst and checker meet the competence requirements, and assure the client that sufficient resources are applied to allow the FEA to be undertaken proficiently. 1,5.1 Academic and Professional Qualifications
The analyst and the checker should be qualified to first degree level in engineering or naval architecture, and have taken at least one full course in structural FEA, Professional Engineer (or equivalent) status is essential for the checker and desirable for the analyst, 1.5.2 Training and Experience
The analyst and checker should have received training in the application of the finite element method, Either of the following is acceptable, in principle, as training:
q
Training provided by various courses offered by educational establishments and software vendors. These courses are only acceptable if they are application oriented. In-house formal or informal training provided by a supervisor capable of satisfying the requirements of a checker, The content of the training should be at least equivalent to a one week application course/s should be documented. oriented training program. The training
The analyst
or checker
must be familiar with the design requirements,
codes of
practice, analysis and design standards have, and the analyst should preferably size and complexity The checker
relating to ship structures. The checker must have, experience with analyses of comparable
as the analysis under assessment, analyst with substantial experience in the
should be an experienced
application of the finite element method, This experience should include working as an analyst on finite element analyses that are comparable in complexity to the analysis the checker will be verifying. The documentation should include a brief outline of previous experiences . The experience requirements for analysts recommended by NAFEMS (NAFEMS, 1990) is summarized in Table 3-1,1, The experience required of the analyst depends on the criticality of the analysis. The criticality failure of the structure being analyzed. category depends on the consequences of
3-5
.-..,. .
Analysis
Category
Engineering Experience Design & Analysis Experience
FE Modelling
and
Problem Solving FE Experience After Formal Training for Each Analysis Type 6 months Relevant Jobs Performed
1. Vital -endanger human life, or property or the environment on a scale of a public disaster 2. Important -Category 1 problem however analysis is not an exclusive part of the integrity demonstration 3. Advisory -All analysis other than the ones covered in Categories 1 and 2
5 years
2 x Category 1 under supervision or 5 x Category 2 properly assessed
2 years
2 months
1 x Category 1 or 2 under supervision or 3 x Category 3 properly assessed
1 year
1 month
Prescribed Benchmarks
For example,
see Part 3 of this report for benchmark Minimum Recommended Experience
problems from NAFEMS, 1990)
TABLE 3-1.1
Levels (adapted
3-6
......
L....
.,
2.0
ENGINEERING
MODEL
CHECKS
The checks recommended in this section are generic in nature, and form part of any of the engineering analysis. The engineering model is a simplified representation physical problem and hence it is crucial that this modelling process is undertaken correctly since the finite element analysis (FEA) cannot improve on a poor engineering model. The aspects covered in this section include type of analysis, problem geometry, material and physical properties, loads, and boundary conditions. The discussion here is restricted to an understanding of the physical problem, Translating these aspects into a finite element Section 2.1 Analysis 3. Type and Assumptions model, in a format recognized by the software program, is covered in
An engineering model is a simplification and idealization of an actual physical structure or component. The contractor should describe the physical problem, and should include, as a minimum, discussion of the following topics: . .
q q
general description purpose of analysis (eg., design, failure investigation, etc.) whether the problem is static or dynamic appropriateness of linear elastic analysis (nonlinear analysis is not addressed document) assumptions and approximations design criteria if appropriate
in this
q q
that have to be made and their likely implications
The underlying assumptions and decisions made in the formulation of the finite element (FE) model should also be described. This description should include the rationale for:
q
.
q q
including and excluding parts of the structure taking advantage of symmetry, antisymmetry, or axisymmetry identification of dominant structural action whether the structure can be modelled with line elements, area elements, elements or a combination of different element types
or volume
Ship structures are usually complex in nature, and can only be analyzed after idealization of the structure, Several simplifying assumptions are made in the idealization process, In order to do this successfully, it is necessary to have a reasonable qualitative understanding of the expected response. This will allow reduction of the complex response of the actual structure to its essentials. The elements that need to be considered in this idealization process are the character loading, the primary loading paths, and the parts of the structure that participate response,
of in the
The loading will be static or dynamic. Many dynamic loads can be treated quasistatically, Where this is not possible, it will be necessary to consider the frequency range over which there is significant energy in the forcing function. This will determine the number of modes to be extracted.
3-7
..,<..~$
Consideration of the likely load paths will help establish the extent of the structure should be modelled, and what boundary conditions might be appropriate.
that
Most real structures are discontinuous and irregular at a local level, For example, it is likely that there will be brackets attached to the structure, openings, access holes, etc. The explicit modelling of these features is not practicable, and not necessary if global response is of interest. All structures are three-dimensional. Depending on the configuration it is often possible
to reduce the number of dimensions 2.2 Geometry Assumptions
to be considered.
One of the first questions to arise during the planning phase of a FEA is how much of the structure needs to be modelled to yield answers of the required accuracy. This is best approached by considering what the influence on the results of interest is of extending or reducing the extent of the model. If the influence is negligible then the extent of the model can be established in advance. However, performing such an exercise on complex structures through intuition alone is difficult. It is recommended that in complex structures the main structural actions should be identified. Once the main structural actions are identified, it is possible to apply simplified structural models to guide the analyst in deciding the extent of the structure to be modelled; Figure 3-2.1 illustrates the concept with simple examples. The following general principles should be borne in mind when using this approach:
q
Drastic changes in stiffness are potential regions to end the model. Figure 3-2.2 presents an example in which the left-hand side of a beam is supported by stiff structure. The bending stiffness of beams is proportional to l/L3 where I and L are the second moment of area and the span respectively. In this example a difference in stiffness of, say, two orders of magnitude would be sufficient to justify the modelling approach shown in the figure. This general approach can be adapted for other more complex structures. Identification of load paths is a good indicator of which parts of the structure are best to model,
The actual extent
of the finite element
model depends on a tradeoff
between
the
resources available for the analysis and the general requirement portions of the structure be model led. The contractor statement should describe and justify the extent of:
that all significant
of the model.
The justification
should include a discussion
3-8
-,
7
3- SPAN BEAM: SPAN - L; W = 1
IF MODELLEDAS 2. SPAN BEAM
(.+
-6
M = 0.063L
w
M
=0
IF MODELLEDAS 1. SPAN BEAM
,. PLATEWITH HOLE
a
STRESSESESSENTIAUY UNIFORM
-4x
IJNE$OF SYMMETRY
DIAMETER d
FIGURE 3-2.1
Examples
of Simple Models that can Indicate 3-9
Extent of Structure
to be Modelled
THIS PART OF STRUCTURE MUCH STIFFER THAN THIS PART
BENDING STIPFNESS / DECK [:::%%!:N
ill
LOADING
-/
w-
SHORT SPAN OEEp GIRDER
I
BULKHEAO
CAN BE MODELLED AS
FIGURE 3-2.2
q q q
Large Changes
in Stiffness
to Indicate
Extent of Model
.
q
. .
all significant structural action captured by model. requirement to accurately predict stresses and/or deflections. region of structure of patlicular interest, whether St. Venants Principle is satisfied obvious changes in structural stiffness that suggest a model boundary very local application of the load to a large uniform structure for large models, can top-down analysis be used?
If the FEA is concerned primarily with local effects then the concepts underlying St. Venants Principle can be helpful in establishing the extent of model. Essentially this principle states that the replacement of a load (which could be caused by a restraint) by a different, but statically equivalent, load causes changes in stress distribution only in regions close to the change. Figure 3-2.3 illustrates the principle. 2.3 Material Properties materials used in the construction of ships are metallic. Other
The most common
materials also used include GRP and wood. The scope of these guidelines is confined to isotropic materials working in the elastic range. However, certain important considerations in modelling material properties of composite materials are discussed in the paragraphs below.
3-1o
,.
<-...
While Poissons ratio for steel is not very sensitive to increases in temperature, Youngs Modulus does reduce significantly when the temperature starts to get above a few hundred degrees Centigrade, Nuclear air blast explosions can cause thermal effects of sufficient magnitude to influence the value of Youngs Modulus. High strain rates can increase the value of the yield and ultimate stresses of the material. However, these strain rates have to be very high to have a significant effect, Examples where structures may be subject to high strain rates include structural response to underwater explosions and nuclear air blast. As a general guide, the effects of strain rate should be considered for strain rates over 0,1 S-l i
DISTRIBUTED SUPIWRT
POINT
SUPPORTS 1 1
----Illustration
&R;E;s&l;;:mEl--
FIGURE 3-2.3
of St. Venants
Principle
2,3.1
Composite
Materials
Modelling the behaviour of composite materials is more complex than modelling isotropic materials such as steel. Composite materials are anisotropic and cannot always be regarded as a continuum, In cases where global response is of interest, it may be reasonable to model composite materials using an anisotropic continuum model. More local analysis requires explicit modelling of the material. Most general purpose FEA software systems include the capability to compute the elastic properties of composite materials. This is done by defining the individual layers that comprise the composite, Alternatively, it is often possible to input the constitutive matrices that define the relationship between generalized forces and moments to generalized strains and curvatures, The failure modes of composite materials are also more complex than those that typically apply to isotropic materials. To check the adequacy of a structure made from composite materials, it is necessary to define the failure criteria that must be applied. Whereas with isotropic materials a single failure criterion (e.g. yield stress) is typically applied, with composite materials failure criteria are generally different for different directions strains, and can be applied to strains, stresses and combinations of stresses and
3-11
There are other modelling issues that are particular to composite materials. Depending on the design of the composite, it may not be possible to apply symmetry conditions even when the loading and the overall geometry are symmetrical about one or more axes, 2.4 Stiffness and Mass Properties in form and the only physical property required is cross on the other hand, are considerably more complex. needed to define beam elements are discussed in the
Truss elements
are the simplest Beam sections, sectional area. The various sectional properties following paragraphs.
The basic sectional properties required to define beam elements are cross sectional area, shear areas in two orthogonal directions normal to the longitudinal axis of the element, torsional constant, and the second moments of area about two orthogonal axes, The axes are usually chosen to coincide with any axes of symmetry that may exist. While this definition of beam properties is complete for the vast majority of cases, there are circumstances in which additional factors need to be considered. The torsional stiffness is based on the torsional constant alone and therefore no
account is taken of warping effects. Warping is most relevant for open sections. The error introduced by ignoring warping is, fotiunately, usually not serious because of the circumstances in which open sections are generally used in structures. However, in situations where the main structural force acting on an open-sectioned beam is torsion this shortcoming should be considered in calculating rotations and torsional stresses. Structures modelled using standard beam elements in most general purpose FEA software would yield incorrect results. Some FEA software does offer beam elements that account for warping effects. Shear flexibility is important for deep short beams. Ignoring shear effects configuration would result in an overestimate of flexural stiffness. The input data required for plate and shell members computer programs can accommodate nonuniform input different thicknesses at each node. 2.4.1 Mass for Dynamic Problems is thickness. thickness for this
Most finite element
and have the facility to
The subject of mass modelling cannot be treated without some preliminary discussion. The discussion concentrates on two main issues. The first matter is the necessity for reducing most dynamic problems to a manageable size. The second concerns two alternative methods for mathematically representing mass. Each is treated in turn. The main difference between static analyses and dynamics analyses is the far greater
computational effort required for the latter compared with the former, Therefore, it is usually not practicable to treat dynamic problems in the same way as static problems except in the most trivial cases. It is usually necessary to reduce the size of the problem by reducing the number of dynamic degrees of freedom (dof), This may be done explicitly or implicitly depending on the algorithm used for extracting eigenvalues
3-12
Certain techniques, such as Subspace Iteration, implicitly reduce the and eigenmodes. size of the problem. The degree of reduction depends on the number of modes that need to be extracted. The reduction process can also be accomplished more directly by a procedure known as condensation and perhaps the best known such technique is Guyan reduction. While the condensation process is generally detrimental to accuracy, the loss of accuracy need not be significant if the appropriate guidelines are followed. There are two alternative methods for mathematically modelling mass. The simpler of
the two methods is the lumped mass method in which concentrated mass is located at nodes, The value of the mass represents the mass of the surrounding structure and equipment. This approach yields mass matrices that are diagonal. Rotational inertias may also be modelled in this fashion, or can be condensed out, Rotational inertias are often ignored when this method is used, The alternative approach is called the consistent mass method. This is a theoretically rigorous method that results in a mass matrix with off-diagonal terms. The presence of these off-diagonal terms in the mass matrix is responsible for making dynamic analysis using consistent mass matrices more computationally demanding than when using lumped mass matrices. For large models there does not appear to be much difference between the two methods in terms of the accuracy attained, at least for lower frequencies. Whatever the technique may be for calculating natural frequencies mass distribution needs to be accurately modelled. Natural frequencies 1. 2. to compare and modes are calculated and modes, the
for one of the following
reasons: of
natural frequencies
and modes of a structure of structural response.
with the frequency/ies
some source of vibration as the first stage in the calculation
In either case it is necessary to anticipate the results to some extent. In the first case the natural frequencies calculated must bracket the frequency of the vibration source. In the second case the spectrum of the forcing function, for example harmonic forces from the propellers or impulse loads from underwater shock, will suggest the range of natural frequencies of the structure that need to be calculated.
The higher the vibration mode, the more detailed the mass distribution needs to be. The general principle is illustrated in Figure 3-2.4. In the actual structure the mass is distributed over the length. Hence, a reasonable number of lumped masses are required to represent the distributed mass. For higher modes a more detailed representation of mass is required because the mode shape is more complex. In the example shown in the figure essentially a single mass is being used to represent the dynamics of one lobe of the third vibration mode. This is in contrast to the five masses used to represent the dynamics of the single lobe in the first mode. 2.4,2 The Influence of Surrounding Fluid
Certain problems in ship structures require that the interaction between the structure and the fluid be considered, The comments made here are limited to cases in which
3-13
.,
fluid displacements structures adjacent
are small. to fluid.
The most common
example
is the vibration
of plated
For vibrations of plated structure adjacent to fluid, the practice is to account for the presence of the fluid by adding masses to the structure to represent the fluid. This mass is usually termed added mass and represents the part of the mass of fluid the There are several sources for data on structure has to accelerate during vibrations. added mass appropriate to plate vibrations (see ISSC, 1991- Report 11.2 for typical sources),
I-IuII
can be treated methods for computing
similarly.
Chalmers
(1 993)
provides guidance
on
approximate
added mass for the hull girder.
The use of added masses to account for fluid-structure effects is generally quite approximate. More rigorous methods require the finite element modelling of the surrounding fluid. Many general purpose FEA systems include fluid elements that allow certain types of acoustics, sloshing and fluid-structure analysis problems to be solved. This is a specialist area, For guidance the reader is referred to finite element texts and the user manuals of the FEA system to be used in the analysis.
BEAM VIBRATIONS
q
1ST MODE
MASSES
ACCEPTABLE
2ND MODE
MARGINAL
3RD MODE
UNACCEPTABLE
FIGURE 3-2.4
Mass Distribution
Required for Accurate
Determination
of Natural
Frequencies
3-14
2.5
Dynamic
Degrees
of Freedom range of interest is decided upon, the mode shape for the highest
Once the frequency
frequency in this range needs to be estimated. This will indicate the number of dynamic Predicting a mode shape in advance is usually dofs required to yield accurate results. very difficult unless the structure is relatively simple. Therefore, it may be necessary to follow an iterative process in which the mass distribution is refined at each iteration. Certain algorithms require any problem size reduction to be undertaken by the analyst. In this case the analyst selects the number of dynamic dofs to be used in the analysis. The selection of the dynamic dofs to be used in the dynamic analysis requires considerable skill except for the simplest structures. The selection of dynamic dofs can be automated. The principle underlying the Guyan reduction process provides a guide on how this should be done, if done manually. The most important dynamic dofs are those that have the largest mass-to-stiffness ratio. This is because such masses are responsible for most of the vibration energy at lower modes. The concept underlying the selection of dynamic dofs is shown in Figure 3-2.5. Viewing a plot of the mode shapes will allow an assessment to be made of the reasonableness of the selection of dynamic dofs,
BEAMVIBRATIONS-
LUMPED MASSES SMALLMASS INCLUDEWITH ADJACENTMASSES
IARGE MASS INCLUDE
MODERATEMASS RtGIDSTRUCTURE IGNORE
IARGE MASS FLEXIBLESTRUCTURE INCLUDE
FIGURE,3-2.5
Selection
of Dynamic
dofs
For most structural dynamics problems translational masses are sufficient to define the problem. However, when components and equipment with large dimensions are being modelled it is prudent to model their rotational inertia, If a single mass element is being used to model the component then three rotational inertias should be input in addition to translational mass data, Alternatively, several masses can be input that approximately simulates the mass distribution, The procedures are summarized in
3-15
Figure 3-2.6. A summary of guidelines to be followed in selected in dynamic dofs is given below:
1, The number of dynamic dofs should be at least three times the highest mode required. For example, if thirty modes are required at least ninety dynamic degrees of freedom should be specified, 2. Dynamic dofs should be located in regions where the highest modal deflections are 3. 4. 5. 6. 7. anticipated. Dynamic dofs should be located where the highest mass-to-stiffness ratios occur on the structure. If a dynamic response computation is to be eventually performed dynamic dofs should be located at points where forces are to be applied, For slender structures, such as masts, only translation dynamic dofs need to be selected. For stiffened plate structures only dynamic dofs at right angles to the plane of the mass does not
structure need be selected. Enough dynamic doffs should be retained such that the modelled differ from the actual mass by more than 10YO,
MODELLED AS MOMENTS OF INERTIA SHOULD
BEINCLUDED lx,If+ 12
/
FIGURE 3-2.6
Modelling
Rotational
Inertia
2.6
Loads and Boundary
Conditions should be described. level of the load. 1984) typically applied in ship structural The description should
All loads that need to be considered include a brief discussion
of the accuracy
Loads (compiled by Giannotti & Associates, analyses include the following:
3-16
\, .... . .
1.
Hull Girder Loads consist of wave induced and still water loads on the hull girder. This load should be considered for longitudinal structure in the main hull, and for interaction of a long continuous deckhouse (superstructure). Hydrostatic Loads are pressure loads due to fluids. The pressure could be either internal or external, Examples of hydrostatic loads are external pressure of the bottom and sides of shell plating, and internal pressure in tanks and on water tight bulkheads, Hydrodynamic Loads consist of liquid sloshing in tanks, shipping of green water on the weather deck and impacting on the house front, and wave slap on all exposed structure and equipment above the waterline, etc. Live Loads consist of uniform deck loading, concentrated loads such as forklift aircraft landing and parking loads, support reactions from stanchions and equipment, cargo container reactions, etc. Dead Loads consist of the weight of the structure.
2,
3.
4.
or
5. 6.
7.
Ship Motion loads consist of inertial forces that act on the entire ship and are important design loads for masts and topside foundations, such as topside cargo attachments. The effect of ship motion loads on the hull girder is to produce vertical and horizontal bending moments and torsion, A lengthy analysis is required to determine these values for a particular ship and service characteristics. Shock Loads consist of displacements, velocities and accelerations in all three directions, This load is important for naval ships in the design of vital equipment and their foundations, and ship structure in the vicinity of these foundations. Missile and Gun Blast Loads consist of a transient pressure and thermal load for all structure within the blast impingement area, usually a static equivalent pressure is used. Nuclear Overpressure consists of transient traveling pressure wave from a nearby nuclear air blast, this is an important consideration in the analysis of deckhouses (superstructures), Vibratory Loads consists of cyclic loading from rotating machinery, especially from propellers, low frequency full girder response from slamming and springing can also be significant, Thermal Loads are caused by heat inputs from: solar radiation exhaust impingement from stack gases combustion engines (important to diesel generator foundations and .
8.
9.
10.
11.
operation of machinery, especially deckhouses and exhaust ducting), condenser 12, Environment foundations
loads consist of wind, snow and ice loads. applied to the model, and the reasons for the The description should include, but not be
A description of the boundary conditions approach adopted, should be described. limited to, a discussion of:
q q q q
model symmetry, antisymmetry and axisymmetry material property changes at the boundary stiffness changes at the boundary assessment of influence on results of assumptions conditions
made concerning
boundary
3-17
3.0
FINITE
ELEMENT
MODEL
CHECKS to ensure that the
The subject of this section is the checks that should be performed physical problem is appropriately provided on various aspects of a element type/s used, the density substructuring and submodelling
translated into the finite element model. Hints are finite element model such as appropriateness of the of finite element mesh used for plated structures, used to optimize the problem size, loads and boundary
conditions, and the solution process. There is also a short subsection on graphical checks using the softwares pre and post processors to scrutinize the finite element model and results. Since access to the software is essential to perform many of these checks, it is the responsibility of the contractor to ensure that these checks are performed. However, documentation, in the form of plots and graphs, should be available for audit. Several examples illustrating finite element modelling practice are presented in Appendix
C. The purpose of these examples is to show the effect of varying certain finite element modelling parameters on the results. The main modelling parameters addressed in this appendix are element type and mesh density. 3.1 Element Types
To some extent all finite element types are specialized and can only simulate a limited number of types of response. An important step in the finite element modelling procedure is choosing the appropriate element/s. The elements best suited to the particular problem should be selected while being aware of the limitations of the element type. A good guide to the suitability of an element type is their performance in other similar situations. Element performance is generally problem dependent, An element or mesh that works well in one situation may not work as well in another situation. An understanding is required of how various elements behave in different situations. The physics of the problem should be understood well enough to make an intelligent choice of element type. As a rough guideline, Cook et al. (1 989) consider elements of intermediate complexity work well for many problems. According to this reference the use of a large number of simple elements or a small number of very complex elements should be avoided. Linear stress field elements are currently the most commonly used. Almost all finite element analysis (FEA) software have families of elements that include elements with linear stress capabilities. For many portions of structures a mesh of linear stress elements can provide a good description of the stress state. In areas of discontinuitie% high thermal gradients, fatigue studies, or nonlinear material problems, where there is an interest of evaluating more than just a linear stress state, linear elements in a relatively fine mesh can give excellent results. Elements with quadratic and higher order stress fields require cubic or higher order displacement functions. These elements have either more nodes per elements and/or more degrees of freedom per node, This make them more expensive in terms of
3-18
L. ,,
computational effort to form the element stiffness matrices, but fewer of them are required than a model using simpler elements to attain the same level of accuracy. Complex structures (eg,, ship deck structure with openings) require relatively fine meshes to model the geometrical discontinuities adequately. According to Kardestuncer (1984) higher order elements are practical only when modelling areas of high stress gradient with a relatively coarse mesh. Even then, the quadratic or higher order fit may over or underestimate the stresses at the free surfaces. The order of the stress function must match the gradient properly, The behaviour of linear stress elements is easy to visualize which is one reason for their popularity. Another limitation higher order elements suffer is the limited availability of companion elements. Lower order element families have a complete range of elements, and therefore it is easier to use these element beams). 3.1.1 Structural types when it is necessary to mix different elements (eg,, plates and
Action
to be Modelled
When a finite element model of a structure is being planned, it is necessary to have a clear concept of the main structural actions. Each element type has limitations and is designed to model a single or limited number of structural actions.
Before modelling a structural problem, it is useful to have a general idea of the anticipated behaviour of the structure. This knowledge serves as a useful guide in several modelling decisions that need to be made in building the model, In an ideal situation the first model will yield adequate results. However, the first model is seldom adequate. Hence, one or more revisions will usually be necessary. In triangulated framed structures, if the members are relatively slender, then the main action is axial with limited bending action. In this case, the use of truss elements would be justified, and the use of beam elements may introduce an unnecessary complication. In certain cases a mixed approach may be appropriate. Consider a lattice mast as shown in Figure 3-3.1. The main legs, which are continuous, should perhaps be modelled using beam elements whereas the bracing members would be better modelled using truss elements. Similarly, deck structure in ships that is subject primarily to in-plane loads, rather than transverse loads, is better modelled using membrane elements rather than plate/shell elements, However, if the analysis of deck structure is local in nature and the loading is transverse, then plate bending elements would be required. In this case transverse shear effects may be significant. Certain element formulations do not account for shear. Some FEA software provide plate bending elements in which the ability to model transverse shear is optional and has to be selected by the analyst. If through elements thickness is prudent. stresses are considered to be important, then the use of solid
3-19
3.2
Mesh
Design
Mesh design, the discretization of a structure into a number of finite elements, is one of the most critical tasks in finite element modelling and often a difficult one. The following parameters need to be considered in designing the layout of elements: mesh density, mesh transitions and the stiffness ratio of adjacent elements. As a general rule, a finer mesh is required in areas of high stress gradient. It is possible, of course, to use a fine mesh over the whole model. This is undesirable on two counts: economy and the greater potential for manipulation errors. Hence, meshes of variable density are usually used, Care is required in transitioning of mesh density. Abrupt transitioning introduces errors of a numerical nature. This subsection provides tips on these aspects of mesh design.
beam elements
truss elements
FIGURE 3-3.1
Typical
Lattice Structure
3,2.1
Mesh
Density
The density of the mesh depends upon the element type used, distribution of applied load and purpose of the analysis. The basic rule is that the mesh is refined most in the regions of steepest stress gradients. Therefore, if such regions can be identified during mesh design, the probability of developing an economical mesh with sufficient refinement is high. In this regard experience plays an important role in striking a balance between economy and adequate mesh density, Analysis of similar structures under similar loading conditions in the past can help in the identification of stress concentrations and regions of rapid changes in stress patterns.
3-20
In cases where experience of a particular configuration is lacking and where it is difficult to anticipate the nature of the stress gradients, an iterative approach is Where stresses show a sharp variation between adjacent elements, the necessary. mesh should be refined and the analysis rerun. If the primary goal of the analysis is to assess deflections, and not stresses, then a comparatively coarse mesh may be used. Mesh density also depends on the type of analysis. A nonlinear or vibration analysis
usually requires a more refined mesh compared to a static stress analysis. Predicting higher frequency modes usually requires a finer mesh than that required for lower frequency modes. Load distribution and load type also have an influence on the mesh density. Nodes at which loads are applied need to be correctly located, and in this situation can drive the mesh design, at least locally, In the case of a uniformly distributed load, such as edge pressures or face pressures, element types that support the particular type of load should be used. Finally, if higher order elements are used with quadratic or cubic stress fields, then a
relatively coarse mesh can be used in the areas of high stress gradients, since the order of the stress function will match the gradient more accurately. For lower order elements with linear or constant stress fields, proper refinement of the mesh is required to obtain accurate 3.2.2 Element Shape results.
Limitations the longest and shortest element
The element aspect ratio is the ratio between dimensions as shown in Figure 3-3.2,
A crude rule of thumb that can be used is to limit the aspect ratio of membrane and bending elements to three for good stress results, and to five for good displacement results. The ideal shape for quadrilateral elements is square and equilateral for triangular elements. Hence, the use of ideally shaped elements is particularly desirable in areas of high stress gradients. In general, higher order elements are less sensitive to departures from the ideal aspect ratio than lower order elements, This observation also applies to solid elements. Since an elements sensitivity to aspect ratio is dependent upon both element formulation and the nature of the problem, general tests and problem dependent may be justified in cases where element performance is not well known, Generally the performance of elements degrades as they become more skewed.
checks
Skewing is defined as the deviation of vertex angles from 90E for quadrilaterally shaped elements, and from 60E for triangularly shaped elements as shown in Figure 3-3,3. For quadrilateral elements, angles greater than 135E and smaller than 45E are not recommended. The limiting range recommended for triangular elements is 45E and 90E. Skewed quadrilateral elements shaped more like parallelograms generally perform better than more irregularly shaped ones.
3-21
0
a
J
b
+3 forsww ~ 5 for displaoamanl
FIGURE 3-3.2
Aspect
Ratio of Plane Elements
When element nodes are not in the same plane, the element is warped as shown in Figure 3-3.3. This is undesirable and the degree to which this impairs the performance of plate elements depends on the element formulation, Hence, the best guidance in regard to limiting levels of warping is contained in the particular FEA programs user manual. high, The use of triangular elements is an option where curvature of the structure is
(a) Skewed Elements
lb) Warped Element
FIGURE 3-3.3
Element
Shape Limitations
3,2,3
Mesh
Transitions
If the mesh is graded, rather than uniform, as is usually the case, the grading should be done in a way that minimizes the difference in size between adjacent elements. Figure 3-3.4 presents several examples of transitions using quadrilateral elements. These examples attempt to keep within the guidelines for element distention discussed in Section 3.2,
3-22
,,
,k --
Another way of viewing good transitioning practice is to minimize large differences in stiffness between adjacent elements. A useful measure of stiffness is the ratio E/Ve, where E and Ve represent the elastic modulus and the element volume respectively. As a working rule, the ratios of E/Ve for adjacent elements should not change by more than a factor of two (Connor and Will, 1969). Sometimes transitions are more easily achieved using triangular elements. Transitions of this type are illustrated in Figure 3-3.5. Most FEA programs will allow two nodes of a quadrilateral element to be defined as a single node in order to collapse the element to a triangular shape.
(a)
(b) a, b) RECTANGULAR
PLATE
(c)
c) CIRCULAR
PLATE
FIGURE 3-3.4
Transitions from Coarse to Fine Meshes
CLOSER APPROXIMATION OF LOAD SINGLILARITV
CLOSER APPROXIMATION OF REALISTIC LOAD
FIGURE 3-3.5
Transitions
Using Triangular Elements
In modern FEA installations
most analysts rely on preprocessors to develop the finite
element mesh. [n general, automatic mesh generators yield adequate meshes. However, in very demanding configurations the mesh generator may produce a poor mesh. In such situations the mesh should be manually improved to meet the guidelines.
3-23
<L
.-
In regular rectangular meshes there are two basic types of transition. One is the change in element density in the direction of the stress gradient, the second is transverse transitioning, which is used between areas with different element size and densities across a transverse plane as shown in Figure 3-3.6.
TRANSITION AREA
(n)
(b)
ELEMENT SIZE CHANGE FIGURE 3-3.6 Mesh Transitions
TRANSVERSE TRANSITIONING
Many rules of thumb for transitioning of elements are based on element strain energy and strain-energy density calculations. The ideal finite element model should have a mesh with constant strain energy in each element. To achieve constant strain energy of elements the volumes must be relatively small in regions of high stress or strain and large in regions of low stress or strain, Transverse transition regions should be used only in areas of low stress gradient deflection, 3,2,4 Stiffness Ratio of Adjacent and never near regions of maximum stress or
Structure
In modelling complex structural assemblies there is a possibility of constructing models where adjacent structural elements have very different stiff nesses. These types of stiffness combinations can cause ill-conditioning of the equilibrium equations which can seriously degrade results, The transitioning guidance given above avoids this problem in models that use two or three-dimensional elements, For truss and frame structures a different approach is required. To prevent large numerical errors in these cases, stiffness ratios of the order of 104 and more between members making up a model should be avoided. This is admittedly a conservative number. More realistic guidance can be obtained by undertaking tests. The problem of stiffness mismatch is most severe in structures where a relatively rigid portion of structure is supported on flexible structure. In such cases the deflections in the rigid portion are due more to rigid-body movement rather than elastic distortion. In these cases it is suggested that the stiff portion be treated explicitly as a rigid body using rigid links, rigid regions, constraints, or combinations of these approaches.
3-24
..
3.2.5
Miscellaneous Improper
Problems between elements of different types can cause errors. Solid
connections
elements types, for example, have only translational nodal degrees of freedom. If solid elements are interconnected with beam or plate/shell type elements, which have rotational degrees of freedom, in addition to translational ones, care must be taken 10 allow for the transfer of moments if that is what is intended, If this is the case then it is best accomplished with linear constraints or multipoint constraints. In case the program does not offer such options, the beam (or plate) can be artificially extended Figure 3-3,7 illustrates the problem and a solution for a through the solid elements. sample problem.
NOMOMENT CONTINUITY
.:
(.
END OF BEAM ELEMENT
MOMENT CONTINUllY PRESERVEI) /
FIGURE 3-3.7
Connecting
Elements
with Different
Nodal Degrees of Freedom
Most flat plate/shell element formulations do not have a shape function for the rotational degree of freedom about a normal to the surface of the element. Hence, inplane rotational stiffness is not modelled, Some programs provide a nominal rotational stiffness to prevent free rotation at the node. Other programs use certain formulations to improve this aspect of performance but at the cost of the presence of spurious modes. The user should be aware of the possible limitations in the program that is being used when modelling situations in which moments are to be transferred into the plane of assemblages solution, is illustrated of flat plate/shell in Figure 3-3.8, elements. The problem, and one possible
3-25
NO ROTATIONAL STIFFNESS
ROTATIONAL STIFFNESS RESTRAINED
RIGIDLINK
FIGURE 3-3.8 3.3 3,3.1
Modelling
in-Plane Rotational
Stiffness
of Membrane
Elements
Substructures Substructuring
and Submodelling
The primary reason for using substructuring is to reduce computational effort in the solution process, However, this saving has to be traded-off against certain other computations that substructuring requires which a normal analysis would not entail. Irons and Ahmed (1 980) identify three circumstances in which substructuring might be attractive: 1. 2. 3, The same substructure is used repeatedly in the structure,
A relatively small portion of a structure may behave nonlinearly, In a major design effort, different teams may be developing different parts of the structure. The use of substructuring would allow substructures of different versions of parts of the structure to be analyzed together. This feature could be very useful during the exploratory and concept design phases of large structures, is becoming of less
Limited computer core capacity as the reason for substructuring concern as the cost of computer memory decreases.
The use of substructuring in the FEA of ships is only likely to be attractive for models involving a substantial portion of the ship. If a general purpose FEA system is used it is essential to have an understanding of the substructuring technique, Even in the case of design-oriented FEA programs it is useful to have an appreciation of the technique. The ease with which substructuring can be undertaken depends on the features
available in the FEA system being used. This section will be confined to a broad description of the steps necessary to undertake successful FEA using substructuring, guidelines in using substructuring techniques, and structural configurations where such techniques might be considered.
3-26
The basic steps in FEA using substructuring 1. Review repeat, 2. 3.
are:
of the global model and identification of portions of the structure that Sketch of the global model indicating substructure boundaries, Design of
mesh in substructures and determination of boundary nodes, Enter input data. Undertake condensation of substructures and develop substructure stiffness and load matrices, Generation of global stiffness matrix which, in general, will require combining the reduced substructure matrices with portions of the structure not modelled as substructures. At this point all the elements of the system equilibrium equations are available. Solve the system equilibrium equations. This run will only yield displacements at substructure boundaries and portions of the model that were modelled in the usual way. The displacements from the global model can be back substituted into the substructure equations, as described below, to yield displacements and stresses within the substructures. This will be repeated for each substructure since, in general, the boundary displacements for identical substructure models will be different,
4.
5.
The following
guidelines
for substructure
analysis are adapted
from Steele (1 989):
1, Substructures can be generated from individual finite elements, from other substructures, or both. 2. Master nodes to be retained must be identified and specified as input when the stiffness matrices for substructures are calculated, Master nodes include boundary nodes and nodes subject to loads, Nodes on substructure boundaries that will be used to connect the substructure to the rest of the global model must be retained as master nodes, Nodes constrained in substructures when substructure stiffness matrices are calculated will be constrained in subsequent stages of the analysis. These constrained nodes cannot be released in later stages. However, master nodes can be restrained during analysis of the global model. For a substructure to be cost-effective it should be used at least three times (i.e., replicated twice).
3. 4.
5.
The following paragraphs contain a description of static condensation, which is a technique fundamental to substructuring. Also discussed is the two-stage analysis technique which has found favour with many analysts. This is followed by a summary of recommendations. 3.3.2 Static Condensation a portion of (dof) the Figure 3-3.9. indeed, is are
In the condensation technique the number of degrees-of-freedom (dofs) in the structure is reduced by condensing out the internal degrees-of-freedom remaining active ones being on the boundary. The process is illustrated in This substructure can be regarded as a special type of finite element, and, sometimes referred to as a superelement. The mathematics of the process relatively simple,
3-27
The equilibrium follows:
equations
of the substructure
with all its dofs intact is partitioned
as
Iuk}=-t}
then be substituted in the upper partition to yield: {fC}
(3.3.1)
in which the subscripts r and c refer to dofs to be retained and condensed out respectively, An expression for i5Ccan be extracted from the lower partition, which
can
( [I(J [km] [Q[kw] ){~r}= [km][k=] {f,}

or in more compact form:
(3.3.2)
[m}= {%}
where [%] and [%-
(3,3,3)
[%1[w%]
Fcl {fr} [%] [Q
{f.} can be solved in the usual way, can be recovered by static Static reduction process, are called If
The equilibrium equations given by Equation (3,3.3) required, displacements internal to the substructure
condensation of Equation (3,3,1) using the Gaussian reduction procedure. condensation amounts to eliminating selected variables using the Gaussian procedure. It is important to note that no approximation is involved in this The condensed out dofs are often called slave dofs and the retained dofs master dofs, 3.3.3 Two-Stage Analysis local mesh refinement is required a two-stage
In cases where
analysis may be justified
The first stage of this technique involves the analysis of a coarsely meshed global model. The local area of particular interest is remeshed using a finer mesh and reanalyses using prescribed displacements at the boundary of the refined model as boundary conditions, The prescribed displacements are taken from the global analysis. The process is illustrated in Figure 3-3.10. The applied loading, i.e., stresses from the global analysis translated into pressure loading for the refined model, can also be used as boundary conditions. Howeverr the use of displacements as boundary conditions is a more common practice since it eliminates the need to provide additional restraints for sufficiently supporting the model.
(see Steele, 1989 for practical aspects of two stage analyses).
3-28
INTERNAL
dofs TO BE
CONDENSED OUT REPEATED SUBSTRUCTURE GIRDERS 6C /
7 BEAMS
.P.
ONLY BOUNDARY dofsTOBERETAINED
FIGURE 3-3.9
Schematic
lllustrationof
The Static Condensation
Process
Design-oriented FEA programs such as MAESTRO, which model the whole or a The displacements from a model substantial part of a ship, suit this technique. developed employing such programs can be used as prescribed boundary conditions a local fine mesh model. In general, there will be several nodes on the boundary
for
of the refined mesh model that
are not modelled in the global model, Therefore, prescribed displacement values are only available for boundary nodes that exist in the global model. The practice is to assume a linear variation in displacement, interpolated from the displacements from the global model, for intermediate nodes. This observation is suggestive of where the appropriate position for the boundary might be, Ideally, boundaries should be placed in areas where gradients in displacement are small, A comparison of unreflected and deflected plots of the global model will yield this information. A finer finite element model is generally more flexible than its coarser equivalent. Hence, there will be a tendency to underpredict the stresses in the refined model when using displacements generated in the global model. R is possible to correct approximately for this tendency using a procedure described by Cook et al, (1 989), The procedure requires the computation of the nodal loads produced by the prescribed The nodal loads for the local area in the global model are boundary displacements. given by:
3-29
EXTRACTRESULTS FROMGLOBAL ANALYSIS 9 / 1= ~ J
H /
DISPLACEMENTSOR F INTERMEDIATE MODES LINEARLY INTERPOLATED

FROM ADJACENT NODES
BP
Ml
/
I
DEVELOPAND ANALYSIS REFINED MODEL
4
\
PRESCRIBED DISPLACEMENTS FROM GLOBALANALYSIS APPLIEDAT ORIGINAL NODES
FIGURE 3-3.10
Two-Stage
Analysis
in which KW 5g, and F~ are the stiffness matrix, displacements, and calculated forces pertaining to the degrees of freedom associated with the nodes on the boundary of the local area. The corresponding expression for the refined model is:
{Fr)=[Kr]~r}
The subscript r refers to the refined model. Note that only the nodes common to both, the local area in the global model and refined model, are included in the above expressions. Once the forces for both cases have been derived, the vector norms for these quantities are calculated. The norm, is a measure of the size of vector, or the size of the nodal loads. There are many types of norms, but for present purposes the following version is recommended:
II~11=
(5I ~,v)%
i-l
where Fi refers to the value of nodal load and n is the number of degrees of freedom the boundary that are common to both the local area of the global model and the
on
3-30
refined model, as follows:
The ratio of the norms for both the cases is calculated
to yield a factor
Factor
-~
r
This factor, which usually exceeds unity, when applied to all stress results from the refined model, approximately corrects for the overstiffness of the global model results. The convenience software 3.4 with which this technique can be applied will depend on the FEA
being used. Conditions appropriate boundary conditions for the model is often
Loads and Boundary The task of selecting
challenging. Generally, the support condition assumed for the degree of freedom concerned is idealized as completely rigid or completely free. In reality the support condition is usually somewhere in between. Several techniques are used to minimize the impact on the analysis of the assumptions made in boundary conditions. The most popular is to develop models large enough such that the area of interest is sufficiently practice to make conservative assumptions bound solutions. remote from the boundary, It is also the so that the results will represent upper
The best guide for determining the extent of structure to model and determining locations for boundaries are natural structural restraints or rigid or stiff supports as: major structural bulkheads, vertical pillars and columns or other structural components such as deep fabricated beams and girders. It is possible to simulate various types of symmetry, antisymmetry
the such
and axisymmetry
by
applying the appropriate boundary conditions. These and other topics related to boundary conditions are discussed in greater detail below, 3.4.1. Minimum Support Conditions
For certain models it is necessary to provide the minimum support for the structure. A good example of this is hull girder modelling in which the structure is, in reality, supported by the pressure distribution on the hull, In FEA modelling a structure with self-equilibrating forces, without any supports, is not admissible. Without proper support the equilibrium equations would be singular and therefore not solvable. Models in a plane have three degrees of freedom, and hence need to have two
translations and a rotation constrained. Care is needed in avoiding the possibility of Models in threerigid body motion. These principles are illustrated in Figure 3-3,11, dimensional space need three translations and three rotations constrained. Examples illustrate minimum support conditions required are provided in Figure 3-3,11.
to
3-31
3,4.2
Boundary
Conditions
for Simulating
Symmetry
Many structures have one or more planes of symmetry, It is possible to take advantage of this in FEA, and model just one portion of the structure. Through various devices it is possible to analyze structures with a plane of symmetry but subject to nonsymmetric loads. Such approaches are used to reduce modelling and computational effort. In engineering applications, the most commonly encountered types of symmetry are: reflective symmetry, rotational symmetry and inversion symmetry as shown in Figure 33.12, In engineering symmetry, problems the characterization of symmetry requires not only geometrical and restraintsi
but also symmetry
with respect to material
properties
When only part of a symmetric structure is modelled, the symmetric or antisymmetric boundary conditions must be applied at artificial boundaries introduced because of symmetry, If the y-z plane is the plane of symmetry, and Ux, Uy, Uz, and Rx, Ry, Rz are assumed as the x, y and z components of displacement and rotation respectively, the following boundary conditions have to be applied to the nodes on the plane of symmetry or antisymmetry: Ux = Ry = Rz = O Rx = Uy = Uz = O - for symmetry - for antisymmetry
In the case of symmetry the points lying in a plane of symmetry can suffer no translation out of the plane and no rotation about the inplane axes. For antisymmetry the complementary set of degrees of freedom are constrained.
The above discussion has been devoted exclusively to static problems, but free vibration problems (eigenvalue problems) can also exploit symmetry. The calculation of all natural frequencies and mode shapes of a symmetric structure would require one modal analysis for each unique combination of symmetric and antisymmetric boundary conditions. symmetry, When only symmetric boundary conditions are applied to the plane of antisymmetric frequencies and mode shapes are not calculated.
The conditions for static problems discussed above apply equally to linear (timeIn addition, if the load is not symmetric or antisymmetric it will be dependent) analysis. necessary to decompose the load into symmetric and antisymmetric components and run the problem twice for each case and combine the results,
3-32
ACCEPTABLE RIGID
NOT
ACCEPTABLE POSSIBLE POINT
BODY MOTION
ABOUT + THIS
EE1 t t
W,u =
W.o free
w tU
2-D problems;
3 independent
conditions
required
A w
U=v=w
oru=O
= o
u plate -
3-D
problems:
6 independent
conditions
required
FIGURE 3-3,11
Minimum
Support
Conditions
for Models 3-33
L
Y (a) Reflective
PLANE OF SYMMETRY
L-
AXIS OF
SYMMETRY
(c) Inversion
(c)
T
3-34
CENTER OF SYMMMETRY
FIGURE 3-3.12
Different
Types of Symmetry
3,4,3
Constraints Constraints are enforced relationships between the dofs of several nodes. There are many situations in which constraints can be useful modelling devices. Various types are discussed below and illustrated using simple examples. The circumstances in which they may be applied, and limitations in their application, are also discussed.
The simplest form of constraint is when certain dofs of different nodes are coupled. Coupling can be used to enforce symmetry and to release forces and moments. A simple example is presented in Figure 3-3.13. During analysis, if the independent node is displaced in the y-direction and/or rotates about the y-axis, the dependent nodes are automatically displaced by the same magnitude conveniently in the same directions. For example, a pin can be
Releases can be introduced
using coupling.
introduced at mid-span in a continuous beam by coupling translational degrees of freedom of two coincident nodes, In certain circumstances coupling can introduce apparent violations of equilibrium. A more powerful and general method for introducing constraints is by using constraint
equations: A constraint equation is a linear equation that relates the displacement or rotational dofs of nodes, These are sometimes referred to as multi-point constraints (MPC). Constraint equations may be used for many purposes such as coupling of nodes by rigid members, rectifying small geometric discrepancies, and coupling adjacent nodes representing locally offset supports and attachments. Rigid regions in structure may be defined using constraint equations, Figure 3-3,14 illustrates the use of constraint equations using the example shown in Figure 3-3,13. In this case the equation ensures that there is no relative movement between Nodes 1 and 2 in the x-direction.
3,4,4
Loads - General Loading in finite element modelling may be applied in a variety of ways, Typical structural loads are forces, pressure load, gravity, body forces and temperatures at nodes and on elements of the model. The load can be applied to: 1. 2. 3. nodes (eg., nodal forces and body forces); element edges or faces (eg., distributed line loads, pressure) the entire model (eg. gravity loads).
applied
Generally the load types and method of its application to the model are specific to a particular FEA software package. However, descriptions of typical load types are provided in the following paragraphs.
3.4,5
Loads - Nodal
Force and Prescribed
Displacement A nodal force
A nodal force is the combination consists of: 1, 2, force magnitude moment
of forces applied to the six nodal dofs.
in X, Y and Z direction;
and elements).
magnitude
about X, Y and Z axes (for structural
3-35
)Y
Node 1 is independent FIGURE 3-3.13 Axis Nodal forces are usually applied in Nodal Coordinate System as shown in Figure 3-3.15. Coupled dof: Nodes 1, 2 and 3 Coupled in the y-Direction and About the y
Applied nodal loads must be compatible with the element type used. For example, a model consisting of only solid elements has no rotational degrees of freedom, Any nodal moment loads would have to be applied in such a case as a force couple with the forces acting at different nodes, Also forced or prescribed nonzero displacement may be input directly to nodes as a load because small changes
case, This displacement should be prescribed with precision, can cause large differences in stress response. 3.4.6 Loads - Nodal Temperature A nodal temperature is a single temperature value illustrated in Figure 3-3.16, A pair of values may surface temperatures. Some programs allow the representing the shell mid-plane temperature and 3.4.7 Loads - Face Pressure
or pair of values applied to a node as represent the shell top and bottom specification of a pair of values a gradient,
A face pressure is a single pressure value applied to selected faces of elements as shown in Figure 3-3.17. The units of pressure value are force per unit area. The pressure is applied to each selected element face across the entire face, and acts in a direction perpendicular to the face. Some FEA programs allow the user to specify pressure at nodal points. A variation of pressure over an element surface can thus be defined. A constant pressure is then a special case corresponding to all element nodes having the same pressure,
3-36
NODE 1- INDEPENDENT
NODE2 - DEPENDENT
MPC: (1)X1-(1)X2 FIGURE3-3.14

Constraint Equation
= 0.0
w t
kY x FIGURE 3-3.15 Definition of Nodal Force 3-37
IL
Y x FIGURE 3-3.16 Definition of Nodal Temperature
FP
FIGURE 3-3.17
Definition
of Face Pressure 3-38
3.4.8
Loads - Edge Loads An edge load is the combination of the forces and moments that can be applied to the edge of an element as shown in Figure 3-3.18. The types of edge loading depend on the type of element, An edge load can be applied to beam elements as: 1. 3. 4. axial force torque bending moment, /
2. shear force
Uniformly distributed loads on beam elements can be handled exactly and no further subdivision of the beam element is required to improve the representation of the load. For membrane elements edge loads can be applied as in-plane forces, and for plate
bending elements both in-plane and out-of-plane bending moments.
forces can be applied along with
3,4,9
Loads - Thermal A beam temperature is the temperature at the centroid of the beams cross section and is applied as temperature, Y axis gradient or Z axis gradient in degrees as shown in Figure 3-3,19, Most programs allow for input of thermal loading directly on elements. Others permit, in addition, specified nodal temperature and temperature-dependent material properties.
FIGURE 3-3.18
Definition
of Edge Pressure
3-39
BEAM
N2
NI
FIGURE 3-3.19
Definition
of Beam Temperature
3,4,10
Gravity
and Acceleration
Inertial loads are generated as a result of the body accelerating. A special case is the self weight of a structure, or body, which is generated by the acceleration due to gravity. Inertial loads are generated 1. 2. 3. translational acceleration angular velocity angular acceleration as a result of one or more of the following:
FEA software systems treat weight data in different ways, It is important therefore, particularly for dynamics problems, to be aware of the way in which the system treats mass, and gravitational forces.
3.5 3,5.1
Solution Static
Options
and Procedures
Analysis the displacements, stresses, strains, and forces in
Static analysis is used to determine
structures due to loads that do not induce significant inertia and damping effects. The loads and the structures response are assumed to vary slowly, if at all, with respect to time, The primary application of FEA in ship structures is in support of design and this usually involves static analyses. These may range from global models encompassing the whole ship, to very detailed local models, Apart from FEA performed in support of design, static analysis is also used in the investigation of certain types of structural failures.
3-40
<: .%. .
3.5.2
Dynamic Dynamic 1,
Analysis analyses in ship structures are usually performed for the following reasons:
2.
To ensure that the natural frequencies of sensitive structures and components do not coincide with those of the hull girder or with the forcing frequencies associated with propellers and other mechanical sources of vibration energy. In preparation for dynamic response computations.
Several quasi-static design procedures have been developed for design against dynamic load conditions, For some of these procedures, for example the Design Response Spectrum Method used for shock analysis, it is often necessary to compute several tens of natural frequencies of the subject structure or component. In complex structures such as masts the natural frequencies and modes can usually only be calculated using FEA. As an alternative to quasi-static procedures, more rigorous dynamic response
calculation may be used. Two methods are available: direct integration of the equations of motion, or the superimposition of modal responses. For nonlinear behaviour, such as that associated with large deflections and/or plasticity, only the former is appropriate. Transient dynamic response analysis is used primarily for computing response to suddently applied loads and/or short duration loads. Examples include forces due to collisions, wave slamming, and shock and blast. In these cases the loading is very uncertain. Various procedures have been developed to compute loads from these types of loading. For example procedures are available to model the shock forces generated as a result of underwater explosions. The procedure models the underwater explosion, the pressure induced on the hull, and finally the transmission of the dynamic forces through the hull structure to the structure or component in question, Many transient dynamic problems involve fluid structure interaction phenomena where the structural response affects the loading on the structure. Sometimes it is possible to treat such phenomena very approximately elements adjacent to the fluid. by adding a certain amount of fluid mass to the
3,5,3
Buckling
Analysis
Depending on the structural element, the estimate of buckling load can be very sensitive to the inevitable presence of discontinuities, imperfections and residual stresses. The application of FEA techniques to solving buckling problem should be approached with caution. The results can be very sensitive to assumptions made in regard to deviations from the ideal, more so than is typical for linear static analysis The usual practice ,in design situations is to adapt classical solutions to the problem.
3-41
4.0
FINITE
ELEMENT
RESULTS
CHECKS analysis (FEA) should always be verified, and
The results obtained
from a finite element
their validity established. To make sure that the results are devoid of any errors in modelling or analysis, it is necessary to perform the checks outlined in this section. These checks ensure that the FEA results are calculated, processed, and presented consistently with the analysis requirements.
4.1
General
Solution
Checks checks can be performed systems, using the graphical such features display features
Many of the following available
with most FEA sollware
Where
are not available,
these checks will have to be performed
by examining
printed results output.
4.1.1
Errors & Warnings Well established finite element software systems to identify poor modelling and analysis practices. generally have several built in checks A warning or an error message is
issued when built in criteria are violated. The correct practice is to resolve any such message is not messages and take the appropriate remedial action, If the warning/error applicable to the analysis, proper justification should be provided. An example could be a warning message for angle between adjacent edges in a quadrilateral shell element. The generally recommended range is between 45 and 1350. If this rule is not followed, valid justification could be that the element in consideration is located well away from the area of interest.
4,1,2
Mass
and Centre
of Gravity of the models
It is good practice to verify the mass of the model and the location
centre of gravity of the model. Several programs provide the mass without the need for a full analysis, If this option is unavailable, the analysis could be run with a 1 G loading (with no other applied loads).
4,1.3
Self-Consistency The results should be checked for self-consistency, For example, displacements at fixed supports should indeed have zero displacements, and any symmetries in the model should be reflected in the stress and deflection results.
4,1.4
Static
Balance
This is a fundamental check. The applied loads should be compared with the reactions. The check should ,include moments where appropriate. This check ensures that the applied loads and reactions are in balance, and ensures that the user specified loading definitions are properly interpreted by the program, When the applied loads and reactions are not in balance this is an indication of a serious error.
3-42
Checking the forces and reactions also ensures that the results are actually for the intended load, In the case of pressure loads, due to possible discrepancies in arriving at nodal forces from pressures, the actual load level could be different from that intended.
4,1.5
Defaults For certain input parameters default All FEA software packages have built-in defaults. values or options are assumed if a value has not been input, or if an option has not been selected. Hence, checks should be performed to ensure that where defaults have been used, they are consistent with the assumptions of the analysis.
4.1.6
Checklist The following is a list of checks to ensure the quality of the FEA, The checklist cover
both prerun and postrun checks. 1. Pre-Run Checks - Graphical: a. b, c, d. e. f, g. h. i. j. k, 1, m. 2. Extremities of model - global dimensions OK Free edges - look for element connectivity Shrunken elements - no missing elements Duplicate nodes Duplicate elements Size of adjacent elements - avoid ill-conditioning Mesh density Mesh transitions Plot material properties by colour Plot physical properties by colour Loads applied to correct elements Direction of loads correct Boundary conditions applied to correct nodes Checks:
Post-Run a, b.
c.
Static balance Comparison i. classical results ii, simple finite element Numerical accuracy i. residuals ii. stiffness ratio Methods
model
4.2
Postprocessing
Methods used for postprocessing of derived quantities from a FEA should be explained. The derived quantities include parameters such as stresses, design margins, factors of safety, etc.
3-43
The need and justification for applying correction factors for FEA results should be explained. The need for applying correction factors may arise due to the necessity to compare FEA results with design codes.
4.3
Displacement
Results
In the design of ship structures the primary result parameter of interest is stress. Most design criteria are expressed as allowable stresses. Although deflection criteria are not as numerous as stress criteria in design codes and standards, they can be just as critical, Stiffness requirements for various components of navigation and combat systems are often quite onerous. Stiffness requirements are often related to dynamic requirements in which the coincidence of equipment operating frequencies and those of the equipment-support structure system is to be avoided. As noted elsewhere, modelling for dynamic analysis is considerably more difficult than modelling for static analysis. This is particularly true for higher modes of vibration. of
In interpreting displacements, it is essential to have an understanding of the accuracy the FEA, how they vary for different response parameters, and the influence on accuracy of modelling decisions made earlier, In general, displacements are more accurately determined by FEA than stress.
The methods used for plotting the displacements of framed structures and certain plated structures in many FEA software packages may understate the actual accuracy. Beams are of-ten plotted as straight lines. In reality the displacement function for beam elements is a cubic polynomial, elements. , The same observation applies to plate bending
In general, displacements in structures composed of beam and truss elements are accurately predicted within the limitations of the engineering model. In terms of the finite element model doubling the number of beam elements in, say, a grillage will not improve the accuracy of the result. The response of two and three-dimensional structures is much more complex and hence, in general, displacement results are sensitive to the fineness of the mesh, Therefore interpreting displacement results in plated and solid models require more care. Gross errors are generally uncovered by the application of intuition and knowledge of previous analyses and physical experiments, More subtle errors are more difficult to uncover. 4.4 Stress Results
As noted earlier, stresses are more difficult to predict accurately than displacements. Limitations in the finite element method are such that stresses are not normally continuous across boundaries between elements. For ease of interpretation of results, most FEA software averages stresses in some fashion before presenting the results. These results are presented attractively as stress contours in colour plots, and the underlying discontinuous nature of the stresses may be obscured as a result of averaging processes, thus engendering a false sense of confidence in the results.
3-44
,.,. ,,
These problems
can be compounded
by misunderstandings
in regard to the type of
stress being plotted. Stress contours provide a good qualitative indication of the adequacy of the density of the mesh, Smoothly changing contours usually indicates that the mesh is suitably fine. Alternatively, stresses in adjacent elements can be compared, It is difficult to give firm qualitative guidance since the accuracy required depends on the nature of the analysis. A change in stress of more than +/- 20A would be regarded as unsatisfactory for design purposes.
4.4.1
Stress
Components and rotations). These and hence the stresses.
The unknowns solved for in FEA are displacements (translations displacements are then used to calculate strains in the element,
For some element types intermediate steps are involved, The nature of inter-element stress discontinuities depends on the element type concerned. In one-dimensional elements such as truss and beam elements, there are no discontinuities because the displacement functions are sufficiently detailed. For example, the standard beam element is based on cubic displacement and hence can represent linear variations of bending moment. Two and three-dimensional lower order elements generally have discontinuities in the stress field at element boundaries unless they are in a constant stress field. For plane and solid elements, stresses depend on displacement plate bending elements. The stress state at a point is defined element type, These are summarized derivatives, and on curvature for
by several stress components in Table 3-4-1.
depending
on the
ELEMENT TYPE
STRESSES
Truss Beam Plane Element Plate Bending Solid TABLE 3-4-1 Stresses
ax ox, TY, T, % ~Y/ TW OX, OY, T,, (Top & Bottom) UX, Ov, ~,, TN, TV,, T= Represented by Element Type
The state of stress in plated and solid structures is generally quite complex, and has to be combined in some way for design situations. Many failure theories have been developed wherein failure is said to have occurred when some equivalent stress exceeds the yield stress. The equivalent stress combines all the stresses acting at a point in the material. The most popular of these is the Von Mises stress which is given by:
3-45
%=
(Oy-q)+ (w)} +6 (fy+fz+fi)l
The use of the equivalent appropriate.
stress for checking the critical buckling stress is not normal stress (OX,OY)and shear stress (Txy), as
For buckling checks,
appropriate, should be used. Generally normal stresses will not be uniform across the panel, Where this is the case, it will be necessary to approximate the stress by a linear distribution for which there are standard buckling formulae. In some cases, the stress state may be biaxial and/or there may be significant shear stresses. To check these situations, it is usual to calculate the ratios of actual stress and critical stress for individual 4.4.2 Average stress states, and combine the effects using interaction formulae.
and Peak Stresses elements, each stress component for each element have been elements may
Except for the one-dimensional
meeting at a node will be different, In FEA programs various techniques developed to average stresses, The stresses in four adjacent membrane look something like the distribution depicted in Figure 3-4.1.
FIGURE 3-4.1
Distribution
of Element Stresses
Stresses can be calculated at any point in the element. It has been shown, however, that depending on the element formulation there are optimal points for computing stresses. In general, stresses are least accurate at corners, more accurate at mid sides, and most accurate at certain interior points. For two and three-dimensional elements based on the isoparametric formulation (by far the most popular) these interior points
3-46
are the so-called Gauss points (integration points), One popular method is to extrapolate the stresses calculated at the Gauss points to the nodes using a more suitable formula than the actual interpolation functions such as, for example, least squares, However, in some FEA software, the values at the Gauss points are copied to the nearest node without extrapolation, unless otherwise instructed, There are yet other methods for estimating nodal stresses. Once the nodal stresses have been calculated for all elements contributing to the node, they can be averaged to yield an average nodal stress. This will be done for all appropriate stress components, Averaged nodal stresses are much more reliable than element nodal stresses, although the extent of the stress discontinuity at the nodes should decrease The different with mesh refinement. used by FEA software systems for extrapolating Gauss point
methods
stresses to the nodes is perhaps the main reason analyses of the identical problem, using different systems, can yield identical displacement results yet differing stress results. One technique used to overcome this problem is to employ dummy line elements in critical regions of structure. In this technique a dummy truss element is included in the model in the area of interest. An example of such a situation is the placement of such an element at the edge of an opening. The stress results from the truss element are directly calculated and are not dependent on extrapolation. The area of the truss element should be small enough to have negligible influence on response. An area of t2/1 00, where t is the thickness of the plate, is a reasonable upper bound. The use of such elements in the interior of plated structure, or indeed any structure, should be undertaken with caution. direction of the axis of the element. the direction of principal stress. Line elements will yield only normal stresses in the In general line elements will not be aligned with
The current popularity of producing smoothed stress fields in stress plots have hidden dangers, It hides large disparities in stress in adjacent elements, Large disparities indicate too coarse a mesh, A more revealing plotting technique is stress contours. These should be smooth and not jagged. It is evident from Figure 3-4.2 that the contours in the coarse mesh are not smooth, This might be regarded as an unacceptably coarse mesh. An even more revealing method with modern postprocessing checkerboard systems is stress isoband plots. These plots will show a type of distribution for unacceptable stress distributions.
The stress results from a FEA undertaken in support of design are often plot-ted in terms of Von Mises stresses, although principal stresses and component stresses are, also sometimes plotted. There are two potential pitfalls that should be guarded against in interpreting stresses: 1. At nodes on boundaries between membrane elements of different thickness stresses, of course, cannot be simply averaged. A check should be made to ensure that the software does not perform averaging blindly in such a configuration. Care should be taken in interpreting stresses at nodes where two-dimensional elements are not in the same plane. Clearly simple averaging is not appropriate.
2,
3-47
!,
\-4,,p.,
FIGURE 3-4.2
Stress Contours
in Coarse and Fine Meshes
4.5 4.5.1
Other Natural
Results Frequencies and Modes
A feature of the finite element method is that the lower vibration modes are more
accurately determined than higher modes. The curvatures in structures in higher modes are more severe than at lower modes, and several masses are required to represent the kinetic energy accurately at higher modes. These features conspire to make the accurate prediction of higher modes difficult.
In assessing the results from a dynamic analysis, a good starting point is the value of frequency, As an approximate guide, the following may be used for the first few modes: 1. 2. 3. 4. Hull Girder Main Mast Superstructure Typical Stiffened 1- 5Hz 5- IOHZ 1O-2OHZ 1O-4OHZ
Plate Decks
The reliability of higher vibration modes can be assessed by considering the number of masses represented in the lobe of a mode shape, Figure 3-4.3 illustrates this idea.
3-48
,\ ~., ..-,,
SIX MASSES
IN LOBE - GOOD REPRESENTATION
TWO
MASSES
IN LOBE - POOR REPRESENTATION
~f~
FIGURE 3-4.3
Assessing Accuracy
of Higher Modes
3-49
5.0
CONCLUSIONS This section
CHECKS and recommendations, of a finite
deals with the final phase, conclusions
element analysis (FEA), It is necessary to perform these loading, strength, and acceptance criteria are considered This is a critical aspect of a finite element analysis since typically be based on recommendations contained in this sections are grouped conclusions. 5.1 FEA Results into five subsections
checks to ensure that the in arriving at the conclusions. engineering decisions will section, The following of FEA
dealing with various aspects
and Acceptance
Criteria quality assessment checks have
A statement confirming that all analysis procedure been executed satisfactorily should be included.
Finite element analysis is an approximate solution technique, and, in spite of careful effort, the results can only be approximations of the real solution, Therefore, the FEA results should always be validated using an alternative method/s. Alternative methods include comparison with experimental data, approximate analytical models, text book and handbook cases, preceding numerical analyses of similar problems, numerical analysis of a related but simpler problem, and results for the same problem predicted by a different program (which could be based on a different numerical method). Many closed-form solutions of structures with simple geometry are available in handbooks and manuals, which could provide a good means for comparison. Numerical analysis using FEA of similar but simpler models could also be used for comparison An example could be the use of a grillage model to check the results of a finite element model of typical deck structure, Despite the remarks made in the previous paragraph the results from alternative solution
methods should also be treated cautiously. Analytical models incorporate idealizations, mistakes may be made in the calculations, textbooks and handbooks may contain errors, numerical solutions are subject to errors in coding and in data preparation, and experiments may be improperly performed and the results misinterpreted. Therefore, when the FEA results do not compare well with alternative methods, the possible reasons should be investigated. The results should be presented so that they can be easily compared with the design/acceptance criteria. Finite element analysis results are identified based on node numbers and element numbers. These should be translated into the actual physical problem. For example, in a lattice mast, the members that do not meet the safety requirements should be highlighted on a figure of the model for easy identification. When the FEA results do not meet the acceptance criteria, possible reasons should be In case of large deviations, further justification regarding the explored and documented. validity of the FEA results should be provided. The results should be assessed based on the knowledge of the physical problem. For analyses of high category of importance, an independent assessment should always be done by a qualified and experienced person.
3-50
[,
.,, ,,
5.2
Load Assessment In case of discrepancies in the results, the loading applied to the model should be reviewed as part of the investigation into the source of the problem. The appropriateness of the types of loads, load cases, magnitudes, directions, load combinations, load factors, boundary conditions, etc., should be reviewed. The loads applied to a finite element model are approximations contractor should provide a general description on the method of the actual loads. used to approximate The the
actual loads, If the load distribution is simplified to a more regular or uniform distribution, this should be justified to ensure that the simplified load distribution closely approximates the actual distribution in magnitude and direction. For example, if concentrated forces, at nodes, are used to approximate a pressure distribution, the calculations used in assigning the values of nodal forces should be explained. When concentrated forces are used to duplicate pressure, it is important that the load is applied such that the resultant acts through the centre of pressure, on
Details on load factors used in the analysis should also be provided, The information whether the loads are based upon serviceability limit states or ultimate limit states should also be provided, Finally, an assessment of the accuracy the results from the analysis. 5.3 Strength/Resistance Assessment
of the applied loads should be used in describing
In design situations using traditional methods the practice is to apply a nominal design load to the structure and compare the computed stress with some allowable stress. The latter is usually some fraction of the yield stress or the theoretical buckling stress, In the modelling process several assumptions are made which may, or may not be, conservative. An assessment of the conservatism, or otherwise, should be made particularly in regard to the underlying assumptions implicit in the design criteria that are being applied. Often design criteria have evolved with design methods based on hand calculation, Different design criteria may be approrpiate if FEA is used to compute stresses. This factor should be included as part the strength/resistance assessment. In making an assessment of the strength/resistance of the structure based on the results of a FEA, appropriate allowances should also be made for factors that were not accounted for in the analysis, Some of these factors include geometric and material imperfections, misalignments, manufacturing tolerance, initial strains, and corrosion. The design criteria being applied may implicitly include an allowance for some, or all, of these factors. 5.4 Accuracy Assessment
In assessing the accuracy of FEA results, factors to be considered include: the level of detail and complexity modelled, type of behaviour modelled, mesh refinements, etc. In deciding the level of detail the analyst would necessarily have omitted some elements
3-51
of the structure.
The effect
of these on the results
should be assessed.
The limitations
of the element type/s used should also be assessed with respect to its capacity to For example, the element type used might model only model the required behaviouri the membrane actions when both membrane and bending behaviour are significant. The joints and connections between members might not be properly detailed in the model, making the model behave in a significantly different way. The effect of the mesh density used on the results should also be assessed. Simple parametric studies on smaller models may sometimes be necessary to assess the accuracy of the mesh used in the model. Performing checks on the numerical accuracy of an FEA is difficult. Generally reliance is placed on a combination of following good modelling practice and on parameters Common parameters output include the ratio of the largest output by the FEA program. and smallest stiffness found in the stiffness matrix, and the so-called residua/s. Unfortunately, satisfactory values for these parameters are necessary, but not sufficient, conditions for satisfactory numerical performance. depends provided
The acceptability, or otherwise, of the ratio of the largest to smallest stiffness on the computer hardware and software and it is suggested that the guidance by the warning and error messages issued by the FEA program are heeded.
5.5
Overall
Assessment
All of the above described factors should be used in conducting an overall assessment of the FEA. The results of this overall assessment should be included as part of the documentation. Deviations, if any, from the actual response should be justified. Recommendations, if any, for future FEA should be clearly stated, anticipated continuation for the project at a later date, information documentation, etc. should be documented, If &here is an on all computer
files,
3-52
PART BENCHMARK PROBLEMS
FOR ASSESSING
FEA SOFTWARE
1.0
INTRODUCTION The assessment methodology presented in Part 2 includes a requirement that suitable
FEA sollware be used. The determination of the suitability of a particular FEA code should involve, among other things, an assessment of its capability to analyze the types of problems that will be applied. This part describes the development and application of a series of standard benchmark test problems that can be used to assess the suitability of new, or significantly modified, FEA software for ship structure analysis, As a means of qualifying FEA sottware, the benchmarks represent a category of test between that of large scale validation efforts and that of smaller scale verification problems. The actual structural behaviour of even the simplest component depends on such a large number of variables of varying complexity, that isolating the response modelled by FEA codes is extremely difficult. As such, large scale validation of FEA software is typically very complex and expensive, of-ten requiring comparison of FEA predictions with physical test results, Although such validation testing may be a requirement for certain critical structure applications, it is not a practical approach for assessing FEA software on a routine basis.
Most FEA software developers perform verification tests as part of their internal quality assurance procedures. For example, the verification test set for the ANSYS FEA program consists of over 5500 test cases at revision 5.1, Some software developers publish and / or make available a subset of the tests in the form of examples or verification manuals. Other developers include text book verification examples in their marketing media. Verification problems of this sort are usually simple and small-scale in character and typically have closed-form theoretical solutions. They are generally designed to test a very specific aspect of the FEA coder such as the numerical performance of a certain type of element in a certain geometry, loading condition and type of analysis. However, the verification problems rarely resemble real life engineering problems involving irregular geometries with large numbers of element types, in various shapes and sizes, combined with several load types and boundary conditions. Thus, while verification problems of the type described above are a necessary step in verifying and validating FEA software, they are not sufficient on their own, The benchmark problems presented here are intended to represent the next step in ensuring that the candidate FEA software is appropriate for the FEA of linear elastic ship structure. The benchmarks are summarized in Table 4.1-1 and cover a range of typical problems and requirements encountered in real life ship structure FEAs. The problems involve simple configurations of a number of representative ship structures, but are detailed enough to retain the key characteristics of the structural assembly or detail. Tha problems typically require that several types of elements, materials, and loads be used in combination. An attempt has been made to design the benchmarks such that, collectively, all key features that determine the quality of FEA packages are
4-1
addressed. The benchmark problems details given in Appendix D.
are described
in Part 4, Section
2,0 with complete
The benchmarks are designed to exercise the FEA software rigorously without making the evaluation process overly demanding. The problem size has been limited to a maximum of 200 nodes to ensure that the process of benchmarking new and modified software is. not onerous, The 200 node limit should also allow, in some cases, for the user to test demonstration or evaluation versions of FEA software. Such versions are usually based on the full versions of the FEA coder but typically have limits on the number of nodes and elements that can be modeled. These are usually available from the FEA software developer at a small nominal fee to allow testing and evaluation prior to making a larger financial commitment, The benchmarks do not have closed form theoretical solutions. Instead, the results from analyzing the benchmark problems using three well known FEA software programs are used to establish the reference benchmark results, The three programs used were ANSYS, MSC / NASTRAN, and ALGOR and are described in Part 4, Section 3.0. Presentation and discussion of the benchmark results is included in Appendix D.
Care has been taken to ensure that the test models for the benchmark problems are sufficiently detailed or refined that the results approach a converged solution, Element formulations, stress averaging / extrapolation algorithms, and other aspects of FEA software performance tend to be optimized for ideal configurations. Testing different FEA software of an ideal configuration (e.g. a rectangular plate with uniform rectangular elements) will tend to give virtually identical results, However once the FEA model deviates from an ideal configuration, as is the case for the benchmarks, differences in the results manifest themselves, In these circumstances the rate of convergence of results from different FEA programs may differ, Ensuring that the results obtained by the test models are near a converged solution should minimize any discrepancies that can be attributed to poor mesh design of the benchmark test models.
New, or significantly modified, FEA software can be evaluated by exercising the software with the benchmark problems and comparing the results obtained with the reference benchmark results. The process by which this should be accomplished is presented in Part 4, Section 4,0,
WARNING The benchmark problems and associated FEA models presented in this document are intended for the express purpose of evaluating FEA software for ship structural analysis applications. While attempts have been made to ensure that the FBI models follow good modelling practice, they should not necessarily be regarded as appropriate for any other purpose than that for which they are intended.
4-2
Benchmark Problem Features

BM-I Reinforced Opening 2D 3D Analysis Types q q q q q q q q q q q q q q q q q q q q q q BM-2 Stiffened Panel BM-3 Isolation System BM-4 Mast BM-5 Bracket Detail
Static Modal Mass Spring Truss / Spar Element Types

Beam Membrane Shell Brick Force Pressure
q q
LoadTypes
Acceleration Displacement Boundary Conditions Displacement Symmetry Displacement Reactions Results Stress Frequency
q q q q q q q q q q q
q q q q
q q q q
q q q
TABLE 4.1-1
Summary
of Ship Structure
FEA Benchmark
Problems
4-3
2.0
THE BENCHMARK The ship structure
PROBLEMS FEA benchmarks Deck Opening include the following problems :
1 - Reinforced 2345-
Stiffened Panel Vibration Isolation System Mast Bracket Detail
Table 4.1-1 summarizes the main modelling and analysis features that the benchmarks are intended to test. The following sections provide a summary description of the benchmark test problems. Complete details of the benchmark problems are presented in Appendix D.
2.1
BM-I
Reinforced
Deck Opening are among the most commonly encountered sources of high
Openings
and penetrations
stress levels in surface ship structures. In most cases, the openings are reinforced by coamings or insert plates to attenuate the resultant stress concentrations. FEA may be required to evaluate the stress levels and the effectiveness of the reinforcement technique. This benchmark tests the capability of FEA packages to analyze this category of ship structure problem and is shown in Figure 4.2-1. The benchmark tests the FEA programs capability to analyze a plane stress concentration problem using either 4-node or 8-node shell elements. However, it goes beyond the classical hole-in-aplate problem by including two plate thicknesses for the deck and the reinforcement insert plate, and by including stiffeners in the plane of the deck.
FIGURE 4.2-1
Benchmark
Problem BM-1
: Reinforced
Deck Opening
4-4
,.
2.2
BM-2
Stiffened
Panel structural component in ships. This benchmark using various plate
Stiffened
panels are the most common
tests the capability of FEA packages to analyze this type of structure and stiffener element modelling techniques. These include : a) 4-node b) 4-node c) 4-node d) 8-node shell elements shell elements shell elements shell elements for plate and in-plane beam elements for plate and off-set beam elements and
for stiffeners. for stiffeners;
for plate and stiffeners; for plate and stiffeners;
Both static and modal analyses are conducted for each model. The static analysis involves surface pressure loading causing out-of-plane panel bending under symmetric boundary conditions (i.e. quarter model). The modal analysis tests the programs capability for calculating natural frequencies and mode shapes under symmetric and antisymmetric boundary conditions.
FIGURE 4.2-2
Benchmark
Problem BM-2
: Stiffened
Panel
4-5
2.3
BM-3
Vibration
Isolation
System
Vibration isolation systems are often required for ships equipment and machinery, FEA analyses may be used to optimize the isolation system and ensure that vibration and shock design criteria are achieved. This benchmark considers a 12 degree of freedom system consisting of a generator which is mounted and isolated on a raft structure which is, in turn, isolated from the foundation structure. The problem is summarized in Figure 4.2-3. Some of the key testing features include of this benchmark include : . .
q
Modal analysis; Point mass including rotational inertia terms (to model generator) Spring elements with stiffness in three directions; and Rigid beam elements connecting generator mass and isolator springs to raft.
,- --
I
-- --
...
- --
-. .._
lm=1800kg : lam= 90 kg m2 ilW=350 kg@ ,l~=370kgm2
+x
Q
I
I I
M
7 d
# nII n
In
l;ll
Ml
Mass Rigid Links \
I
1 I [
Bsams (Saction PfOPnflY 2)
springs
FIGURE 4.2-3
Benchmark
Problem BM-3
: Vibration
Isolation System
4-6
2.4
BM-4
Mast
Structure
loads (wind and Mast structures on ships must be designed to withstand environmental naval ships usually have additional requirements for resisting ship motions). Masts on shock and blast loading. The mast benchmark problem is summarized in Figure 4.2-4 and the key modelling
q
and testing features
include :
Beam elements (with axial and bending stiffness) for main legs and polemast; Axial line elements (spar, truss, rod) for braces; Point mass elements for equipment payloads; Inertial loading in three directions combined with nodal force loading; Two materials (steel and aluminum); Modal analysis.
it can be used to While the benchmark problem is that of a lattice mast structure, assess the FEA programs capabilities for modelling similar frame or truss like structures such as booms and derricks, especially where beam and spar elements are used in
combinations.
FIGURE 4.2-4
Benchmark
Problem BM-4 4-7
: Mast Structure
2.5
BM-5 Welded
Bracket
Connection
Detail loading. Poorly designed or
connection
details on ships are subject to fatigue
constructed details can lead to premature fatigue failure. Finite element methods are frequently used to calculate fatigue stresses and to aid in the development of improved detail geometry and configurations. This benchmark problem is summarized in Figure 4.2-5, Some of the key modelling and testing features of this benchmark include :
q
3-D geometry
with shell elements
of varying thicknesses;
.
q
Axial line elements for bulkheads, deck and flange of bracket; Transition from coarse to fine mesh at the bracket weld; Prescribed non-zero nodal displacement boundary conditions. for a detail
The latter feature was included since in many cases the boundary conditions FEA are obtained from displacements and loads derived from a global FEA.
This particular bracket detail problem is complicated by the existence of a stress In a linear elastic analysis, the stress singularity at the end corner or toe of the bracket. at this point is theoretically infinite. Refining the finite element mesh gives One method which is commonly progressively higher stresses which are meaningless. used to get around this problem is to use the so called hot spot stress, In calculating the hot spot stress no account is taken of the weld geometry, and in an idealised finite element representation (ignoring the weld) the stress is equal to the value at about one plate thickness from the corner (Chalmers, 1993).
FIGURE 4.2-5
Benchmark
Problem BM-5 4-8
: Bracket Detail
, [ -...-,,
3.0
THE BENCHMARK As previously
TEST
FEA PROGRAMS the benchmark problems do not have readily obtainable
mentioned,
theoretical solutions. Instead, the results from analyzing the benchmark problems using three well known FEA software programs are used to establish the reference benchmark results. The three programs used were ANSYS, MSC / NASTRAN, and ALGOR, The ANSYS FEA program is developed~and marketed by ANSYS Inc. of Houstan, PA. ANSYS is a mature, general purpose FEA program that has been commercially available on various computer platforms since 1970. It includes extensive analysis capabilities, a larger comprehensive library of elements, and extensive pre- and post-processing capabilities, The ANSYS Version 5,1 program was run on a DEC 3000 workstation for the benchmark test cases, The MSC / NASTRAN FEA program is developed and marketed by The MacNeal-
Schwendler Corporation, Los Angeles, CA. Traditionally it has been most widely used by the aerospace industry, having evolved from the National Aeronautics and Space Administration (NASA). MSC / NASTRAN is a very comprehensive and mature FEA program that has been commercially available for several decades. It is to some extent regarded, along with ANSYS, as the industry standard. MSC / NASTRAN For Windows 1,0 on an IBM 486 PC was used for the benchmarks, PA. It
The ALGOR FEA program is developed and marketed by ALGOR Inc., Pittsburgh, was one of the first FEA programs to be developed especially for the personal
computer, and has become one of the most popular FEA programs for PC applications, The program features a relatively wide range of modelling and analysis capabilities.
4.0
APPLICATION
OF BENCHMARKS
FOR ASSESSING
FEA SOFTWARE
The intended application of the benchmarks is to provide a methodology for assessing FEA software, This assessment consists of modelling and analyzing the benchmark problem with the FEA software and comparing the results with those obtained by the reference FEA programs as presented in Appendix D. The data files for the benchmark problems in ANSYS, NASTRAN and ALGOR formats may be obtained by contacting the Ship Structure Committee, As was discovered in the benchmark results of the three reference FEA programs, there are liable to be differences between the results obtained by different FEA software packages. The differences may arise from a multitude of factors ranging from the numerical accuracy of the hardware and software platforms, to different element formulations, solution algorithms, and results presentation techniques, to actual errors or limitations in the FEA software. The question that arises is how much variation or deviation from the reference results is acceptable. The authors suggest the following approach be used to judge the acceptability otherwise of the benchmark results for any FEA software : or
4-9
1.
Result differences less than 2%0 with respect to the reference FEA software results for displacements, reaction forces, and lower mode natural frequencies are considered acceptable. The 2% limit is generally within what would normally be the required engineering accuracy for these types of problems.
2,
Result differences between 2% and 5% are probably acceptable for beam and plate element stress results and higher mode natural frequencies. However the user should endeavour to ensure that there are plausible explanations when differences get much past 2%0, This may involve further testing of the problem by, for example, refining the FEA mesh or switching from the defaults used by The FEA program. the analysis options to /
3.
Result differences greater than 5 YO should be considered as abnormal and require an explanation, If a reason cannot be found, the developer of the FEA software should be contacted and requested to investigate the difference. Where no explanation exists, the FEA software should probably be viewed as suspect for the particular type of analysis covered by the benchmark problem.
Particular attention should be paid to ensure that the proper loads and boundary conditions have been applied, and that the stress contours, deformed shape or mode shapes (depending on what is applicable) are consistent with the reference results. The user should also be sure of the default analysis assumptions and solution techniques used by the software. These can be especially impofiant for problems where transverse shear effects need to be considered, or when performing modal and analyses. The user should also be aware of how the FEA software extrapolates averages plate element stress results at nodes, The benchmarks are a necessary but by no means complete method of validating an FEA program, The benchmarks primarily check that a particular FEA code will perform and produce results that are consistent with the three reference FEA codes. However, it is strongly recommended that users of new or significantly modified FEA software become fully aware of all features and limitations of that program for the particular applications involved. This should include testing the software on simplified versions of the main problems of interest in order to build confidence in the modelling approach, choice of elements, mesh densities, etc. as discussed in Part 3, Section 1,3.
I or
4-1o
.. . ,
PART
CONCLUSIONS
AND RECOMMENDATIONS
From a historical perspective the use of finite element analysis (FEA) as a technique for ship structural analysis is relatively new. In contrast to traditional ship structural analysis and design practice, the application of finite element technology to ship structural analysis is not as well established. As a result the body of experience in the application of this technology is limited. In common with most new technologies FEA is relatively unregulated in terms of the tools that are used in its practice, and the qualifications of organizations and individuals who perform the analysis. This presents a special problem for those that are required to evaluate finite element models and results. The work presented in this report seeks to provide guidance to those that are faced with the problem of evaluating the FEA work performed by other parties. As an aid to the evaluation process a comprehensive and systematic assessment methodology is presented in this report. It is designed to be flexible in terms of the level of skill expected of the evaluator, and in terms of the size and complexity of the FEA that the methodology can be applied to.
The methodology is structured in three levels, The first level is essentially an overview checklist of features of a FEA that need to be evaluated. A more detailed checklist, based on the first level, is presented in the second level of the methodology. The third level provides guidance in narrative and illustrative form, and is structured to match the first and second level checklists. Further guidance is provided through a series of illustrative examples which show the influence of varying finite element modelling practice on FEA results. These are intended to help the evaluator in assessing the levels of accuracy that might be attained in the FEA that is being evaluated. The proliferation of FEA software on the market presents a particular problem for the evaluator, and hence quality of the FEA software is considered to be a key element of the evaluation, While well established FEA software houses follow rigorous comprehensive quality procedures their tests tend to concentrate on small problems, particularly those for which closed-form solutions are available. Benchmark problems of the type presented in this report can be regarded as a further level of qualification. These benchmark problems are intended to test the ability of software to provide accurate solutions for structural assemblies typical of ship structures. Unlike the typical verification problem used by software houses benchmark problems consider non-ideal configurations, multiple element types, several load cases etc. FEA codes are large and complex and hence can never be guaranteed to be free of errors. However, it is suggested that FEA software that has been thoroughly tested by the vendor at the verification example level, will, by successfully yielding solutions for the benchmark problems, provide another level of assurance that the software is fit for performing ship structure FEA.
5-1
._,...
Several 1.
recommendations
are presented
below for consideration:
2, 3.
The assessment methodology as presented is entirely new and can certainly be refined. This is best done by seeking feedback from evaluators of FEAs who have used the methodology. The scope could be broadened to include dynamic response computation, nonlinear behaviour, and composite materials. The benchmarks presented in this report might be considered as a starting point for building a library of benchmark problems, These problems could also include high quality well documented experiments on ship structure assemblies. On a broader front consideration should be given to the important question of design criteria for structure analyzed using FEA. Traditional structural design methods have evolved over many decades of use, and the design criteria used implicitly allow for, among other things, uncertainties associated with the structural analysis and design method used, Compared with traditional structural analysis and design methods the finite element method has quite different capabilities, and limitations. The subject of structural design criteria when the analysis is based on FEA should be the subject of investigation and research,
4.
5-2
PART
REFERENCES
CHALMERS,
D, W,, Design of Ships Structures,
HMSO,
London,
1993. Displacement
CON NOR, J.J. and WILL, G .T.,Computer-Aided Method, MIT Report 69, Feb. 1969,
Teaching of Finite Element
COOK, R. D., MALKUS, D. S., and PLESHA, M. A,, Concepts and Applications Ana/ysis, Third Edition, John Wiley & Sons, New York, 1989. GIANNOTTI & ASSOCIATES, prepared for the Department IRONS, B., and AHMAD, UK, 1980. ISSC, 1991, International
of Finite Element
IN C., Structura/ Guidelines for A/umerica/ Ana/ysis, report of the Navy, NAVSEA, Washington, DC., USA, 1984. ot Finite E/ements, Ellis Horwood Limited, Chichester,
S., Techniques
Report of Committee 11.2: Dynamic Load Effects, Proceedings of the 11th Ship and Offshore Structures Congress held in Jiangsu, Peoples Republic of
China, 16-20 September 1991, Volume 1, edited by P.H, Hsu and Y.S. Wu, Elsevier Applied Science, London, UK and New York, 1991. KARDESTUNCER, New York, 1987, H. (Editor in Chief), Finite E/ement Handbook, McGraw-Hill Book Companyr
NAFEMS, Quality System Supplement to ISO 9001 Relating to Finite Element Analysis in The Design and Validation of Engineering Products, Ref: ROO13, National Agency for Finite Element Methods and Standards, East Kilbride, Glasgow, UK, 1990. STEELE, J. E., Applied Finite Element Modelling, Marcel Dekker, Inc., New York, 1989.
6-1
<.,, ..
..
6-2
Appendix
Evaluation Forms for Assessment of Finite Element Models and Results
A-1
1- Prellmlnary Cheeks 1.1 Doaumentatlon

Perform these ohecks to ensurethatthe analysisdocumentation, job spetication, FEA smlware, and
Result Preliminary tieoks acceptable ? Yes No
1.2 Job Spaclflcetlon 1.3 Finks Element Anaiyais Software 1.4 Contmctor / Analyst Quallfloatlona
contractor analystqualification / requirements have beanaddressed.
6
Rasult
1
2- Engineering Model Checks 2.1 Analysl$ Type& Assumptions PerFarm thesechecksto ensurethat the assumptions usedto developtha engineering modelof the problem are reasonable.
2.2 Geometry 2,3 Material Properties 2.4 Stiffness & Mass PmperUes 2.6 Dynamic Degrees of Freedom
2.6 Loade & Boundmy Conditions
Engineering model is a~ptable ? Yes No
3. Finite Element Modal Checks

3.1 Elamerd Types
Result
3.2 Me$h Doslgn Performthesechecksto ensurethat the finiteelementmodelis an adequate 3.3 SUhtructuras q Submodela nd interpretation f the engineering o model. 3.4 FE Load=& Boundary Condiflone 3.6 FE Solution Optlona & Procedures
Finitaelementmodel is acceptable? Yas No
4
4- Flnita Elemant Analyale Reeulta Checks
4.1 Ganeml Solution Cheeks Reauit
Performthesechacksto ensurethat the finiteelementresultsare calculated, rooessed p and presented in a mannerconsistent iththe analysis w requirements.
4.2 Post Proceealng Mathoda 4.3 DisplacementReauhs 4.4 Strea8 Reaulls 4.5 Other Resul&
Finiteelement resultsare
S - Conclusions Chacks
i Raautt i
Performthesechecksto ensurethat adequateconsideration f the loads, o strength, cceptsnmcriteria,FE a model,and resultsamxracy are included arriving the ~nclusions in at fromthe finite elementanalysis.
-1
.
5.5 Overall Aeaea8msmt
53
Yes No
Conclusions f o the analysisare aoosptable?
FIGURE 1 Overall Evaluation Methodology
Chart
FINITE ELEMENT ANAL YSIS ASSESSMENT Project No. Contractor Name: Analvst: ~ Checker: Project 77tle:
I pRELIMINARy
CHECKS
Date:
1.1
Documentation
Requirements Refer to Guideline Section 3-1.1
Finite Element Analysis Assessment Check 1.1.1 a) b) c) d) e) f) 9) h) i) j) k) 1)
Result
Comments
Has the following information been provided in the FEA documentation? Objectives and scope of the analysis, and acceptance criteria.
Description of physical problem. Description of engineering model, Type of analysis. System of units, Coordinate axis systems.
Description of FEA model, Plots of full FEA model and local details. Element types and degrees of freedom per node. Material properties.
m) Element properties (stiffness & mass properties). n) o) P) q) r) FE loads and boundary conditions. Description and presentation Assessment of the FEA results.
of accuracy of the FEA results.
Conclusions of the analysis, List of references. Result
Based on the above checks answer Question 1.1 and enter result in Figure 1.0. 1.1 Is the level of documentation sufficient to perform an assessment of the FEA?
Comments
A-3
1.2
Job Specification
Requirements
Refer To Guideline section F
Finite Element Assessment 1.2.1
Check
Result
Comments
Is the job specification identified and referenced in the analysis documentation? Are the objectives and scope of the analysis clearly stated and are they consistent with those of the iob specification? Are the analysis requirements clearly stated and are they consistent with those of the job specification? If certain requirements of the job specification have not been addressed (such as certain load cases), has adequate justification been given?
1.2.2
3-1.2
1.2.3
3-1.2
1.2.4
3-1.2
1.2.5 Are the design / acceptance criteria clearly stated and are they consistent with those of the job specification? 1.2.6 Is there reasonable justification for this problem? 1.2.7 for using FEA
3-1.2
3-1.2
Has advantage been taken of any previous experimental, analytical, or numerical works that are relevant to this problem?
3-1.2
Based on the above checks answer Question 1.2 and enter result in Figure 1.0. 1,2 Does the analysis address the job specification requirements?
Result
Comments
A-4
1.3
Finite Element
Analysis
Software
Requirements
Check
1%5 I Ies
3-1.3 3-1.4
Comments
1.3.1
Is the FEA software on the list of approved programs for ship structural analysis applications?
If the answer to Check 1.3.1 is Y, you may skip Checks 1.3.2 and 1.3.3. 1.3.2 Are the capabilities and limitations of the FEA software used to perform the required analysis stated in the analysis documentation? Is evidence of this capability documented and available for review (eg. verification manual, results of ship structure FEA benchmark tests, wevious amxoved FEA of similar moblems)? Does the vendor of the FEA software have a quality system to ensure that appropriate standards are maintained in software develoDment and maintenance,
1.3.3
3-1.3
1.3.4
Based on the above checks answer Question 1.3 and enter result in Figure 7.0.
I I
Result
1.3
Is the FEA software qualified to perform the required analysis?
Comments
NOTE: Part 4 of this report presents benchmark problems for the purpose of assessing the quality and suitability of FEA software for performing ship structural analysis. On its own, successful performance of the candidate FEA software in exercising the benchmark problems is not sufficient evidence of the quality and suitability of the software. The assessor should, in addition, be able to answer the other questions in the table above affirmatively.
A-5
.,
.,,
..
1.4
Contractor
/ Personnel
Qualification
Requirements
Finite Element Assessment
Check
Refer To Guideline Section
Result
Comments
1.4.1
Do the contractor personnel have adequate academic training and experience qualifications to perform finite element analysis?
3-1.5
1.4.2 Do the contractor personnel have adequate engineering experience qualifications for performing ship structural design or analysis? 1.4.3 Do the contractor and contractor personnel have adequate professional certification qualifications? Does the contractor have a working system of Quality Assurance (QA) procedures and checks that are satisfactory for the requirement? Do the contractor personnel have adequate experience with the FEA software used for the analysis?
3-1.5
3-1.5
1.4.4
3-1,5 t3-1.5
1.4.5
Based on the above checks answer Question 7.4 and enter result in Figure 1.0. 1.4 Is the contractor adequately qualified for performing ship structure FEA?
m
I
A-6
FINITE ELEMENT ANAL YSIS ASSESSMENT Project No. Contractor Name: Analyst: Checker: I F701ect 17tie:
ENGINEERING
MODEL
CHECKS
Date : <
2.1
Analysis
Type
and Assumptions Refer To Guideline Section 3-2.1
Check
Result
Comments
Does the engineering model employ enough dimensions and freedoms to describe the structural behaviour (egt 1-D, 2-D, or 3-D)? Does the engineering model address the appropriate scale of response for the problem (egi global, intermediate, or local response)? Is the type of analysis appropriate for the type of response and loading of interest (eg. linear, static, dynamic, buckling analysis)? Does the engineering model address all the required results parameters (eg. stress, displacement, frequency, buckling load)? Are all assumptions affecting the choice of engineering model and analysis type justified (watch for non-standard assumptions)? Is the level of detail, accuracy or conservatism of the engineering model appropriate for the criticality of the analysis and type of problem? Does the analysis employ a consistent set of units?
Does the analysis coordinate employ a consistent global
2.1.2
3-2.1
2.1.3
3-2.1
2.1.4
3-2.1
2.1.5
3-2.1
2.1.6
3-2.1
2.1.7
3-2.1
2.1.8
3-2,1
axis system?
Based on the above checks answer Question 2.1 and enter result in Figure 7.0. 2.1 Are the assumptions of the type of analysis and engineering model acceptable?
A-7
... ~,.; \ -.,.

,-
2.2
Geometry
Assumptions
Finite Element
Analysis
Assessment
Check
Result
Comments
2.2.1
Does the extent capture paths,
of the model geometry actions, load parameters
the main structural and response
of interest? 3-2.2
2.2.2
Are correct assumptions used to reduce the extent of model geometry (eg. symmetry, boundary conditions at changes in stiffness)? Will the unmodelled structure (ie. outside the boundaries of the engineering model) have an acceptably small influence on the results? Are the effects of geometric simplifications (such as omitting local details, cut-outs, etc. ) on the accuracy of the analysis acceptable ?
2.2.3
3-2.2
2.2.4
3-2.2
2.2.5
For local detail models, have the aims of St. Venants principle been satisfied?
3-2.2
2.2.6
Do the dimensions defining the engineering model geometry adequately correspond to the dimensions of the structure? For buckling analysis, does the geometry adequately account for discontinuities and imperfections affecting buckling capacity?
3-2.2
2.2.7
3-2,2
Based on the above checks answer Question 2.2 and enter result in Figure 1.0. 2.2 Are the geometry assumptions in the engineering model acceptable?
Result
Comments
A-8
..
2.3
Material
Properties
Result
Comments
2.3.1
Are all materials of structural importance to the problem accounted for in the engineering model? Are the assumed behaviors valid for each material (eg. linear elastic, isotropic, anisotropic, orthotropic)? Are the required material parameters defined for the type of analysis (eg, E, v, etc.)?
3-2.3
2.3.2
3-2,3
2.3.3
3-2.3
2.3.4
Are orthotropic and / or layered properties defined correctly for non-isotropic materials such as wood and composites? Are orthotropic properties defined correctly where material orthotropy is used to simulate structural orthotropy (eg. stiffened panels)? If strain rate effects are expected to be significant for this problem, are they accounted for in the material properties data? Are the values of the materials properties data traceable to an acceptable source or reference (eg. handbook, mill certificate, coupon tests)?
3-2,3
2.3.5
3-2.3
2.3.6
3-2.3
2.3.7
3-2.3
2.3.8
Are the units for the materials properties data consistent with the system of units adopted for other Darts of the analvsis?
3-2.3
Based on the above checks answer Question 2.3 and enter result in Fi!qure 7.0. 2.3 Are the assumptions and data defining the material properties acceptable?
Comments
A-9
L.>,.
2.4
Stiffness
and Mass
Properties Refer To Guideline Section 3-2.4
Result
Comments
2.4.1
Are all components that have significant effect on the stiffness of the structure accounted for in the engineering model ? Are the assumed stiffness behaviors valid for each structural component (eg. linear, membrane, bending, shear, torsion, etc.)? Are the required stiffness parameters defined for each component, eg. : Truss members - A Beams, bars - A, IYY,1,,, other Plates, shells - t (uniform or varying) Springs - K (axial or rotational) Do the section properties of stiffeners (where modelled with beams) include correct allowances for the effective plate widths? If torsion flexibility is expected to be important, are torsion flexibility parameters correctly defined for beam sections? If shear flexibility is expected to be important, are shear flexibility parameters correctly defined for beam and/or plate elements?
2.4.2
3-2.4
2.4.3
3-2,4
2.4.4
3-2,4
2.4.5
3-2.4
2.4.6
3-2,4
If mass or inertial effects are not applicableto this problem, proceed to Check 2.4.13 on the following page. 2,4.8 Are all components that have significant effect on the mass of the structure accounted for in the engineering model? Have material properties data for density been defined (see also Check 2,3,3)? 3-2.4
2.4.9
3-2,4
2,4.10 Has the added mass of entrained water been adequately accounted for with structure partially or totally submerged under water? 2.4.11 Are lumped mass representations of structural mass and / or equipment correctly consolidated and located? If rotational inertia is expected to be important, are mass moments of inertia properties correctly defined for masses?
3-2.4
3-2.4
2.4.12
3-2.4
A-10
finite Element Analysis Assessment Check 2.4.13 Are the values of the stiffness and mass properties data supported by acceptable calculations and / or references? If relevant, has fluid-structure interaction been accounted for? Has the added mass been included in the model? Are the units for the stiffness and mass properties data consistent with the system of units for other parts of the analysis?
Result
Comments
2.4.14
3-2.4
2,4.15
3-2,4
Based on the above checks answer Question 2.4 and enter result in Figure 1.0. 2.4 Are the assumptions and data defining stiffness and mass properties acceptable? Comments
Result
A-1 1
2.5
Dynamic
Degrees
of Freedom
If the analysis is not a reduced dynamic analysis, you may proceed directly to Part 2.6.
Result
Comments
2.5.1
Are dynamic dof defined in enough directions to model the anticipated dynamic response behaviour of the structure? Are the number of dynamic dof at least three times the highest mode required (eg. if 30 modes required, need at least 90 dof)? Are the dynamic dof located where the highest modal displacements are anticipated?
2.5.2
3-2,5
2.5.3
3-2.5
2.5.4
Are the dynamic dof located where the highest mass-to-stiffness ratios occur for the structure? Are dynamic dof located at points where forces or seismic inputs are to be applied for dynamic response analyses? Are the number of dynamic dof such that at least 90% of the structural mass is accounted for in the reduced model in each direction?
3-2.5
2.5.5
3-2.5
2.5.6
3-2.5
Based on the above checks answer Question 2.4 and enter result in Figure 1.0. 2.5 Are the assumptions and data defining dynamic degrees of freedom acceptable? Comments
Result
A-12
2.6
Loads and Boundary
Conditions
Result
Comments
2.6.1
Are all required loadings / load cases accounted for, and has sufficient justification been provided for omitting certain loadings? Are the loading assumptions stated clearly and are they justified?
3-2.6
2.6.2
3-2,6
2.6.3
Has an assessment been made of the accuracy and / or conservatism of the loads?
3-2.6
2.6.4
Are the procedures for combining loads / load cases (eg. superposition) adequately described and are they justified? Have the boundary conditions assumptions been stated clearly and are they justified?
3-2.6
2.6.5
3-2.6
2.6.6
Do the boundary conditions adequately the anticipated structural behaviour?
reflect
3-2.6
2.6.7
Has an assessment been made of the accuracy of the boundary conditions, and if thev movide a lower or urmer bound solution?
3-2,6
Based on the above checks answer Question 2.6 and enter result in Figure 1.0. ) 2.6 Are the assumptions and data defining loads and boundary conditions reasonable? Comments
I I
Result
A-13
.....
FINITE ELEMENT ANAL YSIS ASSESSMENT Project No. Contractor Name: Analyst: I Checker: ~ F!rojectTitle :
I FINITE ELEMENT MODEL CHECKS
Date :
3.1
Element
Types Refer To Guidaline Section

3-3.1
Result
Comments
3.1.1
Are all of the different types of elements used in the FEA model identified and referenced in the analysis documentation? Are tha element types available in the FEA software used appropriate to ship structural analysis? Do tha element types support the kind of analysis, geometry, materials, and loads that are of importance for this problem? If required, do the selected beam element types include capabilities to model transverse shear and / or torsional flexibility behaviour? If required, do the selected beam element types include capabilities to model tapered, off-set or unsymmetric section properties? If required, do the selected beam element types include capabilities for nodal dof end releases (eg. to model partial pinned joints)? If required, do the selected plate element types include capabilities to model out-of-plane loads and bending behaviour? If required, do the selected plate element types include capabilities to model transverse shear behaviour (ie, thick plate behavior)? If the model is 2-D, are the selected element types (or options) correct for plane stress or plane strain (whichever case applies)? If required, can the selected element types model curved surfaces or boundaries to an acceptable level of accuracv?
3.1,2 3.1.3
3-3.1 3-3,1
3.1,4
3-3.1
3.1.5
3-3.1
3.1.6
3-3.1
3.1.7
3-3.1
3.1,8
3-3.1
3.1.9
3-3.1
3.1.10
3-3.1
Basedon theabovechecks answerQuestion 3.1 and enterresultin Egure 1.0. h
!
I
Result
3.1
Are the types of elements used in the FEA model acceptable?
Comments
A-14
3.2
Mesh
Design Refer To
Finite Element Analysis Assessment Check 3.2.1 Does the mesh design adequately reflect the geometry of the problem (eg. overall geometry, stiffener locations, details, etc.)? Does the mesh design adequately reflect the anticipated structural response (eg. stress gradients, deflections, mode shapes)? Are nodes and elements correctly located for applying loads, support and boundary constraints, and connections to other parts? Does the analysis documentation state or show that there are no illegal elements in the model (ie, no element errors or warnings)? Are the element shapes in the areas of interest acceptable for the types element used and degree of accuracy required? Are mesh transitions from coarse regions to areas of refinement acceptably gradual? Are element aspect ratios acceptable, particularly near and at the areas of interest? Are element taper or skew angles acceptable, particularly near and at the areas of interest? If flat shell elements are used to model curved surfaces, are the curve angles < 10 for stresses, or < 15 for displacement results? If flat shell elements are used for double or tapered curve surfaces, is warping avoided (egq small curve angles, use of triangles)? Is the mesh free of unintentional gaps or cracks, overlapping or missing elements? Is proper node continuity maintained between adjacent elements (also continuity between beam and plate elements in stiffened panels)? Are the orientations of the beam element axes correct for the defined section properties?
Guideline
Result
Comments
Section 3-3,2
3.2.2
3-3.2
3.2.3
3-3,2
3.2.4
3-3.2
3.2.5
3-3.2
3.2.6
3-3.2
3.2.7
3-3.2
3.2.8
3-3.2
3.2.9
3-3.2
3.2.10
3-3.2
3.2.11
3-3.2
3.2.12
3-3.2
3.2.13
3-3.2
A-1 5
---
,,
Finite Element
Analysis
Assessment
Check
Refer To Guideline
Result
Comments
3.2.14
Are differences in rotational dof / moment continuity for different element types accounted for (eg. beam joining solid)? Are the outward normals for plate / shell elements of a surface in the same direction?
3-3.2
3.2.15
3-3.2
Based on the above checks answer Question 3.2 and enter result in Figure 1.0. 3.2 Is the design of the finite element mesh acceptable?
Result
Comments
A-16 -..>
3.3
Substructures
and Submodelling Refer To Guideline Section 3-3.3 I
Finite Element Analysis Assessment Check 3.3.1 Is the overall substructure or submodelling scheme or procedure adequately described in the analysis documentation? Are all individual substructure models, global models and refined submodels identified and described in the analysis documentation? Are the master nodes located correctly and are the freedoms compatible for linking the substructures? Are the master nodes located correctly for application of loads and boundary conditions upon assembly of the overall model? Are loads and boundary conditions applied at the substructure level consistent with those of the overall model? Does the boundary of the refined submodel match tha boundary of coarse elements / nodes in the global model at the region of interest? Is the boundary for the submodel at a region of relatively low stress gradient or sufficiently far away from the area of primary interest? Does the refined submodel correctly employ forces and / or displacements from the coarse model as boundary conditions? Does the submodel include all other loads applied to the global model (egi surface pressure, acceleration loads, etc.)? Have stiffness differences between the coarse global mesh and refined submodel mesh been adeauatelv accounted for?
Result
Comments
3.3.2
3-3,3
3.3.3
3-3.3
3.3.4
3-3.3
3.3.5
3-3.3
3.3.6
3-3,3
3.3.7
3-3.3
3.3.8
3.3.9
3-3.3
3,3.10
3-3.3
Based on the above checks answer Question 3.3 and enter result in Figure 1.0. 3.3 Are the substructuring or submodelling procedures acceptable?
G I
Comments
A-17
,.,..
3.4
FE Model
Loads and Boundary
Conditions
Result
Comments
3.4.1
Are point load forces applied at the correct node locations on the structure and are they the correct units, magnitude, and direction? Are distributed loads applied at the correct locations on the structure and are they the correct units, magnitude and direction? Are surface pressure loads applied at the correct locations on the structure and are they the corract units, magnitude and direction? Are translational accelerations in the correct units, and do they have the correct magnitude and direction? Are rotational accelerations the correct units, magnitude and direction and about the correct centre of rotation? Are prescribed displacements applied at the correct locations on the structure and are they the correct units, magnitude and direction. Are the displacement boundary conditions applied at the correct node locations?
3.4.2
3-3,4
3.4.3
3-3.4
3.4.4
3-3.4
3.4.5
3-3.4
3.4.6
3-3,4
3.4.7
3-3.4
Based on the above checks answer Question 3.4 and enter result in Figure 1.0. 3.4 Are the FE loads and boundary conditions applied correctly? Comments
E I
A-18 ..
3.5
Solution
Options
and Procedures
3.5.1
Have any special solution options and procedures been used and, if so, have they been documented? If non-standard options been invoked have they been documented and the reasons for their use been explained? If the problem is a dynamic analysis is the method for eigenvalue and mode extraction appropriate?
3.5.2
3-3.5
3.5.3
3-3.5
Based on the above checks answer Question 3.5 and enter result in Figure 1.0. 3.5 Are the solution o~tions and rwocedures followed for the FEA acceptable?
Result
A-19
FINITE ELEMENT ANAL YSIS ASSESSMENT
FINITE ELEMENT RESULTS CHECKS
Project No. Contractor Name: Analvst:
Project Title: Date: Checker:
4.1
GeneraI
Solution
Checks
Result
Comments
4.1.1
Are all error and warning messages issued by the software reviewed and understood?
3-4.1
4.1.2
Is the magnitude of mass of the finite element model approximately as expected?
3-4.1
4.1.3
Is the location of centre of gravity of the model, as calculated by the program, reasonable? Are the applied forces in equilibrium with the applied reactions?
3-4.1
4.1.4
3-4.1
Based on the above checks answer Question 4.1 and enter result in Figure 1.0. h 4.1 Are the general solution parameters acceptable?
l==
1 I
Comments
A-20
L>,
,,
4.2
Post Processing
Methods
Check
Result
Comments
4.2.1
Are the methods for reducing analysis results described (eg. calculation of safety factors and other parameters calculated by manipulating raw output)? Are the methods for correcting FE results described (eg. correction factors, smoothing factors) ?
4.2.2
3-4.2
Based on the above checks answer Question 4.2 and enter result in Figure 1.0. 4.2 Is the methodology used for post processing the results satisfactory?
l== I
Comments
A-21
4.3
Displacement
Results
Check
Result
Comments
4.3.1
Are the displacement discussed?
results described and
3-4.3
4,3.2
Are plots of the deformed structure (or mode shape) presented?
3-4.3
4.3.3
Are the directions of displacements consistent with the geometry, loading and boundary conditions? Do the. magnitudes sense? of displacements make
3-4.3
4.3.4
3-4.3
4.3.5
Is the deformed shape (or mode shape) smooth and continuous in area of interest?
3-4.3
4.3.6
Are unintentional slits or cuts (indicating elements not connected where they should be) absent?
3-4.3
Based on the above checks answer Question 4.3 and enter result in Figure 1.0. 4,3 Are displacement results consistent with expectations?
l== I
Comments
A-22
4.4
Stress
Results
Finite Element
Analysis
Assessment
Check
Result
Comments
4.4.1
Are the stress results described and discussed? Are stress contour plots presented? In the stress plots are the stress parameters or components defined (eg. u,, OY,Txy, etc.)? Is the method of smoothing stress results, or averaging stress results described (eg. element stresses vs nodal average stresses)? Are the units of stress parameters consistent? Are the magnitudes of stresses consistent with intuition? In cases where there are adjacent plate elements with different thicknesses does the method for averaging stresses account for the differences? Are the stress contours smooth and continuous, particularly in region of primary interest ? Are the stress contours at boundaries consistent with the boundary conditions applied (eg. stress contours perpendicular to boundary if symmetry be)? Are stresses local to the applied loads reasonable? Are there areas in which stresses are above yield (which would invalidate linear elastic analysis)?
3-4.4
4.4.2
3-4.4
4.4.3
3-4.4
4.4.4
3-4.4
4.4.5
3-4.4
4.4,6
3-4.4
4.4.7
3-4.4
4.4.8
3-4.4
4.4.9
3-4.4
4.4.10
3-4.4
Based on the above checks answer Question 4.4 and enter result in Figure 1.0. 4.4 Are stress results consistent with expectations? Comments
Result
A-23
4.5
Other
Results
Check
Result
Comments
4.5.1
Are the frequencies units?
expressed in correct
4.5.2
Are the magnitudes of natural frequencies consistent with the type of structure and mode number? Are the mode shapes smooth?
3-4.5
4.5.3
3-4.5
Based on the above checks answer Question 4.5 and enter result in Figure 1.0. 4.5 Are dynamics results consistent with expectations? Comments
Result
A-24
FINITE ELEMENT ANAL YSIS ASSESSMENT Project No. Contractor Name: Analvst: Checker: Project Title :
CONCLUSIONS
CHECKS
Date :
5.1
FEA Results
and Acceptance
Criteria
Refer To Guideline
Result
Comments
5.1.1
Are the results summarised in a manner that allows comparisons with acceptance criteria, or alternative solutions or data? Are satisfactory explanations provided where the results do not meet acceptance criteria, or where they differ significantly from other com~arable solutions or data?
5.1.2
Based on the above checks answer Question 5.1 and enter result in Figure 1.0. 5.1 Are the results presented in sufficient detail to allow comparison with acceptance criteria?
Comments
A-25
i .,,
k-..
5.2
Load Assessment
Result
Comments
5.2.1
Has an assessment been made of the accuracy or degree of conservatism of the loads used in the FE model with respect to the following aspects :
a)
types of loads / load cases that were included and excluded basis or theory used to derive loads (eg. linear strip theory for sea motion loads, base acceleration vs DRS for shock, drag coefficients for wind loads, etc.) magnitudes of loads
b)
c) d) e) f) g)
loading directions included / excluded load combinations load factors boundary conditions
Based on the above checks answer Question 5.2 and enter result in Fiaure 1.0. 5.2 Are the accuracy and conservatism, understood? Comments or otherwise, of the applied loading modelled
A-26
5.3
Strength
/ Resistance
Assessment
Check
Refer To Guideline
Result
Comments
5.3.1
Has an assessment been made of the accuracy or degree of conservatism of the strength or resistance of the modelled structure with respect to the following aspects : a) failure theory, failure criteria, allowable stresses, safety factors, etc section properties material properties allowances for imperfection, manufacturing tolerances allowances for corrosion misalignment,
b) c) d)
e)
Based on the above checks answer Question 5.3 and enter result in Figure 1.0. 5.3 Has an adequate assessment been made of the capability of the structure?
E
I
Comments
A-27
5.4
Accuracy
Assessment
Check
Refer To Guideline
Section
Result
Comments
5.4.1
Has an assessment been made of the scale of FE model and its level of detail and complexity? Have the types of behaviour modelled and not modelled (egi membrane only instead of membrane plus bending) been assessed? Has the influence of mesh refinement on accuracy been considered?
3-5.4
5.4.2
3-5.4
5,4.3
3-5.4
5.4.4
Has a comparison solutions,
with
other results etc.)
(eg. other
3-5,4
experiment,
been made?
5.4.5
Based on the above of the accuracy made?
has an overall
assessment been
3-5.4
of the relevant
results
Based on the above checks answer Question 5.4 and enter result in Figure 7.0. ) I 5.4 Has an adequate assessment of the accuracy of the analysis been made? Comments
Result I I
A-28
5.5
Overall
Assessment
Check
Result
Comments
5.5.1
Are conclusions from the FEA provided, and are they consistent with the material presented? If appropriate has a way ahead or potential solutions been presented?
5.5.2
3-5.5
5.5.3
Based on consideration of all previous checks is the overall assessment that the FEA is accemable?
3-5.5 i
Based on the above checks answer Question 5.5 and enter result in Figure 1.0. h 5.5 Is the finite element analysis assessed generally satisfactory?
I I
Result
Comments
A-29
.. .
A-30
.,.
.. . ...
Appendix
Example Application of Assessment Methodology
B-1
.,. /.
;
61.0
INTRODUCTION The purpose of this Appendix is to illustrate the application of the FEA assessment methodology and the guidelines presented in Parts 2 and 3 of this document, An example finite element analysis (FEA) of a web frame from an Arctic-going tanker design subject to ice loads is used for this purpose, The approach used to illustrate the assessment
q
methodology
and guidelines
includes : to
a sample report of the Arctic tanker web frame FEA, annotated with references relevent sections of the FEA assessment methodology and guidelines; and completed checklists as required by the assessment methodology. checklists are presented in Annexes
The annotated report and the completed B-4 respectively. 62.0 EXAMPLE FINITE ELEMENT ANALYSIS
B-1 and
The example FEA is adapted from an analysis for an actual designl of an icebreaking tanker. The tanker is double hulled. Transverse strength is provided by a series of closely spaced web frames, and the longitudinal load transfer is achieved through several longitudinal stringers. The design requirements are based on current Canadian rules. The primary interest for this analysis is the behaviour of a typical web frame in response to ice loads, Other loads are ignored as negligible compared with the ice loads. The analysis was performed to ensure that the side structure that directly resists the ice loads responds in the manner expected optimized as possible. This example illustrates structures including:
q
by the designers,
and that the structure
is as
several aspects of finite element
modelling
common
in ship
.
q q q
behaviour of stiffened plate structures openings in structures discontinuities often found in ship structures integrated nature of typical ship structures use of most types of elements commonly used in the FEA of ship structures.
For reasons explained in Annex B-1 it was necessary to make modifications to the original analysis, particularly in regard to the level of ice load, to make it suitable for the purposes of the present work.
1 The design was undertaken by Canarctic Shipping Co, Ltd., Ottawa, to the Transportation Development Centre, Montreal, Quebec, Canada
Ontario, Canada
under contract
B-2
!-
,,
<.._,....
B3.O
ANNOTATED REPORT
Annex B-1 presents a sample report of the Arctic tanker web frame FEA that has been
prepared by a contractor (BB Engineering) and has been subjected to the assessment methodology. For illustrative purposes the report has been annotated with short descriptions identifying the relevant part of the assessment methodology presented in Parts 2 and 3 of this document. Except for the annotations the report is meant to be typical 64.0 of the documentation that an evaluator of FEA might recieve,
CHECKLISTS A sample of completed FEA evaluation presented in Annex B-4. checklists for the report in Annex B-1 are
Acknowledgement The finite element

performed
analysis described
Engineering,
in the following Ottawa,
pages
is adapted
from
an analysis
by MIL Systems
Ontario
for Canarctic
Shipping
Ltd., Centre,
Ottawa, Ontario under a contract Montreal, Quebec.
awarded
by the Transportation
Development
Warning This example is presented solely for the purpose of illustrating the assessment methodology described in Part 2. As such it is not necessarily complete in all details. particularly in regard to parameters such as number of loading types. design criteria, and number of structural responses considered. Furthermore this example should not be construed as representative of the requirements for a finite element analysis of other marine structures.
B-3
..,
~!, ,,-+.
.
Annex
B-1
Finite
Analysis of Arctic Tanker Web Frame

Element
BB Engineering Ltd. 13-1300 Finite Drive Ottawa, Ontario xxx xxx
May B-4
1995
.,. .-,
ANNEX B-1 1.0 2.0 INTRODUCTION PRELIMINARY 281 2.2 2.3 2,4 3.0
TABLE OF CONTENTS
Page No.
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..i. INFORMATION..
..
B. 6-6
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..B-6 . . . . . . . . . . . . . . . . . . . . . . . . . . B-7
Job Specification FEA Software..,,,,,,.. Contractor
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..B-6 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..B-7 . . . . . . . . . . . . . . . . . . . . . . . . . . . ..B-7 . . . . . . . . . . . . . . . . . . . . . . . . . ..B-7 . . . . . . . . . . . . . . . . . . . . . . . . . ..B-8
Rationale for using Finite Element Method and Analyst Qualifications ,,, ,,,
ENGINEERING 3.1 3,2 3.3 3.4 3.5 3,6 3.7 3.8
MODEL.,,,,,.
, . .,, Tanker
Analysis Type and Assumptions Global Geometry of50000DWT Frame Selected Extent of Model interaction Loads
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..B-7
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..B-8 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..B-9 ,,, . . . . . . . . . . . . . . . . . . . .,, . . . . . . . . . . . . . . . ..B-9 . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..B-10 ,,, ,,, ... . . . .. . ..B-12 B.12I2
Material Properties
with Adjacent Structure
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ...6-11 . . . . . . . . . . . . . . . . . . . . . . . . . . . . ,,, . . . . . . . . . . . . . . . . . . . . . . . . ,,, ,,.., . ., .,, , . .,, .,, ,,, . . . . ,,,
Boundary Conditions ...,.
4.0
FINITE ELEMENTMODEL 4.1 4.2 4.3 4.4 4.5 4.6 4.7 General information, Element Selection Mesh Design
. . . . . . . . . ..B-I2
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..B-12 and Spring Constants . . . . . . . . . . . . . . . , . . , . . . B-14
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..B-13 Boundary Conditions . . . . . . . . . . . . . . . . . . . . . . . . ..B-16 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..B-17 . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..B-18 . . . .. . . . . .. B.18 B.. B-18 .B-18
Finite Element Attributes FE Model Loadsand FE Model Checks
FE Solution Option and Procedures RESLILTS ii,,,,,.. Methods,,.
5.0
ANALYSIS 5.1 5.2 5,3
. . . . . . . . . . . . . . . . . . . . . . . . . ., ,, . . . . . . . . . . . . ...-..,,., . ...,,,,...,,,,,,,,,,. ,,, ,,
General Solution Checks Postprocessing Structural Response
B..B-18
. . . . . . ..
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..B-19
6.0 7.0
CONCLUSIONS REFERENCES
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..B-20 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..B-20 . . . . , , , , , , ,., , .,,,,,,,,,,,.,.,,,,, i i . . . . , , . , , . . , . . . B-35 . . . . . . ..AB-35
Annex B-2 Company and Personnel Qualifications B-2.1 B-2.2 Contractor Qualifications,,, Personnel Qualifications
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..B-35 ,,, ,., .,, ,.. .B. .B-36
Annex B-3 FEA Results Verification
. . . . . . . . . . . . . . . . . . . . . ... Methodology Forms
Annex B-4 Sample Completed Assessment
, . . . . . . . . , . . , , . . . . . . . B-37
B-5
~,
,. .,. -. . .
FINITE ELEMENT ANALYSIS SINGLE MIDBODY 1.0

INTRODUCTION Company
OF 50000
DWT TANKER
WEB FRAME
AA Shipping 50000
Limited has developed midbody
a design for a
DWT Arctic tanker.
The focus of the work has been to and bow structures, a
design cost optimized The BE Engineering finite element frame. the response Section
Co Ltd. (BBE) has been tasked to undertake (FEA) of a typical midbody
analysis
web (diaphragm)
The purpose of the FEA reported of the midship structure
in this report is to assess
to ice loads. of the requirements in
2 of this report provides a summary and data on the software The engineering
for the analysis, Section structure model. 5 presents 3,
and the resources model is described of the subject the engineering model, and Section
applied to the problem. This section
includes a discussion
and the assumptions Section
made in developing
4 describes the finite element
the results of the analysis, INFORMATION
2.0
2.1
PRELIMINARY
Job Specification calls for a static, kN, provided in Bow linear elastic, FEA of a web Job Specification Para. 1.2 in the Assessment Methodology
The job specification
frame from the midbody design ice load of 4435 The finite element Arctic Tanker Sections
section of the 50000
DWT tanker at a
model is based on the drawings Evaluation (Ref. 2).
Structural
- Midship Sections,
and Repair Drawings
The acceptance 1. 2. maximum except
criteria for the analysis are as follows: stress not to exceed the material stresses yield stress
Acceptance Criteria Para 1.2.5
as noted in item 2, in excess of yield stress are acceptable
very localized considered
B-6
2.2
Rationale
for using Finite Element
Method Justification for using FEA Para. 1.2.6
The structure concentrations. 2.3 ANSYS element
under investigation particularly
is too complex to be analyzed
by
hand calculation
in regions of high stress
FEA Software finite element by ANSYS software (Version 5.1), here. developed ANSYS and FEA Sotlware Para. 1.3.1
supported established structures
Inc. of Houston, and presented
PA, was used for the finite is a well has a Inc. are have
work performed FEA package and technical updates
has a proven track record in analyzing BBE currently with ANSYS, support contract
of the type under consideration,
maintenance The software reviewed other ANSYS been validated capability 2.4
and error reports received ANSYSS
from ANSYS
by all BBE staff involved documents. against
in FEA, and filed along with shell and beam elements analysis. FEA. designed to test the ANSYS
by BBE for use in ship structural benchmarks to perform
has been evaluated
of the software
ship structural
Contractor
and Analyst
Qualifications the analysts, and Contractor /Personnel Qualification Para. 1.4
Information
on qualifications
of the contractor,
the supervisor,
to perform the required FEA is provided in Annex
B-2 of this document,
3.0
3.1
ENGINEERING Analysis
MODEL
Type and Assumptions are limited to the yield stress the material to be linear. Similarly because behaviour large geometric is assumed to be Analysis Type & Assumptions Para. 2.1
Since the stresses behaviour deflections
is assumed
are not expected
linear as well, The load is assumed strength of the frame. to be static and interest Hence, the dynamic is centred on the behaviour of the frame is behaviour
is not within the scope of this analysis. also not considered considered in this analysis. as part of the design process.
Instability
However,
it should be
B-7
The overall strength analysis, stress concentrations be addressed 3.2 Global
of the frame is the primary focus of this the analysis is not optimized at structural discontinuities to examine such as those
and therefore
that will exist around openings for example.
Again these should
as pati of the normal design process. Geometry of 50000 DWT Tanker length of 242 metres, metres. a Geometry Assumptions Para. 2.2
The 50000 breadth distance 14.6m.
DWT tanker metres
has a waterline
of 34.6 between
and a depth of 18,1 bulkheads dimensions
The vessel Each m x
has seven cargo tanks.
In the cargo tank region of the vessel the is 19.2 metres, of 18 m x 30.6
transverse
cargo tank has approximate
The vessel is double hulled, outer hulls is 2000 4.0 metres Thereforer 11.0 mm. turn of the bilge and connects above baseline. the deck structure metres.
The distance
between
the inner and at a point with of
The bottom
structure
wraps around the connects
to the side shell structure metres above baseline. vertically
The side shell structure
at a point 15,0
the side shell structure The structure spaced at 1000
spans a distance framed
is transversely mm intervals.
by web frames framing bulkheads,
(diaphragms) is provided
Longitudinal
by several stringers spanning is shown
between
The midship section 3.3 Frame
in Figurel 3,1
Selected DWT tanker is approximately positioned, 22 Extent of Model Para. 2.2.1
The ice load for the 50000 metres across a pair of bulkheads. is resisted transverse by the transverse bulkheads. ring), the deck structure,
in length and therefore,
if centrally
spans as a
The ice load applied to side structure frames (each acting essentially structure and by the the bottom
Any transverse vertically and longitudinally
loads applied to the side structure and deck structures through to bulkheads stringers,
are distributed frames,
to the bottom
by transverse
The most severe loading case for a web frame is from ice load
at the end of this document
1 Figures
are presented
B-8
applied to the frame midway disposed load are discussed Extent in Section
between 3.7.
bulkheads
and centrally of the
with respect to the frame.
The characteristics
3.4
of Model of the vessel, between transverse them. bulkheads, is a Extent of Model Para. 2.2.1 An alternative to account influence surrounding adjacent method
The structure stiffened connected members
series of ring frames
comprising
inner and outer hull plating with a These frames are (framing oriented structure
plate diaphragm and plating).
connecting
by all longitudinally
for the of the structure
It is sufficient boundary centreline
to model a single transverse (structure
ring frame if the correct in Section 3.6, Due
conditions
are applied as discussed
would be to model web frames and stringers approximately.
to the symmetry
and load) that exists along the vessel to model one half of the ring frame. of the ship at centreline The width of the mm) and will include plate diaphragm. Figure
it is also sufficient
This ring frame extends
from the bottom
around to the vessel centreline
at the deck,
model needs to be the frame spacing (1 000
the inner and outer shell plating and the stiffened Figure 3.1 illustrates
the midbody frame that was analyzed. for the frame.
3.2 shows the outer dimensions
3.5
Material
Properties that the vessel material and EH36. in the outer shell Table 3.1 lists the 3 for these Material Properties Para. 2.3
Figure 3.2 indicates components relevant steel grades. The Youngs

types. strains were
plating is Grade EH50 material
and that the inner shell and framing as taken from Reference
are Grades DH36 properties
Modulus
was taken as 208,700

as initial imperfections
MPa for all steel

and residual
Parameters
such
not included
in the analysis,
and no allowance are consistent
is
made for corrosion. design criteria,
These assumptions
with the
B-9
.-.
TABLE 3.1:
Steel Mechanical
Properties
Property ~ Yield Stress (min.) (MPa) Tensile Stress (MPa) Elongation

YO
500 610-770 16
355 490-620 21
Youngs
Modulus
(MPa)
I
I
208700
I
I
208700
Poissons
Ratio
0,3
0.3
3.6
Interaction
with
Adjacent
Structure structural system and it above, for. Influence of unmodelled structure Para. 2.2.3
The midbody comprising longitudinal k reasonable
web frame girders.
is part of an integrated
the inner and outer shells, the transverse However,
frames
for the reasons discussed structure is accounted
to isolate a single web frame for analysis provided with adjacent
that the interaction
The primary interaction pattern of interest longitudinal configuration longitudinal this structure. M.lith reference locations: 1. 2. 3. 4. Centreline centreline (vertically); structure. structure
with adjacent A reasonable
structure
(for the load load transfer for this by the the stiffness of via
to this analysis) is through
approximation
is to account
for the support provided
by using springs representing
to Figure 3,1,
springs are required at the following
of Main Deck to account longitudinal
for the deck
girder (vertically); for the inboard side girder for the outboard side girder
On Main Deck to account On Main Deck to account (vertical and horizontal On side shell to account (horizontal);
components); for the upper stringer
B-10
5. 6. 7. 8.
On side shell to account Bottom structure
for the lower stringer at the for the girders (3 locations to account for the shell
top of the turn of the bilge (horizontal); to account structure - vertically); Centreline centreline of bottom girder (vertically); to account (vertically), as the and for the bottom
Bottom structure Iongitudinals
Spring constants member
for the above items have been calculated at the midspan of the longitudinal
inverse of the deflection being evaluated
(list above) due to a unit point load with a midbody If a fixed end web The ends of the longitudinal assumed as pinned. the stiffness member(s)
placed at each of the points of intersection frame along its length, have been conservatively condition structure in reality, Spring constants Section 3.7 B4.4 calculated would
had been assumed,
of the longitudinal resulting in a greater
have been overestimated
load transfer
from the midbody
web frame than would be the case
and used in the FE model are listed in Properties,
Beam Section
Loads of vessel displacement, associated power of the Taking the Loads Para. 2.6 Para. 3.4 Influence of Extent of Model Para. 2.2.1
The ice load2 is a function account a uniform spacing) positioned
vessel, the region of the ship, and the Arctic Class. of the various factors pressure and 2.85 total load applied to the web frame is 4435 of 1 metre width metre height. kN.
with ship parameters
This is applied as to a pressure of
(which equals the web frame the pressure patch is
This translates
1.556 MPa.
As required by the standard
such that 10+ZO its height is above the waterline. of in Figure 3,3.
The load applied is illustrated
The ice loads
are
adapted
from
Ref. 1. The structural design philosoph y of Hence design loads calculated from
this standard is based on plastic design.
this standard will, for a well designed structure, result in extensive yielding. For the purposes of this example FEA, which assumes linear elastic behaviour, the load applied has been arbitrarily halved to ensure the structure remains elastic.
B-1 1
3.8
Boundary
Conditions Boundary Conditions Para. 2.6 Para. 3.4
Symmetry longitudinal conditions longitudinal
is assumed
about a vertical Therefore,
plane through the symmetry boundary along the
axis of the ship.
are applied to all nodes along the outer (longitudinal) restraint restraint axis of the vessel, and rotational about the
edges of the plates. This provides translational other two axes, Symmetrical structure through boundary conditions
are applied to the bottom intersecting the vertical plane In addition, the bottom
and the deck structure the longitudinal
axis of the ship.
shell plating along the centre line is fixed in the vertical translation to avoid rigid body motion
4.0
4.1
FINITE
ELEMENT
MODEL
General
Information the finite element area, moment model. Units Para. 2.1.7
S1 units were used throughout Therefore, Modulus, respectively. The global coordinate Global X axis : Global Y axis : Global Z axis : system the units of length,
of inertia, Youngs
and pressure were mm, mmz, mm4, MPa, and MPa
for the problem is as follows:
Global axes system Para. 2.1.8
athwartship vertical parallel to ship CL Selection (SHELL63) of ANSYS and stiffeners was selected and Element Types Para. 3.1
4.2
Element
The elastic shell element used for modelling
the web frame,
from the bottom edge of the using 3-D of
stringer of the side shell structure main deck. longitudinal (COMBIN14). The SHELL63 element The stiffeners
at the top of the turn of the
bilge to the start of the sloped section on the outboard in other areas were modelled of ANSYS, elastic beam elements (BEAM44)
The stiffness
girders were modelled
using linear spring elements
is well suited for modelling
linear behaviour The
of flat or warped,
thin to moderately
thick, shell structures,
B-12
element
has six degrees
of freedom
at each node: translations about the element using a mixed
in x,
the nodal x, y, and z directions y, and z axes. directions. interpolation properties formed nodes. element, The stiffeners section elements The deformation The out-of-plane of tensorial (if required).
and rotations
shape is linear in the two in-plane The element shaped element is defined by material may be
motion is predicted
components, A triangular
four corner nodes, four thicknesses,
and the orthotropic
by defining the same node numbers for the third and fourth Pressure load may be applied as surface loads on the
in the deck and bottom structure BEAM44
of the mid-body beam with tension, also in the side
have been modelled (BEAM44). torsion,
using 3-D elastic offset is an uniaxial element The stiffeners
compression, structure (SHELL63). To simulate discussed elements. elements elements, translations
and bending capabilities. per node. were modelled
This element
has six degrees of freedom diaphragms
using shell elements
the overall stiffness in Section COMBIN14 2,4, were modelled elements
of the rest of the structure,
as
the connection
points of the frame to
other structure
with linear springs (COMBIN 14) are uniaxial tension-compression at each node: Two sets of and the other
with up to three degrees of freedom in the nodal x, y, and z directions. one for springs in the horizontal direction,
direction
for springs in the vertical 4.3 Mesh Design
were defined.
The response particularly structure the following 1. 2.
of the side shell structure
is of primary interest the frame in
Mesh Design Para. 3.2
in the vicinity areas:
of the loading. Thereforer
has been modelled
with a fine mesh of shell elements
side shell structure
between and
the turn of the bilge and the between the side shell girder,
side shell upper stringer; upper stringer The remainder
outer edge of the deck structure
and the deck angled outboard
of the frame has been modelled is reasonably
using a coarse in an
mesh of shell and beam elements. of this part of the structure
This ensures that the stiffness modelled
B-13
economical
manner. of beam and shell elements, in Figure 4.1. used for the
The mesh, consisting frame consistent stress
analysis is shown
The mesh design is from the finite element a high
with the results expected

with
model, that is, a fine mesh is provided in the regions where
grdierl~expected is
a coarse mesh provided
elsewhere. establishing adopted purpose. openings accuracy
The mesh is most dense around openings which are Since the primary interest is in the mesh density for this overall adequacy of the structure,
sources of stress concentrations. is designed
to yield stresses that are accurate
Based on preliminary of roughly A 5Y0.
analyses the mesh around these of peak stresses with an
should allow the prediction
The finite element and 18131 4.4
model contains
3758
elements,
3578
nodes,
total active degrees of freedom. Attributes and Spring Constants Stiffness and Mass Propetiies Para. 2.4
Finite Element
The attributes 4,1, properties
of the elements
used in the model are listed in Table based on the stiffness are listed in
The spring constants of the adjacent
calculated structure
Table 4.2.
prints a
To avoid ill conditioning warning than 1.0e08. Therefore, used.
in the stiffness
matrix ANSYS stiffness is 4179
if the ratio of largest to smallest The largest stiffness stiffness
value is greater matrix being N/mm. these of the
in the stiffness allowed
4.1 79e + 11, the smallest
springs with stiffness
less than 4179 low stiffness
N/mm were not values,
Because of their relatively
springs will have a negligible effect web frame.
on the overall behaviour
B-14
TABLE
4.1:
Finite Element
Attributes
Item No.
Description
Element Type &

No.
Mat. Type
& No.
EH36
Real Cons. No
Thickness or
122
Iyy
Area
mmlmmz
X106
mm4
XI03
mma
TKZTI
mm
TKYTI
mm
Diaphragms / Web Plating
Shel143
101
16
2
3
Floors - Web Plating
I
2
I ,, I
,, EH36 EH50
102
26 I 12
Deck Transverses - Web 1500xI
103
Deck Plating
Shel143
n
104
14
Outer Shell Plating
105
36
6 7
Bottom Shell Plating
,,
AH36
106
29
Deck Transverses - Flange Inner Deck Plating
Shal[43 .
EH36 ,,
107 108
19 14
Innar Shall Plating
,,
,,
108
16.5
10
11
Inner Shell Plating - Bilge
,,
!,
110
111
17
Tank Top Plating
,,
,,
13
12
Transverse Stiffeners - Diaphragms
Shel143 ,,
EI-136
112
16
13 74
Stringera Transverse Stiffeners - Tank Top
113
16
Beam44
AH36
114 115
5700 15
38.58
190.0
10
142.5
15
Girders - Tank Top
Shel143 Beam44
AH36
16
Deck Transverse Stiffeners
EH36 EH36
116
1575
2.95
14.47
5.25
75
17 18
Side Girdera Deck Plating (with openings)
Shel143
117 118 119
14 9.34 6576 92.56 140.3 205.5
Shel143 Beam44 Beam44 Combinl 4 Combin14
EH36 EH36 EH36 -
19 20
Beam Elements for stiffeners at Beam Elements for the bilge and
8 8
120
6676
92.e6
140.e4
205.5
21 22
Vertical Springs - to account for Horizontal Springs - to account for
see Table 4.2 for spring atiffneas
see Table 4.2 for spring stiffness
B-15
TABLE
4.2
Spring
Stiffness
Calculated
Based on Stiffness
of Adjacent
Structure
Description
Spring Direction Girder Vertical Vertical Vertical Horizontal Horizontal Horizontal Vertical Vertical
Element Type 5 5 5 6 6 6 5 5 5
Real Constant No. 121 122 123 124 125 126 127 128 129
Spring Stiffness N/mm 231 3785 3012 56 7151 7151 6508 5913 3631
Deck Centreline
Inboard Side Girder Outboard Outboard Side Girder Side Girder Stringer
Upper & Centre Lower Stringers Bottom Outboard Bottom Bottom Girder Girders Centre Girder -
Line
Vertical
4.5
FE Model
Loads and Boundary
Conditions Loads and Boundary Conditions Para. 2.6 in mm Para. 3.4
General information
on the applied load is provided in Section 3.7. MPa.
The design ice load was applied as a pressure of 1.556 The finite element Section 3.8. 4.1, in Section model boundary conditions
are as explained system of + 500 or -500 This
Referring to the global co-ordinate all nodes with Z - co-ordinate boundary conditions restraint boundary
described
have symmetry restraints translations Z axes.
along the Z axis.
provides translation line have symmetry
in the Z - axis, and rotational All nodes along the bottom conditions restrained centre along the X - axes, i.e., in the Y &
in the X and Y axes. restrained
in the X and rotations
The nodes along the bottom
centre line for the bottom For the top along boundary conditions
shell plating were also restrained the X - axis,
in the Y direction.
centre line, all nodes have symmetry
B-1 6
4.6
FE Model
Checks model was run, the following on the FE model : prerun Finite Element Model Checks Paraa 3.0
Before the finite element checks were performed . consistent coordinate element boundary The following user interface requested symbols boundary model. nodal coordinates element . . units system
attributes
and real constants and loads using the graphical provides a listing of entities. Also, such as
conditions
prerun checks were conducted provided by ANSYS. for specifically loads, element ANSYS selected
information can be turned conditions,
on/off to view various aspects, connectivity,
etc., of the
of extremities
of model discontinuities orientation skew,
free edge plots to check for structural shape; aspect ratio, taper, shrink plots and element connectivity checks for property property . . element type, etc. element assignment type,
edge plots to check element to elements material type, coordinate - using colour physical
coding based on element plot showing orientation
system to check
for element
true scale 3D plot of beam elements beam size, orientation, boundary condition pressure conditions symbols load magnitude and offsets
to ensure correct
- using model plots with boundary and direction (using arrows) and are or error The
The following performed messages
prerun checks are built into ANSYS, process.
during the data checking
Warning
are issued when the model fails to pass the check, for warning
output from such a data check run were reviewed and/or error messages. nodes not connected . elements missing material missing physical to structure to structure
not connected
properties properties
B-17
element element . element
aspect ratio warping skewness Option and Procedures used were: Solution Options and Procedures Para. 3.5
4.7
FE Solution
The following .
solution options and procedures
New Analysis Static Analysis No Stress Stiffening Small Deflections Store all results for all load steps Print all output to a listing file
5.0
5.1
ANALYSIS
RESULTS
General
Solution post-run
Checks checks were perfornied: to ensure that are included as General Solution Checks Para. 4.1
The following .
comparison Annex B-3)
with simple hand calculations
the results are reasonable equilibrium inspection inspection mesh used All error and warning investigated messages between
(these calculation
the applied load and the reactions shape of the structure in the model of the to ensure
of the displaced
that there were no discontinuities
of stress contours to ensure the adequacy
output by the program were
and resolved, is 4434.9 kN. No forces in
The total applied load in the X direction are applied in the Y and Z directions. the X, Y and Z directions respectively. 5.2 Post Processing graphical Methods post-processor are 4434.9
The summed
reactions
kN, O kN, and O kN
The ANSYS review
was extensively
used to to
Post-processing Methods
stress and displacement
results.
Listings were reviewed
B-18
obtain specific stress contour element extrapolating
magnitudes
for various quantities. was used.
In all of the For the shell by points.
Para. 4.2 Para. 4.3 Para. 4.4
plots nodal averaging
used in the model, the nodal values are calculated from the values at the integration
5.3
Structural
Response is shown in Figure 5.1, where FEA Results and Acceptance Criteria Para. 5.7
The deflected vertical mm. occurred
shape of the structure
the displacements The maximum
are scaled up by a factor of 20, horizontal displacement is 51,08
The maximum mm and small at 1.96
displacement
at the top centre line of the vessel is 124
on the inner shell in the vicinity of the load application. which was relatively between in the diaphragm two stiffeners This possibility the side shell and the This displacement a possible location for using classical
The out of plane displacement, mm, occurred opening, occurred methods,
also in the area of load application. between indicating
shear buckling.
should be checked
The Von Mises stress plot for the area of interest Figure 5.2, indicated yield (500 outer shell. The contours are arranged MPa) stresses MPa) past yield (355 in the outer shell.
is shown in
such that colour orange in all areas except the
Dark red shading is used to indicate stresses past It is clear from the figure that remains elastic, MPa. The outer shell MPa. except for
at the applied load the overall structure yield. The maximum stress recorded
a small area around the openings where the stresses are pasthere is 573
Figure 5.3 shows is in compression MPa.
contours
of bending stress, Sy. compressive
with a maximum
stress of 307
The inner shell has a maximum
tensile stress of 330
High bending stresses, around openings.
past yield stress, were again observed
Clearly the bending stresses in the outer and the
inner shells are below the yield stress. A contour Figure 5.4. 205 plot of shear stresses The maximum in the diaphragm is shown in were Figure
and minimum
stresses recorded
188 and 164 MPa respectively. MPa, the structure an enlarged 5.5 contains concentrations
The yield stress in shear being
remains elastic at the applied load. view of shear stresses around the under the load. The stress
opening which is directly
around the opening are clearly visible in the figure.
B-19
The smoothness probably adequate
of the contours
suggests that the mesh density
is
for the purposes of this study.
6.0
CONCLUSIONS framing section of the 50000 meets the acceptance DWT tanker as criteria. At the towards two at elastic except Overall Assessment Para. 5.5
The midbody designed applied load,
and analyzed
the structure displacement
remains predominantly in the diaphragm, investigation.
in a very localized an out-of-plane stiffeners higher loads.
region around openings.
The tendency between
in the area of an opening, This needs further
could result in instability
7.0
1.
REFERENCES PROPOSED Safety, March EQUIVALENT Canadian STANDARDS FOR THE Dated
CONSTRUCTION (AMNB) 1993.
OF ARCTIC
CLASS SHIPS; Arctic Ship
Coast Guard - Northern;
2.
Arctic Tanker SECTIONS, Shipping
Structural
Requirement
Evaluation
MIDSHIP AA
BOW SECTIONS Limited;
AND REPAIR DRAWINGS; Dated June 1994,
Company REGISTER 1993
3.
LLOYDS TESTING January
- RULES FOR THE MANUFACTURE, OF MATERIALS; Dated
AND CERTIFICATION
B-20
50,000
CDWT
Midship
Fromlno)
Section
(5D0 MPa Sh.lt
PIotino with 355 UP.
:CK PLAllNG 14mm GR. EH36 DEI :CK LONGLS 200.(2 F.B. I&Q. EH36 DE, 5PACEO 750mm MAK. :CK tRANSVERSES 1500*t2mm !4EB/S50.l 9mm F.F. OEI E GRAOI EH36 5FACED 3002mm OECK TRANSVERSE STIFFENERS 150.10.5mm F.B. GRAOE F.F. CENSEE :EO II+US
I
J
nFCK PLATING 14mm GR. FHS6
1 t o
i 0
1 1 1 I
pl t
!
INNER SKIN L-% SPACED 720 mm 7RAM1710N 2SO*14 mm GE!. EH36 i i SHEU PLATE 36 mm ~AOE ELSE SAME AS BELOW FH5J3
St+ELL PLATE 36 mm GRAOE EH50 DIAPHRAGMS SPACEO 1000 mm 16 mm PLATE GRAOE EUX 000.600 CUTOUTS STIFFENERS SPACSO 7S6mm SnFFENERS 411 q16mm F.E. GR. H436 SNIPEO Al ENOS STRINGERS %IACEO 5500 mm 16 mm PLAIE GRADE EH36 mlfioo CUTC421S snFFENERS SPACEO =Eunm snFFEtJERS 261.16mm F.B. cR. HU6
Iml
ll
[
1
1S.5 mm PLATE CR. 0H36
K
iv NOTE:
! !
TANK TOP PLATING 13mm GiS. AtlS6 TMK TOP LO+JGLS 325.19mm F.B. CR. AHS6 5PACE0 750 mm j
IHIHIHIH ::; :;1:::::[ ::l:;;

i 90TTOM SHELL 29 mm GRADE AH36 BOTTOM LGUGS SPACEO 750 mm 562$29 mm F.B. GRAOE AH36 FLAT BAR STIFFENERS SPACEO 750 mm STIFFENERS 100020 mm GRAOE A+t3B ~LOORS SPACEO 3000 mm X mm PLATE CRAOE AJ+36 T.B. STIFFENERS SPACEO 750 mm SnFFEWERS 2U5.20 mm GRAOE AH36 EOO*600 CUTOUTS
slRtNGERS SPACEO 5500 mm 16 mm PLATE CR. 0H36 8cQ*600 CUTOUTS sTIFFENERS SPACEO 500mm STIFFENERS 261.16mm F.B. G% 0H36
CJROERSSPACEO 4500 mm 15 mm PLATE GRAOE AH36 STIFFENERS S?ACEO 750 mm STIFFENERS 433.15 mm GRAOE Ak+= BOO.600 CUTOUTS
AU OiMEtASIONSARE IN MILUMETRES.
_t-er-:7
....--
........-
u-L
l---+
w + o z
FIGURE
3.2
Outer
Dimensions
of Web
Frame
B-22
.-.
. .
.,.. . .
.,.. .
I I
3364
r
I Lctad Footprint 2850 X 1000 Pressura 1.556 MPa
F
lE 00 15 34
2000 I ) t
.
9502
FIGURE
3.3
Characteristics
of Load
FIGURE
4.1
Finite Element
Model
of Web
Frame
B-24
/,
-, . . . .. .
z-$
FIGURE 5.1 Deflected Shape of Web Frame
B-25
B-26
.,, %,.
B-28
., .....\
.,
18, ,,
,.,
., ,,
, ,
,.
..
B-32
L.. ,.,,,
. .
,.
(L.,...
Annex B-2 Company and Personnel Qualifications
B-2
COMPANY AND PERSONNEL Qualifications
QUALIFICATIONS
B-2. 1 Contractor
66 Engineering also certified its staff. BBE performs system, finite element element adequate,
(BBE) is an ISO 9001
compliant
company
with a firm commitment structural engineers
to quality.
[t is on
by the Association
of Professional
Engineers of Ontario.
BBEs primary
business is Ship
Design and Analysis.
It has several qualified
professional
and naval architects
all its finite element 486 PC. software
analysis on either a DecStation ANSYS
5000,
running on Ultrix operating 5000 was used. The software finite finite element
or on a 60 MHz,
For the current analysis the DecStation
used is called ANSYS. [t has been successfully
is a well established
with a large user base. analyses.
used by BBE in several of its ship structure
ANSYS
provides all the required features
for the current task and hence deemed
B-2.2 Analyst
Personnel
Qualifications
Mr. J. S. is the finite element and is registered finite element method
analyst assigned to this task,
He has a Ph.D. in Structural
Engineering,
as a Professional
Engineer in the province of Ontario.
He has taken two courses in in using finite element out of which that analysis problems
analysis at the graduate specific,
level, and has eight years experience Information on request. on specific finite element
as an analysis tool.
JS has a total of five years experience
in using ANSYS,
three years are ship structure JS has worked Checker
on in the past is available
Ms. J, B, is the project engineer element element analysis, analysis, Professional experience has worked
for this project,
and holds the responsibility Engineering, She has taken one graduate in finite element analysis. and has supervised
of checking
the finite as a
JB has a Masters
Degree in Structural
and is registered
Engineer in the province of Ontario. and has six years experience JB has three years experience on request.
level course in finite several finite element on projects that JB
JB has gained ten years
in the design and analysis of ship structures, in the past is available
analysis projects.
in using ANSYS.
Information
B-35
.. -,
Annex B-3 FEA Results Verification
B-3
FEA RESULTS
VERIFICATION with hand calculations, as follows: Two Accuracy Assessment Para. 5.4
The FEA results were compared analyses 1. have been performed
An elastic beam analysis of the frame with a span of 11000 millimetres, ends fixed, openings (9.373*0.8*0.5*0,83), ignored, subjected millimetres for a total to a uniformly equal to 3.112 load of 8869 The structure support distributed MN/m kN, has a bending stress of 550 MPa at the top in the load of length 2850
in the inner hull plating.
Shear stresses
portion of structure This structure approximately 2.
above the load are 195 MPa. first yield (in bending) at a load of
reached 5700
kN. was FE modelled, was
An elastic frame analysis of the structure except that the inner shell and bottom analyzed thickness centreline side sway moments
structure
with a flange width equal to 40 times the plate and the frame was assumed to be fixed on at the deck and at the bottom, of the frame was ignored. calculated were within In this analysis The bending
a few percent of those
found in the first analysis. By comparison shell framing comparison consistent analyses. the FEA predicts first yield, of the inner hull mm portion of the side 4835 kN, This simplified at a load of approximately suggests
plating at the top of the 11000
that the FEA results are broadly
with the results from the approximate
B-36
Annex B-4 Sample Completed Assessment Methodology Forms
EVALUATION
OF FINITE
ELEMENT
MODELS
AND
RESULTS
Project
#:
Xxxx
Project
Title:
Finite Element
Analvsis
o f Arctic
Tamker
Web Frame
Project Description: Linear. stat ic analvsis o f web frame t~ ensure adeauacv o f frame
tce load
Contractor:
BB Enaineerina
Ltd.
Result
of Generallv satisfactory. Final a~wo val subiect to the sumlv of data
Evaluation:
on some d etails of the model
Evaluator:
John Doe
Date:
Mav
7995
B-37
i{
L,-...
1- Prellmlnarv Checks 1.1 Documentation

Perform these checks to
1 R9sult
eneurathat the analysisdocumentation.ob j specification,FEA software,and wnkactor 1analystqualification requirementshave been addressed.
1.2 JobSpacifititicm
1.3 Fin Its Element Analysis

1.4
Sofwara
Contractor/ Analyst Qualifications
Ye*~
Result
R
Yes / / / Result / / / / /
Preliminarychecks are acceptable? No
LN.
2. Engineering Model Chacks 2.1 Analysis Type & Assumptions Perlormthese checksto ensurethat the assumptionsused to developthe engineetirtgmodel of the problemare reasonable. 2.2 Geometry 2.3 Material Propemes 2.4 Stiffness & Mass Properties 2.5 Dynamic Degrees of Freedom 2.6 Loads & Boundary Conditions
Engineeringmodel is acceptable?
Yes~
f40
3. Finite Element Model ChOcka 3.1 Element Typaa PetiiTn these checks to ensure that the flnlteelement model is an adequate interpretation the engineeringmodel. of 9.2 Meah Design 9,3 Subetructums q Submodels nd 3.4 FE Loads& Boundary Conditions 3.5 FE Solutton Options & Pmceduras
Finiteelement model is a~ptable ?
&
4- Flnits Element Analysis Results Checks 4.1 Gonerel Solution Chaaks Result
Periorrn mesechecks toenaura that thefinite element esults r are calculated. Promssad andprasenlnd in a manner onsistent c withtheanalysis requirements.
/ 5
1[
4.2 PostPmwaalng Methods 4.3 DIaplacemant Results 4.4 StraasRaaults 4.SOtherResults
Finite element results are
r
v
--
--
4
Performmese checksto ensurethat adequate mnsideraUonof the loads. strength,accaptanaeciiterla, FE model,and rasultsaccumy are includedin amivingal the conclusions from the finiteelement analysis.
yes~
-No
5- Concluelona Checks 5.1 FE Results& Acceptsncs Criteria 6.2 Loada Aa=aaamont &3 Strength/ Resistance Assaasmnnt 5.4 Accuracy Assessment 5.5 Overall Assessment
Rasuit / / / / / Conclusions of Itw analysisare
Yes-
Lt.Jo
B-38
,.,, ..--
.,.
FINITE ELEMENT ANALYSIS Project No. Contractor Analyst : XXXX Name: JS
ASSESSMENT
I PRELIMINARY
1
CHECKS
Project Title : FEA of Arctic Tanker Web Frame BB Engineering Ltd Checker : Date : JB hlay 1995
1.1
Documentation
Requirements Refer to
Guideline Section
Result
Comments
1.1.1
Has the following information
been
3-1.1
provided in the FEA documentation? a) b) c) d) e) f) g) h) i) j) k) 1) Objectives and scope of the analysis. and acceptance criteria. d d d # d d d d d d V d # # d d d # ~ K Some detail missing *
Description of physical problem. Description of engineering model. Type of analysis, System of units. Coordinate axis systems.
Description of FEA model, Plots of full FEA model and local details. Element types and degrees of freedom ~er node. Material properties,
m) Element properties (stiffness & mass properties). n) o) P) q) r) FE loads and boundary conditions, Description and presentation Assessment of the FEA results,
of accuracy of the FEA results.
Conclusions of the analysis. List of references.
Based on the above checks answer Question 1.1 and enter result in Fiaure 1.0. 1.1 Is the level of documentation sufficient to perform an assessment
of the FEA? I
Comments
*Request
additional
detail on stiffener/web
connection
B-39
1.2
Job Specification
Requirements
Finite Element 1.2.1
Assessment
Check and
Refer To Guideline Section z
Result
Comments
Is the job specification referenced
identified
in the analysis documentation?
1.2.2
Are the objectives clearly stated
and scope of the analysis with
3-1.2
and are they consistent
those of the job specification? 1.2.3 Are the analysis requirements and are they consistent job specification? 1.2.4 If certain as certain justification 1.2.5 requirements load cases), of the job (such has adequate -3-1.2 N/A clearly stated
3-1,2
d
with those of the
specification
have not been addressed been given?
Are the design / acceptance stated the job specification?
criteria clearly with those -of
3-1.2
and are they consistent
1.2.6
Is there reasonable
justification
for using
3-1.2
FEA for this problem?
1.2.7
Has advantage experimental, that are relevant
been taken of any previous analytical, or numerical works to this moblem?
3-1,2 N/A
Based on the above checks answer Question 1.2 and enter result in Figure 1.0. 1.2 Does the analysis address the job specification requirements?
I I
Result d
Comments
B-40
1.3
Finite Element
Analysis
Software
Requirements
Finite Element Analysis
Assessment
Check
Result
Comments
1.3.1
Is the FEA software applications?
on the list of approved
programs for ship structural analysis
If the answer to Check 1.3,1 is Y, you may skip Checks 1.3.2 and 1.3.3. 1.3.2 Are the capabilities and limitations of the FEA software used to perform the required analysis stated in the analysis documentation? 1.3.3 Is evidence of this capability documented available for review (eg, verification and 3-1.3 # 3-1.4 d
manual,
results of ship structure FEA benchmark tests, previous approved FEA of similar problems)? 1.3.4 Does the vendor of the FEA software standards are maintained in software development and maintenance. have a
quality system to ensure that appropriate
Based on the above checks answer Question 1.3 and enter result in F[qure 1.0. 1.3 Is the FEA software qualified to perform the required analysis?
Comments
B-41
1.4
Contractor
/ Personnel
Qualification
Requirements
Refer To Finite Element Assessment Check Guideline Section 1.4.1 Do the contractor personnel have adequate 3-1.5 Result Comments
academic training and experience qualifications to perform finite element analysis? 1.4.2 Do the contractor personnel have adequate qualifications for 3-1.5
engineering experience
performing ship structural design or analysis? 1.4.3 Do the contractor have adequate qualifications? 1.4.4 Does the contractor Quality Assurance that are satisfactory 1.4.5 Do the contractor experience analysis? have a working system of (QA) procedures and checks for the requirement? 3-1,5 3-1.5 Not documented but and contractor personnel 3-1.5
professional certification
using well established software
personnel have adequate used for the
with the FEA software
Based on the above checks answer ) 1.4 Is the contractor adequately
Question
1.4 and enter result in Figure 1.0. ship structure FEA?
I Id
Result
qualified for performing
Comments
B-42
,, .
FINITE ELEMENT ANALYSIS Project No. Contractor Analyst : XXXX Name: JS
ASSESSMENT
I ENGINEERING MODEL CHECKS
Project Title : F&4 of Arctic Tanker Web Frame 66 EngineeringLtd I Checker: Date : J/3
May 1995
2.1
Analysis
Type
and Assumptions
Refer To
Finite Element Analysis Assessment Check 2.1.1 Does the engineering model employ enough dimensions and freedoms to describe the structural behaviour (eg. 1-D, 2-D, or 3-D)? 2,1.2 Does the engineering model address the appropriate scale of response for the problem (eg. global, intermediate, or local response)? 2.1.3 Is the type of analysis appropriate for the type of response and loading of interest (eg. linear, static, dynamic, buckling analysis)?
2.1.4 Does the engineering displacement, 2.1.5 model address all the (eg. stress, buckling load)? the choice of
Guideline Section
Result
Comments
3-2,1
3-2.1
3-2.1
3-2.1
required results parameters frequency,
Are all assumptions engineering (watch
affecting
3-2.1
model and analysis type justified assumptions)? or conservatism for the of 3-2.1 Appears marginal results to complete evaluation - may
for non-standard
2.1.6
Is the level of detail, accuracy the engineering criticality
model appropriate
require more data on
of the analysis and type of problem?
2.1.7
Does the analysis employ a consistent units?
set of
3-2.1
2.1.8
Does the analysis employ a consistent global coordinate axis system?
3-2.1
Basedon the above checksanswerQuestion2.1 and enterresultin Figure 1.0. 2,1

Are the assumptions of the type of analysis and engineering model acceptable?
Result
Comments See above
B-43
B-44
,,.
Appendix
Examples of Variations in FEA Modelling Practices and Results
EusnQ!E
cl
C2
IM
Stiffened Panel Multiple Deck Openings Mast
Paw
c-3 C-17 C-25
C3
c-1
INTRODUCTION The purpose of this Appendix parameters Three typical is to illustrate the effect of varying certain FEA modelling example problems. presented in Section Cl,
on the results using typical ship structure ship structure examples are used.
The first example,
concerns the modelling of stiffened panels. Four different approaches for modelling stiffened panels are considered and the results presented. In the second example, presented in Section C2, the modelling of stress concentrations arising from openings in a deck structure is considered. In the third example, presented in Section C3, variations in the approach to modelling a truss type mast structure are illustrated, A brief introduction is provided for each problem, followed by a pictorial overview of the FEA model and results, A brief discussion of the results is provided at the end of each example. It is not the intention of this Appendix to endorse any particular modelling method, Rather, it represents an effort to illustrate various modelling practices and present the variations in results. This should provide some insight into the consequences of adopting a particular modelling approach. The choice of the appropriate method, for a given problem, depends on the purpose and objectives of the FEA. In all cases the ANSYS . .
q q
program was used,
The following
element
types were used:
four-node four-node
eight-node
membrane shell elements shell elements with bending capabilities shell elements with bending capabilities
. .
two-node 3-D beam elements two-node 3-D truss elements mass elements
In certain cases converged solutions are referred to. These solutions result from very fine mesh models which are known to have converged (by comparison with less fine mesh models).
c-2
C1.O
STIFFENED
PANEL weight in conventional ship structures is stiffened panels that
The majority
of the structural
comprise the shell, decks, bulkheads and superstructure. The panels are stiffened with structural sections that are usually spaced in a regular fashion. The appropriate modelling approach for stiffened panels depends on both the scale of the response (ie, local or global response) and the main structural actions of interest. Two main structural actions typically modelled are 1 ) bending action due to loading normal to the panel surface, and 2) membrane action due to loading in the plane of the panel, The first part of this section deals with bending action and hence focusses on stiffened plate subjected to transverse loading. Membrane action in a stiffened plate as a result of in-plane loads is briefly examined in the second part.
c-3
,,.;
FEA Example No. Problem Description:
Title :
Stiffened
Panel - Transverse Loading
There are various techniques
available for modelling
stiffened
panels.
The choice of a
particular technique depends on the purpose of the analysis. Using a simple stiffened panel structure, the differences in the accuracy of stress and deflection results for some of these techniques are examined. Engineering Model :
t,,oo-j
Stiffeners:
FB 150x
10.5
T 3000
Plate: t=l
Omm 4
Material Properties : E = 207xIOS v = 0.3 MPa
Geometric Properties : Plate Stiffeners t=lomm 150 x 10.5 FB
Loading : PZ = 15000 Pa
Four modelling approaches are considered: Modelling Features : 1. Modelling stiffeners with off-set beams (beam properties defined at beam centroid which is rigidly off-set from plane of plate); 2. Modelling stiffeners with in-plane beams (beam properties includes an effective width plating and are defined at beam centroid which is in the plane of the plate); 3. Explicit modelling of stiffeners using shell elements; and 4. Modelling the plate with orthotropic material properties (in-plane loads / membrane action only)
of
c-4
FEA Example No.

Finite Element Models
1 I :
Title : Stiffened
Panel
- Transverse Loading
A total of 12 FE models, grouped into four sets, were studied. Each set contained three models representing the three modelling techniques. The mesh and element types are as follows : 4x4 element mesh; 4 noded elements Set 1 Set 2 8x8 element mesh; 4 noded elements Set 3 16x1 6 element mesh; 4 noded elements Set 4 16x1 6 element mesh; 8 noded elements All models are fully fixed along the four edges. kN/m2 is applied. For the in-plane
IS
A uniform transverse
pressure
load of 15
beam models the effective
width of plating was assumed to be 40t,
where
the thickness of the plate. The inertia propetiies of the beam were calculated based on stiffener and an effective width of plating. However, for the area, the area of the stiffener alone was input.
Elements
QE91!M of freedom
25
28
150
Example 1 a - Offset Beams
81
486
Example 1 b - Offset Beams
c-5
FEA Example
No.
Title :
Stiffened
Finite Element
Models :
M!L!Es
Elements
Dearees of freedom
289
304
1734
Example 1c - Offset Beams
833
352
4230
Example 1d - Offset Beams
28
150
Examtde 1e - In-rdane Beams
81
486
Example 1f - In-plane Beams
C-6
... ~. ..
.,
FEA Example No.

Finite Element Models
1
:
Title :
Stiffened
N!2r!Es
Elements
Dearees of freedom,
289
304
1734
Example 1g - In-plane Beams
833
352
4230
=xample
1 h - In-plane
Beams
40
28
240
=xample 1 i - All plate elements
648
Sxample lj - All plate elements
c-7
FEA Example No. Finite Element
I
m
Title :
Stiffened
Models : b!@2S EhlM!tS Dearees of freedom
391
352
2346
Example 1 k - All plate elements
1133
352
5886
Example 1 I - All plate elements
C-8
FEA Example No. DISCUSSION
Title :
Stiffened Panel - Transverse Loading
OF RESULTS
Key results are summarized in Table Cl, 1, The maximum vertical deflection is at the centre of the panel (see Figure Cl. 1). The peak stresses reported in the table are at the ends of the central stiffener (at supports) , The three mode shapes associated with the three frequencies are shown in Figure Cl .2. Figure Cl.3 shows the longitudinal stress contours for the plate and the stiffeners. Figure Cl,4 summarizes the deflection results for all twelve models. From Figure Cl,4 it is evident that the deflection solution starts to converge for an 8x8 mesh. Figure Cl,4 also shows the stress results in the stiffener. Some general observations for the three modelling types are : In-Plane Beams: Despite the approximation of 40t as the effective width of plating this method seems to provide the most economical solution for deflection prediction. The same is true even for stress prediction. Deflection decreases with mesh refinement contrary to the expectation that displacement-based FEA model becomes more flexible with more elements. This is probably due to the presence of a spurious moment generated at the ends of the stiffener as a result of two axial forces {in the plate and in the Howeverr with mesh refinement this effect tends to beam) being offset, diminish resulting in reasonable predictions of deflections. In this case the performance an 8x8 mesh, approaches that of the in-plane beam models with
Offset Beams:
All PIElem ents:
All three techniques
predict natural frequencies and mode shapes fairly well.
In modelling stiffeners as in-plane beams, the greatest uncertainty is the choice for the effective breadth of plating. The most important parameter which determines effective breadth of plating is the ratio of actual flange width to the length between points of zero bending moment. The effective breadth of plating can be estimated from charts (see, for example, Hughesl). Another important aspect to note with this technique is that the effective breadth thus used is only effective at the location of maximum ,bending moment. However, for design purposes the stresses at the section of maximum bending moment is of most importance. In conclusion, the approach recommended will depend on the nature of the analysis, [f the plate-stiffener combination is subjected to transverse loading, modelling stiffeners with in-plane beams provides the most economical approach in terms of overall stiffness, and stresses in the stiffener at the location of maximum bending moment. When more detailed stress information is required then the explicit modelling of the stiffener with plate elements appears most appropriate, The use of the offset beam is attractive since there is no approximation required for effective breadth, With a reasonable mesh density (at least 3 elements between stiffeners) this technique should provide reasonable prediction of the overall stiffness of the structure.
1 Owen F. Hughes, Ship Structural Design - A Rationally-Based, Optimization Approach, John Wiley & Sons, New York, 1983.
Computer-Aided,
c-9
L<,.
.. ,.
TABLE Cl. 1 Stiffened Modelling of stiffener SETI: 4x4 Mesh (mm)
Panel FEA - Results In-plane beams Ie 9,51 5.95 45.20 -246,40 45,20 Plate elements Ii 4,48 16.11 -98.31 5,59
Offset beams la
Max, Vertical Deflection
Max. stress in plate (MPa) Max. bending stress in stiffener at ends (MPa)
32,87 -379,90 289,30
24,94
First three natural frequencies (Hz) SET2: 8x8 Mesh (mm) 29.12 38.34 lb If
30.89
34.00 43.54 Ij
30.02
33.93 35.24
Max. Vertical Deflection
7.70
33.87 -339.20 181.80 28.11
6,86
47.69 -259.95 47.69 29,71 32.40 43.96 lg lk
6.64
24,12 -175,58 15,81 30,50 33.93 45.60
Max, stress in plate (MPa) Max. bending stress in stiffener at ends (MPa)
First three natural frequencies (Hz) SET3: 16x16 Mesh (mm) lC
31.89 43.33
6.90
38.96 -307.50
6,69
48.22 -262.88
I
1
6.80
33.15 -226.17
I
Max, stress in plate (MPa) Max. bending stress in stiffener at ends (MP~)
I
1
112.98 29.59
48.22 29.87 32.60 44.64

Ih
26.02 29.84 33.51 45,55

11
First three natural frequencies (Hz)

SET 4: 16 x 16 Mesh (8 node)
33.31 45,29 Id
6.70
(mm)
6,65
6.88
Max. stress in plate (MPa) Max. bending stress in stiffener at ends (MPa)
1
47.26 -289,67 75,37

30.02
I
48.47 -264.25 48.47

29.94
1
50.55 -287.29 41.42

29,58
First three natural frequencies (Hz)
33,73 45.95
32,70 44.93
33.35 45.53
c-lo
FEA Example No. 1
Title :
Stiffened
Figure Cl.1
Deflected Shape
c-1 1
FEA Example No. 1
Title :
Stiffened Panel - Transverse Loading
? / / /
1 /
/ I
/
1
/ I
Mode 1
Mode 2
Mode 3
Figure Cl.2
Mode Shapes
c-1 2
C-14
,.
FEA Example No. 1
Title :
Stiffened
...
..
10
4 4X4 8X8 16X16 Mesh Density
J 8 Noded
I Maximum
400
Stress
in Stiffener
+-
E!El
8x8
0
4X4
16X16 Mesh Density
8 Noded
Figure Cl.4
Summary of Deflection and Stress Results
C-15
..
FEA Example No. 1 I In-Plane Loading :
Title :
Stiffened
Panel - In-Plane Loading
The second part to this example considers the same stiffened panel subjected to in-plane loading. The problem was modelled in two ways : 4; Using ordinary membrane elements but with orthotropic material properties; and Explicit modelling of stiffeners using 4 node membrane elements as per Example lj. :
Description
To model membrane action of stiffened plate structure advantage can be taken ,of the facility, available in most general purpose FEA packages, to model material orthotropy. Using an approach presented below (adapted from Hughes, see Reference on page C-g)[ it is Possible to simulate structural orthotropy by material orthotropy. The appropriate expressions are: A
1
EX=r
E
.--
.-
.._
____
___
____
EY = r E / [r - v2(r-1)1 Gxy=G=E/[2(l+v)l
____
_____
_______
-_
_____
_____
____
. -
--
____
-_.
____
__
.-
___
_____
____
. R
. -AREA
-A, s . 2 * *I s ~x
OF STIFFENER
The value of r is defined in the figure above, With this approach the stiffened plate structure is modelled using ordinary membrane elements but with orthotropic material properties. The expressions given above assume that the stiffeners are aligned in the x direction. The expressions can be altered to reflect stiffener alignment in the y direction. Care must be taken to ensure that the local coordinate system for the element corresponds with that assumed for defining the material properties, A further assumption implicit in the approach is that the stiffeners are assumed to have identical properties and to be equally spaced.
Results : Table Cl.2 presents the results for the two cases investigated under in-plane loading. The case with orthotropic material properties predicts plate stresses and displacement reasonably accurately. It is important to bear in mind that the plate stresses obtained directly from the FEA for the orthotropic plate are incorrect, However, the actual stress can be derived from the predicted stress by factoring it by 1/r, TABLE Cl.2 Comparison of Finite Element Model Results Orthotropic material proper-ties 346.00 Stiffeners modelled explicitly with plate elements
Description Stress in plate (MPa)
350.00 -1.51
7.52 -0.08
Ux
Displacements u, LIZ *
-1.50
7.51 0.00
Obtained by dividing the predicted FEA stress by the factor r
C-1 6 ... ...
..<-
C2.O
Multiple
Deck
Openings
A deck with multiple openings is used as an example to illustrate the influence of mesh density and the element type on deflection and stress results. The mesh density is gradually increased from coarse to fine, Two types of elements, 4-node membrane elements and and 8-node shell elements, were used. The example also illustrates the effect of varying element aspect ratio. The results obtained solution. from the various trials are tabulated and compared with the converged
c-1 7
?EA Example
No.
Title : Multiple
Deck Openings
Problsm Description: & deck with multiple openings is used to illustrate the influence of mesh density, element aspect ratio, and type of element on deflection and stress results. The density of the mesh is qradually increased from coarse to fine. The use of two types of elements, four node linear and In addition, dummy line elements with very small sight node quadratic shells, are illustrated. wea are used along the edge of the opening to extract maximum principle stresses. The latter may be used to overcome errors resulting from extrapolation of stresses from the shell element ntegration points to the nodes along the edge of the opening. Engineering Model :
~~ ~
750
50
MPa
T ~
A :{ 750 -a 600 450 x 45a
50 MPa
50 R
%
300R
b 1350
+
750 ~ C.L.
x+
1 ~o~++
ship Axis
Material
Properties MPa
Geometric Deck Plate Long. Stiff,
Properties
Loading : Uniform M Pa Tension =50
E = 207x103 v = 0,3
t = 6.35 mm 152x 102 Tee

127x 102 Tee mm
Trans. Stiff. Major Access Coaming
50 x 6.35
BC on +/Symmetry Boundaries
FB
Modelling
q q q q q
Features
modelling around stress concentrations selection of element type effect of varying the mesh density use of higher order elements effect of aspect ratio in the area of stress concentrations
C-18 ...
,.
..
FEA Ex;mple No. Finite Element Models :
Title
: Multiple
Deck Openings
Nodes
Qe~r=s Qf freedom
2a : 4-noded membrane shall elements 2e : 8-noded shell elements
214 995
235 465
642 5970
2b : 4-noded membrane shell elements 2f : 8-noded shell alements

lllllllr T1lll[ I I 1 i 1 I 1 I I 1
351
3044
379 1256
1053 18264
1 I
I I I I I I 1
2C : 4-noded membrane shell elements 2g : 8-noded shell elements
1213 4842
1104 1924
3639 29052
2d : 4-noded membrane shell elements 2h : 8-noded shell elements
3186 9368
3272 3540
9558 56208
c-1 9
FEA Ex;mple No. DISCUSSION OF RESULTS
Title: Multiple
Deck
Openings
The analyses revealed peak stresses at the lower left corner of the smaller opening as shown in Figure C2, 1 (the top figure shows stress contours for the full model and the bottom figure provides a close-up view of stress contours around the smaller opening). The stress concentration near the larger opening was relatively insignificant due to the presence of the coaming. When the mesh density around the openings was increased, with the aspect ratio held constant, the results indicate a progressive increase in the magnitude of peak stress. The results listed in Table C2. 1 indicate a converging trend in the magnitude of peak stress with mesh refinement. Although the peak stress always occurs at the same corner, it should be noted that the precise location of the peak stress varies slightly with the refinement of the mesh (number of nodes around the corner radius). Some of the differences in the results may also be due to different mesh transitioning (from areas of coarse mesh density away from the openings to areas of high mesh density at the openings) in the different models. The results in Table C2. 1 indicates the rate of convergencence of the stress results is greater for the line elements (truss or spar elements with only one degree of freedom per node placed along the edge of the openings) than it is for the plate elements. The use of line elements for obtaining stresses also overcomes stress extrapolation errors that arise in shell elements. Note that the stress results for shell elements must be extrapolated from the element integration points to the node locations at the edge of the opening. Parametric studies were conducted to evaluate the effect of aspect ratio in predicting stress concentrations. The mesh density of Example 2d was used as the basis for this investigation, The aspect ratio of elements around the smaller opening was varied from 1,05 to 3.00. The results, Table C2,2, indicate that the best values for stress concentrations are obtained when the aspect ratio is close to one. The difference in the stress results when the aspect ratio is changed from 1,05 to 3.00 is about 8Y0.
C-20
,.. . .,\._ 1 .
,..,
b
M
.r
TABLE C2. 1 FE Results of Mesh Density
Parametric
Studies Peak Stress
No.
Description
pa:;:
Max. Disp. (mm)
Shell Elem. (MPa)
Line Elem. (Mpa)
2a
four noded one element around the radius four noded elements around the radius
1.29
1.8
300
399
2b
two
1.38
1,8
369
453
2C
four noded four elements around the radius four noded elements around the radius
1.37
1.8
502
556
2d
eight
1,37
1.9
572
593
2e
eight noded one element around the radius -eight noded elements around the radius
1.38
1,9
543
557
2f
two
1.37
1,9
570
606
2g
-eight noded four elemr~d~saround
the
1.36
1.9
583
607
2h
eight
-eight noded elements around the radius
1.37
1.9
591
609
Aspect
ratio of elements
near stress concentration
(see figure on following
page)
C-23
...,., .
ELEMENT ASPECT RATIO = a / b
TABLE C2.2
Results from Aspect
Ratio Parametric Peak Stress in Plate Elem. MPa
Studies
Trial No.
Aspect
Ratio*
Relative ** Peak Stress Ratio
3,00 1.98 1.37 1.05

near stress concentration
537 561 572 585
0.92 0,96 0.98 1.00
2 3 4
* Aspect ratio of elements
* * Ratio of peak stress to that for trial No. 4 (plate element
aspect ratio of 1.05,
i.e. 585
MPa)
C-24
C3.O
MAST
A major factor in modelling of lattice masts is the modelling of the connection details, Depending on the type of connection, the joints can be modelled with fully rigidity at the joint, or some or all members can be modeiled as pinned (hinged) joints. A simple truss-type mast structure is used to illustrate both these options. In the case of rigid jointed structure, the mesh density (i. e., the number of elements per member of the mast) was varied to investigate the influence on the results. Both static and dynamic analyses were performed on all these models.
C-25
FEA Example Problem
No.
Title
: Mast
Description:
The truss-type mast structure shown below, consisting of steel pipe sections, is to be analyzed for shock accelerations loading and to calculate frequencies and mode shapes. Engineering Model :
01 Deck Level
1 Deck Level
Material ;=
J =
Properties 03 MPa
Geometric see Table
Properties C3.1
207x1
0.3
Features :
Loading : Base Accelerations: 8g inX 18g inY 8g inZ
Modelling
q q q
pinned and rigid connections model refinement static and dynamic analyses
C-26
./,-.
I
Finite Element Models : The finite element
Title : Mast
models of the mast are as shown
below.
However, if the member is continuous and Example 3a is modelled with all joints pinned. has nodes between the two ends (viz. two or more elements per member) then rotations are restrained at such nodes to simulate the continuity of the member. The following is a list of members that are treated continuous: - Main legs - Horizontal members - One out of the two cross braces at every level - Principal members of the spur frame Examples 3b and 3C are modelled with all rigid joints.
The three-dimensional beam element (BEAM44) of ANSYS is used in modelling mast members, This element has six degrees of freedom per node, and has the option of suppressing rotational degrees of freedom at nodes to simulate pinned connections. The various payloads and other dead loads were represented by mass elements (MASS21 ). The coordinate system used in the finite element model is as follows (also shown in the figures below): X - Athwartship (positive in pott direction) Y - Vertical (positive upwards) Z - Longitudinal (positive in forward direction) The boundary conditions applied to the mast are as follows: UX=UY=UZ=O Ux = Uz=o at 1 deck level at 01 deck level in the X, Y, and Z
Main Legs:
The static analysis consisted directions. The accelerations Case i. Case ii Case iii
of three load cases of base accelerations applied are as follows: > : :8.48
8 g Athwartship Shock (m/s2): 18 g Vertical Shock (m/s2]: 8 g Longitudinal Shock (m/s2):
a~=O
aY = 9.81 av = 186.39 aY = 9.81
a,=O a, = O a, = 78.48
For the dynamic load case, translational master degrees of freedom are selected at the corner nodes of each level and the first 5 natural frequencies and the corresponding mode shapes are extracted.
C-27
.-.
FEA Example
No.
Title
: Mast
65 Nodes
217 370 Elements Degrees of Freedom
A
.
I
x
=xample
3a - Pinned Joints; Typically
one element
per member
65 Nodes
217 370 Elements Degrees of Freedom
A
:xample
I
x
,3b - Rigid Joints;
Typically
one element
per member
C-28 ,,,. --.,,
FEA Example
No.
Title
: Mast
200
352 1180
Nodes Elements Degrees of Freedom
\ \
I
Example 3C - Rigid Joints;
I
Typically
I
two elements per member
C-29
L. ....
FEA Example No.

DISCUSSION
3 I
Title
: Mast
OF RESULTS
The displacements for the three static load cases are summarized in Table C3.2, When the two modelling approaches (pinned joint versus rigid joint models) are compared, the model with pinned joints predicts the most flexible structure with the most displacements for every load case. Also, in some cases, the maximum displacement is predicted at a location different from the one predicted by the rigid joint model. In the second load case (Vertical shock) the displacement in Y direction, although at the same location for all three models, is excessively overpredicted by the pinned joint model. The maximum vertical deflections occur at the centre of the horizontal cross braces. Under vertical shock loading, these members act similar to beams subject to a unform distributed load (ie, inertial loading) for which the maximum deflection in the simply supported case (ie. pinned ends) is five times that for the fixed ends case. Table C3.3 lists peak stresses, As expected, the axial stresses are approximately the same for the two approaches. However, the bending stresses at mid-span of horizontal members and cross braces are significantly more in the pinned joint model. This is again due to the different end conditions in the two modelling methods. The model with simply supported end conditions naturally predicts higher moments at mid-span. Among the two models with fully rigid connections, the predicted maximum stresses are similar. The probable disadvantage with the one element per member model is that the stress at the centre of the member will not be calculated. It is possible that some members might have peak stresses at the centre as opposed to the ends if the members are also subject to local transverse loads (eg, wind loads, high inertial loads, equipment support loads). The natural frequencies and mode shapes for the two approaches are similar (see Table C3.4). Figure C3i 1 shows the first five mode shapes obtained from example 3b. The variations in deflection and some stress results between the pin jointed and rigid jointed models are significant. Hence, extreme care and proper judgement is needed in deciding on the right modelling approach for the problem.
C-30
FEA Example
No.
Title : Mast
Lx
Lx
Figure C3. 1 The first five mode shapes
c-3 1
Table Real Constant Set No. T 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
C3. I:
Geometry
Properties Real Constants
Member or Component Description Main Legs -1 Main Legs -02 Deck to 02 Deck Deck to Level B
Cross Section or Size Area i (10* m2) I10izm4] Iw [10-6 m4) 29.9700 20.7400 15.3100 6.1000 4.0540 3.1600 1.2490 1.9900 1.7000 0.6840 0.6370 1.9900 1.2490 0.6840 0.6370 0.2771 TKZBI (10-3 m) 92.10 90.13 88.90 63.50 61.91 61.91 46.00 50.80 50.80 38.10 36.51 50.80 46.00 38.10 36.51 30.20 TKYBI (10-3 m] 92.10 90.13 88.90 63.50 61.91 61.91 46.00 50.80 50.80 38.10 36.51 50.80 46.00 38.10 36.51 30.20
7.25 OD
7.1 OD
X
6.0 ID
ID ID ID ID
8392.0 5750.0 4236.0 3520.0 2344.() 1780.0 1306.0 1730.0 1450.0 1069.0 t 1100.0 1730.0 ID 1306.0 1069.0 1100.0 t 693.0
29.9700 20.7400 15.3100 6.1000 4.0540 3.1600 1.2490 1.9900 1.7000 0.6840 0,6370 1.9900 1.2490 0.6840 0.6370 0.2771
6.25
Main Legs - Level B to Level D Main Legs - Level D to Level F Main Legs - Level F to Top V Breces -02 Deck to Level D
7.0 OD X 6.375 5.0 OD X 4.25 4.875 4.875 3.625 4.0013
OD X 4.375 OD X 4.5 OD X 3.25 X 0.226 t ID ID ID
V Braces - Level D to Level G Vr Braces - Level G to Top Horizontals - Level A to Level D Horizontals - Level E to Level G Horizontals - Level MG Horizontals - Level MG X Braces - Level A to Level D X Braces - Level E to Level G X Braces - Level MG Platform
iD
4.0 OD X 3.625 3.0 OD X 2.635 2.875
OD X 0.203 t
4.0 OD X 0.226 3.625 OD X 3.25
3.0 OD X 2.635 2.875 2.375 OD X 0.203
ID t
OD X 0.154
TABLE
C3.2
Comparison
of displacements
for the Mast
finite
element
analyses
Max. Displacement Description Example 3b -rigid joints with 1 element per member
Imrn] Example 3C -rigid joints with 2 elements per member Location
Example 3a -pinned joints
Athwartshi~ ~ 6X 6, 6, Vertical (Yl

n
-15.13 -1.602 1.42
-15.07 -1.59 1.40
-15.07 -1.60 1.41
outer tip of spur frame outer tip of spur frame spur frame at main leg junction
LJ
w
ShLIG!l ax 6, 5= Longitudinal ~ 6, 6, 13z 1 2 3 4
-3.46 -74.24 3.093
-0.76 -16.59 3.07
-0.76 -16.75 3.08
middle of horizontal member - level 2 centre of X brace - level 2 horizontal member at mid span (top of mast)
-0.37 3.94 -27.74 is -26.7 is -3.91 is -3.67 is 0.76
-0.37 3.93 -14.56 at the middle of horizontal
-0.374 3.94 -14.56 member
outer tip of spur frame outer tip of spur frame
spur frame at main horizontal at mid span - level 2 - level 2 - level 4
The maximum The maximum The maximum The maximum
at the centre of cross brace member at the middle of horizontal at the middle of V-brace member - level 2
TABLE
C3.3
Comparison
of stresses
for the Mast
finite
element
analyses
Stress (iUIPa} Description Example 3b rigid joints with 1 element per member Example 3C -rigid joints with 2 elements per member Location
Example 3a pinned joints
AthwartshiD (21 shock Axial stress (OX) Bending stress (@ Bending stress (a~z) Vertical (Y} shock Axial stress (OX) Bending stress (a~Y) Bending stress lab.) Longitudinal (2[ shock
Axial stress Bending (uJ
*lt35 A36f A~58
*104 *39 *6 I
*104 *39 *58
Lower V braces Lower V braces at main ieg junction Horizontal members at mid-span
-81,+41 + 2452 * 343
-81, +40 A163 *41
-81,+40 A163 *4 I
Main legs, spur frame diagonals X braces at main leg junction Spur frame at main lag junction
*g8 *174
*88 *31
&87 +31
Lower V braces Spur frame at main
stress (oJ
leg junction
Bending stress (o~,)
+193
A 60
*58
Horizontal members at mid-span
1 2 3 4
Main Legs at level 1 Cross Braces at mid-span Main Legs at mid-span Main Legs at mid-span
TABLE
C3.4
Comparison
of frequencies
for the Mast
finite
element
analyses
Frequency (Hz) Mode Example 3a -pinned joints 13.30 Example 3b -rigid joints with 1 element per member 13.31 Example 3C -rigid joints with 2 elements per member 13.30 13.76 21.53 34.41 38.16 I I Mode Shape
Bending about Z- axis (1st mode) Bending about X-axis Twisting (Ist mode)
21
3 4 5
13.76 21.56 34.51 38.33
13.77 21.53 34.39 38.13
about Y-axis (2nd mode)
Bending about X-axis

Bending about
Z- axis (2nd model
,..
C-36
Appendix
Ship Structure Benchmarks for Assessing FEA Software
Benchmark BM-I -a BM-1-b BM-2-a BM-2-b BM-2-C BM-2-d BM-3 BM-4 BM-5
TilJg Opening With Insert Plate (4-Node Plate Opening With Insert Plate (8-Node Plate Stiffened Panel (in-Plane Beam Elements Stiffened Panel (Off-Set Beam Elements Stiffened Panel (4-Node Plate Elements) Stiffened Vibration Panel (8-Node Plate Elements) Isolation System Elements) Elements) with 4-Node Plate Elements) with 4-Node Plate Elements)
ME
D-2 D-7 D-9 D-15 D-17 D-19 D-21 D-24 D-29
Mast Structure Bracket Detail
WARNING The benchmark problems and associated FEA models presented in this document are intended for the express purpose of evaluating FEA software for ship structural analysis applications. While attempts have been made to ensure that the FEA models folio w good modelling practice, they should not necessarily be regarded as appropriate for any other purpose than that for which they are intended.
D-1
,..,, ,,
3enchmark No. : Nnalysis Type :
BM-1 -a 2D Static
Benchmark Title : Element Type(s) :
Opening with Insert Plate 4-Node Plane Stress 2-Node Line (Axial Stress)
zroblem Description: 14rectangular deck opening with rounded corners is reinforced with insert plates at each corner. 3etermine the maximum von Mises stress in the 20 mm insert plate and the 10 mm deck pIate. Sketch of Benchmark Problem :
2 z
o 0 z II
a) Deck Opening Wtih Insert Plate
DeckPlate t=lOmm T 1000

lz-looo-q
Stiffeners
T 400
300
Insert Plate
t=20mm I
600
600 1
4004
&+ s o 0
1200 ;-b
x
300 R
4
b 9
1200~700+ k -4 b) Detail of Shaded Region of Deck Opening L
Material Properties : E = 207000

v = N/mm2
Geometric Propenies :
Deck Plate Insert Plate t=l Omm
Loading :
P. = 100 (Applied N/mmz as nodal force
0.3
t=20mm
Stiffeners Line Elements
A = 1575 mm2 A = 1 mm2
loading)
D-2
3enchmark No. : Analysis Assumptions
BM-1-a :
Benchmark Title :
Opening with Insert Plate
2ue to symmetry, only one-quarter of the opening is modeled. The deck stiffeners are modelled ~sing axial stress line elements since only in-plane loading is considered. Finite Element Model :
el #12
node
node #l 37
No. of Nodes : No. of Elements 1. 2, 3, 4,

Boundarv :
200
212 120 4-Node 48 4-Node 25 2-Node 19 2-Node Plate Elements t= 10 mm Plate Elements t= 20 mm Line Elements A= 1575 mm2 Line Elements A= 1 mm2 (for stresses at free edge)
Deck Plate Insert Plate Stiffeners Line Elements

Conditions :
Ux
Oat
X=O
Uy=Oat Y= Oand Y= 1600 Uz = O at (X= O;Y= O), (X=O;Y=
1600),
(X= 2600; Y= O), and (X=2600,Y=
1600)
D-3
.,..,.. .
Benchmark No. :
BM-1-a
Benchmark Title :
Opening with Insert Plate Converged Solution 4 (ANSYS 5.1) SHELL93 LINK8
Finite Element Software
Results
ANSYS 5.1 SHELL63 LINK8
MSC I NASTRAN W!ndows 1 CQUAD4 CROD
ALGOR 3.14 TYPE 6 TYPE 1
FEA Software
Element Tv~es :
Maximum Str esses 1, Deck Plate
(MPa) 192.8 198.3 139,8 204.4 193,5 189.2 139.8 203.3 192.3 199,3 139.8 204.4 196.9 206.3 140.3 209,0
ue~v 1 (node # 10) (node #1 63) (el # 129)
2. Insert Plate u,~v 1 3, Stiffeners u, 2
4. Edge Elements u, 3 (el # 205) Maximum Deflections Ux Uy (mm) (node ,#1 37) (node # 1)
1,496. 0.1:57
1.496 0,157
1.496 0.157
1.506 0.157
Comments on Benchmark Results : 1. a,qv is the maximum von Mise$ m equivalent slress reported for the plate elements (section properties 1 and 2) .The values-presented are the nodal averaae d stresses within each group of elements of the same section propertyi The. nodal averaged stresses are obtained by extrapolating stresses at the element integration points to the node locations, and then averaging the values at each notle. Different FEAsoftware may use different extrapolation and averaging methods which can lead to slight differences in the nodal stress results, 2, a, is the maximum axial or direct stress inthe line elements.
3. The benchmark FE model includes line elements of small arbitrary area (section property 4 with A = 1 mm2) which ,are used .to obtain stresses around the free edge of the opening. The maximum axial stress reported in the line elements corresponds approximately to the maximum principal and von Mises stress at the edge of the opening, irrespective of the stress extrapolation method used for the plate elements.
4,
The converged solution for this benchmark was obtained using a more refined model of the same problem consisting of 8 node shell elements with ANSYS 5,1. The stress contour plot for the converged solution is shown on the following page. Note that the plot shows element stresses, ~ nodal averaged stresses, so as to permit presentation of the results for the two plate thicknesses cm the same plot. Although the plot shows slight discontinuities in the stress contours, these are mainly away from the areas of interest. The difference between the maximum element stresses and the nodal averaged stresses is minimal at the two locations reported in the above table. There is a real stress discontinuity at the border between the insert plate and the deck plate due to the abrupt change in plate thickness. The stress contour values are in units of MPa. The MX on the plot signifies the location of maximum stress.
D-4
\L..,-
Benchmark No. : Analysis Type :
BM-1-b 2D Static
Opening with Insert Plate 8-Node Plane Stress 2-Node Line (Axial Stress)
Problem Description: Repeat Benchmark 1-a using a coarser mesh with 8-node elements in place of 4-node elements.
Finite Element Model :
el # 42
node # 19 el # 93
Y L
node #149 .. ~~ :
200
: 103
No.
of Elements
1. Deck Plate 2. Insert Plate 3. Stiffeners 4. Line Elements

Boundarv Conditions :
41 18 22 22
8-Node 8-Node 2-Node 2-Node
Plate Elements t =10 mm Plate Elements t =20 mm Line Elements A= 1575 mm2 Line Elements A= 1 mm2 (for stresses at free edge)
As defined for BM- l-a, Loadinq : As defined for Benchmark 1-a,
D-7
.-._.. -
3enchmark No. :
BM-1-b
Benchmark Title :
Opening with Insert Plate Converged Solution4 (ANSYS 5.1) SHELL93 LINKS
Results
MSC I NASTRAN Windows 1 CQUAD8 CROD
ALGOR
~:
NA*
M!All 1i Deck Plate u,~v 1
(MPa) (node # 30 ) (node #1 72) (el # 42) (@l# 93) (mm) (node #149) (node # 19) 1,505 0.157 1,505 0.157 1,506 0.157 195.6 207,8 140,3 207.8 195.6 204.5 140.3 207,8 196.9 206.3 140.3 209,0
2. Insert Plate a.~vl 3. Stiffeners 0,2
4, Edge Elements us Maximum Deflections Ux Uy
Comments on Benchmark Results : *ALGOR does not include 8-node plate elements for stress analysis. 1. u,, is the maximum von Mises or equivalent stress reported for the plate elements (section properties 1 and 2). The values presented are the nodal averaaed stresses within each group of elements of the same section property. The nodal averaged stresses are obtained by extrapolating stresses at the element integration points to the node locations, and then averaging the values at each node, Different FEA software may use different extrapolation and averaging methods which can lead to slight differences in the nodal stress results. o~ is the maximum axial or direct stress in the line elements. The benchmark FE model includes line elements of small arbitrary area (section property 4 with A = 1 mm2) which are used to obtain stresses around the free edge of the opening. The maximum axial stress reported in the line elements corresponds approximately to the maximum principal and von Mises stress at the edge of the opening, irrespective of the stress extrapolation method used for the plate elements. The converged solution for this benchmark was obtained using a more refined model of the same problem consisting of 8 node shell elements with ANSYS 5.1. The stress contour plot for the converged solution is shown on Page D-5. Refer to the BM-1 -a results for further discussion of the converged solution.
2. 3.
4.
D-8
......
% .,.
3enchmark No. : Analysis Type :
BM-2-a 3D Static 3D Modal
Stiffened
Panel
4-Node Shell 2-Node Beam (In plane of plate)
Jrob[em Description: A rectangular stiffened paneI is subject to a uniform pressure load applied to its surface, 3etermine the maximum deflection, stresses and natural frequencies for the panel. Sketch of Benchmark Problem :
Benchmark Problem 2: Stiffened Panel
Material Properties : E = 207x109 N/m2 v = 0.3 p = 7850 kg/m3
Geometric Properties : Plate Stiffeners t=l Omm 15 OX1O.5FB
Loading : P= = 9810 Pa
D-9
..
I
Benchmark No. : Finite Element Model : BM-2-a
1
Benchmark Title :
Stiffened
Panel
nods#133
No. of Nodes : No. of Ele ments : 1. Panel 2, Stiffeners
143 144 120 24 A= l,, = IYY= IX, =

4-Node 3-D Plate Elements \ nod,# 2
Ax
t=l
Omm
2-Node 3-D Beam Elements 0,001575 m2 53.35 x 10-E m4 ** 10.19 x10-8m4 0.0553 x 108 m4 (Torsion) Y! Y~ Z; Z~ = = = = 0,1352 m 0.0148 m 0.00525 m 0.00525 m
* * In-Plane Beam elements l,, includes 40 t effective Boundarv Conditions : 1.tic Analvsis
plate width.
- All nodes fixed at edges along x=O and along y=O, . Symmetry about YZ plane along edge at x = 2.250 m Symmetry about X2 plane along edge y = 1.500 m - All nodes fixed at edges along x= O and along y= 0, Symmetry about YZ plane along edge at x = 2,250 m Antisymmetry about X2 plane along edge y = 1,500 m
2. Modal Analvsi5*
This benchmark test only requires calculation of the first four natural frequencies for symmetry / antisymmetry boundary conditions, In order to capture all modes of vibration, the modal analysis of the quarter model would also have to consider symmetry / symmetry, antisymmetry / symmetry, and antisymmetry / antisymmetry boundary conditions.
D-10
..
3enchmark No. :
13 M-2-a
Benchmark Title : MSC I NASTRAN Windows 1 CQUAD4 CBAR
Stiffened ALGOR 3,14
Plate Converged Solution (ANSYS 5.1) SHELL93 SHELL93
Results
ANSYS 5.1 SHELL63 BEAM4
~lement Tvms
Plate Stiffeners (MPa) (node # 2) (MPa) (node #1 33) (node #1 44) (mm) (node W 18) 5: (Hz) (Hz) (Hz) (Hz)
TYPE 6 TYPE 2
ylll u tresses I. Plate a,~v 2 2. Stiffeners 0, 3 Tension Compression blaximum Deflections Uz 4 Natural Frequencies 1t Mode 2n~ Mode 3~ Mode 4h Mode
39.3
38.2
36,5
42.1
69.0 -135.8
69.0 -135.8
69.0 -135.0
61.3 -126.5
3,30
3,29
3.29
3.50
36,5 60.9 100.1 110.2
36,5 61,1 100.4 111.4
36.6 61.2 102.4 111.9
35.9 61,0 96,5 106.5
1. The converged solution results were obtained using a refined mesh model with 8-node shell elements on ANSYS 5,1, The von Mises Stress contours for the converged model are shown on Page D-13. The stress contours are in units of Pa (N/m2). 2. The maximum stress in the plate occurs at the middle of the long fixed edges (node 2). Reported are the maximum nodal averaged von Misas stress of the top or bottom surface of the plate elements. Note that different FEA programs may use different conventions for defining the top and bottom surfaces of plate elements, Also, different FEA programs use different extrapolation and averaging techniques for computing plate / shell element stresses which may lead to slight differences (refer to BM-1 -a for discussion). 3. Reported are the maximum stresses in the beam elements (axial stress + bending stress). The maximum tensile stress occurs at the centre of the middle stiffeners (node 133). The maximum compressive stress occurs at the fixed ends of the middle stiffeners (node 144). % The maximum out-of-plane deflection (Uz) occurs at the centre of the panel (noda 11 8). in deflection and stress results relative to the converged model are due mainly to Differences the simplifying assumption of 40 t effective plate width used in defining the beam properties. 5. The frequencies and mode shapes for symmetry / antisymmetry boundary conditions from the The mode shapes predicted by the BM-2-a FEA converged model are shown on Page D-12. models are the same as those for the converged model. The frequencies predicted by the BM-2-a model deviate slightly from those predicted by the converged model, particularity for the 3rd and 4th modes. These are more complex modas involving torsion of the stiffeners for which the beam + plate element model is probably too simplified. However, the plate + beam model gives very good predictions for the first two modes.
D-1 1
Benchmark No. :
BM-2-a
Benchmark Title :
Stiffened
Plate
lS Mode
:35.9
Hz
2nd Mode :61.0
Modal Analysis Results of Converged Model for
EM-2
(ANSYS
5.1) I
D-12 ..-,-----
D-14
..
BM-2-b 3D Static 3D Modal
Stiffened
Panel
4-Node Shell 2-Node Offset Beam
Problem Description: Repeat BM-2-a using 2-node offset beams in place of in-plane beam elements. Finite Element Model :
nod. #1
noti #ha
nada #133
No. of Nodes :
143 144 npds2 # 120 24 4-Node 3-D Plate Elements 2-Node 3-D Beam Elements** Y, = 0.075 m Y~ = 0.075 m Zt = 0.00525 m Z~ = 0.00525 m m in global Z direction. A
No. of Elements : 1. Panel 2. Stiffeners
t = O.OIOm
A = 0,001575 mz IZz = 0.0145 x 10-E m4 Iw = 2.95 x 10E m4 Ixx = 0,0553 x 10-6 m4 (Torsion)
* * Beam element centroid off-set 0.075 Bounda rv Conditions : 1. Static A nalvw
- All nodes fixed at edges along x=0 and along y= O. - Symmetry about YZ plane along edge at x = 2.250 m - Symmetry about X2 plane along edge y = 1.500 m - All nodes fixed at edges along x= O and along y = O. - Symmetry about YZ plane along edge at x = 2.250 m - Antisymmetry about X2 plane along edge y = 1.500 m
2. Modal Ana Ivsis *
This benchmark test only requires calculation of the first four natural frequencies for symmetry / antisymmetry boundary conditions.
D-15
Brmchmark No. :
BM-2-b
Benchmark Title : MSC / NASTRAN Windows 1 CQUAD4 CBEAM
Stiffened
Plate Converged Solution 1
Results
ANSYS 5.1 SHELL63 BEAM44
(ANSYS 5.1)
SHELL93 SHELL93
Element TvDes :
Plate Stiffeners (MPa) (node # 2)

(MPa)
Maximu m Stresses 1. Plate O,qv 2 2, Stiffeners

Tension
42,1
38.2
34.4
42.1
UX 3
Compression Maximum Deflections Uz 4 Natur I Fr ~5: 1 Mode 2nd Mode 3rd Mode 4th Mode
(node #1 33) (node #144) (mm)
70,3 -153.7
70.4 -154.0
70.3 -153.7
61,3 -126.5
(node #1 18)
3.42
3,41
3.41
3.50
(Hz) (Hz) (Hz) (Hz)
36.3 61.1 97.0 107,0
36,3 61,2 95.7 106.8
36.5 61.7 101.9 111.9
35.9 61.0 96.5 106.5
1. The converged solution results were obtained using a refined mesh model with 8-node shell The von Mises Stress contours for the converged model are shown elements on ANSYS 5.1. on Page D-13. 2. The maximum stress in the pIate occurs at the middle of the long fixed edges (node 2). Reported are the maximum nodal averaged von Mises stress of the top or bottom surface of the plate elements. Note that different FEA programs may use different conventions for defining the top and bottom surfaces of plate elements. Also, different FEA programs use different extrapolation and averaging techniques for computing plate / shell element stresses which may lead to slight discrepancies (refer to EM-1-a for discussion). Reported are the maximum stresses in the beam elements (axial stress + bending stress). The maximum tensile stress occurs at the centre of the middle stiffeners (node 133). The maximum compressive stress occurs at the fixed ends of the middle stiffeners (node 144). The off-set beam element introduces an artificial moment into the problem which results in over prediction of the stresses and under prediction of deflections. This effect also influences stress results for the plate elements, Refer to Example 1, Appendix C for further discussion of this effect. The maximum out-of-plane deflection (Uz) occurs at the centre of the panel (node 11 8).
3.
4. 5.
The frequencies and mode shapes for symmetry / antisymmetry boundary conditions from the converged model are shown on Page D-12, The mode shapes predicted by the BM-2-b FEA models are the same as those for the converged model.
D-16
3enchmark No. : rhalysis Type :
BM-2-C 3D Static 3D Modal
Stiffened
Panel
4-Node Plate
Problem Description: depeat BM-2-a using 4-node plate elements to model the stiffeners and plate explicitly.
Finite Element Model :
nodn #
ude 9118
#172
No.
of Nodes
: :
120
No. of Elements
Panel
nod. # 2
4-Node 3-D Plate Elements 4-Node 3-D Plate Elements
t=l
Omm
Stiffeners
48
t = 10.5 mm
3oundarv Conditions : 1. Static Analvsis - All nodes fixed at edges along x=O and along Y=O. - Symmetry about YZ plane along edge at x = 2.250 m . Symmetry about X2 plane along edge y = 1,500 m - All nodes fixed at edges along x=0 and along y = O. - Symmetry about YZ plane along edge at x = 2.250 m - Antisymmetry about XZ plane along edge Y = 1.500 m
2. Modal Analvsis
x.
D-17
L..
Benchmark No. :
BM-2-C
Benchmark Title : MSC I NASTRAN Windows 1

CQUAD4 CQUAD4
Stiffened ALGOR
Plate Converged Solution 1 (ANSYS 5.1)

SHELL93
Firrite Element Software
Results
ANSYS 5.1
SHELL63 SHELL63
3.14
Element
TvDes
Plate Stiffeners
TYPE TYPE
SHELL93
Maximum Stresses 1. Plate u,~v 2 2, Stiffeners 0, 3 Tension Compression
(MPa) (node # 2) (MPa) (node #172) (node #170) (mm) (node #1 18) 6: (Hz) (Hz) (Hz) (Hz)
42,3
41,3
39.3
42.1
68,9 -126.0
69,0 -126,0
68.2 -124.0
61.3 -126.5
Maximum Deflections Uz 4 Natural Frequencies 1t Mode 2d Mode 3rd Mode 4th Mode
3.47
3.43
3.42
3.50
36.1 60,8 95.0 104.9
36.2 61.1 94.9 105.8
36.1 61.2 97.4 106.3
35.9 61.0 96.5 106.5
1. The converged solution results were obtained using a refined mesh model with 8-node shell elements on ANSYS 5.1. The von Mises Stress contours for the converged model are shown on Page D-13. 2. The maximum stress in the plate occurs at the middle of the long fixed edges (node 2). Reported are the maximum nodal averaged von Mises stress of the top or bottom surface of the plate elements. Note that different FEA programs may use different conventions for defining the top and bottom surfaces of plate elements. Also, different FEA programs use different extrapolation and averaging techniques for computing plate / shell element stresses which may lead to slight discrepancies (refer to EM-1-a for discussion), Repor-ted are the maximum nodal averaged stresses, crX, in the stiffener plate elements (maximum of top or bottom surface stress), The maximum tensile stress occurs at the centre of the middle stiffeners (node 172), The maximum compressive stress occurs at the fixed ends of the middle stiffeners (node 170), deflection (Uz) occurs at the centre of the panel (node 11 8).
3.
4. The maximum out-of-plane

5.
The frequencies and mode shapes for symmetry / antisymmetry boundary conditions from the converged model are shown on Page D-12. The frequencies and mode shapes predicted by the EM-2-c FEA models are very similar to those from the converged model,
D-18
.,\-_ ,,.,
lenchmark
No. :
BM-2-d 3D Static 3D Modal
Stiffened
Panel
inalysis Type :
8-Node Plate
)roblem Description: Iepeat BM-2-a using 8-node plate elements to model the stiffeners and plate explicitly. finite Element Model :
node #174
>
node # 176
Y L
No.
of Nodes :
199 56 40 16 8-Node 3-D Plate Elements 8-Node 3-D Plate Elements t=l Omm
Vo. of Elements : Panel Stiffeners Boundarv Co nditions : 1. Static Analwk
t = 10.5mm
- All nodes fixed at edges along x=O and along Y=O. - Symmetry about YZ plane along edge at x = 2.250 m - Symmetry about XZ plane along edge y = 1.500 m
2. Modal Analvsis*
- All nodes fixed at edges along x= O and along y = O. . Symmetry about YZ plane along edge at x = 2.250 m - Antisymmetry about XZ plane along edge y = 1.500 m
D-1 9 ..
!
Benchmark No. :
BM-2-d
Benchmark Title : MSC I NASTRAN Windows 1 CQUAD8 CQUAD8
Stiffened
Plate Converged Solution 1 (ANSYS 5.1) SHELL93 SHELL93
Results
ANSYS 5.1 SHELL93 SHELL93
ALGOR
Element Tv~es :
Plate Stiffeners (MPa) (node # 2) (MPa) (node #1 76) (node #1 74) (mm) (node #1 22) 5: (Hz) (Hz) (Hz) (Hz)
NA*
Maximum Stresses 1. Plate o,~v 2 2. Stiffeners cq 3 Tension Compression
41.7
41.7
42.1
69.9 -143.0
69.9 -143.0
61.3 -126.5
Maximum Defle ctionq Uz 4 Natural Frequencies It Mode 2nd Mode 3rd Mode 4h Mode
3,49
3,49
3,50
36,0 61.0 96.6 105.9
36,0 61,0 96.1 105.6
35,9 61,0 96.5 106.5
*ALGOR does not include 8-node plate elements for stress analysis. 1, The converged solution results were obtained using a refined mesh model with 8-node shell elements on ANSYS 5.1. The von Mises Stress contours for the converged model are shown on Page D-13. 2.
stress in the plate occurs at the middle of the long fixed edges (node 2). The maximum Reported are the maximum nodal averaged von Mises stress of the top or bottom surface of the plate elements. Note that different FEA programs may use different conventions for defining the top and bottom surfaces of plate elements. Also, different FEA programs use different extrapolation and averaging techniques for computing plate / shell element stresses which may lead to slight discrepancies (refer to EM- I -a for discussion),
3.
u,, in the stiffener plate elements (maximum of top or bottom surface stress). The maximum tensile stress occurs at the centre of the middle stiffeners (node 176). The maximum compressive stress occurs at the fixed ends of the middle stiffeners (node 174). The maximum out-of-plane deflection (Uz) occurs at the centre of the panel (node 122).
Reported are the maximum nodal averaged stresses,
4. 5,
The frequencies and mode shapes for symmetry / antisymmetry boundary conditions from the converged model are shown on Page D-12. The frequencies and mode shapes predicted by the BM-2-d FEA models are very similar to those from the converged model, despite the relative coarseness of the mesh of the former,
D-20
BM-3 3D Modal
Machinery Vibration
Isolation System
3D Beams 1 DOF Springs (in X, Y, Z directions) Mass (with Rotational Inertia)
Problem Description: Determine the natural frequencies for this generator vibration isolation system. Sketch of Benchmark Problem :
IsolatorStiffness 1$= 350 kN/m I-$ = 350 kN/m &= 800 kN/m a) Generator VibrationIsolationSeat @~=
++.+05+.5+
T
. 0.7 1 b) Plan VW of Seat Frame Material Properties : Gaomatric 207x103 0.3 7850 kg/m3 Generator elements 2, Rigid Links E = v = 207x104MPa 0,3 centroid. modelled and point as rigid mass at link MPa Propetiies :
IH,
O.015m
u7.5x 10sm4
1~~u IOX 10-5m4 u 17,5x1O-sm4
@~=0,010m2 1==5.0x105m4 lW=7,5x104m4 ln2 = 12.5 x 10-sm4

is
~earn~ Z (Verlical) ~ .... .. n z y
Loading
1. Steel
E = v = p =
Refer
to above
sketch,
Not
Applicable.
p = O kg/m3
D-21
Benchmark No. : Finite Element Model :
EM-3 I
Benchmark Title :
Machinery Vibration
Isolation System
Mass
Rigid Links \,
\ \
[
1
1 i /
B8~ms
(Section Property 21
Beams (Section x
Property 1)
L
No. of Nodes : Nrj 14 5 14 14 14 1 51
81 90 Beams (Section Property 1) Beams (Section Property 2) Springs (X-Direction) Springs (Y-Direction) Springs (Z-Direction) Mass Rigid Links
Boundarv Conditions : Isolator springs fixed at deck seating level.
D-22
... . ...
!.> ,,
Benchmark No. : EM-3
Benchmark Title :
Machinery Vibration Isolation System MSC 1 NASTRAN Windows 1 CBAR CONM2 CROD
Results
ANSYS 5.1 BEAM4 MASS21 COMBIN14
ALGOR 3.14 TYPE 2 TYPE 1 & 7
FEA Software
Element Tv~eq :
Total Mass and C of G Location : Total Mass Cof G x Y z (kg) (m) (m) (m) 2545.7 1,0000 0.3500 0.4066 2!545,7 1,0000 0,3500 0.4066 2545.7 1,0000 0.3500 0.4065
Modes and Frequencies (Hz~ 1 2 3 4 5 6 7 8 9 10 11 12 Translation in Y direction Translation in X direction Translation in Z direction Rotation about Z axis Rotation about Y axis Rotation about X axis Translation in X direction Rotation about Z axis Translation in Y direction Translation in Z direction Rotation about Y axis Rotation about X axis (1 $) (1 ) (1 t) (1 t) (1 ) (1 ) (2) (2) (2nd) (2nd) (2) (2) 2,85 3.60 6.30 6.62 9.61 11.12 14.76 15.28 16.92 21.51 22.86 23.12 2,85 3,60 6.30 6.62 9.61 11,12 14,76 15,28 16.92 21.51 22,86 23.12 2.80 3.66 6.30 6.98 10.04 11,45 ~4.a9 16.61 16.79 21.51 23.60 24.44
Comments on Benchmark Results : Modes 1 to 6 involve vibration modes with the generator and raft masses moving in phase, while the two masses are out-of phase for modes 7 to 12,
D-23
BM-4 3D Static 3D Modal
Mast Structure 3D Beam 3D Spar Mass
Problem Description: Determine the stresses, displacements, natural frequencies and modes under the specified loading conditions for the mast structure shown in the sketch below. Sketch of Benchmark Problem :
1 S* $ ~ n f
.. .........
Material Properties : 1, Steel E = 207xI OgN/m2 v = 0.3 p = 7850 kg / m3 E = 70x10g N/mz v = 0.3 p = 2900 kg / m3
Geometric Properties : Refer to table of section properties.
Loading : Accelerations a, = 5 m/s2 aY = 5 m/s2 aZ = 15 m/s2 N
2. Aluminum (pole mast)
Nodal Forces FX =3000 (Applied on all nodes)
D-24
Benchmark No. :
BM-4
Benchmark Title :
Mast Structure
Member Section Properties Section Description No. 1 2 3 4 5 8 9 10 Main Legs Pole Mast Support Vertical Braces 0.09200 Main Horizontals Pole Mast (Aluminum) Horizontal Braces Platform Braces Platform Chords
O. Dia. (m) 0.12700 0.09200 1,306 0,07620 0.24130 0.07302 0.06040 0.06040
Area (xl 03 m2) 3,520 1,306 1.069 4,887 1.100 0.693 0,693
122 & lyy lxx (xl OG m4 (xl O-Em4) 6,100 1.249 0.684 33.70 12.2 2.50 1.37 67.4
Element Type Beam Beam Spar Beam Beam Spar Spar Beam
No, Elems 32 8 32 32 5 16 10 12
0.2771
0.554
Finite Element Model : The main legs, polemast, main horizontals and platform frame chords are modelled as continuous beams (ie. with full continuity), while the various brace members are modelled as spars with pinned ends,
~: No.
Of
67 Elements : 150
mu ndarv Condition s : UX, UY, & UZ translations of node at base of each leg restrained.
Static An alvsh Loads : Nodal force of 3000 Accelerations N in X direction (Fx) at every node, a, = 5 m/s2, aY = 5 m/s2, a, = 15 m/s2.
D-25
,,
,.
Benchmark No. :
BM-4
Benchmark Title :
Mast Structure
Plot of Finite Element Model Showing Critical Element Numbers :
1-
D-26
%.
Benchmark No. :
EM-4
Benchmark Title :
Mast Structure
Results
ANSYS 5.1 BEAM4 LINKS MASS21
MSC 1 NASTRAN Windows 1 CBAR CROD CONM2 1415.8 0.0336 0.0003 2.3797 12.00 -0.37 -0.62 -190921 7079 21237
FEA Software
Element Tv~es :
Total Mass : Centre of Gravitv:
(kg) (m)
m x Y z UX Uy LIZ FX F, FZ (node #63) (node #63) (node #56)
1415.8 0.0336 0.0003 2.3797 12.00 -0.36 -0,62 -190920 7079 21236
1418.7 0.0335 0,0003 2,3841 12.65 -0.41 -0.65 NA*
M ~
(mm)
Total Reactio n Forces : (N)
Stresses
1. Main
(MPd
Leas
2.
Max. Tensile Max. Compressive Maxi Tensile Max. Compressive Maxi Tensile Max. Compressive Maxi Tensile Max. Compressive Max. Tensile Max. Compressive Max. Tensile Max. Compressive Max, Tensile Max, Compressive Max. Tensile Max. Compressive
(el #1) (el #5) (el #143) (el #1 42) (el #45) (el #61) (el #74) (el #68) (ei #1 36) (el #1 36) (el#lll) (el #1 09) (el #130) (el #1 22) (el #1 16) (el #1 27)
33.70 -36,09 99.42 -108.96 34.94 -35.54 48.41 -38.11 53.53 -53.88 10.77 -4.32 4.60 -15.64 71,90 -73,43
33.67 -36.11 99.41 -108.95 34,94 -35,54 48.40 -38.09 53.54 -53.86 10.77 -4.32 4.61 -15.64 71.92 -73.41
33.72 -31.35 95,85 -97.76 38.15 -37.78 47.81 -39.61 49.98 -50801 10897 -4,29 4.73 -16.40 75.97 -74.85
z pole Mast SUDDort
3. Ve rtical Braces
4. Main Horizontals
5. Pole Mast
5. Horizontal Braces
~ 3
] O. Platform Chor~
D-27
..-. ,,,
Benchmark No. :
EM-4
I
Results
Benchmark Title :
Mast Structure
ANSYS 5.1
MSC I NASTRAN Windows 1
ALGOR 3.14
Modes a nd Freauenc!e s : 3 (Hz) 1
Pole Mast Cantilever Bending

Pole Mast Cantilever Bending Local Bending of Main Horizontals Platforms Bending in X Direction
20,75
20.79 41.13 47.46
20.76
20.80 41.13 47.46
20,72
20.76 41.13 47.45
2
3 4
Comments 1.
on Benchmark Results : The
The maximum deflections in the X and Y directions occur at the top of the polemast. maximum vertical deflection occurs at the starboard spur frame.
2.
The stresses listed are the maximum combined (axial + bending) stress in the beam elements. Note that the maximum combined stresses are calculated assuming -that the beam element section is rectangular, although the members in this particular problem are actually of circular cross section, This is a limitation of most FEA programs, although some programs may allow calculation of element stresses assuming a circular cross section. The simplified assumption of a square or rectangular cross section is conservative. The first two vibration modes involve cantilever bending of the polemast in the X and Y The directions (they are in fact identical modes due to the symmetry of the polemast), third mode involves local bending of the main horizontal members of the mast. The fourth mode involves bending of the two platforms at the top of the mast in the X direction, The ALGOR program requires a separate module to output reaction forces which is not included with the basic solution module and, as such, reaction results were not available. In addition, the ALGOR program does not include mass elements for linear static analysis. Insteadr the inertia loads due to the payload masses were modelled by applying nodal forces at the appropriate locations. The difference in modelling approach and the inability to confirm the total applied loads may explain, in part, the differences in the ALGOR results com~ared to those obtained bv ANSYS and NASTRAN.
3.
4.
D-28
,, .
3enchmark No. : 4nalysis Type :
BM-5 3D Static
Bracket Detail
4-Node Thick Shell

(With Transverse Shear)
%oblem Description: 3etermine the maximum stress for the VLCC Top Bracket detail shown in the sketch below. Sketch of Benchmark Problem :
,n,cr.b~~~ r;
300 2.3 25 u
2A $,00
,
/ End B Ux=l.Omm Uy, Uz, Rx, Ry, Rz = o
100
Deck Longitudinal 300xIo0 T 13 mm Web
E E
gmmange
[:
ii u m # x 5 m a! 2 a! > m != ? jj
UY N w
~ ~
I%+
End c Ux
u
-0.5
mm
Uy=o 600
L-1
Material Properties : E = 207x103 v = 0.3 MPa
Geometric Properties : As defined in above sketch.
Loading : Applied displacement
constraints.
D-29
Ierrchmark No. : inite Element Model :
BM-5
Benchmark Title :
Bracket Detail
hick shell / plate elements with transverse shear flexibility are used to model the bracket, deck mgitudinal, and the web of the bulkhead stiffener, The transverse bulkhead, and upper deck Ire modelled using line elements of 40 t2 section area (9000 mmz for deck, 4850 mm2 for iulkhead). The flange of the bulkhead stiffener is modelled with line elements using the 2250 nm2 area of the flange. The areas of the flange line elements taper down to 923 mm2 at the md of the bracket, Line elements of a small arbitrary area (0.01 mm2) are included at the toe of he bracket for obtaining stresses.
A-J
N. ~: No. of FI ements :
Boundarv Conditions :
199 227
End C m
Translation bulkhead.
in Z direction restrained at nodes representing main deck and transverse
At end A of the model, all nodal degrees of freedom are fixed. At end B of the model, a 1 mm displacement all other nodal degrees of freedom are fixed.
is applied
in the positive X direction and
At end C of the model, a 0,5 mm displacement is applied in the negative X direction and the vertical displacement in the Y direction is constrained to zero. D-30
..
..
,, -,.!.
Benchmark No. :
BM-5
Benchmark Title :
Bracket Detail
Plot showing Critical Element Locations at Toe of Bracket :
\el#71 el #15a
Y I
l?D-31
Benchmark No. :
BM-5
Benchmark Title :
Bracket Detail Converged Solution 1 (ANSYS 5.1) SHELL93 LINKS
Results
MSC I NASTRAN Windows 1 CQUAD4 CROD
ALGOR
Element
TvDes
* NA
Plate Element Str esses a,~v 2 (M Pa) 1. 2, ~ 1. Bracket 2. Deck Long. Web
M ~
Bracket Deck Long, Web
(el # 71) (el # 105) o, (MPa)
209.3 248,9
209.6 247,6
203.5 243.4
(et # 158) (et # 211)

(mm)
119.8 235.5
121.5 236.0
133,0 240.1
xim
Ux Uy Uz
Rea
(node # 86)
1.000
(node #1 85) (node #106) : (N)
-0.339 -0,366
1.000 -0,336 -0.354
1.000 -0.348 -0.388
ction Forces at A Fx Fy Fz
-1194400 -28343 5967
-1194700 -28302 6019
-1191800 -26414 -5064
* ALGOR does not provide a thick shell element with transverse shear, 1, The converged solution results were obtained using a refined mesh model with 8-node shell elements on ANSYS 5,1. The von Mises Stress contours of the toe of the bracket for the converged model are shown on Page D-31. The stress contours are in units of MPa (N/mm2). 2. This particular bracket detail problem is complicated by the existence of a stress singularity at the end corner or toe of the bracket, In a linear elastic analysis, the stress at this point is theoretically infinite. Refining the finite element mesh gives progressively higher stresses which are meaningless. One method which is commonly used to get around this problem is to use the so called hot spot stress. In calculating the hot spot stress no account is taken of the weld geometry, and in an idealised finite element representation (ignoring the weld), the stress is equal to the value at about one plate thickness from the corner (Chalmers, 1993). In this benchmark, the hot spot stress is calculated two ways : a) b) Using element centroidal von Mises stresses at the elements 10 mm from the corner (elements 71 and 105, see figure on Page D-29); and Using line element stresses at 10 mm from the corner (elements 158 and 211).
The advantage of these methods are that they do not rely on the techniques used to extrapolate stresses to the node points which may vary for different FEA programs.
D-32
.-,
.. ,.,,
; ,. ,,,:!.4 .:
!, %...
Project Technical Committee Members

The following persons were members of the committee that represented the Ship Structure Committee to the Contractor as resident subject matter experts. As such they performed technical review of the initial proposals to select the contractor, advised the contractor in cognizant matters pertaining to the contract of which the agencies were aware, and performed technical review of the work in progress and edited the final report.
Chairman
LCDR Stephm Gibson
National efenceHeadquarters, D CANADA
Members
Mr. Rickard Anderson Mr. RichardSonnenschein LT PatLittle Mr. James White Mr. Natale Nappi Mr. JohnAdamchek
Military Sealift Command Maritime Administration U.S. Coast Guard U.S. Coast Guard Naval Sea Systems Command Carderock Division Naval Surface Warfare Center
Mr. Gary Horn Mr. Tom Ingram Stephen Yang Mr. William Siekierka
American Bureau of Shipping American Bureau of Shipping Defence Research Establishment Naval Sea Systems Command, Contracting Officers Techical Representative Atlantic
Dr.Robert Sielski
National Academy of Science, Marine Board Liaison
CDR Steve Sharpe
U.S. Coast Guard, Executive Director Ship Structure Committee

Guideline For Evaluation of Finite Elements and Results

Uploaded by

Copyright:

Available Formats

Guideline For Evaluation of Finite Elements and Results

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Guideline For Evaluation of Finite Elements and Results

Uploaded by

Copyright:

Available Formats

NTIS #PB96-1 53077

SSC-387 GUIDELINE FOR EVALUATION OF

FXEC UTIVE DIRECTOH CDR Stephen E. Sharpe, USCG U, S, Coast Guard

Mr. William J. Siekierka Naval Sea Systems Command

CAPT George Wright Mr. Walter Lincoln Mr. Rubin Sheinberg

AMERICAN Mr. Mr. Mr. Mr.

TRANSPORT Mr. Mr. Mr. Mr.

Glenn Ashe John F, ConIon Phillip G, Rynn William Hanzelek

John Grinstead Ian Bayly David L. Stocks Peter llmonin

NATIONAL ACADEMY ~D Dr. Robert Sielski NATIONAI ACADEMY

Dr. John Landes WELDING RESEARCH Dr. Martin Prager COUNCIL

u. s. NAVAL ACADEMY Dr. Ramswar Bhattacharyya

E OF TE CHNOLOGY CAPT Alan J. Brown

STUDENT MEMBER Mr. Jason Miller Massachusetts Institute of Technology

RECENT SHIP STRUCTURE

SSC-381 SSC-380 SSC-379 SSC-378

SSC-376 SSC-375 SSC-374 SSC-373 SSC-372 SSC-371

by A. Mansour and A. Thayamballi A State of the Art Summary by

Maintenance of Marine Structures: S. Hutchinson and R. Bea 1993

National Academy of Sciences - National Research Council

The COMMllTEE interagency

has technical cognizance

Ship Structure Committees of Tennessee,

research program. Knoxville, TN

John Landes, University

of Michigan, Ann Arbor, Ml Consultant, Houston, TX

Bruce G. Collipp, Marine Engineering NKF Services,

Texas A & M University, Pennsylvania

Barbara A. Shaw, Chairman,

Ovide J. Davis, Pascagoula, Maria Celia Ximenes,

Chevron Shipping Co., San Francisco,

Edison Welding Institute, Columbus, Polymer Sciences, Avon, OH Bethlehem

Steel Corp., Bethlehem, Bethlehem, PA

Bruce R. Somers, Lehigh University,

Ship Structure Committee

Rear Admi al, U.S. Coast Chairman, Ship Structure C{ $//y f; 5

Technical Report Documentation Page

PerForming OrganizationReportNo. SR-1364

10. Work Unit No. (TRAIS)

11. Contmctor Grant No.

13. Type of Reportand PeriodCovered Final 14. SponsoringAgencyCode G-M

18. DistributionStatementD~s t~ibution

Available from: National Technical Information Service Sprirmfield, VA 22161

SecurityClassIf.(of this report) Unclassified

20. SecurityClassification this page) (of Unclassified Reproduction completedpageauthorized of

Form DOT F 1700.7 (8-72)

METRIC CONVERSION CARD

Cofiversions fmm Meiric Measures

Symbol When YouKnow Multiply by

teaspoons ~;p tabl&poons cubic nches i fl oz fluid ounces

centimeters meters meters kilometem

inches inches feet yards miles

~dj ,--~, ---F

muldoli ~] 9/5, degtees Fahrenheit

PART 1 PROJECT 1,0

1-1 1-1 1-1 1-2 1-2 1-3 1-3 1-4 1-4

PART 2 ASSESSMENT 1,0

FOR FINITE ELEMENT