Materials and Cross Sections: Sofistik 2025
Materials and Cross Sections: Sofistik 2025
SOFiSTiK | 2025
AQUA
Materials and Cross Sections
SOFiSTiK AG
This manual is protected by copyright laws. No part of it may be translated, copied or reproduced, in any form or by
any means, without written permission from SOFiSTiK AG. SOFiSTiK reserves the right to modify or to release
new editions of this manual.
The manual and the program have been thoroughly checked for errors. However, SOFiSTiK does not claim that
either one is completely error free. Errors and omissions are corrected as soon as they are detected.
The user of the program is solely responsible for the applications. We strongly encourage the user to test the
correctness of all calculations at least by random sampling.
Contents | AQUA
Contents
Contents i
1 General 1-1
1.1 Task Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-1
1.2 Types of sections . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-1
1.2.1 Static Properties of Cross Sections . . . . . . . . . . . . . . . . . 1-1
1.2.2 Standard Cross Sections . . . . . . . . . . . . . . . . . . . . . . . 1-1
1.2.3 Freely Defined Thin-walled Cross Sections . . . . . . . . . . . . 1-1
1.2.4 Freely Defined Solid Cross Sections . . . . . . . . . . . . . . . . 1-2
1.2.5 Freely Defined FE Cross Sections . . . . . . . . . . . . . . . . . . 1-2
1.2.6 Selection of Section Type . . . . . . . . . . . . . . . . . . . . . . . 1-2
1.3 Creating variants of sections . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1-2
SOFiSTiK 2025 i
AQUA | Contents
ii SOFiSTiK 2025
Contents | AQUA
Bibliography 4-21
iv SOFiSTiK 2025
General | AQUA
1 General
After definition with AQUA, the cross sections can be represented graphically with Re-
sultViewer.
Thus all other sections do not allow neither the calculation of von Mises stresses nor the design
for fatigue.
Section moduli for all stresses are available at all points of the cross section. Torsional moment
of inertia and warping resistance, as well as centre of shear and shear deformation areas, are
determined for open or closed shapes, but they can also be specified explicitly for special
cases. Composite cross sections can be defined.
Section moduli for all stresses are only available at distinct points of the cross section. The
torsional moment of inertia, the centre of shear and the shear deformation areas can be calcu-
lated, or they can be input separately. The warping resistance can not be determined. Com-
posite sections or effective widths of the polygons can be defined.
All sectional values including warping are evaluated. Section moduli for all forces and moments
are available for all element mid points. Composite sections may be defined.
A thin walled section has much in common with a standard frame analysis, while the solid
section requires a continua solution. This means that the simplifications of the thin walled
approach allow a faster, more robust and more extensive solution, but neglects local effects,
which may become visible within a continua solution with integral equations or finite elements.
For example it is to be noted, that for a thin hollow box, the shear stress of the continua solution
is not really constant across the web thickness and may have higher intensities at the corners
which might require a smoothing of the contour. On the other side a thin walled section has
problems to take account of the positive effect of the fillets of a thin rolled steel shape which
has a considerable contribution to the torsional strength.
The checks for the c/t ratio for steel sections are much more easily performed with a thin walled
section and the modelling of discrete dowels is only possible with this type of approach.
In general problems have to be expected if a section is not modelled with the optimum method.
Especially very thin plates with stiffeners modelled as polygons need a very high numerical
effort.
• You may describe the section via CADINP variables within a block, which is then used
multiple times.
• You may interpolate between two sections linearly
• You can define a cross section template consisting of several construction points. Other
points are referenced hierarchically to those original points. You can then generate other
cross sections by changing these points.
• You may describe the position of those construction points by a 3D model with curved
reference lines.
It is also possible for AQUA to update all interpolated or otherwise generated sections with a
single command (INTE).
2 Theoretical Principles
2.1 Materials
Properties of materials must be distinguished according to whether they are to be kept as close
as possible to real values (e.g. for dynamic calculations) or to be used with a safety coefficient
for calculating an ultimate load-bearing capacity.
A small, but subtle contradiction is given by the fact, that many design codes use a factor of
10.0 to convert the density of a material t/ m3 to the weight kN/ m3 . The US customary
units use the pound as base unit, the masses are converted with a nominal value of the gravity
acceleration.
To avoid any confusion and to use all the benefits of the SI-System SOFiSTiK has established
the following rules:
• All specified weights have to be specified in (kN/ m3) or (b/ ƒ t3), but they are saved as
masses to the data base. Thus a direct definition of the mass with a unit of (kg/ m3) or
(bm/ ƒ t3) is also possible.
• Masses for dynamic or densities for thermal analysis will be always determined directly
from that data.
• If the design code or the user allows for more precise loadings, the user may specify the
true value of the gravity acceleration for the reference location with the NORM record.
Somehow confusing, there is an additional density value of the raw material. This value is
relevant for many physical properties of concrete or timber, but its value is less than the value
used for self weight.
Whereas the safety factors were formerly assigned more-or-less at random, sometimes to the
load and sometimes to the material, more recent regulations (Eurocode) provide a clearer
separation between safety factors for the loads and factors for the material.
Since the material safety factors still depend on the nature of the load or the type of design,
AQUA generates and stores only the genuine properties of the material. However, AQUA
accounts for some safety factors which are independent of the particular loading case, such as
long term reduction factors.
Nevertheless, a safety coefficient can be entered in AQUA for each material; this is used in
AQUA for calculating the full plastic section forces and moments, and can be used in AQB for
the strain checks.
• A polygon/circle may create a hole in another area only if it has a different material number
and at least one vertex of the periphery within the other area. The latter is only true provided
that the edges of the two polygons maintain a minimum distance of 0.5 [ mm] apart.
• If a polygon/circle is completely within another area, it will always create a hole, if the areas
overlap partly the sequence of the input will decide: The area defined last will be considered
”in front of” and will create a hole in the area laying behind.
The coordinate system of the section is identical with the local beam coordinate-system, i.e.
the local x-axis is along the beam on the line between the nodes, the y- and the z-axis are right
handed perpendicular to it. The z-axis defines the main bending direction and is in general
oriented downwards in the gravity direction.
For the description of the forces and moments and the support conditions, three points along
the beam have to be distinguished within a section:
• Centre of gravity (G) - This is the physical meaning of centre and is used for the calculation
of the self weight and masses.
• Elastic reference (S) - This is the centre weighted with the elasticity modulus of the gross
section. For homogeneous material it is identical with the centre of gravity, thus often just
named ”centre”. This value may include the basic reinforcements.
• Elastic centre (N) - This is the mechanical point where the normal force is acting. Stresses
induced by a constant strain within the total section will not generate any moments refer-
ence to that point. It differs from the elastic centre by the inclusion of holes and the current
reinforcements. Quite often it is also called the centre of the net or ideal section.
• Effective bending centre (B) - If some parts of the section do not take fully part due to
shear deformations, this may be accounted for bei noneffective areas for bending. This
has the consequence, that the bending strains and stresses are relative to the effective
centre. As a moment can be freely moved around, no additional forces are introduced by
this operation.
• Shear centre (M) - This point may be defined in several ways: A shear force acting at this
point does not create torsion. A section with torsion will rotate about this point. Thus if an
explicit rotation axis should be described, this point has to be defined explicitly.
• Beam axis - The beam axis is the (straight) connection between the beam nodes. The
section of the beam have to be defined at a specific location in the system, to allow the
correct treatment of all eccentricities. The sectional point on the beam axis is preferred to
be the origin of the sectional coordinate system (beam with reference axis). Alternatively
the elastic center may be used instead.
Hint
– Hinges are always defined at the elastic centre or shear centre.
– Enforced rotational axis must coincide with the beam axis.
Differences between these points will create changes in the forces and moments. On the other
side it is possible to describe a complete geometry with any eccentricities and unsymmetrical
haunched beams and construction stages with ease.
S
M
x, y, z Local beam coordinate system, freely selectable, is defined relative to the global
coordinate system by the definition of the beam.
y’, z’ sectional coordinate system for minimum moment of inertia (=coordinate system
shifted to the elastic centre)
For rotations the sign is always defined by the rotation about the local x-axis. This is clockwise
if you look in the direction of the beam and it is counter clockwise if looking on the positive face.
Sections will be saved in the database with the periphery in that same orientation. The sign of
the radius of a circular arc is defined positive if the area is increased compared to the secant,
and it is negative if the area is decreased as in the case of a fillet.
m torsional load
The rest of the parameters are static properties of the cross section (geometrical area mo-
ments). Since it is impractical to incorporate all of the static properties into the calculation,
certain standardisations are normally adopted:
• The axial force refers to the elastic centre of the beam, i.e.
Fy = Fz = 0 (2.2)
Fyz = 0 (2.3)
F = 0 (2.4)
• The torsional moment and the shear forces refer to the centre of shear, i.e.
Conversely, the conditions in 2.2 through 2.5 can be used in determining the elastic centre, the
orientation of the principal axes, the free moduli of warping and the centre of shear.
The determination of the area moments is relatively simple and shall not be described in any
further detail.
The normal stresses of a bar cross section can be described by means of Swain’s expression
and the normalized warping:
My z + Mz yz
N Mz y − My yz Mb
σ = + ·z− ·y+ · (2.6)
A 2
y z − yz 2
y z − yz CM
Effective width are also used for metallic sections of class 4, where parts of the section buckle
and do not allow any longitudinal stress.
AQUA is able to define the non-effective areas by means of rectangular regions with selectable
effectivity for different forces or moments. AQUA then stores the cross section values for the
total cross section as well as for the effective cross section.
Static analysis usually refers to the stiffness of the effective parts, whereas the normal force
(e.g. prestressing) refers to the center of the total cross section.
For the design the stresses from the normal force on the total section are then combined with
stresses from the bending moments on the effective section.
The non effective parts remain effective for shear stress calculations due to many consistency
reasons. But as the longitudinal stress is zero, the shear stress is always constant within an
element.
All the problems of the elasticity theory can be analysed by use of the force method or the
displacement method. While the force method is frequently used in calculations by hand, the
displacement method is better suited for processing with the computer. Both procedures are
implemented in AQUA for solid cross sections. The latter are by default computed by means
of the displacement method and in particular the Finite Element Method (FEM). Alternatively,
the Boundary Element Method (BEM) could be also activated. Certain simplifications of the
following equations can be made in case of thin-walled sections which facilitate a quick solution
of all tasks. These sections are therefore always analysed by the matrix displacement method.
A general formulation for the cross section warping according to the displacement method
conforms to the equilibrium condition
2
∂2
∂ ∂σ
G· + = − (2.7)
∂y 2 ∂z 2 ∂
The right side of 2.7 can be computed for example by 2.6. Assuming constant normal force
and constant cross section properties, one gets:
∂σ Vz z − Vy yz Vy y − Vz yz Mt2
− = · z + · y + · (2.10)
∂ y z − 2
yz y z − 2
yz CM
These equations will be approximated by AQUA either with the Boundary Element Method
(BEM) or the Finite Element Method (FEM).
For the Saint Venant’s torsion problem (∂Θ / ∂ = 1) the right side of 2.7 is identical to zero
As long as AQUA does not solve the differential equation 2.7, only an estimate of the torsional
moment of inertia is possible. The last equation shows that the polar moment of inertia can be
substituted for t in case of warp-free cross sections
t ≤ p = y + z . (2.13)
For all cross sections Eq. 2.13 provides an upper limit, which e.g. is about 10% above the
exact value for a square.
A4
T = (2.14)
4 · π 2 · y + z
This value is exact for circular and elliptical cross sections. For compact solid cross sections
this value provides a good approximation.
For hollow cross sections with more than 30% inner perimeters, an equivalent hollow cross
section based on the external and internal perimeters is used for a more refined estimate.
For composite sections Eq. 2.14 is used for each partial cross section and the components are
added.
can be solved with a variety of methods, which is controlled by the CTRL option STYP:
CTRL STYP FEM Finite Element Method - for all shear related cross-sectional prop-
erties (default)
CTRL STYP FEMX Finite Element Method as described above, but with a secondary
database containing the section mesh in a subfolder
CTRL STYP BEM 1 Boundary Element Method only for t and location of shear centre
ysc ,zsc
CTRL STYP BEM 2 Boundary Element Method for torsion and Force Method for shear
using shear cuts
CTRL STYP BEM 3 Boundary Element Method for torsion, shear deformation areas are
also determined
CTRL STYP FORC Force Method
Hint
The default of version 2014 can be set with CTRL STYP BEM.
In post-cracking (state II) analysis, AQB always employs the force method with proportional
axial force. In case of composite cross sections options BEM 2 and BEM 3 should be used
with caution. The input of explicit shear sections is required as a rule.
Hint
The Finite Element Method is the default method of computation in version 2016 for all
polygonal thick cross-sections.
If necessary, one could switch to any of the remaining analysis methods by means of the CTRL
option prior to the definition of a cross-section or in the cross-section editor.
It is also worth noting that cross-sections defined in version 2014 when reopened in version
2016 will automatically be converted by default to FEM cross-sections with the following con-
version scheme:
General Advantages
The Finite Element Method has several numerical and performance-related advantages over
the methods previously available in AQUA:
• provides a more stable and accurate approach for the elastic computation of torsional,
shear force and warping dependent properties and stresses
• provides a nonlinear incremental approach for the computation of plastic and other
nonlinear properties. FEM yields plausible values for the plastic shear forces and the
plastic torsional moment
• provides a more accurate estimation of the elastic limit forces Ne , My,e , Mz,e , Mb,e for
general cross-sections
• allows the direct computation of arbitrary geometry and complexity of hot sections with
HYDRA (thermal time-dependent analysis)
• provides symmetrical meshing of singly or doubly symmetrical cross-sections, thus en-
suring that symmetry of results is maintained and reduces numerical fluctuations of
derivatives (centroid, shear centre, points of symmetrical warping etc.)
• enables a detailed visualization of stresses and warpings with WINGRAF / Result
Viewer over the entirety of the domain
As with the remaining methods, one can also define none-effective parts, various isotropic /
anisotropic material models, shear cuts, reinforcements and construction stages which shall
be taken into account when constructing and meshing, interpolating, computing and post-
processing the model.
Secondary CDB
Upon activating the extended control option CTRL STYP FEMX the generated cross-sectional
mesh shall be saved into a databank which is located in a subdirectory pertaining to the project,
entitled projectname_sections. For example, Section 1 would then be saved in the following
database: projectname_sections\section_0001.cdb. The created database is an important
feature which is necessary for HYDRA if the section is to be later on used in a thermal analysis
or in a column design. It is not recommended to activate this option unless one is explicitly
interested in the physical output details of the analysis as its creation may demand considerable
additional HDD space and writing time. Out of optimization reasons, the default setting is CTRL
STYP FEM ensuring that the mesh for each cross-section shall be created and kept only for
the period of the analysis.
Databases with meshes could be also imported into AQUA using the SECT FEM definition for
a freely defined cross-section. The user is also allowed to import sectional systems meshed
by SOFiMSHC 2D as a whole or group-wise via the same principle.
The Animator could be used for a general graphical overview too. For most general pur-
poses, the Result Viewer would always extract the boundary warpings/stresses along those
quad edges belonging to the outer contours of the domain and save them in the main project
database, ready for visualization under the tree.
Hint
Explicitly defined stress points and single reinforcement bars shall not necessarily be
meshed in such a way that AQUA will always create a FE node at those locations. There-
fore results at those points are usually obtained via interpolation of results belonging to
the surrounding quads.
Mesh Quality
The quality of the FEM analysis depends entirely on how fine the underlying model is to be
discretized, i.e. the number of quad elements the mesh consists of. By default, AQUA aims at
providing a sufficiently reasonable mesh in order to obtain mechanically reasonably accurate
results for all cross-sectional elastic and plastic values, especially for steel shapes. One has
to however note that the last requirement is very difficult to be fulfilled for any possible input. It
is also a matter of compromise between accuracy, computational time and demand. More ge-
ometrically complex and involving, especially multi-domain / multi-material composite sections
might need further mesh densification in the zones of stress concentration, interfaces and fil-
lets. Same applies to sections which are to be subsequently undergoing thermal analysis with
HYDRA, where acceptable results can only be obtained with finer meshes. Before meshing a
model, AQUA determines for each subregion of the domain a corresponding mesh size which
is to be applied, also termed ”h-value”.
Hint
The cross-sectional FE-mesh provided by AQUA using the default mesh parameters is
in most cases sufficient to obtain mechanically reasonably accurate results of all cross-
sectional elastic and plastic values.
For the cases where a denser or a coarser mesh is required on an overall level, a control option
is provided by means of which the user could :
• factorize the ”h-values” of all geometrical regions by the same predefined factor CTRL
HDIV factor [ -] . For example a factor of 0.5[ -] would mean that each subdomain is to be
scaled down to 50% of its originally defined by AQUA mesh density. Same effect could
be achieved by the alternative input 50 [ %] . A factor greater than 1.0[ -] or 100[ %] can
be also applied.
• a direct input of a global mesh size by CTRL HDIV mesh size [ mm or m] . In this case
AQUA determines the ratio of the input value to the internally computed ”h-value” and
applies the quotient as a factor on that region.
Since the mesh size cannot be uncontrollably rescaled, AQUA has fixed allowable internal
limits in the range [ 0.1 ... 5.0] or [ 10% ... 500%] . Any value inputted beyond this scope forces
the program to rescale it to the nearest extreme and yield a warning.
It is more than recommended to use the first option since an estimate of the ”h-values” is quite
cumbersome and is a function of several complex invariable criteria.
Hint
Reducing the mesh could result in overwhelmingly long analysis time, large databases
and excessive data output.
In order to avoid various meshing problems, one should bear in mind that polygons must not
intersect each other at one or more singular points. It is thoroughly possible for polygons to
intersect each other in various fashions, provided that the intersections are continuous and
generated intersection points are not singular points.
Integration Schemes
When integrating any FE-related properties, AQUA/AQB would always employ a Lobatto-
integration scheme. For some very special cases, further integration schemes are available
using CTRL SINT. The options are :
1 Integration scheme using only quad mid points (useful for fine meshes)
2 Integration scheme using only quad nodes
3 Lobatto-integration scheme with nodes and quad mid-points (default)
Hint
It is always recommended to use the Lobatto-Integration scheme (default setting) in order
to achieve optimum result accuracy at the cost of a negligibly higher computational time.
Nevertheless, the user should be aware that this option also produces the most output as
(Hint continued...)
both FEM quad nodes and quad mid-points are being considered and additionally saved
into the database.
If SECT is given before PROF, the plastic shear values are also taken from section
tables (default CTRL PLAS 0)
SECT 3 MNO 1
PROF 3 TYPE IPE 300 DTYP S
However, with an explicit input of CTRL PLAS 1 the method of the determination of
the plastic shear values can be modified to FE-method.
• Freely defined sections with SECT obtain their plastic values by default according to
the method they are computed with:
Section 7 will be handled by FEM with the enforced condition of using the standard
method PLAS 1 (default if nothing else specified) in combination with the condition of
using elastic shear areas for Vy,p and Vz,p . Section 8 shall be computed by BEM,
therefore none of the PLAS controls related to FEM is relevant here.
If the CTRL PLAS 0 option is explicitly activated for the sections above, section 7 will
obtain only the plastic values of the PROF definition in them, disregarding any further
elements in the section definition which might lead to dangerously small and thus
erroneous plastic values. Users are thus advised to use PLAS 0 only in cases where
a single profile is to be analyzed. Finally, section 8 will remain unaffected.
• Any of the obtained in AQUA plastic characteristic values can be overwritten or
rescaled if needed as per cf. SV.
The example input below illustrates the activation of the advanced algorithm with the
forced condition of overwriting the computed plastic shear forces with the ones obtained
from the elastic shear areas. Example input file: [Link]/english/fem_plastic.dat
The number of elements determines the accuracy of the solution. In the case of a square,
for instance, the unit lateral warping on all the axes of symmetry is zero. A non-vanishing
solution can therefore be obtained only by defining at least four elements per side. AQUA
uses each polygon edge as one element, which can be further subdivided depending on its
size. Duplicate edges are automatically removed. As the results along the edges will vary only
linearly, it is strongly recommended to use the input value CTRL HMIN or the item SMAX in
POLY in order to have a coarse subdivision visible.
Internally a finer subdivision is however needed. Since a finer subdivision increases the com-
putational time with a power of two or three, the subdivision should not be made too fine. The
user can control the mesh size as absolute or relative by CTRL HDIV / SDIV.
The shear force components are then calculated by the classic formula
V S
τ = · (2.15)
b
The separated part of the cross section is defined on the positive side of a shear cut, which is
to the left of the cut’s direction for the positive face of the section. It is extremely important to
input the sections correctly, and especially to maintain their sequence.
For multiple connected cross sections, the difference of the normal stress is in equilibrium
only with the sum of all shear stress. Thus the user must either know the location of zero
shear stress or specify the distribution factors of the shear to the multiple segments of the cut.
Problematic are those where the modulus of the sectional part is very small.
Since in cracked sections the common usage is only based on the force method, these dis-
tribution values FVY and FVZ are needed in any case for reinforced concrete sections. The
displacement method allows however to establish reasonable estimates for many cases. Then
the entering and leaving shear stresses have different signs but nearly the same value. Thus
very large distribution factors are obtained, even if the sum of these still should remain 1.0
Torsional stress analysis is not elementary even for the force method (stress function with soap
film analogy). The resistance areas for the torsional shear stresses are therefore prescribed
by two values per section. The first value defines the shear at mid-area (Bredt’s equivalent
section). The second value defines the increase along the cut.
The recommended width of the equivalent section is in general twice the distance of the re-
inforcement to the outer edge. If the width is small, the kernel area is larger an the required
reinforcement areas are smaller. However the compressive stress in the concrete is larger and
if the allowed limits are exceede, the user should specify a larger width explicitly.
For the design of an equivalent hollow section we have additional integral values:
• A mean torsional shear stress which, after being multiplied by the width of the partial sec-
tion, must be covered by reinforcement. This corresponds to a section modulus for the
shear flow.
• A mean effective transverse shear or an effective width to be multiplied with the shear stress
at the mid of the equivalent hollow section to obtain the design shear flow from shear force.
To obtain these values in AQUA with CTRL STYP FEM/FEMX or BEM 2/3 the elastic shear
stresses are integrated over the total width. For a true hollow sections, this is the design value
directly, for an equivalent hollow section, the moment of the shear components is converted to
a constant shear flux across the width of the equivalent section (Figure 2.2).
beff
btot
Figure 2.2: Elastic and constant design shear in equivalent hollow cross section
Mt
Z
Ms = τs · r · dr = Mt · 1/ Wtm · btot − beƒ ƒ
≤ · btot − beƒ ƒ (2.16)
2Ak
These values will be slightly larger in general as the equivalent hollow section will have smooth
corners in the elastic analysis. If a distinct value of Ak is defined or Ak is given with a torsional
reinforcement in a section, this value will be used to establish an upper limit on the torsional
shear flow.
The principal compressive stress will be calculated with the equivalent width from the design
shear flow.
identification, which consists of three characters. The cut can be defined parallel to an axis or
as a free form polygon line. Every segment has its own material number and it will only cut
through cross section elements with the same material number. Gaps between the segments
will be closed by means of virtual connections. The width of the substitute torsional cross
section is available as a special option for the description of equivalent hollow cross sections
of reinforced and prestressed concrete. Two partial cuts are generated for each section in this
case.
If the user does not supply any input, one or two axis-parallel cuts will be created through the
elastic centre. This is generally not sufficient even for a simple T-Beam, nor for composite
sections, where the reference material number of the section is not necessarily represented at
that location. The user will see a warning for general sections therefore.
CTRL SCUT allows the user to control how many of these standard cuts will be generated
(0/1/2).
The cut can dissect the cross section at several locations creating partial cuts. Each partial cut
has a direction s and three defined points of interest: beginning (A), middle (M) and end (E):
T M
s z
E E A
n
N
A
M
The internal forces perpendicular to the cut M and N act in such way that positive axial forces
cause tensile stresses across the cut, and positive moments cause tensile stresses at the
End-Point.
Several polygons of the same material are united before the sectional check. In case of com-
posite cross sections it may be helpful to pay attention to the cut direction.
In the definition of cuts across several materials the user must take care that each segment of
the cut has the correct material number, because a cut will dissect only parts having the same
material number. It makes a difference for the horizontal shear in a composite flange if a dowel
is before or behind the cut.
Cuts through cross sections with ”open air” between their parts can not be analysed as the
section does not hit any elements. A similar problem occurs if a cut has the wrong material
number. This may happen especially with the standard cuts through the elastic centre.
Some additional advice applies to oblique cuts. Since the shear force at an oblique cut does
not vary significantly compared to the straight cut, however the width of the cut does. Since the
selection of an inappropriate cut direction can result in the analysis of too small shear stresses.
AQUA uses the deformation method in the frame of a simplified FE model for thin-walled sec-
tions in all cases. The CTRL option STYP has in here no effects at all, but you may specify
explicit values via SV.
Materials 999
Materials per cross section 31
Cross sections 9999
Longitudinal reinforcement layers 10
Polygon vertices per polygon 255
Shear sections per cross section 255
Shear reinforcement layers 15
Construction stage per cross section 31
3 Input Description
AQUA allows the user to define general cross sections with arbitrary geometry and materials.
For simple sections and materials you do not need a special license, but for all sections starting
with record SECT you need a license for AQUA.
Before defining a section you have to specify the materials. Materials are addressed by an
arbitrary number. Please note, that by keeping track of construction phases in AQBS, it is
assumed that materials with higher material numbers were added at a later time.
A standard section is defined by just one input record. All sectional values will be calculated
including torsional and shear properties. The maximum components for all stresses are known,
but a detailed analysis at different locations within the section will not take place.
With AQUA cross sections can be redefined at any time during the processing of the project
without affecting other defined sections. However if any material definition is made, all existing
cross sections are deleted. The distributions of reinforcements and stresses are deleted too,
unless otherwise specified with CTRL REST.
Freely defined cross sections always start with the record SECT, which specifies the cross
section number. All subsequent input records describe this one cross section, which may
consist of several partial cross sections (external perimeter, inner perimeter, reinforcement
layout etc.). The input for a cross section is concluded either by the next SECT record or by
two END records.
3.2 Units
SOFiSTiK programs offer the possibility to carry out all input and output of data in engineering
units. A number of unit sets are provided for this purpose, which are preset according to
the design code used in the given project. This default can additionally be changed for each
program run separately using the keyword PAGE. More information about unit sets can be
found in the general SOFiSTiK manual, section ’Units’.
[mm] Explicit unit. Input defaults to the specified unit. Alternatively, an explicit as-
signment of a related unit is possible (eg. 2.5[m] ).
[mm] 1011 Implicit unit. Implicit units are categorised semantically and denoted by a cor-
responding identity number (shown in green). Valid categories referring to the
unit ”length” are, for example, geodetic elevation, section length and thickness.
The default unit for each category is defined by the currently active (design code
specific) unit set. This input default can be overridden as described above. The
specified unit in square brackets corresponds to the default for unit set 5 (Eu-
rocodes, NORM UNIT 5).
For sections the units for all dimesnions are expected in [ mm] in general, the unit sets 0, 3 and
4 expect [ m] , the Unit-Set 1 [ cm] . For reinforcement areas values are expected in [ cm2 ] ,
the unit sets 6 and 7 expect [ mm2 ] .
The default unit set (UNIT) of a corresponding design code is described in the record NORM
at the respective design code in the tables.
Record Items
CTRL OPT VAL VAL2
NORM DC NDC COUN CAT ALT WIND SNOW
SEIS
MATE NO E MUE G K GAM GAMA
ALFA E90 M90 OAL OAF SPM FY
FT TYPE TITL
MAT NO E MUE G K GAM GAMA
Table continued on next page.
Record Items
ALFA EY MXY OAL OAF SPM TITL
FLAY NO T MNO BETA TITL
MLAY NO T0 NR0 T1 NR1 T2 NR2
T3 NR3 T4 NR4 T5 NR5 T6
NR6 T7 NR7 T8 NR8 T9 NR9
TITL
NMAT NO TYPE P1 P2 P3 P4 P5
P6 P7 P8 P9 P10 P11 P12
BMAT NO C CT CRAC YIEL MUE COH
DIL GAMB TYPE MREF H
HMAT NO TYPE TEMP KXX KYY KZZ KXY
KXZ KYZ S NSP A B C
QMAX TK TITL
CONC NO TYPE FCN FC FCT FCTK EC
QC GAM ALFA SCM TYPR FCR ECR
FBD FFAT FCTD FEQR FEQT GMOD KMOD
GC GF MUEC TITL
STEE NO TYPE CLAS FY FT FP ES
QS GAM ALFA scm EPSY EPST REL1
REL2 R K1 FDYN FYC FTC TMAX
GMOD KMOD QS TITL
TIMB NO TYPE CLAS EP G E90 QH
QH90 GAM ALFA SCM FM FT0 FT90
FC0 FC90 FV FVR FVB FM90 OAL
OAF KMOD KMO1 KMO2 KMO3 KMO4 KDEF
TMAX RHO TITL
MASO NO STYP SCLA MCLA E G MUE
GAM ALFA SCM E90 M90 OAL OAF
FCN FC FT FV FHS FTB TITL
SSLA EPS SIG TYPE TEMP EPST TS MUET
MNOC FCTF
MEXT NO EXP TYPE VAL VAL1 VAL2 ...
VAL9
BORE NO X Y Z NX NY NZ
ALF TITL
BLAY S MN0 ES MUE DES VARI PMAX
Table continued on next page.
Record Items
PMAL C PHI
BBAX S1 S2 K0 K1 K2 K3 M0
C0 TANR TAND KSIG D0 D2 CA0
CA2
BBLA S1 S2 K0 K1 K2 K3 P0
P1 P2 P3 PMA1 PMA2 CL0 CL1
CL2 CL3 SM0 SM2
SMAT NO LTYP MTYP P+ P- VY+ VY-
VZ+ VZ- MY+ MY- MZ+ MZ- ALPH
P1 P2 TITL
SFLA NO U F S SH FP TYPE
LEV TITL
LCA NO TYPE PHAS VAL
SVAL NO MNO A AY AZ IT IY
IZ IYZ CM YSC ZSC YMIN YMAX
ZMIN ZMAX WT WVY WVZ NPL VYPL
VZPL MTPL MYPL MZPL BCYZ TITL
SREC NO H B HO BO SO SU
SS MNO MRF MRFL RTYP ASO ASU
DASO DASU DASS A AMIN AMAX ASL
INCL REF YM ZM IT AY AZ
BCYZ SPT INTE CINT MUE BEFF BDL
TITL
EDGB NO H B HO BO SO SU
SS MNO MRF MRFL REF YM ZM
FSYM FLIP ASO ASU DASO DASU DASS
A AMIN AMAX PASU ASL INCL IT
AY AZ SPT CUTZ CUTF INTE CINT
MUE BEFF BDL TITL
SHRW NO B L SS SZ+ SZ- MNO
MRF MRFL A D RTYP LZ+ AZ+
DZ+ RZ+ LZ- AZ- DZ- RZ- TITL
SCIT NO D T SA SI MNO
MRF MRFL RTYP ASA ASI DAS A
ASL IT AY AZ TITL
CABL NO D TYPE INL MNO F K
Table continued on next page.
Record Items
W KE REF GAMR TITL
SECT NO MNO MRF ALPH YM ZM FSYM
BTYP BCY BCZ KTZ TITL
- CS NO TITL ATIL
- SV IT AK YSC ZSC CM CMS AY
AZ AYZ LEVY LEVZ MNO DEFF FACE
FACG AG
- POLY TYPE MNO YM ZM DY DZ SMAX
EXP REFP REFD REFS
- - VERT NO Y Z R PHI TYPE EXP
REFP REFD REFS
- CIRC NO Y Z R MNO EXP
REFP REFD REFS REFR
- NEFF TYPE YMIN ZMIN YMAX ZMAX MNO WIDT
NO SMIN SMAX
REFI RFDI RFSI REFA RFDA RFSA
- CUT NO YB ZB YE ZE NS MS
WTM WTD MNO MRF LAY ASUP
OUT TYPE VYFK VZFK INCL BMAX
BRED BCT MUE SXE TANA REFA
RFDA RFSA REFE RFDE REFS
- PANE NO YB ZB YE ZE T MNO
REFA RFDA RFSA REFE RFDE RFSE
R PHI OUT FIXB FIXE TYPE
AS ASMA LAY MRF TORS DAS
A
- PLAT NO YB ZB YE ZE T MNO
REFA RFDA RFSA REFE RFDE RFSE
R PHI OUT FIXB FIXE TYPE
- WELD NO YB ZB YE ZE T MNO
REFA RFDA RFSA REFE RFDE RFSE
- PROF NO TYPE Z1 Z2 Z3 MNO ALPH
YM ZM REFP REFD REFS REFR
DTYP SYM REF MREF VD VB
VS VT VR1 VR2 VB2 VT2
CW BCYZ WU1 WU2 WU3
Table continued on next page.
Record Items
- SPT NO Y Z WTY WTZ WVY WVZ
SIGY TEFF CDYN SIGC TAUC MNO
FIX REFP REFD REFS
- SFLA NO U F S SH FP TYPE
LEV TITL
- WPAR CS KR ICE TRAF YMIN YMAX ZMIN
ZMAX YREF ZREF CW ALF
REFI RFDI RFSI REFA RFDA RFSA REFP
- WIND ALPH CWY CWZ CWT REF CHYD CLAT
S AG VR0 V0 ... VR19 V19
- RF NO Y Z AS ASMA LAY MRF
TORS D AR SIG TEMP REFP
REFD REFS
- LRF NO YB ZB YE ZE AS ASMA
LAY MRF TORS D A DIST AR
REFA RFDA RFSA REFE RFDE RFSE
R PHI
- CRF NO Y Z R PHI AS ASMA
LAY MRF TORS D A DIST AR
REFP REFD REFS REFR
- CURF H EXP AS ASMA LAY MRF
TORS D A DIST AR CENT
TVAR NAME VAL SCOP CMNT
INTE NO NS0 NS1 S NREF ICS ...
ICS9
IMPO MAT SECT FROM
EXPO MAT SECT TO PASS
ECHO OPT VAL VAL2
The records HEAD, END and PAGE are described in the general manual SOFiSTiK: ’Basics’.
CONC
SSLA
Definitions of materials (can be
skipped, if not all cross sec-
STEE
tions are redefined)
SSLA
TIMB
Type of cross section
SV SV
POLY PANE
VERT PLAT
CURF PROF
PROF WELD
CIRC
SPT SPT
CUT
RF RF
LRF LRF
CRF
END
CTRL
The CTRL options may be defined at any location within the input data. However if they are
intended to be different for individual sections it is mandatory to define them before the sections
intended to use them. Within a restart CTRL REST 3 all defined options will be effective for all
reanalyzed sections.
CTRL REST controls what AQUA should do with existing data in the database. As default
AQUA will erase everything if materials are defined, and only the basic reinforcements, limit
stresses and beam stiffness if only sections are defined. This is usually the best choice to
avoid unforeseeable results. In some cases though it is desirable to process these results
further. This can take place without problems only, if the assignment of the layers and the use
of the material numbers in the individual cross sections are not changed.
0 deletes all materials, cross sections and all beam design results of the following pro-
grams in the database (default)
1 keeps all materials and cross sections and all beam design results of the following
programs in the database
2 re-calculation of all cross sections and materials only when a material is input and
also retention the beam design results in the database. When no material is new
defined only the new defined cross sections are calculated and retention the beam
design results in the database.
CTRL REST 2 nn also forces the calculation of the cross-section nn.
3 re-calculation of all cross sections and materials even when no material is new de-
fined and retention the beam design results in the database
CTRL FACE defines the standard view on the section. While the physical orientation in space
is only specified by the orientation of the local beam coordinate system, the possible values for
graphical views on the section are defined that the literal POS and all positive numbers define
a view on the positive face (i.e. in the inverse direction of the beam or axis) while the literal
NEG and all negative values define a view on the negative face (ie. along the direction of the
beam or axis) :
CTRL RFCS controls whether the basic reinforcement should be considered in the calculation
of the cross section values:
0 do not consider
1 do not consider for Boundary Element Method ("CTRL STYP BEM"), otherwise con-
sider for sections with more than one material (default)
2 consider for all sections (also for circular sections SCIT and edge beam section
EDGB, but not for rectangular sections SREC)
3 consider also effect on dead load (the weight of the concrete should be reduced in
that case)
HMIN defines a maximum allowed length for linear or circular polygon edges. (Default: no limit)
HTOL defines the maximum allowed stitch (error) of an approximation of a circular arc by a
polygon. This is applicable to fillets and arcs. Default is: 0.0035. There are two alternative
ways of defining an allowed stitch:
• HTOL factor [ -] - define the stitch in the way that all circles will use the same angular
resolution, which is nicer for visualization. But for a small radius below 50 mm the angle
of the segments will be enlarged up to a factor of 6.
In both cases the angles of the segments will not be taken less than 1 degree or larger than 30
degrees.
HDIV defines a maximum element length for the boundary elements (BEM). In case of FE
cross-sections, HDIV shall be interpreted as mesh density scaling. There are two alternative
ways of applying a scale factor on the mesh density:
• HDIV factor [ -] - for example a factor of 0.5 [ -] would mean that each subdomain
is to be scaled down to 50% of its originally defined by AQUA mesh density. Same
effect could be achieved by the alternative input 50[ %] . A factor greater than 1.0 [ -] or
100[ %] can be also applied.
• HDIV mesh size [ mm or m] - a direct input of a global mesh size.
In this case AQUA determines the ratio of the input value to the internally computed
”region element size” and applies the quotient as a factor on that region.
Since the mesh size cannot be uncontrollably rescaled, AQUA has fixed allowable internal
limits in the range [ 0.1 ... 5.0] or [ 10% ... 500 %] . Any value inputted beyond this scope
forces the program to rescale it to the nearest extreme and yield a warning.
SDIV defines a relative desired element length per edge for the boundary elements. For finite
elements this should be kept to the default of 0 and if not so, the value will be ignored.
0 no explicit subdivision
1 maximum 1/2
2 maximum 1/4
3 maximum 1/8
4 maximum 1/16
5 maximum 1/32 (Default BEM)
A snap distance (always in m) for the detection of cross section parts connecting together can
be defined additionally at item VAL2. The value SDIV 4 0.001 defines 1 mm, as snap measure.
The detailed background of STYP and SCUT are explained in Section 2.8.
STYP FEM Finite Element Method - for all shear related cross-sectional properties
(default)
STYP FEMX Finite Element Method as described above, but with a secondary
database containing the section mesh in a subfolder
STYP BEM Boundary Element Method with default sub option: 1 for concrete, 3 oth-
erwise)
STYP BEM 1 Boundary Element Method only for It and location of shear centre ysc ,zsc
(Default for concrete)
STYP BEM 2 Boundary Element Method for torsion and Force Method for shear using
shear cuts
STYP BEM 3 Boundary Element Method for torsion, shear deformation areas are also
determined
STYP FORC Force Method
Hint
The default of version 2014 can be set with CTRL STYP BEM. However, for rectangular,
T-Beam and circular explicit concrete sections, the values Ay and Az must be set to 0.0
explicitly if required.
For FE solid sections a FE-method of scaling the elastic shear stresses flux is
used.
+1*16 Replace current method’s computed Vy,p and Vz,p by plastic forces based on
the scaled elastic shear areas
For more information the user is kindly advised to consult Sections 2.8 and 2.10
of the AQUA manual for a detailed description of the control options and their
peculiarities.
By specifying the additional literal BUCK at item VAL2 it is possible to activate the
considering of buckling (reduction of cross section area) of thin walled elements
for a negative normal force -N (compression) or in the compression zone of the
bending moment My or Mz. Considering of buckling is not the default setting
because bracings are normally not applied in the design phase.
By specifying the additional literal GAM at item VAL3 it is possible to define the
material safety factor for the design forces via item VAL4:
- Material safety for cross sectional design (γM0 of steel) (default)
M0 Material safety for cross sectional design (γM0 of steel)
SCM Material safety for stability check
The material safety of single materials can be modified via the material records,
e.g. for γM0 of a steel with ’STEE SCM0 ...’.
Input example for negative limit forces considering buckling with material safety for stability
check (γM1 ):
For a composite section according to EN 1994-1-1, ch. 6.7. the concrete is used with 0.85*fcd
in the compression zone of Npl and Mpl (ch. [Link], [Link]). There are similar requirements
in the codes SIA, AASHTO, DECRETO MINISTERIALE and so on. A composite section is
detected if there are only steel and concrete materials and if the reference material is a steel
material.
A concrete section is detected if the reference material is a concrete material. Then no plastic
limit forces are computed, except of in case of a ’hot’ design (temperatures from program
HYDRA), even if steel materials exists.
Thin walled sections are always calculated with a FEM-type method. For solid sections the
more advanced algorithms may create high stress values at singularities. To disable these
values for the design, it is possible to select STYP with option 1, defining that for all shear cuts
the simplified values according to the force method are used. If options 2 or 3 are used, the
integral mean values will be used for this case and all other stress points will get the calculated
values.
SCUT controls the generation of the two standard cuts parallel to the coordinate axis and the
standard stress points for solid cross-sections:
As the input of references is best done with relative offsets (default: CTRL REFD 1), the
export however with the actual absolute coordinates, the input CTRL REFD 0 allows to change
the default to the definition in absolute coordinates. In that case it is also allowed to make
a reference to an element defined later in the input stream, but all coordinates will be taken
exactly as they have been specified. Reference IDs should be unique within every element
type. As long as the sections are only defined via CADINP in specified sequence, same Ids
may be used however. With CTRL REFD 2+xx the algorithm may be changed to the old purely
sequential scheme for test purpose.
3.5 Materials
SOFiSTiK supports a large number of different material descriptions. All will be addressed by
a unique material number and should in general be usable everywhere. The default for the
material type and any national deviant parameters is established dependent on the selected
design code NORM. In some cases parameters may be also changed by the selection of a
country code. All values may be modified explicitly if required.
For any material CONC, STEE, TIMB and MASO a unique literal for the material subtype and
a class value are required. The class is always a string with integer numbers and/or optional
characters. In general these two values define all material properties (Some require also a
density or a thickness) It is then possible to change any value individually, although there are
a few implicit dependencies.
These records are mutually exclusive but may be enhanced by other records:
Input of material is possible in all parts of the program system. However, it is self-evident that
not all parameters are used for all types of analysis or system. Each material has a standard
name given by its classification, which might be extended by the user. If the user wants to
replace the standard completely, he has to start his own text with an exclamation mark (e.g.
’!my own Text’) or to quote it a second time (eg. " ’my own Text’ ").
Properties of materials must be distinguished according to whether they are properties which
are close to the realistic behaviour (e.g. for dynamic calculations) or to which have some
lower or upper limit to be multiplied with a safety factor for the calculation of an ultimate load-
bearing capacity. Whereas the safety factors were formerly assigned more-or-less at random,
sometimes to the load and sometimes to the material, more recent regulations (Eurocode)
provide a clearer separation between safety factors for the loads and factors for the material.
However, since the material safety factors still depend on the nature of the load or the type of
design, it will not be possible to define all safety factors with the material itself.
If some design codes (DIN 18800, DIN 1045-1) apply additional safety-factors to the mean
values, this may be defined with the stress-strain relation via SSLA. The safety factor defined
with the material will thus be used only for the full plastic forces in AQUA.
Note: The following pages are valid in all details only for AQUA, for other programs
(SOFiMSHB) deviations are possible due to older versions with missing or changed items.
Note: Hints for material properties of strange materials may be found on the internet at www.
[Link] (The A to Z of Materials).
NORM
Many defaults for materials, superposition and design are selected according to the selected
design code and an optional country code and all the other data provided with this record. It is
therefore requested to specify this data with the beginning of the project.
A redefinition of the design code after the definition of actions or load cases have been
defined or the editing of the INI-File to include ”missing” materials does not comply with
the provisions of SOFiSTiK for a proper use of the software.
It is possible to redefine the design code NORM temporarily for the design (eg. concrete / steel)
if the parameters of the actions remain the same. As this might have some special risque, the
user should use this option thoroughly.
Although there are still explicit code fragments in the software unavoidable, many of the defaults
are specified in so called INI-Files located in the SOFiSTiK directory. The name of the matching
INI-file is derived from the given data as DC_NDC.INI.
Some properties (e.g. Eurocode) are dependant on national variants (boxed values). Corre-
sponding INI-files to EN 1992-2004 and EN 1993-2005 and the country code may be used to
select those values, as far as we have got notice of them. The country code for example is valid
for deviations in Hong Kong to the British Standard 8110, 5400, 5950 or similar. Examples for
the country code:
00 General EN
HK / 852 Hong Kong
The ”boxed values” can be modified manually with the record TVAR.
Some codes require or allow the selection of a category or class. This can then be specified
with CAT. The possible items are given in the INI-File. For evaluations with historic design
codes not available with an INI-file, the definition of CAT USER allows to specify any name of
the design code.
In the case of a subsequent modification of these classes or categories the input which
depends on the design code has to be checked and adapted if necessary.
The extend to which the specified latitude, altitude, wind/snow or earthquake zone definitions
are accounted for is described in the program manuals of the modules using those values. The
user should never assume that all regulations of the design codes are automatically fulfilled
when selecting such a value. The possible items and defaults are given in the matching INI-
File. The resultant values which result from the altitude or the wind/snow/earthquake zones
have to be checked in the corresponding programs in the case of a subsequent modification.
E.g. for some design codes the combination coefficients of the snow depend on the altitude.
In the case of a modification of the altitude combination coefficients have to be adapted by the
user if necessary.
If the user wants to suppress such a value completely he may specify it with ”NONE”. The
load case of self weight is always accounted for following the definitions of the design codes,
an explicit definition of G has to be done only in those cases where the true local value should
be applied.
The item UNIT or d LANG will be processed only in AQUA or TEMPLATE. With a definition of
UNIT a set of units will be selected globally for all input and output data in all other modules.
The default is specified in the INI-file. Definitions with record PAGE will be active only within
the current module.
The following design codes are available as INI-Files and/or special program code has been
created to cope with special regulations. The marks A and B indicate if this code has been
implemented in AQB and BEMESS. For more detailed information, especially which provisions
of the codes have been implemented, please check the manuals and the HTM-files of the
design programs. In many cases it is possible to add some clauses within short time within the
program or with CADINP.
EN - Eurocodes
For the old design codes OEN 4200, OEN 4250, OEN 4253 no INI files exist. The program
AQB is so programmed that the appropriate design is done with input of the design code. As
BS - British Standard
SP - Russian Standard
IS - Indian Standard
AS - Australian Standard
SS - Swedish Standard
DS - Danish Standard
DS EN 1991-1-x/NA:200x
including materials according to:
DS EN 1995-1-1 Timber
DS EN 1996-1-1 Brickwork
DS EN 1999-1-1 Aluminium
CAT LE low safety & extended control
CAT NE normal safety & extended control
CAT HE high safety & extended control
CAT LN low safety & normal control
CAT NN normal safety & normal control (de-
fault)
CAT HN high safety & normal control
CAT LR low safety & reduced control
CAT NR normal safety & reduced control
CAT HR high safety & reduced control
safety according to DS EN 1990/NA:2007
control according to DS EN 1992-1-
1/NA:2007 and DS EN 1993-1-1/NA:2007
inspection level
DS EN199X-200X-BRIDGE National Annexes DS bridge construction 5 A,B
DS EN 1992-2/NA:2015 Design of concrete
structures
DS EN 1993-2/NA:2015 Design of steel
structures
DS EN 1994-2/NA:2005 Design of compos-
ite steel and concrete structures
including actions and combinations accord-
ing to:
DS EN 1990/A1/NA:2017
DS EN 1991-2/NA:2017
including materials according to:
DS EN 1995-1-1 Timber
DS EN 1996-1-1 Brickwork
DS EN 1999-1-1 Aluminium
for road bridges:
CAT BNE normal safety & extended control
CAT BHE high safety & extended control
Table continued on next page.
NS - Norsk Standard
GB - Chinese Standard
ZA TMH7 Code of Practice for the Design of Highway Bridges and Cul- 6 A
verts in South Africa (1989)
CAT 0 without prestress
CAT 1/2/3 prestress for class 1/2/3
MATE
NO Material number − 1
E Elastic modulus [MP] 1090 *
MUE Poisson’s ratio (between 0.0 and 0.49) − 0
automatic default for nonlinear materials NMAT
G Shear modulus [MP] 1090 *
K Bulk modulus [MP] 1090 *
kN/ m3 1091
GAM Specific weight 25
kN/ m3 1091
GAMA Specific weight under buoyancy *
ALFA Thermal expansion coefficient [1/ K] 1083 E-5
kN/ m3
E90 Anisotropic elastic modulus 1090 E
M90 Anisotropic Poisson’s ratio − MUE
OAL Meridian angle of anisotropy about the local x axis deg 0
OAF Descent angle of anisotropy about the local x axis deg 0
MATE defines Materials which cannot be defined by CONC, STEE, TIMB or BRIC otherwise.
The number of these materials must not be used for other materials.
With the definition of a literal at TYPE from the following list, default values will be selected:
The mechanical properties of those materials are not always known with a distinct value or
even suitable for a linear material description at all!
For geotechnical design the TYPE may be specified by the Unified Soil Classification System
(USCS) :
For fine grained soils the effect of the liquidity index will be defined by the uniaxial strength FY,
which may be preset by appending a literal to the type:
Further TYPE may be used to preset fluid material constants (depending strongly on the tem-
perature) for:
where MUE is the kinematic viscosity ν [ m2 /sec] . Sometimes it is more convenient to define
the elastic constants by other values than the Elasticity modulus, the Shear modulus and the
Poisson ratio. You may transform your values by the following formulas:
E Elastic modulus
Es subgrade modulus (horizontally constrained)
K Bulk modulus
G Shear modulus
μ Poisson’s ratio
E E
K = G = (3.1)
3(1 − 2μ) 2(1 + μ)
9·K ·G 3K − 2G
E = μ = (3.2)
(3K + G) 6K + 2G
(1 − μ)
Es = E · (3.3)
(1 + μ)(1 − 2μ)
3 · K · E 3 · K · (1 − 2μ)
G = G = (3.4)
9 · K−E 2 · (1 + μ)
If not specified, missing values will be calculated according to these formulas. It is, however,
possible to define non-consistent constants. If no values are given, E defaults to 30000 MPa
and MUE to 0.0.
Orthotropy may be defined via material and thickness of QUAD elements. (see record QUAD
in SOFiMSHA and manuals to ASE and TALPA).
The material law of a transversal orthotropy according to Lechnitzky has one direction that
has different properties, while the description in the plane perpendicular to this direction re-
mains isotropic. This covers most practical problems like timber and rock. Unfortunately the
designation of the general material constants for timber materials is not compatible with that
convention, so we strongly recommend to use TIMB for the description of timber materials.
The principal values E and μ are related to the isotropic plane, while E90 (MAT still old literal
EY) acts normal to that plane and μ90 and G90 (attention: input parameter G) describe the
transverse straining behaviour between normal direction and isotropic plane. With z being the
normal direction we have:
σ σy σz
ε = −μ· − μ90 · (3.5a)
E E E90
σy σ σz
εy = −μ· − μ90 · (3.5b)
E E E90
σz (σ + σy )
εz = − μ⋆90 · (3.5c)
E90 E
E
μ⋆90 = μ90 · (3.5d)
E90
It should be noted, that poisson’s ratios μ90 and μ90∗ are no longer limited by an upper bound
of 0.5 (this would hold for the isotropic case). According to the formula given above, their
relation is determined by the ratio of the elasticity moduli; this preserves symmetry of the
stress strain matrix, which is required to fulfill elemental equilibrium.
The order of the indices of stress and strain components for subsequent equations is defined
as:
With axial symmetry x denotes the axis of rotation while y represents the radial and z the
tangential direction.
Furthermore holds:
E1
E1 = E , E2 = E90 , μ1 = μ , μ2 = μ90 , G1 = , G2 = G (3.6)
2(1 + μ1 )
μ1 +n·μ2
1−n·μ2 μ
E ·
1 (1+μ1 )·m
2 2
E1 · (1+μ )·m E1 · m2 0 0 0
1
E1 · μ1 +n·μ2
2 1−n·μ2
μ
(1+μ1 )·m
2
E1 · (1+μ )·m E1 · m2 0 0 0
1
μ μ 1−μ1
E1 · m2 E1 · m2 E2 · 0 0 0
D=
m
(3.7)
0 0 0 G1 0 0
0 0 0 0 G2 0
0 0 0 0 0 G2
E1
n = , m = 1 − μ1 − 2 · n · μ2
2
(3.8)
E2
E1
n = , m = 1 − μ1 − 2 · n · μ2
2
(3.10)
E2
Axial symmetry:
A general case of anisotropy does not need to be considered since axial symmetry would be
impossible to achieve under such circumstances. A case of interest in practice is that of a
stratified material in which the rotational axis x is normal to the plane of isotropy. For such a
case the material stiffness matrix reads:
1−μ21
n μ2 (1 + μ1 ) 0 μ2 (1 + μ1 )
μ2 (1 + μ1 )
1 − nμ2 0 μ1 + nμ2
D=A· 2 2
(3.12)
G90
0 0 0
A
μ2 (1 + μ1 ) μ1 + nμ2
2 0 1 − nμ2
2
E2 · n E1
A= , n= , m = 1 − μ1 − 2 · n · μ2
2
(3.13)
(1 + μ1 )·m E2
Skew orthotropy:
Consideration of ’skew’ orthotropy is also possible. In geological terms, the three-dimensional
orientation of the isotropic plane is defined by means of the meridian and descent angle. They
describe the deviation of the steepest descent to the north direction and the inclination of the
layers. Mathematically, the angles are equivalent to the first and third of the Eulerian angles.
The transformation is defined by two rotations, the north axis (N) corresponding to the element
y-direction and the G-axis corresponding to the element z-direction. Axes K, N and G form a
right handed Cartesian coordinate system.
Apart from 3D continuum elements these transformation rules apply to shells and plates, as
well.
For planar systems (TALPA) the value OAL defines the slope of the stratification, i.e. the angle
between the element x-direction and the stratification direction. Input for OAF is not evaluated
for the plane case.
For axial symmetry input of OAF and OAL is not evaluated (see above: axial symmetry).
MAT
NO Material number − 1
kN/ m2
E Elastic modulus 1088 *
MUE Poisson’s ratio (between 0.0 and 0.49) − 0
automatic default for nonlinear materials NMAT
kN/ m2
G Shear modulus 1088 *
kN/ m2
K Bulk modulus 1088 *
kN/ m3
GAM Specific weight 1091 25
kN/ m3
GAMA Specific weight under buoyancy 1091 *
ALFA Thermal expansion coefficient [1/ K] 1083 E-5
kN/ m2 1088
EY Anisotropic elastic modulus Ey E
MXY Anisotropic Poisson’s ratio m-xy − MUE
OAL Meridian angle of anisotropy about the local x axis deg 0
OAF Descent angle of anisotropy about the local x axis deg 0
SPM Material safety factor − 1.0
TITL Material name Lt32 -
Materials which can be used for SVAL or QUAD and BRIC elements may be defined with the
record MAT and MATE. The number of the material must not be used for other materials.
The differences between the two records are mainly the used dimensions. MATE is analogue to
CONC,STEE etc. (MPa) and has additional strength values, while MAT uses (kN/m2 ) analogue
to NMAT. MAT has older item names for the orthotropic parameters.
FLAY
With FLAY you may define a composite layered material for QUAD elements with up to 99
layers. Each layer may be defined with a absolute thickness (positive) or a relative (negative)
one. The total thickness of the element is calibrated to the sum of the thicknesses of the
material definition. If some layers have negative thickness, only these layers are adapted.
Otherwise a uniform scaling will take place. The fibres direction of layer may be defined via
BETA. The direction of layers material (eg meridian angle OAL) is ignored. If NO is not defined,
last valid input for NO is used. So it is sufficient to input NO only at first layer. Same for
description TITL.
For example, if you have a sandwich element with two outer laminates with a given thickness:
then this data is applied to match two QUAD elements with a total thickness of 0.10 or 0.15 as
follows:
MLAY
With MLAY you may define for QUAD elements a composite layered material with up to 10
layers. Each layer may be defined with a positive absolute thickness or a negative relative
one. The total thickness of the element will be calibrated to the sum of the thicknesses of the
material definition. If some layers have negative thickness only these layers will be adapted.
Otherwise a uniform scaling will take place.
If you have a sandwich element with two outer laminates with a given thickness for example:
then this data will be applied to match two QUAD elements with a total thickness of 0.10 or
0.15 as follows:
NMAT
NO Material number − 1
TYPE Kind of material law LT !
LINE Linear material
VMIS von Mise law, optional viscoplastic exten-
sion
DRUC Drucker-Prager law, optional viscoplastic
extension
MOHR Mohr Coulomb law
GRAN Hardening Plasticity Soil
SWEL Swelling
FAUL Faults in rock material
ROCK Rock material
MISE Mise / Drucker Prager law
GUDE Gudehus law
LADE Lade law
MEMB Textile membrane
USP1 to USP8 and USD1 to USD8 reserved for
user defined material models
UNDR Undrained Stress Analysis
The types of the implemented material laws and the meaning of their parameters can be found
in the following pages.
In a linear analysis the yield function for the non-linear material is merely evaluated and output.
This enables an estimation of the non-linear regions for a subsequent non-linear analysis.
• The first principal invariant of the Cauchy stress tensor σ is defined as:
1 = σ + σy + σz (3.14)
1
• The deviatoric stress tensor can then be established as s = σ − 3 1:
1
s = σ − (3.15a)
3
1
sy = σy − (3.15b)
3
1
sz = σz − (3.15c)
3
with its second and third principal invariants:
1
J2 = (s2 + s2 + s2 2
) + τy 2
+ τyz 2
+ τz (3.16a)
2 y z
2 2 2
J3 = s sy sz + 2τy τyz τz − s τyz − sy τz − sz τy (3.16b)
A related set of quantities is frequently adopted for the description of yield surfaces of cohesive
frictional materials:
1
p= (3.17a)
3
v
t1
u h i
Æ 2 2
q= 3 · J2 = σ − σy + σy − σz + (σ − σz )2 + 3 τy
2 + τ2 + τ2
yz z
(3.17b)
2
1/ 3
1
r=3 J3 (3.17c)
2
A third set of widely used invariants (ξ, ρ, θ) describe a cylindrical ccordinate system, the
Haigh-Westergaard stress space:
1
ξ= (3.18a)
sqrt3
Æ
ρ = 2 · J2 (3.18b)
p
3 3 J3 π
cos (3θh ) = with 0 ≤ θh ≤ + (3.18c)
2 3/ 2 3
J2
(3.18d)
In this documentation, unless stated otherwise, for the Lode angle θ an alternative definition is
adopted:
p
3 3 J3 π π
sn (3θ) = − with − ≤θ≤+ (3.19)
2 3/ 2 6 6
J2
Angle of friction
The angle of friction is zero for most fine grained cohesive soils under undrained conditions.
Friction angles larger than 40 degrees are encountered rarely.
Note: A slope without cohesion world cannot be steeper than the material’s ultimate friction
angle.
Cohesion
The cohesion as well as the friction have to be clearly distinguished for drained and undrained
conditions. For fine grained soils a pore pressure is created for sudden loading which de-
creases the possible friction considerably. As time passes, the water will leave the soil, the
friction increases, but the cohesion will be reduced by a factor up to 10.
Dilatancy
Dilatancy denotes the volumetric plastic deformation behavior of a material sample under
shear. The quantity which is used to describe the dilatancy effects is called the dilatancy
angle ψ:
For soil materials, the plastic deformation behavior depends on the material’s effective density,
which in turn changes with the material’s loading state - therefore, the dilatancy angle is in fact
not a constant quantity. This coherence is described by the well-established stress dilatancy
theory (Rowe [4]), which links the mobilized dilatancy angle ψm to the actual shear strain-
ing level, the latter being characterized by the mobilized friction angle φ′m (cf. section NMAT
Hardening Plasticity Soil Model - GRAN).
In soil mechanics usual way to qualitatively illustrate the dilatancy effects of normally (NC) and
overconsolidated soil (OC) is to idealize the soil with an assembly of spherical particles (Fig.
3.3). In case of a loosely packed soil any shear deformation will result in a denser packing of
the particles, followed by a volume decrease, also know as the contraction (Fig. 3.3a). This
process is characterized by the negative angle of dilatancy ψ. On the other hand, shearing
of the initially densely packed soil is only possible if the volume increases (Fig. 3.3b). This
process is also know as dilatation and it is characterized by the positive dilatancy angle.
ψ>0
ψ<0
σ′ σ′
τ τ
ΔV < 0 ΔV > 0
Figure 3.3: Model of soil as spherical particles assembly: Contraction and dilatation under shear
The easiest way to understand the mechanical interpretation of the dilatancy angle is to look
at the behavior of the soil under a plane-strain conditions subjected to shear τy (Fig. 3.4).
′
σyy
τy
ϵ̇pyy ψ
γ̇py
1
Figure 3.4: Mechanical interpretation of the dilatancy angle ψ for a soil under a plane-strain condi-
tions
p p
where ϵ̇yy is the plastic strain rate in y-direction describing the change in volume and γ̇y is
the plastic shear strain rate.
The internal friction of a densely packed granular material is composed of a sliding friction
between the grains increased by an additional resistance as a result of the granular interlock-
ing. The effect of the dilatancy on the mechanical behavior of the soil, or more precisely on
the friction angle can be illustrated by the sawtooth model shown in Figur0e 3.5 (Maksimović
[5]). It is assumed that the friction angle of the non-compacted soil, when the volume gradient
is zero (constant volume), is a constant value and it describes only the friction resistance of
the grains. This angle is called critical state friction angle and usually denoted as φ′c . For
the densely packed soil as a result of the interlocking between the grains, the sheering does
not develop along a straight horizontal plane s − s, but along the sawtooth planes t − t with
the angle of inclination ψm , called the mobilized dilatancy angle. The friction angle mobilized
along the imaginary horizontal shear plane s − s, φ′m , can now be determined as the sum of
the friction angle along the inclined shear planes t − t , φ′c , and mobilized dilatancy angle, ψm ,
defining the inclination of t − t planes, i.e.
σ′
τ
t φ′m
ψm t
s s s s
ψm
t t φ′cs
ψm
For the triaxial stress state the dilatant soil behavior is more complicated than the one de-
scribed by the sawtooth mechanism but the overall behavior is qualitatively similar. In case of
the triaxial stress state, the relationship between the mobilized dilatancy angle and the plastic
strain rates can be described as follows
p
ϵ̇
sin ψm = , (3.22)
γ̇p
p p p p
were ϵ represents the plastic volumetric strain and γp = ϵ1 − ϵ2 − ϵ3 is the measure of the
plastic shear strain.
Figure 3.6a schematically shows the differences in the behavior of the normally consolidated
(NC) and heavily overconsolidated (OC) soils in drained triaxial experiment. As aforemen-
tioned, under shear normally consolidated soils exhibit contractive behavior, while the over-
consolidated soils have a tendency to expand (dilatation) after some small initial range of con-
traction. For NC soils deviatoric stress - axial strain curve will monotonically grow until the
shear strength qcs is reached. On the other hand, OC soils reach their peek shear strength
qpek at a maximal dilatation (maximal inclination of the ϵ − ε curve), followed by the reduc-
tion of the shear stiffness up to the remaining level qcs corresponding to the zero volumetric
strain gradient. This state is also know as the critical state.
q φ′m
qpek φ′pek
OC φ′cs
qcs
NC φ′ƒ
φ′μ
ϵ ϵ
ϵ ϵ
dϵ
ψm → 0
dϵ m
OC
contraction dilatation
max. dilatancy
ψpek
ϵ ϵ
ψmn
ψm = 0
NC
(a) Deviatoric stress q vs. axial strain ϵ and vol- (b) Mobilized friction angle φ′m vs. axial strain ε
umetric ϵ vs. axial strain ε for NC and OC soils and volumetric ϵ vs. axial strain ε for OC soils
Figure 3.6b shows the dependency of the mobilized friction angle ψ′m on the level of axial shear
strains for the OC soils. The value φ′ƒ marks the transition between the contractive and dilatant
behavior (a.k.a. phase transformation). Above this level of straining, the behavior of the soil
is purely dilatant. The value of φ′ƒ is in between the critical state friction angle φ′cs and the
grain-to-grain friction angle φ′μ . The friction angle range [φ′μ , φ′cs ] is not larger than a couple
of degrees (Wehnert [6]). For dense soils, the value of the φ′ƒ will be closer to the φ′μ . For
loose soils as well as for the case of the plain-strain conditions the value of φ′ƒ will be closer
to the upper limit, i.e. to the value φ′cs . Peak friction angle φ′pek corresponds to the maximal
dilatancy, i.e. ψpek .
Evidently the dilatancy is not a constant value - it depends on the stress level. Most advanced
soil models recognize this fact by incorporating models which can take into consideration stress
dependent dilatancy angle. One of the most famous dilatancy models is a model according to
Rowe [4]. This model establishes the relationship between the mobilized dilatancy angle ψm
Mobilized friction angle φ′m represents the measure of the stress and it can be computed as
σ1′ − σ3′
sin φ′m = . (3.24)
2c · cot φ − σ1′ − σ3′
ψm [ ◦ ]
ψpek
0 φ′m [ ◦ ]
φ′ƒ φ′pek
ψmn
Figure 3.7: Mobilized dilatancy angle - mobilized friction angle relationship according to Rowe [4]
Classical elasto-plastic material models (e.g. MOHR, DRUC) adopt a constant dilatancy angle.
If the dilatancy has considerable effects (e.g. due to arch action) the predominant loading
situation should be assessed carefully. As a coarse rule of thumb, for dense soils a value of
+φ/ 2, for middle dense soils a value of 0.0 and for loose soils a value of −φ/ 2 can be adopted.
In contrast, the advanced theoretical setting of the Hardening Plasticity Soil model directly
incorporates a loading state dependent variation of the dilatancy angle according to the above
mentioned stress dilatancy theory (cf. subsection 3.11.7).
Plastification number
Value of the corresponding yield function for the uncorrected (=linear elastic) stress state, pos-
sibly scaled to stress units. If >0 the material undergoes plastification. The value is computed
for each loading step anew. Therefore, regions that possibly have plastified previously, still can
get values < 0 in a subsequent loading step. For GRAN the value of the MOHR yield function
is computed, here.
For the ’Hardening Plasticity Soil’ (GRAN) material model, the plastification number is an iden-
tifier that holds more detailed information about the current state of loading (see section ’Hard-
ening Plasticity Soil Model’) - instead of storing the current value of the yield function, only.
rates. It reflects the volume neutral (shearing) portion of the plastic deformation:
v v
t2t 1
u u
ϵ̇p,de = ė2
p,
+ ė2
p,y
+ ė2
p,z
+ γ̇2
p,y
+ γ̇2
p,yz
+ γ̇2
p,z
(3.26a)
3 2
v
t2 1
u h
2 2 2 i
= ϵ̇p, − ϵ̇p,y + ϵ̇p,y − ϵ̇p,z + ϵ̇p, − ϵ̇p,z + γ̇2
p,y
+ γ̇2
p,yz
+ γ̇2
p,z
9 3
where
ϵp,
ėp, ≡ ϵ̇p, − , = , y, z (3.27a)
3
X 1h 2 2 2 i
⇒ ė2
p,
= ϵ̇p, − ϵ̇p,y + ϵ̇p,y − ϵ̇p,z + ϵ̇p, − ϵ̇p,z (3.27b)
=,y,z
3
Hint
For triaxial conditions with ϵ̇p,1 ≡ ϵ̇p,2 we have
2 1
ϵ̇p,de = ϵ̇p,1 − ϵ̇p,3 = ϵ̇p, − 3ϵ̇p,3 (3.28)
3 3
Assuming, that the plastic volumetric strain portion is negligible, this can be further re-
duced to
σ1′ − σ3′
sin φm = (3.30)
2c · cot φ − σ1′ − σ3′
Utilisation level
Ratio ≡ φm / φnp ≤ 1 , where the material input parameter φnp marks the maximum (ultimate)
friction angle (np - input) .
Application range
Metals and other materials without friction. Simulation of creep effects.
Parameters
kN/ m2
P1 Yield stress 1089 !
kN/ m2
P2 Hardening modulus (tangent modulus) 1088 0.0
P10 Type of creep law (overstress function) − 0
0 no viscous effects, elasto-plastic
1 Perzyna model
P11 Creep parameter, exponent m ≥ 1.0 − 1.0
kNs/ m2
P12 Viscosity η ≥ 0.0 1220 0.0
Formulation of the viscoplastic material behaviour is based on the Perzyna model. Accordingly,
the viscoplastic strains are defined by
Z t
p
Δϵ = ϵ̇p dτ
t0
Z t
∂g(σ, κ)
= λ̇ · dτ (3.32)
t0 ∂σ
Z t
〈(ƒ (σ, κ))〉 ∂g(σ, κ)
= · dτ
t0 η ∂σ
In case of an associative flow-rule (e.g. von Mise material) the plastic potential g equals the
yield function. The overstress function reads
ƒ (σ, κ)m , ƒ ≥0
〈〉 = (3.33)
0 , ƒ <0
This frequently used form can be calibrated to reproduce a wide range of time-dependent
material phenomena adequately.
References
The implementation of this material model essentially adopts concepts and strategies from
Chrisfield [7], Chrisfield [8], Zienkiewicz and Taylor [9] and Zienkiewicz and Cormeau [10].
Application range
Soil and rock with friction and/ or cohesion. Modelling of time-dependent effects (e.g. short
term strength)
Parameters
P1 Friction angle φ [ ◦] !
<0 inner cone
≥0 outer cone
kN/ m2
P2 Cohesion c 1089 !
kN/ m2
P3 Tensile strength ƒt 1089 0.0
P4 Dilatancy angle ψ [ ◦] 0.0
P5 unused - -
P6 Plastic ultimate strain ϵp,de,
h
1081 0.0
P7 Ultimate friction angle φ [ ◦] P1
kN/ m2
P8 Ultimate cohesion c 1089 P2
P10 Type of creep law (overstress function) − 0
0 no viscous effects, elasto-plastic
1 Perzyna model
P11 Creep parameter, exponent m ≥ 1.0 − 1.0
kNs/ m2
P12 Viscosity η ≥ 0.0 1220 0.0
2 sin φ Æ 6c cos φ
ƒ=p · 1 + J2 − p ≤0 (3.34a)
3 (3 ± sin φ) 3 (3 ± sin φ)
2 sin ψ Æ
g= p · 1 + J2 (3.34b)
3 (3 ± sin ψ)
This formulation describes a cone in principal stress space that either embraces the MOHR
yield surface (- sign) or is embraced by and tangential to it (+ sign).
For a description of the material’s viscoplastic extension, please refer to NMAT VMIS.
Tension cut-off
By default, the Drucker-Prager yield surface is accomplished by a tensile stress limit of R ANK -
INE style
ƒ = σ ,, − ƒt ≤ 0 (3.35)
where σ ,, denotes the vector of principal stresses and ƒt the entered tensile strength.
2 sin φ 6 cos φ
ƒ= · 3p + q − ·c ≤0 (3.36)
(3 − sin φ) (3 − sin φ)
The DRUCKER/ PRAGER yield surface features two strength parameters φ and c which can
be derived from two failure conditions, one for the uniaxial and one for the biaxial compression
limit state, respectively.
Given the uniaxial compressive strength ƒc , the corresponding failure stress state renders
−ƒc
1
σ ,m = 0 ⇒ p=− ƒc ; q = ƒc (3.37)
3
0
For many materials, in particular concrete, an increased strength ƒc2c = · ƒc can be observed
in biaxial compression tests, cf. Kupfer and Gerstle [11]. In this notation we have introduced
the amplification factor . The corresponding biaxial stress state at failure thus renders
− · ƒc
2
σ b,m =
− · ƒc
⇒ p=− · ƒc ; q = · ƒc (3.38)
3
0
! !
The corresponding failure conditions ƒ σ ,m = 0 and ƒ σ b,m = 0 constitute a set of two
equations, which can be solved for the two material parameters φ and c to be determined.
sin φ (3 − sin φ)
!
− + · ƒc − c =0 (3.39a)
3 cos φ 6 cos φ
sin φ (3 − sin φ)
!
−2 + · · ƒc − c = 0 (3.39b)
3 cos φ 6 cos φ
Hint
In the absence of experimental data, Fédération internationale du Béton [12] proposes
the following relation for the biaxial strength ƒc2c
ƒc [MP]
ƒc2c = 1.2 − ·ƒc (3.41)
1000
| {z }
Assuming normal strength concrete, ≈ 1.16. Consequently, from eq. 3.40a and eq.
3.40b:
References
The implementation of this material model essentially adopts concepts and strategies from
Zienkiewicz and Taylor [9] and Chrisfield [8].
Application range
Soil and rock with friction and/ or cohesion. Modelling of time-dependent effects (e.g. short
term strength).
Parameters
P1 Friction angle φ [ ◦] !
Table continued on next page.
kN/ m2
P2 Cohesion c 1089 !
kN/ m2
P3 Tensile strength ƒt 1089 0.0
P4 Dilatancy angle ψ [ ◦] 0.0
P5 unused - -
P6 Plastic ultimate strain ϵp,de,
h
1081 0.0
P7 Ultimate friction angle φ [ ◦] P1
kN/ m2
P8 Ultimate cohesion c 1089 P2
P10 Type of creep law (overstress function) − 0
0 no viscous effects, elasto-plastic
1 Perzyna model
P11 Creep parameter, exponent m ≥ 1.0 − 1.0
kNs/ m2
P12 Viscosity η ≥ 0.0 1220 0.0
1 sin θ sin φ
Æ
ƒ = 1 sin φ + J2 cos θ − p − c cos φ ≤0 (3.43a)
3 3
2 sin ψ Æ
g= p (3 ± sin ψ) · 1 + J2 (3.43b)
3
Special comments
The following expressions are better suited for checking the yield criterion:
1 − sin φ 2c cos φ
ƒ = σ1 − · σ3 − (3.44)
1 + sin φ 1 + sin φ
For a description of the material’s viscoplastic extension, please refer to NMAT VMIS.
References
The implementation of this material model essentially adopts concepts and strategies from
Zienkiewicz and Taylor [9], Chrisfield [8] and Zienkiewicz and Taylor [9].
Application range
Realistic stiffness and hardening behavior of soil, settlement analysis.
Parameters
P1 Friction angle φ [ ◦] !
kN/ m2
P2 Cohesion c 1089 !
kN/ m2
P3 Tensile strength ƒt 1089 0.0
P4 Dilatancy angle ψ [ ◦] 0.0
kN/ m2
P5 Stiffness modulus Es,reƒ (GRAN-extended) 1088 *
P6 Lateral earth pressure coefficient K0nc (GRAN-extended) − 1 − sin φ
kN/ m2
P9 Modulus for primary loading E50,reƒ 1088 !
P10 Exponent m ≥ 0.0 − 0.7
P11 Failure factor 0.5 < Rƒ < 1.0 − 0.9
kN/ m2
P12 Reference pressure preƒ > 0.0 1089 100.0
The extended version of the GRAN-model (two-surface model, double hardening) is activated
by specification of the oedometric stiffness modulus Es,reƒ (P5) - only in this case the lateral
earth pressure coefficient K0nc (P6) takes effect. In case no input of Es,reƒ is provided, the basis
version of the GRAN material model (single-surface model, single hardening) is adopted.
The hardening rule is based on the hyperbolic stress-strain relationship proposed by Kondner
and Zelasko [13], which was derived from triaxial testing. Hardening is limited by the mate-
rial’s strength, represented by the classic MOHR/COULOMB failure criterion. Additionally, the
model accounts for the stress dependent stiffness according to equations (4-6). A further es-
sential feature is the model’s ability to capture the loading state and can therefore automatically
account for the different stiffness in primary loading and un-/reloading paths.
In the subsequent notation, compression and contraction are defined as negative; for the prin-
cipal stresses the relation σ1 ≥ σ2 ≥ σ3 holds. Accordingly, for the triaxial state index 3 denotes
the axial and index 1 the lateral direction.
q = σ1 − σ3 [kN/ m2 ]
-1200
-1000
σ1 -800
← MOHR
-600
-400
-200
σ3
-1000 -800 -600 -400 -200 0
Mohr-Coulomb
σ2 p [kN/ m2 ]
q = σ1 − σ3 [kN/ m2 ]
-1200
-1000
σ1 -800
← MOHR
-600
-400
-200
σ3
-1000 -800 -600 -400 -200 0
Mohr-Coulomb
σ2 p [kN/ m2 ]
The extended version (GRAN-extended) enhances the model by an additional hardening two-
parameter cap surface. An appropriate calibration of the cap’s hardening and shape param-
eters is done automatically, based on the input of physically sound input parameters - and,
hence, allows for
• a realistic modelling of the contractant behaviour and stiffness during primary compression
(oedometric testing) ⇒ plastic straining
σter
• preservation of a realistic stress ratio K0nc = σ , e.g. according to Jaky as K0nc = 1− sin φ
parameter : Es,reƒ ; K0nc ; (m; preƒ )
−ϵ3
q = σ1 − σ3 = , (3.45)
b − · ϵ3
with
1
= E ≃ 2 · E50 , (3.46a)
b
1 qƒ
= q = , (3.46b)
Rƒ
where q is the asymptotic value of the deviatoric stress and E is the initial stiffness. Duncan
and Chang [14] extended this formulation by introducing the limit value to the deviatoric stress
q in form of the ultimate deviatoric stress qƒ = q · Rƒ . The ratio between qƒ and q is given by
the failure factor Rƒ , which is smaller than 1.0.
1
E =
1 b
q =
E50
qƒ
σ1 − σ3
E
ε3
σ3 [kP]
σ reƒ Ereƒ
s
1.0
ε3 [−]
m
|σ3 | · sin φ + c · cos φ
Es = Es,reƒ · (3.47)
preƒ · sin φ + c · cos φ
In contrast to oedometric test conditions, lateral expansion is not constrained under triaxial
conditions. Due to the changed boundary conditions the triaxial modulus’ stiffness E50 deviates
from the stiffness modulus. E50 is defined as secant stiffness that corresponds to a 50-percent
mobilisation of the maximum shear capacity (figure 1). Choosing the smaller compressive
stress σ1 as reference stress, a relation anlogous to equation (4) can be established for the
stiffness evolution of the triaxial modulus E50 (Kondner & Zelasko 1963, Duncan & Chang
1970), which is then used in the model equations (1) to (3).
m
|σ1 | · sin φ + c · cos φ
E50 = E50,reƒ · (3.48)
preƒ · sin φ + c · cos φ
For the Hardening Plasticity Soil model, the so-called mobilized dilatancy angle ψm is defined
from the well-established stress dilatancy theory (Rowe [4]) as (Fig. 3.12)
see also subsubsection 3.11.2. Therein, the critical state friction angle φcs marks the transition
between contractive (small stress ratios with φm < φcs ) and dilatant (higher stress ratios with
φm > φcs ) plastic flow. The mobilized friction angle φm in Eq. 3.51 is computed according to
σ1′ − σ3′
sin φm = (3.52)
2c · cot φ − σ1′ − σ3′
At failure, when φm ≡ φ , also the dilatancy angle reaches its final value ψm ≡ ψ . Accordingly,
from Eq. 3.51 the critical state friction angle can be derived as
sin φ − sin ψ
sin φcs = (3.53)
1 − sin φ sin ψ
SOFiSTiK performs the computation of the critical state friction angle φcs automatically on
basis of the user specification for the final angles φ and ψ .
It has been recognized that in some cases the Rowe’s model for dilatancy angles (Eq. 3.51)
can overestimate the contractive behavior of the soil at low mobilized friction angles, φm < φcs .
As a remedy, several researchers have developed modified formulations based on the original
Rowe’s model. Some of these models which are implemented in SOFiSTiK are described
below.
One of the models which do not require additional input parameters is the model according to
Søreide [15] which modifies the Rowe’s formulation by using the scaling factor sin φm / sin φ
(Fig. 3.12)
This method is activated with the TALPA/ASE commands TALPA: CTRL MSTE EMAX 2 and
ASE: CTRL MSTE W4 2.
Wehnert [6] proposed a model based on a lower cut-off value ψ0 for the modification of the
Rowe’s formulation from Eq. 3.51 at low mobilized friction angles (Fig. 3.12)
sin ψ0 ; 0 < ψm ≤ ψRoe
m
sin ψm = sin φm − sin φcs . (3.55)
; ψRoe
m
< ψm ≤ ψ
1 − sin φm sin φcs
ψm [ ◦ ]
Rowe
Soreide
20 Wehnert, ψ0 = −3◦
Constant
ψ = 10
φcs
0 φm [ ◦ ]
5 15 25 φ = 35
ψ0
-10
-20
Figure 3.12: Comparison of models for mobilized dilatancy angle ψm implemented in SOFiSTiK for
φ = 35◦ and ψ = 10◦
Hint
It should be noted that for saturated soils under undrained conditions, the evolution of ex-
cess pore water pressures shows a pronounced sensitivity regarding the volumetric defor-
mation behavior of the soil phase. As a consequence, when carrying out corresponding
simulations using the 2-phase model, the dilatancy model parameters should be carefully
selected. Ideally, model parameters should be calibrated according to laboratory tests.
For more information regarding the effect of different modeling approaches on the
volumetric deformation behavior of soil under undrained conditions, also refer to
the documented benchmark “Triaxial Consolidated Undrained (CU) Test” in VERIFICA-
TION_MANUAL_MECHANICAL: BE48.
Ratio ≡ φm / φnp ≤ 1 , where φnp is the maximum (ultimate) friction angle provided as
material input parameter.
• Isotropic pre-consolidation stress
Hydrostatic stress pc marking the highest state of compression that was reached in loading
history (cap); the hydrostatic stress being defined as p ≡ σ + σy + σz / 3 .
Special comments
The model can easily be calibrated according to triaxial/oedometric test data. Therefore, de-
formation behaviour of the material prior to failure can be captured with a good accuracy. This
feature, combined with the consideration of specific stiffnesses for primary and un-/reloading,
respectively, constitutes a significant progress when compared to the behaviour of classic
elasto-plastic soil material models. Consequently, GRAN is particularly suited for tasks that
require more precise settlement predictions.
If no precise data is available, then the following estimations may be used for an approximation
of the properties of normally consolidated soil:
References
The implementation of this material model essentially adopts concepts and strategies from
Schanz [16], Benz [17], Kondner and Zelasko [13], Desai and Christian [18] and Duncan and
Chang [14].
Parameters
P1 Swelling modulus Kq
h
1081 3.3
kN/ m2
P2 Swelling limit stress (absolute value) |σc | > 0 1089 10.0
(magnitude of smallest compressive stress below which no
further increase of swelling occurs)
kN/ m2
P3 Historical swelling equilibrium stress (absolute value) from 1089 2000.0
oedometer testing |σ0,hst | > |σc |
P4 Viscous extension: retardation time η ≥ 0.0 h 0.0
where
1..3
σ principal normal stresses
σ0 equilibrium state of stress w.r.t. swelling (initial condition), transformed to the
direction of the principal normal stresses σ
Special comments
Swelling of soils is a complex phemomena that is influenced by various factors. There are two
swelling mechanisms of practical importance that can be distinguished - for both processes
the presence of (pore-) water is a common prerequisite. The first mechanism is termed as the
”osmotic swelling” of clay minerals, which basically is initiated by unloading of clayey sedimen-
tary rock. The second mechanism takes place in sulphate-laden rock with anhydride content.
In this case the swelling effects are due to the chemical transformation of anhydride to gypsum-
which goes along with a large increase in volume (61%).
For both described mechanisms a principal dependency between the increase in volume,
caused by swelling, and the state of stress was observed both in laboratory and in in-situ
experiments. The formula employed represents a generalization of the 1-dimensional stress-
strain relationship that HUDER and AMBERG derived from oedometer tests for the final state.
Here the time dependent evolution of the swelling process is not considered.
The equilibrium stress state with respect to swelling σ0 is defined by means of the GRP record.
For this we use the option PLQ in order to reference a (previously calculated) load case as
”primary state for swelling”. This state is regarded as an equilibrium state with respect to
swelling (normally in-situ soil prior to construction work). I.e. swelling strain increments caused
by an eventual ”unloading” from the historical equilibrium state σ0,hst to this new ”primary
state” σ0 have already occurred. Swelling strain increments in the course of construction work
are only due to unloading related to the new ”primary state for swelling” σ0 :
Upper Limit
Test Data
Grob (Eq.2)
1.0
σc = −2 kN/ m2
0.9
0.8
0.7
0.6
0.5 εq [%]
0.4
kq = 0.33 %
0.3
0.2
0.1
-800 -700 -600 -500 -400 -300 -200 -100 0 100 200
σ [kN/ m2 ]
q q q
(Δ) ϵ∞ = ϵ,tot − ϵ,hst
σ σ0
§ ª § ª
= −p1 · log − −p1 · log
σ0,hst σ0,hst (3.58)
σ
= −p1 · log
σ0
The constitutive equation reproduced above is limited to the final (stationary) state, i.e. it relates
the evolved swelling strains to the stress state that is present at time t = ∞ . To account for time
dependent behaviour, the relation is extended to the time scale by a formal viscous approach.
Correspondingly, the rate of swelling strains is defined as
q
ϵ∞ (σ)−ϵq
ϵ̇q = (3.59)
η
with the retardation time η as a viscosity parameter and ϵq denoting the swelling strains that
have developed at the considered time t . In rheological terms this approach can be interpreted
as a parallel coupling of a ’swelling’ and a dashpot device.
The time dependent response can be calibrated via the retardation time η (P4) - the greater
η the more accentuated is the retardation in the evolution of swelling strains. For η = 0 the
response is instantaneous, identical with the non-viscous (instationary) case. Furthermore,
for t → ∞ the model’s response converges to the instationary solution - independent of the
adjusted retardation time η . This property enables application of the viscous model also for
stabilisation of the solution process, even if one is not explicitly interested in modelling time
effects.
The SWEL record is specified in addition to a linear elastic or elastoplastic basic material.
References
The implementation of this material model essentially adopts concepts and strategies from:
Wittke-Gattermann [19], Wittke [20], Wittke and Rissler [21], Baugrund-Institut [22], and
Zienkiewicz and Taylor [9].
Parameters
The shear stress over the designated shear plane and corresponding plastic strain develop-
ment is controlled by means of the yield condition ƒ1 and plastic flow rule g1 according to
ƒ1 = tn φs · σ − cs + τ ≤ 0 (3.60a)
g1 = tn ψs · σ + τ (3.60b)
where:
The tensile stress normal to the plane is limited by the additional cut-off criterion
ƒ2 = σ − ƒt,s ≤ 0 (3.61)
where again:
Special comments
This material law may be specified up to three times in addition to the material law of the base
material (elastic, MOHR, DRUC). This allows for the consideration of different distinct fault
directions. Increasing the number of specified shear planes per material, also increases the
number of possible equilibrium states for a material point - this may possibly affect the stability
of the overall equilibrium iteration process.
Specification of meridian angle OAL and descent angle OAF follows the instructions given
in the descriptions for input records MAT/MATE. For planar systems the value OAL directly
defines the slope of the stratification, i.e. the angle between the local x direction and the global
X direction. Input for OAF is not evaluated for the plane case.
For P9>0 a scalar damage model with exponential softening of the tensile strength is applied.
The softening obeys
§ ƒt,s ª
ƒt,s = ƒt,s · ep − · (3.62)
Gƒ
where denotes the crack opening. In this context, the tensile fracture energy Gƒ represents
an objective material parameter. In order to minimize discretization dependent spurious side
effects, a characteristic element size is incorporated into the softening formulation. This re-
quires, however, a sufficiently fine finite element discretization in the corresponding system
domains.
In case of P9=0 a tension cut-off with respect to ƒt,s without consideration of softening is
executed.
References
The implementation of this material model essentially adopts concepts and strategies from
Chrisfield [8] und Zienkiewicz and Taylor [9].
Parameters
kN/ m2
P8 Tensile capacity ƒt,r , rock 1089 0.0
P9 Dilatancy angle ψr , rock [ ◦] 0.0
Shear plane:
where:
Special comments
The mpdel allows to specify the strength of the rock as well as the capacity of a predominant
shear plane, which is are defined by the angle P5. The flow rule of the hear failure is non-
associated if P4 is different from P1.
References
The implementation of this material model essentially adopts concepts and strategies from
Chrisfield [8] and Zienkiewicz and Taylor [9].
Application range
For soils saturated with water with low permeability and high loading rate and when the short
term behavior has to be assessed.
Parameters
Undrained Soil
Soil is a porous material consisting of solid particles which together represent the grain skele-
ton. The pores can be filled with fluid (commonly water) and air. Hence we talk about soil as a
three phase system (Figure 3.14).
V Air n
Vn n
V Water n
V 1
Solid
Vs 1−n
Particles
The solid particles occupy the volume Vs . An important basic parameter of soil is the porosity
n, which is defined as the ratio of the volume of the pores Vn and the total volume of the soil V
(Figure 3.14b), i.e.
Vn
n= . (3.65)
Vs
As aforementioned, a part of pore volume is usually occupied by water V while the rest is filled
with air V . When the volume of air in pores approaches zero, i.e. V → 0, the water completely
occupies the pore space (Vn = V ) and we speak of soil fully 1 saturated with water.
For the following, unless stated otherwise, we restrict ourselves to the soils saturated with
water.
Let us consider an element of soil saturated with water subjected to some external total force
P acting on a plane with area A (Figure 3.15). These external forces are transmitted by a
pore water pressure p and by stresses in the particles. The stresses in the particles are
a result of the the concentrated forces acting on the contact points between the grains and
the pressure in water encompassing the grains (Figure 3.15a). Assuming that the particle
stiffness is much higher than the stiffness of the soil skeleton as a whole, it follows that the
soil element can deform only by sliding and rolling in the contact points between the particles.
Hence the deformation of the soil skeleton depends solely on the concentrated forces acting
1 As will be illustrated below, a complete saturation with V = 0 is a theoretical concept rather than a realistic
option. When we use the terminus saturated soil, it should be thought of as quasi-saturated.
P
P
s s As ≈ A
s s
p
N′
Having aforementioned in mind, inside the soil element we can construct a surface s − s which
crosses trough the points of contact between the soil particles (Figure 3.15a). Due to the
relatively small size of soil particles, the s − s surface does not differ much from the plane and
hence it has an area As which is approximately equal to A. The contact forces acting on s − s
can be resolved into tangential and normal component N , acting in the direction of the force P.
Assuming that the contact areas between the soil particles are practically points (their area is
negligible compared to A), the pressure in the water p will also act on the surface s − s over
the entire area As ≈ A. We can now write the equilibrium equations for the surface s − s in the
vertical direction (Figure 3.15b) as follows:
X
P= N′ + p · A . (3.66)
P ′
Dividing Eq. 3.66 with A and introducing the notation σ = P/ A and σ ′ = N / A we get:
σ = σ ′ + p . (3.67)
Eq. 3.67 represents the basis for one of the fundamental principles of soil mechanics pos-
tulated by Karl Terzaghi - the principle of effective stress. This principle states that the total
normal stress σ comprises of effective normal stress σ ′ and pore water pressure p . The
effective normal stress σ ′ is a measure of contact forces of granular material and hence, under
the assumptions stated above, it governs the deformation of soils.
Hint
The Therzaghi’s effective stress principle is based on three main assumptions:
For standard saturated soil these assumptions are valid; however for rock some of these
assumptions may not hold. In this case the compressibility of rock must be taken into
account and therefore Eq. 3.67 should be modified to take this fact into an account.
The Terzaghi’s principle generalized for a three-dimensional body reads: the total normal
stresses in a saturated soil consist of two parts: effective normal stresses and pore (water
or fluid) pressures (Terzaghi [23]). The effective shear stresses are equal to the total shear
stresses since the resistance of pore water (fluid) to shear is neglegible. This can be described
by the following equations:
where:
Writing the Eqs. 3.68 in matrix form, the total stress increments can be expressed as
Δσ = Δσ ′ + Δσ , (3.69)
where:
h iT
Δσ ′ = Δσ′ Δσy′ Δσz′ Δτy Δτz Δτyz , (3.70a)
h iT
Δσ = Δp Δp Δp 0 0 0 . (3.70b)
Δσ ′ = D′ · Δϵ , (3.71)
where D′ denotes the elasticity matrix of the soil skeleton, as described in MATE.
Pore stress comprises of two parts: steady-state pore stress ps and excess pore stress pe :
The rate of change of the steady-state part of the pore stress Δps is equal to zero and hence
the rate of change of pore pressure Δp is equal to the rate of change of excess pore pressure
Δpe , i.e.
In the undrained conditions the relative movement between the water (fluid) and soil skeleton
is negligible and therefore the strains are the same for each phase of the soil. Hence, for the
water stress increments we can write:
Δσ = D · Δϵ , (3.74)
with D representing the pore water stiffness matrix which is given by:
1 1 1 0 0 0
1 1 1 0 0 0
1 1 1 0 0 0
D = Ke · . (3.75)
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
Ke is a constant and represents the equivalent bulk modulus of the pore fluid (water). Sub-
stituting Eq. 3.75 into Eq. 3.74, and having in mind Eqs. 3.70b and 3.76, we can write the
relationship between the rate of excess pore water pressure Δpe and the rate of volumetric
strain Δϵ
The modulus Ke is related to the bulk modulus of the pore water K . This relationship is
examined next.
In an elementary volume V of saturated soil under undrained conditions the pore water occu-
pies a volume V = n · V and the solid soil particles occupy a volume Vs = (1 − n) · V , where n
is the porosity (Figure 3.14) . The total volume change in the soil ΔV is equal to the change of
volume in the soil particles ΔVs plus the change of volume in the pore water ΔV , i.e.
Since the soil particles are assumed to be incompressible (ΔVs ≈ 0), the change of the volume
V can take place only on the account of the change of porosity. Since the soil is saturated with
water and under undrained conditions, the change of porosity is only possible if the water in
the pores changes its volume. Taking this into an account and dividing Eq. 3.77 with V , the
total volume change per unit volume of soil Δϵ reduces to:
ΔV ΔV
Δϵ ≈ = ·n . (3.78)
V V
The change of volume in the pore water per unit area of pore water, ΔV / V , is directly pro-
portional to the change of pore water pressure Δpe and inversely proportional to the bulk
modulus of pore water K . Now the volumetric strain of soil can be written as
Δpe
Δϵ ≈ ·n . (3.79)
K
Comparing Eqs. 3.76 and 3.79 the relationship between the equivalent bulk modulus of the
pore water Ke and the bulk modulus of the pore water K can finally be established Potts and
Zdravković [24]:
K
Ke ≈ . (3.80)
n
Writing the effective stress principle given by the Eq. 3.69 in terms of the isotropic stress
p = (σ1 + σ2 + σ3 )/ 3, and taking into consideration Eqs. 3.71 and 3.74, we get
K K
Δp = Δp′ + Δpe = K ′ Δϵ + Δϵ ⇒ K = K ′ + , (3.81)
n n
where K ′ = Δp′ / Δϵ is the bulk modulus of soil skeleton, while K = Δp/ Δϵ represents the
undrained total bulk modulus of soil.
The equivalent pore water bulk stiffness K / n can be defined by the bulk modulus of the
soil skeleton K ′ and the Skempton’s B-parameter, which is convenient since the Skempton’s
B-parameter can be measured experimentally in the undrained triaxial test. It is for this rea-
son that the Skempton’s B-parameter is chosen as an input parameter (parameter P1) for the
undrained material properties 2 .
In the undrained triaxial compression test where the soil sample is subjected to the isotropic
compression (Δσ1 = Δσ2 = Δσ3 ), Skempton’s B-parameter represents the relationship be-
tween the change of pore pressure Δp and the change of isotropic total stress Δp =
(Δσ1 + Δσ2 + Δσ3 )/ 3, i.e.
Combining Eqs. 3.82 and 3.81 we can finally obtain the relationship between the equivalent
bulk modulus of the pore water and the bulk modulus of the soil skeleton K ′ trough Skempton’s
B-parameter3
K B
Ke = = · K′ . (3.83)
n 1−B
Theoretically, for soils fully saturated with water, K is equal to the bulk modulus of pure water
K0 . In many applications the water is assumed to be an incompressible fluid (K 0 = ∞), which
would then imply that the saturated soil is also incompressible (K = ∞). This corresponds to
the value B = 1, which also marks the upper limit for the Skempton’s B-parameter. This limit
value describing fully incompressible behavior of the pore water ( K / n → ∞, when B → 1; Eq.
3.83) and saturated soil cannot be reached within a FE calculation due to the singularity of the
stiffness matrix which then occurs.
2 Skempton [25] has expressed the pore pressure change Δpe occurring under changes in total principle
stresses Δσ1 and Δσ3 by the equation: Δpe = B · [Δσ3 + A · (Δσ1 − Δσ3 )]. He termed the coefficients A and B
as “Pore-pressure Coefficients“, also know as “Skempton’s A and B Coefficients“.
3K / n = B · K .
However, fully incompressible behavior of the water is not realistic - water has a very large but
finite bulk modulus K 0 = 2 · 106 kN/ m2 , and hence the bulk stiffness of the pore water K / n
is always finite, even if the soil is fully saturated. It follows that Skempton’s B-parameter in real
applications is always smaller than one.
Hint
Moreover, full saturation is also not feasible in reality. Practically, it is not possible to
squeeze out the pore air volume V completely. Since the ’stiffness’ of air is considerably
smaller than the bulk stiffness of pure water K 0 , even a very small fraction of remaining
pore air volume has a significant impact; the resulting equivalent bulk modulus of the pore
fluid K / n shows very high sensitivity with respect to the degree of saturation S = V / Vn .
This relationship can be approximated by the following expression (Verruijt [26]):
0 ·p
K
1 S 1−S K 0 1
= + ⇒ Ke = = · , (3.84)
K 0
K p0 n (1 − 0
S) · K + S · p0 n
where p0 is the absolute pressure in the water, considered with respect to vacuum (p0 =
100kN/ m2 ).
For illustration: despite an almost complete saturation of 99.9%, i.e. S = 0.999, this
0 / K ≈ 20 – the compressibility of the pore fluid is about 20
relation indicates a ratio of K
times higher than the one of pure water!
To avoid potential numerical problems associated with very high stiffness values in FE calcula-
tions, the following upper limits are used in the program:
K / n 500
≤ 500 ⇒ B≤ ≈ 0.998 . (3.85)
K′ 501
On the other hand, in order to have the realistic computational results (B ≈ 1), the pore water
stiffness K / n must be substantially larger than the stiffness of the soil skeleton K ′ . Therefore
the lower limits are set to
K / n 10
≥ 10 ⇒ B≥ ≈ 0.909 . (3.86)
K′ 11
Default values for bulk stiffness of pore water and Skempton’s B-parameter are
K / n 45
= 45 ⇒ B= ≈ 0.978 . (3.87)
K′ 46
Very often, the stiffness of the pore water K / n is expressed in terms of the undrained Pois-
son’s ratio ν , the bulk modulus of the soil skeleton K ′ and the effective Poisson’s ratio ν ′ :
4
K 3(ν − ν ′ )
= · K′ . (3.88)
n (1 − 2ν )(1 + ν ′ )
It follows from Eqs. 3.88 and 3.83 that the undrained Poisson’s ratio and the Skempton’s B-
parameter are not independent of each other. They are connected by the following expressions
4E
= G · (1 + ν ) , K = E / [3(1 − 2ν )]
3ν ′ + B(1 − 2ν ′ ) 3(ν − ν ′ )
ν = ⇔ B= . (3.89)
3 − B(1 − 2ν ′ ) (1 + ν )(1 − 2ν ′ )
Once the bulk stiffness of pore water K / n is know so is the pore water stiffness matrix D ,
and the total stresses increments can be related to strains:
Δσ = (D′ + D ) · Δϵ . (3.90)
Application range
All materials with friction including rock and concrete (module ASE, 3D solid elements)
Parameters
P1 Parameter η1 − !
P2 Exponent m − !
kN/ m2
P3 Uniaxial tensile strength ƒt 1089 0.0
P4 Parameter η2 for flow rule − P1/10
kN/ m2
P5 Compressive strength (cap) 1089 -
P6 Tensile failure strain ϵt
h
1081 -
p m
ƒ= 13 − 27 + η1 · · 3 ≤ 0 (3.91a)
1
p m
g = 13 − 27 + η2 · · 3 (3.91b)
1
with
Special comments
Due to their non-physical nature, calibration of the LADE model parameters may not be obvious
at first sight. For this reason, the basic procedure for a material with known uniaxial tensile and
compressive strength (e.g. concrete) is described in the following. Of particular interest is
the section of the three-dimensional yield surface with one of the principal planes (-> ”Kupfer
Curve”).
Parameter m (exponent) affects the curvature (convexity) of the yield surface towards the hy-
drostatic axis - the larger m the stronger the curvature. In this manner m determines the
shape of the intersection curve. For most types of concrete a value of m between 1.0 and 2.0
is reasonable.
Using the known quantities of uniaxial tensile and compressive strength and the chosen pa-
rameter m, parameter η1 can now be determined from the condition: For the stress state
corresponding to the uniaxial compressive stress limit the yield condition must be fulfilled.
13 |1 | m
!
η1 = −27 · (3.92)
3 p
σ = σ = −ƒt
σ = − (ƒt + ƒc )
(3.93)
1 = σ + σ + σ
3 = σ · σ · σ
Where ƒt (= P3 ) and ƒc are the magnitudes of the uniaxial tensile and compressive strength,
respectively, 1 and 3 the required invariants for this stress state. Substituing into the rewritten
yield function yields the yet unknown parameter η1 .
The following table contains exemplary parameters for selected concrete types, derived from
the procedure described above (classification according to EC2, Ultimate Limit State).
By specification of parameter P5 the model can optionally be extended by a spherical cap (in
principal stress space) that limits the volumetric compressive stress to a maximum value. This
can be meaningful in particular for mainly hydrostatic compression. The cap is defined by:
q q
ƒ= σ1 2 + σ2 2 + σ3 2 − P5 2 + P5 2 + P5 2 ≤ 0 (3.94)
The optional parameter P6 defines a linear decay of the material’s tensile capacity after the
initial tensile strength ƒt has been exceeded, such that ƒt (ϵt ) = 0. If not specified, the behavior
is ideally plastic, the tensile strength remains at the initial value ƒt (tension cut-off). The tension
model is suitable for monotonic loading, only.
References
The implementation of this material model essentially adopts concepts and strategies from
Lade [27].
P1 and P2 should not be used for standard membranes. P1 can be used to limit the maximum
tension force, e.g. for geo textiles.
The failure of a membrane under pressure is generally activated in ASE via SYST...NMAT JES.
For linear elastic membrane materials only NMAT MEMB must be set in AQUA. Examples for
lineare membrane materials siehe [Link]\english\membranes\[Link].
With P3 and P4 the nonlinear yarn warp-fill material law according to Galliot and Luchsinger
[28] can be activated. P3 (ΔE ) and P4 (ΔEƒ ) define the nonlinear variation of the stiffness
in warp and fill direction depending on the stress ratio σ and σƒ (warp fill) . Example see
[Link]\english\membranes\nonlinear_warp_fill_behaviour.dat.
Please notice that in input MAT the E modulus must be input in kN/m2 but the values NMAT
MEMB P1+P3+P4 are in kN/m, see [Link].
σ σƒ
γ = q γƒ = q (3.95)
σ 2 + σƒ 2 σ 2 + σƒ 2
Stress-strain relation:
1 −νƒ
ϵ E (γ ) E (γ ) σ
= −νƒ
(3.96)
1
ϵƒ E (γ ) Eƒ (γƒ ) σƒ
1
Eƒ (γƒ ) = ΔEƒ γƒ − p + E1:1 ƒ (3.98)
2
For the advanced user the modules TALPA (for QUAD-elements) and ASE (for BRIC-elements)
offer the possibility to plug in self-developed non-linear material models via an interface (cur-
rently only for WINDOWS-OS). The following paragraphs describe the interface in detail.
The user-defined material models have to be provided in a Dynamic Link Library (DLL) with
arbitrary name. The variable SOFISTIK_USERMATDLL must be set with the name of this DLL.
This can either be done by specification of
SET SOFISTIK_USERMATDLL=my_material
SOFISTIK_USERMATDLL=my_material
into the [Link] file. In both cases the user defined material models, in the DLL
my_material.dll, are loaded at run-time. The interface routine itself reads:
Input parameters:
Return values:
BMAT
NO Material number − 1
kN/ m3
C Elastic constant normal to surface Cs 1097 0.
kN/ m3
CT Elastic constant tangential to surface Ct 1097 0.
kN/ m2
CRAC Tensile cracking strength of interface 1089 0.
kN/ m2
YIEL Maximum stress of interface 1089 -
MUE Friction coefficient of interface − -
kN/ m2
COH Cohesion of interface 1089 -
DIL Dilatancy coefficient − 0.
t/ m2
GAMB Equivalent mass distribution 1184 0
BMAT defines for an existing material (e.g. MATE/CONC/BRIC) properties for elastic support.
For a QUAD element it is thus possible to select the properties of the plate and the soil by a
single material number. The support stiffness may be calculated from the elasticity constants
of MREF, a geometric dimension H and a specific geometry rule. The bedding approach works
according to the subgrade modulus theory (Winkler, Zimmermann/Pasternak). It facilitates the
definition of elastic supports by an engineering trick which, among others, ignores the shear
deformations of the supporting medium. The bedding effect may be attached to beam or plate
elements, but in general it will be used as an own element. (see SPRI, BOUN, BEAM or QUAD
and the more general description of BORE profiles)
The determination of a reasonable value for the foundation modulus often presents consider-
able difficulty, since this value depends not only on the material parameters but also on the
geometry and the loading. One must always keep this dependance in mind, when assessing
the accuracy of the results of an analysis using this theory.
The subgrade parameters C and CT will be used for bedding of QUAD elements or for the de-
scription of support or interface conditions. A QUAD element of a slab foundation will thus have
a concrete material and via BMAT the soil properties attached to the same material number.
The value C is than acting in the main direction perpendicular to the QUAD surface in the local
z-direction, while CT is acting in any shear direction in the QUAD plane.
If subgrade parameters are assigned to the material of a geometric edge (GLN), spring ele-
ments will be generated along that edge based on the width and the distance of the support
nodes.
• Unspecific:
E G
Cs = Ct = (3.99)
H H
• Planar layer with horizontal constraints e.g. for modeling elastic support by columns and
supporting walls (plane stress condition):
E 1 G
Cs = · Ct = (3.100)
H (1 + μ)(1 − μ) H
• Planar layer with horizontal constraints for settlements of soil strata (plane strain condition):
E (1 − μ) 2·G
Cs = · Ct = (3.101)
H (1 + μ)(1 − 2μ) H
• Circular hole with radius R in infinite disk with plane strain conditions (bedded pipes or
piles):
E 1
Cs = · Ct = Cs (3.103)
R (1 + μ)(1 − 2μ)
Including a dilatancy factor describing the volume change induced by shear deformations, we
have for the bedding stresses the following equations depending on the normal and transverse
displacements:
Cracking: Upon reaching the failure stress, the interface fails in both the axial
and the lateral direction. The failure load is always a tensile stress.
If the bedding reaction is applied to a QUAD element, a deformation
in the direction of the local z-axis will create compressive (negative)
stresses.
Yield load: Upon reaching the yield stress, the principal deformation component
of the interface increases without an increase of the stress.
Friction/cohesion: Defining a friction and/or a cohesion coefficient, the lateral shear
stress can not become larger than:
Friction coefficient * normal stress + Cohesion
Please note, that before reaching this limit the stiff-ness CT will pro-
duce the shear stress only if a deformation is present.
If the principal interface has failed (CRAC), then the lateral bedding acts only if 0.0 has been
entered for both friction- coefficient and cohesion.
The non-linear effects can only be taken into account by a non-linear analysis. The friction is
an effect of the lateral bedding, while all other effects act upon the principal direction.
Hint
For the Rigid Frictional Interface Element (I2D), an explicit provision of stiffness parame-
ters is normally not required; instead, the program establishes appropriate stiffness val-
ues automatically in a mesh adaptive manner. Explicit provision of stiffness parameters
overrides this mechanism, the provided values are adopted unmodified.
Interface strength is determined from the reference material by means of the strength
reduction factor:
For the rigid frictional interface, the strength reduction factor defaults to RFAC = 1.0.
(Hint continued...)
Dilatancy coefficient and equivalent mass distribution have no effect on the interface be-
havior.
HMAT
NO Material number − !
TYPE Type of material law LT FOUR
DARC Seepage (Darcy, linear)
FORC Seepage (Forchheimer, nonlin.)
MISS Seepage (Missbach, nonlin.)
FLOW General diffusion (e.g. vapor)
FOUR Heatflow (Fourier, linear)
Hydration of concrete:
JONA Jonasson model
HSCM ”Shrinkage-Core” model
WESC Danish model acc. to Wesche
HMAT allows to specify three individual sets of constants for three tasks diffusion, seepage
and heat flow for a material. Material properties have to be selected according to literature or
experimental data. But some rough estimates (without warranty) are given here.
For any material there might be up to 15 different sets of materials for different temperatures
[ ◦ Celsius] or pore water pressures [ kPa = kN/m2 ] . With a nonlinear analysis the material
values will then be interpolated between those values.
Large jumps (e.g. by more than 10 times) of the material values K and S in adjacent regions
can lead to numerical problems in a non-linear calculation, i.e. to a poor iteration process.
These problems can be solved by material transition regions, smaller and more iteration steps
and/or a coarser iteration tolerance.
Positive pressures represent saturated flow regions while negative values describe unsaturated
soils. Conductivities and Capacities will be interpolated. Free surface problems also use a
variation of the porosity to account for the effective capacity of the free surface.
Material values may be defined isotropic or anisotropic depending on a stress state from the
database. This is performed by additional data given with the same material number and the
types FVOL, FSIN or FSIT. The difference between FSIN and FSIT is given by the fact that
FSIN modifies the values across the crack, while FSIT does this perpendicular to the crack. It
is recommended to use a linear stress field for that purpose and not a plasticity field.
Within HYDRA the user has the possibility to define material properties for the elements either
via a given material number (engineering constructions) or via the nodes (element material
number = 0), especially for ground water models. Both methods may be used together within
the same system.
For a general diffusion according to potential theory, the unit for conductivity (KXX, KYY, KZZ)
is [ sec/m] , for storage coefficient (S) [ kg/m3/Pa] =[ sec2/m] .
e.g.:
If the materials have been defined with CONC, STEE and TIMB, the constants will be preset or
calculated as default for FOUR according to EN 12524 or DIN EN ISO 10456 and the design
codes EN 199x-1-2. For concrete and timber a definition of NSP will select the moisture ratio
with default of 3% for concrete and 10% for timber. For timber, the storage coefficient S is
determined taking into account the ratio of density to dry density (EN 1995-1-2, Table B2).
With explicit definition of a TEMP-value it is possible to create and modify tabulated values.
The default is to use the formulas defined in the Eurocode. The thermal conductivity of the
4500
1.800
4000
1.600
3500
1.400
3000
1.200
1.000 2500 S [kJ/m3/K]
0.800 2000
K [W/m/k]
0.600 1500
0.400 1000
0.200 500
0.00 [◦ C] 0.00 [◦ C]
200.00
300.00
400.00
500.00
600.00
700.00
800.00
900.00
200.00
300.00
400.00
500.00
600.00
700.00
800.00
900.00
100.00
100.00
1000.00
1100.00
1200.00
1000.00
1100.00
1200.00
0.00
0.00
concrete is given by an upper (A=1.0) and a lower (A=0.0) limit. The special effect for S is the
evaporation of pore water.
.0 50.00 18000
18000
.0 45.00
16000 16000
.0 40.00
14000 14000
.0 35.00
12000 12000
.0 30.00 K [W/m/k]
k [W/m/K] 10000 10000
.0 25.00
8000 8000
.0 20.00
6000 6000 S [kJ/m3/K]
.0 15.00 S [kJ
.0 10.00 4000 4000
0.00
100.00
300.00
400.00
500.00
700.00
1000.00
1200.00
1100.00
0.00
100.00
200.00
300.00
500.00
600.00
700.00
800.00
900.00
1000.00
1200.00
1100.00
100.00
200.00
300.00
400.00
500.00
600.00
700.00
800.00
900.00
100.00
200.00
300.00
400.00
500.00
600.00
700.00
800.00
900.00
1000.00
1100.00
1200.00
1000.00
1100.00
1200.00
0.00
0.00
Figure 3.17: Conductivity and Capacity of Structural Steel
1.50 7.00
1.40 6.50
1.30 6.00
1.20 K [W/m/k] 5.50
1.10 5.00
1.00
4.50
0.90
4.00
0.80
0.70 3.50
0.60 3.00
0.50 2.50
0.40 2.00
0.30 1.50
0.20 1.00
0.10 0.50 S [kJ/m3/K]
[◦ C] [◦ C]
100.00
200.00
300.00
400.00
500.00
600.00
700.00
800.00
900.00
100.00
200.00
300.00
400.00
500.00
600.00
700.00
800.00
900.00
1000.00
1100.00
1200.00
1000.00
1100.00
1200.00
0.00
0.00
TYPE JONA
Function of Jonasson, an extension to the Byfors definition:
τ
α = ep b n 1 + (3.107)
τk
Examples of those constants a,b and τk can be found in Appendix A of Heft 512 of the Ger-
man DAfStB, printed with the input record HMAT but the general rule is that you need tests!
Unfortunately there are publications with exchanged parameters a and b.
The values have to be determined from experiments, values for the total heat in the literature
are often defined in [ kJ/kg] . However the following values might give a rough idea:
Q-max [ kJ/m3] a b τ k [ h]
B25 (Heft 512 DAfStB) 167700 -0.925 -0.42 41.82
B35WU (Heft 512 DAfStB) 135750 -1.185 -1.00 17.53
B35 (Heft 512 DAfStB) 105940 -1.605 -1.00 37.46
TYPE HSCM
”Shrinkage-Core” model:
· (τ − τk )
α = τ > τk (3.108)
1 + · (τ − τk )
τk is a period with very low chemical reaction, named d in the original formula. Values for a
mass concrete are given by Dussinger:
TYPE WESC
Danish model according to Wesche:
τk b
α = ep − (3.110)
τ
ƒcc (t) 0
β = = · ep c · t −0.55 c = c1 · (3.111)
ƒcc (28d) z
The water cement ratio w0 /z has to be incorporated in the input value c. Please note that
parameter a is used in the original with two different meanings. The following parameters for
the Wesche model can be found in the literature.
Q-max [ kJ/kg] a c1
w0 /z 0.35 0.55 0.75
Z 55 / Z 45 F 380 1.10 1.15 1.20 -1.50
Z 45 L / Z 35 F 295 1.20 1.30 1.40 -2.80
Z 35 L 216 1.30 1.50 1.70 -4.40
Z 25 1.50 1.90 2.30 -7.10
τk = 24.87[h]
b = 0.84[−]
0 / z = 0.68[−] (3.112)
= 1.63[−]
c1 = −4.4[−] ⇒ c = −4.4 · 0.68 = −2.99
Hint: If the maturity function according to Saul is used, the exponent s there might be specified
with item S.
CONC
N/ mm2
FC Compressiv strength of concrete 1092 *
N/ mm2
FCT Tensile strength of concrete 1092 *
N/ mm2
FCTK Lower fractile strength value fct =fctk,0.05 1092 *
N/ mm2
EC Elastic modulus 1090 *
MUE Poisson’s ratio or shear modulus − 0.2
kN/ m3
GAM Unit weight 1091 25
ALFA Thermal expansion coefficient [1/ K] 1083 1E-5
SCM Typical material safety factor − *
TYPR Type of service state line LT *
LINE constant elastic modulus
A,B short time lines (EN 1992)
N/ mm2
FCR Strength for non-linear analysis 1092 *
N/ mm2
ECR Elastic modulus for servicability 1092 *
N/ mm2
FBD Design bond strength 1093 *
N/ mm2
FFAT Fatigue strength fcd,ƒ t 1092 *
N/ mm2
FCTD Design tensile strength 1092 *
N/ mm2
FEQR Equiv. tensile strength after cracking 1092 0.0
N/ mm2
FEQT Ultimate tensile strength 1092 0.0
N/ mm2
GMOD Shear modulus 1090 *
N/ mm2 1090
KMOD Bulk modulus *
kg/ m3 1189
RHO Density *
GF Fracture energy (MC 2010, [Link]) N/ mm *
MUEC Friction in cracks − *
Only if a value is given for FEQR or FEQT the tensile strength of the concrete is applied for the
design and nonlinear analysis.
C = regular concrete
LC = light-weight concrete
As the German DIN Fachbericht anticipates the design rules of the upcoming Eurocode, but
references cite explicitly the old DIN 1045-1, it is recommended to specify the correct material
type with C_EN or CDIN unambiguously.
The cylindrical strength is always to be input for FCN. The default value is 20.
Some properties are dependent on so called boxed values or other national regulations. The
definition of NORM is used to select those variants. The possible values for example are given
in table 3.1 and 3.2 of EN, resp. table 9 and 10 of DIN 1045-1 resp. table 4 of OEN B 4700
and not repeated here in detail.
FCN = ƒck 12 16 20 25 30 35 40 45 50
ƒck,cbe 15 20 25 30 37 45 50 55 60
FCN 55 60 70 80 90 100
ƒck,cbe 67 75 85 95 105 115
The default values for strength and elastic modulus are derived as follows:
FC = 0.85 · ƒck
2/ 3
FCT = 0.3 · ƒck (ƒck < 55)
= 2.12n((ƒck + (ƒck > 55)
8)/ 10 + 1)
EC = 9500(ƒck + 8)0.3 (EN 1992-1992)
EC = 22000(ƒcm / 10)0.3 (EN 1992-2004)
FBD = 2.25 · ƒct,0.05 / γ (Tbl. 5.3.)
The coefficients αcc and αc for the long term strength effects are defined in the national an-
nexes. The Eurocode suggests values between 0.8 and 1.0 and recommends as default the
value of 1.0. However SOFiSTiK uses a default on the safe side of 0.85, if not specified ex-
plicitly in the INI-file according to the national annex. An explicit definition of FC = αcc ·FCN
is possible of cause. However, for the fatigue, bond or tensile strength values all corrective
factors and the safety should be included in the input data.
For the elasticity modules we have to distinguish between a secant modulus Ecm (Input value
EC of this record) for elastic deformations, especially constraining forces and a tangential mod-
ulus Ec0,28 = 1.05 Ecm or Ec0,28 = Ecm / α , used for creep and nonlinear analysis (input item
The value ECR may be also defined with an explicit unit [-] as factor to the default value, as
recommended in Heft 600 of DAfStB with 0.85 for non lnear analysis.
For light-weight concrete (LC) according to density RHO has to be defined explicitly in kg/ m3 .
GAM and EC will then be defined appropriately.
For light-weight concrete, the tensile strength and bond values and limit strains will be scaled by
a factorη1 . For the ultimate limit stress strain law the bilinear version is selected. The different
coefficients for natural sand and other components are selected by AQUA automatically based
on the defined weight and strength.
The fatigue strength may be specified with item FFAT. The formula given in the EN 1992-1-1
Gl. 6.76 is dependent on a some parameters not fully known here: The number of load cycles,
the kind of the cement, the age at the time of the design and its own safety factor. Thus the
user has to specify the final value directly:
ƒck
ƒcd,ƒ t = k1 · βcc (t0 ) · ƒcd · 1− ; k1 = 0.85 / 1.0
250 (3.115)
p
βcc (t0 ) = es(1− 28/ t0 ) ; s = 0.20(R) / 0.25(N) / 0.38(S)
For detailed analysis of concrete according to appendix 1 you need to know the kind of cement.
You may specify this by appending a literal to the concrete grade
In case of a fire design it is required to distinguish between quartzitic and calcareous aggre-
gates. For the second case an additonal character C may be appended to the cement type:
NC, SC and RC.
The usual stress-strain curve of the C types is the parabolic-rectangular stress-strain diagram.
For non-linear analysis or deformation analysis, there are other types A/B available, following
the expression:
σ k · n − n2
= (3.116)
ƒc 1 + (k − 2) · n
with
n = ϵ/ ϵc1
(3.117)
k = Ec0 · ϵc1 / ƒc
For fc the value fck +8 is used for the curves A and B . The maximum strain is limited according
to the strength. The A line does not possess a descending branch, and it is thus possibly more
stable numerically. The C line has its stress values even for very large strains and will be the
most robust case.
The safety factors SCM are preset to 1.5 (in Italy to 1.6). However, they should be selected at
the design explicitly, because they are dependent on the loading combinations. For concrete
with high strength the factor will be increased by γ’, which will be incorporated in the strain-
stress laws immediately, to allow a global safety factor to be used for the design.
For non-linear analysis with a unified safety factor according to DIN 1045-1 the strength of the
concrete will be reduced to a value of 0.85α fck , while those of the reinforcements will be raised.
These non linear analysis stress-strain laws are generated automatically for DIN 1045-1. For
DIN EN1992 these alternate curves have to be selected with SSLA CALC 1.3 for CONC and
STEE explicitly.
For steel fibre concrete according to DAfStb-Richtlinie (Juni 2021) it is allowed to use higher
concrete tensile bending strength values for elements with concretes up to C50/60. If desired
the user has to enter that value of FCTD explicitly. Values are defined as follows:
ƒ ƒ
FCTD = ƒ,ctd = ƒ,ctk,ƒ · αct / γct
ƒ ƒ ƒ ƒ
FEQR = ƒeq,ctd, = ƒeq,ctk, · αct · κF · κG / γct
ƒ ƒ ƒ ƒ
FEQT = ƒeq,ctd, = ƒeq,ctk, · αct · κF · κG / γct
ƒ
αc = 0.85 (C) / 0.75 (LC)
ƒ
κF = factor for fibre direction
ƒ
κG = factor for member size
ƒ
γct = 1.25
With these values the stress strain laws according to pictures R.1 or R.2. of the DAfStb paper
are created.
The increased saftey factor according [Link] EN 1992 (2004) for cast in place piles with kf =
1.1 has to be specified by the user explicitly.
The default FCN is 25 for B and LB, and 45 for SB. FCT is defined by:
Defaults in accordance with old DIN 1045 / DIN 4227 / DIN 18806:
FCN 10 15 25 35 45 55
FC:
B (DIN 1045) 7 10.5 17.5 23 30
B (DIN 4227) - - 15.0 21 33
FBD: -
B (DIN 1045) 1.4 1.8 2.2 3.0
EC 22000 26000 30000 34000 37000 39000
The density has to be specified in case of light-weight concrete. The default for GAM and EC
then complies with DIN 1055. A bilinear stress-strain curve is usually employed for light-weight
concrete.
For detailed analysis of creep and shrinkage according to DIN 4227 you need the kind of ce-
ment and the consistency. You may specify this by appending a Literal to the class of concrete
DIN 4227 has some contradictions about the bond stress. Chapter 13 gives values which
correspond quite well to the ratios given in table 7 of appendix A1. But these values do not
match those given at DIN 1045 Table 19. Thus we have decided to enlarge the FBD values
for concrete SB by a factor of 1.43. With that amendment the value may be used for the bond
design according to chapt. 13 of DIN 4227.
For standard concrete a parabola-rectangular stress-strain diagram will be selected. SCM will
default to 1.00. If you analyse composite sections you might want to change the value. High
strength concrete will have lesser ultimate strains.
The design strength is calculated according section [Link] and [Link] with two factors ηƒ c
and ηt . The first factor is applied for the definition of FCN automatically. As the second factor
is depending on the design situation it can be modified similar to the Eurocode with a boxed
value ALF-CC beforehand or later.
The elastic moduli are calculated based on the mean strength. For light-weight concrete the
values are corrected depending on the specific weight. The default stress-strain diagram is
always according to the deformation stress strain law of EC-2, even for the design. SCM will
be preset with 1.5.
FCN = ƒck 12 16 20 25 30 35 40 45 50
ƒck,cbe C 15 20 25 30 37 45 50 55 60
ƒck,cbe LC 13 18 22 28 33 38 44 50 55
This code allows a characteristic strength of the concrete depending on the age of the concrete
and a calculation strength depending on the duration of the loading t , to be defined by the user
explicitly. Further we have:
FC = 0.85/Θ · ƒck
Θ = 1.00 (t > 24 h)
Θ = 0.90 (24h > t > 1 h)
Θ = 0.85 (1h > t )
FCT = 0.6 + 0.06 · ƒck (ƒck < 60)
= 0.275 · ƒck 2/ 3 (ƒck ≥ 60)
1/ 3
EC = 11000 · ƒck
FBD = 0.60 · Ψs 2 · ƒct,0.05 / γ ; Ψs = 1.5
There is also a class DUCT for the UHPRFC Ductal FM. For the stress strain law the values
GF is taken as w0.3/lc.
C = 16, 20, 25, 28, 30, 35, 40, 45, 50, 54, 55, 58, 60
LC = 8, 12, 16, 20, 25, 30, 35, 40, 45, 50, 55, 60
As the BBK does not state any details about the transition zone in the work law, the user has
the possibility to influence with FCR between 0.6FCK and 1.0FCK the shape of this curve.
For E-modulus and tensile strength specific formulas are provided in the design codes. For the
design it is allowed to use a rectangular stress block, but the default is the Parabula-Rectangle-
Diagram.
CAN = 2008 / 2018 regular concrete with cube strength fk (γ=1.50)
CAL = 2008 / 2018 light weight concrete with cylindrical strength fck (γ=1.50)
CAN = 1996 regular concrete with cube strength Rck (γ=1.60)
CAP = 1996 prestressed concrete with cube strength Rck (γ=1.50)
The default values for design strength and elastic modulus are as follows:
The nominal strength FCN is the cube strength. The design strength is obtained by
British Standard employ a parabolic rectangle curve, starting from a design cube strength β =
p p
FC/0.67 with 0.24 β strain at full plasticity and an initial stiffness of 5.5 β according to Figure
2.1. The safety factor SCM is preset to 1.5. The bond strength will be set to the non-physical
p
maximum value for table 3.28 of BS 8110 of ƒc .
For Hong Kong slight modifications to British Standards are selected with the country code
p
852. The initial stiffness will then be 5.0 β.
The test values of the cylindrical strength have to exceed the class value by a certain amount
based on the standard deviation. Chapter 5.3 specifies default values for this required distance
as 7.0, 8.5 and 10.0 MPa for class values of fc ’ up to 21, until 35 and above.
Æ
As the value of ƒc ′ should not exceed the value of 25/3 MPa in general and different reduc-
tions
Æ have to be applied for lightweight concrete, we use the tensile stress to define the value of
ƒc ′ . The modulus of rupture fr is the upper fractile value of the tension strength. ACI [Link]
defines:
q
ƒr = 0.75 · ƒc ′ < 0.75 · 25/ 3 (3.120)
The bond strength will be set to the non-physical maximum value for chapter 12.2 of ACI 318
p
to a value of ƒc .
The modulus of elasticity is given by Ec = 5600 · fck 1/ 2 and Ecs = 0.85Ec . The design curve is
a Parabola-Rectangle with 0.85 · fcd, as tensile strength we have fct,m = 0.3 · ƒck 2/ 3 ; fct,knƒ =
0.7 · fct,m ; fct,sp = 1.3 · fct,m . Bond strength fbd = η1 · η2 · η3 · fctk,nƒ /γc ; η1 =2.25.
The nominal strength FCN (15 to 80) and the the design strength are taken from table 4.1.3./4.
Youngs modulus is derived from 4.1.5.
The compressive and tensile strength is selected according to tables 6.7/6.8. The elasticity
modulus is taken from table 6.11.
TYPE=B,C,LC,SIA TYPE=LB
ƒc ƒc
ε [%◦ ] ε [%◦ ]
−2.0 −3.5 −2.0 −3.5
ƒc ƒc σ
E=
ε
ε [%◦ ] ε [%◦ ]
ε −3.5 −2.0 −3.5
TYPR=B,C,R,RS TYPR=A
σ[N/ mm2 ] σ[N/ mm2 ]
ƒc ƒc
ε [%◦ ] ε [%◦ ]
εc1 εc −2.2 εc
STEE
N/ mm2
FY Yield strength (f0.01 or f0.02 ) 1092 *
N/ mm2
FT Tensile strength 1092 *
N/ mm2
FP Elastic limit 1092 *
N/ mm2
ES Elastic modulus 1090 *
MUE Poisson’s ratio − 0.3
kN/ m3
GAM Unit weight 1091 *
ALFA Thermal expansion coefficient [1/ K] 1083 *
scm Typical material safety factor obsolete − *
use SCM0,SCM1,SCM2
There are some, but not very rigorous checks about the usage. While structural steel types
(e.g. S 235, ST, AL etc) can be used only for cross sections, prestressing steel is only allowed
for reinforcements, cables and tendons.
For structural steel that should get the material values of the product standard EN 10025-2 till
-6 the plate thickness TMAX has to be input. E.g. the input
defines in the first case the limit values of the stiffness for the structural steel of the quality
355 according to EN 10025-2 for the plate thickness of 63 mm < t ≤80 mm and in the second
case the limit values of the stiffness for the structural steel of the quality 355M according to EN
10025-3 for the plate thickness of 80 mm < t ≤100 mm.
Note: In the National Annexes the material values may be specified differently.
The specifications which are a part of the steel quality and are printed in italics (e.g. T, T8,
4) describe the thicknesses. They have to be input by the user only in the case of a CADINP
input.
The strength values especially for the high strength steels may vary depending on the manu-
facturer and the alloy considerably, a check of the assumed values against the actual values is
strongly recommended.
The maximum allowed plate thickness may be specified via TMAX, reducing the strength val-
ues according to the design codes. As different strength values within the section may lead to
consistency problems for some design tasks, we allow only for a unified approach based on
the maximum thickness. It will be checked during the generation of sections. For many grades
one may append an identifier for a maximum plate thickness for a direct definition.
Attention: The maximum thickness will be also used to control if the design of a composite
section with DESI in AQB will allow compressive strains beyond the yield limit. For sections of
classes 3 and 4, this limit will be observed by default. For sections of class 1 or 2, the strain
becomes unlimited by definition of TMAX 0.0.
The safety factor SCM0 is the value for the resistance of sections (γM0 = 1.0 according EN
1993-1-1). It will be preset according to the material and the design code and will be used to
calculate the design forces for steel and composite sections. For the design by member checks
the input ’CTRL FORC VAL3 GAM SCM’ can be used which is switching to the old safety factor
scm (γM1 according EN 1993-1-1).
For the Russian design steel the defaults are nominal values according to SB 16.13330.2011
any conversions, round off or higher values for rolled steel have to be specified by the user
explicitly.
FY FT EPST BC np TMAX
AWP 3004 180. 220 10 B 23 6.0
AWP 3005 150. 170 10 B 38 6.0
AWP 3103 120. 140 20 B 31 25.0
AWP 5005 95. 125 20 B 25 12.5
AWP 5052 160. 210 40 B 17 40.0
AWP 5049 190. 240 30 B 20 25.0
AWP 5454 220. 270 20 B 22 25.0
AWP 5754 190. 240 30 B 20 25.0
AWP 5083 250. 305 30 B 22 40.0
AWP 6061T4 110. 250 120 A 23 12.5
AWP 6061T6 240. 290 60 B 23 12.5 / 80.0
AWP 6082T4 110. 205 120 A 25 12.5
AWP 6082T6 240. 295 60 B 25 6.0 / 12.5 /100.0
AWP 7020T6 280. 350 90 A 19 12.5
AWP 8011 110. 125 20 B 37 12.5
FY FT EPST BC np TMAX
AW 6063T5 110. 160 70 B 13 3.0 / 25.0
AW 6063T6 160. 195 80 A 24 25.0
AW 6063TX 180. 225 80 A 21 10.0 / 25.0
AW 6063DT 190. 220 100 A 31 20.0
AW 6005T6 200. 250 80 A 20 5.0 / 10.0 / 25.0
AW 6106T6 200. 250 80 A 20 10.0
AW 6082T4 110. 205 140 B 8 25.0
AW 6082T5 230. 270 80 B 28 25.0
AW 6082T6 260. 310 100 A 25 5.0 / 15.0
AW 6082ER 240. 295 100 A 17 15.0 / 40.0
AW 6082DT 240. 310 100 A 17 5.0 / 20.0
AW 7020T6 275. 350 100 A 19 15.0 / 40.0
AW 7020DT 280. 350 100 A 18 20.0
FY FT EPST BC np TMAX
Eurocode
AC 42100 147. 203 20.
AC 42200 168. 224 15.
AC 43000 63. 126 12.5
AC 43300 147. 203 20.
AC 44200 56. 119 30.
AC 51300 70. 126 20.
EC 9 and the new DIN 4113 (2002) use the American system for classification of aluminium
alloys. As there are more than 300 different materials available, with significant differences of
properties, the user should check the thickness limit and strength parameters thoroughly. For
• For plates we use the type AWP to distinguish them properly from the profile and tubes.
For untreated alloys only the number of the alloy is necessary.
• For tubes and profiles the type AW is used and either ER or DT (drawn tubes) or the
important criteria of the heat treatment as T4, T5, T6 or TX for T66 has to be appended to
the alloy as characters.
• For the castings the case of a cocille and temper F or T6 has been selected as default.
As the reduction of the strength in the HAZ is depending on the welding process and the
thickness, the user has to define a separte Materialnumber for those regions with a explicitly
reduced values for FY and FT.
DIN 4113 requires the stress for aluminum to be reduced for creep effects according to chapter
6.3 with a factor c between 0.8 and 1.0. As the exact evaluation would be rather complex, this
is accounted for by the global factor of safety 1/c. The default on the safe side is a value of
1.25. Better values have to be specified explicitly by the user.
A reduction for generally higher temperatures has to be specified explicitly, but the reduction
for the fire case is available via the stress strain law.
The values ES and GAM are for all classes with 70000 [ MPa] and 28.0 [ kN/m3] preset, the
values FP and EPSY are selected according to the data of table 10 of the DIN 4113.
FY FT FP FYC ES GAM γs
SP
52-101-2003
SNIA 240 235 400 190 235 - 200000 78.5 1.10
SNIA 300 300 500 235 300 - 200000 78.5 1.10
SNIA 400 400 500 320 400 - 200000 78.5 1.13
SNIA 500 500 600 345 460 - 200000 78.5 1.15
SNIB 500 410 500 290 360 - 200000 78.5 1.20
FY FT FP FYC ES GAM γs
RK 5.03-
33-2005
SNIA 600 600 800 360 540 - 200000 78.5 1.15
SNIA 800 800 1000 360 575 - 200000 78.5 1.15
SNIA 1000 1000 1000 360 600 - 200000 78.5 1.20
FY FT FP FYC ES GAM γs
SNIP
2.03.01:
SNIA 240 210 400 210 210 - 210000 78.5 1.05
SNIA 400 350 500 280 350 - 210000 78.5 1.05
SNIA 500 435 600 300 435 - 200000 78.5 1.10
SNIA 600 520 800 300 470 - 190000 78.5 1.15
SNIA 800 695 980 300 500 - 190000 78.5 1.20
FY FT FP FYC ES GAM γs
SNIA 1000 870 1000 300 500 - 190000 78.5 1.20
SNIB 500 415 1000 290 390 - 170000 78.5 1.20
SNIB 1200 1050 1200 785 500 - 200000 78.5 1.20
SNIB 1300 1130 1300 835 500 - 200000 78.5 1.20
SNIB 1400 1210 1400 940 500 - 200000 78.5 1.20
SNIB 1500 1300 1500 1000 500 - 200000 78.5 1.20
SNIB 1600 1390 1600 1000 500 - 200000 78.5 1.20
SNIB 1700 1475 1700 1000 500 - 200000 78.5 1.20
For the steel type B it is possible to attach to the grade extra characters defining the ductility:
For the German steel type BST the characters SA/SB resp. MA/MB for bars resp. meshes
allow also to select between new DIN 1045-1 and old DIN 1045.
The safety factor is preset for most reinforcing and prestressing steels to 1.15 and 1.05 (BS)
respectively. The safety factor becomes effective immediately for the calculation of the full
plastic internal forces of steel and composite sections.
For non-linear analysis with a constant safety factor according to DIN 1045-1 the strength of the
concrete will be reduced, while those of the steel will be raised. For this a special serviceability
stress-strain law is generated with a safety factor of 1.3.
The Russian SNIP has a reduced strength for shear links and inclined bars. This value is taken
from the value FP. In the very general case, it might be necessary to use a separate material
with a reduced strength.
Attention:
Some material parameters may depend on other parameters not known to the program. E.g.
the dynamic stress range of the reinforcements is not only depending on the diameter TMAX
of the bars but also on the curvature and the type of the material of the duct (steel, plastic). In
all those cases it is necessary to use different material numbers and to specify the deviating
values explicitly.
Table 3.1 of EN 1993-1-11 provides different elastic moduli for different prestressing steels
depending on the selected cable sections (group B tension components).
The material safety is preset with 1.0. Further information is available in the description of the
record CABL.
3.15.5 Relaxation
Relaxation of tendons is implemented in AQUA/AQB as a product. While the time factor is
specified in AQB, AQUA defines the stress dependant factor for the reference time of 1000 h.
This may be accomplished either by a linear relation established by two values at 0.55·fpk and
0.70·fpk or via selected literals for item REL2 as quadratic function according to CEB / EN1992
or the general function according to BPEL annexe 2 or AS 3600.
In Germany the values of the general technical approval. The values of a typical approval are
provided as default, but may be redefined with MEXT if required.
If different bond properties should be applied, different materials have to be specified. The
bond value is needed for
• Reinforcements
• Crack width
• Fatigue
• Limitation of stress increase for tendons
The relative bond properties are mainly intended in the design codes to account for the reduced
values of tendons. But it should also account for different diameters and properties of passive
reinforcements. For the crack width a mean value has to be used where the diameter and the
The second coefficient is used for those design codes using special values not clearly de-
ductible from the relative bond coefficients alone. This is especially the base version of the
Eurocode, but also the russian SNIP:
The defaults are given in principal in the following table. However, there are some deviations
depending on the design code and the strength possible:
R K1
Reinforcing steel 1.0 0.8
Prestressing steel 0.5 1.6
Depending on the steel type and grade the values EPSY and EPST as well as FP will be
preset. With explicit definitions you may suppress:
In general the stress-strain laws are identical for serviceability and ultimate limit design. How-
ever, for reinforcing steel according to EC2 and DIN 1045-1 there are numerous explicit
changes.
As the safety factor concept will not generate an affine curve for the ultimate state, the safety
factor will be applied immediately.
Although the tensile strength for reinforcement steel with standard ductility will be reached at
25 o/oo, it is not allowed to use this in the design according to DIN 1045-1. The stress strain
laws for design and non linear analysis differ therefore for those materials.
ƒy ƒt
ƒy
ϵ(%) ϵ(%)
2.0(−EPSY)
S (Eurocode 2) PS (Eurocode 2)
ƒt ƒt
0.8ƒt
ƒy / 0.9ƒ t
ϵ(%) ϵ(%)
(EPST) 5.0(−EPSY)
TIMB
N/ mm2
EP Elastic modulus parallel to fibre 1090 *
N/ mm2
G Shear modulus 1090 *
N/ mm2
E90 Elastic modulus normal to fibre 1090 *
QH Poisson’s ratio yz (polywood panels) − *
QH90 Poisson’s ratio xy / xz (solid wood) − *
kN/ m3
GAM Unit weight 1091 *
ALFA Temperature elongation coefficient [1/ K] 1083 0.0
SCM Material safety factor − 1.3/*
N/ mm2
FM Bending strength 1092 *
N/ mm2
FT0 Tensile strength parallel to the fibre 1092 *
N/ mm2
FT90 Tensile strength normal to the fibre 1092 *
N/ mm2
FC0 Compressive strength parallel to fibres 1092 *
N/ mm2
FC90 Compressive strength normal to fibres 1092 *
N/ mm2
FV Shear strength at center (shear force) 1092 *
N/ mm2
FVR Shear strength at the edge (torsion) 1092 *
N/ mm2
FVB Shear strength for plate bending 1092 *
N/ mm2
FM90 Bending strength normal to fibres 1092 *
N/ mm2
G90 Shear modulus for plate bending 1090 *
OAL Meridian angle of anisotropy degree 0.0
OAF Descent angle of anisotropy degree 0.0
TIMB allows the definition of all materials with a preferred fibre orientation. As type you have the
timber classes of Eurocode, respective DIN, OENORM and SIA and the German compound
fibre types. As the EN 1995 does not specify any strength values, those values have to be
specified for all other countries explicitly.
With the EN 1995 (EC5) and the derived design codes correction factors kmod have been
introduced for the permissible stresses and kdeƒ for the deformations. The distinct value is
depending not only on the material but also on the service class and the duration of the loading.
The following table shows the values for solid timber, however there are much more values
available within the program and the user may change theses values explicitly.
Service Class
Class of load duration
1 2 3
permanent 0.60 0.60 0.50
long term 0.70 0.70 0.55
medium term 0.80 0.80 0.65
short term 0.90 0.90 0.70
very short term 1.10 1.10 0.90
k-def 0.60 0.80 2.00
The service class may be either specified in general with the definition of the design code
NORM or be appended to the class definition with a colon. The input of TIMB C 30:2 selects a
solid soft wood of strength class 30 and the service class 2.
The base values of the strength are given in EN 338 resp. EN 1194. But the transverse shear
strength FV has to be defined reduced as kcr · ƒ accounting for cracks. In Germany this
coefficient has been selected to provide the same shear strength for all strength classes.
There are many composite materials in timber constructions. Although a precise treating is
possible with composite sections or FLAY, the design codes provide equivalent materials for
that purpose. As the strength is no strongly dependant on the thickness of the construction
part, the definition of that value with TMAX is mandatory.
The description of a transverse orthotropy material law has one direction that has different
properties (fibre direction), while the description in the plane perpendicular to this remains
isotropic. The law defined with TIMB is formal equivalent, however not identical with that de-
fined via MATE. If x is chosen as this special direction it holds:
σ (σy + σz )
ε = − μ90 · (3.122)
E E90
σy σz σ
εy = − μ· − μ⋆90 · (3.123)
E90 E90 E
σz σy σ
εz = − μ· − μ⋆90 · (3.124)
E90 E90 E
E
μ⋆90 = μ90 · (3.125)
E90
It should be noticed, that the poisson’s ratios μ90 and μ90 ∗ are no longer bound to 0.5 and are
strongly connected to the ratio of the elasticity moduli, as the resulting stress-strain matrix has
to be symmetric.
The order of the indices of stress and strain components notation is defined as follows:
For the axial symmetric case x denotes the axis of rotation while y represents the radial and z
the tangential direction.
Furthermore holds:
1−n·μ
2
2 μ1 +n·μ2 μ
E1 · (1+μ )·m 2
E1 · (1+μ )·m E1 · m2 0 0 0
1 1
μ1 +n·μ22 1−n·μ22 μ
E1 · (1+μ )·m E1 · (1+μ )·m E1 · m2 0 0 0
1 1
μ2 μ2 1−μ1
E1 · E1 · E2 · 0 0 0
D= m m m (3.127)
0 0 0 G1 0 0
0 0 0 0 G2 0
0 0 0 0 0 G2
E1
n = , m = 1 − μ1 − 2 · n · μ2
2
(3.128)
E2
Plane strain conditions: The x direction is defined as the fibre direction (=normal to the isotropic
1−μ1 μ
E2 · m E1 · m2 0
μ 1−n·μ2 2
E1 · m2 E1 · (1+μ )·m 0
D = 1
(3.129)
0 0 G2
1μ +n·μ2
2 μ
E1 · (1+μ E1 · m2 0
1 )·m
E1
n= , m = 1 − μ1 − 2 · n · μ2
2
(3.130)
E2
Plane stress conditions: The x direction is defined as the fibre direction (=normal to the
isotropic plane). The material stiffness matrix is obtained via inversion of the reduced strain-
stress matrix and reads:
E2 E1 ·μ2
0
1−n·μ22 1−n·μ22
E1
D = E1 ·μ2
, n = E
E1 (3.131)
0
1−n·μ22 1−n·μ22
2
0 0 G2
Axial symmetry:
The general case of anisotropy does not need to be considered since axial symmetry would
be impossible to achieve under such circumstances. A case of interest in practice is that of
the fibre direction parallel to the rotational axis x, i.e. the x direction is normal to the plane of
isotropy. For such a case the material stiffness matrix reads:
1−μ21
n μ2 (1 + μ1 ) 0 μ2 (1 + μ1 )
μ2 (1 + μ1 ) 1 − nμ2 0 μ1 + nμ2
D = A · 2 2 (3.132)
G90
0 0 A 0
μ2 (1 + μ1 ) μ1 + nμ2
2 0 1 − nμ2
2
E2 · n E1
A = , n = , m = 1 − μ1 − 2 · n · μ2
2
(3.133)
(1 + μ1 ) · m E2
For beams the fibre direction is identical with that of the beam axis.
For planar systems (TALPA) the value OAF is the angle between the fibre direction and the
element x-direction. The values E90 and μ90 then hold within the isotropic plane whose normal
is given by the (skew) fibre direction.
For shells and plates it might be possible (eg. plywood) that there are fibres in both x and y
direction. The anisotropy effects thus reduces to different shear moduli for in plane membrane
shear (Gm=0.5E90 /(1+μ)) and the transverse shear force directions (G). This may be accom-
plished either with an explicit definition E90 == E or with a layered material (FLAY). For the case
of vertical boards tied together, you may use the orthotropy factors of thickness description (→
QUIAD) or use a 3D model.
In three-dimensional systems (continuum elements) the default fibre direction is the element z-
direction. Other fibre directions can be specified by defining the three-dimensional orientation
of the isotropic plane (=plane, whose normal direction is the fibre direction) via meridian and
descent angle, known from geology (compare MATE).
MASO
N/ mm2
FK Compressive strength fk 1092 *
N/ mm2
FB Compressive brick strength fb 1092 *
N/ mm2
FM Mortar strength fm 1092 *
N/ mm2
FVK0 Adhesional shear strength fk,0 1092 *
N/ mm2
FV Maximum shear strength 1092 *
N/ mm2
FT Tensile strength 1092 *
N/ mm2
FTB tensile brick strength 1092 *
SCM Material safety factor − *
N/ mm2
E Elastic modulus 1090 *
N/ mm2
G Shear modulus 1090 0.4E
MUE Poisson’s ratio − 0.25
kN/ m3
GAM Unit weight 1091 *
ALFA Temperature elongation coefficient [1/ K] 1083 *
The parameters follow EN 1996-1-1. It is to be defined in the national annex if the properties
are based on experiments (i) or on calculatoric values (ii) specified with table 3.3 and formulas
3.1 to 3.2. In both cases there are groups 1 to 4 depending on the holes in the brick, which
have to be defined appended to STYP with a colon (eg. SB:3)
In all those cases where the design code describes a reduction of the K-value, this may be
defined with a definition like 80[%] for FK.
AQUA supports the German national annexe. Table Na.4 for HLzA/B etc. is selectable with
group 2, table NA.5 for HLzW etc. is selectable with group 3. Table NA.9 is selectable for solid
bricks (Vn, Vbn, Vm, Vmb) with group 1, solid bricks (Vbl S und Vbl SW) with group 2, solid
bricks (V und Vbl) with group 3 and hollow bricks (Hbl S und Hbn) with group 4.
For masonry according to BS 5628-1/2 the Literals "‘BS-1"’ resp. "‘BS-2"’ with the group iden-
tifier A to D appended with a colon has to be given for the STYP definition. Selectable are the
mortar designations I to IV. FT is the tensile strength for bending according Table 3 ”parallel to
bed joints”, FV is the vertical shear strength according pict. 2 and clause (25, part 2), FV0K is
the basic shear value according clause 25 part 1, FBT is the bending tensile strength according
Table 3 ”perpendicular to bed joints”. FT and FBT vary considerably and should therefore be
specified.
SSLA
Stress-Strain curves define the dependency of the stress from the strains. They are required
for the design and nonlinear analysis. There are three different types to be distinguished:
For the materials according to the several design codes CONC, STEE, TIMB and MASO there
For reinforcing and prestressing steel the EN 1992-1-1 specifies two variants for the stress-
strain curves:
An exception to this is the National Annex DIN, which only explicitly prescribes variant a. For
EN 1992-1-1 and all available National Annexes the variant a is the default. A switch to variant
b is made by input 0 for the ultimate strain of the reinforcing or prestressing steel in the record
STEE item EPST.
For a nonlinear analysis it is most reliable and most economical to use the method of double
book keeping (Quast) where the true deformations are evaluated based on the mean values
and the ultimate capacity is evaluated separately with the design values. This is congruent to
the classic approach of a linear analysis and a design as a post processing step. EN 1992-
1 covers this in section 5.8.(3) with the very first sentence or the more detailed hints in the
German NA to this subject. This is controlled in AQB/ASE/STAR2 with an input of:
For stability problems it might be requested to apply a safety factor on the stiffness too. If this
factor is included in the stiffness or at the load is a matter of taste, but the nonlinear nature
of the problem does not allow to use the same factor on either side of the equations. This is
controlled in AQB/ASE/STAR2 with an input of:
DESI ULTI KSV ULD ULD ; NSTR ULTI KSV SLD SLD
While the mean values of the concrete are specified in the design codes, there is no definition
for the reinforcing steel. A common assumption is to increas the strength by a factor of 1.1. As
this is comparable to the safety factor and a yielding will decrease the stiffness considerably it
is also common practice to use the option SLD SL instead.
However, it might be desirable to perform design and nonlinear analysis within a single run.
Then special curves are needed based on design values of the strength, but a stiffness based
on the mean values with a special saftey factor. This is controlled in AQB with an input of:
EN 1992-1 covers this in section 5.8.(3) with the second and third sentence. As the German
NA provides an alternate version in section 5.7 with a global system safety of 1.3, the user has
to decide between those two methods by the definition of the safety factor for the CALC curve.
The default is according to 5.7. It has to be pointed out, that this method does not include any
safety on the stiffness. Thus there are three possibilities to include that:
• Reduction of the ECM according Heft 600 DAfStB with 0.85, to be specified with TVAR
ALF-CE 0.85 SCOP DESI
• Division of the total curves by a factor 1.3, which is not provided here
• Ultimate load analysis with a factor of 1.3 and stress-strain curves without a safety factor
If the default stress-strain curves are not applicable, stress-strain curves must be defined im-
mediately after the input of the material. For special applicationes each set may have multiple
temperature levels, to be defined in ascending order of temperature. A stress-strain curve
starts with one of the possible headers:
The safety factors are predefined as specified with the material, but may be changed if needed.
In particular it is possible to modify a standard stress-strain law with an own safety factor, by
defining a header record only. The design codes use the safety factors quite differently, some
materials will be divided by the safety factor in total, others reduce only the maximum stress
value and keep the elasticity modulus. A positive safety factor will select the first case, while
a negative value will select the second one. For the standard design tasks the material safety
factors are chosen by AQB depending on the loading condition and design code.
For the CALC curves we have a slightly other rule to account for the German NA. A negative
safety factor will activate the curve according to chapter 5.7, while a positive safety factor
will trigger the rule according to chapter 5.8 with a safety factor γce for Ecm . This factor is
predefined in the INI-file, but may be overwritten as boxed value with a TVAR GAM-CE. A
safety factor of 0.0 will deactivate the calculatoric curves.
The stress-strain curve follows. Each consists of several data points in an ordered sequence.
For each data point it is specified whether it should behave as a vertex (linear polygon line) or
it should be part of a smooth curve (quadratic or cubic parabola).
σ σ
POL SPL
ϵ ϵ
The user must make sure that a sufficiently large strain range gets covered and that the zero
point constitutes a data point of its own. Strains outside the defined range will have for TYPE
EXT the last defined stress value and will use the tangent at the last point, provided this has a
positive elasticity module. If the TYPE is specified as LIM, stresses outside the defined range
will become zero, which however lead to trouble in numerical behaviour of non-linear iterations.
On the other hand the stresses outside the defined range are extrapolated by an input of TYPE
EXT. Default is EXT, however for ultimate limit state of concrete it is LIM.
For concrete without explicit data points and without temperature it is possible to define with
EPSS a factor for the strains. This may be used to account for creep effects as specified in
the EN 1992 with 1+ϕ. It should be noted however, that in this case the creep becomes also
effective for variable actions.
For a fiber or a FE-section the general analysis method will account for temperatures according
to the Eurocodes EN 1992 to 1999 automatically. However for special cases it is also possible
to define between 2 and 15 discrete temperature levels TEMP, to be interpolated. EPST defines
the thermal elongation for that temperature, the default will be taken from the Eurocode.
For a zonal method, the section will be subdivided in several zones (polygons) by different
material numbers, where every zone has a constant average temperature to be specified for
each of these materials with TEMP. One may shift the stress strain curve wih a definition of
EPSS SHIF by the thermal strain allowing to activate all eigenstresses directly without the
need to define a temperature at the section itself in detail.
The contribution of the concrete between the cracks (Tension Stiffening) may be taken into
account by a modification of the stress strain law of the concrete or the steel. With a single
record SSLA SERV a modification of the reinforcement stress strain law for that purpose is
possible. However a solution is only possible if the ratio of reinforcement is large enough to
avoid the complete rupture of the reinforcement, which is equivalent to the requirement that
the initial crack stress in the reinforcement must not exceed the yield limit.
The reinforcement ratio MUET is only determined in advance for an annular section (or those
with a similar evenly distributed reinforcement). For all other cases, the design task should
adopt the stress strain law accordingly.
MEXT
NO Number of material − 1
EXP Name of an (exposure) class Lt4 -
TYPE Type of constant LT !
VAL Value of material constant ∗ -
VAL1 First additional material value ∗ -
VAL2 Second additional material value ∗ -
VAL3 Third additional material value ∗ -
VAL4 4th additional material value ∗ -
VAL5 5th additional material value ∗ -
VAL6 6th additional material value ∗ -
VAL7 7th additional material value ∗ -
VAL8 8th additional material value ∗ -
VAL9 9th additional material value ∗ -
VA10 10th additional material value ∗ -
...
VA20 20th additional material value ∗ -
With MEXT you may define special material values for any type of material. The values may
be assigned to different regions (edges) of the material with separate values. This is defined
by a freely selectable literal EXP. Defaults for those values may be specified in the INI-file. It
has to be checked individually to which extend the defined data is really used in the analysis
modules.
You can find further explanations and examples for this workflow in our Online Help -> Tutorials
[Link]
cross-section_air-[Link]
Warning
There are no practical experiences for the procedure according to EN 1992-2 B.104.
The use is recommended only with caution in exceptional cases. The feature can be
activated by adding "H" next to the concrete class, e.g. C 60H.
Explicit coefficients for tendon relaxation according to EN 1992-1 is required as a table. The
default for ENC2 in Germany is a typical table for low relaxation steel from a general technical
approval (abZ). The definition of an explicit table starts with a headder without EXP with up
to 9 descending time values. The following records contain on ascending order the losses in
percent for selected ratios of Ri/Rm. As EXP the literal "‘Rnn"’ has to be specified, where nn
ist the ration of Ri/Rm in percent.
MEXT 12 TYPE EIGE EXP - VAL 1000000 500000 5000 1000 200 10 1
MEXT 12 TYPE EIGE EXP 'R55' VAL 1.2 1.0 0.0 0.0 0.0 0.0 0.0
MEXT 12 TYPE EIGE EXP 'R60' VAL 2.8 2.5 1.2 0.0 0.0 0.0 0.0
MEXT 12 TYPE EIGE EXP 'R65' VAL 5.0 4.5 2.0 1.3 0.0 0.0 0.0
MEXT 12 TYPE EIGE EXP 'R70' VAL 7.0 6.5 3.0 2.0 1.0 0.0 0.0
MEXT 12 TYPE EIGE EXP 'R75' VAL 10.0 9.0 4.5 3.0 2.5 1.2 0.0
MEXT 12 TYPE EIGE EXP 'R80' VAL 14.0 13.0 6.5 5.0 4.0 2.0 1.0
Multiple curves for creep are allowed with different loading start times of the creep curve. For
any creep interval the curve nearest before the current start time will be selected and shifted
to the true start point.
Example for a definition of three explicit creep curves and an explicit shrinkage curve with
name CURV for a concrete no. 1:
Several EXP may be used to define sets of variant curves for different environmental conditions
or nominal effective thickness.
Hint: In table 4 of DIN 1055 part 4 slightly larger values are defined for k. The value VAL1 is
the friction coefficient to be used for forces from flow along a surface, especially along a beam
axis.
With BORE a soil or bore profile is described defining material layers along an axis. The use
of which is different.
For simple regular cases it is also possible to define design stresses for the soil pressure with
record BTAB. These values are than applicable without reference to the strata data.
Hint
PHI, GAM and GAMA are not supported at the moment for the HASE calculation.
BLAY is used for the definition of the soil layers of a corresponding BORE profile. This data is
then used in program HASE for the determination of the stiffness and resistance properties of
soil and piles.
BORE NO #nb X #x Y #y Z #z
BLAY S #s1 ES #Es1 VARI para PMAX #P1
BLAY S #s2 ES #Es2 VARI line DES #dEs2 PMAX #P2
BLAY S #s3 ES #Es3 VARI cons PMAX #P3
BLAY S #s4
s1 Es1 s1 P1
L1 L1
s2 Es2 s2 P2
dEs2
L2 L2
s3 Es3 s3 P3
L3 L3
s4 s4
• First layer L1 starts at the depth #s1 and ends at the next defined depth #s2. L1 is assigned
a parabolic stiffness distribution. Since there is no explicit stiffness increase #dEs1 defined,
a continuous distribution is realized - the concluding stiffness value will be equal to the
stiffness value of the subsequently defined BLAY (#Es2).
• Second layer L2 has a linear stiffness distribution. This time, #dEs2 is defined, so the
concluding stiffness value will be #Es2+#dEs2.
• Third layer L3 has a constant stiffness distribution with the value #Es3.
If there is only one BLAY defined, then the ending depth is by default 999m. If the number
of BLAY record is larger than 1, the last BLAY defines the ending depth S while the other
properties of this BLAY are ignored.
Apart from stiffness modulus ES, all other properties of the BLAY record (MNO, MUE, PMAX,
PMAL, C, PHI) are constant within a layer.
Stiffness modulus ES and Poisson’s ration MUE can alternatively be set by reference of a
material number, where corresponding elastic material properties are defined. If within the
same BLAY record MNO and ES and/or MUE are defined, the values defined within BLAY
have precedence over those defined by material MNO. If the Poisson’s ratio MUE is defined
and larger than 0.0, then ES is interpreted as the elastic modulus, and a Boussinesq method
is used in program HASE (see HASE manual for more details).
By default (no input) the non-linear resistance properties (PMAX, PMAL and C) are switched
off, meaning that the x-pile nodal forces are not limited (elastic analysis). As in the default
case, the input of a resistance property smaller or equal zero will result in the x-pile contact
forces to be unlimited. Any input of the resistance properties larger than zero will activate the
Hint
• As a prerequisite for a meaningful soil-layer interpolation, all defined bore profiles
(BORE) within a grid must have the same number of layers (BLAY).
• For a consistent input of the BLAY-resistance properties it is required that within a
layer of each of the used bore profiles the resistance properties be either defined
or undefined. Otherwise the interpolated properties can have unpredictable values.
In other words, the input where the resistance properties within a layer for some bore
profiles are defined, while for the others are not defined, will most likely yield undesired
results.
BBAX
Hint
The ordinates S1 and S2 are measured along the axis of the bore profile, defined with the
BORE command.
The axial bedding describes the skin friction of the pile in dependence from the deformation
and from the lateral bedding force of the pile. A positive value of the load M0 acts on the pile
in the direction of the pile head.
skin friction
CO + TAND · SIG
k
−M0 displacement
Further explanations for the axial beddings are contained in the record BBLA.
BBLA
Hint
The ordinates S1 and S2 are measured along the axis of the bore profile, defined with the
BORE command.
Elastic supports have many related parameters. Therefore those values are combined to spe-
cial property elements for a geometric line.
All the corresponding BBAX and BBLA records follow the BORE record in the order defined by
the s ordinate. All data for the s ordinate refer to the parametric system of coordinates. The
default is the global z axis.
The pile bedding at the beginning of the section is K0, and the one at its end is K0+K1+K2+K3.
The individual values correspond to constant, parabolic, linear and quadratic distributions.
The default value for S1 is the latest S2 value. The initial default is -999.99.
The factors for the variation along the periphery are effective in the four quadrants (angle of 0,
90, 180 and 270 degrees). The angle refers to the local z axis. For linear analyses the factor
(P0+P2)/2 is used for the principal bending (MY,VZ), while (P1+P3)/2 is used for the transverse
bending (MZ,VY).
P2
P1 y P3 P1 y
P3
P0
P0
z z
P2 P2
y y
P3 P3
P1 P1
P0 P0
z z
The form factor is generated from the fact that the acting bedding force pL per length is given
by the following simplified relation from the deformation uL :
Es
ϖpL = D · C · L = D · · L (3.135)
Deƒ ƒ
According DIN we have e.g. Deƒ ƒ = min(D,1.0). For more complex cases where the bedding
stress is not uniform but more like a cosine, there are of course other values possible. There-
fore several design codes recommend to use any value between 0.5 and 2.0 to get the most
unfavourable results. Thus SOFiSTiK will not change the prescribed values in any kind.
For the bedding in axial direction a similar form factor may be defined based on a shear mod-
Thus the factors cancel each other in general and it is sufficient to use the stiffness modulus
Es for the axial bedding as well. In most cases the maximum skin friction is the more essential
part of the relation. However some value has to be specified, otherwise there would be no skin
friction at all.
Further there is a a rather sophisticated approach for the interaction of both directions available.
The friction has very different causes:
The first part of the pressure is described by the vertical earth pressure and the horizontal
pressure coefficient. The second part is given by the elastic constants which consist of a
stiffness and a dilatation.
BTAB allows the definition of design soil pressures which might be specified instead of the
geotechnical design for simple regular cases.
• The user defines only one value for the soil pressure V1 without specifying any options for
OPT
• The user defines for OPT explicitly a literal for the table in DIN 1054: A6.1, A6.2, A6.5,
A6.6, A6.7 or A6.8 with optional appendix for the conistency (only for tables A6.6 - A6.8).
E.g.:
or
• The user defines a Header-Record with OPT B and the desired width values V1 to V10,
followed by an arbitrary number of soil pressures (OPT S) per depth D for the selected
width values. E.g.:
BTAB D OPT V1 V2 V3 V4 V5
- B 0.1 0.2 0.3 0.4 0.5
1.0 S 210 210 210 210 210
2.0 S 280 420 560 700 700
3.0 S 380 520 660 800 800
4.0 S 480 620 760 900 900
General
SMAT defines a link material with the corresponding properties (header record). The defini-
tion of further reaction specific characteristics including work laws is followed directly after by
the command SFLA; SMAT and SFLA together provide a complete definition of link material
properties.
LTYP
The following link types are available:
STD Standard spring material with the reaction types P, PT, M (cf. SARB TYPE).
IHNG Material for the definition of the force-displacement relationship of the beam end
reactions (Implicit Hinge) (cf. SARB TYPE).
PMM Special material definition for the consideration of the P-My-Mz interaction. Available
are the reactions analogous to IHNG; the reactions My and Mz are additionally
coupled by the interaction exponent α (see ALPH).
MTYP
The following material types are available:
PLAS Anisotropic elasto-plastic material law. For the initial loading the material reaction
is developed according to the defined work law (plastic hardening). Unloading and
reloading on the other hand are linear elastic actions and are defined by the unload-
ing/ reloading modulus Er of the work law. The anisotropic model distinguishes
tension and compression, with mutually independent hardening behavior in the re-
spective directions. The provided work law may be non-symmetric, featuring differ-
ent branches in tension and compression.
PISO Elasto-plastic material law with isotropic hardening. Similar to PLAS, except that the
evolution of plastic hardening strains is not differentiated according to tension and
compression; plastic strains developed during tension are effective with respect to
compressive hardening and vice versa.
This type of material is usually accompanied by a work law with symmetric branches
for tensile and compressive response – its functionality, however, is not restricted to
the symmetric case.
PKIN 5 Elasto-plastic material law with kinematic hardening. For the kinematic hardening
rule, the initial elastic range is preserved during plastification – in contrast to the
isotropic and anisotropic hardening cases. In this model, the hardening effect is
represented by a shift rather than by an extension of the elastic domain.
This type of material is usually accompanied by a work law with symmetric branches
for tensile and compressive response exhibiting positive final gradients – its function-
ality, however, is not restricted to the symmetric case.
HYPE Hyperelastic material law. The material reaction develops according to the defined
work law, both for loading and unloading processes. Since there is no differentiation
between loading and unloading paths, there is no plastic deformation – the material
behaves elastic.
The material types PISO and PKIN are not available for the link types IHNG and PMM.
For further information regarding the respective material type characteristics, please refer to
subsection 3.25.2.
5 Kinematic hardening is supported by ASE, only.
Hint
This definition alone does not cause the actual limitation of the corresponding permissible
reaction force. The actual reaction, and with it the possible limitation of the permissible
force, must be defined together with the corresponding work law.
The values for the tension P+ and compression capacity P- play a special role. With the
definition of these capacities it is possible to formulate other reactions R with respect to the
axial load of the link element. The input of the normal force dependent family of work law
curves (see SFLA LEV) for the reaction R describes an interaction between the actual acting
normal force and the reaction R ( P-R interaction). For P ≥ P+ or P ≤ P− due to the interaction
the reaction is R ≡ 0.0; for the normal force levels in between, the permissible reaction R is
described by the user specified family of the work laws.
The capacity values are taken into consideration for the IHNG and PMM link types.
ALPH
The interaction exponent α ≥ 1.0 describes the interaction between the My and Mz reactions
according to
α α
My Mz
+ ≤1 (3.139)
My,m Mz,m
For example, with ALPH = 1.0 a linear interpolation between My and Mz is carried out. With
increasing exponent ALPH, the fullness of the interaction curve increases as well; with ALPH =
2.0 the interaction curve corresponds to an ellipse.
The default value ALPH = 0.0 defines a decoupled behavior of the My and Mz reactions (
mathematically this corresponds to α → ∞ in Eq. 3.139)
The interaction exponent ALPH is only active for the PMM link type.
[Link]/english/spring/spring_law_1_plas.dat
[Link]/english/spring/spring_law_2_piso.dat
[Link]/english/spring/spring_law_3_pkin.dat
[Link]/english/spring/spring_law_4_hype.dat
use the same SMAT work laws, as shown in Figure 3.28. The only difference is that the work
law for PKIN has an increasing end-tangent, see remarks to PKIN in Figure 3.28. Normal
plastic work laws (MTYP PLAS+PISO+HYPE) should have a horizontal end-tangent to clearly
define the final strength (Figure 3.28).
The four material types PLAS, PISO, PKIN and HYPE result in a different response under a
single cyclic loading. Resulting spring force-displacement curves are illustrated in Figure 3.28.
p p y
P[kN]
P [kN]
tensile and compressive behavior, e.g.,
P[kN]
500
and compression.
7
200
200.0
(
(5)
0.0
-400.0
-200.0
0.0
200.0
400.0
0.0
200.0
400.0
[mm] [mm]
0.0
−400
−200
−400
−200
200
400
200
400
−200
-200.0
P
P P
P
−500
MTYP PLAS + PISO + HYPE MTYP PKIN
compressive strength 6 . The yield function reflects an according differentiation between tensile
and compressive behavior:
|σ| − ƒ
y,c εp,c for σ ≤ 0
ƒ :=
(3.140)
σ − ƒy,t εp,t for σ > 0
In the example (Figure 3.28), the initial tensile hardening (1) is followed by an unloading step
to σ = 0 (2). So far, no compressive hardening has been triggered. Therefore, the continued
unloading path follows the initial loading path defined by the compressive branch of the work
law (3).
Isotropic Hardening
In contrast to the anisotropic approach, plastic strains are effective for hardening irrespective
of their direction (compression/ tension). Tensile plastification, e.g., also induces compres-
sive hardening and vice versa.7 Consequently, the yield function adopts the same hardening
variable tensile and compressive behavior:
|σ| − ƒ
y,c −κp for σ ≤ 0
ƒ :=
(3.141)
σ − ƒy,t κp for σ > 0
κp := m εp,c , εp,t (3.142)
In the example (Figure 3.28), the initial tensile hardening (1) is followed by an unloading step
to σ = 0 (2). For compression, the so far developed hardening strain is also effective, causing
a correpondingly extended elastic domain (4).
tensile and compressive behavior, e.g., concrete. The work laws in this case are usually non-symmetric for tension
and compression.
7 Work laws for isotropic hardening are usually symmetric for tension and compression – functionality, however,
N[kN] N[kN]
(1)
spring element 10 spring element 10
300 300
200 200
100 100
(2)
0.0 [mm] 0.0 [mm]
0.0
0.0
−200
−150
−100
−50
−200
−150
−100
−50
50
100
150
50
100
150
200
−100 −100
−200 −200
−300
−300
(3) (4)
0.0
−200
−100
−100
−150
100
200
50
100
150
−100
(5′ )
−200
−200
(1′ )
(c) MTYP PKIN (d) MTYP HYPE
aims at describing this effect and defines plastic hardening as a translation of the initial yield
strength ƒy0 without changing the extension of the initial elastic domain. Mathematically, this
notion is adopted by introducing the so-called back stress q as a function of the plastic defor-
mation εp . 8
ƒ := σ − q εp − ƒy0 (3.143)
N[kN]
feder element 1001
400
200
0.0
[mm]
−100
0.0
100
200
300
400
−400
−300
−200
−200
−400
• The allowable stress space is limited by the final tangents of the defined work law.
• For sufficient plastification in positive or negative direction, the end tangent is reached.
Possibly, a stretching of the nonlinear portion of the work law beyond ƒy0 is performed,
siehe Figure 3.31.
Subsequently, we describe the response based on the provided work law according to ??.
Important is the first point in the work law (5) that defines the length of the unloading part.
The unloading part from point (1) does not end at σ = 0 but at point (6). From there first the
linear part is attached to reach (5’). Then the nonlinear part of the worklaw (5)-(1) is streched
(5’)-(1’) and added to reach the lower end tangent of the negative part of the worklaw at (1’)
(end tangent = green line in Figure 3.31).
For multiple cyclic loading with increasing load factor PKIN shows the typical kinematic hard-
ening force-displacement curve (see Figure 3.29).
Figure 3.30 shows the kinematic hardening for a curved work law without linear part (fy=0).
The stretching of the nonlinear part of the work law for a curved work law with linear part
(fy=60) is demonstrated in figure 3.31. In a first unloading the new zero point of the work law
is reached (red circle). If at that point the original work law curve would be added, you would
never reach the lower end tangent (two upper pictures). Therefore the nonlinear part of the
work law is streched to reach the lower end tangent (lower picture).
N[kN]
feder element 1001
300
200
100
0.0
[mm]
0.0
−200
−100
100
200
300
400
−100
−200
−300
HYPE Hyperelastic
The response is nonlinear elastic, without plastic deformation. Upon load reversal, HYPE
follows the same path as in the loading phase (7) (Figure 3.28).
P[kN]
400
200
0.0
0.0 [mm]
200
400
600
−600
−400
−200
−200
−400
P
P[kN]
400
200
0.0
[mm]
0.0
200
400
600
−600
−400
−200
−200
−400
P
P[kN]
400
200
0.0
[mm]
0.0
−600
−400
−200
200
400
600
−200
−400
P
Standard spring
TYPE Type of reaction LT4 !
P Axial force reaction
PT Transversal force reaction
M Axial moment reaction
General link
TYPE Type of reaction LT4 !
N Normal force
Vy Shear force Vy
Vz Shear force Vz
Mt Torsional moment Mt
My Bending moment My
Mz Bending moment Mz
Performance limits
PERF Performance point
General
The SFLA command complements the SMAT command and allows specification of nonlinear
characteristics and work laws for the individual reactions of the link material. The work law
describes the loading path under uni-axial monotonic loading (backbone curve). Material re-
action under unloading and reloading conditions takes place according to the chosen material
type (→ SMAT MTYP).
A work law consists of up to 20 (U,F) tuples, identifying points on the curve. These points
are entered by consecutive SFLA commands. Segments inbetween these points are linear
(polygonal shape); alternatively they can be chosen to be smoothened by cubic splines. The
input of of the points must be made in such a way that U is arranged in monotonically ascending
order.
Partial material safety factors are not considered automatically; they must be incorporated in
the specified work laws explicitly.
Hint
For strains beyond the defined range, a linear continuation using the tangent of the bound-
ing defined curve point is assumed (linear perturbation). This holds both for tensile and
(Hint continued...)
compressive branches of the curve; a continuation from the tension domain into the com-
pression domain and vice versa is not performed, however.
TYPE
The available reaction types TYPE depend on the chosen link type (→ SMAT LTYP).
Hint
The reaction type TYPE must be specified with the first command line, i.e. the one which
identifies the first point of the curve (work law).This definition is retained for subsequent
points, i.e., subsequent command lines. A subsequent command line with changed TYPE
marks the beginning of a new work law or the beginning of a new family of work laws (see
LEV)), respectively.
3.26.1 Work laws and nonlinear properties for standard spring elements
(→ SMAT LTYP STD)
kNm/ m2
M Axial moment 1156
Work law
U Displacement [mm] 1003 or rotation [mrd] 1004 ∗ !
F Area related force or moment, units according to TYPE ∗ !
(see above)
S Type of supporting point ∗ POL
POL Polygon point, sharp bend
SPL C2 -continuous point (cubic spline)
or stiffness kN/ m2 / m 1097 , kNm/ m2 / rd 1100
∗ -
Hardening modulus kN/ m2 / m 1097 , kNm/ m2 / rd 1100
SH ∗ -
FP Elastic (proportional) limit, units according to TYPE (see ∗ -
above)
Additional nonlinearities
kN/ m2
CRAC Optional tensile strength (cracking) 1155 -
kNm/ m2
1156
kN/ m2
COH Optional cohesion (for TYPE PT) 1089 -
Additional nonlinearities
Additional nonlinear characteristics can be specified for respective reaction types. These fea-
tures are described in subsequent subsections.
The following table provides an overview which additional nonlinear features are available de-
pending on reaction type, material type and shape of the work law.
Hint
For material type MTYP PKIN, additional nonlinearities are not feasible.
Initial slip is differentiated in positive (tension) and negative (compression) directions. The
magnitude of maximum slip is not affected by possible plastic deformation during the response
history, it remains at its initial value. Upon load reversal, slip may be activated again (cf.
Figure 3.32).
9 Combination of reaction type PT and PLAS is supported by FEABENCH, only.
10 Supported by FEABENCH, only.
11 For reaction type PT, GAPN ≡ GAP is always implied.
N[kN]
300
springelement 10
200
100
0.0
0.0
−100
[mm]
100
150
200
250
300
−50
−100 50
−200
−300
Figure 3.32: Spring work law with GAP > 0 and GAPN = GAP
Crack model
An optional tensile strength can be provided independent of the work law via definition of a
corresponding CRAC property for the respective reaction.
• upon reaching the specified cracking strength, the transferable tensile force drops to zero.
• continued loading is accompanied by corresponding crack-growth.
• upon load reversal, the crack-opening must be completely reversed before compressive
forces can be transfered again.
Hint
Coulomb friction model and crack model are mutually exclusive, i.e., it is not feasible to
(Hint continued...)
simultaneously define cracking strength and friction parameters for reaction type PT.
Work laws are entered in a similar fashion as for standard springs, including the possibility of
refering to special predefined shapes (cf. subsubsection 3.26.1).
LEV
The formulation of a normal force-dependent reaction R (P-R interaction) is possible by defining
a family of work laws; each curve within this family being characterized by a specifically defined
normal force reference level LEV12 . The input of the reference level is carried out relative to the
link capacity SMAT P+/ P- (in [%] or as factor [−] ).
1 Inthe case that corresponding link capacities have been defined via SMAT, a relative input of F in [%] oder
[−] is recommended. Doing so, a value < 0.0 refers to the MIN-values while a value > 0.0 refers to the MAX-
values of the capacities corresponding to the respective type of reaction.
12 The reference level LEV must be specified at the input of the first point of each family of curves. This definition
is retained for the subsequent points of the curve, until a change of LEV occurs. This change marks the beginning
of a new curve for the family of curves for the present TYPE. A change in the reaction type TYPE completes the
definition of the family of curves and marks the beginning of a new curve (family) for the new TYPE.
Hint
The input for the definition of a family of curves must be arranged in a monotonically
increasing order for normal force reference level LEV.
The pictures 3.33 and 3.34 show for example the definition of a family of curves for the reaction
My and the resulting P − My interaction diagram.
MyM[kNm] S U
400.0
S U
S U
S U
200.0
S U
S U
-U -S 0 S U
0.0
-U -S
phi-yM[mrad]
-20.0
-10.0
0.0
10.0
20.0
-U -S
-U -S -200.0
-U -S
-U -S
-U -S
-400.0
P-level=M -1000.0[kN] P-level=M -980.0[kN]
P-level=M -950.0[kN] P-level=M -700.0[kN]
P-level=M -546.0[kN] P-level=M -200.0[kN]
P-level=M -70.2[kN] P-level=M 0.0[kN]
P-level=M 100.0[kN]
P [kN]
-400.0
-200.0
200.0
400.0
0.0
My [kNm]
0.0
-500.0
-1000.0
phi-y,pl= 0.0[mrad] phi-y,pl= 0.8[mrad]
phi-y,pl= 2.3[mrad] phi-y,pl= 5.4[mrad]
phi-y,pl= 11.6[mrad] phi-y,pl= 24.0[mrad]
Performance point
U Displacement [mm] 1003 or rotation [mrd] 1004 ∗ !
F Performance criterion label LT4 -
Table continued on next page.
To categorize and simplify assessment of the link reaction (Performance), a work law along the
deformation axis can optionally be divided in arbitrary sections, also known as Performance
Intervals. The subdivision takes place by the definition of Performance Limits, which mark the
interval boundaries.
A set of performance limits is associated with the preceding work law; it is performed by enter-
ing the corresponding tuples of U and F values in consecutive commands:
Just as for the work laws themselves, input of the performance limits within a set must be
arranged according to U in monotonically ascending order. A new set of performance limits
must be marked by a corresponding TYPE PERF in the first command line.
The following input example sets 5 performance limits for the previously defined My work law13
(see also Figure 3.33).
!Link material:
SMAT NO 112 LTYP IHNG P+ 100 P- -1000 MY+ 400 TITL "Plastic Hinge"
...
!Family of curves for reaction My:
...
!Curve for P-level 0.2*(P-) = -200kN:
SFLA S 400000 SH 20000 FP 0.8[-] F 1.0[-] LEV -0.2[-] TYPE MY
• Generally, for each work law of a family of curves corresponding performance limits can be
specified.
• Within one family of curves, the consistency of the performance limits must be ensured, i.e.
– The number of limits with U < 0.0 must be the same for all work laws of the family of
curves.
– The number of limits with U > 0.0 must be the same for all work laws of the family of
curves.
– The limit identifiers (labels) are the same for all the work laws of the family of curves
(the first set of performance limits defines the identifiers for the entire family of curves).
LCA
The definitions for the life cycle assessment can be done for individual materials with the record
LCA. Currently, it is possible to enter the global warming potential, which is defined via the CO2
equivalent per material.
The corresponding phases of the life cycle of a building are defined with PHAS. It is possible
to input the values individually for each phase or as a sum for several phases.
The basic unit for the global warming potenial (GWP) is kgCO2e = kilogramm CO2 equivalent.
The values for the individual materials can be defined per mass or per volume. The default for
all steel materials is kgCO2e/kg or kgCO2e/lbm and for concrete kgCO2e/m3 bzw. kgCO2e/ft3.
With a restart CTRL REST 2 it is possible to define the values of the life cycle assessment in
a separate AQUA run for already existing materials. Please note that only the last calculated
input is valid. This means that existing values are overwritten at restart.
With the following input variants it is possible to delete the material values of the life cycle
assessment from the CDB:
NO Cross-section number − 1
MNO Material number or preferred beamtype −/ LT 1
CENT centric beam
BEAM excentr. beam (Reference axis)
TRUS only truss (no bending)
CABL only cables
m2
A Cross section area 1012 1.0
m2
AY Shear area for y 1012 -
m2
AZ Shear area for z 1012 -
m4
IT Torsional moment of inertia 1014 *
m4 A3 /12
IY Moment of inertia y 1014
m4
IZ Moment of inertia z 1014 IY
m4
IYZ Moment of inertia yz 1014 0
m6
CM Warping modulus 1016 0
YSC Coordinates of shear centre [mm] 1011 0
ZSC relative to the elastic centre [mm] 1011 0
This record allows the input of cross sections without the corresponding geometric data, which
are necessary of course in detailed stress analysis, yield zone theory or reinforced concrete
dimensioning. Thus they can be used only for the static analysis.
With NO and a Literal for MNO you may also subsequently specify which element type should
be selected for elements with automatic type selection with that section. This definition can be
redefined at any time for any existing section. All other input values will be ignored in that case.
The default for YS is 0.0. If the ordinates YMIN and YMAX are defined, YS is in the middle
between YMIN and YMAX. The default of ZS is defined analogously.
If IT is defined as zero, special attention should be paid so that the torsional degree of freedom
does not lead to undefined rotation capability during the assembly of the total static system
(Error message: Parts of the system can move freely.).
The default for IY is equivalent to a rectangular section with a width of 1 m and the given area
A.
In accordance with Saint Venant’s estimate, the default value for the torsional moment of inertia
is
A4
T = (3.145)
4 · π 2 · (y + z )
The defaults for ymin up to zmax assume as a first guess a rectangular cross section and apply
then appropriate corrections from the radius of gyration.
SVAL can also be used for defining a modified cross section. This can be done either by
using a negative NO to modify an already defined cross section, or by making a copy of an
existing cross section by means of a negative MNO. The values A through CM may be then
then defined with a unit of [−] ,[0/ 0] or [0/ 00] as factors for the corresponding values, and
are thus preset to 1.0[−] . The new cross section has no geometric properties any more.
Cross section 1 receives 50% of the torsional moment of inertia. The geometry of the cross
section gets erased. Cross section 3 has 50% of the torsional moment of inertia of cross
section 2 and a slightly higher Iz but no geometry. Cross section 2 was not modified. Cross
section 4 is a cross section with IT reduced by half and with complete geometry. (Only possible
with AQUA license)
Depending on the definition of values one of the following section types is generated:
In the case of a T-beam, the total width of the plate BDL can be defined additionally in order to
design it as a single beam (without T-beam philosophy). The increased dead weight is taken
into account here and the greater width can be used for load applications. Only the widths BO
or B are used for the static values in order to obtain realistic shear areas.
BO BO
SO
HO
HO
H
SU
B B B
When nothing is input for REF, the origin of the coordinate system of the cross section is
assumed to be at the elastic centre. This is of relevance for beams with a reference axis. The
required dimensions of the cross section can be calculated by AQB. For this task, B or H can
be input negative when only that dimension should be changed.
The distribution of the reinforcement is controlled with the option RTYP, the diameter DASO/U,
the distance A and some entries for maximum distances between bars of the selected INI-File
of the design code. A definition of the total numbers of bars is possible by the usage of the unit
[-] with ASO resp ASU.
Columns:
CU ASO as layer 0 at all 4 sides with bar spacing (circumferential)
SYM ASO as layer 0 each at upper and lower side, however if the distance is
greater than MaxBarDistanceC (300 mm), intermediate bars with layer 3 will
be added
SZM ASO as layer 0 each at right and left side, however if the distance is greater
than MaxBarDistanceC (300 mm), intermediate bars with layer 3 will be
added
CORN ASO as layer 0 concentrated in the corners, bars between according to the
selected design code
CORN:n as CORN, but bars in corners with multiplicity of n (1≤n≤7)
CORN:nZ as CORN:n, but oriented left/right
p
CORN:nB as CORN:n, but with bundles (bars with a larger equivalent diameter D · n)
In all cases the diameter will be preset to the smallest allowed value from the INI-File (MinBar-
DiameterC = 12mm), and the absolute basic reinforcement according to the number of bars
(total minimum of 4) with that diameter.
As the reinforcement is saved as single bars, there will be some maximum allowable distances
AMAX established in the design codes for columns, bending members and torsional members.
The bars in the corners are positioned with the minimum distance AMIN, along the layouts with
A. Unless AMAX is defined quite large, intermediate constructive bars will be inserted with that
maximum distance in between.
Beams:
ASYM two main layers ASU at lower (1) and ASO at upper (2) side, optional inter-
mediate bars in layer (3) at the sides
ASZM similar to [ASYM] but right (1) and left (2)
If the distance of the reinforcement layers (1 below) and (2 above) becomes greater than the
limit MaxBarDistanceT (350 mm), additional reinforcement at the sides of the web will be intro-
duced with layer (3). Although this is mandatory for torsion, you might suppress this by entering
a zero value for IT or a large value for AMAX, but then these sections will not be designable
for torsion any more. This layer (3) will be fixed unless a "‘:F"’ is appended to the RTYP literal,
then it might be increased for biaxial bending design.
The cover of the reinforcement from the side edge is preset to the minimum cover from the
upper or lower edge, but not larger than one-fourth of the width. It may be changed with item
SS.
Please note, that DIN 1045-1 uses the cover of the compressive reinforcement (effective dis-
tance - D/2) as a limit for the lever arm during the shear design of the cracked section. Thus
the diameter has always an effect on the shear design.
For ASL a minimum value for the shear links may be specified. The default is zero for sym-
metric reinforcements (compression members) and taken from the design code otherwise. As
some design tasks require to distinguish between the minimum and the provided shear rein-
forcements, the definition should specify only the minimum shear reinforcements. The value
may be specified in cm2 / m 1021 or cm2 / m2 1022 or as ratio in [-].
MRF = 0 must be specified for unreinforced concrete sections. The input of MRF or MRFL is
not allowed for steel or timber cross sections.
The torsional inertia and elastic torsional stress is always calculated based on the distribution
of the torsional moment on the web with the full height and the plate with a reduced width. To
avoid constraint torsion for concrete constructions the effective torsional moment of inertia can
be reduced by IT explicitly or with the unit [%] to a relative value. A value of 0.0 is allowed but
may lead to kinematic systems.
For the equivalent hollow cross section used in computing the torsion reinforcement only the
web is used. The thickness of the equivalent hollow cross section BEFF is set according to the
selected design code:
• EN 1992-1-1 and National Annexes: The width is taken as the double value of the rein-
forcement cover, but not shorter than A/U in accordance with EN 1992-1-1 6.3.2 (1).
• DIN EN 1992-1-1/NA, DIN FB 102: The width is taken as the double value of the reinforce-
ment cover in accordance with DIN EN 1992-1-1/NA 6.3.2, DIN FB 102.
• OEN EN 1992-1-1/NA: The width is taken with a sixth of the diameter of an inscribed circle
in accordance with F. Leonhadt und E. Mönning: Vorlesungen über Massivbau, Erster Teil:
Grundlagen zur Bemessung im Stahlbetonbau, Kapitel 9.3.3, Seite 233-235, but not shorter
than the double value of the reinforcement cover in accordance with OEN EN 1992-1-1/NA.
• SIA 262: The width is taken to at least an eighth of the diameter of an inscribed circle in
accordance with SIA 262.
If the compressive stress in the concrete is exceeded, a larger value can be specified with
BEFF.
If SPT is defined with a value +8, and the reinforcement in the plate is specified within the outer
third, then a second layer is created in the inner third and the torsion is distributed on web and
flanges.
By default, shear deformation areas are always computed – independent of the material type –
and accounted for consistently for beam analysis according to Timoshenko-theory. The shear
correction factor for the shear deformation areas is 5/6 of the area for rectangular cross sec-
tions. For T beam cross sections, a numerical integration is occured according to the Boundary
Elements Method (CTRL STYP BEM 3). This behavior can be overruled by explicit specifica-
tion of the respective shear deformation areas Ay and AZ with a value in m2 or by a factor
with the explicit unit of [%] . Specifying 0.0 deactivates the effect of shear deformation areas
for the beam analysis.
The design of the transverse shear takes place at the most unfavourable location (neutral axis
or connection to the plate). A definition of SPTS 3 selects one of the flanges, SPTS 5 or 7 both
flanges for a design. If a definition of the interface is made, the design of the interface will be
done always.
cm2 /[−]
ASO Basic reinforcement position -Z or number of bars 1020 *
cm2 /[−]
ASU Basic reinforcement position +Z or number of bars 1020 *
DASO Diameter of top reinforcement [mm] 1023 *
DASU Diameter of bottom reinforcement [mm] 1023 DASO
DASS Diameter of side reinforcement [mm] 1023 DASO
A Distance of reinforcements [mm] 1011 *
FIXS Fixing of the side reinforcement LT YES
YES No increase of the side reinforcement dur-
ing design.
NONE Possible increase of the side reinforcment
during design
PASU Additional lower (or upper) plate reinforcement for tor- LT NONE
sion
YES With
NONE Without
cm2 / m
ASL Area of shear link reinforcements 1021 *
INCL Inclination of shear links cot|deg 0
Similar to the standard T-beam in record SREC, an edge beam section can be easily defined
with EDGB. However, this is stored in the database as a polygonal cross section and calculated
with CTRL STYP BEM 3 for the shear stresses.
When nothing is input for REF, the origin of the coordinate system of the cross section is
assumed to be at the elastic centre. This is of relevance for beams with a reference axis.
The distribution of the reinforcement is controlled with the diameter DASO/U, the distance A
and some entries for maximum distances between bars of the selected INI-File of the design
code. A definition of the total numbers of bars is possible by the usage of the unit [-] with ASO
resp ASU. If the distance of the reinforcement layers (1 below) and (2 above) becomes greater
than the limit MaxBarDistanceT (350 mm), additional reinforcement at the sides of the web will
be introduced with layer (3). Although this is mandatory for torsion, you might suppress this by
entering a zero value for IT, but then these sections will not be designable for torsion any more.
This layer (3) for the side reinforcement will be fixed. With FIXS NO it might be increased for
biaxial bending design.
The cover of the reinforcement from the side edge is preset to the minimum cover from the
upper or lower edge, but not larger than one-fourth of the width. It may be changed with item
SS.
For ASL a minimum value for the shear links may be specified. The default is zero for sym-
metric reinforcements (compression members) and taken from the design code otherwise. As
some design tasks require to distinguish between the minimum and the provided shear rein-
forcements, the definition should specify only the minimum shear reinforcements. The value
may be specified in cm2 / m 1021 or cm2 / m2 1022 or as ratio in [-].
MRF = 0 must be specified for unreinforced concrete sections. The input of MRF or MRFL is
To avoid constraint torsion for concrete constructions the effective torsional moment of inertia
can be reduced by IT explicitly or with the unit [%] to a relative value. A value of 0.0 is allowed
but may lead to kinematic systems.
For the equivalent hollow cross section used in computing the torsion reinforcement only the
web is used. The thickness of the equivalent hollow cross section BEFF is set according to the
selected design code:
• EN 1992-1-1 and National Annexes: The width is taken as the double value of the rein-
forcement cover, but not shorter than A/U in accordance with EN 1992-1-1 6.3.2 (1).
• DIN EN 1992-1-1/NA, DIN FB 102: The width is taken as the double value of the reinforce-
ment cover in accordance with DIN EN 1992-1-1/NA 6.3.2, DIN FB 102.
• OEN EN 1992-1-1/NA: The width is taken with a sixth of the diameter of an inscribed circle
in accordance with F. Leonhadt und E. Mönning: Vorlesungen über Massivbau, Erster Teil:
Grundlagen zur Bemessung im Stahlbetonbau, Kapitel 9.3.3, Seite 233-235, but not shorter
than the double value of the reinforcement cover in accordance with OEN EN 1992-1-1/NA.
• SIA 262: The width is taken to at least an eighth of the diameter of an inscribed circle in
accordance with SIA 262.
If the compressive stress in the concrete is exceeded, a larger value can be specified with
BEFF.
If PASU YES is defined and the reinforcement in the plate is specified within the outer third,
then a second layer is created in the inner third and the torsion is distributed on web and
flanges.
By default, shear deformation areas are always computed – independent of the material type –
and accounted for consistently for beam analysis according to Timoshenko-theory. The shear
correction factor for the shear deformation areas is 5/6 of the area for rectangular cross sec-
tions. For T beam cross sections, a numerical integration is occured according to the Boundary
Elements Method (CTRL STYP BEM 3). This behavior can be overruled by explicit specifica-
tion of the respective shear deformation areas Ay and AZ with a value in m2 or by a factor
with the explicit unit of [%] . AYZ is also determined automatically, unless AZ and / or AZ was
not specified manually. The shear area AYZ can also be negative under certain circumstances.
An explicit input is only possible here in the unit m2 . Specifying 0.0 deactivates the effect of
In order to check the shear stresses due to shear forces, a shear cut is created in the web
of the edge beam section with default of CUTZ H/2 for the z position. The position of this
shear cut can be modified by the user with an input of CUTZ and can be adapted to the most
unfavourable position of the shear stess. With CUTF YES a shear cur in the flange is defined
for the check of flange connection. If a definition of the interface between web and plate is
made, the design of the interface will be done always.
For the standard edge beam, the control options STEU FACE and STEU BEWQ can also be
used.
If further modifications are required for this cross section, it can be exploded in the program
SOFiPLUS, i.e. it can be divided into the components of the corresponding polygonal cross
section. The polygonal cross section can then be graphically edited and modified using the
Cross Section Editor.
Base reinforcement
A Distance of reinforcements [mm] 1011 150
D Diameter of base reinforcement [mm] 1023 10
RTYP Reinforcement of the edges LT SYM
SYM symmetric reinforcement
ASYM asymmetric reinforcement
explantion below
Edge reinforcement Z+
LZ+ Length of edge reinforcement Z+ [mm] 1011 -
AZ+ Number of bars Z+ (one side ) - 3
DZ+ Diameter of edge reinforcement Z+ [mm] 1023 D
cm2 1020
RZ+ Maximum total edge reinforcement Z+ -
Edge reinforcement Z-
LZ- Length of edge reinforcement Z- [mm] 1011 LZ+
AZ- Number of bars Z- (one side ) - AZ+
DZ- Diameter of edge reinforcement Z- [mm] 1023 DZ+
cm2 1020
RZ- Maximum total edge reinforcement Z- RZ+
The record SHRW is used to create simple cross sections as shear walls for stability and
earthquake checks including design.
According to the input, the cross section can receive only a base reinforcement in layer M0 or
also additional edge reinforcements. The reinforcement distribution of the edges is determined
with the literal for RTYP:
Apart from the material number for the shear reinforcement, no further information on shear
design (e.g. minimum shear reinforcement) is given. Torsion is not taken into account.
The location of the local coordinate system lies in the elastic centre.
SA
D T
S
AS ASA
Record QC is intended for the definition of reinforced or timber sections. It is not suited to
define steel sections. These should be defined with PROF and the types BAR. COLD.
The distance SI is taken relativ to the inner radius for an annular section, but on the outer radius
for a solid section. The basic reinforcement may be also specified with explicit unit in [cm2 /m].
There will be at least four locations for the reinforcements. The default will provide six bars, the
number is controlled with item A. If that value is defined negative or with the explicit unit [-], the
given value is taken as the number of bars directly.
Bending members do not have a minimum longitudinal reinforcement, while columns do not
have a minimum shear link area. If a bending member with a longitudinal basic reinforcement
is required, an explicit value for A will trigger this. As the provisions of most design codes are
not covering circular sections, any minimum shear reinforcement has to be selected explicitly
via ASL in cm2 / m 1021 The default of diameter of reinforcement DAS is depending on chosen
design code.
The definition of CTRL RFCS in AQUA is also effective for the SCIT section.
Values for IT, AY and AZ are either absolute values or if defined with unit [%] as factors to the
theoretical values.
NO Section number − 1
D Nominal diameter [mm] -
TYPE Type of cable section (see remarks) LT -
INL Type of Inlay LT FE
FE, FEN, FEC = Fiber Inlays
SE, SES, SEL = Steel Inlays
MNO Material number of a prestressing steel for cables − 1
record STEE TYPE YC
The ultimate forces of the cables are obtained only based on the tensile strength values given
above)
FE = Fibre inlay
FEN = Natural fibre inlay
FEC = Chemical fibre inlay
SE = Steel inlay
SES = Steel cable inlay
SEL = Steel strand inlay
In the Eurocode EN 12385-4 appendix A to C the definitions of the factors are given:
In DIN 3051 part 3 the factors have been defined in a different way as:
2
Metallic cable section qm = ƒ · d 4· π
Minimum ultimate force Fmn = k · Fr = k · qm · σz
Weight per length G = qm ·
The loss factor ke taking into account the type of fixing of the cable endings is defined else-
where e.g. in EN 1993-1-11 Tab. 6.2 or DIN 18800. The formula given in EN 12835-4 appendix
A.2 for the ultimate force for cables with diameters larger than 60 mm is not applied. It does
not contain a strength and is not applicable in general cases.
The user has to check all factors in detail, as they depend on the type of anticorrosive lining and
the intended usage (e.g. for cableways). Cable sections may only be used for cable elements
or automatic elements, which derive their type from the cable section type.
The safety factor 1.5 according to EN 1993-1-11 equation (6.2) is considered automatically.
The partial factor γR according to EN 1993-1-11 Tab. 6.2 is defined at GAMR. The determined
tension resistance Frd is the minimum value of the equation (6.2) EN 1993-1-1.
See also: CS, INTE, SV, POLY, CIRC, PANE, PLAT, PROF, Reinforce-
ment, CUT, SPT, SFLA, WIND, WPAR
SECT
Freely defined cross sections always begin with the record SECT, which defines the cross sec-
tion number. All subsequent input records describe this one cross section, which may consist
of several partial cross sections (external outline, inner perimeter, reinforcement arrangement
etc.). The input for each cross section ends either by the next SECT record, which can be
empty, or by the END record.
Freely defined cross sections are divided into three groups (see 1.1):
A coordinate system y - z is established for every SECT definition, the origin of which is in
general on the reference axis defined by the two nodes of a beam. All coordinate data of the
input records which follow a SECT record refer to this coordinate system. The directions of the
axes are identical to those of the local beam element and defined in accordance with Chapter
2.1 (y to the left, z downward). The local coordinate axes y’ and z’ of a centric beam are only
shifted parallel to the cross section coordinate system, so that the origin is at the elastic centre
of the cross section.
Prestressing Tendons (AQBS, GEOS) always refer to the input coordinate system. It may be
appropriate to take this into account when selecting the zero point.
If desired, however, it is also possible to rotate the cross section by any angle or into the
direction of the principal axes (ALPH ’S’) as well as to translate it (items ALPH, YM and ZM).
For all rotations, it is not the reference coordinate system but the cross section elements that
are rotated. The definition of ALPH will force AQB to do the stress analysis for the rotated
principal axis system. The input of CTRL AXIS -2 within AQB has then no effect on this section
any more.
On the other hand a definition of FSYM YES will suppress the rotation of the principal axis
completely. This is intended for sections describing only half of a full symmetric section, but
may have very dangerous effects if applied to general sections. The value of Iyz has consider-
able effects on deformations and also forces within constrained systems. This option enforces
also uni-axial bending (Vy=0, Mz=0) within AQB.
With BTYP you may specify your preferred beam element type for that section. This info may
be used for the mesh generation and the specification of minimum reinforcements.
For the design of a reinforced concrete cross section with AQB, it is always necessary, to
specify the location of the intended reinforcement - single, linear, or perimetric reinforcement-
by means of the records RF, LRF, CRF or CURF. Freely defined cross sections cannot be
dimensioned with STB/STBA. MRF defines the material number of reinforcement. Only if the
link reinforcement has a different quality of material a combined value is to be entered. Then
this combined value represents the default for the records CUT.
The material number should, in general, be specified by SECT. The declaration of a material
number with individual cross section elements is only appropriate for composite cross sections.
In case of composite sections, ideal cross section values are calculated, based on the material
(A · E)
A = (3.146)
Ereƒ
All cross section elements are addressed with an arbitrarily selected identification number,
which has up to four characters in general. In AQB and ResultViewer you may specify a mask
to select specific elements for the output. You might for example then select all elements with
a zero at the end.
The buckling strain curve can be input for stability checks. The permissible input values are 0
(none), a, b, c or d. Appropriate technical knowledge is required for making this choice in case
of general cross sections. AQUA, attempts to model most cases with the following defaults:
strong weak
axis axis
Profiles without welding joints
Annular and SH-shapes a a
U, L and solid circle shapes c c
Double T-shapes h/b > 1.2 t ≤ 40mm a b
t > 80mm d d
others b c
all others
I-strong > 1.67 · I-weak and
I-t < 0.50 · I-weak
tmax ≤ 40 mm b c
tmax > 40 mm c d
I-strong < 1.67 · I-weak or
I-t > 0.50 · I-weak
tmax ≤ 40 mm b b
tmax > 40 mm c c
The buckling curve ”e” may be input to select the old AISC-curve with the Euler hyperbola for
λ > 0.5 and a quadratic parabola for the plastic region. The safety factors have to be large
enough for this curve!
For ultimate plastic design according EN 1993-1-1 the interaction of the bending moments
requires to exponents α and β. As these factors depend strongly on the type of the section, it
is possible to define them explicitly or to chose from the two predefined formulas.
Hint
For some records (e.g. TVAR and PROF) it makes a difference if those records are
(Hint continued...)
defined within a section or separately. A definiton of SECT 0 will allow to terminate the
current section.
Primary solution for that task are formula expressions to be defined for any coordinate or per-
manent radius (CIRC, CRF) as well as the thicknesses of thin-walled elements and their shear
connections with up to 256 characters in the form of ”=formula”. These formulas will be saved
with the section and may be reevaluated for any section with different values along an axis (see
GAX/GAXP) or with explicit definitions locally with TVAR commands.
Formulas should use variables with a designated unit, only then the interpretation of #B + 20
can be handled in a unique way. Simple cases may be also defined in the default sectional
units. As an alternate possibility the unit of the formula may be appended, but in that case all
the used variables must not have assigned a unit.
Then we may define coordinates relative to up to three other reference points. As a reference
you may use every stress point, polygon vertex, panel or circular element. The reference
is done via this explicit 4-character identifier. The specified coordinates are relative to that
reference. However if the coordinates are not a formula and either CTRL REFD 0 has been
specified or the reference is preceeded by an @ the coordinates are taken as absolute, which
is used for the export of sections.
It is also possible to specify the position of the reference point by an axis, for that case the
reference has to be specified as the ID of that axis with a colon prefix, e.g. ”:AX_0”. The local
coordinates (y,z) are then given by the 3D-distance to the same parameter on the reference
axis projected in the plane of the section.
If points with multiple identifiers are present, then only the first occurrence of that point is used
and the others are neglected. If a reference point is not part of the section itself, you should use
the material number 0 for it. However, you may have an arbitrary number of nested references,
i.e. a reference point may use a reference itself. For thin walled elements start and end point
may be addressed by a reference with an index.
As the references may be used on single points, start points and end points, the following
examples use a generic description. Thus REF holds for REFP resp. REFA or REFD.
PT PTZ
PT z
z
PT0 PTY
y y
• It is also possible to specify negative references. The coordinates will then be used with an
alternate sign, allowing easy description of [Link] the coordinates itself, the double
of the coordinate values of the mirroring center or line have to be specified. (yne =2•ymrr -
yorg )
PT0 PT0
The point PTD (at item REFD, RFDA or RFDE) defines the direction of the reference relative
to the PT0 point. Instead of a point it is also possible to specify the name of a variable or a
formula containing the angle of the direction in radians: ”=#VARNAME”
The y-coordinate is then measured in the radial direction along, while the z-coordinate is per-
pendicular and positive to the left. The input is done via:
For a variant construction the (+) will maintain the ratio of distances in the radial direction and
maintain the value perpendicular to that.
The point PTD is defined as geometry point e.g. G1 outside the polygon and is referenced with
∼G1. The variable POL1 at the y coordinate of the polygon point P varies the length from y =
250 to y’ = 500 mm during the interpolation. The z coordinate of the geometry point G1 as well
as the point PT0 can also input with a variable TVAR along an axis (carthesian reference).
SECT 100 MNO 1 MRF 2 FSYM NO BTYP BEAM TITL "master 100"
TVAR 'POL1' VAL 0[mm]
POLY TYPE O MNO 1
VERT '0100' Y -500 Z 0 EXP 1
VERT '0101' Y 500 Z 0 EXP 1
VERT '0102' Y "=#POL1,250" Z 1000 EXP 1 REFP '@0101' REFD '@~G1'
VERT '0103' Y "=#POL1,-250" Z 1000 EXP 1 REFP '@0100' REFD '@~G1'
SPT 'G1' Y 0 Z 2000 MNO 0
Example for polar reference with affine scaling only along direction:
For the variant with affine scaling only along direction the point PTD has to be defined here
as geometry point G1 which is variable along an axis. In this example the geometry point
G1 moves upwards in z direction. The referencing is done with +G1. For the extension in
longitudinal direction the y coordinate of the polygon point P has to get an information for the
extension, here with the variable POL3 with 250 mm at the start and 500 mm at the end.
The length of the polygon edge y changes to y’ according to the distance ratio:
y/(PT0-PTD) of the first cross section to y’/(PT0-PTD’) of the last cross section = 0.5 resulting
from the ratio variable POL3 of the first cross section to last cross section 250/500 = 0.5.
SECT 300 MNO 1 MRF 2 FSYM NO BTYP BEAM TITL "master 300"
TVAR 'POL3' VAL 0[mm]
TVAR 'G13' VAL 0[mm]
POLY TYPE O MNO 1
VERT '0100' Y -500 Z 0 EXP 1
VERT '0101' Y 500 Z 0 EXP 1
VERT '0102' Y "=#POL3,250" Z 1000 EXP 1 REFP '@0101' REFD '@+G1'
QP '0103' Y "=#POL3,-250" Z 1000 EXP 1 REFP '@0100' REFD '@+G1'
SPT 'G1' Y 0 Z "=#G13,2000" MNO 0
As in the previous cross section the geometry point is variable along the axis and the variable
POL4 is defined for the extension in longitudinal direction. Additionally, for the extension in
transverse direction the variable POLZ is defined. This is input for the z coordinate of the
polygon point P. The reference to the geometry point G1 is done with *G1 for the extension of
both directions.
As result the width z is changing to z’ according the distance ratio of the variable POLZ of the
first cross section to the last one 100/200 = 0.5 additionally to the length of the polygon edge y
to y’ as described for the previous cross section.
SECT 400 MNO 1 MRF 2 FSYM NO BTYP BEAM TITL "master 400"
TVAR 'G14' VAL 0[mm]
TVAR 'POL4' VAL 0[mm]
TVAR 'POLZ' VAL 0[mm]
POLY TYPE O MNO 1
VERT '0100' Y -500 Z 0 EXP 1
VERT '0101' Y 500 Z 0 EXP 1
VERT '0102' Y "=#POL4,250" Z "=#POLZ,1000" EXP 1 REFP '@0101' REFD '@*G1'
VERT '0103' Y "=#POL4,-250" Z "=-#POLZ,1000" EXP 1 REFP '@0100' REFD '@*G1'
SPT 'G1' Y 0 Z "=#G14,2000" MNO 0
PT0 PT
If three points are given, then the third point may be provided with a prefix operator defining
the distance or elevation to be taken from that point and searching the corresponding point on
the line REF-RFD:
If for the first two cases the third point is specified as the point itself, this will create points with
the same selected coordinate, which is needed for example for points with a fixed distance, but
a height depending on the cross inclination.
For a circle you may specify an additional point. The distance of this point to the centre will then
specify the radius of the circle. For steel shapes the angle of orientation is used in a similar
way.
Example for a construction reference with distance (y) - cut with a vertical line:
The point PTD is defined as geometry point e.g. G1 outside the polygon and is referenced
with >G1. The value y of the polygon point P should be constant and is input with 250 mm.
However, it is also possible to define the coordinate y of the polygon point P variable along an
axis with TVAR.
Since the geometry point G1 shifts upwards with a variable height as a result of the Carte-
sian reference the cross section becomes smaller in height, while the width of the lower edge
remains the same.
SECT MNO 1 MRF 2 FSYM NONE BTYP BEAM TITL "master 100 vertical"
TVAR 'G11' VAL 0[mm]
POLY TYPE O MNO 1
VERT '0100' Y "=250,250" Z 1000 EXP 1 REFP '@0103' REFD '@>G1'
VERT '0101' Y "=-250,-250" Z 1000 EXP 1 REFP '@0102' REFD '@>G1'
VERT '0102' Y -500 Z 0 EXP 1
VERT '0103' Y 500 Z 0 EXP 1
SPT 'G1' Y 0 Z "=#G11,2000" MNO 0
Example of a construction reference with elavation (z) - cut with horizontal line:
The point PTD is defined as geometry point G1 outside the polygon and is referenced with
ˆ G1. The value z of the polygon point P should be constant and is input with 1000 mm.
However, it is also possible to define the coordinate z of the polygon point P variable along an
axis with TVAR.
Since the geometry point G1 shifts upwards with a variable height as a result of the Cartesian
reference the lower edge of the cross section becomes smaller, while the heigth of the cross
section remains the same.
SECT MNO 1 MRF 2 FSYM NONE BTYP BEAM TITL "master 200 horizontal"
TVAR 'G12' VAL 0[mm]
POLY TYPE O MNO 1
VERT '0100' Y 250 Z "=1000,1000" EXP 1 REFP '@0103' REFD '@^G1'
VERT '0101' Y -250 Z "=1000,1000" EXP 1 REFP '@0102' REFD '@^G1'
VERT '0102' Y -500 Z 0 EXP 1
The point PTD is defined as geometry point G1 outside the polygon and is referenced with G1.
For the definition of the perpendicular, a third point - in this case the opposite polygon point -
has to be selected for REFS. In this example asymmetric cross sections should be generated
with the interpolation. Therefore only one polygon point is referenced.
Since the geometry point G1 shifts upwards with a variable height as a result of the Cartesian
reference the line PT0 - PTD changes. Furthermore, since the default is that the line PT3 - P
is the perpendicular to PT0 - PTD, the polygon point P shifts in the interpolated cross sections.
SECT MNO 1 MRF 2 FSYM NONE BTYP BEAM TITL "master 300 lot"
TVAR 'G13' VAL 0[mm]
POLY TYPE O MNO 1
VERT '0100' Y -250 Z 1000 EXP 1
VERT '0101' Y -500 Z 0 EXP 1
VERT '0102' Y 500 Z 0 EXP 1
VERT '0103' Y 250 Z 1000 EXP 1 REFP '@0102' REFD '@G1' REFS '0100'
SPT 'G1' Y 0 Z "=#G13,2000" MNO 0
In addition of the examples explained here, more examples showing most of these features are
given with [Link] to [Link] in the example directory [Link].
Hint
The offset definitions YM and ZM for SECT and POLY cannot be used as parameters.
3.34.2 FE-Sections
For advanced design tasks like a hot design it is not sufficient to describe the section just by
its outer contour. For such cases it is necessary to us a finite element mesh of the section.
These meshes may be generated by AQUA or may be imported from a secondary database.
The secondary database will contain in general also temperature distribution fields. It is also
possible to clone a section from the same database, which is especially useful for sections
with different reinforcement layouts based on the same thermal master analysis, performed
only once. The defintion of the name ”SECT:nnn” will use a FE-mesh of this section nnn
directly. In that case an import of reinforcements via RF is also possible.
The section will be converted to single integration points given by the center of all QUAD ele-
ments. These may be imported group wise with POLY to specify any construction sequences
or in total if no such definition is made.
As this feature allows the evaluation of all sectional values for any type of section (e.g. sec-
ondary torsion for a solid section) it is the most general type of section definition. All additional
elements (Stress points, reinforcements, shear cuts) are defined as usual.
With the definition of LTAU the unit warping and the shear stress distributions will be saved to
the original database. So it becomes possible to view those results with WING/WINGRAF.
With CS you may specify a cross section for (up to 31) construction stages. All elements
following this record will be added within that stage. The current construction stage will also
contain all the elements of the previous construction stages.
If part of the section is active only temporarily, a value for ATIL may be specified. This is then
the last construction stage where this part is active.
Hint
The omission of cross section elements with ATIL is not recommended because this work-
flow is continuously not supported by all programmes.
The CS-stage records have to be given in monotonic ascending sequence of NO, but for this
case multiple records with the same NO value may be given. However the designation should
be given with the first record of such a CS block.
10 20 21 22 24 30 40 ∞
The picture above defines a construction process. As phase 21 has no impact on the section
itself it does not need to be defined. The same holds for any other construction phase like
pre-stress stages. The input scheme is given by:
NO is the total section for the final phase 40, QGroup "‘C040"’
NO.1 is the 1st construction stage, phases 10-19, QGroup "‘C010"’
NO.2 is the 2nd construction stage, phases 20-21, QGroup "‘C020"’
NO.3 is the 3rd construction stage, phases 22-23, QGroup "‘C022"’
NO.4 is the 4th construction stage, phases 24-29, QGroup "‘C024"’
NO.5 is the 5th construction stage, phases 30-39, QGroup "‘C030"’
The definitions of ATIL for phases 20 and 22 will be extended automatically to the next defined
construction stage, only for phase 24 the end value has to be specified explicitly. The transition
to the final stage may be defined either by an explicit construction stage 40 or by the latest
removal of any other phase. There is no generation of intermediate values, a CS 20 ATIL 33
will not generate a construction stage at 34.
The specified QGroups are the primary and secondary groups assigned to the elements in the
secondary database of the section (CTRL STYP FEMX only)
The construction stages are assigned to the individual elements with the group definition of the
analysis program, a construction stage number defined there is inserted between the defined
numbers here. The number NO will be incremented by default. You may want to use larger
gaps to allow prestressing stages to be mixed in. With the construction stages defined in AQUA
as 10, 20 and 30 the selection of stage 25 would use the 20-section and all tendons up to stage
25.
For every construction stage it is also possible to specify with SV for every material a factor
for the elasticity and the shear modulus after the CS record. Further it is possible to specify a
construction phase as hot design phase, if the section is imported with a HYDRA analysis.
m4 1014
IT Torsional moment of inertia -1.0
m2 1012
AK Area of Bredt’s box *
YSC Coordinates of the shear center in the LT/ ∗ *
ZSC reference coordinate system LT/ ∗ *
m6 1016
CM Warping modulus *
m4 1014
CMS Warping shear modulus *
m2 1012
AY Shear deformation area y *
m2 1012
AZ Shear deformation area z *
AYZ Shear deformation area yz *
LEVY Minimum lever arm for VY [mm] 1011 -
LEVZ Minimum lever arm for VZ [mm] 1011 -
In general, these cross section properties will be computed automatically. You only have to
specify explicit deviations.
AK defines the area of the equivalent hollow cross section according to Bredt. This value
is used for determining the longitudinal reinforcement and the link reinforcement due to tor-
sion. However, in general AK should be defined implicitly by specifying the torsion resisting
reinforcements.
An absolute value (0. or positive) can be input for IT. If a scalar value is input, the moment of
inertia computed by the program is multiplied by a factor. The value 0.5[−] or 50[%] results
in a cross section with 50 percent of the torsional moment of inertia.
YSC and ZSC can define a mandatory centre of rotation or, in case of only one value, a
restraining plane. The literal ’C’ can be used in both cases to set this ordinate to the value of
the center of the elastic centre.
If shear deformations are not to be taken into account, despite a detailed shear stress evalu-
ation, AY, AZ , AYZ and CMS should be set to ”0.” . The value 0.5[−] or 50[%] results in a
cross section with 50 percent of the calculated shear area.
For shear force dimensioning in state II, minimum lever arms can be specified with LEVY and
LEVZ, in order to obtain more economical results (e.g. 0.90d) or cover extreme cases (e.g.
moment=0). Values are definable as absolute values or as factors to the height or width of the
section.
For the analysis of creep and shrinkage effects an effective depth deƒ ƒ is required. This value
is defined by the area of the section A and the length of the periphery U which has air contact
by the formula 2A/U. The air contact ratio may be defined for the vertices of a polygon and for
circles. If not otherwise stated, outer peripheries will have a ratio of 1.0 and inner peripheries
of 0.0. However you may specify the value for every material within a section via record SV
explicitly.
For construction stages it is possible to define with SV MNO two factors for the elasticity and
the shear modulus. They will be mainly used for the evaluation of the sectional values. A
consistent treating in AQB for all effects is still under investigation.
In many practical cases it is recommended that the user defines his own plastic forces and
moments or limits the ones computed by AQUA if special considerations are to be taken into
account. The latter could easily be achieved either by a direct input of the plastic value that
shall directly overwrite the one computed by AQUA, or by defining a fractional part by means of
[ -] resp. [ %] in the frame of the CadInp language, thus scaling the default one. For example,
setting MBPL to 0.8[ -] would result in AQUA using 80% of the originally computed plastic
bimoment. The default for each entry is 1.0[ -] or 100%, i.e. unless explicitly defined by the
user, AQUA would not change the value of a plastic force or a moment.
It is worth noting that only characteristic plastic forces and moments could be defined or altered
in such a way. If so, AQUA will produce a corresponding design value by dividing the input value
by the material safety factor.
Unless a REC*-type has been selected, the record POLY must be followed by the input of the
polygon vertices with VERT . The sequence of the polygon (clockwise or counter-clockwise)
has no effect. The polygon will be closed by the program automatically. In case of symmetry,
the polygon is extended by mirroring before being closed.
The use of variables and formulas for interpolation of the cross-section is only possible for
polygons with geometry definitions in record VERT. Variables and formulas are not evaluated,
if symmetry or dimensions of the rectangles DY,DZ are specified.
A hole is automatically created when polygons or circles overlap, a true hole is thus defined by
a polygon or circle with material number 0. The definition of a special inner polygon with the
same material number or the formerly common method of combining several polygons (outer
perimeter and inner perimeters) into one single polygon by making two passes along the same
edge should be avoided.
The input of YM and ZM causes a corresponding shift of the given polygon. The coordinates
of VERT or DVER refer then to the shifted (by YM, ZM) coordinate system. The symmetry data
also refers to the shifted coordinate system. Thu, similar openings can be generated easily.
A rectangle can easily be defined by means of a special type. These define a rectangle with
sides DY and DZ, with its centre or upper or lower midside point at the coordinates YM,ZM.
AQB prints stresses only in the defined polygon or stress points. Thus a coarse sub division
of longer edges might be useful. This can be done with a value SMAX or as an alternate way
with the value PHI at VERT. However in this case it is recommended for symmetrical polygons
to define additional vertices on the symmetry axis to allow a subdivision of the closing edges.
For the REC*-types SMAX is preset to DZ/4.
The exposition class EXP allows the definition of special material parameters (MEXT) to indi-
vidual poygon edges. The value specified here becomes the default for that polygon.
NO is used for identification during any output of stresses. If nothing is input, AQUA generates
internal numbers in sequence.
The distances between adjacent polygon vertices must be at least 0.0001 m. The polygon is
defined by the sequence of the vertices, not by their numbers. The number of points is limited
to 255 per polygon.
For tangential points defined with TYPE TP a fillet is created at that vertex with the given radius
R. If the second radius R2 is specified, the fillet is elliptical. Without the radius R2 it is a circle.
If the radius R is defined negative, an isosceles chamfer is created instead with the value of R
used as a distance along the edges (see picture below). If R2 is also defined negative, then an
unequal chamfer is generated.
In the variant with two radii for the definition of an ellipse, the radii have to be specified accord-
ing to the direction of rotation of the polygon. In this case, the radius R is the radius before the
polygon point and the radius R2 is the radius after the polygon point.
Only this input option of a fillet is possible for a cross section interpolation with references along
an axis.
Without thew input TYPE TP a tangential intersection point (with a ’T’ at beginning of the name)
is inserted between this (’E’) and the previously defined vertex (’A’) (see picture below). With
an angle > 120 degree two tangential points (’R’ and ’T’) and an additional midpoint (’S’) are
inserted. An additional start point (’Q’) of rounding is inserted if the radius is less than the
half distance between the original polygon points. The inserted points get the reference of the
original polygon points. For this a new geometry point ’G’ with the original reference is inserted.
Then all new points are referenced to this geometry point.
In both cases the radius of the inserted points will be adjusted to achieve exactly the same
area as the circular arc respective as the elliptical arc.
Thus, a full circle section should be defined via CIRC or an explicit loop with trigonometric
functions.
Non effective areas can be defined by means of NEFF-areas, where parts of the polygons
within those areas become non effective for selected forces by the internal introduction of
deductional polygons.
The sectional values of the total sections are only used for the calculation of the area as well
as the torsional stress and the integral equation solution. All geometrical moments of inertia
are computed based on the effective parts only (refer to the AQB manual). It is to be noted that
the effective width is actually dependent on the load case and the on the purpose of design.
The input definition EXP allows the definition of special material parameters (e.g. air contact,
temperature of environment, cover of reinforcement, ... see MEXT) to individual poygon edges.
The input definition EXP is set at the first point of the [Link] can find further explanations
and examples in our Online Help -> Tutorials
[Link]
cross-sections-general/cross-section-air-mext/cross-section_air-[Link]
Instead of a literal at EXP, only the degree of air contact can be set. For the degree of air
contact of the whole edge the geometric mean value is used. I.e. if one of the two vertices of
an edge has the degree 0.0, the total edge will have this value. The values ’0.0’ and ’1.0’ are
defaults for concrete materials which are redefined in AQUA. Inner edges and holes are get
the value ’0.0’, outer edges are get the value ’1.0’.
R
y, z
+ϕ
Internal sequential numbering takes place when nothing is input for NO. NO can be selected
arbitrarily. When CIRC is overlapping another material, a whole of that material is automatically
calculated.
With a CUT one specifies a part of the section to be used for shear design and/or the minimum
reinforcement or crack width of a partial section. The input of shear shear cuts is necessary
for the checking of shear stresses in concrete cross sections. Without input for CUT, up to two
axis parallel shear cuts are defined through the elastic centre (see CTRL STYP). As these do
not necessarily pass through the smallest width of the tensile zone, a warning is issued when
the user enters complicated polygons without shear cuts. It is the responsibility of the user to
define enough shear cuts to account for all unfavourable locations within a section.
The TYPE of the cut defines some important properties for the design. The following values
are available:
• ACT Only for the crack width design, cuts off the cracked part
• WEB Web of a cross section
WRED Web with reduced allowed strength
• FLAN Flange of a cross section
FFUL Flange with enhanced strength allowance
• INTE Construction interface
INDE indented shear interface
ROUG rough shear interface
EVEN even shear interface
SMOO smooth (very even) shear interface
The coefficients CINT and MUE are defined in EN 1992 or DIN 1045-1 (2008). For design
according to the old DIN 1045-1 (ch.10.3.6) the value CINT has to be specified with the
value of βct > 1.
• SLAB For special cases (e.g. hollow plates) it is possible to design a cut of any section
like a plate.
In general the definition of WRED / FLAN is expected to describe structured sections, while
WEB / FFUL is defining a compact section. There are some design tasks requiring the knowl-
edge of the minimum web thickness, which is found by scanning the shear cuts of type WEB.
For DIN 1045 (1978) this info is used to allow the shear region 3. However in this case, the
For DIN EN1992-1 this info is used to distinguish for the superposition of shear. WEB and
FFUL use a quadratic superposition for open and compact sections as specified in the national
annexe. WRED and FLAN use a linear superposition for closed sections according to original
equation 6.29.
AASHTO 2005 distinguishes for the shear design (5.8.3.) between two alternates depending on
the presence of sufficient minimum shear reinforcement. If less minimum reinforcement should
be provided, a ”crack spacing parameter” has to be specified depending on the maximum
aggregate size ag [ mm] and a maximum distance between longitudinal crack reinforcements
s :
35
se = s < 2000mm (3.147)
g + 16
The definition of the sign of the shear cut, its internal forces and torsional resistances can be
seen in the following figure. The sign of INCL applies to a rotation from the bar axis towards
the normal direction n.
T M
s z
E E A
n
N
A
M
A cut can be defined parallel to the axis by the input of YB or ZB only. In such cases the literal
S can be used for describing the location of the elastic centre.
However, the cut can also be defined by the input of multiple segments each with same identi-
fier NO and two points (YB,ZB) and (YE,ZE). Then every segment may treat a different material
and shear link properties. Intermediate segments not cutting any part of the section are not to
be defined.
• The separated part is on the right of the cut segment for CTRL FACE NEG, where right is
taken again relative to the direction A-E.
• It is recommended to separate the simpler part of a section. Holes in the separated parts
as well as any additional other parts should be avoided. They may be properly detected
only if defined in the correct sequence.
• If there is biaxial bending there are two options, see further below!
Figures 3.56, 3.57 and 3.58 show sensible and adequate definitions of shear cuts when there
is predominantly uniaxial bending (load from My + Vz ). In these cases CUT ... VYFK 0 could
also be defined explicitly in shear section definition in order to avoid numerical problems (due
to little biaxial bending).
y y
S S
z z
y y
S S
z z
y y
S S
z z
y y
S S
z z
y y
S S
z z
1. By manually defining one of the two factors VYFK 0 or VZFK 0 a shear cut is clearly
defined for only one load direction:
Like that numerical problems, which would arise from an area moment of inertia of almost
0 (see next point 2!) are prevented in a simple way. This simplification can be used
therefore with a shear cut definition according to Figure 3.56 (shear cuts in the web) and
3.58 (shear cuts in the central flange) with CUT. . . VYFK 0.
2. The shear cuts have to be defined more elaborately because the user must assure that
the area moment of inertia never becomes "0" (for this neither the lever arm in y-direction
nor in z-direction must become "0"). This is the case when separated parts are always
entirely within one quadrant, which are defined by the main axes.
FACE POS
y
z
S
CUT NO 2 CUT NO 1
>0
>0
S2 S1
CUT NO 2
CUT NO 1
>0 >0
Figure 3.59: Here the distances between the elastic centers of the separated areas perpendicular
to the main axes are always > 0 (lever arms to determine the area moments of inertia).
Figures 3.60 and 3.61 show sensible definitions of shear cuts for biaxial bending. The blue and
red shear cuts always consist of two partial cuts with the same number / same name. However,
the layer numbers are defined in such a way that they are the same in all partial shear cuts in
the bottom flange or in the left or right web. In addition, TYPE (here WRED or FLAN) can be
defined explicitly for each partial cut.
z z
CUT NO 2 CUT NO 1 CUT NO 1 CUT NO 2
S S
LAY 3 LAY 1 LAY 1 LAY 3
CUT NO 2 LAY 2
CUT NO 1 LAY 2
CUT NO 1 LAY 2
CUT NO 2 LAY 2
y y
z z
S S
CUT NO 4 CUT NO 3 CUT NO 3 CUT NO 4
LAY 3 LAY 1 LAY 1 LAY 3
CUT NO 4 LAY 2
CUT NO 3 LAY 2
CUT NO 3 LAY 2
CUT NO 4 LAY 2
In Figure 3.60 the question about necessary stirrup reinforcement in the two webs is primarily
interesting. The vertical partial cuts through the bottom flange directly next to the main vertical
axis are only required to keep each separated part within a quadrant (the results/shear flows
in these partial cuts would only be decisive if a load Mz + Vy would predominate). The two red
partial cuts are defined with CUT NO 1. . . LAY 1 (in right web, 1st partial cut) and CUT NO
1. . . LAY 2 (in bottom flange, 2nd partial cut). The two blue partial cuts are defined with CUT
NO 2. . . LAY 2 (in bottom flange, 1st partial cut) and CUT NO 2. . . LAY 3 (in left web, 2nd
partial cut). Since the relevant point for the maximum shear reinforcement is not clear (since
it depends on the load), two shear cuts through each web are examined here as an example.
Similarly for an examination of the web above the main axes further seperated parts might
have to be defined, which would have to lie in the upper quadrants.
z z
S S
CUT NO 6 LAY 2
CUT NO 6 LAY 2
CUT NO 5 LAY 2
CUT NO 5 LAY 2
CUT NO 5 LAY 2
CUT NO 5 LAY 2
CUT NO 6 LAY 2
CUT NO 6 LAY 2
Abbildung 3.54
Figure 3.61: Shear cuts in bottom flange
For shear sections 5 and 6 in Figure 3.61 the question about required reinforcement between
web and bottom flange is interesting. And here as well each separated part always remain
within a quadrant.
Hint
The possibilities shown here for the inner flange can treat a transverse shear Vy only with
the complete methods (CTRL STYP FEM as default or BEM 3).
Additionally for each edge, not only special torsional resistances can be defined for each edge,
but a proportional factor for the shear force as well. This makes it possible to describe outer
dowel joints, reduced web widths, and cuts in multiply connected sections. These factors
specify the portion of the total shear force V·S/I applied to the partial cut. The integral equation
algorithm (CTRL STYP FEM or BEM 3) will establish these factors based on the integrals of
the shear stress along the cuts automatically. For a manual definition the sign is very important
and it has to be assured that the sum of the factors is 1.0. If the center of the separated part
has the same ordinates that the center of the section, no shear can be calculated any more
based on the force method.
BMAX defines whether an equivalent hollow cross section should be used (e.g. for reinforced
concrete cross sections). If the half cut width is greater than BMAX, the found cuts are auto-
matically subdivided into two partial cuts, each of them processing the external width of the
equivalent hollow cross section. For cuts of type WEB or WRED, the automatic division into
two partial cuts is done even if no BMAX has been specified and:
• it is a simple cut with one input part on a simple cross section without any holes
• or the shear flow at both ends of the cut has an opposite direction.
The width of the equivalent hollow section is then set to half the web width.
The default of the width BMAX is derived from the defined data (material type, peripheral
reinforcements, Bredt-Area etc.) and some literals can be specified too:
EC The width is taken as the double value of the reinforcement cover, but not shorter
than A/U or the effective wall thickness for hollow sections in accordance with EN
1992-1-1 6.3.2 (1).
DIN The width is taken as the double value of the reinforcement cover, but not larger than
the effective wall thickness for hollow sections in accordance with DIN EN 1992-1-
1/NA 6.3.2, DIN FB 102.
LEON The width is taken with a sixth of the diameter of an inscribed circle in accor-
dance with F. Leonhadt und E. Mönning: Vorlesungen über Massivbau, Erster Teil:
Grundlagen zur Bemessung im Stahlbetonbau, Kapitel 9.3.3, Seite 233-235, but not
shorter than the double value of the reinforcement cover or not larger than the ef-
fective wall thickness for hollow sections in accordance with OEN EN 1992-1-1/NA.
SIA The width is taken to at least an eighth of the diameter of an inscribed circle in
accordance with SIA 262.
The position of the actually inserted torsion reinforcement is checked in this order:
• First AQUA try to determines the distance of the polygon describing the Bredt’s equivalent
cross section to the outer edge of the material.
• If this Ak polygon does not exist, the lines LRF and circles CRF of the torsionally active
longitudinal reinforcement are checked.
• If a perimeter reinforcement CURF is present, the defined value for the cover is taken into
account.
If the value BMAX is not specified explicitly, the program tries to determine for unstructured
shear cuts the effective width iteratively according to the specifications of the design code
using the selected literals EC, DIN or SIA.
Since the engineer is free to choose the thickness BMAX according to EN 1992-1-1 or the
National Annex DIN EN 1992-1-1/NA, provided that the reinforcement is inserted only in the
centre, the upper limit A/U or alternatively d/6 should be observed in any case. Analog a
lower limit of 0.4·A/U is reasonable. An average value such as min(d ; 0.75·A/U) seems ideal,
since the choice of the equivalent hollow section with presumed reinforcement arrangements
involves numerous uncertainties.
BMAX will not be preset for steel or timber sections and multiple cuts from several CUT records.
An increase in this value results in a smaller equivalent hollow cross section, thus leading to
greater shear reinforcement but smaller shear stress.
But the width can be defined directly if needed, which might be necessary for very thin open
sections.
Design of the shear reinforcement only takes place for cuts when a genuine or an equivalent
hollow cross section has been defined.
The elastic stress will be calculated completely with the selected analysis method via CTRL
STYP. For the design of concrete the resulting shear flux will be obtained by the integral of the
stresses.
For the torsion the user may prescribe an explicit value with a definition of WTM explicitly. For
the force method or CTRL STYP BEM 1, this value will be preselected via Bredt’s area as
1/(2·AK·tef). Flanges which should not be designed for torsion have to be defined with WTM
0.0.
The maximum stress in a shear cut is only informative in most cases. However, if requested, it
may be prescribed explicitly with WTD.
The effective width for the shear force can be reduced by BRED. This is appropriate, for exam-
ple, when ungrouted ducts weaken the web and the width must be reduced accordingly.
The forces N and M perpendicular to the cut are only used to describe the stress. There will
be no design for bending and normal force with the links as longitudinal reinforcement. Thus a
tensile force may lower the required link area while the total steel area becomes higher.
Attention!
In cracked state the reinforcement is used in the shear checks only when it is situated inside the
separated polygon. In particular, a partial area without reinforcement leads to a shear stress
0.0 if it is in the tensile zone! Therefore we use a minimum shear stress which is taken from
the uncracked state and a reduction factor based on the lever arms.
The individual segments of one cut are treated separately during the design. The extreme
values of reinforcement of every shear cut, which are stored separately for each rank of links,
are always per physical section, the following rules apply:
• If a cut intersects multiple webs of an articulated section, the obtained reinforcement has
to be placed in every web, The total shear link area is n times the printed value for n webs.
• If a cut in a compact section intersects the equivalent hollow section twice, the reinforce-
ment is the double of the most unfavorable segment. This assures, that the total reinforce-
ment is computed by superposition of the shear and torsion components as required by
some design codes.
• Cuts having the same link rank enter the maximum value to the resulting reinforcement.
This value has to be place then in every segment. A typical case are multiple cuts for the
same web referencing the same shear link.
• Cuts having different link ranks enter the results separately for each rank with the corre-
sponding maximum. Different link ranks should thus always be used when different rein-
forcement is to be placed in the individual parts of the cut, e.g. for multiple webs of a hollow
section.
The minimum reinforcements ASUP may be specified in cm2 / m 1021 which is referred to the
length or with or in cm2 / m2 1022 or as ratio in [-] which is referred to the length and the width
of the cut as specified in most design codes. If not defined, the default according to the design
code will be selected. If ASUP is defined as zero, it will be tried not to use any shear links (slabs
and beams of minor importance). A special minimum reinforcement enlarged by a factor of 1.6
according the Nachrechnungsrichtline resp. DIN_EN 1992-1-1 9.5b or DIN 1045-1 13.2.3 (5)
as ”articulated sections with pretensioned tensile flanges” requires the definition of Literal
’PFLA’ for ASUP. As some design tasks require to distinguish between the minimum and the
provided reinforcements, the definition should specify only the minimum reinforcements.
R Radius mm -
PHI Maximum sector angle degree 15
y , z
ye , ze
Uniform normal and shear stresses are generally assumed across the thickness of thin-walled
elements. Therefore the moment of inertia about the longitudinal axis (B-E) is zero. The shear
stresses due to torsion, however, are distributed linearly across the thickness. The thickness
is considered likewise for the maximum stress. A mixture of the cross section elements POLY
or CIRC inside a cross section is not permitted.
The transmission of shear is only possible at interconnected elements. The program does not
recognise any penetrations of thin- walled elements. In the determination of shear stresses,
elements are considered to be connected with each other when their coordinates are no further
than 0.5 mm (or snap defined by CTRL SDIV VAL2) apart.
If a radius is specified, additional straight elements are generated in order to simulate a circular
arc (180◦ maximum). The aperture angle is defined as less than PHI per segment (Default de-
fined by CTRL HMIN/HTOL). The orientation of the arc is defined by the sign of PHI respective
R. Positive values describe an arc rotating about the positive x-axis.
For a buckling design it is necessary to define the boundary conditions (free, built-in) and to
combine several plates to an integral field. To allow this we have the convention that plates will
be combined if:
The user may define the built-in face with an explicit value measured from the beginning (FIXB)
or the end (FIXE) respectively, which may be defined by the following rules:
If the thin-walled elements are cut by non-effective areas "NEFF NYZ" the resulting edges are
assumed to be free ends.
A c/t calculation via "NEFF SIG" does not change the location of the clamped edges.
Uniform normal and shear stresses are generally assumed across the thickness of thin-walled
elements. Therefore the moment of inertia about the longitudinal axis (A-E) is zero. The shear
stresses due to torsion, however, are distributed linearly across the thickness. The thickness
is considered likewise for the maximum stress. A mixture with section elements POLY or CIRC
is not permitted. A mixture with section elements PROF or WELD is permited, but PROF must
be defined as thin-walled.
The transmission of shear is only possible at interconnected elements. The program does not
recognize any penetrations of thin- walled elements. In the determination of shear stresses,
elements are considered to be connected with each other when their coordinates are no further
than 0.5 mm (or snap defined by CTRL SDIV VAL2) apart.
In the case of a variable plate thickness a variable TVAR can be defined for T in order to use it
for the interpolation with a master cross section.
If a radius is specified, additional straight elements are generated in order to simulate a circular
arc (180◦ maximum). The aperture angle is defined as less than PHI per segment (Default de-
fined by CTRL HMIN/HTOL). The orientation of the arc is defined by the sign of PHI respective
R. Positive values describe an arc rotating about the positive x-axis.
yb , zb
ye , ze
For a buckling design (c/t, section class 4) it is necessary to know the effective length c. This
value is defined for every plate by the location of the supports (FIXB, FIXE), which may be
defined by the following rules:
1001
1021
1022
1013
SC 1012 1011
1002
20
20
40 SC 420 40
20 20
If the user does not specify distances explicitly, it will be tried to guess the support condition
and to combine consecutive plates to a single field if the following conditions hold:
For the automatic determination of the locations of the clamped edge of panel sequences the
panels be aligned the same, if necessary. With different input of the start end end points thus
individual panels are changed in the direction.
This element connects thin-walled section elements PLAT and PROF shear-resistant. It is not
used for the sectional values for bending and normal force. Thus one can describe:
- Longitudinal seams of welded joints (T > 0), with contributing to the torsional moment,
the shear forces and shear deformation areas
- Dowels in composite sections (T < 0)
- Trussed walls (T < 0)
For real welds the equivalent seam thickness is to be used for T. The element is then used for
the determination of the shear stresses in welds.
For bracing walls the effective thickness T is given through the shear stiffness Sd of the frame
divided by the product of the WELD-shear modulus and the length of the shear connection in
the section plane. T = Sd /(G·L)
In the case of a variable effective thickness a variable TVAR can be defined for T in order to
use it for the interpolation with a master cross section.
The coordinates of the end points must be placed exactly within 0.5 mm (or snap defined by
CTRL SDIV VAL2) to the end points of the corresponding elements.
E E
PLAT 2 PLAT 2
DTYP Representation − S
T thin-walled
S solid cross section
TP T using positive z ordinates only
SP S using positive z ordinates only
(not available for L, T, Z and SH shapes)
SYM Symmetric options LT -
o-o duplicate symmetrically to origin
y-y duplicate symmetrically to y-axis
z-z duplicate symmetrically to z-axis
QUAD four items (all three options)
YFLP just mirror at the y-axis
ZFLP just mirror at the z-axis
REF Location of shape‘s reference point LT *
C elastic center
SC shear center
Y+Z-/Y-Z-/Z- top right/left/middle
Y+/Y-/M right/left/middle
Y+Z+/Y-Z+/Z+ lower right/left/middle
MATI Material of a filling for hollow sections − 0
Table continued on next page.
Record PROF may be entered without a preceding SECT. In that case a section with the given
shape-number NO will be generated. The plastic shear values are determined as the table
values. A faster calculation for a pre-design can be done via input ’CTRL STYP FORC’ which
is generating tabulated values only. These values are not suitable for any really calculation. In
this case do not use this feature.
If the record SECT with a section number is defined before PROF, this behaviour can be mod-
ified for a solid steel profil with CTRL PLAS.
By default profiles are oriented with their legs in the direction of the y and z coordinate axes,
so that the y axis shows the larger moment of inertia. Channels (U-sections) are oriented with
their opening to the right. Angles (L-sections) stand like the letter L.
The value Z1 - Z3 denote the title of a profile as specified in the profile tables. Not all identifiers
Z1 - Z3 always apply depending on the profile type. Only Z1 is generally used for I-profiles,
e.g. for HEB 300 with Z1 300. For L-profiles, Z1 stands for the height, Z2 for the width and
Z3 for the thickness, e.g. L 200x100x10 with Z1 200, Z2 100 and Z3 10. The identifiers to be
specified are marked with an * in the tables below.
Hint
The values Z1 - Z3 are independent of the units and cannot be changed. Deviating
geometries have to be specified with the values VD - VT2.
The actual list of shapes can be printed by an incomplete input. If nothing is selected you will
get a list of possible shapes, if you have selected a shape type you will get the list of all shapes
of that type, if you have selected only partial identifiers you will get the set of matching shapes.
Hint
For all profiles of the AISC in each case two sets of identifiers (values Z1 - Z3) are avail-
able, e.g. identifiers for in. x lb/ft and mm x kg/m
Japanese Shapes
Indian Shapes
The geometry of the shapes can be modelled only roughly. The most important data is the
weight of a single pile Z3 in kg/m and a few more dimensions. AQUA will create the locks
as equivalent solid sections such that the total area (based on a weight of 7850 kg/m3) is
achieved. Input of R1 is the half diameter of the lock (U-shapes only), R2 is the angle of the
web against the horizontal length of the wall and B2 may be used to specify the width of the
”upper” flanges. For the sheet piles the input for SYM defines:
Cold formed sections (CDS) may be defined as U, Z, C or OMEGA shapes by defining the
values of the width B and/or the grps T positive or negative:
B
B B<0
T<0 T>0
T2 > 0
B2 B2 > 0
B2
This input always defines polygons or thin-walled section elements and not any fixed cross
section properties. Therefore the analysis of thin-walled section types is performed with small
deviations from the tabulated cross section values. For torsional values with very thin shapes
with significant fillets, more severe deviations may occur. The default for DTYP is changed to
T, if beforehand a thin walled element in that cross section has been defined. If the profile
is included (oartly or complete) in concrete (also core concrete) a hole in the surrounding
concrete is created automatically unless MREF is specified with zero.
ALPH can rotate the cross section about the x axis. For the rotation of profiles with ALPH, with
the exception of L-, hollow and circular profiles, it is to be noted that the interaction conditions
of the design codes (e.g. Eurocode) refer to standard rolled sections or welded beams with the
same geometry. This means that the possibility of profile rotation means that these proofs are
no longer valid in every case. If it is only a matter of aligning the profiles correctly in the system
(on the beams), it is better to use the rotation of the local coordinate system of the beams.
When nothing is input for REF, the reference point of the shape (YM,ZM) is located at mid-
height for I- and SH-sections, left outside at mid-height for Channels, bottom left outside for
Angles and top middle for T-shapes. For centric beams REF C should be applied as with
the GUI default. The coordinates of these reference points and the angle of rotation may be
defined with reference to other points, but the shape and its size itself may of course not be
influenced by other reference points of.
The shapes can be coupled with other cross section elements, but then the type of element
must match. Thus thin walled elements (PANE, PLAT, WELD) may be combined only with
DTYP D profiles, while thickwalled elements are allowed with DTYP V profiles, overlapping
definitions with other materials will generate automatic holes. And it must be kept in mind
that the transfer of shear forces for thin walled profiles is only possible at ends, vertices, or at
the centre point (I-sections and Channels). In case of solid cross sections the edges must lie
exactly on top of each other. As can be seen in the following example, the bisected shapes
must therefore be positioned on the outer edge of the other profile’s web.
All solid I profiles and hollow sections (SH, CHS, HSS, HSSr) have preset edge properties that
are necessary for a thermal analysis.
SECT 1
PROF 101 HEB 400 ALPH 0 DTYP T
PROF 102 HEB 400 ALPH 90 DTYP T
SECT 1
PROF 101 HEA 300 ZM -145 DTYP TP ALPH 180.
PROF 102 HEA 300 ZM 145 DTYP TP
WELD 200 150 -7 150 7 14
200 -150 -7 -150 7 =
SECT 1
PLAT 101 120 -106.2 120 106.2 3.6 R -106.2
PLAT 101 120 106.2 120 -106.2 3.6 R -106.2
PLAT 102 -120 -106.2 -120 106.2 3.6 R -106.2
PLAT 102 -120 106.2 -120 -106.2 3.6 R -106.2
!Connections
WELD 300 120 110 120 106.2 6.0
300 -120 110 -120 106.2 =
300 120 -110 120 -106.2 =
300 -120 -110 -120 -106.2 =
For the combination of solid polygon items, one has to check that the polygon components
are perfectly aligned with each other along their edges, so that the shear connection can be
identified. For the definition of the following section one has to stipulate the exact coordinates
and the exact height of the second divided HEM 1000:
133
132
103134
102
152 112
123 163
104 124 164
S M
154 114
153 113
122 162
105
144106
142
143
A defined profile is mirrored three times with SYM QUAD. With the input
Additional output points for normal and shear stresses or arbitrary geometry points may be
defined with SPT. The design for fatigue is only possible with those points. A stress point
will create its own stress results in the database visible with WinGRAF along a beam or may
be addressed in MAXIMA for the superposition of stresses at that point. For most cases the
evaluation of maximum stresses might be sufficiently performed with all polygon vertices and
the intermediate points of the thin walled elements.
If the user does not select stress points, automatic points may be created controlled by CTRL
SCUT. For construction stages it should be noted, that these points will be generated for the
final section. If those points will be in a non effective part of the section for a single construction
stages, no stresses will be obtained for those points.
The normal stresses can be calculated from the X and Y coordinates using Swain’s formula.
The stress coefficients will be determined completely for thin walled sections if the stress point
is within a sectional element. If the Integral equations are active for the shear stress (CTRL
STYP BEM) you will get unit shear stress for all polygon vertices and stress points inside a
polygon. For some design tasks an effective width is needed, which will only be provided
automatically for thin walled sections.
For the fatigue design according to DIN 15018/4132 you may also specify literals for item CDYN
indicating the working conditions and notch types (B1W0 to B6K4) as well as explicit allowed
stress values SIGC and TAUC.
For the fatigue design according to DS 804 / DS 805 the notch groups (WI, WII or WIII as well
as KII, KII, KIV, KV, KVI, KVII, KIIX(!), KIX and KX) may be selected only for the design codes
DIN 18800, DIN 1045-2008 and older. SIGC and TAUC are as stress sways the double value
of those of the first row of the tables of appendix 6. Precise defaults are available for S 235
and S 355 (DS 804), smaller class values will be treated according to the formulas of DS 805.
Specifying only SIGC and TAUC will select the check of the absolute stress range as required
for solid sections.
For the stress superposition in MAXIMA, the corners and edge-mid points of the encased
rectangle will be available for all sections without the need of a SPT definition.
For a plate buckling design via a c/t ratio it is necessary to define the effective thickness TEFF
and two stress points with identical identifier describing the end points of the plate. The points
have to be defined at the mid line at the location of the built in effect. In the case of a vari-
able effective thickness a variable TVAR can be defined for TEFF in order to use it for the
interpolation with a master cross section. Following possibilities for FIX are available:
With NEFF it is possible to define parts of polygons or an FE section or a thin walled section
to be non effective for different types of forces or moments. The default NYZ, deactivating all
longitudinal stress components is the correct choice for elements which are not effective at all
(e.g. local buckling). If the non uniform distribution of longitudinal stresses is to be modeled,
the normal force should remain effective. AQB uses the total section for the normal force due
to prestress and all creep and shrinkage evaluations. For the shear stiffness the total section
is used in all cases, but the non effective parts have a constant shear stress due to the non
existing longitudinal stress.
If no information is given for the reinforcement at TYPE, it is active in the non effective part (=
+R). An input +/-R for the reinforcement does not work with th especification DIST NONE for
the distribution of the bars.
For the bending parts it has to be distinguished if the Moment My or the lever z becomes zero
in Swain’s formula. SOFiSTiK uses always the lever, thus MY and MZ are just synonyms for Z
and Y.
In general the implemented method will generate multiple deductional areas if the windows
overlap, so this must be avoided.
For the definition of the NEFF-area it is possible to define the corner points (ymin,zmin) and
(ymax,zmax) of a rectangular window or with a specification of WIDT a rectangle along a line
from (ymin,zmin) to (ymax,zmax).
The NEFF window will create additional polygon vertices for the entries and exits points. This
option may be deactivated with an extra character "‘V"’ appended to the TYPE. For NEFFs
defined within a construction stages this option is always enforced.
For sections of class 4 the procedure in the design code deactivates the non effective areas
beforehand.
357.14
-357.14
This can be facilitated by defining non effective areas in AQUA beforehand allowing a simple
stress check or even an interaction of ultimate forces in AQB.
This reduction (e.g. EN 1993-1-5) is dependent on the stresses and thus in general an iterative
process. AQUA does not use sectional forces to be used for a not yet existing section but an
explicit stress distribution The most simple case is to assume a maximum uniform compressive
stress, just by selecting SMIN without any coordinates. (EN 1993-1-5 4.3 (3))
For uni-axial bending the stress is defined by just two points on the principal axis hitting the
elastic center of the section.
But even for bi-axial bending the stress plane may be defined by just two points if we require
that the stress does not vary perpendicular to the connecting stress defining line. The exact
values may be obtained from an iterative analysis with AQB NSTR.
SMIN(YMIN,ZMIN)
Figure 3.73: Stress definition with two points on the stress line perpendicular to the neutral line
σ = E · [ε0 + (z − zs ) · ky − (y − ys ) · kz
(3.149)
The stress for two arbitrary points (not required to be within the section) symmetric to the
centroid along the stress line may be taken as
σmn (ys − · kz , zs + · ky ) = E · ε0 + · ky2 + · kz2 (3.150a)
σm (ys + · kz , zs − · ky ) = E · ε0 − · ky2 − · kz2 (3.150b)
r
= h/ ky2 + kz2 (3.150c)
Accounting for eccentricity acc. EN 1993-1-1 [Link] is done automatically for beams with a
reference axis.
These NEFF are effective for thin walled sections and solid section with defined c/t-stress
points.
CS Construction stage no − 0
KR Absolute roughness [mm] *
default according to material (MEXT)
ICE Thickness of ice cover or marine growth [mm] 1011 0
TRAF Height of additional area of wind attack due to traffic [mm] 1011 0
For the analysis of wind loading you may specify different parameters for construction stages.
This is mainly the wind attack area due to ice, traffic and building extensions. The effectiveness
of these definitions has to be aligned to the selection of reference length in record WIND.
Specifying CW enables default values of the drag coefficients based on literature for some
specific shapes. If neither WPAR nor WIND has been defined and no PROF has been defined,
a rectangular shape will be assumed. For standard circular cross sections with record SCIT or
for a circular cross section with record CIRC, the values for a circular cross section are used
without input.
The reference point of the wind forces will be by default in the geometric midpoint of the section.
Its location relative to the shear center of the section may create extra torsional moments. In
special cases you may even redefine the complete set of wind coefficients. The selected CS-
number is valid for all consecutive stages selected in SOFiLOAD with a group.
The force coefficients are needed for the calculation of wind or wave loads in dependence of the
flow direction. The values CWY, CWZ and CWT may be defined as formulas resp. variables
along a reference axis. For rectangular sections and some steel shapes these values are
known to a certain extend, For circular sections very detailed values are known depending on
the Reynolds number and roughness and available within the programs. But there are some
explicit definitions possible (e.g. cables always with Cd = 1.2). As there are no directional
dependencies in that case, a single record is sufficient and the resulting force is given by the
equation from Morrison:
1 ∂U
Fd = ρ · Cd · B · U2 + ρ · Chyd · A · (3.151)
2 ∂t
These force components act always in the flow direction. If a value for cz is specified, this force
is calculated with the absolute value of the cosine, allowing to describe simple cases also with
a single record.
For other section types more complex force coefficients may be needed, defining forces trans-
verse to the flow direction or torsional moments depending on the angle of attack of the flow.
Than CWY to CWT have to be defined depending on an angle. These values have to be stipu-
lated according to an angle sequence. Curves are extended for angles outside the range 0 to
90 or 0 to -180, if they were not defined explicitly for this range. More than a total of 99 values
will not be permitted however. For rectangular sections and standard steel shapes again de-
tailed values are known and available as default. If no values are defined a rectangular section
with the outer bounds is assumed. Defining ECHO WIND FULL will visualize the hydrodynamic
coefficients for all sections.
It is especially important to consider the exact definition of the coefficients, because the liter-
ature often uses very different coefficients. In wind engineering the generally used coordinate
system is to align the y-axis in wind direction and the z-axis vertical to the top:
• The angle for load application is the angle against the (lateral) local y-axis of the cross
section. 0 degrees correspond to an angle from left to right in negative direction of the
y-axis, +90 degrees from below to above in direction of negative z-axis. This definition cor-
responds to the normal designations, yielding forces in negative direction with the positive
signs of the force coefficients.
+ cwz
α < 0° + cwt
y
WI N D α = 0° + cwy
α > 0°
z
Figure 3.74: Angle of attack of wind
• The references WB, WH, WHB and WBH are similar to their counterparts without the W, but
the coefficients CWY and CWZ are in wind direction and transverse upwards as measured
in a wind tunnel. They are converted to the values of the schema mentioned above.
+ cl + cl
+ cm + cm
α < 0° α > 0°
y y
WI N D + cd + cd
z z
• The coefficients are be to determined from wind-tunnel tests or numerical flow simula-
tions. In order to expand the band width with regard to smaller modifications (e.g. ice
arrangement), the flow areas for lateral and vertical flows are recommended to determine
separately. For torsion both dimensions are then used. Thus for the sectional coordinate
System and REF HB it is valid:
These definitions seem to be strange, with an input to REF one can also stipulate either the
height or the width as reference for all coefficients. It has to be stated however, that this will
hamper an interpolation and some of the WPAR effects. If the values are all taken on the width,
a change of the height will not change any wind force. Many values specified in the literature
may be generalized easily if defined with the REF HB or BH. For the dervativa with HB or BH
and a rotation of WHB or WBH the value B2H is required.
For extended analysis it is also possible to define the derivatives according to Scanlan. For
every defined angle up to eight consecutive records with up to 27 data pairs may be defined.
1 z′ Bα ′ z y′ y
L
Fm = ρU2 B KH1 + KH2 + K 2 H3 α + K 2 H4 + KH5 + K 2 H6 (3.153)
2 U U B U B
1 z′ Bα ′ z y′ y
M 2 2 2 2 2
Fm = ρU B KA1 + KA2 + K A3 α + K A4 + KA5 + K A6 (3.154)
2 U U B U B
1 z′ Bα ′ z y′ y
D 2 2 2 2 2
Fm = ρU B KP1 + KP2 + K P3 α + K P4 + KP5 + K P6 (3.155)
2 U U B U B
All these values are thus referenced on the total chord width b. When importing data it has
to be taken care or the fact, that Scanlan used the half chord B and has used other indices
with the P (P1 exchanged with P5 and P4 with P6). Then he used the heave downwards, so
the coefficients H2, H3, A1 and A4 are with opposite signs. Thus some of the values might be
corrected, adjusted or shifted.
The quasi static wind coeffecients and their derivatives may be used to create asymptotic
values of the curves for static loading (K=0). The user may select this asymptotic point by a
last point onyl with VRi.
The printout of the coefficients is done with ECHO WIND FULL in the notation of the Polytec-
nico Milano which provides a better visual control. With ECHO WIND EXTR it is also provided
in the Scanlan notation.
Longitudinal reinforcements are organized in groups by Layers. For the design all elements of
a single layer will be activated in total to achieve a useful distribution for all layers [Link]
achieve this in a consistent manner, all elements of layer must have the same material and
type and should be activated in the same embedding material and construction stage. This is
achieved easily by specifying a LAY before defining the elements. The definition of a short text
(TID) may be used to clarify the location of the layer for the print. The default for TYPE is MIN
for Layer 0 and 9 for all other layers.
In general, all reinforcements are saved as single bars to allow unlimited usage of the sections
for non linear design. Several possibilities for the distribution are available via item DIST.
EVEN generates a bar in the middle of all segments, thus all areas have the same
weight.
FULL generates bars at the start and end point and between
The following scheme shows the distribution of the reinforcement bars e.g. in the case of the
partitioning of the line reinforcement in four segments:
B E
EVEN
FULL
INS,INSC
ADJA
ADJE
FULL INSC
In all other cases of DIST the value of A is preset with the maximum reinforcement bar distance
specified in the INI-file, respective a maximum distance of 400 mm. Negative definitions of A
will be taken directly as number of bars.
For the reinforcement area and especially to distinguish the behaviour with an interpolated
section three cases are distinguished, controlled primarily by the definition of a specific unit.
mm2
defines the total reinforcement area
mm2 / m
defines the grade of reinforcements
[−] defines a factor on the single bar with Diameter D
(Units [o/ o] and [o/ oo] have the same effect)
If AS is not defined, the total area of reinforcements will be obtained by the count on the
generated single bars.
The three most important variants of the input of a line reinforcement with reinforcement AS and
the maximum reinforcement ASMA are realized in the cross section editor. The corresponding
inputs for AQUA are:
LRF 'L1' YB 220. ZB 470. YE -220. ZE 470. LAY1 M1 DIST EVEN D 20 A 4[-]
Over a line length of 440 mm 4 reinforcing bars with the diameter 20 are generated respec-
tively in the middle of the four segments. The reinforcement AS which is to be considered
is used automatically with 12.56 cm2 . The maximum reinforcement ASMA can be input
with cm2 .
LRF 'L2' YB 220. ZB 470. YE -220. ZE 470. LAY M2 DIST EVEN D 20 A 100. $$
AS 28.56[cm2/m]
Over a line length of 440 mm the reinforcement bars with diameter 20 should be consid-
ered respectively in the middle of the corresponding segment using the reference length
A 100 mm. In this case the program determines four segments with a length of 110 mm.
Without an explicit input the program generates automatically 3.14 cm2 per reinforce-
ment bar for the reinforcement AS. This corresponds to a grade of reinforcement AS of
28.56 cm2 / m . The maximum reinforcement can be input with cm2 / m .
3. Input of the reference lenght A and a total reinforcement area for AS,ASMA
LRF 'L2' YB 220. ZB 470. YE -220. ZE 470. LAY M2 DIST EVEN D 20 A 100. $$
AS 12.56[cm2]
Over a line length of 440 mm with a reference length A of 100 mm a total reinforcement
area of 12.56 cm2 (corresponds to 4 diameter 20) is given for a diameter 20. In this
case the progam determines four segments of 110 mm with a reinforcing bar in the middle
of a segment with 3.14 cm2 . The maximum reinforcement can be defined also in cm2 .
For the definition of perimetric reinforcements the first variant of the input of the reinforcement
with the quantity of the reinforcing bars is not possible. In the cross section editor the variants
2 and 3 are available for the definition of the reinforcement in analog mode.
If, however, processing in the order of the layer numbers is desired, the layer numbers S0 - S9
should be used as a special case. S-layers cannot be used in combination with M- or Z-layers.
As an exception to this rule, however, a minimum (basic) reinforcement can be defined for the
lowest layer by M0.
Each layer is allowed to have only one material number. For the further processing to be
consistent, it is also useful that every one layer lies only in one material number since only
then the right deduction areas or equivalent cross section values can be determined.
Reinforcement with the same layer number is always laid in proportion to the input AS-values.
In general each layer can only be increased to the value where its first reinforcement reaches
its maximum value.
- symmetric reinforcement
Compression and tension reinforcement are defined symmetric with respect to the elas-
tic centre of the cross section and with the same layer number. Minimum (basic) rein-
forcement requires layer type M.
- non-symmetric reinforcement
Compression and tension reinforcement are assigned different layer numbers.
Crack Widths
The item AR can be used for each reinforcement element in order to introduce additional
properties for checking the crack width. D specifies already the diameter for which the crack
width must be maintained. AR defines the reference area for a single check of the crack width,
as required for instance by DIN 4227 10.2 Section 3. Here AR defines the reinforcement ratio
(AS · ƒ ctor)
μ−z = (3.156)
AR
If a FEM section with a temperature field has been selected, the temperatures TEMP will be
interpolated from that mesh. If the master section is within the same data base, a definition
of RF without coordinates will import all the reinforcements matching the selection criteria
specified wit NO (e.g. "A?" or "*" ).
The default for YB and ZB are the last defined values of YE and ZE of the previous line rein-
forcement.
b b
If a radius is input, then single reinforcement points are defined on the arc ( < 180 degree)
above the defined chord. For each segment the aperture angle is set smaller than PHI (default
from CTRL HMIN/HTOL), the orientation of the arc is defined by the sign of PHI or R.
The points are always arranged in the middle of the considered sectors. Hence, the beginning
and the end point of the arc are not reinforcement points. If an angle of 180 degrees is subdi-
vided into 30-degree segments, the single points lie at angles 15, 45, 75 as well as 105, 135
and 165 degrees.
PHI is required to shift the generated bars to a different angle. It defines the base direction
(PHI: 0 = at the z axis, -90 = at the y axis) With a definition of DPHI an arc of reinforcements is
generated from PHI-DPHI/2 to PHI+DPHI/2. A definition of DPHI 0 will generate a single point
only.
CURF can be used to define circumferential reinforcement for the last defined polygon. The
distance may be defined by two components, a constant offset H and a variant cover defined
for the material with MEXT CNOM. Any input for EXP will add the value CNOM+D/2 for all
edges. If no CNOM is defined only H is used as true cover.
For the exposition class a match is defined using the wild card rules. Thus if all edges should
be treated a value of "‘*"’ or equivalent has to be specified. Any other definition will create a
reinforcement only for those edges matching the selection criteria. Thus it is possible to define
different reinforcements for different edges. The definitions for the distribution DIST will be
applied for any connected number of edges.
The value of A specifies the maximum distance of the generated single reinforcements. How-
ever, unless DIST INS, INSC or EVEN is specified, at least one reinforcement bar is placed at
each corner. The meaning of AS, ASMA may be pre length, in total or by bar and are controlled
by the unit (eg. mm2 or cm2 / m). If definition of A as number of bars is not possible here. The
use of the line reinforcement LRF with the radius R and the angle PHI is recommended for
polygonal circles instead of the perimetric reinforcement.
Because it is useful for the design method that the centre of reinforcement coincides with the
centre of the cross section, AQUA attempts to change the reinforcement distribution with a
least-square method so that this aim is fulfilled. The sum of the perimetric reinforcement is
kept unchanged, however the reinforcement densities are increased or decreased at the single
edges. With CENT 0.0 or definition of AS as per bar the method can be deactivated.
SOFiSTiK variables, defined via STO# are globally valid. The variables LET# are only valid
locally. For more complex tasks like template section generation, it becomes necessary to
define variables valid only within a scope. TVAR allows a very general definition of those
numeric variables within a freely selectable scope. The definition of literals or the definition of
tables or functions is not possible.
TVAR without the definition of a scope saves the value for the current section. Thus the use of
this variable becomes possible for template formulas even if that variable is not defined globally.
Multiple definitions should be used for very special cases only.
• Highest Rank have variables defined along an axis during interpolation along the axis.
• Secondary rank have all variables defined in public scopes (0 to 99999), either with CADINP
and STO# or with TVAR and an explicit scope.
• Finally for all variables the default-value will become effective, which is the value which has
been defined at the time when the section has been created. If no other scope has been
used, this will be the value defined with TVAR for the section. These values are saved with
the section in the database and will be updated with every INTE command.
Variables with TVAR can also be used in formulas, e.g. when interpolating a master cross-
section along an axis. It should be noted that mixing TVAR with other variables can lead
to problems, as the evaluation and interpolation of the TVAR only takes place after the
SOFiMSHC run in a separate AQUA-INTE run. LET# variables must not be used here, as
they only work locally in a program run and are not saved in the database. It is recommended
to input all fixed values of a formula as a numerical value or with a mathematical literal, e.g.
SIN for sine, next to the TVAR.
The input of a formula in the record TVAR is not possible. Formulas have to be defined at the
corresponding cross sectional elements. Example for a variable coordinate z of an polygonal
point with TVAR H:
For the master cross section in the AQUA run with the cross section definition the coordinate
z is used with 1200 mm = value after the comma. The different heights TVAR H in mm are
defined in SOFiMSHC at GAXV along the axis. In the following interpolation, H+100 mm is
taken into account for each of the interpolated cross-sections.
There is a range of Variables in the Euro-codes allowing to adopt national factors called boxed
values. They are defined in the INI files in general. A survey of the accepted names and
subjects (DESI,COMB etc.) is given in the file [Link].
For special cases, these values may be defined with TVAR. If for example a reduction of the
elasticity modulus for the CALC curve is required, it is possible to specify factors ALF-CE and
GAM-CE :
All sections marked as to be interpolated in the database, will be established if defining INTE
ALL (or INTE 0). If NREF is specified, only sections along this structure line or axis or sec-
ondary axis are treated. The additional specification of S selects only the beam section at this
position for the interpolation. If NREF is not defined, sections are treated along all available
axes.
It is not possible to mix both variants. In the case of inconsistent definitions AQUA prints a
warning and tries to interpolate with the master cross section.
When performing this action AQUA will use free section numbers above 100 and above the
cross section template for the new sections. These sections are deleted and reassigned if this
procedure is repeated and will be printed in detail only if ECHO SDEF EXTR is given. In the
case of an interpolation with a master cross section, the numbers of the interpolated cross
sections can be influenced as follows:
E.g. if the master cross section receives the number 200, then the interpolated cross sections
start with number 201.
INTE will change the coordinates either by interpolation with a second section or by variables
along an axis. In a second step all defined reference rules of the master section will be applied.
All non geometric data is taken from the master section. Thus the interpolation between two
sections will not necessarily be identical if the sections are exchanged.
AQUA will check the generated sections against duplicates and may use a section multiple
times. A definition of INTE ALL 1 will suppress this behaviour.
The construction stage numbers of the generated sections may be changed. The definition of
a variable CSI_SELECTi (i=1...9) defines a new construction stage number for the ith construc-
tion phase. All following construction phase will be shifted by the same value if not explicitly
defined with another CSI_SELECTi definition.
Individual sections may be also directly interpolated or derived or inserted in the current sec-
tion. For the elements to be inserted there is:
• Just copy a section definition. For a cross section template currently defined variables
(TVAR) may be evaluated. This method is selected by entering NS0 only.
• You may use linear interpolation between two sections having an identical layout. For this
method two cross section numbers and an interpolation factor S (0 for NS0, 1 for NS1) have
to be stipulated.
• Last but not least you may extrude a cross section template NS0 with reference points along
a general curved axis. This is selected by entering a reference NREF of the structural line
on the axis, the common station value S for all the axes and the reference cross section
NS0.
When inserting typed standard sections, only the geometric definitions will be inserted into the
current section.
IMPO
With the record IMPO you may import materials and sections from a database to the current
project database. The imported sections may obtain a new section number and the first con-
struction stage may be shifted to a different number.
During an import used material data not defined in the target will be added automatically. A
replacement of material data is only done if a MAT definition has been made explicitly.
Imported sections will be always reanalyzed. An import will upgrade material and sectional
data as far as possible.
EXPO
With the record EXPO you may export the materials and sections in the database to an input
file for AQUA. This may be useful in special cases. If MAT or SECT is defined negative, then
the export of materials and cross sections is deactivated respectively.
If the filename is not specified the data will be appended to the most recently defined file or a
file with the name project_AQU.DAT is generated.
The units of the values will be set to the current setting of UNIE from record PAGE. The lan-
guage of the new file will be the same as the current CADINP input file.
ECHO
In case of no ECHO input all options are set to YES (ECHO MAT, SYST and LCA has NO).
The input of the option alone is therefore sufficient for increasing the value to FULL. The record
name ECHO must be entered for every record.
MAT
NO Only name of design code and materials
YES Material constants
FULL Plus stress-strain curves of the materials
EXTR Plus thermal and hydraulic material constants and Life Cycle As-
sessment
The ECHO-options may be defined anywhere in the data record, the last definition will be ap-
plied to all sections. However it is possible with a definition of ECHO SNO n1 n2 to remember
the currently active values for the print out of sections n1 to n2.
ECHO SECT allows to select individual items for FULL and EXTR, if these options are not
additional, only the selected elements will be shown.
SECT
YES Overview of cross section values only (SV,PLA 1,LAY)
FULL Plus the most important values for each cross section
(SV,PAR,ADD,PLA 2,LAY,TYPE)
EXTR Plus the individual elements of the cross section added
(SPT,PPT,CIR,PAN,WEL,CUT 2,NEFF,PRF 1,LRF)
SV sectional values of total section
PAR sectional values of partial sections
ADD additional sectional values
PLA Limit forces (1=single, 2=combined)
LAY Reinforcement summary
TYPE Typed section definitions
SPT stress points
PPT polygon vertices
CIR circles
PAN thin-walled elements
WEL welds and dowels
CUT shear cuts (1=cut, 2=points)
NEFF non effective areas
PRF Point reinforcements (1=user, 2=generated)
LRF Linear and circular reinforcements
REF Reference points (always additional)
FEM FEM-Points (always additional)
GEN all generated sectional elements (always additional)
GPRF all generated reinforcement elements (always additional)
SDEF
YES The cross sections which have been input in this calculation run
only
FULL plus the unmodified cross sections in the database
EXTR plus all interpolated sections
SYST
YES Statistics of total sum of sections and masses in the system (only
available for restart)
PICT
NO No pictures to be included
YES Nice pictures with shading
4 Description of Output
Default design code is ... The used design code is output here.
class The class input at NORM CAT is printed here
(e.g. safety class according to the design code).
Altitude above sea level
Wind zone
Snow load zone
Earthquake zone
Concrete Material
Steel Material
Timber Material
Masonry / Brickwork
Nonlinear Material
Cohesion
Tensile strength ft
Dilatancy angle
Compressive strength
Ultim. plast. strain Ultimate plastic strain
ultimate frict. angle Ultimate friction angle
ultimate cohesion
Fault dilatation
Nordic angle
Friction angle
Cohesion
Tensile strength ft
Dilatancy angle
Compressive strength
The defined stress-strain curves are output by ECHO MAT FULL. With ECHO MAT EXTR, the
standard curves are output as well:
The tangential elastic modulus is given in each case for the following range of the stress-strain
curve.
For the creep and shrinkage curves explicitly defined with the record MEXT the corresponding
curves are shown graphically.
In case a description was defined (which is the default for standard cross sections), then it is
appended after each cross section number.
After the calculation of the system a summary of the cross section types can be printed then
via a restart and ECHO SYST YES:
Following a repetition of the cross section properties from the overview, the principal moments
of inertia and the locations of the principal axes are output.
Cross section properties are also released separately for each material in the case of compos-
ite sections. For cross sections with effective width, the total cross sectional properties of the
unrestricted effective cross section are also released.
Some analysis methods, eg. a 2nd order theory of a linear calculation of steel structures (see
manual of program STAR2, record GRP, item SECT, page 3-6), require that the stiffness to
be reduced by the safety factors referenced to the reference material. In case the materials
have safety factors, this would result in totally different values for composite cross sections.
Therefore an extra table for the design sectional values is output.
For the detailed output of cross sections it is important to know, if it refers to a restart. This
is because in such a case only the newly defined cross sections get printed out. Even in the
case, however, one can output all the cross sections with ECHO SDEF FULL.
The following data are output in the table of additional cross section properties:
The table of the effective static properties and table of the design values of cross section are
printed then:
The fully plastic internal forces are output for steel or composite sections with STEE reference.
The first line contains the single forces and moments (Points A and B of the interaction dia-
gram).
Np,Rd A
Npm,Rd C
Npm,Rd
D
2
B
Mp,Rd Mm,Rd
Figure 4.1
The plastic forces are then followed by the combined values of point D marked as COMB. (For
most composite cross sections point C of the interaction curve is given by double the value of
the axial force in point D and the same moment as in Point B.)
As long as the tensile and compressive strengths of the material do not match, a second set of
values will be denoted with an inverse sign of the normal force or the moments, whereas shear
force and torsional moments are based only on the different shear strength values.
If prestressed reinforcement has been defined, the internal forces due to prestress are output.
An additional table includes the output of the design data. The values which are printed at
thet-p, thet-y, thet-z and thet-yz are with masses multiplied moments of inertia (= rotational
masses).
thet-z ρ·Iz
thet-yz ρ·Iyz
If reinforcement has been defined, then the output for each layer includes the sum of the
input steel areas, the upper and lower limits of the reinforcement, and the elastic centre of the
reinforcement.
Rectangular cross-section/T-beam
Ra Outer radius
Ri Inner radius
Rsa Radius of the outer reinforcement
Rsi Radius of the inner reinforcement
Asa Outer reinforcement
Asi Inner reinforcement
D Diameter of the reinforcement
a Distance of bars
D Nominal diameter
T Wall thickness for tube
Type Cable type
strands Number of strands
wire Number of wires per strand
W*100 Weight factor
C Sectional or fill factor
K Rupture or cable factor
ke Loss factor
Zr,k/Zr,d Failure load
Circular element
Id Designation
Mat Material number
ym Distance of the strand to the overall elastic centre
zm of the cable
R Radius of the strand
Id Section number
Type Type of section
WEB/WRED Web with/without shear region 3
FLAN/FFUL Flange with/without shear region 3
Mat Material number of the partial section
beta Parameter of friction in construction interfaces
mue Friction coefficient of construction interfaces
y-A/z-A Coordinates of the partial section
y-E/z-E
b0/bt Width of section / effective width
1/WTm/1/WTd Reciprocal torsional stress moduli
FVy/FVz Proportion factors for shear force
Ns/Ms Lateral bending internal forces
MoR Material number of the link reinforcement
Lay Layer of the link reinforcement
Asl Minimum link reinforcement
beta inclination of links w.r.t. bar axis
Id Point designation
Mat Material number
y Coordinates of the point
z
1/WT Shear stress due to torsional moment Mt=1.0
1/WVy/1/WVz Shear stress due to shear force Vy or Vz=1.0
Solid: first row = τy , second row = τz
Thin: first row = 0, second row = τ
sig-t Stress in transverse direction
W0 Unit warping
Id Point designation
Mat Material number
T Temperatur
y Coordinates of the point
z
A Effective area
t Thickness
sigma-t Prestress
1/WMy/1/WMz Inverses of the section moduli for bending
My and Mz (Swain’s formula)
1/WT Shear stress due to torsional moment Mt=1.0
1/WVy/1/WVz Shear stress due to shear force Vy or Vz=1.0
Solid: first row = τy , second row = τz
Thin: first row = 0, second row = τ
W0 Unit warping
sig/tau-dyn Permissible range of steel stresses or notch type
Id Designation
Neff Non effective flags:
N = normal force
Y = y-ordinate, bending Mz
Z = z-ordinate, bending My
Mat Material number
y-A/z-A Coordinates of beginning
y-E/z-E Coordinates of end
t Thickness
Id Designation
Mat Material number of reinforcement
y Coordinates of the reinforcement
z
T Temperature
sigma Prestress
As Base value or minimum reinforcement
As-max Maximum reinforcement
Lay Layer number
D Diameter
a Distance of bars
Ar Reference area for cracked widths
Neff Non effective flags:
N = normal force
Y = y-ordinate, bending Mz
Z = z-ordinate, bending My
Id Designation
Mat Material number of reinforcement
No Identification number
ya/za Beginning point of the linear reinforcement
ye/ze End point of the linear reinforcement
As Base value or minimum reinforcement
As-max Maximum reinforcement
Lay Layer number
D Diameter
a Distance of bars
Ar Reference area for cracked widths
Dist Distribution of bars
Circular reinforcement
Id Designation
Mat Material number of reinforcement
ym/zm Coordinates of reinforcement centre
R Radius of the circular reinforcement
As Base value or minimum reinforcement
As-max Maximum reinforcement
Lay Layer number
D Diameter
a Distance of bars
Ar Reference area for cracked widths
Dist Distribution of bars
The table of wind coefficients is output for steel profiles (record PROF) or in case of an explicit
input with WIND:
Reg Region
edge Boundary number
node-a Start node
node-b End node
M area number
MNo Material number
Conn. Hint to a connecting edge
YA, ZA Coordinates at beginning
YE, ZE Coordinates at end
Bibliography
[1] U.J. Wienecke. Zur wirklichkeitsnahen Berechnung von Stahlbeton- und Spannbeton-
stäben nach einer konsequenten Theorie [Link] unter allgemeiner Belastung. Dis-
sertation Technische Hochschule Darmstadt, 1985.
[2] C. Katz. Self-Adaptive Boundary Elements for the Shear Stress in Beams. BETECH 86,
Boundary Element Technology Conference 1986 Massachusetts Institute of Technology,
Cambridge U.S.A., 1986.
[3] D. Schade. Zur Berechnung von Querschnittswerten und Spannungsverteilungen für Tor-
sion und Profilverformungen von prismatischen Stäben mit dünnwandigen Querschnit-
ten. Z. [Link] 11 , 167-173., 1987.
[4] P.W. Rowe. “The stress-dilatancy relation for static equilibrium of an assembly of particles
in contact”. In: Proceedings of the Royal Society of London. Series A. Mathematical and
Physical Sciences 269.1339 (1962), pp. 500–527.
[5] M.M. Maksimović. Mehanika Tla. Gradjevinska knjiga, 1995.
[6] M. Wehnert. “Ein Beitrag zur drainierten und undrainierten Analyse in der Geotechnik”.
PhD thesis. Universität Stuttgart, Institut für Geotechnik, 2006.
[7] M.A. Chrisfield. Non-linear Finite Element Analysis of Solids and Structures. Volume I,
Essentials. Wiley & Sons, 1991.
[8] M.A. Chrisfield. Non-linear Finite Element Analysis of Solids and Structures. Volume II,
Advanced Topics. Wiley & Sons, 1997.
[9] O.C. Zienkiewicz and R.L. Taylor. The Finite Element Method, volume 2. McGraw Hill,
London., 1991.
[10] O.C. Zienkiewicz and I.C. Cormeau. Visco-Plasticity - Plasticity and Creep in Elastic
Solids - a Unified Numerical Solution Approach. In International Journal for Numerical
Methods in Engineering, volume 8, 1974, pp. 821–845.
[11] H.B. Kupfer and K.H. Gerstle. “Behavior of concrete under biaxial stresses.” In: Journal
of the Engineering Mechanics Division 4 (1973), pp. 853–866.
[12] Fédération internationale du Béton. fib Model Code for Concrete Structures 2010. Wil-
helm Ernst & Sohn, Berlin, 2013. ISBN: 978-3-433-03061-5.
[13] R.L. Kondner and J.S. Zelasko. A hyperbolic stress strain relation for sands. Proc. 2nd
Pan. Am. I-COSFE Brazil 1, 1963, pp. 289–394.
[14] J.M. Duncan and C.Y. Chang. Nonlinear analysis of stress and strain in soil. J. Soil Mech.
Found. Div. ASCE 96, 1970, pp. 1629–1653.
[15] O.K. Søreide. “Mixed hardening models for frictional soils”. PhD thesis. Norwegian Uni-
versity of Science and Technology, 2003.
[16] T. Schanz. Zur Modellierung des mechanischen Verhaltens von Reibungsmaterialien.
Habilitationsschrift, Institut für Geotechnik der Universität Stuttgart, 1998.
[17] T. Benz. “Small-strain Stiffness of Soils and Its Numerical Consequences”. PhD thesis.
Universität Stuttgart, Institut für Geotechnik, 2006.
[18] C.S. Desai and J.T. Christian. Numerical Methods in Geotechnical Engineerin. Chapter
2, McGraw-Hill Book Company, 1973.
[19] P. Wittke-Gattermann. Verfahren zur Berechnung von Tunnels in quellfähigem Gebirge
und Kalibrierung an einem Versuchsbauwerk. Dissertation RWTH-Aachen, Verlag Glück-
auf, 1998.
[20] W. Wittke. Grundlagen für die Bemessung und Ausführung von Tunnels in quellendem
Gebirge und ihre Anwendung beim Bau der Wendeschleife der S-Bahn Stuttgart. Veröf-
Freely defined thin-walled cross sections contain thin elements where stress variation over the thickness is negligible, resulting in properties like negligible moment of inertia about the weak axis and availability of section moduli for all points. They're suited for faster solutions but ignore local effects such as shear stress concentration at corners. Solid sections, composed of circles or polygons, have section moduli available only at distinct points and require a continua solution, which considers these local effects .
Finite element mesh-based cross sections enable detailed local stress and strain analysis by leveraging a robust computational grid that can import temperature fields and reinforcement structures, allowing high precision in simulations involving complex load interactions. In contrast, standard cross sections rely on tabulated properties, offering less precision in capturing localized effects, making them less suitable for cases where detailed stress distribution is crucial .
The von Mises viscoplastic material law is applied to metals and other frictionless materials, modeling creep effects using an elastoplastic framework with optional viscoplastic extension, typically suited for simulations involving metals. Conversely, the Drucker/Prager model is used for soil and rock, featuring friction, cohesion, and a conical yield surface, emphasizing time-dependent effects like short-term strength, making it more applicable to geotechnical engineering .
Shear deformation areas account for the distribution and magnitude of shear stresses across a cross section. They play a vital role in evaluating shearing deformation affecting overall structural integrity, especially in non-uniform or composite sections. However, their applicability is limited in capturing local stress variations in complex geometries without detailed element definition, as in freely defined FE cross sections .
Specifying reinforcement details, such as reference areas and precise distribution, ensures the reinforcement covers stress-critical zones adequately, enhancing structural stability against cracking, tension, and shear forces. This specificity is crucial in optimizing material use while meeting safety standards, ultimately supporting the durability and serviceability of reinforced elements under expected loading conditions .
Standard cross sections are advantageous when simplicity, uniformity, and availability of pre-computed properties suffice for design, such as in rolled steel shapes where torsional properties are critical. These sections provide more available data and allow skipping detailed analyses, making them preferable in projects with time constraints or when the detailed local stress effects are not critical. However, they limit flexibility compared to freely defined sections .
The mobilized friction angle measures the degree of shearing strain by evaluating the ratio between mobilized and ultimate friction angles under the Mohr-Coulomb criterion. It facilitates the assessment of the material's shearing behavior, ensuring it remains within safe limits by comparing to a calculated ultimate friction angle, a critical component for managing material stability, particularly in soil mechanics .
Linear calculations apply straightforward stress-strain relationships, which, combined with safety factors, adjust design specifications to address real-world inconsistencies and variabilities in materials. By reducing sectional stiffness in composite sections, safety factors ensure structural resilience against unforeseen loads or manufacturing defects, enhancing safety margins. This adjustment is crucial for effective design integrity, despite leading to conservative designs potentially offsetting efficiency .
The Perzyna model integrates creep parameters, allowing for time-dependent strain evaluations by applying overstress functions. It affects long-term structural performance by enabling simulations of slow deformation under constant load, crucial for assessing durability and predicting potential failure due to prolonged stress exposure. Accurately tuning these parameters ensures safety and reliability in structures subjected to sustained loads .
For open cross-sectional shapes, warping resistance and torsional moment of inertia can vary and are less effective due to an increased tendency for warping. Closed shapes, on the other hand, exhibit better torsional rigidity and warping resistance as the closed geometry inherently resists deformation and distributes stresses more uniformly, making them preferable in applications requiring high torsional stability .