PULS User Manual
PULS User Manual
PULS User Manual
USER MANUAL
PULS
PULS
JANUARY 2009
All rights reserved. No part of this book may be reproduced, in any form or by any means, without
permission in writing from the publisher.
Published by:
Telephone: +47 67 57 76 50
Fax: +47 67 57 72 72
E-mail, sales: [email protected]
E-mail, support: [email protected]
Website: www.dnv.com/software
If any person suffers loss or damage which is proved to have been caused by any negligent act or ommision of Det Norske Veritas, then Det norske Veritas shall pay compensation to such person
for his proved direct loss or damage. However, the compensation shall not exceed an amount equal to ten times the fee charged for the service in question, provided that the maximum
compensation shall never exceed USD 2 million. In this provision "Det Norske Veritas" shall mean the Foundation Det Norske Veritas as well as all its subsidiaries, directors, officers, employees,
agents and any other acting on behalf of Det Norske Veritas.
Nauticus Hull User Manual PULS
ProgramVersion 11.0 / PULS Version 2.0.9 January 2009 1
Contents
1 INTRODUCTION ....................................................................................................... 5
1.1 General ....................................................................................................................... 5
1.2 Approach...................................................................................................................... 7
1.3 Revision history ........................................................................................................... 7
1.4 Buckling of hull elements - problem identifications.................................................. 10
8 REFERENCES ........................................................................................................ 85
1 Introduction
1.1 GENERAL
The present report gives a User’s guide to the Windows program PULS 2.0 with a brief review
of theoretical background, design principles and purpose. Main emphasis is on program features
and functionality, input description, results and example illustration.
PULS is a computerised buckling code for thin-walled plate constructions. It assess the elastic
buckling stresses and ultimate load bearing capacities under combined loads of stiffened and
unstiffened plates used as building blocks in larger plated constructions such as ships and
offshore constructions. The PULS element library is illustrated in Fig.1 below
PULS 2.0
ELEMENT LIBRARY
U3 S3 T1
Unstiffened plate Stiffened plate Stiffened plate
Non-linear theory Regular geometry Non-regular geometry
Non-linear theory Linear theory
The present PULS 2.0 version has three elements, which are briefly summarized in the
following. A detailed description of each element is given in separate chapters. For more
theoretical details see publications given in the reference list at the end of this report.
1.2 Approach
The PULS buckling models for the S3 and U3 elements apply the non-linear large deflection
plate theory of Marguerre and von Karman. Discretizations of the buckling displacements follow
the Rayleigh-Ritz method using Fourier series expansions across the plate and stiffener surfaces.
Energy principles are used for establishing the algebraic non-linear equilibrium equations and
incremental perturbation techniques are used for solving these. Hot spot stress control using the
redistributed membrane stresses in a selection of critical locations determines the ultimate
strength allowing for overcritical strength.
For stiffened panels (S3 element) an orthotropic version of the same theory is used for the global
buckling mode. The local buckling model treats the plate and stiffeners as discrete elements
including all relevant effects such as buckling of the plate between the stiffeners, buckling of the
stiffener web plate as well the rotational restraints between plate and stiffeners. The interaction
between local and global (lateral) stiffener buckling is coped with using modified orthotropic
stiffness coefficients (reduced stiffness/modulus coefficients).
For the T1 element (plates with stiffeners in arbitrary direction), linearized plate theory is used
(not overcritical strength). This means that the buckling capacity predicted for this element is
maximum the minimum eigenvalue. The incremental approach and hot spot stress control is used
also for this element for assessing the strength of geometrical imperfect plates.
The PULS code is programmed in a Windows environment using Visual Basic (VB) tools.
Subroutines are programmed in Fortran 2000.
Two different user-interfaces are available
i) Advanced Viewer (AV). Provides 3D graphical presentation of buckling modes,
membrane stress redistributions, load interaction curves etc.
ii) PULS Excel version with compatible input/output data files with the AV version.
1.3 May New local buckling model -Local eigenvalues for Aluminium
2002 implemented for stiffened panels: pure shear and shear in option locked
S2. comb. with normal ( alu can be
stresses shown. Shown specified in
-More correct local buckling
in: steel mode but
assessment of slender stiffener webs,
no HAZ
rotational restraint effects from ..detail result window,
effects
stiffeners
..in local eigenmode included.)
-Local shear buckling assessment window,
with rotational restraint effects from
..in capacity curves
stiffeners
Stiffened panels
Typically the main load components acting on a local stiffened panel/plate are
i) In-plane load in the direction of the stiffener, compression or tension
ii) In-plane load in the direction perpendicular to the stiffener, compression or tension
iii) In-plane shear
iv) Lateral pressure from sea or cargo
Buckling and ultimate strength of plates depends on the nature of the locally applied loads and
the boundary conditions enforced from the surrounding structure. The boundary conditions can
be categorised in two groups, i.e. out-of-plane support and in-plane support.
Out-of-plane support: In most codes the out-of-plane support along the outer plate edges is
assumed to be rigid in the lateral direction while free to rotate. This corresponds to the classical
simply supported boundary conditions. This is also the default for the present PULS elements,
but alternative boundary conditions may be specified. See also description for each element type
for more details.
In-plane support: The in-plane (membrane) support is also important particularly with respect to
the elements ability to carry loads beyond the elastic buckling load (LEB and GEB eigenvalues)
level, i.e. the nature of in-plane support influences the postbuckling and ultimate capacity
behaviour (UC) of thin plates.
The PULS elements U3 and S3 are by default assumed to be integrated elements, Fig.2.
Integrated elements mean that the plate edges can transfer second order normal membrane
stresses to neighbouring elements, compressive as well as tension (i.e. tension fields). In a
theoretical language this means that the plate edges are constrained to remain straight but free to
move in-plane. This is a constraint that is very much used both in the analytical and numerical
published literature on plate buckling. In practise this constraint means that the present PULS
elements are relevant for dimensioning and strength control of “internal” plates and panels such
as in bottom and deck structures, bulkheads, shipsides etc. Weaker “free to pull in edges”
membrane conditions may be specified as an alternative being more relevant for plate girders etc.
For some type of structures and type of loadings (e.g. serviceability loads, SLS) it may be that
elastic buckling and thereby large deflections are not acceptable. This philosophy implemented
in a design code means that the plate thickness and stiffener proportions are to be increased as
compared to a ULS philosophy. For girders in particular such a SLS philosophy will be very
reasonable and ensure robust designs with extra margins to take “additional redistributed” loads
coming from accepting elastic buckling of surrounding structures, i.e. typically inner and outer
bottom plating in ship hulls.
For all elements the present PULS code gives a buckling control (BS) and ultimate capacity
control (UC) under a given load combination (nominal applied stresses, load control) as
specified by the user. This corresponds to the standard mode of program operation. The results
are presented as elastic buckling (eigenvalues) and ultimate strength nominal capacities and
summarized in a detailed result table. As a single parameter result, the safety margin in the form
of usage factor is given, both related to buckling (BS) and ultimate capacity (UC). The usage
factor provides a measure of the difference between the user specified loads and the
corresponding ultimate capacity (UC) or buckling strength (BS).
The PULS AV program also provides capacity curves under combined loads. The capacity
curves are illustrated in two-dimensional load-spaces. They are to be understood as limit
boundaries covering the load-space selected by the user. They inform about the strength of the
plates in the different load directions and under any load combination. In this mode the usage
factor is not calculated as it is not defined. The term capacity curve is demonstrated in Section
5.7.
2.1 General
The U3 unstiffened plate element applies Marguerre’s non-linear plate theory (geometrical non-
linearity). The elastic buckling (eigenvalues) and non-linear postbuckling problem are solved
using a multiple degree of freedom model in terms of Fourier expansions.
The plate can be subjected to combined load situations, and the numerical procedure scale the
applied loads up to collapse. Stress control criteria describe the onset of material yielding in the
highest loaded position along the plate edges using the redistributed membrane stress
distribution. The redistributed membrane stresses consist of the external applied nominal stresses
added to the second order stress distribution arising due to elastic buckling and due to presence
of geometrical imperfections from production.
The values of the proportionally scaled loads, at the onset of first edge membrane yield, is taken
as representative for the UC values. UC values based on such first yield criteria for thin-walled
designs are close to the real UC values, and on the conservative side. In addition to the UC value,
the ideal elastic buckling stress (eigenvalue) is calculated (Linear Elastic Buckling).
Using default tolerance settings and boundary conditions, and characteristic yield strength as
specified in the rules, the code predicts UC strength values representative for integrated plates in
larger flat plate constructions. Alternatively, the in-plane boundary conditions may be specified
as free, so that an isolated unstiffened plate is represented.
User specified tolerance input for the max amplitude is optionally. The default tolerance
amplitude is specified in terms of a maximum amplitude (delta = s/200). The tolerance shape is
automatically taken to harmonise with the minimum eigenmode with some added trigger modes
for safe UC strength assessment. The latter is called imperfection model and is not possible to
control by the user.
Aspect ratio limit: L1/L2 < 20 for L1 > L2 (or equivalent L2/L1 < 20 for L1 < L2)
Plate slenderness ratio: Li/tp < 200 (Li = minimum of L1 and L2)
By ensuring the maximum membrane stresses along the plate edges to stay below the yield stress
condition (von Mises), excessive permanent sets and buckles are prevented ( principle ii) above).
In plates second order bending stresses adds to the second order membrane stresses. The
resulting surface stress is generally accepted to exceed the yield condition in local areas. This is
not considered critical for the ULS strength of plates and it will not induce permanent sets
beyond what is normally accepted in ship designs.
In some cases, it will from serviceability/functional reasons (SLS), not be acceptable that elastic
buckling takes place. In such cases the local elastic buckling stress (LEB cut off) can be used as
the upper limit of allowable load application. In effect a SLS type of strength assessment allows
no elastic buckling and may prove to be a useful approach in particular for girder webs, webs
frames i fore and aft ship etc. and in general for designs that are not accepted to buckle
elastically nor plastically.
σ1,1 = Λ σ10,1
σ1, 2 = Λ σ10, 2
σ 2,1 = Λ σ 20,1
σ 2,1 = Λ σ 20, 2
σ 3 = Λ σ 30
where Λ is the proportional load factor automatically controlled by the program. A subscript 0 on
the nominal external stresses indicates the values of the input values, i.e. they corresponds to Λ=
1. This correspond to a proportional load history control up to elastic buckling (LEB) and
beyond to final collapse (UC).
Lateral pressure can also be specified. However, it should be noted that by default the plate is
assumed to have simply supported boundary conditions along all four edges, i.e. a type of single
span model in both directions. This is not relevant for analysing cases where the lateral pressure
is acting across several spans (i.e. several stiffener and girder spans) since then symmetric
(clamped) deflection modes needs to be included on top. In order to analyse this case more
realistic a multi span plate model is needed. Multi span unstiffened plate model is not available
in the current PULS version. Alternatively, the boundary conditions should be prescribed so that
all four edges are clamped if this can be documented to be realistic boundary conditions.
In practical cases the stiffeners will be designed to carry the lateral pressure and thus the S3
model will ensure all relevant failure modes.
For strength assessment the lateral pressure p is kept fixed, equal to the input value, while the in-
plane loads are scaled proportionally until elastic buckling and UC strength is identified.
Alternatively, a plate subjected to lateral pressure alone may be analysed. In this case, stress and
deflection results are presented for the specified pressure magnitude, and no UC assessment is
performed.
The considered external loads typically take the form as illustrated on Fig.3.
σ3
σ2,2
L1
x1
x2
σ3 σ2,1
σ1,2 σ1,1
L2
Graphical results available are plots of the minimum eigenmode, the imperfection mode, and the
ultimate limit state deflection mode. Membrane stress plots for the ultimate limit state are also
available. For pure lateral pressure computation, it is possible to request bending stress output,
rather than membrane stress output.
For a given prescribed load combination, the code calculates the value ΛE of the load factor at
ideal elastic buckling. The elastic buckling stresses follow as
σ1E ,1 = Λ E σ10,1
σ1E , 2 = Λ E σ10, 2
σ 2 E ,1 = Λ E σ 20,1
σ 2 E , 2 = Λ E σ 20, 2
σ 3E = Λ E σ 30
σ1u ,1 = Λ u σ10,1
σ1u , 2 = Λ u σ10, 2
σ 2 u ,1 = Λ u σ 20,1
σ 2 u , 2 = Λ u σ 20, 2
σ 3u = Λ u σ 30
(σ1u ,1 ) 2 + (σ1u , 2 ) 2 + (σ 2 u ,1 ) 2 + (σ 2 u , 2 ) 2 + (σ 3u ) 2
η=
(σ10,1 ) 2 + (σ10, 2 ) 2 + (σ 20,1 ) 2 + (σ 20, 2 ) 2 + (σ 30 ) 2
The calculated usage factor η is to be measured against the allowable ηallow (or called ηmax) as
specified in the relevant Rules and Standards.
i) Bi-axial load space σ1-σ2 for constant shear say τ12 = 10, 20, 30 MPa and for fixed p
= 0.0, 0.1 MPa etc.
ii) Axial-shear load space σ1-τ12 for constant transverse compression say σ2 = 10, 20, 30
MPa and for fixed p = 0.0, 0.1 MPa etc.
iii) Transverse-shear load space σ2-τ12 for constant axial compression say σ1 = 10, 20, 30
MPa and for fixed p = 0.0, 0.1 MPa etc.
t
p f = 2σ F ( ) 2
s
which corresponds to first material yielding in extreme fibre along the long edges due to pure
bending stress across the plate thickness.
3.1 General
The S3 element code applies Marguerre’s non-linear plate theory (geometrical non-linearity) in
combination with stress control criteria for ultimate strength assessment.
The stiffened panel may be subjected to combined loads. The model calculates the ideal elastic
buckling loads (eigenvalues) as a separate procedure. A non-linear model is used for assessing
the elastic postbuckling strength, while stress control criteria covers the inelastic response and is
used for determining the ultimate strength. The stress control criteria describes the onset of
material yielding in a selection of critical positions in the panel (hard corners) and are called
limit state functions.
The stresses in the hard corners are calculated as the sum of the direct applied membrane stresses
added to the second order membrane stress due to buckling. The second order membrane stresses
have contributions from the local buckling of the plate between stiffeners-sideways/torsional
buckling of the stiffeners and global buckling of the stiffeners (out-of plate bending of
stiffeners). Membrane stresses are mid-plane stresses of each component plate in the cross-
section. Bending stresses across any component plate thickness are not included in the limit state
yield criteria. The limit states solved explicitly gives the ultimate strength.
Using default tolerance settings and characteristic yield strength as specified in the rules, the
code predicts UC strength values as being representative for integrated plates in larger flat plate
constructions consistently including redistribution of stresses between primary stiffeners and the
plate.
User specified tolerance for the maximum amplitude input is optionally. The tolerance shape is
always (automatically) taken to harmonise with some critical modes being most influential for
the ultimate strength. The latter is called the imperfection model and can not be controlled by the
user.
Validity limits: Plate slenderness: s/tp < 200 (s, shortest edge between primary and sec stiffeners)
Max aspect ratio of plate between prim/sec. stiffeners; 0.17 < L1/s < 20
Primary stiffeners:
Stiffener web slenderness for flat bar stiffeners, h/tw < 35
Stiffener web slenderness Angle and T profiles; h/tw < 90
Free flange slenderness for Angle and T profiles; ff/tf < 15
Minimum width of flange for Angle and T profiles bf/hw > 0.22
Secondary stiffeners(AV):
Stiffener web slenderness for flat bar stiffeners, h/tw < 15
Stiffener web slenderness Angle and T profiles; hw/tw < 33
Free flange slenderness for Angle and T profiles; ff/tf < 15
Geometry: Uni-axially stiffened plate with open continuous profiles welded to plating. The
stiffeners may also be specified as sniped. The stiffeners may have a tilt angle
relative to the plate normal.
Profiles: Open profile type with standard angle shape, eccentric angle, symmetric T or flat
bar shape. Standard profile table included. Profiles run continuously across
several spans. Bulb profiles are modelled as equivalent angle profile.
Sec. stiffeners: Secondary stiffeners perpendicular to continuous stiffeners; option for sniped
stiffeners, continuous stiffeners or tripping brackets.
Loads: Uniform membrane stresses in axial σ10 direction, linearly varying stresses in
transverse σ 20 direction, and constant shear stress σ 30 . All combinations of bi-
axial compression/tension and shear can be analysed. Uniform lateral pressure
acting across several bays is included.
Edge Support: All four edges are assumed to be straight, free to move in-plane and free to rotate
out-of-plane(simply-supported edges) and constrained to be straight and stay in
the original plate plane. This implies full utilisation of normal compressive and
tension fields. Alternatively, two opposite edges may be specified as free to pull
in, for analysis of web girders, stringers etc.
In order not to exceed the range of validity of the theory used in the S3 sub-element, the
following slenderness requirements of component plate elements in a cross-section are
formulated:
i) General:
Web slenderness for flat bar stiffeners: hw / t w < 35
* NB! Note that the “real” validity limit for the plate between stiffeners may well exceed this recommended limit if
the actual load condition analysed enforces a wave pattern sufficiently described by less than 20+1=21 half-waves in
the axial direction (x1-direction). The local minimum eigenmode can give a hint of the “real” validity limit of the
ULS capacity calculated by PULS.
These validity limits are checked by the program. An error message is given when limits are
exceeded. A summary is given on the status bar shown in the bottom of the PULS AV window.
i) Elastic local buckling of any of the component plates in a panel section is accepted.
Local buckling is classified as all modes where the stiffener/plate junction lines act as nodal lines
in the buckling pattern.
For open profiles local elastic buckling means buckling of the plating between stiffeners,
sideways/torsional buckling of stiffeners, stiffener web plate buckling and interactions between
these modes.
iii) Global (overall) buckling of the panel is not accepted (GEB cut-off).
This principle ensures the panel as a whole to have sufficient out-of plane bending stiffness to
avoid global buckling (overall stiffener buckling, see Fig.8 for illustration). Sufficient overall
bending stiffness of the stiffeners ensures lateral support to the component plates, which is a
reasonable requirement for accepting local elastic buckling of these (principle i).
In sum the present ULS philosophy predicts an ultimate strength value accepting elastic local
buckling deflections of plates and stiffeners while preventing excessive permanent damages.
ULS design for ships corresponds to an extreme load condition typically the probably largest
loading experienced in a 20 years period.
These ULS design principles are established in order to constrain the panel designs to have some
minimum stiffness properties for efficient in-plane load transfer. The principles are consistent
with the present DNV rules and guidelines, even though in the latter they are not explicitly stated
and not consistently included for combined load situations.
In some cases, it will from serviceability/functional reasons (SLS), not be acceptable that local
elastic buckling takes place. In such cases the local elastic buckling boundary (LEB and GEB
cut-off) can be used as the upper limit of allowable load application.
The ultimate limit state calculation procedure for the S3 element can be split into three levels:
Open profiles;
T, L, bulb, flat-bar
Profile type: n
The considered external loads typically takes the form as illustrated on Fig.5.
Stiffeners σ3
σ2
L1
x1
x2
σ3 σ2
σ1 σ1
L2
i) and ii). The eigenvalue and the corresponding eigenmode are ideal elastic buckling stresses
with associated buckling shape for a stiffened panel with perfect flat geometry. They are
categorised into local (LEB) and global (GEB) modes. They represent an idealized reference
state at which the panel will start buckling (deflecting out of plane).
iii) The ultimate stress (UC) is the maximum nominal stress the panel can carry for the defined
proportional load history.
iv) The buckling load is defined as the minimum of the eigenvalues and the ultimate load.
Applicable if functional requirements are to be imposed (SLS) with the purpose of avoiding
elastic buckling deflections of plates and stiffeners.
v) The UC usage factor describes the margin between the applied loads and the corresponding
ultimate capacity stresses. The usage factor BS is measured against the buckling load.
vi) Orthotropic macro coefficients represent the in-plane and out-of-plane stiffness of the panel
in an unloaded and loaded state. They are reduced compared to the linear smeared macro
coefficients taking into account the non-linear effect of local elastic buckling of plates and
stiffeners. They can be used as reduced efficiency elements in linear FE models or as reduced
stiffness elements in simple hull girder models.
Graphical results available are plots of minimum eigenmodes, imperfection modes, and ultimate
limit state deflection modes. In addition, membrane stress plots are produced.
The calculated parameters are more thoroughly explained in the following sections.
For all buckling and deflection modes the stiffener flanges are considered locally strong and they
bend and twist like a beam constrained to follow the stiffener web deflections.
Interactions between these mode categories are typical for stiffened panels. Fig.6 illustrates the
three main categories as assessed by PULS
ii) torsional stiffener buckling (typical for tall profiles and flat bar profiles in particular); LEB
iii) stiffener web plate buckling interaction with plate buckling and some torsional stiffener buckling; LEB
The values for the in-plane stresses at the instant where local elastic buckling starts are defined
as
σ1LE = Λ LE σ10
σ 2,1LE = Λ LE σ 20,1
σ 2, 2 LE = Λ LE σ 20, 2
σ3LE = Λ LE σ30
where Λ LE is the eigenvalue of the load proportionality parameter Λ calculated by the program.
The values σ 1LE , σ 2,1LE , σ 2, 2 LE , σ 3 LE are called the local elastic buckling stresses under a
combined load situation.
The local elastic buckling stresses are not critical with respect to the ultimate load bearing
capacity of integrated panels. Higher loads can be carried due to a positive postbuckling effect.
In an ultimate strength context the local buckling stress (LEB) can be viewed as a reference state,
beyond which there exist extra load bearing capacity (typically for thin plates or high local
slenderness in general).
The local elastic buckling stress can be useful as an upper limit for panels constrained to follow
strict functional requirements (SLS), i.e. for design where elastic buckling deflections are not
accepted.
According to non-linear plate theory this macro material law takes the form of an incremental
relation between the in-plane loads ( N 1 , N 2 , N 3 ) and moments (M 1 , M 2 , M 3 ) , and the
corresponding strains (ε1 , ε 2 , ε 3 ) and curvatures ( κ1 , κ 2 , κ 3 ) of the continuous plating. In
mathematical terms the orthotropic macro material law takes the following form
The symbol Δ indicates incremental properties. By accepting local elastic buckling in stiffened
panels under extreme loads, the panel will behave in a non-linear and more flexible way than
under the standard linear response hypothesis. This local panel flexibility is assessed in the
PULS code as a set of reduced orthotropic macro stiffness coefficients defined as
C αβ ≡ C αβ + C αβ
L N
D αβ ≡ D αβ + D αβ α, β = 1, 2, 3
L N
Q αβ ≡ Q αβ + Q αβ
L N
Each of the coefficients in the stiffness matrix has two contributions, i.e. a linear part and a non-
linear part. The linear part represents coefficients with the full stiffener rigidities smeared out
N N N
over the total plate area. The non-linear corrections C αβ , D αβ , Q αβ are assessed using a
numerical procedure. This and more details about the theory is explained in the Theory Manual,
Ref.[3].
⎡N⎤ ⎡ C Q⎤ ⎡ ε ⎤ ⎡ε ⎤ ⎡M S⎤⎡ N ⎤
⎢M ⎥ = ⎢Q T D ⎥⎦ ⎢⎣ κ ⎥⎦
⇔ ⎢ κ ⎥ = ⎢S T F ⎥⎦ ⎢⎣M ⎥⎦
⎣ ⎦ ⎣ ⎣ ⎦ ⎣
The bending/extension coupling coefficients Q αβ , (Q) enter due to the eccentricity effect of the
one-sided welded stiffeners. This coupling effect is eliminated using a simple neutralisation
procedure. The resulting uncoupled moment-curvature relation is
~
M = Dκ
~
[
D = D − (Q) T (C) −1 Q ] (= F −1 )
~
The PULS S3 element calculates the macro model submatrices C, Q, D, D for the linear state, for
the initial unloaded geometrically imperfect state and for an averaged state (secant) being
representative for the ultimate strength assessment.
a) b)
The corresponding global elastic buckling level (GEB eigenvalue) is assessed using the classical
orthotropic plate theory, though as explained in chapter 3.6.3, with modified orthotropic macro
material coefficients accounting for the local buckling effects.
With reference to the DNV rules, Ref.[4], Ref.[5], the lateral buckling mode (or column
buckling) corresponds to the global buckling mode in PULS. However, in the DNV rules, and
most other international codes, the lateral buckling mode neglects the plate effect altogether (i.e.
the column approach assumes lateral support only at two opposite edges typically being at
transverse frame supports). This is the reason while the column approach is not suited for
analysing integrated stiffened panels in ship and offshore structures addressing the problem of
combined in-plane bi-axial and in-plane shear loads.
According to design principle iii), Section 3.2, the present S3 element takes the global
eigenvalue as the upper limit of the buckling capacity of the panel. This is refereed to as GEB
cut-off in the present code and is consistent the Perry approach much referenced in the literature.
In the standard program mode the global eigenvalues are found by scaling the simultaneously
combined loads σ 10 , σ 20,1 , σ 20, 2 , σ 30 in proportion until global buckling takes place. The global
buckling loads are accordingly
σ1GE = Λ GE σ10
σ 2,1GE = Λ GE σ 20,1
σ 2, 2GE = Λ GE σ 20, 2
σ 3GE = Λ GE σ 30
where Λ GE is the global eigenvalue of the load parameter Λ found by the program. The nominal
stresses σ 1GE , σ 2,1GE , σ 2, 2GE , σ 3GE are called the global elastic buckling stresses under a combined
load situation.
(K − ΛK g ) q = 0
K is the small displacement stiffness matrix, Kg is the geometrical stiffness matrix and q is the
eigenvector.
3
1
Fig. 9 Stress control points in critical positions in a panel defining the ultimate limit states
The six limit states fi’s are stress controls in the following positions:
i = 1; Plate criterion: Stress control along plate edges – based on max edge stresses along
supported edges (typical: transverse load when local buckling dominates)
i = 2; Stiffener tension criterion: Stress control in stiffener; at midspan x1 = L1/2 ; in stiffener
flange for global panel deflecting towards stiffener flange, tension criterion - rare for
compressive loads, but kicks in for tension loads (will also kick in for transverse
compressive loads for panel with small stiffeners, i.e. large global effects)
i = 3; Plate compression criterion: Stress control in plate; at midspan x1 = L1/2; in plating
for global panel deflecting towards stiffener flange, compression criterion (PI collapse)
i = 4; Stiffener compression criterion: Stiffener criterion Stress control in stiffener; at
midspan x1 = -L1/2: in stiffener flange for global panel deflecting towards plating,
compression criterion (SI collapse) (typical for pure axial load)
i = 5; Plate tension criterion: Stress control in plate; at midspan x1 = -L1/2; in plating for
global panel deflecting towards plating, tension criterion – rare for compressive loads,
but kicks inn for tension loads
(Note that the limit state criteria i = 2-5 is not always evaluated at midspan. Maximum
curvature in x1 – direction, and thereby the highest bending stress, could be closer to the
ends for certain geometrical proportions of stiffened panels. Typical are cases with small
stiffeners for which the panel behaves more as a plate than a “column”, with a global
buckling mode pattern flattening in the mid-regions.)
i = 6; Stiffener bending stress criterion at support: Stress and capacity control at support x1
= 0; compressive or tension criterion, kicks in for cases with lateral pressure. This limit
state is used to control the bending and shear capacity of the stiffeners under the
influence of combined lateral load and in-plane loads. Yielding in the stiffener flange at
the transverse frames is accepted, since stiffeners have significant strength reserves after
first yield when subjected to lateral pressure. The panel is loaded until the plastic capacity
of the stiffeners is reached. Two criteria are used for this limit state. The first is the
capacity of the top and bottom flanges to carry the combined axial force and bending
moment resulting from the applied loads, and the second is the capacity of the web to
carry the shear force and axial force due to the applied loads.
The value for the four independent in-plane stresses at the point of ultimate strength is defined as
σ1u = Λ u σ10
σ 2,1u = Λ u σ 20,1
σ 2, 2 u = Λ u σ 20, 2
σ 3u = Λ u σ 30
The value of the Λ u factor is calculated as the minimum explicit solution of the six limit states
functions. The solutions are found using a numerical procedure. There exist one Λ u for each
limit state and the minimum of these is used as representative for the ultimate strength. It follows
per definition that it is the inverse of the UC usage factor η as defined in Chapter 3.9.
3.10.2 STEEL
The steel limit states follow the six main stress control points as given above.
Different material yield stresses can be specified in the plate and stiffeners (same in all
stiffeners) while only one set of values for the Young’s modulus E and Poisson ratio ν is
possible.
The geometrical imperfection size and shape effect, and the residual stress effects, are implicitly
considered in the ULS values when using default tolerance values. The imperfection values used
as default values are typical for welded and fabricated steel plates used in ships and offshore
constructions.
3.10.3 ALUMINIUM
A special option for analysis of aluminium panels, including HAZ effects for NV 5082 T6, is
implemented in PULS 2.0.
A first membrane yield criterion in the HAZ zone is added to the limit states used for steel. This
first HAZ yield criterion will be conservative, but reasonable as a design limitation since limited
knowledge is available for soft HAZ zones strained beyond the material yield level. Crack
initiations and fracture can be the result if the HAZ zones around welds frequently are loaded
beyond the first yield level.
L u = (σ1u + σ 2 u + σ 3u )
2 2 2
L 0 = (σ10 + σ 20 + σ 30 )
2 2 2
Applied loads
L0
LU
Usage factor:
η= L0/LU
Fig. 10 Definition of safety margin/usage factor; capacity curve example for bi-axial loading
For the single load cases the definition of usage factor becomes
η = σ i 0 / σ iu i = 1, 2, 3
as is the familiar form used in the DNV steel ship rules. For the present PULS code approach the
final rule acceptance criterion will be in the form
L0
< η allow i.e. η < ηallow
Lu
where η allow (or also called η max ) is to be specified in the rules and η is calculated by the PULS
code.
In the offshore marked the LRFD format is used, i.e. the acceptance criterion is on the form
Sd < R d
where Sd is the load effect and Rd is the design resistance. The design resistance is related to the
characteristic resistance as
Rd = Rk / γm
For offshore application of the PULS code, the following definitions applies:
Characteristic resistance:
R k = (σ1u + σ 2 u + σ 3u ) (= L u )
2 2 2
Load effect:
S d = (σ10 + σ 20 + σ 30 ) (= L 0 )
2 2 2
i) Bi-axial load space σ1- σ2 for fixed shear σ3 (τ12) and fixed p (MPa)
ii) Axial-shear load space σ1- τ12 for fixed transverse stress σ2 and fixed p (MPa).
iii) Transverse-shear load space σ2 - τ12 for fixed axial stress σ1 and fixed p (MPa)
Limits for acceptable lateral pressures are set in the program and capacity curves generation will
be aborted if these limits are exceeded. The pressure limits are summarised in the detailed result
menu under standard parameters.
SNx1 Sx
N1
p
Fig. 11. Deflection of a stiffener in the simply supported mode (SS asymmetric modes, top)
and the clamped mode (CS symmetric modes, bottom)
Maximum pressure criteria are introduced representing practical design boundaries for the
stiffeners. They kick in and present themselves in dialog boxes with stated limits not to be
exceeded. Two different pressure limits are specified; one for controlling the bending stiffness of
the stiffener/plate unit and one for controlling the shear area of the stiffener web. They are based
on linear beam theory.
The bending stiffness pressure limit pFs is
Wmin
p Fs = 12 σ F 2
s L1
where Wmin is the minimum section modulus of the stiffener/plate unit normally being at the
stiffener flange position
I
Wmin =
h w + 0.5 t f + 0.5 t p − z g
and zg is the neutral axis measured from the plate middle-plane. This limit correspond to first
bending stress yield at support for a clamped stiffener with span L1. The pFs value is tabulated in
the PULS output detailed result/standard parameter list for info.
The lateral pressure which gives shear stress yielding in the web is
2V
p τs =
sL1
σ Fs t w I
V =
3S p
where the moment of area about the neutral axis, where the maximum shear stress occurs, is
S p = st p z g + 0.5t w ( z g − 0.5t p ) 2
This pressure limit correspond to first pure shear yield in the stiffener web as calculated for a
clamped stiffener with span L1. The pτs value is tabulated in the PULS output detailed
result/standard parameter list for info.
Another pressure limit of practical interest is the first onset of surface yield due to pure local
bending of the plate between stiffeners. This limit is simply
t
p F = 2σ F ( ) 2
s
and is based on a clamped plate unit strip formulation. This limit is not included as design limit
in the present PULS 2.0 version as it will be too strict for practical applications. For ship bottom
designs, pressure level significantly beyond this pF is normal. It is also documented using
advanced elasto-plastic FE analysis, that pressure levels well beyond pF can be carried for plates
supported by typical solid stiffener/frame structures. It is also well documented that the in-plane
uls capacity for stiffened panels is not very influenced by local surface yielding along
stiffener/plate supports. The pF value is tabulated in PULS output detailed result/standard
parameter list for info.
Conclusion: The present S3 lateral model has practical pressure limits implemented
reflecting reasonable design requirements to the stiffener/plate unit. These are based on the
philosophy that the stiffeners shall carry the pressure and transfer these to the supporting
transverse frame/bulkhead structures without significant plastic deformations. Input
pressures beyond these limits are not accepted by the program and warnings are given in
dialog boxes.
4.1 General
The T1 stiffened plate element applies a linearized version of Marguerre’s non-linear plate
theory. Linearized in this context is means that the theory do not assess the non-linear
postbuckling strength (no overcritical strength), i.e. loads beyond the ideal elastic buckling level
(eigenvalue) is not possible. The elastic buckling (eigenvalue) is found using a multiple degree
of freedom model in terms of Fourier expansions of lateral deflections. Stress control criteria
checking the onset of material yielding in hot spot stress locations along the plate edges and in
the stiffeners are used as method for the buckling strength assessment.
The stiffened plate can be subjected to combined load situations and the numerical procedures
implemented scale the applied loads up to buckling. The results are presented in terms of a single
parameter called the usage factor. The numerical algorithms are based on the linearized load-
deflection solution using the minimum eigenvalue in the amplification factor.
The element can be used for assessing the buckling strength of panels where the stiffener can
have variable geometry and arbitrary orientations. This means it is suitable for non-regular
stiffener arrangement and triangular plates.
More details concerning the T1 element can be found in ref.(18).
Model - Automatically set default model tolerances consistent with as welded panels
imperfections:
Output: - Buckling Strength (BS)
- Usage factor BS (ratio: applied loads/ BS)
- Minimum Eigenvalue
The theory for stiffeners is simplified in sense that they are not considered to carry any axial
stress. This means that buckling of stiffener web is not checked and thus a certain max
slenderness of the stiffener web has to be ensured. The T1 stiffener model is similar to the
secondary stiffener option in the S3 element, and the web proportions should fulfil the table 2
requirements at least if the stiffeners are continuous across the boundaries. For sniped stiffeners
the stiffener web requirement can be relaxed. As an upper max slenderness limit the following
ratios are implemented:
i) General:
Web slenderness for flat bar stiffeners: hw / t w < 35
Web slenderness for L or T profiles: h / t w < 90
Free flange for L or T profiles: f f / t f < 15
These validity limits are checked by the program. An error message is given when limits are
exceeded. A summary is given on the status bar shown in the bottom of the PULS window.
Another alternative option for modelling stiffeners is also available using directly a user defined
effective moment of inertia and rotational spring stiffness. This option enables the user to specify
his own preferred stiffener efficiency properties.
By ensuring the maximum membrane stresses along the plate edges and the plate/stiffener
connections to stay below the yield stress condition (von Mises), excessive permanent sets and
buckles are prevented ( principle ii) above). In plates second order bending stresses adds to the
second order membrane stresses. The resulting surface stress is generally accepted to exceed the
yield condition in local areas. This is not considered critical for the strength of plates and it will
not induce permanent sets beyond what is normally accepted in ship designs.
Since linear theory is used, the ULS strength calculated will always be below the local elastic
buckling load.
where Λ is the load factor automatically controlled by the program. A subscript 0 on the nominal
external stresses indicates the values of the input values, i.e. they corresponds to Λ=1.
i) The eigenvalue and the corresponding eigenmode are ideal elastic buckling stresses with
associated buckling shape for a stiffened panel with perfect flat geometry. They represent an
idealized reference state at which the panel will start buckling (deflecting out of plane).
ii) The buckling load is defined as the load at the onset of yield. Since linear theory is used, the
buckling load is always lower than the eigenvalue.
iii) The usage factor describes the margin between the applied loads and the corresponding
buckling load.
Graphical results available are plots of the minimum eigenmode, the model imperfection, and the
ultimate limit state deflection mode. In addition, membrane stress plots are produced.
(σ1u ,1 ) 2 + (σ1u , 2 ) 2 + (σ 2 u ,1 ) 2 + (σ 2 u , 2 ) 2 + (σ 3u ) 2
η=
(σ10,1 ) 2 + (σ10, 2 ) 2 + (σ 20,1 ) 2 + (σ 20, 2 ) 2 + (σ 30 ) 2
The calculated usage factor η is to be measured against the allowable ηallow (or called ηmax) as
specified in the relevant Rules and Standards.
The U3 element is assumed to have uniform plate thickness. The S3 element is assumed to have
constant plate thickness and stiffener proportions across the panel. The T1 element is assumed to
have uniform plate thickness, but may have stiffeners with different proportions.
The PULS code concept isolates a rectangular panel and prescribes a set of external loads across
the panel surface. For the S3 element the prescribed external nominal loads (stresses) have to be
uniformly distributed in the axial direction, but may be linearly varying in the transverse
direction.
σ30
σ20
Fig. 12 Rectangular panel; definition of supported edges and applied loads
Limitations and fields of applications of the PULS 2.0 code are summarised in Table 1 in Section
5.2. Usage hints are given in Chapter 6.
5.2 Applications
For overview, typical ship hull areas of application for the U3 and S3 elements are tabulated
below.
Table 1. PULS application matrix. ULS OK means that local elastic buckling is accepted; i.e. a GEB cut-off
is implemented ensuring sufficient stiffener proportions, but plate between stiffeners can buckle elastically.
SLS OK means that not any type of elastic buckling is accepted; i.e. both GEB and LEB cut-off are
implemented.
The PULS program has a multi document (MDI) layout for easy input and output. All
parameters are given in cells with text identification and unit assignment.
The program has two levels of organising data. The highest level is the “Project” typically being
a specific ship identification or similar. The level below is the “Panel” level typically being a
specific panel in the ship. Any number of “Panels” can be specified, but it is recommended not
to exceed 15 especially when there are generated graphics for most panels.
The “Project”, with the corresponding “Panels” organised as subfolders, are saved in a file with
a format notation of *.pbp. The Project name will automatically be updated to be the same as the
name given to the file using the Save as… option in the File menu. The Panel can be given any
name by directly editing in the Tree view control menu or using the identification text box.
The program features and the project control is facilitated using the
Menubar, Toolbar, Tree view control menu and the Window
Delete panel
Save project as…; any name can be used on file, preferred equal to project
Print Word report or data Windows; directly or to file Recalculate all panels
Edit curtain:
N/A
N/A
Result curtain:
Window curtain:
New panel; same as “new Open project; same as Direct button for
panel” from File menu “open project” from File Capacity curves
Save file; same as Print report in Word format to Solve Buttons for
“Save project” on printer or to file. Input button; 3D graphic
File menu parameters and main results same as operations.
“Solve Alternative:
Panel” in right hand
Panel menu, mouse click
requires when in 3D
values on graphic
input loads! window
The tree view control menu operates equal as to Windows Explorer. It has the following layout
Explicit geometrical model tolerances can be specified. For the unstiffened element (U3) there is
only one model tolerance amplitude available, i.e. the maximum initial plate amplitude δp0. For
the stiffened panel (S3) three independent imperfection amplitudes can be specified, i.e. see
figure below; δP0, δS0, and δT0. Default model tolerances correspond to ULS code strength values
according to normal fabrication standards of welded integrated structures used in the
shipbuilding and offshore industry.
The model imperfection shape is not operated by the user, only the imperfection amplitude.
Imperfection models combining eigenmodes for constructing relevant imperfection modes are
pre-set in the program. They can be viewed in the tree control menu after a calculation has been
carried out.
Both amplitude and shape is defined by the program for the T1 element.
Opening an existing file is done by clicking the Open button on the toolbar or by using the Open
option in the File menu.
A new blank input page is generated by clicking the New file button on the tool bar or using the
New option in the File menu.
After an analysis is run the results are saved by clicking the Save button or using the Save or
Save as… in the File menu. Save as…will provide a dialog box for setting of file name, which is
automatically set to be the same as the “Project” name.
The “Panel” identification names can be edited and given any name by clicking on the
corresponding labels in the Tree view control menu.
The Input Window for stiffened panels takes the following form
Input;
Allowable
usage
factor is
defined in
the rules
Output/result:
Calculated actual
usage factor
Input: Plate dimensions and number Tab strips for material, load and
of continuous stiffeners Ns. boundary condition specification-
Secondary transverse stiffener option: see next page
own dialog if different from zero.
Profile table for flat bars, bulbs, angles and T-bars (same as in DNV Nauticus Hull)
Input: Material
parameters.
NB! G not input.
Only pressure
computation: Linear
computation of
deflections for
specified pressure. No
ULS computation.
NB!! Note that all input numbers entered has to be followed by a tab entry for proper registration.
A tab reminder is given on the status bar at bottom whenever a new number is entered into a cell.
Unstiffened plate figure - automatically adjusted with correct aspect ratio proportion
Input;
Allowable
usage
factor
defined in
the rules
Specification of
material properties.
Aluminium
properties also
available.
Specification of
applied loads.
NB!! Note that all input numbers entered has to be followed by a tab entry for proper registration.
A tab reminder is given on the status bar at bottom whenever a new number is entered into a cell.
Change between
moment of inertia
input or stiffener
scantlings input.
Defining the end Buttons for defining and Profile data for
points of each deleting stiffeners. each stiffener
stiffener. Stiffeners can be defined or individually.
deleted with buttons or
directly in picture.
Material input.
Rotational spring
support for plate from
surrounding structure
For pre-stress selection of the “ULS collection” in the dialog box (see next page) will put all
ULS curves into one diagram. The S3 and U3 elements assess capacity curves for preloading
w.r.t. a set of specified lateral pressure levels and stresses. In girder mode the prestress option is
limited. The T1 element supports prestresses, but not lateral pressure.
The dialog box for capacity curves is given under Capacity Curves on the menu (or separate
button) and it takes the form
2D load space
alternatives
Fixed/prestress loads:
Generate curves
For the S3 element the Curve settings for the 2D curves are given in the following dialog box
(for the U3 and T1 elements the dialog box is similar but no distinction between global and local
elastic buckling is necessary)
A set of capacity curves in the 2D load space (according to i) classification above) is shown in
Fig.13 for a constructed S3 element example. The red curve indicates the global elastic buckling
boundary (GEB), which always will be the very outer curve according to the adopted design
principles (see chap.3.3.. design principle iii)). The blue curve indicates local elastic buckling
(LEB) and the violet curve is the ULS curve. Note that the local elastic buckling boundary (LEB)
is below the ULS boundary in a limited region in load space for this example. This is typical for
cases where the plate is thin. For thicker plates the ULS boundary curve will be the inner curve
in the whole load space. See Fig.14 for illustration of such a case.
Fig. 13 Capacity Curves in 2D load space; bi-axial loading, all quadrants, thin plate
Fig. 14 Capacity Curves in 2D load space; bi-axial loading, all quadrants, thick plate
Similar capacity curves are available in the other two load spaces where shear stresses are one of
the load components. Fig. 15 shows an example including shear. Note the symmetry in the shear
load τ12 (= σ3).
Fig. 15 Capacity Curves in 2D load space; transverse load σ2 and shear σ3(τ12), all quadrants.
Fig. 16 Capacity Curves in 2D load space for U3 element (unstiffened plate); axial stress σ1 and shear stress
σ3(τ12), all quadrants.
An example of 2D Capacity curves for the U3 element is given in Fig.16 for combined axial
nominal stress σ1 and nominal shear stress σ3 (= τ12). As can be seen, the example is for
relatively thin plate for which the UC strength is beyond the elastic buckling stress for most of
the load combinations.
In the standard mode of program operation 3D plots of the eigenmode, model imperfections,
pressure plots, ULS mode and ULS membrane stresses is available. These data can be found
under the Result menu.
a)
b)
c)
In particular, the ULS membrane stresses show the non-linear redistributed stress distribution at
the first edge membrane yield condition. This is then, as per definition in the PULS approach, the
ULS state. Some more capacity may be expected for thin plates involving spread of plasticity,
and consequently involving permanent sets, but the extra capacity is marginal.
Fig. 18 Capacity Curves in 2D load space with associated 3D visualization of ULS modes for two different
points in load space
Detail result summary sheet: In the standard mode a detailed results window is available by
selecting Detailed results from the Result menu. The output parameters are commented in figure
below.
Ultimate loads with Minimum local
GEB cut-off = do Buckling loads, eigenvalue; LEB
Applied loads, not accept global GEB and LEB cut-
mapped from stiffener buckling off Minimum global
input window eigenvalue; GEB
Real ULS; accepts
global stiffener
buckling
Calculated
usage factor
Tree
view
control
menu
Ultimate stress and usage factor for each Table showing the maximum
limit state, minimum ultimate stress displacement selected locations
given in detail result summary above. in panel.
Standard Parameters sheet: The stiffened panel element has also a Standard parameter list. Most
of these parameters are not a part of the PULS buckling models and they are only meant as
general user info.
They are categorised into
- Cross-sectional data; cross sectional data for plate and stiffener, x1 = const.
- Cross-sectional bending data; valid for one stiffener/plate unit, full plate width s is included
- Cross-sectional capacity parameters; standard moment and shear load capacities
- Numerical parameters; number of increments in PULS procedure
- Design lateral pressure parameters; Stresses in single stiffener/plate unit according to linear
beam theory for the prescribed pressure
- Maximum lateral pressure parameters - one stiffener/plate unit; pressure limits according to
linear beam theory for a single plate/stiffener unit
- Maximum lateral pressure parameters - plate; reference pressure limit pF based on first edge
bending stress limit for a clamped plate unit with length s.
In the standard mode a Detail results Window is available by selecting Detailed results from the
Result menu. The output parameters are commented on figure below.
Press Tab key - a new number has entered a cell, tab to next field to register the
number.
Recalculation necessary A new set of data has been registered, recalculation is necessary
for updated buckling and ULS strength.
The usage hints are given here for overview. It refers both to program operation and how to
address special problems and applications cases.
General:
i) User specified model tolerances can be specified for the plate and stiffener (U3 and S3),
see under Panel menu on tool bar. In current version only model imperfection amplitudes
not imperfection forms can be specified. Default settings correspond to normal
fabrication standards of steel structures.
ii) Excel Report can be printed, see under File menu. .
Always remember to press the solve button before printing in order to avoid mismatch between
input data and results. This is controlled also by the status bar information at the bottom of the
PULS Windows.
Stretching/compressing axes: Press the shift key on the keyboard and drag the mouse with left
button pressed.
Moving origin: Press the shift key down and drag the mouse with the right button pressed.
iii) The text box in the capacity curve figures can be resized and moved in the window
iv) Graphical visualisation of buckling (or ULS mode) mode from any point in load space laying
on the capacity curves.
Procedure: A point in load space is selected by first select the whole curve (all points visualised
by coloured marks) and then click once more on the point of interest (only this point is to be
marked by colour symbol). Then double click on the selected point and the buckling mode
window pops up. The corresponding co-ordinates in load space are given in the inserted frame
above (always given in MPa). Several points may be “called up” and moved around inside the
window for comparison purposes, illustrating different type of failure modes as function of the
load combination.
U3 element:
i) Detailed results summary
The applied nominal stresses are listed together with the ultimate stresses, buckling loads, the
eigenvalues and the usage factor.
The results are based on a proportional load path for the in-plane stresses meaning the same ratio
between each nominal stress component; applied, at the point of elastic buckling (eigenvalues)
and ultimate state. The lateral pressure is as specified in the input, i.e. it is kept fixed.
ii) Eigenmode
The eigenmode shows the minimum ideal elastic buckling mode according to linear elastic
buckling theory. The deflections are normalised to unity.
iv) Pressure
The pressure plot shows the displacements in the plate caused by the applied pressure. The plot
generated by a regular solution is non-linear while the plot generated by the only pressure mode
is linear. The shown stresses are bending stresses or membrane stresses depending on the users
choice in the load tab strip.
v) UC mode
Shows the deflected form in the UC state.
The UC deflection shown in the figure is the total deflection (wtot) including the geometrical
imperfection and the deflections due to pressure; wtot = wimp + wload + wpressure.
By clicking on the UC figure a yellow frame pops up showing the co-ordinates of the point and
the corresponding UC out-of-plane deflection in [mm].
S3 element:
i) Detailed results summary
The applied nominal stresses are listed together with the ultimate stresses, the eigenvalues in
local and global mode and the usage factor.
The results are based on a proportional load path for the in-plane stresses meaning the same ratio
between each nominal stress component; applied, at the point of elastic buckling (eigenvalues)
and ultimate state. The lateral pressure is as specified in the input, i.e. it is kept fixed.
The results for each of the five limit states are listed. The critical limit state is picked out and
given in the detailed summary list.
A set of orthotropic stiffness coefficients are given indicating the degree of reduced stiffness
(increased flexibility) of the panel due to local buckling from the plate between stiffener and
sideways buckling of stiffeners. The stiffness coefficients given in PULS 1.5 represent linear
values (sub/superscript L) and tangent values for zero load (sub/superscript I) and at the secant
values as representative for the ultimate state (sub/superscript U). A superscript N symbolises
neutral relevant for the bending stiffness about the instaneous neutral axis of plate and stiffener
combination.
ii) Eigenmodes
Shows the minimum ideal elastic buckling mode according to linear elastic buckling theory in
the local and global mode separately.
The global eigenvalue is calculated based on an orthotropic plate theory with the reduced
bending stiffness coefficients in the ultimate state.
The eigenmodes are normalised to unity.
iv) Pressure
The pressure plot shows the displacements in the plate caused by the applied pressure. Both the
solution for regular and only pressure mode is linear.
iv) UC mode
Shows the deflected form in the UC state separated into local and global modes.
The ULS deflection shown in the figure is the total deflection (wtot) including the geometrical
imperfection and the displacements due to pressure; wtot = wimp + wload + wpressure.
v) UC membrane stresses
The membrane stresses is shown at the local UC state. This is not exactly as the same as the real
UC state, but for qualitative illustration of non-linearly redistributed stresses the 3D graphics
may be useful. The graphical presentation is only an added feature, it has no influence on the
UC results which consider interaction between local and global buckling effects.
7.1 Introduction
The PULS Excel spreadsheet and the PULS GUI use the same computational routines. The
spreadsheet offers easy input of a large number of panels, and therefore makes parameter studies
easy to perform. The spreadsheet is organized in input sheets and output sheets for the S3 and U3
elements. The input and output parameters are the same as in the Puls GUI. The Puls spreadsheet
is able to read and write pbp-files which are compatible with the Puls GUI.
as corrugated)
Numerical values can also be
specified.
2) Delete/save old results: Delete or save the results already written to the output sheet.
8 References
1 Brush and Almroth. “Buckling of Bars, Plates and Shells”, McGraw-Hill 1975
4 DNV, Rules for Classification of Ship, Part 3, Chapter 1, Sec.14, January 2000
5 DNV, Classification Note No.30.1, July 95
6 Steen, E., Østvold, T.K and Valsgård, S. “A new design model for ultimate and
buckling strength assessment of stiffened plates”, PRADS 2001, Shanghai, China.
7 Byklum, E. and Amdahl, J. “Nonlinear Buckling Analysis and Ultimate Strength
Prediction of Stiffened Steel and Aluminium Panels”, The Second International
Conference on Advances in Structural Engineering and Mechanics,
Pusan, 21-23 August 2002.
8 Steen, E. “Application of the perturbation method to plate buckling problems”,
Research Report in Mechanics, No. 98-1, ISBN 82-553-1149-1, University of Oslo,
Department of Mathematics, Mechanics Division, 1998.
9 Steen, E. “Buckling of Stiffened Plates using a Shanley Model Approach”, Research
Report in Mechanics, No. 99-1, ISBN 82-553-1182-3, University of Oslo, Department
of Mathematics, Mechanics Division, 1999.
10 Steen, E and Østvold, T. K. “Basis for a new buckling model for strength assessment of
stiffened panels”, DNV Seminar: Buckling and Ultimate Strength of Ship structures, 20
September 2000, Høvik, Norway.
11 Østvold, T. K and Steen, E. “Windows program implementation of new buckling
code”, DNV Seminar: Buckling and Ultimate Strength of Ship structures, 20 September
2000, Høvik, Norway.
12 Steen, E and Østvold, T. K. “Simplified treatment of buckling effects in global linear
FE ship models using an orthotropic macro material modelling technique”, DNV
Seminar: Buckling and Ultimate Strength of Ship structures, 20 September 2000.
Høvik, Norway.
13 Byklum, E. and Amdahl, J. “A simplified method for elastic large deflection analysis
of plates and stiffened panels due to local buckling”, Thin Walled Structures, vol. 40,
no. 11, pp. 923-951, 2002.
14 Byklum, E. “Ultimate strength analysis of stiffened steel and aluminium panels using
semi-analytical methods", Dr.ing. thesis, Department of Marine Technology,
Norwegian University of Science and Technology, 2002.
15 Byklum, E., Steen, E. and Amdahl, J. “A semi-analytical for global buckling and
postbuckling analysis of stiffened plates”, Thin Walled Structures, vol. 42, no. 5, pp.
701-717, 2004.
16 Steen, E. , Byklum, E. and Vilming, Kjetil G. “Computer efficient non-linear buckling
models for capacity assessments of stiffened panels subjected to combined loads”,
ICTWS 2204, Fourth International Conference on Thin-Walled structures, 22-24 June,
Loughorough University, UK.
17 Steen, E. , Byklum, E., Vilming, K. G. and Østvold, T. K.“Computerized Buckling
Models for Ultimate Strength Assessment of Stiffened Ship Hull Panels”, PRADS
2204, Lubeck/Travemunde, Germany
18 Brubak, L., Hellesland, J., Steen, E., Byklum, E. “Aprroximate buckling strength
analysis of plates with arbitrarly oriented stiffeners”, NSCM 2004, Nordic Seminar on
Computational Mechanics, KTH, Stockholm, Sweden
NOTATION
σ1 - Nominal uniform stress in stiffener direction, i.e. along x1 – axis, compression positive –
tension negative, unit N/mm2
σ2 - Nominal uniform stress in perpendicular to stiffener direction, i.e. along x2 – axis,
compression positive – tension negative, unit N/mm2
σ3 - Nominal uniform in-plane shear stress, positive/negative, unit N/mm2
τ12 Shear stress; same as σ 3
Λ Load proportionality factor, unit load factor (Λ=1) corresponds to applied loads
σ10 , σ 20 , σ 30
ΛE Load proportionality factor at ideal elastic buckling for U3, eigenvalue
ΛGE Load prop. factor at ideal elastic buckling in global mode, global eigenvalue GEB
ΛLE Load prop. factor at ideal elastic buckling in local mode, local eigenvalue LEB
Λu Load proportionality factor at ideal ultimate load
L0 Load effect (= Sd)
Lu Characteristic resistance (= Rk)
Sd Load effect (= L0)
Rk Characteristic resistance (= Lu)
γm Material factor, offshore rule notation
Superscripts or subscripts
L Linear properties
N Non-linear properties
U Ultimate limit state
I Initial properties for zero load
T Transposed
0 Applied nominal stresses/loads, reference nominal stresses/loads
- o0o -