ATENA Input File Format
ATENA Input File Format
Na Hrebenkach 55
150 00 Prague
Czech Republic
Phone: +420 220 610 018
E-mail: [email protected]
Web: http://www.cervenka.cz
Written by
Jan Červenka and Libor Jendele
Trademarks:
ATENA is registered trademark of Vladimir Cervenka.
Microsoft and Microsoft Windows are registered trademarks of Microsoft Corporation.
Other names may be trademarks of their respective owners.
CONTENTS
1 INTRODUCTION AND SCOPE OF THE DOCUMENT 9
2 PROGRAM EXECUTION 9
3 INPUT COMMANDS 13
3.6.14 Material Types for Simplified Nonlinear Analysis Using CCBeam Element 170
3.6.14.1 Sub-command &BEAM_MASONRY_MATERIAL 170
3.6.14.2 Sub-command &BEAM_RC_MATERIAL 174
3.6.14.3 Sub-command &BEAM_REINF_BAR_MATERIAL 178
2 PROGRAM EXECUTION
There are several methods how to execute the finite element module ATENA. The heart of
the analysis module is contained in a dynamically linked library ATENADLL.DLL. This
module can be utilized either via its COM object interface CCCoAtena or from the command
console by executing either AtenaConsole.exe or AtenaWin.exe or ATENAStudio.exe
program. The CCCoAtena is used by AtenaGUI graphical pre and postprocessor and its use
is described in a separate part of ATENA documentation. Here, the starting the analysis using
AtenaConsole, AtenaWin, and ATENAStudio executables is described. The programs are
executed as follows:
[input_file | project_file]
AtenaConsole front-end is aimed for batch analyses. Hence, it works only with input and
output files, produces no graphics and does not need any user interaction. On the other hand
AtenaWin is a windows based application. On start it creates an editable window for each of
ATENA’s window. The user can use these windows to edit content of the files, inspect
ATENA’s output during the analysis etc. Of course, similar windows can be used for editing
any other text file. It also provides graphical windows post processing and windows for 2D
plots, which are useful for example for assessing load-displacement diagram of analyzed
structure. Note that all windows in AtenaWin are updated already during the analysis.
/D path = specifies path to the working directory where input and output files will be stored.
/P = this option forces the program to request manual specification of input and output files.
/M module_name = name of main DLL library used for execution. By default, it is assumed
CCStructures. The CCStructuresCreep DLL is needed for creep analysis.
/O = specifies overwrite flag for output_file, message_file and error_file files. This means
that during execution, (or re-execution within AtenaWin) the files are created or overwritten.
By default the files are appended, i.e. output of the new analysis is added at the end of the
files.
input_file = name of a file with Atena input commands. If not specified, standard input from
keyboard is assumed.
output_file = name of a file for Atena output. If output_file doesn’t exist, it is created.
Otherwise it is appended. If output_file is not specified in the command line, then standard
output to the screen is assumed.
message_file = name of a file for Atena message output. The message file contains compact
information on Atena execution as for instance: a log of the execution start and end,
convergence performances, severe warning and error messages during execution etc. If
message_file doesn’t exits, it is created. Otherwise it is appended. If message_file is not
specified in the command line, then standard output to the screen is assumed.
error_file = name of a file for Atena error output. The error file contains full information on
Atena execution as for instance: a log of the execution start and end, convergence
performances, all warning and error messages during execution (incl. their place of
invocation) etc. If error_file doesn’t exits, it is created. Otherwise it is appended. If error_file
is not specified in the command line, then standard output to the screen is assumed.
[/translate_ids] = this option is only for internal use for debugging. Don’t use it.
[/extend_int_output_width] [/extend_real_output_width] = double number of digits used to
output integer or real numbers, respectively.
[/catch_fp_instructs] = flag to catch, (i.e. unmask) floating point exceptions during the
execution. Upon occurrence of such exception it will get caught, reported and the
execution will be terminated. By default, floating point exceptions are ignored.
[/demo-mode] = flag for trial execution. All features are available in trial mode, however,
there apply some restrictions towards size of the analyzed problem.
ATENA Input File Format 11
[/batch_execute] = option which forces AtenaWin automatically execute the given problem
without any user intervention. After the execution all output data are saved and
AtenaWin gets terminated. Use this option for batch execution.
[/execute] = option which forces AtenaWin automatically execute the given problem without
any user intervention. After the execution the AtenaWin session remains running,
thereby enabling a subsequent interactive postprocessing
[/silent] = flag that forces AtenaWin to output eventual error messages into message_file and
error_file. By default, they are output to a message box on the screen. Use this
option for batch execution.
[/num_threads=n] [/threads n] = use n threads during the execution. By default all available
processor's cores are used.
[/num_unused_threads=m]= same as the above but Atena will use number of processor's
available threads minus m. The parameter [/num_threads=n] has higher priority.
[/num_iters_per_thread=m] = chunk size for dynamic schedule, =0 for static load distribution;
default=0
[/inp] = precedes the Input File name. ATENA Studio derives the .out, .msg, and .err
filenames from the .inp filename by replacing the extension.
[input_file | project_file] = opens the Input file (.inp) or the Project file (.ccs).
If the Project file is given, it opens the existing project.
If the Input file is given there are two possibilities. If the directory, where Input file
is located, contains one or more project files, which link on this Input file, ATENA
Studio allows to choose, whether to open one of corresponding projects or create
new one. Otherwise it creates a new project.
Table 1: Environmental variables for AtenaConsole, AtenaWin 1, and ATENA Studio execution
Command
AtenaConsole 32-bit execution
%AtenaConsole% Basic AtenaConsole command, by default executes
statics module
%AtenaConsoleD% AtenaConsole execution for dynamics analysis
%AtenaConsoleC% AtenaConsole execution for creep analysis
%AtenaConsoleT% AtenaConsole execution for transport analysis
AtenaConsole 64-bit execution
%AtenaConsole64% Basic AtenaConsole 64-bit execution, by default
executes statics module
%AtenaConsoleD64% AtenaConsole 64-bit execution for dynamics analysis
1
AtenaWin program can be used for runtime visualization of the analysis progress and postprocessing. Starting
from ATENA version 5, AtenaWin program is replaced by ATENA Studio. Please, check the program
documentation of these programs for more details.
12
3 INPUT COMMANDS
C++ style where it is assumed that everything following “//” (i.e. two slash) characters
up to the end of line is considered to be a comment.
ATENA Input File Format 15
Parameter Description
NAME “task name“ Task name.
E.g.: NAME “task name“
TITLE “title“ Title of the analysis.
TITLE “title“
DIMENSION n Problem dimension. n equals 2 or 3 for two or three-
dimensional analysis. Note that setting of DIMENSION sets
also SPACE type. If DIMENSION is 2, then 2D SPACE type is
expected. Once DIMENSION type is set, it cannot be changed
elsewhere.
SPACE Set type of space approximation. It can be 2D, 3D or
AXISYMMETRIC, i.e. 2D in axis x and y symmetric with
respect to axis y. (Radius of rotation corresponds to axis x).
Note that setting of SPACE type sets also problem
DIMENSION. Once SPACE type is set, it cannot be changed
elsewhere.
Note: This command should be the first input, as it specifies dimension several entities read
later, i.e. nodal coordinates.
&BREAK:
BREAK {[AT] [MODULE module_name] ID break_id [IGNORE_HITS n_hits]} |
{[" user’s string "]}
Break Atena execution at a particular break point break_id at module module_name after
number_of_hits_to_ignore hits. The parameter module_name can be CCFEModel,
CCStructures, CCFEModelGenerate.... If no MODULE is specified, the execution terminates
at the given break point break_id at any module. If the parameter IGNORE_HITS
number_of_hits_to_ignore is not specified, the execution is terminated at the first approach of
the specified break point. Several break point ids are recognized, but break point ids 1 and 2
are probably the most important. The former one is placed at entry of a main execution
routine of each Atena's modul. Similarly, the latter one is located at the exit of that routine.
Alternatevily this command terminates the input commands stream, (i.e. no ID break_id was
input), thereby terminating the execution and optionally displays user’s string.
If the execution is run from a GUI window, (e.g. AtenaWin), a dialog is displaied before the
actual termination/break action that gives the user choice to either accept the break or ignore
it and continue the analysis. Batch analyses are broken unconditionally, see the /batch_exec
command line switch.
The commands BREAK and TERMINATE behave identically, the latter one supported only
for input compatibility reasons.
Examples:
BREAK "Joints' coordinates read"
BREAK ID 1
BREAK AT MODULE CCFEModel ID 2 IGNORE_HITS 3
&LABEL:
LABEL “string with label name"
The first command instructs Atena to ignore all subsequent input data until the second
command is found. Thereafter, the input commands are processed in the usual way. Several
&JUMP/&LABEL commands can be used in the same file providing they have unique “string
with label name". Note that &LABEL commands are ignored, unless a &JUMP command is
being processed.
ATENA Input File Format 19
Examples:
Set parameters for multihreading. num_threads sets number of threads to be used for Atena
execution. By default all available processor's cores are used. num_unused_threads is the
same same as the above but Atena will use number of processor's available threads minus
num_unused_threads. The parameter [/num_threads=n] has higher priority.
number_iters_per_thread sets chunk size for dynamic schedule. By default 0 value is
assumed, i.e. static load distribution.
Operators:
& | << >>
= <> < > <= >=
+-
* / % ||
^
!
ATENA Input File Format 21
Functions:
Abs, Exp, Sign, Sqrt, Log, Log10
Sin, Cos, Tan, ASin, ACos, ATan
Factorial. Erf, ErfInv, Atan2, Pow,
SOLVE_QUADRATIC_EQN, SOLVE_CUBIC_EQN,
Rand, SRand
Variables:
Pi, e
You can also define your own variable/parameter. e.g.
EVALUATE =cc=10" ->10
EVALUATE =cc+5 -> 15
Other:
Scientific notation supported, e.g. -1.24e-1
Error handling supported
EVALUATE =Poisson=0.3
EVALUATE =Dens=0.023
EVALUATE =FCompr=42.5
EVALUATE "FTens=0.3*pow(FCompr,0.6666)"
EVALUATE "Fract=0.000001*73.*pow(FCompr,0.18)"
EVALUATE =Young=21500.*0.9*pow(FCompr/10.,0.3333)
Parameter Description
„destination_name“ Name of the created or modified selection list.
CLEAR Clear current content of the list but doesn’t remove the
selection itself
REMOVE_DUPLICATES Remove duplicate entries in the selection
{ SOFT | FORCE } If SOFT is defined, then ATENA tolerates non-
existent source selection(s). Otherwise a input error
exception is generated.
{COMBINE | SEPARATE} “list1” Combines two or three selection lists into one list or
“list2” [“list3”] split one list into two or three selection lists. Used to
convert multi_selection lists into ordinary selection list
and vice versa.
RENAME “source_name” Rename selection “source_name” to
„destination_name“
{ FROM | AT} from_id Set interval for entity ids to be generated.
[TO to_id [BY by_id]] } They are generated for recursive formula
id1 from _ id
id n id n 1 by _ id up to id n to _ id
By default
to _ id from _ id , by _ id 1
Example:
LIST AT 1 AT 10 FROM 100 TO 150 BY 10
LIST id Entity to be added into the selection , e.g.
LIST 23 26 100
INSERT „selection_name“ Insert entities from the selection_name selection into
INCLUDE „selection_name“ the selection destination_name. Source entities, which
are already present in the selection destination_name,
are not inserted, thus avoiding entities‘ duplication.
EXCLUDE „selection_name“ Remove entities defined in the selection_name
selection from the selection destination_name. Source
entities, which are already not present in the selection
destination_name, are skipped.
CONNECT “selection_name” Connect the source selection "selection_name" with
destination selection "destination_name". This is done
in the following way: Loop from the first to the last
entry of the source selection. For each such entry loop
from the last to the first entry of the destination
selection. If the current source and destination entries
24
GENERATE { [NODES] | Data for the selection list generation. The list will
[ELEMENT] [OF] [ {GROUP | include either all nodes or all elements of the group
GROUP_FROM} group_id} <group_id.... group_id_to> from within distance x
[GROUP_TO group_id_to] ] } with respect to the point defined by the macro nodes i1.
[WITHIN] DISTANCE x FROM If <group_id_....group_id_to> is specified, elements
POINT [MACRO] [NODES] i1 are generated, otherwise nodes are generated. The
[EXECUTE] EXECUTE keyword forces to carry out the generation
immediately. Otherwise it is done prior a first step
execution.
GENERATE { [NODES] | Data for the selection list generation. The list will
[ELEMENT] [OF] [ include either all nodes or all elements of the group
{GROUP|GROUP_FROM} <group_id.... group_id_to> from within distance x
group_id} [GROUP_TO with respect to the line defined by the macro nodes i1
group_id_to]] } [WITHIN] and i2. If group_id is specified, elements are
DISTANCE x FROM LINE generated, otherwise nodes are generated. The
[MACRO] [NODES] i1 i2 EXECUTE keyword forces to carry out the generation
[EXECUTE] [INSIDE] immediately. Otherwise it is done prior a first step
execution. If the keyword INSIDE is used, the
generation is reestricted only to entities with a node
located between the macro node i1 i2.
GENERATE { [NODES] | Data for the selection list generation. The list will
[ELEMENT] [OF] [ include either all nodes or all elements of the group
{GROUP|GROUP_FROM} <group_id.... group_id_to> from within distance x
group_id} [GROUP_TO with respect to the plane defined by the macro nodes I,
group_id_to] ] } [WITHIN] i2 and i3. If group_id is specified, elements are
DISTANCE x FROM PLANE generated, otherwise nodes are generated. The
[MACRO] [NODES] i1 i2 i3 EXECUTE keyword forces to carry out the generation
[EXECUTE] [INSIDE] immediately. Otherwise it is done prior a first step
execution. If the keyword INSIDE is used, the
generation is reestricted only to entities with a node
located between the macro node i1 i2, i3.
GENERATE { [NODES] | Data for the selection list generation. The list will
[ELEMENT] [OF] include the nearest node or element of the group
[{GROUP|GROUP_FROM} <group_id.... group_id_to> with respect to the i1. If
group_id} [GROUP_TO group_id is specified, an element is included,
group_id_to]] } NEAREST otherwise a node is added. The EXECUTE keyword
[MACRO] [NODES] i1 forces to carry out the generation immediately.
[EXECUTE] Otherwise it is done prior a first step execution.
[{{IP | Generated a multiselection that includes integrated
IPS}|{ENODE|ENODES}|{GNOD points (or element nodes ) instead of global nodes. Use
E|GNODES}}] {GNODE|GNODES} to generate selection with global
nodes, where each entry must be incidented by a
element with group_id>=grouip_id_from and
group_id<=group_id_to.
{ ACTIVE | INACTIVE } GROUP Make active or inactive all elements contained in the
group_id selection list that belongs to the group group_id
26
Example:
Parameter Description
„destination_name“ Name of the created or modified selection list.
CLEAR Clear current content of the list but doesn’t remove the
selection itself
{ SOFT | FORCE } If SOFT is defined, then ATENA tolerates non-
existent source selection(s). Otherwise a input error
exception is generated.
RENAME “source_name” Rename selection “source_name” to
„destination_name“
{ FROM | AT} from_id Set interval for entity ids to be generated.
[TO to_id [BY by_id]] } They are generated for recursive formula
id1 from _ id
id n id n 1 by _ id up to id n to _ id
By default
to _ id from _ id , by _ id 1
Example:
LIST AT 1 AT 10 FROM 100 TO 150 BY 10
LIST id Entity to be added into the selection , e.g.
LIST 23 26 100
INSERT „selection_name“ Insert entities from the selection_name selection into
INCLUDE „selection_name“ the selection destination_name. Source entities, which
are already present in the selection destination_name,
are not inserted, thus avoiding entities‘ duplication.
EXCLUDE „selection_name“ Remove entities defined in the selection_name
selection from the selection destination_name. Source
entities, which are already not present in the selection
destination_name, are skipped.
REMOVE Remove the modified selection list.
ATENA Input File Format 29
Parameter Description
&ANALYSIS_TYPE Set what type of analysis is executed, i.e. static,
transient etc.
&LINEAR_SOLVER_TYPE Use direct or iterative solver (and set some vital
parameters for the iterative solver).
&CONVERGENCE_CRITERIA Convergence criteria during iteration process
within each load step.
&SOLUTION_METHOD Choose solution method for the analysis.
&ARC_LENGTH_PARAMS Set parameters for Arc Length method.
&LINE_SEARCH_PARAMS Set parameters for Line Search method.
&PREDICTOR_TYPE Set which type of predictor should be used for
building stiffness matrix, (i.e. elastic, tangential
or secant).
&UPDATE_DISPLS_STRATEGY Strategy for updating displacements during
iterations, either each iteration or each load step.
2
Not available in ATENA version 4.3.1 and older.
30
&ANALYSIS_TYPE:
{ STATIC | &TRANSIENT | &EIGENVALUES }
Parameter Description
STATIC Specify static analysis. There are no additional parameters
&TRANSIENT Set transient analysis and set some parameters for it.
&EIGENVALUES Set some parametyers for eigenvalues analysis.
&TRANSIENT:
TRANSIENT { [TIME] CURRENT x | [TIME] INCREMENT x |
TIME_INTEGRATION { {CRANK_NICHOLSON | THETA x }+ |
ADAMS_BASHFORTH } | NEWMARK BETA x | NEWMARK GAMMA x |
HUGHES ALPHA x | WILSON THETA x | DAMPING { STIFFNESS
[COEFFICIENT] x | MASS [COEFFICIENT] x | ®RESSION_DATA }+
®RESSION_DATA:
REGRESSION { MODE mode_id | OMEGA omega_val | KSI ksi_val | WEIGHT
weight_val }+ CALCULATE
Parameter Description
[TIME] CURRENT x Sets current time.
[TIME] INCREMENT x Sets time increment in steps.
TIME_INTEGRATION Set type of temporal integration scheme. If this parameter is not
input, then Newmark integration will be used.
CRANK_NICHOLSON Use linear trapezoidal integration.
THETA x parameter for trapezoidal integration. By default = 0.5.
Several other linear temporal integration may be utilized
depending on the , e.g. implicit Newton integration for = 1.,
explicit integration for = 0 etc. For good compromise between
convergence and possibility of oscillations values about =
0.85 is recommended.
ADAMS_BASHFORTH Adams – Bashforth quadratic temporal integration.
NEWMARK BETA x Defines the Newmark’s parameter.
ATENA Input File Format 33
&LINEAR_SOLVER_TYPE:
{ SOLVER { LU | DSS_LLT | DSS_LDLT | JAC | GS | ILUR | DCG | ICCG | DCGN | LUCN
| DBCG | LUBC | DCGS | LUCS | DOMN | LUOM | DGMR | LUGM | PARDISO
} | SLAP_ITERATION [LIMIT] n | SLAP_SAVED_VECTOR [LIMIT] n |
SOLVER_BLOCK_SIZE n | EXTEND_ACCURACY_FACTOR x |
PARDISO_REQUIRED_ACCURACY y | MIN_LHS_BCS_MASTER_NODE_COEFF n
}+
Parameter Description
{ SOLVER { LU | Type of solver for computing linear problem Ax=y. It can be
DSS_LLT | DSS_LDLT | either a direct skyline storage solver, (i.e. LU), or direct sparse
JAC | GS | ILUR | DCG | storage solver, (i.e. DSS_LLT, DSS_LDLT), or iterative sparse
ICCG | DCGN | LUCN | storage solver (i.e. the remaining types). Alternatively, it can be
DBCG | LUBC | DCGS | parallel direct sparse solver PARDISO from the MKL provided
LUCS | DOMN | LUOM | by Intel Visual Fortran. The skyline and sparse (SLAP) storage
DGMR | LUGM } schemes are described in the Theoretical Manual for Atena
software. The direct sparse solvers DSS_LLT and DSS_LDLT
34
In the above:
n is number of degree of freedom of the problem. nel is the number of nonzeroes in the lower
triangle of the problem matrix (including the diagonal). nl and nu is the number of nonzeroes
in the lower resp. upper triangle of the matrix (excluding the diagonal).
&CONVERGENCE_CRITERIA:
{ ABSOLUTE [ ERROR] | RELATIVE [ERROR] } | RESIDUAL ERROR x |
DISPLACEMENT ERROR x | ENERGY ERROR x | STEP_STOP_RESIDUAL
ERROR FACTOR x | STEP_STOP_DISPLACEMENT ERROR FACTOR x |
STEP_STOP_ENERGY ERROR FACTOR x | ITER_STOP_RESIDUAL
ERROR FACTOR x | ITER_STOP_DISPLACEMENT ERROR FACTOR x |
ITER_STOP_ENERGY ERROR FACTOR x | NEGLIGIBLE_RESIDUAL x |
NEGLIGIBLE _DISPLACEMENT x | NEGLIGIBLE_SIZE x |
NEGLIGIBLE_TIME_FRACTION x | ITERATION [ LIMIT ] n }+
Parameter Description
38
ABSOLUTE [ERROR] The convergence criteria values are computed using the absolute
norm that is using the maximal element of an array in its absolute
value. The error is then computed by dividing an iterative value
with the value cumulated within the whole step. Note that this
keyword can be used also in conjugation with the input
NEGLIGIBLE _SIZE n, in which case it has slightly different
meaning, see below.
RELATIVE [ERROR] The convergence criteria values are computed using the
Euclidean norm. The error is then computed by dividing an
iterative value with the value cumulated within the whole step.
Note that this keyword can be used also in conjugation with the
input NEGLIGIBLE _SIZE n, in which case it has slightly
different meaning, see below.
RESIDUAL ERROR x Convergence limit for absolute value of residual forces. Default
value is 0.01.
E.g. RESIDUAL ERROR x
DISPLACEMENT ERROR Convergence limit for absolute value of displacement increments.
x Default value is 0.01.
E.g. DISPLACEMENT ERROR x
ENERGY ERROR x Convergence limit for value of residual energy, i.e. norm of
displacement increment multiplied by norm of residual forces.
Not used in transport analysis.
Default value is 0.01.
E.g. RESIDUAL ERROR x
STEP_STOP_RESIDUAL Factors for appropriate convergence criterion value. If a
ERROR FACTOR x | convergence criterion value multiplied by the appropriate factor
STEP_STOP_DISPLACEME exceeds the related calculated analysis error, then the execution is
NT ERROR FACTOR x | immediately killed. They are two sets of factors: the first one for
STEP_STOP_ENERGY checking each iteration and the other one to be exercised at the
ERROR FACTOR x end of each step. The default value for iteration related factors is
| ITER_STOP_RESIDUAL 1000, whilst the default value for step related factors is 10.
ERROR FACTOR x | E.g.
ITER_STOP_DISPLACEME
NT ERROR FACTOR x | SET
ITER_STOP_ENERGY Absolute stop_displacement error factor 15.
ERROR FACTOR x Step_stop_displacement error factor 10.
Step_stop_residual error factor 53
Iter_stop_displacement error factor 201
Iter_stop_residual error factor 203
SET Relative
Step_stop_displacement error factor 54
Step_stop_energy error factor 55
Step_stop_residual error factor 56
Iter_stop_displacement error factor 204
Iter_stop_energy error factor 205
ATENA Input File Format 39
&SOLUTION_METHOD
{ LINEAR | NEWTON-RAPHSON | NEWTON-RAPHSON_AND_ LINE-SEARCH |
ARC-LENGTH_AND_ LINE-SEARCH | { MODIFIED_NR | FULL_NR } }+
}
Parameter Description
NEWTON-RAPHSON Use Newton Raphson nonlinear solver.
ARC-LENGTH Use Arc Length nonlinear solver.
Recommended for force loading up to
peak load or behind, can scale (reduce)
the load.
Only for static analysis, i.e., not for
40
&PREDICTOR_TYPE:
{ ELASTIC_PREDICTOR | TANGENTIAL_PREDICTOR | SECANT_PREDICTOR }
Parameter Description
ELASTIC_PREDICTOR Elastic stiffness matrix shall be used to predict displacement
increments from structural unbalanced forces. There are no
additional parameters for this command. This is option is
set by default
TANGENTIAL_PREDICTOR Tangential stiffness matrix shall be used to predict
displacement increments from structural unbalanced forces.
There are no additional parameters for this command. By
default elastic stiffness matrix is used.
&UPDATE_DISPLS_STRATEGY:
{ UPDATE_IP_EACH_STEP | UPDATE_IP_EACH_ITERATION }
Parameter Description
UPDATE_IP_EACH_STEP Specify that material points, (i.e. integration points)
ATENA Input File Format 41
&ARC_LENGTH_PARAMS:
{ &ARC_LENGTH_TYPE | &CONSTRAINT_LENGTH_CONTROL |
&LOAD_DISPLACEMENT_RATIO | &LOCATION_PARAMS
Parameter Description
&ARC_LENGTH_TYPE Set type of Arc Length method and associated
constrain.
&CONSTRAINT_LENGTH_CONTROL Set several parameters that control Arc Length
method
&LOAD_DISPLACEMENT_RATIO Control load – displacement scale for
calculating Arc Length constrain.
&LOCATION_PARAMS Set location where the Arc Length step_length
and/or Line Search energy criterion should be
calculated.
&ARC_LENGTH_TYPE:
{ CRISFIELD | NORMAL_UPDATE | CONSISTENTLY_LINEARISED |
EXPLICIT_ORTHOGONAL}
Parameter Description
CRISFIELD Crisfield variant of constant step length (including
loading space) is to be used.
NORMAL_UPDATE Updates of displacements within iteration kept normal
42
&CONSTRAINT_LENGTH_CONTROL:
{ &ARC_LENGTH_BASE_STEP_LENGTH | &ARC_LENGTH_OPTIMISATION }
Parameter Description
&ARC_LENGTH_BASE_STEP_LENGTH Set base step_length .
&ARC_LENGTH_OPTIMISATION Set the way how to optimize step_length in the
current step based on base step_length and
convergence performance in the previous step.
The base step_length is defined by
&ARC_LENGTH_BASE_STEP_LENGTH
and by default it corresponds to step_length in
the previous step.
&ARC_LENGTH_BASE_STEP_LENGTH
{ARC_LENGTH_PREVIOUS_STEP_LENGTH |
ARC_LENGTH_RESET_STEP_LENGTH | STEP_LENGTH x |
STEP_LENGTH_ONCE x | REL_STEP_LENGTH x |
REL_STEP_LENGTH_ONCE x | REL_REF_STEP_LENGTH x |
REL_REF_STEP_LENGTH_ONCE x | DLAMBDA_MIN x | DLAMBDA_MAX
x | REF_DLAMBDA_MIN x | REF_DLAMBDA_MAX x |
MIN_STEP_LENGTH x | MAX_STEP_LENGTH x |
MIN_REL_STEP_LENGTH x | MAX_REL_STEP_LENGTH x |
MIN_REL_REF_STEP_LENGTH x | MAX_REL_REF_STEP_LENGTH x }
Parameter Description
ARC_LENGTH_ For the current step use base step_length (for
PREVIOUS_STEP_LENGTH possible optimization by
&ARC_LENGTH_OPTIMISATION) from
the previous step. In case of the 1st step, it
ATENA Input File Format 43
acts according to
ARC_LENGTH_RESET_STEP_LENGTH.
ARC_LENGTH_RESET_STEP_LENGTH For the current step reset base step_length.
The actual step_length is step_length
resulting from applied load in the 1st iteration
of the current step (for = 1). It is always
calculated for the 1st step, 1st iteration.
STEP_LENGTH x Set directly required step length to x. By
default, it is initiated based on load
increment, see
ARC_LENGTH_RESET_STEP_LENGTH.
STEP_LENGTH_ONCE x Same as the above but it is appkued only
once.
REL_STEP_LENGTH x | Allows direct setting of in the next step
REL_STEP_LENGTH_ONCE x | relative to previous or reference step length.
REL_REF_STEP_LENGTH x | It can be set only "ONCE", i.e. only in the
REL_REF_STEP_LENGTH_ONCE x next subsequent step or in all subsequent
steps until a new relevant input. If x==-1, this
input is ignored. By default, all these input
valus are set to -1, i.e. they are ignored.
MIN_STEP_LENGTH x | Set minimum and/or maximum value step
MAX_STEP_LENGTH x length. If the x value is negative, this check is
ignored. By default, x=-1. This input can
overwrite DLAMBDA_MIN,
DLAMBDA_MAX
MIN_REL_STEP_LENGTH x | Set minimum and/or maximum value of
MAX_REL_STEP_LENGTH x current step length related to the step length
in the previous step. If the x value is
negative, this check is ignored. By default,
x=-1
MIN_REL_REF_STEP_LENGTH x | Set minimum and/or maximum value of
MAX_REL_REF_STEP_LENGTH x current step length related to the step length
in firrst previous Arc-Length/ Line Srearch
step. If the x value is negative, this check is
ignored. By default, x=-1
DLAMBDA_MIN x | DLAMBDA_MAX x Set minimum and/or maximum value of delta
step load increment factor. If the x value is
negative, this check is ignored. By default,
x=-1. This input can be overwritten by
MIN_STEP_LENGTH and
MAX_STEP_LENGTH
REF_DLAMBDA_MIN x | Set minimum and/or maximum value of delta
REF_DLAMBDA_MAX x step load increment factor with respepect
to reference load. If the x value is negative,
this check is ignored. By default, x=-1. This
44
&ARC_LENGTH_OPTIMISATION:
{ { ARC_LENGTH_CONSTANT |
ARC_LENGTH_VARIABLE_CONSERVATIVE_1/2 |
ARC_LENGTH_VARIABLE_CONSERVATIVE_1/4 |
ARC_LENGTH_VARIABLE_PROGRESSIVE } |
REFERENCE_NUMBER_OF_ITERATIONS }+
Parameter Description
ARC_LENGTH_CONSTANT For the current step use step_length
unchanged from the previous step.
ARC_LENGTH_VARIABLE_CONSERVATIVE_1/2 Adjusts step_length for each load
step based on the previous structural
behavior:
step_length_new=
pow(reference_number_of_iteration
/last_number_of_iteration,1/2)
ARC_LENGTH_VARIABLE_CONSERVATIVE_1/4 Adjusts step_length for each load
step based on the previous structural
behavior:
step_length_new=
pow(reference_number_of_iteration
last_number_of_iteration,1/4)
ARC_LENGTH_VARIABLE_PROGRESSIVE Adjusts step_length for each load
step based on the previous structural
behavior:
step_length_new=pow(last_number
_of_iteration/
reference_number_of_iteration,1/2)
REFERENCE_NUMBER_OF_ITERATIONS n Set optimum number of iterations
per load step to n. This value is used
in Arc Length optimization of
step_length. By default it is set to
n=5.
ATENA Input File Format 45
&LOAD_DISPLACEMENT_RATIO:
{ LOAD_DISPLACEMENT_RATIO x |
LOADING_DISPLACEMENT_RATIO_CONSTANT |
LOADING_DISPLACEMENT_SCALE_CONSTANT |
LOADING_DISPLACEMENT_BERGAN_CONSTANT }
Parameter Description
LOAD_DISPLACEMENT_ Sets the parameter ratio to x. By default, it is 0.2.
RATIO x
LOADING_DISPLACEMENT_ The SW first (i.e. in the 1st load increment) calculates
RATIO_CONSTANT scaling factor =ratio||displacements|| , where
=1 and displacements is derived from the loading
increment. The calculated is afterwards kept constant.
The ratio ||displacements||/ is called bergan
coefficient.
LOADING_DISPLACEMENT_ Adjusts (see the previous option) for each new load
SCALE_CONSTANT step as follows
=ratioberganlast
This strategy tries to keep the same impact of changes
happening in loading and geometric space.
LOADING_DISPLACEMENT_ Adjusts (see the previous option) for each new load
BERGAN_CONSTANT step as follows
=last berganold / berganlast
Subscript old stands for one before the last results. This
strategy tries to keep the same ratio of influence of
loading and geometric space.
&LOCATION_PARAMS:
LOCATION { NODE { AT n | FROM n1 [TO n2 [BY n3] ] }+ DOF { AT n | FROM n1
[ TO n2 [BY n3] ]}+ COEFF x | REMOVE}
Parameter Description
LOCATION Specifies list of domains. Each from these domains contains list
of structural DOFs and their coefficients used for calculation Arc-
length step length.
REMOVE It destroys list of domains and in the subsequent steps all
structural DOFs will be accounted for.
46
&LINE_SEARCH_PARAMS:
{ &LINE_SEARCH_ITERATION_CONTROL | &LIMIT_ETA_CONTROL |
REFERENCE_ETA x | UNBALANCED_ENERGY_LIMIT x |
&LOCATION_PARAMS }+
Parameter Description
&LINE_SEARCH_ Control several parameters for Line Search iteration process.
ITERATION_CONTROL
&LIMIT_ETA_ Set minimum and maximum value for parameters etc.
CONTROL
REFERENCE_ETA x Resets to x.
UNBALANCED_ Limit for relative work of out-of balanced forces within the
ENERGY_LIMIT x “main” iteration. When satisfied, it stops Line search internal
iteration loops. By default it is set to x=0.8 It says that Line
search has by default reduce work of out-of balanced forces
by 20%.
&LINE_SEARCH_ITERATION_CONTROL:
{ LINE_SEARCH_WITHOUT_ITERATIONS | {
LINE_SEARCH_WITH_ITERATIONS | LINE_SEARCH_ITERATION_LIMIT
n }+ }
Parameter Description
LINE_SEARCH_WITHOUT_ Do not carry internal Line search iteration loop within each
ITERATIONS “main” iteration.
ITERATION_LIMIT n
&LIMIT_ETA_CONTROL:
{ LIMIT_ETA | MINIMUM_ETA x | MAXIMUM_ETA x }+
Parameter Description
LIMIT_ETA Apply limit value for = min ... max. Only multiple of
coordinate changes are applied to the next iteration.
It is set automatically when issuing either of the commands
MINIMUM_ETA x and/or MAXIMUM_ETA x.
MINIMUM_ETA x Sets min = x. By default it is set to x=0.1
MAXIMUM_ETA x Sets max = x By default it is set to x=10.
&OPTIMIZE_PARAMS:
OPTIMIZE [BAND] WIDTH {SLOAN | GIBBS-POOLE | NONE}
Parameter Description
BAND Dummy keyword.
WIDTH Activates bandwidth minimisation and set default method to
SLOAN.
SLOAN Use Sloan’s algorithm for optimization process
GIBBS-POOLE Use Gibbs-Poole’s algorithm for optimization process
NONE Don’t optimize band-width. This is default setting.
&SERIALIZE_PARAMS:
SERIALIZE [MODEL] [STATE] { { BASICS | [AND] NODAL |
[AND] ELEMENT | ALL }+ | { DEEP | STANDARD } }+
Parameter Description
MODEL Dummy keyword
BASICS Stores just basic information about the model like number of
nodes, materials etc.
AND Dummy keyword
NODAL Stores data related to nodes of the model, (e.g. displacements)
ELEMENT Stores data related to elements of the model, (e.g. strains)
48
&FATIGUE_PARAMS:
{ FATIGUE_TASK f_task | FATIGUE_CYCLES f_cycles |
FATIGUE_MAX_FRACT_STRAIN_MULT f_mult |
FATIGUE_COD_LOAD_COEFF f_codcoeff }+
These parameters only have influence on materials that support fatigue, see the
description of the CC3DNonLinCementitious2Fatigue material.
Parameter Description
FATIGUE_TASK f_task The FATIGUE_TASK parameter determines the operation
(fatigue calculation phase) for the analysis step.
0 - nothing to do with fatigue
1 - store base stress
2 - reset FATIGUE_MAX_FRACT_STRAIN
4 - calculate fatigue damage induced by
FATIGUE_CYCLES load cycles. The calculated damage is
added to FATIGUE_MAX_FRACT_STRAIN.
8 - apply the fatigue damage stored in
FATIGUE_MAX_FRACT_STRAIN, multiplied by
FATIGUE_MAX_FRACT_STRAIN_MULT
To combine operations in one analysis step, the values are
added together (combined by binary or), e.g. storing base stress
and resetting fatigue max.fract.strain are requested by the value
3.
&CREEP_ANALYSIS_PARAMS:
{ SAMPLE_TIMES_PER_DECADE ndecl | RETARD_TIMES_PER_DECADE
ndecl_retard | STOP_TIME execution_stop_time | FIRST_DTIME
first_time_increment | {MP_METHOD | CS_METHOD} }+
Parameter Description
SAMPLE_TIMES_PER_ Number of integration times per log10 of time span. Note that
DECADE ndecl this command affects generation of integration (sample) times
by the &CREEP_STEP_DEFINITION sub-command. Hence,
the ndecl parameter must be set before the
&CREEP_STEP_DEFINITION sub-command.
This parameter defines the number of time steps, the program
will use to integrate the structural behavior. Creep or other
nonlinear effects will cause a redistribution of stresses inside the
structure. In order to properly capture such processes a
sufficiently small time steps are needed. This time spacing is
defined by the number of sample times. Its definition depends
on the type of the analyzed structure as well as on the choice of
time units. For typical reinforced concrete structures and for the
time unit being a day, it is recommended to set this parameter
to 2. This will mean that for each load interval longer then 1
day, two sub-steps will be added. For a load that is interval
longer then 10 days, 4 sub-steps will be added. For an interval
longer than 100 days, it will be 6 sub-steps.
Default value: 2.
RETARD_TIMES_PER_ Number of retardation times per log10 of time span. Note that
DECADE ndecl_retard this command affects generation of retardation times by the
&RETARDATION command and hence it must be set
beforehand. Alternatively, this value can be set directly in
&RETARDATION.
Example: If number of retardation times is set to 2, the creep
law will be approximated by two points for each time unit in the
logarithmic scale. This means two approximation points will be
50
used for the time interval between 0 - 1 day, two points for the
interval 1 - 10 days, then two points for 10 - 100 days, etc.
So the proper values will depend on the choice of time units. If
the time unit is a day, the recommended value is 1 - 2.
Default value: 1.
FIRST_DTIME Time increment for the first step of steps' serie that is
first_time_increment automatically generated after issuing a STEP definition
command, (without the option "RESUME_AT" option).
Default value: 0.1 day
Exanple SET FIRST_DTIME 0.08 ;
STOP_TIME Time at which the execution should stop [days]. This value must
execution_stop_time be input at leatest (or anywhere earlier) just before executing a
step that should by stopped by this command. If it has not been
specified, ATENA assumes STOP_TIME equal to time_end
from the &retardation times command. The inputted value of
STOP_TIME gets inserted in (automaticly generated) series of
sample times but the higher sample times are not modified.
Default value: 0 [days]
MP_METHOD | Creep analysis method. CS_METHOD uses simplified
CS_METHOD approach, in which temperature and humidity in a material point
depend only on cross sectional shape and average exterior
temperature and humidity. The MP_METHOD uses accurate
temperature and humidity at each structural material point and
therefore it need additional analysis of moisture and heat
transfer. Currently only CS_METHOD is supported.
Default value: CS_METHOD.
&DYNAMIC_ANALYSIS_PARAMS:
{ STOP_TIME execution_stop_time | LAST_TIME last_time |
{NEWMARK_METHOD | HUGHES_ALPHA_METHOD |
WILSON_THETA_METHOD | MODIFIED_ WILSON_THETA_METHOD } }+
Parameter Description
STOP_TIME execution_stop_time Time at which the execution should stop. If it is not
defined, (i.e. execution_stop_time=0), then it is
assumed execution_stop_time=last_time.
Default value: 0
LAST_TIME last_time Last time of the whole analysis.
Default value: 0
ATENA Input File Format 51
&MAX_REF_IDS:
MAX_REF_ID { { MACRO_NODES_SMART_IDS_MAP |
MACRO_ELEMENTS_SMART_IDS_MAP |
MATERIALS_SMART_IDS_MAP | LOAD_CASES_SMART_IDS_MAP |
STEPS_SMART_IDS_MAP | FUNCTIONS_SMART_IDS_MAP |
GEOMETRIES_SMART_IDS_MAP | ELEMENT_TYPES_SMART_IDS_MAP
| NODES_SMART_IDS_MAP | ELEMENT_GROUPS_SMART_IDS_MAP |
ELEMENTS_SMART_IDS_MAP [FOR] [GROUP] group_id } max_ref_id }+
Parameter Description
{MACRO_NODES_SMART Set maximum reference id for a specified data entity.
_IDS_MAP | ….. | The given value is typically used by the internal ATENA
ELEMENTS_SMART_IDS_ generator, when a request for next reference id is
MAP [FOR] [GROUP] processed. Note that if it is specified max. ref_id for
group_id } max_ref_id }+
elements, i.e. the command
ELEMENTS_SMART_IDS_MAP [FOR] [GROUP]
group_id } max_ref_id , then the group id must be id of
an already input element group. Any “forwards”
specification is not allowed here.
Default value: 50000 (for all queues).
&FORCE_UNITS:
FORCE { N | kN | MN }
&TEMPERATURE_UNITS:
TEMPERATURE { oC | oF | oK | C | F | K }
&LENGTH_UNITS:
LENGTH { MM | M | IN }
&MASS_UNITS:
MASS { KG | TON | LB}
52
&TIME_UNITS:
TIME { sec | day }
In some parts of the manual, the default values of certain material parameters are specified. If
the parameter is not specified in the input manual, the default value is used. The used default
value depends of coarse on the selected unit set. This means that the program converts the
default value to the selected unit set. The conversion is done with the help of the following
factors, whose value depends on the selected units.
Table 36: Value of factor f F for the conversion of force default values
Jednotka Faktor f F
N 1 000 000
KN 1 000
MN 1
lbf 224809.024733489
Table 37: Value of factor fl for the conversion of length default values
Jednotka Faktor fl
mm 1 000
cm 100
m 1
ATENA Input File Format 53
in 39.3700787401575
Jednotka Faktor f s
Pa 1 000 000
kPa 1000
MPa 1
psi 145.037680789469
ksi 0.145037680789469
&COORDINATES_SPEC:
COORDINATES { [ID] n [NCOORDS] ncoords [X] { x }ncoords }+
This command is used to set model joint coordinates. Each joint coordinate should be on a
separate line, e.g.
[ID] n [X] x1 x2 x3
Syntax:
54
&LOCAL:
LOCAL DOFS JOINTS { n }+
Parameter Description
LOCAL DOFS JOINTS List of nodes with local degree of freedom.
{ n }+ E.g. LOCAL DOFS JOINTS n1, n2, n3, ..... nk
Parameter Description
ID Geometry identification,
e.g. ID n
NAME User defined geometry name in quotes, also for identification.
E.g.: NAME „geometry name“
TYPE Geometry type in quotes and other geometry type dependent
parameters, see &GEOMETRY_SPECIFICATION.
&GEOMETRY_SPEC:
{ &2D_GEOMETRY_SPEC | &3D_GEOMETRY_SPEC |
&TRUSS_GEOMETRY_SPEC | &SPRING_GEOMETRY_SPEC |
&EXTERNAL_CABLE_GEOMETRY_SPEC | &BEAM_GEOMETRY_SPEC |
&LAYERED_SHELL_GEOMETRY_SPEC | &BEAM_3D_SPEC |
&BEAM_1D_SPEC }+
&2D_GEOMETRY_SPEC:
{ “2D” THICKNESS x | { REF_V1_IDS node1 node2 | REF_V1_VECTOR x y [ z] |
[DIAMETER dia]}+
Parameter Description
THICKNESS Thickness of the two-dimensional object.
E.g.: THICKNESS x
REF_V1_IDS node1 Define position of an arbitrary vector v1 used throughout
node2 definition of local coordinate system for plane 3D and 2.5D
elements. The vector is set by coordinates of finite element
ATENA Input File Format 55
&3D_GEOMETRY_SPEC:
“3D”
Parameter Description
none No parameters needed.
&TRUSS_GEOMETRY_SPEC:
“Truss” AREA x
Parameter Description
AREA Cross sectional area of a truss object.
E.g.: AREA x
&SPRING_GEOMETRY_SPEC:
“Spring” { { AREA | THICKNESS} x | { LOCAL | GLOBAL } [ SPRING]
DIRECTION { x }ncoords }2
Parameter Description
AREA | THICKNESS Cross-sectional area or spring “thickness” of a point spring or
line spring object respectively. Default = 1.0.
E.g.: AREA x
[ { LOCAL | GLOBAL } ] Spring direction in local or global coordinate system. Local
[ SPRING] DIRECTION coordinate system is applicable only for line or plane springs.
By default, global coordinate system is assumed. ncoords
coordinates defines direction vector, (ncoords equals to
problem dimension from &TASK). The direction vector
represents not only spring direction, but also its length that is
significant in case of geometrically nonlinear analyses.
56
&EXTERNAL_CABLE_GEOMETRY_SPEC:
“Cable“ { AREA x | [FRICTION] COEFFICIENT x | [FRICTION] CONSTANT x |
RADIUS x | BOND_COHESION bond_cohesion | WOBBLE_COEFFICIENT
wobble_coeff | SLIP_UNLOAD_COEFF slip_coeff_unload | FUNCTION [SLIP]
slip_function_id | FUNCTION LOCATION location_function_id |
{FIXED|PRESTRESSED} [START | END | BOTH] | PERIMETER x |
PRECISION [FACTOR] x | DAMPING [FACTOR] x | RESET_SLIPS}11
Parameter Description
AREA Cross-sectional area or spring “thickness” of a point spring or
line spring object respectively. Default = 1.0.
E.g.: AREA x
COEFFICIENT frictionlin Parameters defining friction force at a deviator.
CONSTANT frictionconst Ffrict =( ( 1.- a ) max( Fright, Fleft ) + b), where
RADIUS radius
Example:
BOND_COHESION bond_cohesion
FIXED [START | END | If specified, the starting node and/or the end node of the
BOTH] reinforcement bar is fixed with respect to the concrete, i.e. it
cannot slip. By default, if FIXED command is not used, it can
slip everywhere.
PRESTRESSED [START | Similar info as that above. PRRESTRESSED START means
END | BOTH] the same as FIXED LEFT etc.
FUNCTION [SLIP] Id of a function, by which all the coefficients are multiplied,
slip_function_id i.e. frictionlin , frictionconst. If not specified, no multiplication
occurs. The functional argument is current (total) deviator slip.
FUNCTION LOCATION Id of a function, by which all the coefficients are multiplied,
location_function_id i.e. frictionlin , frictionconst. If not specified, no multiplication
occurs. The functional argument is distance between the 1st
node and the current node, for which the slip parameters are
calculated.
For cables, the two current friction values are calculated
frictionconst_current = frictionconst fs(s) fd(dist),
and
58
&BEAM_GEOMETRY_SPEC:
“Beam” { AREA x | [MOMENT] INERTIA_Y x | [MOMENT] INERTIA_Z x |
[MOMENT] POLAR x | [MOMENT] TORGUE x | [MOMENT] SHEAR_Y x |
[MOMENT] SHEAR_Z x | [WINKLER] [COEFFICIENT] C_1_X x | [WINKLER]
[COEFFICIENT] C_1_Y x | [WINKLER] [COEFFICIENT] C_1_Z x | [PASTERNAK]
[COEFFICIENT] C_2_X x | [PASTERNAK] [COEFFICIENT] C_2_Y x | [PASTERNAK]
[COEFFICIENT] C_2_Z x | [LOCAL] [Z] [AXIS] DIR_X x | DIR_Y x | DIR_Z x |
[{SIZE_LOCAL_Y | WIDTH} x ] | [{SIZE_LOCAL_Z}|{HEIGHT} x] |
[{KIRCHHOFF}|{MINDLIN}|{TIMOSHENKO}|{TIMOSHENKO_CSF}] |
[REDUCE_TM_STIFF ] | [REDUCE_MT_STIFF] | [RO_N x] | [EFF_WIDTH_FACTOR x]
| [EFF_HEIGHT_FACTOR x] | [UPDATE_BEAM_DIR] |
[MAX_NUMBER_OF_ITERATIONS_FOR_REDUCE_FORCES n] |
[MAX_ERROR_FOR_REDUCE_FORCES x] | S_MIN s_min S_MAX s_max T_MIN t_min
T_MAX t_max| [BARS NUMBER n { MATERIAL n BAR_AREA x BAR_LOCAL_Y x
BAR_LOCAL_Z x}n] }
Parameter Description
AREA Cross-sectional area of a beam object. Default = 1.0.
E.g.: AREA x
INERTIA_Y Cross-sectional inertia moment of a beam object with respect
to local Y-axis.
E.g.: INERTIA_Y x
INERTIA_Z Cross-sectional inertia moment of a beam object with respect
to local Z-axis.
E.g.: INERTIA_Z x
POLAR Cross-sectional polar moment of a beam object with respect to
local X-axis.
E.g.: POLAR x
TORGUE Cross-sectional moment of a beam object in torque.
E.g.: TORGUE x
SHEAR_Y Cross-sectional shear moment of a beam object with respect to
local Y-axis.
E.g.: SHEAR_Y x
SHEAR_Z Cross-sectional shear moment of a beam object with respect to
local Y-axis.
E.g.: SHEAR_Z x
C_1_X Winkler (or C1 Pasternak) coefficient with respect to local X-
axis.
E.g.: C_1_X x
60
default it is 1.
E.g. RO_N 0.5
EFF_WIDTH_FACTOR x Coefficient for buckling widtf of comperessed columns’ cross
section. By default it is 1.
E.g. EFF_WIDTH_FACTOR 0.5
EFF_HEIGHT_FACTOR Coefficient for buckling height of comperessed columns’ cross
x section. By default it is 1.
E.g. EFF_HEIGHT_FACTOR 0.5
UPDATE_BEAM_DIR Flag for updating beam’s position already during iterations
with a load step. By default it is updated only at e ach step.
MAX_NUMBER_OF_ITE Maximum number of iterations for establishing force/moment
RATIONS_FOR_REDUC equilibrium. Such procedure is needed typically after any of
E_FORCES n beam’s nodal forces/moments have been reduced due to
material nonlinearity. By default 30 iterations are allowed.
MAX_ERROR_FOR_RE Acceptable relative error for the iteration process described
DUCE_FORCES x above. By default the value 0.01 is used.
S_MIN s_min S_MAX Definition of reinforcement bars in the cross section. First
s_max T_MIN t_min number of bars is read and then for each bar its material, area
T_MAX t_max BARS and coordinates are inputed. Note that all the values are
NUMBER n { specified in isoparametric coordinate system, i.e. in
MATERIAL n coordinates <s_min..s_max>, (for direction of the cross
BAR_AREA x sectional width) and <t_min…t_max>,m(for height). By
BAR_LOCAL_Y x default, these intervals are set to <-1..1>, which corresponds
BAR_LOCAL_Z x}n] to isoparametric coordtinates. If the intervals <0..width>,
<0..height> are use. the the bar areas and coordinates are input
in real coordina system with origin in the left bottom corner.
&LAYERED_SHELL_GEOMETRY_SPEC:
“LayeredShell” { DETECT_DEPTH {DETECT_VECTOR x1 x2 x3 } | { REF_V1_IDS
node1 node2 | REF_V1_VECTOR x y z } | INTERFACE interface_nodes_list | [ SOLID |
REINFORCEMENT ] LAYER n [ { [MATERIAL mat_id ] [THICKNESS thick ]
[POSITION pos ] DIAMETER dia]} | {SAME_AS layer_id }| REF_THICK x |
IGNORE_REF_THICKNESS | { REDUCE_TAU_XZ_YZ | REDUCE_TAU_XY |
FULL_TAU} }+ THICKNESS_EQN " eqn_string"
[REDUCE_TAU_XY] Reduce shears by the factor 0.85.
[REDUCE_TAU_XZ]
[FULL_TAU]
Parameter Description
SOLID | The data that follow specify a solid, (i.e. concrete) or
REINFORCEMENT reinforcement, (i.e. steel) layer.
62
SAME_AS layer_id Specifies that the layer n has the same properties as a
previously defined layer layer_id
DETECT_DEPTH Detect depth of shell elements and reorder element’s
{DETECT_VECTOR x1 x2 incidences. If DETECT_VECTOR is not specified, the depth
x3 } is chosen to comply with the smallest dimension of the
element. Otherwise it is chosen to have the smallest angle
with the given vector {x1, x2, x3}.
REF_V1_IDS node1 node2 Define position of an arbitrary vector v1 used throughout
definition of a shell local coordinated system, see the Atena
Theory Manual. The vector is set by coordinates of finite
element nodes node1 (tail) and node2 (head). By default, this
input need not be specified. In such a case, Atena kernel will
construct v3 using the default definition from the Atena
Theoretical Manual .
REF_V1_VECTOR x y z Same as tha above, but the arbitrary vector is input directly.
REF_THICK x Reference thickness used to transform normalized layer
coordinates to real coordinates. By default, this value is not
specified and in this case actual shell thicknesses at
integration points are used instead. This input is particularly
useful, if a reinforcement layer is placed at constant distance
from the shell bottom or top surface, whereby the shell real
thickness is variable.
IGNORE_REF_THICKNE This flag can be input for each reinforcement layer. If it is
SS specified, the consequences implied by inputting above
REF_THICK x are for this particular reinforcement layer
ignored.
INTERFACE Name of list that includes nodal ids, where all 6 shell DOFs
interface_nodes_list should be retained. Use this feature to connect shell elements
with other solid elements, e.g. bricks.
{ REDUCE_TAU_XZ_YZ | Reduce the specified shear(s) by 1/6 of its original value to
REDUCE_TAU_XY | compensate for constant shear strain thru cross section. By
FULL_TAU}+ default, no reduction is carried out, (recommended). .
(Ahmad elements use always full shear strains without any
reduction).
THICKNESS_EQN " String containing equation to caculate shell's thickness. It can
eqn_string conation placeholders "x", "y", "z" that are replaces by actual
shell coordinates.
Example:
THICKNESS_EQN "0.2+x*0.001+y*0.002"
64
&BEAM_3D_GEOMETRY_SPEC:
“Beam3D” [DETECT_AXIS [DETECT_AXIS_VECTOR x1 x2 x3 ]]
[DETECT_HEIGHT [DETECT_HEIGHT_VECTOR x1 x2 x3 ]]
[NUMBER_OF_IPS_IN_R n] [SOLID] HEIGHTS NUMBER n VALUES val1,
val2 .. val_n WIDTHS NUMBER n VALUES val1, val2 .. val_n DOMAINS
NUMBER n MATERIAL {n|0} QUAD_IDS {FROM n [TO n [BY [n]]] | AT n |
LIST i1,i2…} [[REINFORCEMENT] BARS NUMBER n {MATERIAL mat_id
ST_AREA a S_COORD s T_COORD t [IGNORE_REF_HEIGHT]
[IGNORE_REF_WIDTH] ]n [REDUCE_TAU_XY] [REDUCE_TAU_XZ]
[FULL_TAU] [REF_HEIGHT x ] [REF_WIDTH x]
Parameter Description
SOLID | The data that follow specify a solid, (i.e. concrete) or
REINFORCEMENT reinforcement, (i.e. steel) layer.
HEIGHTS NUMBER n Total number of solid heights, i.e. number of rows of the s, t
VALUES val1, val2 .. val_n raster. It is followed of actual height values. Isoparametric
coordinates are used. Otherwise, the input heights are scaled
so that their sum will equal to 2.
WIDTHS NUMBER n Ditto for widths.
VALUES val1, val2 .. val_n
DOMAINS NUMBER n Definition of material domains. The quad_ids are counted
MATERIAL {n|0} rowvise starting from the bottom left corner. If material_id is
QUAD_IDS {FROM n [TO zero, a hole is assumed.
n [BY [n]]] | AT n | LIST
i1,i2…}
[REINFORCEMENT] Number of reinforcement “bars”, i.e. quads, where
BARS NUMBER n reinforcement is assumed
MATERIAL mat_id For n bars specify its material id, area and position via s, t
ST_AREA a coordinates. Isoparametric coordinates are used, otherwise
S_COORD s the scaling factors are applied. The factors are those used for
T_COORD t scaling solid heights and widths.
[IGNORE_REF_HEI
GHT]
[IGNORE_REF_WID
TH]
DETECT_AXIS Detect axis of beam elements and reorder element’s
{DETECT_AXIS_VECTO incidences. If DETECT_AXIS_VECTOR is not specified,
R x1 x2 x3 } the axial direction is chosen to comply with the biggest
dimension of the element. Otherwise it is chosen to have the
smallest angle with the given vector {x1, x2, x3}.
DETECT_HEIGHT Detect height of beam elements and reorder element’s
{DETECT_HEIGHT_VEC incidences. If DETECT_HEIGHT_VECTOR is not
TOR x1 x2 x3 } specified, direction of the beam’s height is chosen to comply
with the bigger dimension of the element’s cross section.
ATENA Input File Format 65
&BEAM_1D_GEOMETRY_SPEC:
“Beam1D” CS_WIDTH_EQN "eqn_expression" CS_HEIGHT_EQN
"eqn_expression" VT_X_EQN "eqn_expression" VT_Y_EQN "eqn_expression"
VT_ZEQN "eqn_expression" [NUMBER_OF_IPS_IN_R n] [SOLID]
HEIGHTS NUMBER n VALUES val1, val2 .. val_n WIDTHS NUMBER n
VALUES val1, val2 .. val_n DOMAINS NUMBER n MATERIAL {n|0}
QUAD_IDS {FROM n [TO n [BY [n]]] | AT n | LIST i1,i2…}
[[REINFORCEMENT] BARS NUMBER n {MATERIAL mat_id ST_AREA a
S_COORD s T_COORD t [REF_HEIGHT x ] [REF_WIDTH x] ]n ]
[REDUCE_TAU_XY] [REDUCE_TAU_XZ] [FULL_TAU]
Parameter Description
SOLID | The data that follow specify a solid, (i.e. concrete) or
REINFORCEMENT reinforcement, (i.e. steel) layer.
HEIGHTS NUMBER n Total number of solid heights, i.e. number of rows of the s, t
VALUES val1, val2 .. val_n raster. It is followed of actual height values. Isoparametric
coordinates are used. Otherwise, the input heights are scaled
so that their sum will equal to 2.
WIDTHS NUMBER n Ditto for widths.
VALUES val1, val2 .. val_n
DOMAINS NUMBER n Definition of material domains. The quad_ids are counted
MATERIAL {n|0} rowvise starting from the bottom left corner. If material_id is
QUAD_IDS {FROM n [TO zero, a hole is assumed.
n [BY [n]]] | AT n | LIST
66
i1,i2…}
[REINFORCEMENT] Number of reinforcement “bars”, i.e. quads, where
BARS NUMBER n reinforcement is assumed
MATERIAL mat_id For n bars specify its material id, area and position via s, t
ST_AREA a coordinates. Isoparametric coordinates are used, otherwise
S_COORD s the scaling factors are applied. The factors are those used for
T_COORD t scaling solid heights and widths.
[REF_HEIGHT x ]
[REF_WIDTH x]
CS_WIDTH_EQN Width and height of beam's cross section. Both are given in
"eqn_expression" terms of algebraic expression f(x,y,z), in which the
CS_HEIGHT_EQN parameters x,y,z, (i.e. coordinates) are substituted
"eqn_expression" automatically based on location a beam using this geometry.
Example:
CS_WIDTH_EQN "0.5+0.1*x" CS_HEIGHT_EQN "0.1"
VT_X_EQN Algebraic expressions for x,y,z coordinates of the vector vt.
"eqn_expression" Theey are input in similar way to the above cross section's
VT_Y_EQN dimensions.
"eqn_expression" Example:
VT_ZEQN VT_X_EQN "0" VT_Y_EQN "0" VT_Z_EQN "0.3"
"eqn_expression"
[NUMBER_OF_IPS_IN_R Number of integration points in beam’s longitudinal axis. By
n default 2 IPs are used, however especially in case of heavy
material nonlinearity, more IPs may yield more accurate
results, as the beam can better locate a material failure. Max.
value is 6.
[REDUCE_TAU_XY] Reduce shears by the factor 0.85.
[REDUCE_TAU_XZ]
[FULL_TAU]
[REF_HEIGHT x ] Reference height and width of the beam's cross section used
[REF_WIDTH x] to transform normalized cross sectional local s,t coordinates
to real coordinates, (unless they are ignored). By default,
theses values are not specified and in this case the actual
beam's dimensions at integration points are used. This input
is particularly useful, if a reinforcement cell is placed at
constant distance from surface of the beam, whereby the
beam's cross section dimensions are variable.
Parameter Description
&ELEMENT_GROUP This sub-command begins the definition of a new element
group. This command should be followed by the definition
of element connectivity by using the sub-command
ELEMENT INCIDENCES
&ELEMENT_ TYPE Define a new element type. This element type is later
referred to by the sub-command &ELEMENT_GROUP to
specify an element type/formulation for an element group.
&ELEMENT_INCIDENCES This sub-command should follow the command
&ELEMENT_GROUP. It is used to define element
connectivities.
&ELEMENT_MATERIALS This sub-command should follow the command
&ELEMENT_GROUP. It sets material types individually
for each material point of the element. If not specified,
default material type from &ELEMENT_GROUP is used.
&ELEMENT_GROUP:
GROUP { ID n [NAME = “element group name”] TYPE n MATERIAL n
GEOMETRY n | DELETE | ACTIVE | INACTIVE | ASSOC_LC_ID lc_id}+
Parameter Description
ID n Element group identification
E.g.: ID n
NAME “element group Element group name in quotes, also for identification
name” E.g. NAME “element group name“
TYPE n Element type identification.
E.g.: TYPE n
MATERIAL n Identification number of material to be used for this element
group.
E.g.: MATERIAL n
GEOMETRY n Identification number of geometry to be used for this element
group.
E.g.: GEOMETRY n
DELETE Resets content of the element group to default, i.e. removes its
all-previous input data.
ACTIVE | INACTIVE Marks all elements within the group as active or inactive. Active
elements are included in the analysis, whereas inactive elements
are ignored.
68
ASSOC_LC_ID lc_id Associated load case id. This input is generated automatically,
however in some cases it allows to manually specify load case
id associated with this group. For example, if discrete
reinforcement bars are input manually, i.e. not generated, the
lc_id says, which load case is used to bind the bar with the
surrounding solids.
&ELEMENT_TYPE:
TYPE { ID n | NAME “element type name” | { LINEAR | NONLINEAR |
IGNORE_NEGATIVE_JACOBIAN |
IGNORE_ELEMENT_TYPE_EXCEPTION | SEMINONLINEAR } TYPE
“element_type” | GAMMA_REF x | GAMMA_COEFF x |
PREPARE_CALCULATION |
[ "DEFAULT_PROCESSING" | "INITIAL_STRAIN_ONLY_INTO_SOLID" |
"INITIAL_STRESS_ONLY_INTO_SOLID" |
"INITIAL_STRAIN_ONLY_INTO_REINF" |
"INITIAL_STRESS_ONLY_INTO_REINF" ] }+
Parameter Description
ID n Element type identification
E.g.: ID n
NAME Element group name in quotes, also for identification
“element type name” E.g.: NAME “My_CCIsoBrick“
LINEAR Forces to ignore all terms due to geometrical non-linearity. Material
linearity still may exist.
NONLINEAR Forces to account for all terms due to geometrical non-linearity. This
is the default setting.
SEMINONLINEAR Linear in the 1st iteration, nonlinear in the next iterations. This option
is sometimes advantageous, if the structure is loaded by
deformations.
IGNORE_NEGATI By default, if an element at an IP has negative jacobian, (i.e. there is
VE_JACOBIAN something wrong with its geometry), then the appropriate exception
is thrown and the execution is terminated. If this keyword is
specified, the negative value of the jacobian is ignored and the
execution continues (, although with disputable results).
IGNORE_ELEMEN By default, if an exception occures within element calculation, then
T_TYPE_EXCEPTI the appropriate exception is thrown and the execution is terminated.
ON If this keyword is specified, execution of the offending element is
skipped and remaining elements get proceeded.
TYPE Element type in quotes.
“element_type” E.g.: TYPE “element_type”, where “element_type” adopts form
name<xx_x…>, where x and _ characters in the <> brackets indicate
ATENA Input File Format 69
CCCircumferentialTruss
CCPlane-Spring
CCIso-Triangle
CCLine-Spring
CCQ10Sbeta
CCIso-Truss
CCIso-Brick
CCIso-Quad
CCIso-Gap
CCSpring
CCQ10
CC1DElastIsotropic X X X X X
(*)
CCPlaneStressElastIsotropic (*) X X X X X X X X X
CCPlaneStrainElastIsotropic (*) X X X X X X X X X
CC3DElastIsotropic X X X X X X X X X X
(*)
CCASymElastIsotropic X X X X X X X X
(*)
CC3DBiLinearSteelVonMises (*) X X X X X X X X X
CC3DCementitious X X X X X X X X X
CC3DNonLinCementitious X X X X X X X X X
CC3DNonLinCementitious2 (*) X X X X X X X X X
CC3DNonLinCementitious2User (*) X X X X X X X X X
CC3DNonLinCementitious2Variable X X X X X X X X X
CCSBETAMaterial X X X X X X X X X
CC2DInterface X
ATENA Input File Format 73
CC3DInterface X
CCReinforcement X X X X X
CCCyclingReinforcement X X X X X
CCSmearedReinf X X X X X
CCCircumferentialSmearedReinf X X X X
CCSpringMaterial X X X
CC3DDruckerPragerPlasticity X X X X X X X X X
CCMaterialWithVariableProperties X X X X X X X X X
CCMaterialWithTempDepProperties X X X X X X X X X
CCMaterialWithRandomFields X X X X X X X X X
CCCombinedMaterial X X X X X X X X X
Table 56 : Element Type and Material Compatibility, (beam and shell elements)
CCIsoShellTriangle
CCIsoShellWedge
CCIsoBeamBrick
CCIsoShellBrick
CCIsoShellQuad
CCIsoBeamBar
CCBeamNL
CCBeam
Ahmad
CC1DElastIsotropic X3 X3 X3 X3 X3 X3 X3 X3
(*)
CCPlaneStressElastIsotropic (*)
CCPlaneStrainElastIsotropic (*)
CC3DElastIsotropic X X X X X X X X
(*)
CCASymElastIsotropic
(*)
CC3DBiLinearSteelVonMises (*)
CC3DCementitious X X X X X X X X
CC3DNonLinCementitious X X X X X X X X
CC3DNonLinCementitious2 (*) X X X X X X X X
CC3DNonLinCementitious2User (*) X X X X X X X X
CC3DNonLinCementitious2Variable X X X X X X X X
CCSBETAMaterial
CC2DInterface
CC3DInterface
CCReinforcement X3 X3 X3 X3 X3 X3 X3 X3
CCCyclingReinforcement X3 X3 X3 X3 X3 X3 X3 X3
CCSmearedReinf X3 X3 X3 X3 X3 X3 X3 X3
CCCircumferentialSmearedReinf X3 X3 X3 X3 X3 X3 X3 X3
3
For reinforcement.
74
CCSpringMaterial
CC3DDruckerPragerPlasticity X X X X X X X X
CCMaterialWithVariableProperties X X X X X X X X
CCMaterialWithTempDepProperties X X X X X X X X
CCMaterialWithRandomFields X X X X X X X X
CCCombinedMaterial X X X X X X X X
CCBeamMasonryMaterial X
CCBeamRCMaterial X
Table 57 : Beam and shell elements and their element idealisation, material idealisation and
geometry type
CCIsoShellTriangle
CCIsoShellWedge
CCIsoBeamBrick
CCIsoShellBrick
CCIsoShellQuad
CCIsoBeamBar
CCBeamNL
CCBeam
Ahmad
Element
geometry type
LAYEREDSHELL
LAYEREDSHELL
LAYEREDSHELL
LAYEREDSHELL
LAYEREDSHELL
BEAM_ 3D
BEAM_3D
BEAM_3D
BEAM
Element
idealisation 4 LAYERED_SHELL_2D
LAYERED_SHELL_2D
LAYERED_SHELL
BEAM_NL_1D
BEAM_NL_3D
BEAM_3D
BRICK
BRICK
BRICK
Material
BEAM_3D
BEAM_3D
THREE_D
THREE_D
THREE_D
idealisation
ONE_D
SHELL
SHELL
SHELL
4
Defined by a finite element that is used.
ATENA Input File Format 75
Element
SHAPE_SHELL_WEDGE_3D
SHAPE_QUADRILATERAL
SHAPE_SHELL_BRICK_3D
shape4
SHAPE_TRIANGLE
SHAPE_BRICK
SHAPE_BRICK
SHAPE_BRIC
SHAPE_BAR
SHAPE_BAR
The above tables apply in full for static and dynamic analysis. As far as creep analysis is
concerned, it uses time independent and time dependent materials:
Time independent material (as indicated by the name) does not change its behaviour with age.
Such a material is e.g. used for reinforcement. Any material from the above table can be used
as time independent material for creep analysis.
On the other hand, concrete is known to change its properties with time and therefore (within
a creep analysis) it must be modelled by a time dependent material &CREEP_MATERIAL.
Only materials marked with "*" (from the above table) can be used as the parameter
“short_term_material_type", (refering to the definition of &CREEP_MATERIAL).
Transport analysis uses completely different element types and element material models.
They are described in Section 3.14. Any transport element type can be used in conjugation
with any transport material model.
&ELEMENT_INCIDENCES:
[NNODES num_nodes ]
id1 { n }number_nodes_1
id2 { n }number_nodes_2
…
idm { n }number_nodes_m
Parameter Description
[NNODES num_nodes ] Optional number of element incidences. If not defined,
num_nodes is derived from the element’s element type.
id Element id.
E.g.: n
{ n }number_nodes Element incidences, i.e. ids of nodes incidenting with the
element. number_nodes integer numbers is expected, where
number number_nodes is number of element nodes for the
particular element type
E.g.: n1 n2 [n3 ] .... [n number_nodes ]
76
Note:
This command has to follow the command ELEMENT GROUP.
Each element incidences data must be input on a separate line.
&ELEMENT_MATERIALS:
id1 { n }number_of_material_points
id2 { n }number_of_material_points
…
idm { n }number_of_material_points
Parameter Description
id Element id.
E.g.: n
{ n }number_of_material_points Material type at element’s material point. By default, a positive
integer value is expected for each material point of the element.
If the input value n is zero, it indicates that this and all
remaining material points use the default material type.
If the input value n is negative, it indicates that this and all
remaining material points are of type (-n).
If the element uses the same material types in all its material
points, the &ELEMENT_MATERIALS command can be
omitted and a default material type specified in
&ELEMENT_GROUP is adopted.
E.g.: 10 20 30 40
E.g. 10 -20
Note:
This command has to follow the command ELEMENT GROUP.
Each element material type’s data must be input on a separate line.
Syntax:
&NODAL_IMPERFECTIONS :
ATENA Input File Format 77
NODAL_IMPERFECTIONS [SETTINGS] {
&MANUAL_IMPEREFECTIONS_ENTRY |
&GENERATED_IMPEREFECTIONS_ENTRY }
&MANUAL_IMPEREFECTIONS_ENTRY:
{ NODE n { TOTAL | INCREMENT | INCREMENTAL} {VALUE | VALUES } val_x
val_y [val_z] }
Sub-Command Description
NODE n Set initial conditions for node n.
{VALUE | VALUES} val_x Specify initial nodal imperfections in direction of global
val_y [val_z] coordinates. 3D problems need 3 values, 2D problems only
two values..
{TOTAL | INCREMENT | Set input for total or incremental (with respect to the reference
INCREMENTAL } coordinates) values of the imperfect structural geometry.
&GENERATED_IMPEREFECTIONS_ENTRY:
Sub-Command Description
SELECTION Name of selection, for which the generation is requested.
"selection_name"
{GENERATE Keyword for entities to be generated. The values in global
GENERATE_VEL } structural directions are generated as linear combination:
CONST const_vector
COEFF_X coeff_x_vector | valuex const (1) x coeff x (1) y coeff y (1) z coeff z (1)
COEFF_Y coeff_y_vector |
COEFF_Z coeff_z_vecor valuey const (2) x coeff x (2) y coeff y (2) z coeff z (2)
valuez const (3) x coeff x (3) y coeff y (3) z coeff z (3)
{TOTAL | INCREMENT | Set input for total or incremental (with respect to the reference
INCREMENTAL } coordinates) values of the imperfect structural geometry.
Example:
NODAL_IMPEFECTIONS SETTINGS // 3D
NODE 2 TOTAL VALUES 0. 0. 0.001
NODE 3 INCREMENT VALUES 0. 0. 0.0015
…….
NODAL_IMPEFECTIONS SETTINGS // 2D
NODE 2 TOTAL VALUES 0. 0.001
NODE 3 INCREMENTAL VALUES 0. 0.0015
…….
Parameter Description
ID n Material identification
E.g.: ID 1
NAME “material_name” Material name in quotes, also for identification
E.g.: NAME „my_material“
&MATERIAL_TYPE_PARAMS Material type and type specific parameters
&MATERIAL_TYPE_PARAMS:
{ &LINEAR_ELASTIC_ISOTROPIC | &3DCEMENTITIOUS |
&3DNONLINCEMENTITIOUS | &3DNONLINCEMENTITIOUS2 |
&3DNONLINCEMENTITIOUS2VARIABLE |
&3DNONLINCEMENTITIOUS2USER |
&3DNONLINCEMENTITIOUS2SHCC |
&3DNONLINCEMENTITIOUS2SFATIGUE |
&3DNONLINCEMENTITIOUS3 | &SBETAMATERIAL |
ATENA Input File Format 79
&VON_MISES_PLASTICITY | &USER_MATERIAL |
&INTERFACE_MATERIAL | &REINFORCEMENT |
&REINFORCEMENT_WITH_CYCLING_BEHAVIOR |
&SMEARED_REINFORCEMENT | &SPRING |
&DRUCKER_PRAGER_PLASTICITY | &MICROPLANE |
&CREEP_MATERIAL | &COMBINED_MATERIAL |
&VARIABLE_MATERIAL |
&MATERIAL_WITH_TEMP_DEP_PROPERTIES |
&MATERIAL_WITH_RANDOM_FIELDS
&BEAM_MASONRY_MATERIAL | &BEAM_RC_MATERIAL |
&BEAM_REINF_BAR_MATERIAL}
Parameter Description
&LINEAR_ELASTIC_ISOTROPIC Linear elastic isotropic materials for 1D, Plane
Stress, Plane Strain, Axisymmetric and 3D
analyses
&3DCEMENTITIOUS Material suitable for rock or concrete like
materials.
&3DNONLINCEMENTITIOUS Materials suitable for rock or concrete like
materials. Enhanced &3DCEMENTITIOUS
material.
&3DNONLINCEMENTITIOUS2 Materials suitable for rock or concrete like
materials. This material is identical to
3DNONLINCEMENTITIOUS except that this
model is fully incremental.
&3DNONLINCEMENTITIOUS2VARI Materials suitable for rock or concrete like
ABLE materials. This material is identical to
3DNONLINCEMENTITIOUS2 except that
selected material parameters can be defined using
a time or load step function.
&3DNONLINCEMENTITIOUS2USE Materials suitable for rock or concrete like
R materials. This material is identical to
3DNONLINCEMENTITIOUS2 except that
selected material laws can be defined by user
curves.
&3DNONLINCEMENTITIOUS2SHC Strain Hardening Cementitious Composite
C material. Material suitable for fibre reinforced
concrete, such as SHCC and HPFRCC materials.
&3DNONLINCEMENTITIOUS2FATI Based on the 3DNONLINCEMENTITIOUS2
GUE material, suitable for fatigue analysis of rock or
concrete like materials.
&3DNONLINCEMENTITIOUS3 Materials suitable for rock or concrete like
materials. This material is an advanced version of
80
Parameter Description
Basic properties
Ex Elastic modulus.
Units: F/(l2)
Acceptable range: (0; maximal real number>
Default value: 210 x 103 f F / f l 2
idealisation is to be used.
Parameter Description
Basic properties
Ex Elastic modulus.
Units: F/(l2)
Acceptable range: (0; maximal real number>
Default value: 30 x 103 f F / f l 2
Units: F/(l2)
Acceptable range: <minimal real number; -2 FT)
Default value: -30 f F / f l 2
Parameter Description
Basic properties
Ex Elastic modulus.
Units: F/(l2)
Acceptable range: (0; maximal real number>
Default value: 30 x 103 f F / f l 2
Units: l
Acceptable range: <0; maximal real number>
Default value: 0.0
TENSION_STIFF x Tension stiffening
Units: none
Acceptable range: <0; 1>
Default value: 0.0
Compressive properties
EPS_CP x Plastic strain at compressive strength.
Units: none
Acceptable range: <minimal real number; 0>
Default value: -0.001
Generation formula: FC/E
88
Units: l
Acceptable range: <0; maximal real number>
Default value: 0.0
Miscellaneous properties
EXC x Eccentricity, defining the shape of the failure surface
Units:
Acceptable range: <0.5; 1.0>
Default value: 0.52
BETA x Multiplier for the direction of the plastic flow.
Units:
Acceptable range: <minimal real number; maximal real
number>
Recommended range: (-2; 2)
Default value: 0.0
RHO x Material density.
Units: M/l3
Acceptable range: <0; maximal real number>
Default value: 0.0023 f M / fl 3
ATENA Input File Format 89
Parameter Description
Basic properties
Ex Elastic modulus.
Units: F/(l2)
Acceptable range: (0; maximal real number>
Default value: 30 x 103 f F / f l 2
DAMPING_STIFF xK element group. They overwrite the same factor set for the whole
structure by SET command .
Parameter Description
Basic properties
Ex Elastic modulus.
Units: F/(l2)
Acceptable range: (0; maximal real number>
Default value: 30 x 103 f F / f l 2
Table 69: Parameters for MATERIAL TYPE „CC3DNonLinCementitious2User“. Suitable for rock
or concrete like materials
Parameter Description
Basic properties
E Elastic modulus.
Format: E x
Units: F/(l2)
Acceptable range: (0; maximal real number>
Default value: 30 x 103 f F / f l 2
Units: F/(l2)
Acceptable range: (0; -FC/2)
Default value: 3 f F / f l 2
2
Format: TENSION_SOFT_HARD_FUNCTION n
Units: none
Acceptable range: <1;maximal int number>
Default value: none
Generation formula: default function should have the following
points.
( 0.000 ; 1.00 )
0.75 GF
( ; 0.25 )
0.03 f t '
5 GF
( ; 0.00 )
0.03 f t '
where: GF 0.000025 FT
CHAR_SIZE_TENSION Characteristic size for which the various tensile functions are
100
valid.
Format: CHAR_SIZE_TENSION x
Units: l
Acceptable range: (0;maximal real number>
Default value: 0.03 fl
Generation formula: none
X_LOC_TENSION Strain value after, which the softening/hardening becomes
localized, and therefore adjustment based on element size is
needed.
Format: X_LOC_TENSION x
Units: none
Acceptable range: <0;maximal real number>
Default value: 0.0
Generation formula: none
CRACK_SPACING x Crack spacing – average distance between cracks after
localization. If zero crack spacing is assumed to be equal to
finite element size.
Units: l
Acceptable range: <0; maximal real number>
Default value: 0.0
TENSION_STIFF x Tension stiffening
Units: none
Acceptable range: <0; 1>
Default value: 0.0
ATENA Input File Format 101
Compressive properties
COMP_SOFT_ Index of the function defining the tensile hardening/softening
HARD_FUNCTION law. The horizontal axis represents strains and vertical axis
compressive strength, which should be normalized with respect
to f c' .
Format: COMP_SOFT_HARD_FUNCTION n
Units: none
Acceptable range: <1;maximal int number>
Default value: none
Generation formula: default function should have the following
points.
( 0.000 ; 0.25 )
( 0.5*FC/E ; 0.80 )
( FC/E ; 1.00 )
( FC/E - 0.005 ; 0.00 )
Format: CHAR_SIZE_COMP x
Units: l
Acceptable range: (0;maximal real number>
Default value: 0.10 fl
Generation formula: none
X_LOC_COMP Strain value after, which the softening/hardening becomes
localized, and therefore adjustment based on element size is
needed.
Format: X_LOC_COMP x
Units: none
Acceptable range: <0;-maximal real number>
Default value: -0.001
Generation formula: FC/E
102
Format: FC_REDUCTION_FUNCTION n
Units: none
Acceptable range: <1;maximal int number>
Default value: none
Generation formula: default function should have the following
points.
( 0.000 ; 1.0 )
( 0.001 ; 1.0 )
( 0.005 ; 0.6 )
( 0.01 ; 0.4 )
( 0.015 ; 0.3 )
( 0.05 ; 0.1 )
( 0.30 ; 0.02 )
ATENA Input File Format 103
Shear properties
SHEAR_STIFF_ Index of the function defining the shear retention factor
FUNCTION evolution based on tensile strain in the crack direction. The
horizontal axis represents strains and the vertical axis the
relative reduction of the shear modulus.
Format: SHEAR_STIFF_FUNCTION n
Units: none
Acceptable range: <1;maximal int number>
Default value: none
Generation formula: default function should have the following
points.
( 0.00000 ; 1.00 )
( 1.e-7 ; 1.00 )
( 1.e-6 ; 0.79 )
( 1.e-5 ; 0.58 )
( 0.00010 ; 0.36 )
( 0.001 ; 0.15 )
( 0.005 ; 0.01 )
Format: X_LOC_SHEAR x
Units: none
Acceptable range: <0;maximal real number>
Default value: 0.0
Generation formula: none
104
Format: SHEAR_STRENGTH_FUNCTION n
Units: none
Acceptable range: <1;maximal int number>
Default value: none
Generation formula: default function should have the following
points.
( 0.00000 ; 1.10 )
( 0.0001 ; 0.87 )
( 0.0005 ; 0.51 )
( 0.0010 ; 0.34 )
( 0.002 ; 0.20 )
( 0.003 ; 0.15 )
( 0.005 ; 0.09 )
( 0.010 ; 0.05 )
Tension-compression interaction
TENSILE_STRENGTH_ Index of the function defining the tensile strength reduction
RED_FUNCTION based on the compressive stress in other material directions.
The horizontal axis represents relative compressive stress
normalized with respect to f c' and the vertical axis the relative
reduction of the tensile strength with respect to ft ' .
Format: TENSILE_STRENGTH_RED_FUNCTION n
Units: none
Acceptable range: <1;maximal int number>
Default value: none
Generation formula: default function should have the following
points.
( 0.000 ; 1.00 )
( 1.000 ; 0.20 )
ATENA Input File Format 105
Miscellaneous properties
EXC Excentricity, defining the shape of the failure surface
Format: EXC x
Units:
Acceptable range: <0.5; 1.0>
Default value: 0.52
BETA Multiplier for the direction of the plastic flow.
Format: BETA x
Units:
Acceptable range: <minimal real number; maximal real
number>
Recommended range: (-2; 2)
Default value: 0.0
RHO Specific material density.
Format: RHO x
Units: M/l3
Acceptable range: <0; maximal real number>
Default value: 0.0023 f M / fl 3
Parameter Description
Basic properties
E Elastic modulus.
Format: E x
Units: F/(l2)
Acceptable range: (0; maximal real number>
Default value: 27 x 103 f F / f l 2
Fiber reinforcement
FIBER_VOLUME_FRA Volume fraction of the fibers.
CTION Format: FIBER_VOLUME_FRACTION x
Units: none
Acceptable range: <0;1>
108
Tensile properties
TENSION_SOFT_ Index of the function defining the tensile hardening/softening
HARD_FUNCTION law. The horizontal axis represents strains and vertical axis
tensile strength, which should be normalized with respect to ft ' .
Format: TENSION_SOFT_HARD_FUNCTION n
Units: none
Acceptable range: <1;maximal int number>
Default value: none
Generation formula: default function should have the following
points.
( 0.000 ; 1.00 )
( 0.040 ; 1.25 )
ATENA Input File Format 109
( 0.115 ; 0.00 )
CHAR_SIZE_TENSION Characteristic size for which the various tensile functions are
valid.
Format: CHAR_SIZE_TENSION x
Units: l
Acceptable range: (0;maximal real number>
Default value: 0.08 fl
Generation formula: none
X_LOC_TENSION Strain value after, which the softening/hardening becomes
localized, and therefore adjustment based on element size is
needed.
Format: X_LOC_TENSION x
Units: none
Acceptable range: <0;maximal real number>
Default value: 0.04
Generation formula: none
CRACK_SPACING x Crack spacing – average distance between cracks after
localization. If zero crack spacing is assumed to be equal to
finite element size.
Units: l
Acceptable range: <0; maximal real number>
Default value: 0.0
TENSION_STIFF x Tension stiffening
Units: none
Acceptable range: <0; 1>
Default value: 0.0
110
Compressive properties
COMP_SOFT_ Index of the function defining the tensile hardening/softening
HARD_FUNCTION law. The horizontal axis represents strains and vertical axis
compressive strength, which should be normalized with respect
to f c' .
Format: COMP_SOFT_HARD_FUNCTION n
Units: none
Acceptable range: <1;maximal int number>
Default value: none
Generation formula: default function should have the following
points.
( 0.000 ; 0.25 )
( 0.5*FC/E ; 0.80 )
( FC/E ; 1.00 )
( FC/E - 0.005 ; 0.00 )
Note: the x-values should be negative.
CHAR_SIZE_COMP Characteristic size for which the various compressive functions
are valid.
Format: CHAR_SIZE_COMP x
Units: l
Acceptable range: (0;maximal real number>
Default value: 0.15 fl
Tension-compression interaction
TENSILE_STRENGTH_ Index of the function defining the tensile strength reduction
RED_FUNCTION based on the compressive stress in other material directions.
The horizontal axis represents relative compressive stress
normalized with respect to f c' and the vertical axis the relative
reduction of the tensile strength with respect to ft ' .
Format: TENSILE_STRENGTH_RED_FUNCTION n
Units: none
Acceptable range: <1;maximal int number>
Default value: none
Generation formula: default function should have the following
points.
( 0.000 ; 1.00 )
( 1.000 ; 0.20 )
Miscellaneous properties
EXC Excentricity, defining the shape of the failure surface
Format: EXC x
Units:
Acceptable range: <0.5; 1.0>
Default value: 0.52
BETA Multiplier for the direction of the plastic flow.
Format: BETA x
Units:
Acceptable range: <minimal real number; maximal real
number>
Recommended range: (-2; 2)
Default value: 0.0
RHO Specific material density.
Format: RHO x
Units: M/l3
Acceptable range: <0; maximal real number>
Default value: 0.0023 f M / fl 3
Parameter Description
Basic properties
Ex Elastic modulus.
Units: F/(l2)
Acceptable range: (0; maximal real number>
Default value: 30 x 103 f F / f l 2
Tensile properties
GF x Specific fracture energy
Units: F/l
Acceptable range: (0; maximal real number>
Default value: 0.0001 f F / f l
Generation formula: GF 0.000025 FT
CRACK_SPACING x Crack spacing – average distance between cracks after
localization. If zero crack spacing is assumed to be equal
to finite element size.
Units: l
Acceptable range: <0; maximal real number>
Default value: 0.0
TENSION_STIFF x Tension stiffening
Units: none
Acceptable range: <0; 1>
Default value: 0.0
Compressive properties
EPS_CP x Plastic strain at compressive strength.
Units: none
Acceptable range: <minimal real number; 0>
Default value: -0.001
Generation formula: FC/E
ATENA Input File Format 115
Units: none
Acceptable range: <0; +>
Default value: 2.1
SHEAR_FACTOR x Shear factor that is used for the calculation of cracking
shear stiffness. It is calculated as a multiple of the
corresponding minimal normal crack stiffness that is based
on the tensile softening law.
Units: none
Acceptable range: <0; +>
Default value: 20
UNLOADING x Unloading factor, which controls crack closure stiffness.
Acceptable range: <0; 1)
0 - unloading to origin (default)
0.9 - unloading direction parallel to the initial elastic
stiffness
BETA_FATIGUE x Exponent for fatigue calculation.
Units: none
Acceptable range: (0; +>
Default value: 0.06
KSI_FATIGUE x Factor for fatigue damage calculation based on crack
opening and closing (COD).
When set to -1, it activates the trilinear stress based
damage evolution (and no COD based damage is
calculated).
Units: none
Acceptable range: (0; 1> or -1
Default value: 0.0001
IDEALISATION Defines the idealisation if material model with higher
dimension is to be used in a finite element with lower
dimension. For instance in case a 3D model is to be used
in 2D configuration.
Units: none
Acceptable range: { 1D, PLANE_STRESS,
PLANE_STRAIN, AXISYMMETRIC, 3D }
Default value: program tries to determine a suitable
idealisation based on the dimension of the material model
and the dimension of the finite element where it is used.
So in most cases it is not needed to use this command. In
ATENA Input File Format 117
Parameter Description
Basic properties
Ex Elastic modulus
Units: MPa
Acceptable range: (0; maximal real number>
Recommended value : From table below
{ MU | POISSON | NY } x Poisson’s ratio (ν)
Units: none
Acceptable range: <-1; 0.5)
Recommended value : From table below
{ FT | RT | F_T | R_T} x Tensile strength (ft)
Units: MPa
Acceptable range: (0; -FC/2)
Recommended value : From table below
118
Units: none
Acceptable range: <minimal real number; 0>
Recommended value : From table below Generation
formula: (FC/E) * (1-2*MU)
{ FC0 | F_C0 | RC0 | R_C0 } x Onset of non-linear behavior in compression (f )
co
Units: MPa
Acceptable range: (FC,-2 FT)
Recommended value : From table below
SOFT_T x Slope of softening curve t
Units: none
Acceptable range: <0; maximal real number>
Recommended value : From table below
ATENA Input File Format 119
Parameter Description
Basic
Ex Elastic modulus.
Units: F/(l2)
122
Tension
ISOFT x Type of tension softening.
Units: none
Acceptable range: <1.0;5.0>
1.0=Exponential
2.0=Linear
3.0=Local strain
4.0=SFRC
5.0=SFRC local strain
Default value: 1.0
Case ISOFT = 1.0 (Exponential)
GF x Specific fracture energy
Units: F/l
ATENA Input File Format 123
CS x Tension-compression interaction.
Units: none
Acceptable values: 0.2, 0.4, 0.6
0.6=Linear
0.4=Hyperbola A
0.2=Hyperbola B
Default value: 0.6 (Linear)
ROTATED CRACKS Activates rotated crack model. If not used fixed crack model is
considered.
Units: none
Acceptable range: none
Default value: not used
Miscellaneous
DAMPING_MASS xM Mass and stiffness damping factors specified for indiviual
DAMPING_STIFF xK element group. They overwrite the same factor set for the whole
structure by SET command .
RHO x Material density.
Units: M/l3
Acceptable range: <0; maximal real number>
Default value: 0.0023 f M / fl 3
Parameter Description
Basic properties
Ex Elastic modulus.
Units: F/(l2)
Acceptable range: (0; maximal real number>
Default value: 210 x 103 f F / f l 2
K2 Bauschinger memory
Format: K2 x
Units: none
Acceptable range: (0; maximal real number)
Default value: 1000
Miscellaneous properties
RHO x Material density.
Units: M/l3
Acceptable range: <0; maximal real number>
Default value: 0.00785 f M / fl 3
Syntax:
&DRUCKER_PRAGER_PLASTICITY:
TYPE “CC3DDruckerPragerPlasticity” { E x | { MU | POISSON | NY } x K x |
ALPHA_DP x | WD x | BETA x | RHO x | ALPHA x | IDEALISATION
{ 1D, PLANE_STRESS, PLANE_STRAIN, AXISYMMETRIC, 3D} |
DAMPING_MASS xM DAMPING_STIFF xK}+
ATENA Input File Format 129
The parameters for this material model can be generated based on compressive and tensile
strength of the material R c and R t (see Table 75). These values should be specified in MPa
and then transformed to the current units.
Parameter Description
Basic properties
Ex Elastic modulus.
Units: F/(l2)
Acceptable range: (0; maximal real number>
Default value: 30 x 103 f F / f l 2
1
Generation formula: K Rc ALPHA_DP f F / f l 2
3
Compressive properties
WD x Critical compressive displacement
Units: l
Acceptable range: <0; maximal real number>
Default value: -0.0005 fl
130
Miscellaneous properties
BETA x Multiplier for the direction of the plastic flow.
Units: none
Acceptable range: <minimal real number; maximal real
number>
Recommended range: (-2; 2)
Default value: 0.0
RHO x Material density.
Units: M/l3
Acceptable range: <0; maximal real number>
Default value: 0.0023 f M / fl 3
Parameter Description
Basic properties (inherited from elastic material)
Ex Elastic modulus.
Units: F/(l2)
Acceptable range: (0; maximal real number>
Default value: 210 x 103 f F / f l 2
Parameter Description
Basic properties
{K_NN | KNN} x Normal stiffness. Units: F/l3
Acceptable range: (0; maximal real number>
Default value: 200 x 106 f F / f l 3
5
Available starting from ATENA version 4.3.1.
ATENA Input File Format 133
Default value: 0 f F / f l 2
COHESION x Cohesion.
Units: F/(l2)
Acceptable range: <0; maximal real number>
Default value: 0.0 f F / f l 2
Parameter Description
Basic properties
Ex Elastic modulus.
Units: F/(l2)
Acceptable range: (0; maximal real number>
Default value: 210 x 103 MPa
FUNCTION a Function which defines uniaxial stress-strain relationship.
Relationship should be defined as a set of points starting from
(0, 0) and only positive values should be specified. Same
relationship will be used in compression.
Units: none
Acceptable range: (1; maximal integer>
Default value: none, see command &FUNCTION.
RHO x Material density.
Units: M/l3
Acceptable range: <0; maximal real number>
Default value: 0.00785 f M / fl 3
Parameter Description
Basic properties
Ex Elastic modulus.
Units: F/(l2)
Acceptable range: (0; maximal real number>
Default value: 210 x 103 MPa
FUNCTION n Function which defines uniaxial stress-strain relationship.
Relationship should be defined as a set of points starting from
(0, 0) and only positive values should be specified. Same
relationship will be used in compression.
Units: none
Acceptable range: (1; maximal integer>
Default value: none, see command &FUNCTION
Rx Bauschinger effect exponent of Menegotto-Pinto model.
Units: none
Acceptable range: (0; maximal real number>
Default value: 20
C1 x Menegotto-Pinto model parameter
Units: none
Acceptable range: (0; 1)
Default value: 0.925
ATENA Input File Format 137
&SMEARED_REINFORCEMENT
TYPE “CCSmearedReinf” { E x | FUNCTION n | RATIO x | DIRECTION x1 x2 [x3] | RHO x
| ALPHA x | F_MULTIP x }+
Parameter Description
Basic properties
Ex Elastic modulus.
Units: F/(l2)
Acceptable range: (0; maximal real number>
Default value: 210 x 103 MPa
FUNCTION a Function which defines uniaxial stress-strain relationship.
Relationship should be defined as a set of points starting from
(0, 0) and only positive values should be specified. Same
relationship will be used in compression.
Units: none
Acceptable range: (1; maximal integer>
Default value: none, see command &FUNCTION.
RATIO x Cross-sectional area ratio of the smeared reinforcement with
respect to the base material.
Units: none
Acceptable range: <0;1>
Default value: 0.01
138
DIRECTION x1 x2 [x3] Unit vector defining the smeared reinforcement direction. The
third component x3 is required in case of 3D analysis.
Units: l
Acceptable range: <minimal real; maximal real number>
Default value: 1 0 [0]
RHO x Material density.
Units: M/l3
Acceptable range: <0; maximal real number>
Default value: 0.00785 f M / fl 3
Parameter Description
Basic properties
Ex Elastic modulus.
Units: F/(l2)
Acceptable range: (0; maximal real number>
Default value: 210 x 103 MPa
ATENA Input File Format 139
Parameter Description
Basic properties
140
Kx Initial stiffness.
Units: F/l
Acceptable range: (0; maximal real number>
Default value: 1000.0
DAMPING_MASS xM Mass and stiffness damping factors specified for indiviual
DAMPING_STIFF xK element group. They overwrite the same factor set for the whole
structure by SET command .
FUNCTION n Function which defines uniaxial spring relationship.
Relationship should be defined as a set of points starting in
compression passing through (0, 0) and extending into tension. a
Units: none
Acceptable range: (1; maximal integer>
Default value: none, see command &FUNCTION
&MICROPLANE4
ATENA Input File Format 141
Parameter Description
Basic properties
Ex Elastic modulus.
Units: F/(l2)
Acceptable range: (0; maximal real number>
Default value: 30 x 103 f F / f l 2
Default value: 15
K4 x Microplane parameter k4.
Units: None
Acceptable range: <0; maximal real number>
Default value: 150
BAND x Crack band size.
Units: l
Acceptable range: <0; maximal real number>
Default value: 0.003 fl
Miscellaneous properties
RHO x Material density.
Units: M/l3
Acceptable range: <0; maximal real number>
Default value: 0.00785 f M / fl 3
Cx x
Default values: (see theory manual for details)
c1 =6.20e-1 Normal bound. param.
c2 =2.76 Normal bound. param.
c3 =4.00 Normal plasticity (EB_N)
c4 =70.00 Strain ratio: normal/vol
c5 =2.50 Tens. dev. bound. param.
c6 =1.30 Comp. dev. bound. yield
c7 =50.00 Deviatoric plasticity (EB_D)
c8 =8.00 Compressive strength (FCP)
c9 =1.30 Dev. bound. param.
c10 =7.30e-1 Fric.b. initial slope
c11 =2.00e-1 Fric.b.\sig_N inter.@\sig_V=0
c12 =7.00e+3 Fric.b.Speed\sig_N goes zero
c13 =2.30e-1 Tensile vol.b. vert. scalar
c14 =8.00e-1 Tensile vol.b. slope
c15 =1.00 Tensile vol.b. horiz. yield
c16 =2.00e-2 Unl. volumetric coeff.
c17 =1.00e-2 Unl. volumetric coeff.
c18 =1.000 Tensile vol.b. unload.coeff
c19 =0.40 Unloading slope interpolator
c20 =14.00e-2 Residual strength
c21 =0.990 Unloading slope Int. in tens
&CCM4:
TYPE “CCM4” {&CCM4Params}+
&CCM4Params:
{ E x | Nplane n | K1 x | K2 x | K3 x | K4 x | ES0 x | VA x | FC x | TSH x | PSI x |
ETA_V x | ETA_D x | ETA_N x | MY_U1 x | IDEALISATION { 1D,
PLANE_STRESS, PLANE_STRAIN, AXISYMMETRIC, 3D}
Parameter Description
Basic properties
Ex Elastic modulus.
Units: F/(l2)
Acceptable range: (0; maximal real number>
144
Units: 1/T
Acceptable range: <0; maximal real number>
Default value: 0.000012
IDEALISATION Defines the idealisation if material model with higher dimension
is to be used in a finite element with lower dimension. For
instance in case a 3D model is to be used in 2D configuration.
Units: none
Acceptable range: { 1D, PLANE_STRESS,
PLANE_STRAIN, AXISYMMETRIC, 3D }
Default value: program tries to determine a suitable idealisation
based on the dimension of the material model and the dimension
of the finite element where it is used. So in most cases it is not
needed to use this command. In certain cases however the
program cannot determine correctly the idealisation to use, such
a case is for instance if a 3D model is to be used in 2D element.
Then it is necessary to directly specify if plane stress or strain
idealisation is to be used.
&CCM4R:
TYPE “CCM4R” { &CCM4RParams | &CCM4Params }+
&CCM4RParams:
{ REF_TEMPER x | QR x | CR0 x | CR2 x }
Parameter Description
REF_TEMPER x Reference temperature.
Units: 0C
Default value: 25 0C
QR x Activation energy constant.
Units: 0 K
Default value: 1000 0 K
CR0 x Boundary rate shape CR0 constant.
1
Units:
sec
Default value: 106 sec1 6.4
ATENA Input File Format 147
&CCM4RC:
TYPE “CCM4R” { &CCM4RCParams | &CCM4RParams | &CCM4Params }+
&CCM4RCParams:
{ TIME0 x | HUMIDITY0 x | TEMPERATURE0 | TAU1 x |
NUMBER_MAXWELL n | Q1 x | Q2 x | Q3 x | Q4
x | WC x | CC x | AC x | C x | C1 x | CREEP_DEGREE x | VOLUME_POW x |
LAMBDA0 x }
Parameter Description
TIME0 x Initial time.
Units: Days
Default value: 28 days
TEMPERATURE Material initial temperature
Units: 0C
Default value: 25 0C
HUMIDITY Material initial humidity.
Units: None
Default value: 0.94
TAU1 x Te smallest relaxation time.
Units: days
Default value: 1.E-6 days
NUMBER_MAXWELL Number of Maxwell or Kelvin units
n
Units: None
Default value: 14
Q1 x Creep parameter Q1, (refer to Bazant & Baweja Model B3). If
negative, the parameter is estimated according to the above
mentioned creep model.
1
Units:
MPa
Default value: -1
148
The parameter BASE contains material type to be used for the short term material model. See
Table 55 for more information about the available material models for this parameter. After
that the parameters of the short term material will follow.
150
&CCModelB3_DATA
CCModelB3 { CONCRETE concrete_type | THICKNESS thick | FCYL28 fcyl28 | E28
e28 | HUMIDITY humidity | DENSITY density | AC ac | WC wc | [SHAPE]
FACTOR sfactor | {WATER | AIR | STEAM} [CURING] | [END] [OF]
[CURING] TIME endcuring | TOTAL_LOSS total_water_loss | {LOAD |
CURRENT} [TIME] time | {LOSS | SHRINKAGE | COMPLIANCE}
measured_val }+
Parameter Description
CONCRETE Type of concrete. Only type 1 and 3 are supported for static and
concrete_type types 1-4 for transport analysis. More information available in
the Atena Theory Manual.
Default value: 1
THICKNESS thick Ratio volume [m3] / surface area [m2] of cross section. For long
elements it is approximately cross sectional area [m2] / perimeter
[m].
Default value: 0.0767 [m].
FCYL28 fcyl28 Cylindrical material strength in compression [kPa].
Default value: 35100 [kPa].
E28 e28 Short-term material Young modulus at 28 days, i.e. inverse
compliance at 28.01 days loaded at 28 days [kPa].
Default value: calculated from fcyl28.
HUMIDITY humidity Ambient relative humidity (0.3...1).
Default value: 0.780
DENSITY density Concrete density [kg/m3].
Default value: 2125. [kg/m3].
AC ac Total aggregate/cement ratio.
Default value: 7.04
WC wc Water/cement ratio.
Default value: 0.63
[SHAPE] FACTOR Cross section shape factor. It should be 1, 1.15, 1.25, 1.3, 1.55
sfactor for slab, cylinder, square prism, sphere, cube, respectively.
Default value 1.25
{WATER | AIR | Curing conditions, either under in water or air under normal
STEAM} [CURING] temperature conditions (WATER | AIR) or steam condition
(=STEAM).
ATENA Input File Format 151
&CCModelB3Improved_Data
CCModelB3Improved { CONCRETE concrete_type | THICKNESS thick | FCYL28
f cyl ,28 f f G
| E28 E28 | FCYL0_28 cyl 0,28 | FT28 t ,28 | GF28 f ,28 | ALPHA |
HUMIDITY humidity | DENSITY density | AC ac | WC wc | [SHAPE]
FACTOR sfactor | {WATER | AIR | STEAM} [CURING] | [END] [OF]
[CURING] TIME endcuring | EPS_A_INF a , | TAU_A a | TIME_S ts |
h
H_A_INF a , | TOTAL_LOSS total_water_loss | {LOAD | CURRENT} [TIME]
time | {LOSS | SHRINKAGE | COMPLIANCE} measured_val { HISTORY
[TIME] time [HUMIDITY] humid [TEMPERATURE] temper }+ }+
Parameter Description
CONCRETE Type of concrete. Only type 1 and 3 are supported.
concrete_type
Default value: 1
THICKNESS thick Ratio volume [m3] / surface area [m2] of cross section. For long
elements it is approximately cross sectional area [m2] / perimeter
[m].
Default value: 0.0767 [m].
FCYL28 f cyl ,28 Cylindrical material strength in compression f cyl (28 days ) . This
value is crucial for the creep model's prediction, i.e. prediction
of material compliance (t , t ') and cylindrical compression
strength f cyl (t ) , shrinkage etc. The ratio of f cyl (t ) / f cyl (28 days )
152
{WATER | AIR | Curing conditions, either under in water or air under normal
STEAM} [CURING] temperature conditions (WATER | AIR) or steam condition
(=STEAM).
Default value: AIR
[END] [OF] [CURING] Time at beginning of drying, i.e. end of curing. [days].
TIME endcuring Default value: 7 [days].
TOTAL_LOSS Total water loss (at zero humidity and infinite time).
total_water_loss
Default: 0 [kg]
{LOAD | CURRENT} Current or load time for the subsequent measured value.
[TIME] time Default: 0 [days]
{LOSS | SHRINKAGE | Measured water loss (at current humidity) | shrinkage | material
COMPLIANCE} compliance measured_val for previously specified load and
measured_val current time. Units of water loss must correspond to units of
total_water_loss , shrinkage is dimension-less and compliance is
input in kPa-1.
{ HISTORY [TIME] For each entry of material history the data time , temper and
time [HUMIDITY] humid must be input. If the data keywords are used, then it
humid doesn’t matter in which order the 3 data are input. Otherwise the
[TEMPERATURE] indicated order is assumed. The units are days, degrees Celsius
temper }+ and dimension less humidity (in interval 0.3..1).
EPS_A_INF a , Autogenous shrinkage at infinity time, (typically negative!).
Default value= -0.
TAU_A a Half-time of autogenous shrinkage. Default value =30 days
&CCModelFIB_MC2010_DATA
CCModelFIB_MC2010 { CEMENT_CLASS { 32.5N | 32.5R | 42.5N | 42.5R | 52.5N |
52.5R } | AGGREAGETE { BASALTDENSELIMESTONE | QUARTZITE |
LIMESTONE | SANDSTONE | LIGHTWEIGHTSANDSTONE } |
THICKNESS thick | FCYL28 f cyl ,28 | E28 E28 | FCYL0_28 f cyl 0,28 | FT28 ft ,28 |
GF28 G f ,28 | ALPHA | HUMIDITY humidity | DENSITY density | [END]
[OF] [CURING] TIME endcuring | {LOAD | CURRENT} [TIME] time | {
SHRINKAGE | COMPLIANCE} measured_val { HISTORY [TIME] time
[HUMIDITY] humid [TEMPERATURE] temper }+ }+
154
Parameter Description
CEMENT_CLASS { Type of cement, see e.g. http://www.cis.org.rs/en/cms/about-
32.5N | 32.5R | 42.5N | cement/standardization-of-cement :
42.5R | 52.5N | 52.5R } Strength classes of cement
Cements are according to standard strength grouped into three
classes, they being:
• Class 32,5
• Class 42,5
• Class 52,5
Three classes of early strength are defined for each class of
standard strength:
• Class with ordinary early strength – N
• Class with high early strength – R
• Class with low early strength – L
Class L can be applied only on CEM III cements.
remains unchanged.
Default value: 0 [MPa]
GF28 G f ,28 The parameter fracture energy G f (28 days ) . If specified by a
positive value, this value is to calculate G f (t ) and overide the
corresponding value in the base material. . If it is specified as
any negative value, G f (28 days ) is calculated by FIB_MC2010
based on f cyl (28 days) . Othewise, the value in the base material
remains unchanged.
Default value: 0 [MPa]
FT28 ft ,28 The parameter tensile strength ft (28 days ) . If specified by a
positive value, this value is used to calculate ft (t ) and overide
the corresponding value in the base material. If it is specified as
any negative value, ft (28 days ) is calculated by FIB_MC2010
based on f cyl (28 days) . Othewise, the value in the base material
remains unchanged.
Default value: 0 [MPa]
E28 E28 Short-term material Young modulus at 28 days, i.e. inverse
compliance at 28.01 days loaded at 28 days [kPa]. It is used by
the creep model to predict material compliance (t , t ') . If
unspecified, the model calculates its value based on fcyl28.
Default value: calculated from fcyl28.
ALPHA Coefficient of thermal expansion to be used for calculation
t (T ) within the creep material.
Default value: 0
HUMIDITY humidity Ambient relative humidity (0.3…1).
Default value: 0.780
DENSITY density Concrete density [kg/m3].
Default value: 2125. [kg/m3].
[END] [OF] [CURING] Time at beginning of drying, i.e. end of curing. [days].
TIME endcuring Default value: 7 [days].
{LOAD | CURRENT} Current or load time for the subsequent measured value.
[TIME] time Default: 0 [days]
{ SHRINKAGE | Measured shrinkage | material compliance measured_val for
COMPLIANCE} previously specified load and current time. Units of water loss
measured_val must correspond to units of total_water_loss , shrinkage is
dimension-less and compliance is input in kPa-1.
{ HISTORY [TIME] For each entry of material history the data time , temper and
time [HUMIDITY] humid must be input. If the data keywords are used, then it
156
humid doesn’t matter in which order the 3 data are input. Otherwise the
[TEMPERATURE] indicated order is assumed. The units are days, degrees Celsius
temper }+ and dimension less humidity (in interval 0.3..1).
&CCModelEN1992_DATA
CCModel EN1992 { CEMENT_CLASS { 32.5N | 32.5R | 42.5N | 42.5R | 52.5N |
52.5R } | AGGREAGETE { BASALTDENSELIMESTONE | QUARTZITE |
LIMESTONE | SANDSTONE | LIGHTWEIGHTSANDSTONE } |
THICKNESS thick | FCYL28 f cyl ,28 | E28 E28 | FCYL0_28 f cyl 0,28 | FT28 ft ,28 |
GF28 G f ,28 | ALPHA | HUMIDITY humidity | DENSITY density | [END]
[OF] [CURING] TIME endcuring | {LOAD | CURRENT} [TIME] time | {
SHRINKAGE | COMPLIANCE} measured_val { HISTORY [TIME] time
[HUMIDITY] humid [TEMPERATURE] temper }+ }+
Parameter Description
CEMENT_CLASS { Type of cement, see e.g. http://www.cis.org.rs/en/cms/about-
32.5N | 32.5R | 42.5N | cement/standardization-of-cement :
42.5R | 52.5N | 52.5R } Strength classes of cement
Cements are according to standard strength grouped into three
classes, they being:
• Class 32,5
• Class 42,5
• Class 52,5
Three classes of early strength are defined for each class of
standard strength:
• Class with ordinary early strength – N
• Class with high early strength – R
• Class with low early strength – L
Class L can be applied only on CEM III cements.
Default value: 0
HUMIDITY humidity Ambient relative humidity (0.3…1).
Default value: 0.780
DENSITY density Concrete density [kg/m3].
Default value: 2125. [kg/m3].
[END] [OF] [CURING] Time at beginning of drying, i.e. end of curing. [days].
TIME endcuring Default value: 7 [days].
{LOAD | CURRENT} Current or load time for the subsequent measured value.
[TIME] time Default: 0 [days]
{ SHRINKAGE | Measured shrinkage | material compliance measured_val for
COMPLIANCE} previously specified load and current time. Units of water loss
measured_val must correspond to units of total_water_loss , shrinkage is
dimension-less and compliance is input in kPa-1.
{ HISTORY [TIME] For each entry of material history the data time , temper and
time [HUMIDITY] humid must be input. If the data keywords are used, then it
humid doesn’t matter in which order the 3 data are input. Otherwise the
[TEMPERATURE] indicated order is assumed. The units are days, degrees Celsius
temper }+ and dimension less humidity (in interval 0.3..1).
&CCModelBP_KX_DATA
CCModelBP_KX { CONCRETE concrete_type | THICKNESS thick | FCYL28 fcyl28
| E28 e28 | HUMIDITY humidity | DENSITY density | AC ac | WC wc |
[SHAPE] FACTOR sfactor | {WATER | AIR | STEAM} [CURING] | [END]
[OF] [CURING] TIME endcuring | {LOAD | CURRENT} [TIME] time |
{SHRINKAGE | COMPLIANCE} measured_val { HISTORY [TIME] time
[HUMIDITY] humid [TEMPERATURE] temper }+ }+
Parameter Description
CONCRETE Type of concrete. Only type 1 and 3 are supported.
concrete_type
Default value: 1
THICKNESS thick Ratio volume [m3] / surface area [m2] of cross section. For long
elements it is approximately cross sectional area [m2] / perimeter
[m].
Default value: 0.0767 [m].
FCYL28 fcyl28 Cylindrical material strength in compression [kPa].
Default value: 35100 [kPa].
ATENA Input File Format 159
&CCModelACI78_DATA
160
Parameter Description
CONCRETE Type of concrete. Only type 1 and 3 are supported.
concrete_type
Default value: 1
THICKNESS thick Ratio volume [m3] / surface area [m2] of cross section. For long
elements it is approximately cross sectional area [m2] / perimeter
[m].
Default value: 0.0767 [m].
FCYL28 fcyl28 Cylindrical material strength in compression [kPa].
Default value: 35100 [kPa].
HUMIDITY humidity Ambient relative humidity (0.3...1).
Default value: 0.780
DENSITY density Concrete density [kg/m3].
Default value: 2125. [kg/m3].
AC ac Total aggregate/cement ratio.
Default value: 7.04
WC wc Water/cement ratio.
Default value: 0.63
AS as Total aggregate/find sand ratio.
Default value 2.8
SLUMP slump Slump value [m].
Default value: 0.012m
AIR_CONTENT air Air content [%]:
Default value: 5%.
{WATER | AIR | Curing conditions, either under in water or air under normal
STEAM} [CURING] temperature conditions (WATER | AIR) or steam condition
(=STEAM).
Default value: AIR
[END] [OF] [CURING] Time at beginning of drying, i.e. end of curing. [days].
TIME endcuring Default value: 7 [days].
{LOAD | CURRENT} Current or load time for the subsequent measured value.
ATENA Input File Format 161
&CCModelCEB_FIP78_DATA
CCModelCEB_FIP78 { THICKNESS thick | FCYL28 fcyl28 | E28 e28 | HUMIDITY
humidity | [END] [OF] [CURING] TIME endcuring | {LOAD | CURRENT}
[TIME] time | SHRINKAGE measured_val }+
Parameter Description
THICKNESS thick Ratio volume [m3] / surface area [m2] of cross section. For long
elements it is approximately cross sectional area [m2] / perimeter
[m].
Default value: 0.0767 [m].
FCYL28 fcyl28 Cylindrical material strength in compression [kPa].
Default value: 35100 [kPa].
E28 e28 Short-term material Young modulus at 28 days, i.e. inverse
compliance at 28.01 days loaded at 28 days [kPa].
Default value: calculated from fcyl28.
HUMIDITY humidity Ambient relative humidity (0.3..1).
Default value: 0.780
[END] [OF] [CURING] Time at beginning of drying, i.e. end of curing. [days].
TIME endcuring Default value: 7 [days].
{LOAD | CURRENT} Current or load time for the subsequent measured value.
[TIME] time Default: 0 [days]
SHRINKAGE Measured (at current humidity) shrinkage measured_val for
measured_val previously specified load and current time. Unit of shrinkage is
dimension-less.
&CCModelCSN731202_DATA
CCModelCSN731202 { CONCRETE concrete_type, THICKNESS thick | FCYL28
fcyl28 | E28 e28 | HUMIDITY humidity | [END] [OF] [CURING] TIME
162
Parameter Description
CONCRETE Type of concrete. Only type 1 and 3 are supported.
concrete_type
Default value: 1
THICKNESS thick Ratio volume [m3] / surface area [m2] of cross section. For long
elements it is approximately cross sectional area [m2] / perimeter
[m].
Default value: 0.0767 [m].
FCYL28 fcyl28 Cylindrical material strength in compression [kPa].
Default value: 35100 [kPa].
E28 e28 Short-term material Young modulus at 28 days, i.e. inverse
compliance at 28.01 days loaded at 28 days [kPa].
Default value: calculated from fcyl28.
HUMIDITY humidity Ambient relative humidity (0.3...1).
Default value: 0.780
[END] [OF] [CURING] Time at beginning of drying, i.e. end of curing. [days].
TIME endcuring Default value: 7 [days].
{ HISTORY [TIME] For each entry of material history the data time , temper and
time [HUMIDITY] humid must be input. If the data keywords are used, then it
humid doesn’t matter in which order the 3 data are input. Otherwise the
[TEMPERATURE] indicated order is assumed. The units are days, degrees Celsius
temper }+ and dimension less humidity (in interval 0.3..1).
{LOAD | CURRENT} Current or load time for the subsequent measured value.
[TIME] time Default: 0 [days]
SHRINKAGE Measured shrinkage measured_val for previously specified load
measured_val and current time. Unit of shrinkage is dimension-less.
&CCModelBP1_DATA
CCModelBP1 { CONCRETE concrete_type | THICKNESS thick | FCYL28 fcyl28 |
HUMIDITY humidity | AC ac | WC wc | GS gs | SC sc | SA sa | CEMENT
[MASS] cement_mass | [SHAPE] FACTOR sf | {STEAM | WATER |AIR}
[CURING] | [END] [OF] [CURING] TIME endcuring | { LOAD | CURRENT }
[TIME] time SHRINKAGE measured_val }+
ATENA Input File Format 163
Parameter Description
CONCRETE Type of concrete. Only type 1 and 3 are supported.
concrete_type
Default value: 1
THICKNESS thick Ratio volume [m3] / surface area [m2] of cross section. For long
elements it is approximately cross sectional area [m2] / perimeter
[m].
Default value: 0.0767 [m].
FCYL28 fcyl28 Cylindrical material strength in compression [kPa].
Default value: 35100 [kPa].
HUMIDITY humidity Ambient relative humidity (0.3...1).
Default value: 0.780
AC ac Total aggregate/cement ratio.
Default value: 7.04
WC wc Water/cement ratio.
Default value: 0.63
GS gs Coarse/fine aggregate ratio.
Default value: 1.3
SC sc Fine aggregate/cement ratio.
Default value: 1.8
SA sa Fine/total aggregate ratio.
Default value: 0.4
CEMENT [MASS] Cement content.
cement_mass
Default value: 320. kg/m3
[SHAPE] FACTOR sf Cross section shape factor. It should be 1, 1.15, 1.25, 1.3, 1.55
for slab, cylinder, square prism, sphere, cube, respectively.
Default value 1.25
{STEAM | WATER Curing conditions, either under in water or air under normal
|AIR} [CURING] temperature conditions (WATER | AIR) or steam condition
(=STEAM).
Default value: AIR
[END] [OF] [CURING] Time at beginning of drying, i.e. end of curing. [days].
TIME endcuring Default value: 7 [days].
{LOAD | CURRENT} Current or load time for the subsequent measured value.
[TIME] time Default: 0 [days]
164
&CCModelBP2_DATA
CCModelBP2 { CONCRETE concrete_type | THICKNESS thick | FCYL28 fcyl28 |
HUMIDITY humidity | AC ac | WC wc | GS gs | SC sc | SA sa | [SHAPE]
FACTOR sf | {STEAM | WATER |AIR} [CURING] | [END] [OF] [CURING]
TIME time | { LOAD | CURRENT } [TIME] xx SHRINKAGE measured_val }+
Parameter Description
CONCRETE Type of concrete. Only type 1 and 3 are supported.
concrete_type
Default value: 1
THICKNESS thick Ratio volume [m3] / surface area [m2] of cross section. For long
elements it is approximately cross sectional area [m2] / perimeter
[m].
Default value: 0.0767 [m].
FCYL28 fcyl28 Cylindrical material strength in compression [kPa].
Default value: 35100 [kPa].
HUMIDITY humidity Ambient relative humidity (0.3...1).
Default value: 0.780
AC ac Total aggregate/cement ratio.
Default value: 7.04
WC wc Water/cement ratio.
Default value: 0.63
GS gs Coarse/fine aggregate ratio.
Default value: 1.3
SC sc Fine aggregate/cement ratio.
Default value: 1.8
SA sa Fine/total aggregate ratio.
Default value: 0.4
[SHAPE] FACTOR sf Cross section shape factor. It should be 1, 1.15, 1.25, 1.3, 1.55
for slab, cylinder, square prism, sphere, cube, respectively.
ATENA Input File Format 165
&CCModelGeneral_Data
CCModelGeneral { T’ t' | T t | FI fi | EPS eps | FCYL fcyl }+
Parameter Description
T’ t' Set effective loading time t' for following data.
Default value: none
Units: t.
T t Set effective observation time t for following data, i.e. a time,
when an input value is measured.
Default value: none
Units: t.
FI fi Value of material compliance fi(t.t’) for times t,t’ .
Default value: none.
Units: 1/S
EPS eps Material shrinkage eps(t) at time of observation t.
Default value: none
Units: none
FCYL fcyl Current cylindrical strength in compression fcyl(t’) pertinent for
loading time t’ . Note that the value is input as positive value!
Default value: none
Units: S
166
Parameter Description
Basic properties
COMPONENT id Id of the previously defined material, which is to be used a one
component of the combined/composite material.
Units: none
Acceptable range: (1; maximal integer>
Default value: none
RATIO x Relative contribution of this material to the overall behavior of
the combined composite material.
Units: none
Acceptable range: <0; maximal real number>
Default value: 1.0
Parameter Description
Basic properties
BASE id Id of the previously defined base material, whose parameters
ATENA Input File Format 167
Parameter Description
Basic properties
BASE id Id of the previously defined base material, whose parameters
will be modified based on the thermal loading and the provided
function. Only the following materials should be used as a base
material:
CC3DNonLinCementitious2,
CC1DElastIstotropic,
CCPlaneStressElastIsotropic,
CCPlaneStrainElastIsotropic,
CC3DelastIsotropic,
CCASymElastIsotropic,
CC3DDruckerPragerPlasticity,
CC3DBiLinearSteelVonMises,
CCReinforcement,
CCSmearedReinf
Units: none
Acceptable range: (1; maximal integer>
Default value: none
PARAM | Parameter name from the base material whose values will
PARAMETER “name” change based on the thermal loading and provided function. The
original value of this parameter in the base material is
overwritten by the values in the function. The base material
should not be used in any other combined material as well as
a stand alone material. Otherwise results are unpredictable.
Units: none
Acceptable range: any string
Default value: none
ATENA Input File Format 169
Parameter Description
Basic properties
BASE id Id of the previously defined base material, whose parameters
will be modified based on the thermal loading and the provided
function. Only the following materials should be used as a base
material:
CC3DNonLinCementitious2,
CC1DElastIstotropic,
CCPlaneStressElastIsotropic,
CCPlaneStrainElastIsotropic,
CC3DelastIsotropic,
CCASymElastIsotropic,
CC3DDruckerPragerPlasticity,
CC3DBiLinearSteelVonMises,
CCReinforcement,
CCSmearedReinf
Units: none
Acceptable range: (1; maximal integer>
Default value: none
FILENAME “name” File name containing the spatial distribution of material
parameters.
Units: none
Acceptable range: any string
Default value: none
Parameter Description
Ex Young modulus.
Units: stresses
Default value: 0
MU x Poisson ratio
Units: none
Default value: 0
RHO x Material density
Units: mass/volume
Default value: 0
ALPHA x Coefficient of thermal expansion
Units: 1/T
cceptable range: <0; maximal real number>
Default value: 0.000012
DAMPING_MASS xM Mass and stiffness damping factors specified for indiviual
DAMPING_STIFF xK element group. They overwrite the same factor set for the whole
structure by SET command .
F_K x Characteristic material compressive strength, (negative). This
input is not used, if the corresponding design value is given.
Units: stresses
Default value: 0
F_VK0 x Characteristic material initial shear strength, (positive). This
input is not used, if the corresponding design value is given.
Units: stresses
Default value: 0
COEFF_F_VK x Coefficient for normall stress to calculate F_VK.
Units: none
Default value: 0.4
F_VLT x Characteristic material limit shear strength – constant part,
F_VLT_CONST x (positive). Final value is calculated as f vlt f vlt ,const f vlt ,coeff d ,
172
Units: none
Default value: -0.0035
EPS_M x Maximum compressive strain at the centre of cross section,
(negative)
Units: none
Default value: -0.002
LAMBDA x Coefficient to reduce compressed masonry area
Units: none
Default value: 1.
ETA x Coefficient to apply for F_D
Units: none
Default value: 0.8
REL_TOL x Relative acceptable error in moments/forces
Units: none
Default value: 0.001
ITER_MAX n Maximum number of iterations for zeroizing of lateral bending
moment. Note that the moments are calculated in a coordinate
system, whose Y‘ axis is parallel to the resultant moment from
My and Mz load. Therefore, moment along Z‘ must be equal
zero.
Units: none
Default value:30
EPS_SMALL x Strain value already assumed neglibable
Units: none
Default value: 0.001
ALPHA_STEP x Angle step (for resultant moment load) at which the M-N
diagram of cross section is cached. For zero or negative value
nthing is cached and the appropriate M-N diragram is calculated
on run-time basis.
Units: none
Default value:
60
ALPHA_TOL x Angle difference (for resultatnt moment load) thas is assumed
negligible.
Units: none
Default value:
360
174
Syntax:
& BEAM_RC_MATERIAL :
TYPE “CCBeamRCMaterial” {[E x] | [MU x] | [RHO x] | [ALPHA x] | [F_CK x
] | [F_CVK x] | [F_CTK_INPLANE x ] | [{F_CTK_OUTPLANE}|{F_CTK} x ]
| [GAMMA_M x] | [F_CD x ] | [F_CVD x] | [F_CTD_INPLANE x ] |
[{F_CTD_OUTPLANE}|{F_CTD} x ] |[ EPS_CU x ]|{EPS_C x}|{ LAMBDA x ]
| [ ETA x ] | [ REL_TOL x ] | [ITER_MAX n ] | [ EPS_SMALL x ] | [
ALPHA_STEP x ] | [ ALPHA_TOL x ] | [ FLEX_DRIFT_COEFF x ] |
[SHEAR_DRIFT_COEFF x ] | [ STIRRUPS_SPACING x ] | [STIRRUPS_AREA
x ] | [STIRRUPS_MATERIAL n] | [STIRRUPS_K_I x ] | [STIRRUPS_NI_1 x ] |
[STIRRUPS_EFFECTIVE_DEPTH x ] | [STIRRUPS_C_RD_C x] | [
STIRRUPS_NI_MIN x ] | DAMPING_MASS xM DAMPING_STIFF xK}
ATENA Input File Format 175
Parameter Description
Ex Young modulus.
Units: stresses
Default value: 0
MU x Poisson ratio
Units: none
Default value: 0
RHO x Material density
Units: mass/volume
Default value: 0
ALPHA x Coefficient of thermal expansion
Units: 1/T
cceptable range: <0; maximal real number>
Default value: 0.000012
DAMPING_MASS xM Mass and stiffness damping factors specified for indiviual
DAMPING_STIFF xK element group. They overwrite the same factor set for the whole
structure by SET command .
F_CK x Characteristic material compressive strength, (negative). This
input is not used, if the corresponding design value is given.
Units: stresses
Default value: 0
F_CVK x Characteristic material shear strength, (positive). This input is
not used, if the corresponding design value is given.
Units: stresses
Default value: 0
F_CTK_INPLANE x Characteristic material in-plane tensile strength in bending,
(positive). This input is not used, if the corresponding design
value is given.
Units: stresses
Default value: 0
[{F_CTK_OUTPLANE}| Characteristic material out-of-plane tensile strength in bending,
{F_CTK}x (positive). This input is not used, if the corresponding design
value is given.
Units: stresses
Default value: 0
176
Units: none
Default value: 0.15
STIRRUPS_NI_1 x Coefficient of compressive strut strength. Typically no change is
needed.
Units: none
Default value: based on f ck .
STIRRUPS_EFFECTIV Effective depth of the section, typically distance between the
E_DEPTH x centre of the longitudinal reinforcement and the top edge.
Typically no change is needed.
Units: length
Default value: calculated automatically.
STIRRUPS_C_RD_C x Coefficient based on National annex. Typically no change is
needed.
Units: none
Default value: 0.18
c .
STIRRUPS_NI_MIN x Minimal shear strength. Typically no change is needed.
Default value min 0.035 k 2 f ck2
3 1
Syntax:
& BEAM_REINF_BAR_MATERIAL :
TYPE “CCReinfBarMaterial” {[E x] | [MU x] | [RHO x] | [ALPHA x] | [F_YK x ]
| [F_YVK x] | [GAMMA_M x] | [F_YD x ] | [F_YVD x] | [E_YD_HARD x] |
[EPS_YD_MAX x] DAMPING_MASS xM DAMPING_STIFF xK}
ATENA Input File Format 179
Parameter Description
Ex Young modulus.
Units: stresses
Default value: 0
MU x Poisson ratio
Units: none
Default value: 0
DAMPING_MASS xM Mass and stiffness damping factors specified for indiviual
DAMPING_STIFF xK element group. They overwrite the same factor set for the whole
structure by SET command .
RHO x Material density
Units: mass/volume
Default value: 0
ALPHA_TOL x Angle difference (for resultatnt moment load) thas is assumed
negligible.
Units: none
Default value:
360
F_YK x Characteristic material compressive strength, (negative). This
input is not used, if the corresponding design value is given.
Units: stresses
Default value: 0
F_YVK x Characteristic material shear strength, (positive). This input is
not used, if the corresponding design value is given.
Units: stresses
Default value: 0
GAMMA_M x Partial factor of safety
Units: none
Default value: 1
F_YD x Material strength, (positive)
Units: stresses
Default value: 0
F_YVD x Material shear strength, (positive)
Units: stresses
Default value: 0
180
Sub-Command Description
&LOAD_DISPLACEMENT Prescribed nodal displacement (i.e. Dirichlet boundary
condition), either &SIMPLE_LOAD_DISPLACEMENT, or
&COMPLEX_LOAD_DISPLACEMENT
&LOAD_FORCES Prescribed nodal forces (i.e. Neumann boundary condition),
either &SIMPLE_LOAD_FORCE or
&COMPLEX_LOAD_FORCE
&LOAD_MASTER_SLAVE Master slave node pairs – prescribed displacement as a linear
_NODES combination of other displacements and constant value, (i.e.
Cauchy boundary condition).
&ELEMENT_LOAD Element loads, either &BODY_ELEMENT_LOAD or
&ELEMENT_BOUNDARY_LOAD or
&TEMPERATURE_ELEMENT_LOAD or
&ELEMENT_INITIAL_STRAIN_LOAD or
&ELEMENT_INITIAL_STRESS_LOAD or
&LOAD_FUNCTION or &MASS_ACCELERATIONS or
&ELEMENT_INITIAL_GAP_LOAD or &CHLORIDES or
&CARBONATION or &ASR
&LOAD_FUNCTION Time function id, i.e. id of time (or step id) function defining
coefficient for the applied load. See &FUNCTION for the
function definition.
&SPRING_DEFINITION Spring support boundary condition.
&RIGID_BODY, Definition of rigid body and/or inverse rigid body constrains
&INVERSE_RIGID_BODY
&LOAD_SHELL_TO_SOLI A special boundary condition similar to
D &LOAD_MASTER_SLAVE_NODES that is dsigned to
connect nodes of shell 2D elements to solid elements.
ATENA Input File Format 181
&LOAD_DISPLACEMENT:
SUPPORT [&DISPLACEMENT_TYPE] &LOAD_FUNCTION ]
{&COMPLEX_LOAD_DISPLACEMENT |
&SIMPLE_LOAD_DISPLACEMENT | &SPRING_DEFINITION}+
&DISPLACEMENT_TYPE:
TYPE {DISPLACEMENT | VELOCITY | ACCELERATION}
Note that displacements boundary conditions, (i.e type = “DISPLACEMENT”), are treated as
incremental displacements load, whilst in case of velocities and/or accelerations, (i.e. type =
“VELOCITY” or “ACCELERATION”), the input values are considered to be total load, not
incremental load. Hence, “VELOCITY” and/or “ACCELERATION” BCs (because of its
"total" character) must be specified in the group of "fixed" load within the dynamic load step
definition. On the other hand, “DISPLACEMENT" type BCs are typically input within
"increment" loads of the stepd definition.
&COMPLEX_LOAD_DISPLACEMENT:
{COMPLEX {&MASTER_NODES | &SLAVE_NODES | &LOAD_VALUE |
RELAX}+ }+ [PROCESS_FLAG {REFERENCE_COORDS |
USE_CURRENT_COORDS | COPY_DEFORMATION |
COPY_DEFORMATION_ONCE | COPY_NO_DEFORMATION} ]
182
This type of Dirichlet boundary condition sets the following general boundary condition:
N
ui x u j f j , where i j
j 1
&SIMPLE_LOAD_DISPLACEMENT:
{SIMPLE { SAVE_SUPPORT |REMOVE_SAVED_SUPPORT } { &LOAD_PLACE |
&LOAD_VALUE }+ }+
ATENA Input File Format 183
This type of Dirichlet boundary condition sets the following general boundary condition:
u value
It is the simplest way to define prescribed deformation at a specified node and degree of
freedom (defined in &LOAD_PLACE).
Location of the boundary condition is specified by id of supported node and its supported
degree of freedom. Alternatively, the boundary condition can be set for all nodes (and the
specified supported degree of freedom), whose ids are stored in a list of ids, see command
&SELECTION . In this case, the BC’s value is calculated as follow:
u const x coeff _ x y coeff _ y z coeff _ z , see &LOAD_VALUE command fragment.
In the above x,y,z are coordinates of node id from the list. This way it is possible to
prescribe variable load that depends of coordinates of a node, to which it is applied. Typical
example of such a load may by lateral (hydrostatic) pressure applied to a vertical wall of a
water tank.
If idof >0, then the specified displacement applied at DOF idof . Otherwise current
displacement at -idof is applied.
The optional flags SAVE_SUPPORT |REMOVE_SAVED_SUPPORT are used to save or
remove the current load for later use. It is particularly useful for the case of idof >0, (e.g.
for modelling of a tunnel excavation.
&LOAD_FORCE:
LOAD TYPE {CONCENTRATED_LOAD | LUMPED_MASS }
[&LOAD_FUNCTION] { &COMPLEX_LOAD_FORCE |
&SIMPLE_LOAD_FORCE }+
&COMPLEX_LOAD_FORCE:
{COMPLEX { &SLAVE_NODES | &LOAD_VALUE }+ }+
&SIMPLE_LOAD_FORCE:
{SIMPLE { &LOAD_PLACE | &LOAD_VALUE }+ }+
&LOAD_MASTER_SLAVE_NODES:
{ {MASTER | SOLID_TO_SOLID} { &MS_PAIRS | &MS_GROUPS |
&MS_SELECTION } [ &MS_PROCESS_FLAGS ] }+
&MS_PAIRS:
[SLAVE] [NODAL] PAIRS [ACCEPT_OUTSIDE_ELEMENT] [DISTANCE dx] {
ni [{REPLACE | REPLACES}] ii }+
&MS_GROUPS:
[SLAVE] [NODAL] GROUPS [ACCEPT_OUTSIDE_ELEMENT] [DISTANCE dx]
[SHAPE shape] { { ni }+ { REPLACE | REPLACES } ii }+
&MS_SELECTION:
{ { SELECTIONS | LISTS } list_of_masters list_of_slaves [DISTANCE dx] }
SEARCH_RADIUS r | { LIST_MASTER_NODES list_of_masters
LIST_SLAVE_NODES list_of_slaves [DISTANCE dx] } | {
LIST_SLAVE_NODES list_of_slaves LIST_MASTER_NODES list_of_masters
[DISTANCE dx] }
&MS_PROCESS_FLAGS:
[PROCESS_FLAG {REFERENCE_COORDS | USE_CURRENT_COORDS} |
{COPY_DEFORMATION | COPY_DEFORMATION_ONCE |
COPY_NO_DEFORMATION} ] | [ SKIP_DOFS_MASK skip_mask ] }+
The PROCESS_FLAG input can be used to specify a special way of master-slave boundary
conditions generation. These constrains can be generated using either current or reference
coordinate system. Another option is to copy during the generation displacements from
master points to slave points. It is useful in modeling of construction process. For a
complete description of the PROCESS_FLAG options, see Table 107. skip_mask allows
for definition of DOFs that are skipped, i.e. not connected. If skip_mask is not defined, all
nodal DOFs are linked.
The SKIP_DOFS_MASK skip_mask is used to code, which nodal dofs should be skipped,
i.e. which dofs should not be affected by the current master-slave condition. Displacement
x,y, .. rotation z corresponds to 0b1, 0b10 ...0b100000. For example, let us want to constrain
only displacements x,y and rotation y of nodes with 6 dofs, (3 displacements and three
rotations). Using binary biwise notation, we need to constrain dofs 0b010011. The
skip_mask is complement of 0b010011, i.e. 0b101100. Hence you must input skip_mask
as integer number 44. (0b101100=0x2C=44).
SEARCH_RADIUS r set radius, where to look for surrounding elements to fix the master
slave node.
&LOAD_SHELL_2D_TO_SOLID:
186
&LOAD_BEAM_1D_TO_SOLID:
{ LOAD_BEAM_1D_TO_SOLID LIST_BEAM_NODES " list_of_slaves " |
LIST_SOLID_NODES " list_of_masters " | MAX_DISTANCE x |
SEARCH_RADIUS r | { FIX_ROTATION_BM | FIX_ROTATION_MT |
FIX_ROTATION_BT } | { FIX_ROTATION_LM | FIX_ROTATION_MR |
FIX_ROTATION_LR } | BEAM_GROUP beam_group_id | EMBEDDED
SOLIDS_GROUPS solid_group_ids_interval | ROTATION_ARMS n }+
&LOAD_VALUE:
{ [VALUE value ] | { [CONST const ] | [COEFF_X coeff_x] | [COEFF_Y coeff_y ] |
[COEFF_Z coeff_z ] } }
This command can be used to define a general spatial distribution of loads in the form:
f ( x, y, z ) (const x coeff _ x y coeff _ y z coeff _ z ) value
&SLAVE_NODES
SLAVE { [NODE] ni [DOF] ii }+
&MASTER_NODES
MASTER { [NODE] ni [DOF] ii [*] xi }+
&LOAD_PLACE
{ NODE node | SELECTION “list_name” DOF idof
&LOAD_FUNCTION:
{ [ INCREMENT | TOTAL ] FUNCTION i }2
where cti is load multiplier, ftot (t ), fincr (t ) are values of the total and increment load
functions at time t, ti and ti-1 is time at current and previous step, respectively. The above
formula is applicable for loads that have incremental character. For loads with total
character the load multiplier is calculated by:
cti fincr (ti ) f tot (ti ) .
&ELEMENT_LOAD
LOAD [&LOAD_FUNCTION ] { &LOAD_FUNCTION | [INITIAL]
&BODY_ELEMENT_LOAD | &BOUNDARY_ELEMENT_LOAD |
&TEMPERATURE_ELEMENT_LOAD |
&ELEMENT_INITIAL_STRAIN_LOAD |
&ELEMENT_INITIAL_STRESS_LOAD | &PRESTRESSING_LOAD
&FIXED_PRESTRESSING_LOAD |
&FIXED_PRESTRAINING_LOAD|&MASS_ACCELERATIONS_LOAD |
&ELEMENT_INITIAL_GAP_LOAD | | &DURABILITY_ELEMENT_LOAD }
&LOADED_ELEMS:
[ GROUP { group_id [ TO group_id_to [ BY group_id_by]] } | { SELECTION
list_name } [ ELEMENT { element_id [ TO element_id_to [ BY
element_id_by]] } | { SELECTION list_name } ] ]
&LOAD_COEFF :
ATENA Input File Format 189
&BODY_ELEMENT_LOAD:
BODY [ &LOADED_ELEMS ] [ &LOAD_COEF] [ {LOCAL | GLOBAL}] { { X |
Y | Z | ROT_X | ROT_Y | ROT_Z | DOF idof } [VALUE] x }+ | [
MNODE_IDS "mnode_ids" MNODE_LOADS "mnode_loads"
MNODE_DOFS_MASK mask [MNODE_MAX_DISTANCE abs_max_dist ]
[MNODE_DIM1 n1] [MNODE_DIM2 n2] [MNODE_DIM3 n3] ]
&BOUNDARY_ELEMENT_LOAD:
BOUNDARY [ &LOADED_ELEMS ] [ &LOAD_COEF] [{LOCAL | GLOBAL} ]
[{ANY 6 | SURFACE| EDGE | EDGE_NO_DUPLICATES}] | [ MULTIPLE
{YES|NO} ] | [ NODES “loaded_nodes” ] { { X | Y | Z | ROT_X | ROT_Y |
ROT_Z | DOF idof} [VALUE] x }+ [MERGE [ MERGE_STRING str ]
[NO_ELEM_OUTPUT ] | [ MNODE_IDS "mnode_ids" MNODE_LOADS
"mnode_loads" MNODE_DOFS_MASK mask [MNODE_MAX_DISTANCE
abs_max_dist ] [MNODE_DIM1 n1] [MNODE_DIM2 n2] ]
&TEMPERATURE_ELEMENT_LOAD
TEMPERATURE [ &LOADED_ELEMS ] [ &LOAD_COEF] { REFERENCE
[TIME] t_ref TARGET [TIME] t_target [ IMPORT GEOMETRY
geometry_filename] IMPORT [HISTORY] RESULTS results_filename } | {
[VALUE x] | [REF_VALUE ref_x] | [NODE_ID node_id NODE_VALUE
node_value | REF_NODE_VALUE ref_node_value |
{AUTOMATIC|MANUAL} | TIME_UNITS "time_units"}
&ELEMENT_INITIAL_STRAIN_LOAD:
[INITIAL] STRAIN [ &LOADED_ELEMS ] [ &LOAD_COEF] [IP ip_id] { X | Y |
Z | XY | YX | YZ | ZY | XZ | ZX } [VALUE] x_element_initial_strain }+
&ELEMENT_INITIAL_STRESS_LOAD:
[INITIAL] STRESS [ &LOADED_ELEMS ] [ &LOAD_COEF] [IP ip_id] { X | Y |
Z | XY | YX | YZ | ZY | XZ | ZX } [VALUE] x_element_initial_stress }+
6
The option ANY is only available in 4.3.1 and older; starting 4.3.2, the default is "SURFACE" for 3D
problems and "BOUNDARY" for 2D and axisymmetric problems.
190
&PRESTRESSING_LOAD:
PRESTRESSING [ &LOADED_ELEMS ] [ &LOAD_COEF] [VALUE]
{START_NODE | END_NODE | START_AND_END_NODE } prestres_val
&FIXED_PRESTRESSING_LOAD :
FIXED_PRESTRESSING [ &LOADED_ELEMS ] [ &LOAD_COEF]
[ DIRECTION ] { START_TO_END | END_TO_START } ] { [VALUE |
VALUES] { s_coord value_at_s }+ | VALUE_FNC i }
&MASS_ACCELERATIONS_ELEMENT_LOAD:
MASS_ACCELERATIONS [ &LOADED_ELEMS ] [ &LOAD_COEF] {LOCAL |
GLOBAL} { { X | Y | Z | DOF idof } [VALUE] x }+
&ELEMENT_INITIAL_GAP_LOAD:
[INITIAL] GAP [ &LOADED_ELEMS ] INIT_STEP_ID n
Example:
Use the above command structure to define loads applied to finite element(s).
If element a group and element ids range/selection is specified, then all
(group_id,element_id) combinations are processed except for the case, when loaded
elements are specified by two selections, one with group ids and another with element ids
(of the loaded elements). If ngroup and nelem is number of entries in the appropriate
selections, then in this case min(ngroup , nelem ) elements get loaded.
If no group and element ids range/selection is specified, all elements in the structure are
loaded.
Currently the supported types are:
Volumetric (mass or body) load in a general direction (defined as a vector in
ATENA Input File Format 191
In addition, boundary and body element load can also include a load at arbitrary
location(s). It can take a form of a concentrated load at specific points, or a load
applied along one or more lines, areas and volumes, (the last option is available
only for element body load). Units of the load change respectively, i.e. force,
force/length etc. The control location of the load, i.e. points, lines... need not
coinside with the finite element mesh and can overlap several elements, (and/or an
element can be loaded by more such loads). The load must be applied to the above
defined element edge or surface, (the case of boudary load) or to a body, (the case of
body load). It is input as MNODE_IDS "mnode_ids" MNODE_LOADS
"mnode_loads" MNODE_DOFS_MASK mask [MNODE_MAX_DISTANCE
abs_max_dist ] [MNODE_DIM1 n1] [MNODE_DIM2 n2] [MNODE_DIM3 n3]
], where "mnode_ids" is selection with macro node ids, where the concentrated
load is applied, "mnode_loads" is a selection containing corresponding loads, (i.e.
real numbers), defined by SELECTION_REAL command. The loads are sorted by
loaded dofs and mnode_ids. mask is bitwise map indicating, which dofs are to be
loaded.. For example mask=6 means loading of dir_y(=bit 2=2) + dir_z(=bit3=4),
i.e. "mnode_loads" must comprise FY, FZ at mnode_1, followed by FY, FZ at
mnode_2 etc. The parameter abs_max_dist sets an accuracy, with which ATENA
checks, if a macro node is or is not located at the processed element surface. In the
latter case the surface is not loaded. Finaly, the parameters n1, n2, n3 specify type of
loaded loacation. For n1=n2=n3=1 the load is applied to points.This is default
behaviour. If just one of nx>1, then the load is applied to a line. If just two of nx>1,
then the load is applied to an quadrilateral and similarly if all three of nx>1, then the
7
The option ANY is only available in 4.3.1 and older
192
load is applied to a hexahedral, (applicable for element body load only). Layout of
the controle load objects resembles the layout used for finite element of the same
shape and the actual value of nx specify number of points, into which the load
should be distributed, (n1,n2,n1 for directions in isoparametric r,s,t coordinates).
The MERGE flag is used, if the current boundary load should be merged with a
previous boundary load within the same load case. MERGE_STRING str allows
merging only boundary loads with the same MERGE_STRING str. The merging is
successful, if the current and the other boundary load are of the same type,
(edge/surface) and have the same values. Other parameters, (e.g. function_id,
coeff_x etc.) are not tested and values from the other boundary load are adopted. If
the merging is not successful, then the current boundary load is processed in the
same way as it would without the MERGE flag. The NO_ELEM_OUTPUT flag
suppress element boundary related output at element level. Note that only single
element surface or edge can be loaded within single boundary load. Hence, use
MERGE option with caution.
TYPE_STRING str is used only for output data aggregation.
Element temperature load, &TEMPERATURE_ELEMENT_LOAD that
corresponds to element initial strain load, where initial strains are calculated based
on material expansion coefficient and specified temperature. The temperature
history can also be imported from the associated CCStructuresTransport analysis. In
this case one has to input IMPORT subcommand. If results_file_name is specified
without geometry_filename_name , it means that imported and current models are
identical. If geometry_filename_name is specified, an interpolation between the two
models is executed. Note that the IMPORT HISTORY option should be used only,
if target and reference times are given, (see REFERENCE [TIME] t_ref | TARGET
[TIME] t_target. This is because any loading in ATENA is assumed to be of
incremental character. Hence, the TEMPERATURE_LOAD is imported as
temperature increments betveen the structural conditions at target and reference
time. Alternatively, temperature load increments at element nodes can be input
directly using syntax { NODE_ID node_id NODE_VALUE node_value }. Note
that element node related input is always added to average element temperature
load, see [VALUE] x. Some material laws are temperature denpent and thus they
need info about absolute temperatures, rather then temperature increments (used e.g.
for element load due the material thermal expansion). These are input thru
REF_VALUE ref_x and REF_NODE_VALUE ref_node_value in the similar way
as temperature increments via VALUE x and NODE_VALUE node_value. Note
that from the transport analysis, i.e. using the IMPORT command, they are
imported automatically. The reference temperatures ignores any load coefficient
coming from function definition, load case multiplier etc. The AUTOMATIC option
causes Atena to automatically update TARGET and REFERENCE TIME according
to time at the current and previous step. It is usefull particularly for element
tremperature load during creep analysis. If AUTOMATIC, the load is imported
from history files and no additional load is acceptable, (such as via VALUE and
NODE_VALUE). By default, MANUAL regime is assumed. }. The TIME_UNITS
"time_units" allows to specify, which time units were used to calculate and write the
transpored analysis results in the file results_file_name. It is specified in the same
way as in the Unit command. By default no time unit conversion is made.
ATENA Input File Format 193
CHLORIDES, CARBONATION and ASR element load does not represent a real
load and are described in Section 0. It merely forces Atena to calculate degradation
of reinforced concrete elements due to progression of carbonation and/or chlorides
from their outside surfaces.
function), it is assumed that this function defines total accelerations in a time (and
not load increments, as it is usual in most other load types). The corresponding load
increment at time t t is then calculated as a ( f t t ) f (t ) , where f (t ) is the
acceleration function and a is constant acceleration in a particular direction having
been input within this load specification. This load is meanigful in dynamic analysis
only and because of its "total" character, it must be specified in the group of "fixed"
load within the dynamic load step definition, (i.e. not among "increment" loads!
The element load is aplied to element groups specified by GROUP group_id [ TO
group_id_to [ BY group_id_by]] command tokens. Otherwise all element groups are
loaded. For each element group it is possible to load only some elements. Their list is input
in ELEMENT SELECTION list_name command tokens. If the list contains a non-existing
element, the corresponding entry is ignored. Alternatively, the loaded elements can be input
in form of interval ELEMENT element_id [ TO element_id_to [ BY element_id_by]]. In
this case, however, one have to be cautious. element_id [ TO element_id_to must exist in
the group group_id. For the remaining element groups, i.e. up to group_id_to, internal
element numbering is used. E.g. let group group_id has elements 100, 105, 108, 110, 120,
130 and element_id=105, element_id_to=110. Then the remaining loaded element groups
(i.e. groups up to group_id_to) receive the load into their second, third and forth element.
(The elements within each group are sorted according to their element_id). As usuallly, by
default all elements of the group are loaded.
In addition, it is possible to use linear spatial interpolation based on the element’s centre-
point coordinates and COEFF_X coeff_x] [COEFF_Y coeff_y ] [COEFF_Z coeff_z ] see
Table 113. By default, coeff_x=0, coeff_y=0, coeff_z=0 and const=1. If only GROUP
group_id is given (and [ELEMENT element_id] is omitted), then the load applies to all
element of the specified element group. An exception to that is prestressing of external
cable. This load is always applied in element_id=1 (and only once, if element_id is not
specified).
Different values of element initial stress and strain can be applied at each material (i.e.
integration) point, see IP ip_id input. If ip_id=0, the element load is applied into all material
points. Hence, with ip_id=0 the user can specify “uniform” portion of a load (across the
element) and then he can define the load deviation at a particular material point ip_id.
By default ip_id=0.
&DURABILITY_ELEMENT_LOAD:
[ &CARBONATION ] | [ &CHLORIDES ] | [ &ASR]
&CARBONATION 8:
CARBONATION { TIME_FNC_ID id | WATER_MASS x | CEMENT_MASS x |
SCM_MASS x | CONCRETE_COVER x | K_CO2 x | CO2 x | RH x | NODES
“loaded_nodes” ” | A1 | A2 | A3 | F_T_CH | W_D | B | R_CORR |
DX_CORR_DT_SPALLING | CEMENT_MASS | MAX_REINF_DEPTH |
8
development/testing implementation of CARBONATION, CHLORIDES, and ASR in version 5.3.x and older
ATENA Input File Format 195
&CHLORIDES 9:
CHLORIDES { TIME_FNC_ID id | D_REF x | TIME_D_REF x | M_COEFF x |
TIME_M_COEFF x | CONCRETE_COVER x | CS x | CL_CRIT x | NODES
“loaded_nodes” | A1 | A2 | A3 | F_T_CH | W_D | B | R_CORR |
DX_CORR_DT_SPALLING | CEMENT_MASS | MAX_REINF_DEPTH |
T_AVER_OFFSET | MAX_CORR_DT | PSI }+ [TYPE_STRING str ] [MERGE
[ MERGE_STRING str ]] [NO_ELEM_OUTPUT ] }+
&ASR 10:
ASR { TIME_FNC_ID id | U_C x | U_L x | TAU_C_0 x | TAU_L_0 x | THETA_0 x |
SAND_MASS x | REQUIRED_ALKALI_PER_REACTIVE_SILICA x |
PROPORTION_REACTIVE_SILICA x |
PROPORTION_REACTIVE_PARTICLES_IN_SAND x |
ASR_MEASSURED_ASR_STRAIN x | THRESHOLD_ALKALI_IN_CONCRETE x |
TOTAL_ALKALI_IN_MORTAR x | STEP_TIME_INCR x | MAX_NUMBER_OF_ITERS n
| T_AVER_OFFSET x | S_L x | S_U x | EPS_CR x | EPS_U x | BETA_E x | BETA_FT x |
BETA_GF x | H_MIN | H_AVER | BETA_SHAPE_FACTOR | ASR_MODE {
APPLY_NOTHING | ASR_MODE APPLY_INIT_STRAINS | APPLY_REDUCTION_ALL
| APPLY_INIT_STRESSES}+ }+
This command extends the ELEMENT_LOAD command specifically for various durability
analysis options. Use the above command structure to define loads applied to finite
element(s) for durability load types.
CHLORIDES, CARBONATION and ASR element load does not represent a real load. It
merely forces Atena to calculate degradation of reinforced concrete elements due to
progression of carbonation and/or chlorides from their outside surfaces. The input data
resembles &BODY_ELEMENT_LOAD. It applies to the parameters NODES
“loaded_nodes” “loaded_nodes”, MERGE, MERGE_STRING str and
NO_ELEM_OUTPUT. The remaining parameters are:
- WATER_MASS, CEMENT_MASS and SCM_MASS - mass of water, cement and non-
active suplementary cementitious material, SCM per 1m3, [weight/volume], default
200kg/m3, 400 kg/m3, 50kg/m3, respectively,
- CONCRETE_COVER : thickness of concrete cover layer, [length], default value 0.02m,
- K_CO2 : efficiency factor, [-], with typical values 0.3 for silica fume, 0.5 for low- calcium
fly ash, 0.7 for high-calcium fly ash, effective only for concrete with SCM_MASS >0, i.e.
not for Portland cement, default value 0.5,
- CO2 : content CO2 in the ambient air, [-], default 0.00036,
- RH : relative humidity of ambient air RH, [-], default 0.6,
- CL_CRIT : critical mass of chlorides per mass of SCM+cement for initialisation of
9
development/testing implementation of CARBONATION, CHLORIDES, and ASR in version 5.3.x and older
10
development/testing implementation of CARBONATION, CHLORIDES, and ASR in version 5.3.x and older
196
Parameter Description
DIRECTION { x }ncoords Spring direction.
E.g. DIRECTION x1 x2 [x3]
Component x3 is valid only in 3D problems. Positive internal
spring force acts in the direction given by this vector.
NODE n Node number, in which the spring is applied.
MATERIAL n Spring stiffness material id.
198
Parameter Description
ID n Load case identification.
NAME “load case name” Load case name in quotes, also for identification.
E.g.: NAME “load case name”
MASTER { [NODE] ni List of master nodes, their degrees of freedom and
[DOF] ii [*] xi }+ multipliers.
E.g.:
MASTER NODE n1 DOF i1 * f1 .... NODE nk DOF ik * fk
SLAVE { [NODE] ni [DOF] List of slave nodes and their degrees of freedom. They are
ii }+ ordered according to MASTER
E.g.:
SLAVE NODE n1 DOF d1 .... NODE nk DOF dk
VALUE x Prescribed nodal value, either displacement or force
depending on context.
E.g.: VALUE x
MASTER [SLAVE] Ids of master-slave nodal pairs.
[NODAL] [PAIRS] { ni ii }+ E.g.: MASTER [SLAVE] [NODAL] [PAIRS] n i , n i ,
1 1 2 2
n3 i3 ... ni ii
NODE n DOF n Node and its DOF specifying a place, where the simple
boundary condition is applied.
FUNCTION n Id of time function applied atop of a specified boundary
condition.
E.g.: FUNCTION n
{ X | Y | Z | DOF idof } Element body load components in reference coordinate
[VALUE] system, (in force per volume unit). If DOF idof is used, the
specified value applies to a DOF idof .
E.g. X [VALUE] x Y [VALUE] x Z [VALUE] x
TEMPERATURE Element temperature, (in deg).
STRAIN { X | Y | Z | XY | Component of element initial strain components in
YX | YZ | ZY | XZ | ZX } reference coordinates system.
[VALUE]
STRESS { X | Y | Z | XY | Component of element initial stress components in
YX | YZ | ZY | XZ | ZX } reference coordinates system.
[VALUE]
GROUP, ELEMENT Group and element ids, where the ELEMENT_LOAD is
applied.
ATENA Input File Format 199
&INVERSE_RIGID_BODY
INVERSE_RIGID_BODY SLAVE_ID n MASTER_SELECTION list_of_masters
FIX_DOFS dofs_mask MASTER_WEIGHTS (w1,w2…)
&BEAM_NL_CONNECTION
BEAM_NL_CONNECTION LIST_OF_NODES list_of_nodes SKIP_DOFS_MASK
skip_mask MAX_COORDS_TOL max_tol
Parameter Description
ID n1 [ TO n2 [ BY n3 ] ]
Steps interval that would be executed by EXECUTE
subcommand. By default n3=1, n2= n1.
&STEP_TYPE_AND_DA Type and data for a particular load step. Currently STATIC,
TA EXECUTE TRANSIENT, CREEP and DYNAMIC type are available.
&STATIC_STEP_DEFINITION:
[TYPE] STATIC { NAME “step name“ | ID n }+ {[LOAD] [CASE] ni * xi }+
Parameter Description
STATIC Static load step.
NAME “step name“ Step name in quotes that is going to be defined.
ID a Integral identification of the step “step name“.
[LOAD] [CASE] { ni * xi }+ Linear combination of load cases for step “step name”,
which are to be used in this step.
E.g.: LOAD CASE 1 * 1.5 2 * 0.8
&TRANSIENT_STEP_DEFINITION:
[TYPE] TRANSIENT { NAME “step name“ | ID n }+ {[LOAD] [CASE] ni * xi }+
Parameter Description
TRANSIENT Transport analysis load step.
ATENA Input File Format 201
&CREEP_STEP_DEFINITION
TYPE CREEP { NAME “step name“ | ID n | {AT|RESUME_AT} time | [ {FIXED |
INCREMENT } ] [LOAD] [CASE] ni * xi }+
Parameter Description
TYPE CREEP Creep load step. As creep analysis involve numerical time
integration, the creep step consists typically of several
“static like” integration steps, one for each sample time. It
starts at creep step time of the current creep step and stops
at min(time of the next creep step, execution_stop_time)
(see &CREEP_ANALYSIS_PARAMS.) The analysis
cannot exceed time_end, see &RETARDATION.
NAME “step name“ Step name in quotes that is going to be defined.
ID a Integral identification of the step “step name“.
{AT|RESUME_AT} time Time at the beginning of the current creep step, in [days]. If
"AT" label is used, ATENA assumes that an additional
loading is applied in this step and therefore it automatically
refines time integration, (i.e. it resets step time incerements
dt to 0.1 days). If "RESUME_AT" label is used, no
additional loading is assumed and thus, no special time
refinement is carried out. This option can be used for
getting user control and produce some print outs, figures
etc. during execution of creep analyses.
[LOAD] [CASE] | Linear combination of load cases for step “step name”,
[{FIXED | INCREMENT }] which are to be used in this step. The FIXED type of load is
{ ni * xi }+ evenly distributed into all applied integration time sub-steps
of the current creep step, whilst the INCREMENT type is
used only in the 1st integration sub-step. In the remaining
sub-steps they are applied, but load values are a priori
zeroised. Typically loads are specified as of INCREMENT
type and LHS boundary conditions as of FIXED type. By
default the FIXED type is assumed.
E.g.: LOAD CASE FIXED 1 * 1.5 2 * 0.8 INCRENENT
3 * 1.3 4 * 10.8
&DYNAMIC_STEP_DEFINITION
202
Parameter Description
TYPE DYNAMIC Dynamic analysis related load step. As dynamic analysis
involve numerical time integration, the dynamic step
consists typically of several “static like” integration steps,
one for each sample time. It starts at time of the current step
and stops at min( step time of the next dynamic step,
execution_stop_time). It behaves similarly to creep analysis,
however, dynamic analysis uses equal size sub-step time
lenghts.
NAME “step name“ Step name in quotes that is going to be defined.
ID a Integral identification of the step “step name“.
AT time Time at the beginning of the current dynamic step, in
[days].
If the step's id is defined in form of an interval, the value of
time´ is incremente based on current time increment dt.
[LOAD] [CASE] | Linear combination of load cases for step “step name”,
[{FIXED | INCREMENT }] which are to be used in this step. The FIXED type of load is
{ ni * xi }+ evenly distributed into all applied integration time sub-steps
of the current dynamic step, whilst the INCREMENT type
is used only in the 1st integration sub-step. In the remaining
sub-steps they are applied, but load values are a priori
zeroised. Typically loads are specified as of INCREMENT
type and LHS boundary conditions as of FIXED type. By
default the FIXED type is assumed.
E.g.: LOAD CASE FIXED 1 * 1.5 2 * 0.8 INCRENENT
3 * 1.3 4 * 10.8
&EXPORT_IMPORT | &SUPLEMENT_MONITOR
[PRESERVE_OUTPUT_OPTIONS] | | REMOVE | { FILE_OPEN | FILE_CLOSE }
“file_name” “ [ { overwite | append } ] | [MAXIMUM |
MINIMUM|SUMMATION|AVERAGE]
[RECORD] LENGTH x | &LOCATION | | {TRACK | RECORD} &DATA |
TRACE { OFF | ON } | RECOVERY {LUMPED | VARIATIONAL|
NEAREST_IP } }+
&OUTPUT_TYPE:
{ STANDARD | { MONITOR | MONITOR_1 | MONITOR_2 | MONITORS| PLOT |
PLOT_1 | PLOT_2 } [ EACH {ITERATION | STEP} ] }
&EXPORT_IMPORT:
{ INTERPOLATE { FULL | NONE | STEP } EXPORT {DATA|CMDS} TO
"filename" | IMPORT {DATA|CMDS} FROM "filename_1”, "filename_2”…
"filename_n” }
& SUPLEMENT_MONITOR :
SUPLEMENT FROM n ARCHIVES "filename_1”, "filename_2”… "filename_n”"
&FORMAT_CSV:
{FORMAT__CSV_ON | FORMAT_CSV_OFF}
®EX:
{REGEX_ON | REGEX_OFF}
&LOCATION:
LOCATION { ELEMENT_IPS | ELEMENT_ NODES | NODES | GLOBAL |
LOAD_CASES | ELEMENT_TYPES | MATERIALS | GEOMETRIES |
OUTPUT_DATA } &LOCATION_LIST
&LOCATION_LIST:
{{ GROUP[S] &INTERVAL [ ELEMENT[S] &INTERVAL [ IP[S] &INTERVAL ] ]
| GROUP[S] &INTERVAL [ ELEMENT[S] &INTERVAL [ ENODE[S]
&INTERVAL ] ] | NODE[S] &INTERVAL | ID[S] &INTERVAL | LOC_1
&INTERVAL [ LOC_2 &INTERVAL [ LOC_3 &INTERVAL ] ] }+ }|{
MULTI_SELECTION AT [SELECTION] multi_selection_list }
&INTERVAL:
{ AT { n | FROM n [TO n [BY n]] | SELECTION selection_list} }+
&DATA:
DATA { ALL | { ITEM n [TO n [BY n]] LIST {“output_keyword” [RECALCULATE]
{ AT n1 FROM ITEM n1 [TO n2 [BY n3]] }+ END }+ }
Parameter Description
204
OUTPUT | PRINT Two output channels are supported: the primary and secondary.
The primary channel is activated using the keyword OUTPUT and
it is managed in the same way as all the other standard I/O streams
in Atena. The secondary channel is activated by use of the keyword
PRINT. It is aimed for exporting Atena output to a third party
application via a disc file. Both the output channels can operate
simultaneously as required.
MONITOR | Adds output set “set_name“ into monitor output requests.
MONITOR_1 | Output format is set to produce output data records versus time, in
MONITOR_2 | which all output data (for a particular step or iteration, i.e. for a
MONITORS EACH particular time) are written into one line. The first word of such line
{ITERATION | contains “set_name“, followed by current step id , iteration id and
STEP } time, and then all output items are sequentially printed one after
another. Use “grep set_name” or similar to extract output lines
corresponding to “set_name” output data for their import into a
thirty-party post-processing package like spreadsheets etc.
The specified output command is processed after completing of
every iteration or step.
If the keyword MONITOR is specified, the MONITOR_1 set is
used.
Two output sets are available, one called MONITOR_1 and the
other MONITOR_2. Both of them can be used for monitoring
output data per iteration or per step, however, it is not
recommended to mix ouput monitors per iteration with monitors
per step into the same monitor set. (It would result in a table with
data delivered by iterations with empty slot for data monitored per
step, when convergence was not reached yet.). Hence, one of the
monitors is typically used for monitoring output at each iteration
and the other for output at each step.
Two output sets are particularly useful, if AtenaWin/ Atena Studio
is used for execution of the ATENA analysis. This is because
AtenaWin/AtenaStudio can directly plot all the data from the
monitors into 2D plots without need of any thirty-party SW.
However, in this case it is recommended to use the set
MONITOR_1 for output monitors per iteration and the set
MONITOR_2 for monitors per step, because AtenaWin /
AtenaStudio automatically allocates a monitor with information
about analysis convergence called "ConvergenceMonitor" into the
set MONITOR_1 and it produces convergence information per
iteration. The monitor MONITOR_1 is thus pre-selected" to output
monitors per iteration and MONITOR_2 remains free for step
monitors.
The option “MONITORS” is used for export/import data or
command from/to the both monitors, i.e. it operates on both sets
MONITOR_1 and MONITOR_2. It has nothing to do with
definition of a particular output data monitoring.
ATENA Input File Format 205
PLOT | PLOT_1 | The way of using the keywords PLOT | PLOT_1 | PLOT_2 is
PLOT_2 EACH nearly the same as the use of the keyword MONITOR |
{ITERATION | MONITOR_1 | MONITOR_2. When specified, it (also) creates a
STEP } set of data that can be printed or drawn in 2D plots. The following
table points out the differences:
UNSPLIT_MONIT "NODAL_DISPLACEMENT_AT_NODE_1",
OR_DATA_BY_LO "NODAL_DISPLACEMENT_AT_NODE_2" ...
CATION "NODAL_DISPLACEMENT_AT_NODE_n".
ELEMENT_NODE, ELEMENT_IPS AND ELEMENT Location's
data are splitted by elements, e.g.
"FORCES_AT_GROUP_20_ELEMENT_4". (The level 3 is not
accounted for).
STANDARD Output format is set to “table” oriented form, i.e. items are printed
in separate tables. Each line of such a table presents results for one
location.
Output command request is processed immediately after its issuing.
NAME “set_name” Name of monitor output set.
INTERPOLATE { Export/Import data from/to specified monitors. The “export” is
FULL | NONE | always for the current step, i.e. time. The import is for time saved
STEP } EXPORT in import archives. When importing, linear interpolation of
{DATA|CMDS} TO monitored output data can be requested. If “INTERPOLATE
"filename" | STEP” is specified, the imported output data are smoothly
IMPORT connected to the data from the recent step. If “INTERPOLATE
DATA|CMDS} FULL” is input, the imported data get connected to the lastly
FROM entered value, e.g. typically value for a last previous step, where
"filename_1”, the data were monitored for the last time. “INTERPOLATE
"filename_2”… NONE” suppresses any interpolation. "filename" is binary file into
"filename_n” which the data are exported. "filename_1”, "filename_2”…
"filename_n” are filenames of previously exported data that should
be now imported.
The “DATA” and “CMDS” options are used to export/import
actual output data/monitor output command requests.
SUPLEMENT Force Atena to automatically add the output data history into the
FROM n both monitors, (regardless of MONITOR_1/MONITOR_2 option).
ARCHIVES For each of the specified archive files it restores that file, (i.e.
"filename_1”, state), executes current output monitor requests and exports all
"filename_2”… results. After that, it restores back the current state and imports all
"filename_n” the exported data, thereby adding output data history, (i.e.
[PRESERVE_OUTP monitors) from the specified archives. This command is useful, if at
UT_OPTIONS] a later time it is needed to add some monitored data from previous
times, (i.e. from previous archives).
PRESERVE_OUTPUT_OPTIONS causes to use for the
supplemented monitor data current settings of the output data
conditions, (such as recovery type etc.) rather then the settings,
which were in use during the original execution.
REMOVE Removes output set “set_name” from monitor output requests.
{ FILE_OPEN | This command specifies the name of the output file to be opened or
FILE_CLOSE } closed. All output following this command will be redirected to this
ATENA Input File Format 207
“file_name” [ { file or the default stream. The open file is overwrited or appended.
overwite | append } ]
TRACE {OFF | ON} Flag for tracing results during iterations. By default, data (e.g. at
element IPs) can be traced even during iterations; (either by
OUTPUT MONITOR EACH ITERATION ... or from ATENA
GUI). As this extra output service costs not-negligible resources
(mainly RAM), the user may find reasonable to switch off this
service in case of extensive analyses (e.g. at areas being not critical
for structural over-all behavior). This output is available only for
the location ELEMENTS.
RECOVERY Method for recovering output data akin stress, strain etc. from IPs
{LUMPED | to element nodes. It can be either VARIATIONAL, in which case
VARIATIONAL | an energy based is used to do the recovery, or a simplified
NEAREST_IP } LUMPED method. The former one is more accurate and
theoretically thorough, however it is also more costly in terms of
CPU requirements. By default, the LUMPED approach is used.
Note that LUMPED is usually preferred for linear elements, whilst
VARIATIONAL is the best choice for nonlinear elements. The
third option, i.e. NEAREST_IP set values in element nodes to be
equal to those at the nearest integration point.
This output is available only for the location ELEMENTS.
MAXIMUM | Output only maximum, minimum, sum or average of all values
MINIMUM|SUMM over the printed domain, incl. loop over specified data
ATION|AVERAGE items(=components).
{TRACK | This flag is significant only for MONITOR output. If TRACK is
RECORD} used, the monitored output data are stored for later output and they
are also printed immediately. The keyword RECORD inhibits the
immediate output and the data are only stored for later use.
Default value: TRACK
Table 128: Output-type keywords understood by the command &OUTPUT for the location
type OUTPUT_DATA
Table 129: Output-type keywords understood by the command &OUTPUT for the location
type GLOBAL
Table 130: Output-type keywords understood by the command &OUTPUT for the location
type LOAD_CASES
Table 131: Output-type keywords understood by the command &OUTPUT for the location
type ELEMENTS
Table 132: Output-type keywords understood by the command &OUTPUT for the location
type ELEMENT_IPS
Table 133: Output-type keywords understood by the command &OUTPUT for the location
type ELEMENT_NODES
Table 134: Output-type keywords understood by the command &OUTPUT for the location
type NODES
Table 135: Output-type keywords understood by the command &OUTPUT for the location
type GEOMETRIES
Table 136: Output-type keywords understood by the command &OUTPUT for the location
type ELEMENT_TYPES
Table 137: Output-type keywords understood by the command &OUTPUT for the location
type MATERIALS
Table 138: Output-type keywords understood by the command &OUTPUT for the location
type MACRO_ELEMENTS
Table 139: Output-type keywords understood by the command &OUTPUT for the location
type MACRO_NODES
Examples:
OUTPUT LOCATION OUTPUT_DATA DATA LIST "CURRENT_SORTED_LHS_BC"
END
OUTPUT NAME "displ" MONITOR_1 EACH ITERATION
LOCATION NODES NODE AT 132 DATA LIST "DISPLACEMENTS" ITEM AT 3
END
OUTPUT NAME "s_coord" PLOT_2 LOCATION NODES NODE AT SELECTION
"border_nodes" DATA LIST "REFERENCE_BORDER_COORDINATE" END
ITEM FROM 1 TO 1 ;
Parameter Description
results_filename Name of binary file with the history. It must be the same
as that specified for HISTORY EXPORT command in the
CCStructuresTransport module. It should be enclosed in
double quote character (“).
geometry_filename Name of binary file with geometry of the imported model.
It must be the same as that specified for HISTORY
EXPORT command in the CCStructuresTransport module.
It should be enclosed in double quote character (“). If
omitted, identical imported and current models are
ATENA Input File Format 219
assumed.
num_int_hum Number of intervals into which nodal humidities at each
time step should be sorted. By default num_int_hum=1.
num_int_temp Number of intervals into which nodal temperatures at each
time step should be sorted. By default num_int_temp=1.
HUMIDITY_ABS_MAX_ER Relative and absolute humidity and temperature “errors”
ROR err1 | that are considered as negligible. The values are used
HUMIDITY_REL_MAX_ERR during mapping of moisture and humidity histories at
OR err2 | structural material points. If the tested and master values
TEMPERATURE_ABS_MAX differ less than as it is required by these maximum
_ERROR err3 | “errors”, than no new history is created and the tested
material point is mapped towards the master material
TEMPERATURE_REL_MAX point. By default, these “errors” are set to 0.1.
_ERROR err4
TIME_UNITS "time_units" The TIME_UNITS "time_units" allows to specify, which
time units were used to calculate and write the transpored
analysis results in the file results_file_name. It is specified
in the same way as in the Unit command.
By default no time unit conversion is made.
Syntax:
&DYNAMIC_INITIAL_CONDITIONS:
NODAL {ACCEL_VEL | VEL_ACCEL | ACCELERATION | VELOCITY }
[SETTINGS] { &MANUAL_INITIAL_VALUES_ENTRY |
&GENERATED_INITIAL_VALUES }
&MANUAL_INITIAL_VALUES_ENTRY:
{ NODE n VEL vel_x vel_y [vel_z] | ACCEL accel_x accel_y [accel_z] }
Sub-Command Description
NODE n Set initial conditions for node n.
VEL vel_x vel_y [vel_z] Specify initial nodal velocities in direction of global
coordinates. 3D problems need 3 values, 2D problems only
two values..
ACCEL accel_x accel_y Input initial nodal acceleration in similar way as the above
[accel_z] initial velocities input.
&GENERATED_INITIAL_VALUES:
Sub-Command Description
SELECTION Name of selection, for which the generation is requested.
"selection_name"
{GENERATE_ACCEL | Keyword for entities to be generated. The values in global
GENERATE_VEL } structural directions are generated as linear combination:
CONST const_vector
COEFF_X coeff_x_vector |
COEFF_Y coeff_y_vector |
ATENA Input File Format 221
COEFF_Z coeff_z_vecor valuex const (1) x coeff x (1) y coeff y (1) z coeff z (1)
EQN "eqn" valuey const (2) x coeff x (2) y coeff y (2) z coeff z (2)
valuez const (3) x coeff x (3) y coeff y (3) z coeff z (3)
Example:
NODAL VEL_ACCEL SETTING NODE 1 VEL 0.0030 0. 0. ACCEL -0.005370861556 0.
0.
NODAL VEL_ACCEL SELECTION “my_selection”
CONST 0.0030 0. 0. COEFF_X 0. 0. 0. COEFF_Y -0.6523648649 0. 0. COEFF_Z 0. 0. 0.1
GENERATE_VEL
CONST -0.005370861556 0. 0. COEFF_X 0. 0. 0. COEFF_Y 0. 0.1 0. COEFF_Z 0. 0. 0.
GENERATE_ ACCEL.
Parameter Description
&TRANSIENT Set transient analysis and set some parameters for it.
Syntax:
&TRANSIENT:
TRANSIENT { [TIME] CURRENT x | [TIME] INCREMENT x | STOP_TIME
execution_stop_time | LAST_TIME last_time | NEWMARK_METHOD |
HUGHES_ALPHA_METHOD | } | NEWMARK BETA x |
NEWMARK_GAMMA x | HUGHES_ALPHA x | DAMPING STIFFNESS
[COEFFICIENT] x | DAMPING MASS [COEFFICIENT] x }+
222
Parameter Description
[TIME] CURRENT x Sets current time.
[TIME] INCREMENT x Sets time increment in steps.
STOP_TIME Time at which the execution should stop.
execution_stop_time
There are few new SET &EIGENVALUES parameters as described below, see &SET,
subparameter &ANALYSIS_TYPE
Parameter Description
&EIGENVALUES Set some parametyers for eigenvalues analysis.
Syntax:
& EIGENVALUES:
{ NUMBER_OF_EIGENVALS n | MAX_EIGENVAL_ERROR r |
MAX_NUMBER_OF_SSPACE_ITERATIONS n |
REQUEST_STURM_SEQUENCE_CHECK {YES | NO} |
MAX_NUMBER_OF_JACOBI_ITERATIONS n | NUMBER_OF_PROJ_VECS
n | SHIFT_EIGENVALUES shift, }+
Parameter Description
NUMBER_OF_EIGENV Sets number of the lowest eigenmodes that should be calculated.
ALS n Default value: 10
MAX_EIGENVAL_ERR Maximum eigenvalues error that is tolerated.
OR r Default value: 1.E-6
MAX_NUMBER_OF_S Max. number of subspace iterations.
SPACE_ITERATIONS n Default value: 16
STURM_SEQUENCE_C Flag for requesting Sturm check that no eigenvalue got missed
HECK {YES | NO} during the solution. This check is supported only by the direct
skyline solver. Using of a sparse matrix solver will turn down
eventual request for the Sturm check.
MAX_NUMBER_OF_J Max. number of iteration within Jacobi. The Jacobi procedure
ACOBI_ITERATIONS n computes eigenmodes of the projected global eigenvalues
problem via minimization of Rayleigh quotient. Hence, within
each (“main”) iteration of inverse subspace iteration method
another iterating process is executed in Jacobi. The value of n
sets maximum number of these iterations that are allowed.
Default value: 12
NUMBER_OF_PROJ_V Defines number of projection vector used by Rayleigh quotient
ECS n method. It must be equal or bigger than the number of required
eigenvalues.
Default value: min(2*n_eigenvals, eigenvals+8), where
n_eigenvals is the number of required eigenvalues.
224
//----------------------------------------------------------------------------\\
// Material definition \\
ATENA Input File Format 225
//----------------------------------------------------------------------------\\
MATERIAL
id 1
name "Spring"
type "CC3DElastIsotropic"
E 30
Mu 0.00
Rho 0.000000000001
Alpha 1.200E-05
MATERIAL
id 2
name "Spring"
type "CC3DElastIsotropic"
E 30000000
Mu 0.00
Rho 156.
Alpha 1.200E-05
//----------------------------------------------------------------------------\\
// Element type definition \\
//----------------------------------------------------------------------------\\
ELEMENT TYPE
id 1
name "1D Truss"
type "CCAhmadElement33L9"
//----------------------------------------------------------------------------\\
// Geometry definition \\
//----------------------------------------------------------------------------\\
//----------------------------------------------------------------------------\\
// Joint coordinates definition \\
//----------------------------------------------------------------------------\\
JOINT COORDINATES
//----------------------------------------------------------------------------\\
// Element group definition \\
//----------------------------------------------------------------------------\\
ELEMENT GROUP
id 1
name "Spring"
type 1
material 1
geometry 1
ELEMENT INCIDENCES
1 1 13 15 3 6 18 20 8 9 14 10 2 11 19 12 7 4 16 17
5
2
13 25 27 15 18 30 32 20 21 26 22 14 23 31 24 19 16 28
29 17
3 25 37 39 27 30 42 44 32 33 38 34 26 35 43 36 31 28 40
41 29
ELEMENT GROUP
id 2
name "Mass"
type 1
material 2
geometry 1
ELEMENT INCIDENCES
1
37 49 51 39 42 54 56 44 45 50 46 38 47 55 48 43 40 52
53 41
LOAD CASE
id 2
name "Concetrated force"
LOAD SIMPLE node 49 dof 1 value 0.25
LOAD SIMPLE node 51 dof 1 value 0.25
LOAD SIMPLE node 54 dof 1 value 0.25
LOAD SIMPLE node 56 dof 1 value 0.25
NODAL SETTING
node 49 vel 0.0030 0. 0. accel -0.005370861556 0. 0.
node 50 vel 0.0030 0. 0. accel -0.005370861556 0. 0.
node 51 vel 0.0030 0. 0. accel -0.005370861556 0. 0.
node 52 vel 0.0030 0. 0. accel -0.005370861556 0. 0.
node 55 vel 0.0030 0. 0. accel -0.005370861556 0. 0.
node 54 vel 0.0030 0. 0. accel -0.005370861556 0. 0.
node 55 vel 0.0030 0. 0. accel -0.005370861556 0. 0.
node 56 vel 0.0030 0. 0. accel -0.005370861556 0. 0.
//----------------------------------------------------------------------------\\
// Options and switches \\
//----------------------------------------------------------------------------\\
SET TRANSIENT HUGHES BETA 0.2505 GAMMA 0.5 ALPHA -0.05 DAMPING MASS
COEFFICIENT 1.789 STIFFNESS COEFFICIENT 0.
//SET TRANSIENT HUGHES BETA 0.2505 GAMMA 0.5 ALPHA -0.05 DAMPING MASS
COEFFICIENT 0. STIFFNESS COEFFICIENT 0.1396
SET HUGHES_ALPHA_METHOD
//----------------------------------------------------------------------------\\
232
// Executing \\
//----------------------------------------------------------------------------\\
STEP id 1 TYPE DYNAMIC name "Load No. 1" AT 0.0 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * 0.001094800003
STEP id 2 TYPE DYNAMIC name "Load No. 2" AT 0.1 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * 0.001077716015
STEP id 3 TYPE DYNAMIC name "Load No. 3" AT 0.2 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * 0.001043814628
STEP id 4 TYPE DYNAMIC name "Load No. 4" AT 0.3 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * 0.000993624865
STEP id 5 TYPE DYNAMIC name "Load No. 5" AT 0.4 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * 0.000927929917
STEP id 6 TYPE DYNAMIC name "Load No. 6" AT 0.5 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * 0.847754933E-3
STEP id 7 TYPE DYNAMIC name "Load No. 7" AT 0.6 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * 0.754351018E-3
STEP id 8 TYPE DYNAMIC name "Load No. 8" AT 0.7 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * 0.649175706E-3
STEP id 9 TYPE DYNAMIC name "Load No. 9" AT 0.8 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * 0.533870226E-3
STEP id 10 TYPE DYNAMIC name "Load No. 10" AT 0.9 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * 0.410233878E-3
STEP id 11 TYPE DYNAMIC name "Load No. 11" AT 1.0 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * 0.280195968E-3
STEP id 12 TYPE DYNAMIC name "Load No. 12" AT 1.1 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * 0.145785694E-3
STEP id 13 TYPE DYNAMIC name "Load No. 13" AT 1.2 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * 0.9100483E-5
STEP id 14 TYPE DYNAMIC name "Load No. 14" AT 1.3 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.127726738E-3
STEP id 15 TYPE DYNAMIC name "Load No. 15" AT 1.4 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.262560826E-3
STEP id 16 TYPE DYNAMIC name "Load No. 16" AT 1.5 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.393297741E-3
STEP id 17 TYPE DYNAMIC name "Load No. 17" AT 1.6 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.517897375E-3
STEP id 18 TYPE DYNAMIC name "Load No. 18" AT 1.7 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.634415394E-3
STEP id 19 TYPE DYNAMIC name "Load No. 19" AT 1.8 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.741033573E-3
ATENA Input File Format 233
STEP id 20 TYPE DYNAMIC name "Load No. 20" AT 1.9 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.836088172E-3
STEP id 21 TYPE DYNAMIC name "Load No. 21" AT 2.0 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.918095893E-3
STEP id 22 TYPE DYNAMIC name "Load No. 22" AT 2.1 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.985777035E-3
STEP id 23 TYPE DYNAMIC name "Load No. 23" AT 2.2 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.1038075457E-2
STEP id 24 TYPE DYNAMIC name "Load No. 24" AT 2.3 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.1074175059E-2
STEP id 25 TYPE DYNAMIC name "Load No. 25" AT 2.4 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.1093512517E-2
STEP id 26 TYPE DYNAMIC name "Load No. 26" AT 2.5 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.1095786078E-2
STEP id 27 TYPE DYNAMIC name "Load No. 27" AT 2.6 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.1080960265E-2
STEP id 28 TYPE DYNAMIC name "Load No. 28" AT 2.7 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.1049266428E-2
STEP id 29 TYPE DYNAMIC name "Load No. 29" AT 2.8 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.1001199139E-2
STEP id 30 TYPE DYNAMIC name "Load No. 30" AT 2.9 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.937508473E-3
STEP id 31 TYPE DYNAMIC name "Load No. 31" AT 3.0 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.859188300E-3
STEP id 32 TYPE DYNAMIC name "Load No. 32" AT 3.1 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.767460782E-3
STEP id 33 TYPE DYNAMIC name "Load No. 33" AT 3.2 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.663757294E-3
STEP id 34 TYPE DYNAMIC name "Load No. 34" AT 3.3 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.549696096E-3
STEP id 35 TYPE DYNAMIC name "Load No. 35" AT 3.4 LOAD CASE FIXED 1 * 1.0
INCREMENT 2 * -0.427057074E-3
step id 1 execute
step id 2 execute
step id 3 execute
step id 4 execute
step id 5 execute
step id 6 execute
step id 7 execute
234
step id 8 execute
step id 9 execute
step id 10 execute
step id 11 execute
step id 12 execute
step id 13 execute
step id 14 execute
step id 15 execute
step id 16 execute
step id 17 execute
step id 18 execute
step id 19 execute
step id 20 execute
step id 21 execute
step id 22 execute
step id 23 execute
step id 24 execute
step id 25 execute
step id 26 execute
step id 27 execute
step id 28 execute
step id 29 execute
step id 30 execute
step id 31 execute
step id 32 execute
step id 33 execute
step id 34 execute
step id 35 execute
/* end of file */
// Eigenvalue analysis
//
// A cantilever modelled by 4 nonlinear shells
// Cross sectional dimension width=height=1; length=40
//
// Exact solution: (see
c:\AtenaExamples\Examples\Dynamics\SpringWithLumpedMass\Eigenvalues\cantilever.mw
s)
//
// f1=0.0443Hz
// f2=0.278Hz
// f3=0.775Hz
//
// Calculated:
//
// f1=0.0445Hz
// f2=0.299Hz
// f3=0.945Hz
//----------------------------------------------------------------------------\\
// Material definition \\
//----------------------------------------------------------------------------\\
MATERIAL
id 1
name "Spring"
type "CC3DElastIsotropic"
E 30000000
Mu 0.00
Rho 156.
Alpha 1.200E-05
236
//----------------------------------------------------------------------------\\
// Element type definition \\
//----------------------------------------------------------------------------\\
ELEMENT TYPE
id 1
name "1D Truss"
type "CCAhmadElement33L9"
//----------------------------------------------------------------------------\\
// Geometry definition \\
//----------------------------------------------------------------------------\\
//----------------------------------------------------------------------------\\
// Joint coordinates definition \\
//----------------------------------------------------------------------------\\
JOINT COORDINATES
//----------------------------------------------------------------------------\\
// Element group definition \\
//----------------------------------------------------------------------------\\
ELEMENT GROUP
id 1
name "Spring"
type 1
material 1
geometry 1
ELEMENT INCIDENCES
1 1 13 15 3 6 18 20 8 9 14 10 2 11 19 12 7 4 16 17
5
ATENA Input File Format 239
2
13 25 27 15 18 30 32 20 21 26 22 14 23 31 24 19 16 28
29 17
3 25 37 39 27 30 42 44 32 33 38 34 26 35 43 36 31 28 40
41 29
4
37 49 51 39 42 54 56 44 45 50 46 38 47 55 48 43 40 52
53 41
//----------------------------------------------------------------------------\\
// Options and switches \\
240
//----------------------------------------------------------------------------\\
SET NUMBER_OF_EIGENVALS 5
SET MAX_EIGENVAL_ERROR 0.0001
SET MAX_NUMBER_OF_SSPACE_ITERATIONS 14
SET REQUEST_STURM_SEQUENCE_CHECK YES
SET MAX_NUMBER_OF_JACOBI_ITERATIONS 10
SET NUMBER_OF_PROJ_VECS 15
//----------------------------------------------------------------------------\\
// Executing \\
//----------------------------------------------------------------------------\\
// STEP ID 1 STATIC NAME "BCs and load" LOAD CASE 1 * 1.0 EXECUTE
/* end of file */
OUTPUT_FFT_Y | OUTPUT_FFT_ERR_Y |
OUTPUT_FFT_SPECT_A | OUTPUT_FFT_SPECT_T |
OUTPUT_FFT_SPECT_O |
OUTPUT_FILTERED_FFT_Y | OUTPUT_FILTERED _FFT_ERR_Y |
OUTPUT_FILTERED _FFT_SPECT_A | OUTPUT_FILTERED_FFT_SPECT_T
| OUTPUT_FILTERED _FFT_SPECT_O |
OUTPUT_SIEVED_X | OUTPUT_SIEVED_Y] +
&MULTILINEAR_FUNCTION_DATA:
“CCMultiLinearFunction” [ALLOC_POINTS dim]
{ XVALUES { { xi }+ | RETRIEVE "output_data_X" }
YVALUES { { yi }+ | RETRIEVE "output_data_Y" } }|
GENERATE &ANALYTIC_FUNCTION_DATA }
[ { REMOVE_ALL | REMOVE_ALL_FILTERS | REMOVE_ALL_GENERATED |
{ ADD_OMEGA_FILTER omega_min omega_max ] }n |
SIEVE_Y | SIEVE_FFT_Y | SIEVE_FILTERED_FFT_Y |
COPY_FFT_Y | COPY_FILTERED_FFT_Y | COPY_SIEVED |
FFT | X_DOWN xl | X_UP xh | N_MAX nmax | DX_APPROX dx |
PRINT_Y PRINT_FFT_Y | PRINT_FILTERED_FFT_Y |
PRINT_COEFFS |PRINT_FILTERED_COEFFS |PRINT_SIEVED_POINTS } ] +
[ALLOC_POINTS dim] allocates space for dim additional points. It does not effect the
actual number of input points. It merely speeds up input of the subsequent dim points by
preparing beforehand a cumulative storage they will need.
&ANALYTIC_FUNCTION_DATA:
"CCAnalyticFunction" Y_EQN "y_string" [X_MIN x_min] [X_MAX x_max]
[DX dx]
Note that in order to visualize these plots, (using Atena's Graph Series dialog) don't forget to
check the "Values'profile for fixed time" checkbox and set horizontal and vertical fixed time
to zero, see description of the PLOT output option.
Parameter Description
ELEMENT GROUP n Element group id in which the pre-defined crack is to be
inserted.
ELEMENT n Element id in which the pre defined crack is to be inserted.
INTEGRATION Integration point id in which the pre defined crack is to be
[JOINT] n inserted. This is an optional parameter, if it is not specified,
crack is inserted into all integration points.
NORMAL x1 x2 [x3] Crack normal direction.
Parameter Description
ELEMENT GROUP Delete element group from the model or a single element from
the specified element group.
E.g. ELEMENT GROUP 3 [ELEMENT 4]
ELEMENT TYPE Delete element type from the list of element type definitions.
E.g. ELEMENT TYPE 2
GEOMETRY Delete geometry from the model.
E.g. GEOMETRY 6
JOINT Delete joint from the model
E.g. JOINT 3
ATENA Input File Format 245
Parameter Description
MACRO_ELEMENT Delete macro element from the model.
E.g. MACRO_ELEMENT 3
MACRO_NODE Delete macro node from the model
E.g. MACRO_NODE 3
ENFORCED If not specified, all references to a deleted entity remain valid
even after the deletion, thereby it is possible later to re-input the
entity with new data. Otherwise, the entity and all references to
it get unconditionally removed.
The command specifies the name of the input file to be opened or closed. Following this
command the ATENA input stream will be redirected into this file or the default stream.
This command specifies the name of the message file to be opened or closed. All
messages following this command will be redirected to this file or the default stream. The
open file is overwrited or appended.
This command specifies the name of the error file to be opened or closed. All errors
following this command will be redirected to this file or the default stream. The open file
is overwrited or appended.
This command reads the finite element model state from the given binary file name. The
content of the finite element model is overwritten by the file contents.
This command writes the finite element model state to a binary file. It can write
immediately, e.g. STORE TO “file name“, or it can autimatically serialize each n-th,
e.g. STORE TO “file name“ EACH n STEPS, or it can carry out the serialization each
step and m-th substeps, e.g. STORE TO “file name“ EACH m SUBSTEPS, (for dynamic
and creep analyses only). In the case of automatic serialization by steps the filename is
appended by ".step_id". The serialization by substeps appends the file name by
"_substep_id.step_id".
If n==0, then it the automatic serialization is stopped.
PO_STOP_IF_ULS_AND_DLS_FLAG n | PO_STOP_ONLY_IF_UNSTABLE_FLAG n |
STOREY_NODES_IDS { n }+ | VERTICAL_AXIS_ID n | HORIZONTAL_AXIS_ID n |
STOREY_DLS_COEFF x | EXECUTE}n
Units: none
Default: 2.5
SOIL_FACTOR x Soil factor from Eurocode, (called S in
Eurocode)
Units: time
Default: 0
ACCEL_GROUND x Ground acceleration, (ULS), (called a_g in
Eurocode)
Units: length/time2
Default: 0
ACCEL_GROUND_D x Ground acceleration, (DLS), (called a_Dg in
Eurocode)
Units: length/time2
Default: 0
P_D x Relative displacement stopping value, (called
p_d in Eurocode)
Units: none
Default: 1.5
P_F x Relative force drop down coefficient to
violate PO ULS criterion, (called p_f in
Eurocode).
Units: none
Default: 0.8
EXT_P_F x Relative force drop down coefficient to
declare the analysis unstable and stop the
execution.
Units: none.
Default: 0.2
PO_STOP_ONLY_IF_UNSTABLE_FLA If n=1 the analysis continues until the stability
Gn criterion is failed (irrespective of the pushover
analysis status).
If n=0, the pushover analysis is completed
based on the pushoover analysis status and the
flag PO_STOP_IF_ULS_AND_DLS_FLAG.
Default: n=0
PO_STOP_IF_ULS_AND_DLS_FLAG n If n=1, the pushover analysis is completed
after both ULS and DLS criteria are met.
If n=0, to complete the analysis it suffices to
fulfill only the ULS critera.
ATENA Input File Format 251
Default: n=0
Units: none
Default: none
Example: 0 249 693
VERTICAL_AXIS_ID n Id of model axis to be considered vertical, i.e.
axis, where gravity load applies.
Units: none
Default: 3, (i.e. Z axis)
HORIZONTAL_AXIS_ID n Id of model axis, where the ground
acceleration is applied.
Units: none
Default: 1, (i.e. X axis)
STOREY_DLS_COEFF x Coefficient coeff DLS to calculate maximum
interstory drift: d r coeff DLS h . h is height of
store and d r is relative storey drift.
Units: none
Default: 0.005
EXECUTE Carry out pushover analysis immediately. (By
default, this command is not needed, as the
analysis is calculated automatically at the end
of execution of each load step).
252
Syntax:
&STATIC_INITIAL_CONDITIONS:
NODAL TEMPERATURE [SETTINGS]
{ &STATIC_MANUAL_INITIAL_VALUES_ENTRY |
&STATIC_GENERATED_INITIAL_VALUES }+
& STATIC_MANUAL_INITIAL_VALUES_ENTRY:
{ [BASE_TEMPERATURE base_temp ] [NODE n TEMPERATURE nodal_temp] }
Sub-Command Description
NODE n Set initial conditions for node n.
TEMPERATURE base_temp Specify initial nodal temperature for node n. This value is
added to the base temperature below.
Units: T
Default: 0.
BASE_TEMPERATURE Initial base temperature. This value is used for all nodes of the
nodal_temp structure.
Units: [T]
Default: 0.
& STATIC_GENERATED_INITIAL_VALUES:
NODAL [SETTING] SELECTION "selection_name" | CONST const | COEFF_X
coeff_x | COEFF_Y coeff_y | COEFF_Z coeff | GENERATE_TEMP }
Sub-Command Description
SELECTION Name of selection, for which the generation is requested.
"selection_name"
ATENA Input File Format 253
CONST const COEFF_X Generate reference temperature for nodes in the selection
coeff_x | COEFF_Y coeff_y | "selection_name" . The values are generated as linear
COEFF_Z coeff_z combination:
GENERATE_TEMP
temperature base _ temp const x coeff x y coeff y z coeff z
where x,y,z are coordinates of nodes of nodes in the selection.
Units: COEFF_F, COEFF_M, COEFF_Z: [T/L]
CONST: [T]
Default: all constants are set to zero.
Note that initial reference temperatures can be set also by applying element temperature load
that import temperature history from a previous transport analysis of the structure. In this case
the reference nodal tepleratures corresponds to structural conditions at reference time of the
first applied element temperature load. As such values typically represent actual real
temperatures in the structure, the input described in this paragraph is not needed, (actually
temperatures from element temperature load would be added to temperatures from the
command &STATIC_INITIAL_CONDITIONS).
Example:
NODAL SETTING
NODE i TEMPERATURE temp
NODAL SETTING
BASE_TEMPERATURE base_temp // this value is added to specific node temperature,
NODEPROP ‘nodeprop’
ELEMPROP ‘elemprop’
In similar way, the command VERTEX can additionally include:
NODEPROP ‘nodeprop’
The parameter NODEPROP and/or ELEMPROP is used to generate the above mentioned
selection lists. Such a list is given name ‘nodeprop’ resp. ‘elemprop’ (notice use of single
quote ‘ instead of usual double quote “ !) and it will contain identification ids of all internal
FE nodes, resp. elements that were used to generate the T3D entity with the additional
parameters. Specify the parameters NODEPROP and ELEMPROP also for boundary entities,
(such as for surfaces of T3D region), if the generated list should include also boundary nodes
and elements of the T3D entity.
Use of the above new T3D commands and subcommands is demonstrated in the enclosed
sample AtenaWin analyses.
Example:
surface 11 curve 102 100 103 12 equidistant
The subcommand EQUIDISTANT can also be used for unstructured meshes. In this case,
however, no curve with the EQUIDISTANT property is allowed to have fixed vertices and
splitting of a copied curve (as described above) will help.
256
Note also, that the EQUIDISTANT is not always 100% accurate, especially in case of a
higher order meshes.
Syntax:
Only vertices with nearly the same coordinates get connected. The “same” property is judged
based on 1% octree mesh size. Octree is a special technique by which the 3D space around
the model is subdivided into brick shaped regions in order to facilitate faster searching
methods. It works for both structured and unstructured meshes. An error message is produced
and the generation is terminated, if for a SLAVE node no master node is found.
The regular selection lists includes only nodes or elements within the entity and outside its
boundary. They are created automatically during the mesh generation by T3D and they are
using an actual setting of &T3D_EXPAND_SETTINGS during the generation. The expanded
selection lists are regular selection lists expanded by adding nodes and elements on
boundaries of the appropriate entity. They are created by commands
&T3D_EXPAND_SETTINGS after the T3D mesh generation, i.e. in time, when the regular
lists are available.
Syntax:
&T3D_EXPAND:
T3D_EXPAND_SELECTIONS { [ &T3D_EXPAND_SETTINGS ] }+ {
[&T3D_EXPAND_ENTITY] }+
&T3D_EXPAND_SETTINGS :
&T3D_EXPAND_ENTITY:
PROP_GENERATION=NONE and
PROP_GENERATION=SEMIAUTOMATIC. If
PROP_GENERATION equals to AUTOMATIC,
then the nodeprop is ignored, (or reserved) and
DEF_VERTEX_FMT_FOR_NODES "vertex_fmt"
definition is used instead.
Default: "$N$V%i"
Example: “$N$Vertex%i”.
In this case, e.g. all finite nodes associated with a
vertex 13 will be listed in a selection list that calls
$N$Vertex13.
DEF_MNODE_FMT_FOR_NODES The same definition as the above for
"mnode_fmt" DEF_VERTEX_FMT_FOR_NODES, however, it
applies for macro nodes.
Default: "$N$MN%i"
Example: “$N$MacroNode%i”.
In this case, e.g. all finite nodes associated with a
macro node 13 will be listed in a selection list that
calls $N$MacroNode13.
DEF_CURVE_FMT_FOR_NODES The same definition as the above for
"curve_fmt" DEF_VERTEX_FMT_FOR_NODES, however, it
applies for curves.
Default: "$N$C%i"
Example: “$N$Curve%i”.
In this case, e.g. all finite nodes associated with a
curve 13 will be listed in a selection list that calls
$N$Curve13.
DEF_PATCH_FMT_FOR_NODES The same definition as the above for
"patch_fmt" DEF_VERTEX_FMT_FOR_NODES, however, it
applies for patches.
Default: "$N$P%i"
Example: “$N$Patch%i”.
In this case, e.g. all finite nodes associated with a
patch 13 will be listed in a selection list that calls
$N$Patch13.
DEF_SURFACE_FMT_FOR_NODE The same definition as the above for
S "surface_fmt" DEF_VERTEX_FMT_FOR_NODES, however, it
applies for surfaces.
Default: "$N$S%i"
Example: “$N$Surface%i”.
In this case, e.g. all finite nodes associated with a
260
DEF_MELEMENT_FMT_FOR_ELE
MENTS "melement_fmt"
DEF_BAR_REINFORCEMENT_FM
T_FOR_ELEMENTS "rc_fmt"
&COORDINATES_SPEC:
COORDINATES { [ID] n [NCOORDS] ncoords [X] { x }ncoords }+
This command adds new macro joints to the model. The joints are used for example for
reinforcement bar generation. Each macro joint coordinate should be on a separate line, e.g.
[ID] n [X] x1 x2 x3
If ncoords is not specified, it is by default equal to problem dimension, see &TASK.
This command adds new macro joints to the model or deletes the existing one. The joints
are used for example for reinforcement bar generation. Each macro joint coordinate should
be on a separate line, e.g.
[ID] n [X] x1 x2 x3
If ncoords is not specified, it is by default equal to problem dimension, see &TASK. The
“ENFORCED” keyword has the same meaning as in “DELETE” command.
262
Syntax:
&MACRO_ELEMENT
MACRO_ELEMENT melem_id { &GENERATE_SPEC | &UPDATE_SPEC |
&DELETE_SPEC )
&GENERATE_SPEC:
GENERATE TYPE "type_str" { [THROUGH] NODES { mnode_id }+ |
GROUP group_id | COUNTER [{BASE | ELEMENT_BASE | NODAL_BASE}]
base_id | NAME "melem_name" |
ELEMPROP "elem_prop" | NODEPROR "node_prop" { ID id}+ |
MACRO_ELEM_DATA_SPEC | EXECUTE }+
&DELETE_SPEC:
{ ENFORCED DELETE } | { DELETE }
Example:
MACRO_ELEMENT 1000 GENERATE TYPE
"CCIsoMacroElement<xxxxxxxx_x_x>" THROUGH NODES 201 202 204 203
101 102 104 103 205 206
GROUP 1 COUNTER ELEMENT_BASE 1 NODAL_BASE 1 NAME "Macro
block 1"
ELEMPROP "Block_1"
NODEPROP "N1" ID 1
NODEPROP "N2" ID 2
NODEPROP "N3" ID 3
NODEPROP "N4" ID 4
NODEPROP "N5" ID 5
NODEPROP "N6" ID 6
NODEPROP "N7" ID 7
NODEPROP "N8" ID 8
QUADRATIC
SHAPE HEXA DIR 1 DIVISION 4 DIR 2 DIVISION 3 DIR 3
DIVISION 2 DR 0.2 0.2
EXECUTE
Syntax:
SOURCE_NODES {id}3 | {id}4 | SOURCE_ELEMPROP "elemprop" |
SOURCE_GROUP id | SOURCE_NODEPROP "nodeprop" | ACCOMPLISH
count | [TIMES] }+
Parameter Description
SOURCE_NODES {id}3 Defines ids of source macronodes, whose coordinates should be
| {id}4 transformed into destination coordinates of nodes {
[THROUGH] NODES { mnode_id }+ . Note that this input data
only defines transformation of the model and no actual
macronodes will be copied. 2D resp. 3D problem needs 3 resp. 4
of such nodal source-destination nodal pairs.
SOURCE_ELEMPROP All elements defined in the selection "elemprop" will be copied.
"elemprop" If this parameter is not specified, then all (current) elements
from SOURCE_GROUP are copied.
SOURCE_NODEPROP Selection list "nodeprop" of source nodes, whose copy should
"nodeprop" be included in a new node selection. Name of the selection will
be concatenation of destination “elemprop” and "nodeprop" . If
more copies are generated, (see ACCOMPLISH count TIMES
data), the name is appended by “$n”, where n is number of
additional copy. The same applies for destination “elemprop” .
SOURCE_GROUP id Id of element group that contains the elements
SOURCE_ELEMPROP "elemprop". By default, GROUP
group_id is used.
ACCOMPLISH count | Specifies number of copies to be generated. By default one copy
[TIMES] is created, i.e. count=1.
KEEP_ELEM_IDS Copy source element id to target element id. Otherwise a new
target id is generated, which is behaviour by default. This
function is only available for the first copy of the elements, (i.e.
count=1) and at the same time target and source group ids must
be different.
Example:
MACRO_ELEMENT 1001 GENERATE TYPE "CCCopyElementSelection" THROUGH
NODES 102 107 104 202
GROUP 1 NAME "Macro block 2"
ELEMPROP "Block_2"
SOURCE_NODES 101 102 103 201 SOURCE_ELEMPROP "Block_1"
SOURCE_NODEPROP "N1N4N5N8" "N5N6N7N8" "N5N8" "N5N6"
EXECUTE
ATENA Input File Format 267
Parameter Description
SOURCE_NODES {id}3 Defines ids of source macronodes, whose coordinates should be
| {id}4 transformed into destination coordinates of nodes {
[THROUGH] NODES { mnode_id }+ . Note that this input data
only defines transformation of the model and no actual
macronodes will be copied. 2D resp. 3D problem needs 3 resp. 4
of such nodal source-destination nodal pairs.
INPUT_NODEPROP All nodes defined in the selection "prop_in" will be copied.
"prop_in"
Parameter Description
SOURCE_ELEMPROP All elements defined in the selection "elemprop" will be copied.
"elemprop" If this parameter is not specified, then all (current) elements
from SOURCE_GROUP are copied.
SOURCE_GROUP id Id of element group that contains the elements
SOURCE_ELEMPROP "elemprop". By default, GROUP
group_id is used.
ACCOMPLISH count | Specifies number of copies to be generated. By default one copy
[TIMES] is created, i.e. count=1.
KEEP_ELEM_IDS Copy source element id to target element id. Otherwise a new
target id is generated, which is behaviour by default. This
function is only available for the first copy of the elements, (i.e.
count=1) and at the same time target and source group ids must
ATENA Input File Format 269
be different.
Example:
MACRO_ELEMENT 1001 GENERATE TYPE "CCCopyElementSelection" THROUGH
NODES 102 107 104 202
GROUP 1 NAME "Macro block 2"
ELEMPROP "Block_2"
SOURCE_NODES 101 102 103 201 SOURCE_ELEMPROP "Block_1"
SOURCE_NODEPROP "N1N4N5N8" "N5N6N7N8" "N5N8" "N5N6"
EXECUTE
"Block_3_Block_2_N2N3N6N7<-"
"Block_3_Block_2_N2N3N6N7->"
"Block_3$1_Block_2_N2N3N6N7<-"
"Block_3$1_Block_2_N2N3N6N7->"
"Block_3$2_Block_2_N2N3N6N7<-"
"Block_3$3_Block_2_N2N3N6N7->"
Syntax:
SOURCE_GROUP id | SOURCE_NODE id | SOURCE_ELEMPROP "elemprop" |
SOURCE_NODEPROP "nodeprop" | ACCOMPLISH count | [TIMES] }+
Parameter Description
SOURCE_NODE id Defines id of a bottom macronode for the extrusion vector. The
top node is defined by NODE id.
SOURCE_ELEMPROP All elements defined in the selection "elemprop" with nodes
"elemprop" defined in SOURCE_NODEPROP "nodeprop" will be used as
270
Parameter Description
[EMBEDDED] [IN] Interval of element groups defining the “master” material, i.e.
[SOLID] [SOLIDS] { AT | solids ids, where the bar should be generated. In other words,
FROM} solid_group_id_1 the bar will be embedded in the specified material groups.
[TO solid_group_id_2]
{ NORMAL | TINY If TINY size is defined, then the algorithm used to generate
[SIZE] } elements of the bar works correctly even in the case, that more
neighboring NODES are located with the same elements. If it
is not the case, use of NORMAL size is preferable, as it results
ATENA Input File Format 271
11
Not available in ATENA version 4.3.1 and older.
272
Parameter Description
{PLANE n THROUGH Specify 3 or 4 macronodes ids defining triangular or
NODES { n1 n2 n3 n4 | quadrilateral reinforcement plane.
n1 n2 n3 }}+
[EMBEDDED] [IN] Interval of element groups defining the “master” material, i.e.
[SOLID] [SOLIDS] { AT solids ids, where the bar should be generated. In other words,
| FROM} the bar will be embedded in the specified material groups.
solid_group_id_1 [TO
solid_group_id_2]
NORMAL | TINY If TINE size is defined, then the algorithm used to generate
[SIZE] elements of the smeared reinforcement planes works correctly
even in the case, that more neighboring NODES are located
with the same elements. If it is not the case, use of NORMAL
size is preferable, as it results in much faster element generation.
Default value: NORMAL
[SIZE] MINIMUM x Minimum distance between nodes of generated element. If not
satisfied, newly generated node is ignored.
Default value: 0 [length units]
Example:
MACRO_ELEMENT 1001 GENERATE TYPE "CCDiscretePlaneReinforcementME"
PLANE 1 THROUGH NODES 1001 1005 1006 1004
PLANE 2 THROUGH NODES 1005 1002 1003
PLANE 3 THROUGH NODES 1005 1003 1006
NAME "Bottom reinforcement"
MINIMUM 0.
GROUP 10 EMBEDDED AT 1
ELEMPROP "Plame_1"
NODEPROP "N1" ID 1
ATENA Input File Format 273
NODEPROP "N2" ID 2
NODEPROP "N3" ID 3
NODEPROP "N4" ID 4
EXECUTE
/M CCStructuresTransport:
The CCStructuresTransport module is an extension of CCFEModel, (being the base for all
engineering modules in ATENA) and hence most input command for the transport analysis
are the same as those e.g. for static analysis of structures. This section describes additional
commands that are relevant only for the transport analysis.
There are a few input commands that are meaningful only for transport analysis. These are
commands:
related to temporal time integration, &Transport Set parameters (and problem’s time
step marching execution as it is),
needed for definition of transport finite element, &Transport finite elements,
specifying transport constitutive material model, &Transport constitutive material,
inputting structural initial state conditions, &Transport initial value of state variables,
&History export related commands
&Transport analysis additional output data.
Note also that only Modified Newton-Raphson or Full Newton-Raphson execution method
can be used.
0
[C_TEMP_W CTw ] | [C_H_T CTt0 ] | [K_TEMP_H_FNC_TEMP_ID f KTTh ] |
[K_TEMP_TEMP_FNC_TEMP_ID f KTTT ] | [K_TEMP_W_FNC_TEMP_ID
f KTTw ] | [K_TEMP_GRAV_FNC_TEMP_ID f KTTgrav ] |
[C_TEMP_H_FNC_TEMP_ID fCTTh ] | [C_TEMP_TEMP_FNC_TEMP_ID fCTTT ] |
[C_TEMP_W_FNC_TEMP_ID fCTTw ] | [C_TEMP_T_FNC_TEMP_ID f CTTt ] |
[K_TEMP_H_FNC_H_ID f KhTh ] | [K_TEMP_TEMP_FNC_H_ID f KhTT ] |
[K_TEMP_W_FNC_H_ID f KhTw ] | [K_TEMP_GRAV_FNC_H_ID f KhTgrav ] |
[C_TEMP_H_FNC_H_ID f ChTh ] | [C_TEMP_TEMP_FNC_H_ID fChTT ] |
[C_TEMP_W_FNC_H_ID fChTw ] | [C_TEMP_T_FNC_H_ID fChTt ] |
[K_TEMP_H_FNC_T_ID f Kt Th ] | [K_TEMP_TEMP_FNC_T_ID f Kt TT ] |
[K_TEMP_W_FNC_T_ID f Kt Tw ] | [K_TEMP_GRAV_FNC_T_ID f KhTgrav ] |
[C_TEMP_H_FNC_T_ID fCtTh ] | [C_TEMP_TEMP_FNC_T_ID fCtTT ] |
[C_TEMP_W_FNC_T_ID fCtTw ] | [C_TEMP_T_FNC_T_ID fCtTt ] }+ ]
[WATER
0 0 0 0
{ [D_H_H Dwh ] | [D_H_TEMP DwT ] | [D_H_W Dww ] | [D_H_GRAV Dwgrav ]|
0 0 0
[C_H_H Cwh ] | [C_H_TEMP CwT ] | [C_H_W Cww ] [C_H_T Cwt0 ] |
[D_H_H_FNC_H_ID f Dhwh ] | [D_H_TEMP_FNC_H_ID f DhwT ] |
[D_H_W_FNC_H_ID f Dhww ] | [D_H_GRAV_FNC_H_ID f Dhwgrav ] |
[C_H_H_FNC_H_ID fChwh ] | [C_H_TEMP_FNC_H_ID fChwT ] |
[C_H_W_FNC_H_ID f Chww ] | [C_H_T_FNC_H_ID fChwt ] |
[D_H_H_FNC_TEMP_ID f DTwh ] | [D_H_TEMP_FNC_TEMP_ID f DTwT ] |
[D_H_W_FNC_TEMP_ID f DTww ] | [D_H_GRAV_FNC_TEMP_ID f DTwgrav ] |
[C_H_H_FNC_TEMP_ID fCTwh ] | [C_H_TEMP_FNC_TEMP_ID fCTwT ] |
[C_H_W_FNC_TEMP_ID f CTww ] | [C_H_T_FNC_TEMP_ID fCTwt ] |
[D_H_H_FNC_T_ID f Dt wh ] | [D_H_TEMP_FNC_T_ID f Dt wT ] |
[D_H_W_FNC_T_ID f Dt ww ] | [D_H_GRAV_FNC_T_ID f Dt wgrav ] |
[C_H_H_FNC_T_ID fCtwh ] | [C_H_TEMP_FNC_T_ID f CtwT ] |
[C_H_W_FNC_T_ID fCtww ] | [C_H_T_FNC_T_ID fCtwt ] }+ ]
Table 169: &Parameters of the & CCModelBaXi94 within the transport analysis
Parameter Description
CONCRETE TYPE n_type Type of concrete, resp. type of cement. n_type = <1..4>,
n_type =1 for Portland cement etc.
Default value: 1
RATIO_WC ratio Water cement ratio. The allowed range is <0.3..0.7.
ATENA Input File Format 277
Table 170: &Parameters of the & CCTransportMaterial within the transport analysis
Input parameters for user-defined constitutive law for flow governing equations
278
Heat :
W
div( qw )
t
h T w
CTh CTT CTw CTt div KTh grad (h) KTT grad (T ) KTw grad ( w) KTgrav
t t t
Moisture :
Q
div( qT )
t
h T w
Cwh CwT Cww Cwt div Dwh grad (h) DwT grad (T ) Dww grad ( w) Dwgrav
t t t
W , Q states for total amount of moisture per unit volume, [kg/m3 ] and total amount of energy
per unit volume, [J/m3 ]. Note that positive value of CTt , Cht causes consumption, so that e.g.
hydration heat must be input as negative number. Input always a label followed by an
associated real value, (for constant parameter) or integer id of a previously defined function,
(for a function definition). If a parameter is skipped, it is assumed either zero or the
associated function is assumed to have value 1, i.e. neglected. The T subscript for
temperature related parameters is replaced by TEMP string. The subscripts for humidity,
water content and time, i.e. sink related terms remain unchanged, i.e. H, W, t respectively.
For example CTT is entered as C_TEMP_TEMP etc. All functions are defined separately.
Each such a definition is refered by its id, i.e. a integer number. This integer is then specified
as a value following the appropriate label. For example the function fCtTT (t ) is defined with id
k. Then, the material data input would read C_TEMP_TEMP_FNC_ID k.
Significance of the parameters is as follows:
ATENA Input File Format 279
Default values: All functions are constant and equal to one, i.e. they are disregarded. All
other parameters are by default zero with the following exceptions:
kg kg
Chh 225 3
, Dhh 1.5 E 6
m sm
J J
CTT 2.55 E 6 3 , KTT 2.1
mC s mC
Parameter Description
DOH_FNC_ID id Id of degree of hydration DoH(time) function.
It prevails input of DOH25_FNC_ID and
analytical calculation of DoH(time) using B1,
B2, ALPHAINF and ETA.
DOH25_FNC_ID id Id of degree of hydration DoH25(time)
function, i.e. DoH function for reference
temperature 25oC and relative humidity 1. It is
280
Units: [mass]
Default value: 2086 kg
FILLER_MASS val Filler mass in concrete m filler .
Units: [mass]
Default value: 69 kg
CEMENT_DENSITY val Cement density.
Units: [mass/length3]
Default value: 3220 kg/m3
WATER_DENSITY val Water density.
Units: [mass/length3]
282
Units: [mass/length3]
Default value: 127 kg/m3
H80 val Relative humidity h80 for w80 .
Units: [-]
Default value: 0.8
W80 val Water saturation w80 for h80 .
Units: [mass/length3]
Default value: 40 kg/m3
TEMP0 val Parameter T0 to calculate saturaturated water
vapour pressure psat for temperatures
T 0 0C .
Units: [temperature]
Default value: 234.18 0C.
A_WV val Parameter a to calculate saturated water
vapour pressure psat for temperatures
T 0 0C .
Units: [-]
Default value: 17.08
Units: [temperature]
Default value: 272.44 0C.
A_WV_ICE val Parameter a to calculate saturated water
vapour pressure psat for temperatures
T 0 0C .
Units: [-]
Default value: 22.44
EA val Acxtivation energy Ea
Units: [energy/mol]
Default value: 38300 J/mol
All remaining input data in the sections They are the same as those for
TEMPERATURE and WATER: &CCTransportMaterial_PARAMS, except by
Chh 0.0, Dhh 0.
default
CTT 0., KTT 0.
IsoTriangle6_2D CCIsoTriangle6_2D
IsoTriangle3_AS Axisymmetric triangular isoparametric elements CCIsoTriangle3_ASy
ym m
... ...
IsoTriangle6_AS CCIsoTriangle6_ASy
ym m
IsoBrick8_3D Hexahedral isoparametric elements CCIsoBrick8_3D
... ...
IsoBrick20_3D CCIsoBrick8_3D
IsoWedge6_3D Wedge isoparametric elements CCIsoWedge6_3D
.... ....
IsoWedge15_3D CCIsoWedge15_3D
IsoTetra4_3D Tetrahedral isoparametric elements CCIsoTetra4_3D
.... ....
IsoTetra10_3D CCIsoTetra10_3D
IsoTruss2_2D, Truss isoparametric elements, 2D, 3D and CCIsoTruss2_2D,
IsoTruss3_2D, axisymmetric. CCIsoTruss3_2D,
IsoTruss2_3D, CCIsoTruss2_3D,
IsoTruss3_3D, CCIsoTruss3_3D,
IsoTruss2_ASym, CCIsoTruss2_ASym,
IsoTruss3_ASym, CCIsoTruss3_ASym,
3_2D ...
... CCIsoTriangle6_2D
IsoTriangle Gauss
6_2D
IsoTriangle Gauss Axisymmetric triangular isoparametric elements CCIsoTriangle3_ASy
3_ASym m
... ...
IsoTriangle Gauss CCIsoTriangle6_ASy
6_ASym m
IsoBrick Gauss Hexahedral isoparametric elements CCIsoBrick8_3D
8_3D
...
...
CCIsoBrick8_3D
IsoBrick Gauss
20_3D
IsoWedge Gauss Wedge isoparametric elements CCIsoWedge6_3D
6_3D
....
....
CCIsoWedge15_3D
IsoWedge Gauss
15_3D
IsoTetra Gauss Tetrahedral isoparametric elements CCIsoTetra4_3D
4_3D
....
....
CCIsoTetra10_3D
IsoTetra Gauss
10_3D
Syntax:
&INITIAL_CONDITIONS:
NODAL { MAT_H_TEMP | MAT_TEMP_H| H_TEMP_MAT | TEMP_H_MAT |
TEMPERATURE | HUMIDITY | MATERIAL} [SETTINGS] {
&MANUAL_INITIAL_VALUES_ENTRY |
&GENERATED_INITIAL_VALUES }
&MANUAL_INITIAL_VALUES_ENTRY:
{ NODE n | TYPE type | H h | W w | TEMP temp }+
ATENA Input File Format 287
Sub-Command Description
NODE n Set initial conditions for node n.
TYPE type Specify type of material used in node n. Note that transport
analysis is integrated in finite nodes rather than integration
nodes in finite elements and hence material model is related to
finite nodes (and not finite elements).
Hh|Ww Set initial condition for relative humidity h. Moisture
conditions can be equivalently also set by setting the amount
of water content w, see the ATENA Theoretical manual for
definition of w(.
TEMP temperature Set initial temperature in the node [Kelvin]
&GENERATED_INITIAL_VALUES:
Sub-Command Description
SELECTION Name of selection, for which the generation is requested.
"selection_name"
TYPE type Specify type of material used in nodes in the selection.
{GENERATE_H | Keyword for entities to be generated. The value is generated
GENERATE_W | as linear combination:
GENERATE_TEMP }1
CONST const | COEFF_X
value const x coeff x y coeff y z coeff z
coeff_x | COEFF_Y coeff_y |
COEFF_Z coeff_z
x,y,z are coordinates of nodes, where the generation is
processed
Example:
NODAL MAT_H_TEMP SETTING NODE 1 MATERIAL TYPE 1 H 1. TEMP 20
NODAL SELECTION “my_selection” GENERATE_TYPE 1
CONST 0.5 COEFF_X 0. COEFF_Y -0.6523648649 COEFF_Z 0. GENERATE_H
CONST -10. COEFF_X 0. COEFF_Y 0. COEFF_Z 0 GENERATE_T
288
Parameter Description
&TRANSIENT Set transient analysis and set some parameters for it.
&CONVERGENCE_CRIT Convergence criteria for the transport analysis
ERIA
&TRANSIENT:
TRANSIENT { [TIME] CURRENT x | [TIME] INCREMENT x |
TIME_INTEGRATION { {CRANK_NICHOLSON | THETA x }+ |
ADAMS_BASHFORTH } | REFERENCE_ETA eta}+
Parameter Description
[TIME] CURRENT x Sets current time.
[TIME] INCREMENT x Sets time increment in steps.
TIME_INTEGRATION Set type of temporal integration scheme. If this parameter is not
input, then CRANK_NICHOLSON integration will be used.
CRANK_NICHOLSON Use linear trapezoidal integration.
THETA x parameter for trapezoidal integration. By default = 0.5.
Several other linear temporal integration may be utilized
depending on the, e.g. implicit Newton integration for = 1,
explicit integration for = 0 etc. For good compromise between
convergence and possibility of oscillations values about =
0.85 is recommended.
ADAMS_BASHFORTH Adams – Bashforth quadratic temporal integration.
REFERENCE_ETA eta Damping factor. t dt t t dt . 0.3..1 ; 1 set
totally un-damped analysis.
Default: 1
&CONVERGENCE_CRITERIA:
{ ABSOLUTE [ ERROR] | RELATIVE [ERROR] } | TEMPERATURE ERROR x |
HUMIDITY ERROR x | STEP_STOP_TEMPERATURE ERROR FACTOR x |
STEP_STOP_HUMIDITY ERROR FACTOR x |
ITER_STOP_TEMPERATURE ERROR FACTOR x | ITER_STOP_HUMIDITY
ATENA Input File Format 289
Parameter Description
ABSOLUTE [ERROR] The convergence criteria values are computed using the absolute
norm that is using the maximal element of an array in its absolute
value. The error is then computed by dividing an iterative value
with the value cumulated within the whole step.
RELATIVE [ERROR] The convergence criteria values are computed using the
Euclidean norm. The error is then computed by dividing an
iterative value with the value cumulated within the whole step.
TEMPERATURE ERROR x Convergence limit for absolute value of temperature increments.
Default value is 0.01.
E.g. TEMPERATURE ERROR x
HUMIDITY ERROR x Convergence limit for absolute value of humidity increments.
Default value is 0.01.
E.g. HUMIDITY ERROR x
STEP_STOP_TEMPERATU Factors for appropriate convergence criterion value. If a
RE ERROR FACTOR x | convergence criterion value multiplied by the appropriate factor
STEP_STOP_HUMIDITY exceeds the related calculated analysis error, then the execution is
ERROR FACTOR x immediately killed. They are two sets of factors: the first one for
| checking each iteration and the other one to be exercised at the
ITER_STOP_TEMPERATU end of each step. The default value for iteration related factors is
RE ERROR FACTOR x | 1000, whilst the default value for step related factors is 10.
ITER_STOP_HUMIDITY E.g.
ERROR FACTOR x
SET Absolute
Step_stop_humidity error factor 15.
Step_stop_ temperature error factor 53
Iter_stop_ humidity error factor 201
Iter_stop_ temperature error factor 203
SET Relative
Step_stop_ humidity error factor 54
Step_stop_ temperature error factor 56
Iter_stop_ humidity error factor 204
Iter_stop_ temperature error factor 206
NEGLIGIBLE_TEMPERAT Values that are negligible, i.e. that can be ignored. By default
URE x | NEGLIGIBLE they are set to 1.E-11.
_HUMIDITY x E.g.
SET
Absolute error Negligible_temperature 0.1
Relative error Negligible_temperature 0.2
290
Syntax:
&HISTORY_EXPORT:
HISTORY [{APPEND | OVERWRITE}] [EXPORT] [TO] [GEOMETRY
geometry_filename] | [RESULTS] results_filename] 2
Parameter Description
results_filename Name of binary file with the history. It must be the same as that
specified for HISTORY IMPORT command in the
CCStructuresCreep module. It should be enclosed in double
quote character (“).
geometry_filename Name of binary file with geometry of the exported model. It
must be the same as that specified for HISTORY IMPORT
command in the CCStructuresCreep module. It should be
enclosed in double quote character (“). If omitted, identical
imported and current models are assumed.
[{APPEND | Open option for the file. By default, the file gets during
OVERWRITE }] execution overwritten.
[EXPORT] [TO] Ignored keywords.
&BOUNDARY_ELEMENT_LOAD
&BODY_ELEMENT_LOAD
&FIRE_BOUNDARY LOAD
&MOIST_TEMP_BOUNDARY_LOAD
&FIRE_BOUNDARY _LOAD:
FIRE_BOUNDARY [GROUP group_id [ TO group_id_to [ BY group_id_by]]
[ELEMENT { element_id [ TO element_id_to [ BY element_id_by]] |
ATENA Input File Format 291
Important: Note that unlike other types of static loads (that are input in incremenental
manner), the fire boundary load has character of a load potential and thus it must be input in
total form. Therefore the load describes (total) fire load conditions !
Parameter Description
[FIRE] [TYPE] { GENERIC | Type of fire load to be applied.
NOMINAL_HC, MODIFIED_HC
USER_CURVE
[CONVECTION hc Convection heat transfer coefficient [W/m2/K].
Default value 50 W/m2K.
EMISSIVITY r Emissivity parameter.
Default value 0.56.
TEMPERATURE_MAX Tg ,ref Max. temperature parameter.
[TIME_FUNCTION time_id Id of an user-defined time dependent function. It acts
as an extra multiplier of the generated or directly
inputed fire boundary load.
TEMPERATURE_MIN Tg ,min Ambient temperature prior the fire broke up. (Any
generated temperature cannot fall below this value).
NODES "boundary_nodes_list" List of boundary load that are load.
({EDGE | Type of boundary load, that is applicable for the
EDGE_NO_DUPLICATES } | given fire load. For more explanation see
SURFACE } &BOUNDARY_ELEMENT_LOAD.
&MOIST_TEMP_BOUNDARY_LOAD:
MOIST_TEMP_BOUNDARY &ELEM_LOAD_DATA
&MOISTURE_FLUX_DUE_TO_RELATIVE_HUMIDITY_GRADIENT
&MOISTURE_ FLUX_DUE_TO_HUMIDITY_RATIO_GRADIENT
&MOISTURE_ FLUX_DUE_TO_CEMSTONE_CALC
&HEAT_ FLUX_DUE_TO _TEMPERATURE_GRADIENT
&HEAT_ FLUX_DUE_TO_EVAPORATED_MOISTURE
&COMMON_ MOIST_TEMP _BC_DATA [MERGE [ MERGE_STRING str ]]
[NO_ELEM_OUTPUT ] { CONSTANT_IN_STEP | VARIABLE_IN_STEP |
SEMIVARIABLE_IN_STEP | [MULTIPLE {YES|NO}] }
&MOISTURE_FLUX_DUE_TO_RELATIVE_HUMIDITY_GRADIENT:
[{ACCOUNT|NEGLECT} [GRADIENT] [OF] RELATIVE_HUMIDITY]
[CONVECTION_W hcw ]
&MOISTURE_ FLUX_DUE_TO_HUMIDITY_RATIO_GRADIENT:
[{ACCOUNT|NEGLECT} [GRADIENT] [OF] HUMIDITY_RATIO]
[EVAPORATION_MOISTURE ] [AIR_PRESSURE p] [AIR_VELOCITY
v ] [AIR_VELOCITY_FUNCTION air_velocity_fnc_id]
ATENA Input File Format 293
&MOISTURE_ FLUX_DUE_TO_CEMSTONE_CALC:
[{ACCOUNT|NEGLECT} [GRADIENT] [OF]
HUMIDITY_CEMSTONE_CALC]
&HEAT_ FLUX_DUE_TO_EVAPORATED_MOISTURE:
[{ACCOUNT|NEGLECT} [GRADIENT] [OF] EVAPORATED_MOISTURE]
[EVAPORATION_HEAT hwe ]
Important: Note that unlike other types of static loads (that are input in incremenental
manner), the moisture-temperaturee boundary load has character of a load potential and thus
it must be input in total form. Therefore the load describes (total) moisture-temperature load
conditions !
Parameter Description
AMBIENT_HUMIDITY hg Ambient air relative humidity, [-].
Default value: 0.6
Ambient temperature, [oC].
AMBIENT_TEMPERATURE Tg
Default: 20 oC
CONVECTION_W hcw Convection moisture transfer coefficient [kg/s/m2].
Default value 0. kg/s/m2
EVAPORATION_MOISTURE Evaporation moisture transfer coefficient [kg/(m2s)].
Default value (25+19*v_)/(3600.) kg/s/m2, where v_
is air velocity in ms-2.
AIR_PRESSURE p Total (absolute) ambient air pressure, [Pa], (=sum of
partial dry air pressure and partial water vapour
pressure).
Default: normal pressure 101325Pa
AIR_VELOCITY v Average ambient air velocity, [m/s],
Default =0.m/s
Convection heat transfer coefficient [W/m2/K].
CONVECTION_T hcT
Default value 20 W/m2K.
EMISSIVITY_T rT Heat emissivity parameter, [-]
Default value 0.85.
294
Table 182: Transport analysis related Output-type keywords understood by the command
&OUTPUT for the location type NODES
Table 183: Transport analysis related Output-type keywords understood by the command
&OUTPUT for the location type NODES
Sample analysis:
-----------------
Analysis of a simple 2D wall comprising quadrilateral and
triangle
elements, subject to displacement load at nodes 600 and 700.
|y
|
|
| 300
|400 800 700
-x-> ------------------------------------| -> 3.33e-6 *
| | /| function 20
| | / |
| | 20 / |
| | / |
| 10 | / |
^| | / |
|| | / 15 |
x| | / | 3.33e-6 *
|| |/ | -> function 20
-x-> ------------------------------------------------------
-
100 200 600
(50) 500
x
*/
// Coordinate definition
JOINT COORDINATES
50 0. 0. // dummy object for deletion checking
298
100 0. 0.
200 0.1 0.
700 0.2 0.1
300 0.1 0.1
500 0.1 0.
400 0. 0.1
600 0.2 0.
800 0.1 0.1
// Material definition
MATERIAL ID 71 NAME "Steel" TYPE "CCPlaneStressElastIsotropic"
E 210000 mu 0.2 rho 0.0023 alpha 1.2e-5
MATERIAL ID 70 NAME "Steel" TYPE "CCPlaneStressElastIsotropic"
E 210000 mu 0.2 rho 0.0023 alpha 1.2e-5 // dummy object
for deletion checking
// Geometry definition
GEOMETRY ID 81 Name "Steel thickness" TYPE "2D" thickness 0.1
GEOMETRY ID 80 Name "Steel thickness" TYPE "2D" thickness 0.1
// dummy object for deletion checking
// Testing of deletion
DELETE ELEMENT GROUP 500
DELETE JOINT 50
DELETE ELEMENT GROUP 2000 ELEMENT 15
DELETE GEOMETRY 80
DELETE ELEMENT TYPE 90
DELETE MATERIAL 70
DELETE LOAD CASE ID 60
DELETE FUNCTION 10
Testing input data format - LHS and RHS boundary conditions; their values and sign.
(for 3D version see transp2_bricks_test.inp)
Structure:
Initial condition: dT/dy= -20/1.= -20; dT/dx=0; dh/... irrelevant, h=fixed everywhere
Flux: qy = K_TEMP_TEMP * dT/dy = 103680 * -20 = -2073600
External forces: sum(Q) = qy*width*thick= -2073600 * 0.15 * 10. = 3110400
Individual force: Q = sum(Q)/2 = 3110400/2 = 1555200
Internal forces: positive value corresponds to the flow in direction of outwards normal to the
boundary surface
External load: positive value corresponds to the flow in direction of inwards normal to the
boundary surface
top surface (nodes 9,10), i.e. y=1 .... internal forces negative, i.e. -1555200; external load
positive, i.e. 1555200
bottom surface (nodes 1,2), i.e. y=0 .... internal forces positive, i.e. 1555200; external load
negative, i.e. -1555200
ALL EXTERNAL LOADS as well as NON_ZERO LHS BCs (i.e. fixing psi, h) HAVE
INCREMENTAL CHARACTER.
This means that e.g. LOAD SIMPLE SELECTION "all9-10" dof 2 const 1555200. applied to
all steps
will produce external forces 1555200. in the 1st step, 3110400. in the 2nd step.... The same
applies to nonzero SUPPORT SIMPLE ..... specification.
NODAL SETTING
NODE 1 MATERIAL TYPE 1 H 1. TEMP 20
NODE 2 MATERIAL TYPE 1 H 1. TEMP 20
NODE 3 MATERIAL TYPE 1 H 1. TEMP 25
NODE 4 MATERIAL TYPE 1 H 1. TEMP 25
NODE 5 MATERIAL TYPE 1 H 1. TEMP 30
NODE 6 MATERIAL TYPE 1 H 1. TEMP 30
NODE 7 MATERIAL TYPE 1 H 1. TEMP 35
NODE 8 MATERIAL TYPE 1 H 1. TEMP 35
NODE 9 MATERIAL TYPE 1 H 1. TEMP 40
NODE 10 MATERIAL TYPE 1 H 1. TEMP 40
Boundary conditions:
*/
// Coordinate definition
JOINT COORDINATES // 4 elements 0.15*0.25 placed vertically
ATENA Input File Format 303
1 0. 0.
2 0.15 0.
3 0. 0.25
4 0.15 0.25
5 0. 0.5
6 0.15 0.5
7 0. 0.75
8 0.15 0.75
9 0. 1.
10 0.15 1.
// Material definition
MATERIAL ID 1 NAME "Baxant-Xi"
TYPE "CCModelBaXi94"
CONCRETE
CONCRETE TYPE 1
RATIO_WC 0.5
CEMENT_WEIGHT 0.27
TEMPERATURE
K_TEMP_TEMP 103680
C_TEMP_TEMP 0.000008
// Geometry definition
GEOMETRY ID 1 Name "Concrete column" TYPE "2D" thickness 10.
ELEMENT INCIDENCES
1 1 2 4 3
2 3 4 6 5
3 5 6 8 7
4 7 8 10 9
// initialisation
STEP ID 1 STATIC NAME "BCs and load" LOAD CASE 1 * 1.0 EXECUTE
OUTPUT LOCATION NODES DATA LIST "Q_T" "CURRENT_PSI_VALUES"
"EXTERNAL_FORCES" "INTERNAL_FORCES" "REACTIONS" END
ATENA Input File Format 305
// load alternative 2 - additional temperature increment induced by dT/dy and dQ at the top
and bottom
LOAD CASE ID 3 NAME "LC-2 -additional temperature increment" // total dT/dx=-40, i.e.
increment at fixed nodes -20, (incr. of Q already in LC-1)
SUPPORT SIMPLE SELECTION "all" dof 1 const 0. // fix h
SUPPORT SIMPLE NODE 3 DOF 2 VALUE 5
SUPPORT SIMPLE NODE 4 DOF 2 VALUE 5
SUPPORT SIMPLE NODE 5 DOF 2 VALUE 10
SUPPORT SIMPLE NODE 6 DOF 2 VALUE 10
SUPPORT SIMPLE NODE 7 DOF 2 VALUE 15
SUPPORT SIMPLE NODE 8 DOF 2 VALUE 15
LOAD SIMPLE SELECTION "all9-10" dof 2 const 1555200. // fix T
LOAD SIMPLE SELECTION "all1-2" dof 2 const -1555200. ; // fix T
// load alternative 3 - additional temperature increment induced by dT/dy and dQ at the top
and bottom
LOAD CASE ID 4 NAME "LC-2 -additional temperature increment" // total dT/dx=-40, i.e.
increment at fixed nodes -20, (incr. of Q already in LC-1)
SUPPORT SIMPLE SELECTION "all" dof 1 const 0. // fix h
306
//STEP ID 2 STATIC NAME "BCs and load" LOAD CASE 2 * 1.0 EXECUTE // step
execute command for the load alternative 1
//STEP ID 2 STATIC NAME "BCs and load" LOAD CASE 3 * 1.0 EXECUTE // step
execute command for the load alternative 2
STEP ID 2 STATIC NAME "BCs and load" LOAD CASE 4 * 1.0 EXECUTE // step
execute command for the load alternative 3
/* End of File */
ATENA Input File Format 307
3D......................................................................81
B
3DNONLINCEMENTITIOUS2FATIGUE....113
B 117, 119, 121, 126, 270
A BAND .............................................. 47, 141, 142
A 32, 33, 37, 81, 117, 119, 121, 126, 140, 209, BEAM_3D ......................................................... 81
221, 222, 235 BETA32, 83, 84, 86, 88, 90, 92, 94, 96, 98, 106,
ABSOLUTE..........................33, 37, 38, 288, 289 113, 115, 128, 130, 221, 222
AC147, 148, 150, 151, 152, 158, 159, 160, 162, BETA_FATIGUE .................................. 113, 116
163, 164 BODY ............ 180, 188, 189, 190, 195, 290, 292
AIR150, 151, 153, 158, 159, 160, 162, 163, 164,
165 C
ALL.47, 48, 49, 50, 203, 207, 255, 299, 300, 301
C 99, 107
ALPHA81, 82, 83, 84, 86, 89, 90, 92, 94, 96, 98,
106, 113, 115, 121, 126, 128, 129, 130, 131, C_1_X .............................................................. 59
135, 137, 138, 139, 142, 145 C_1_Y ........................................................ 59, 60
ALPHA_DP ............................................128, 129 C_1_Z......................................................... 59, 60
AND....................39, 40, 47, 49, 50, 51, 190, 200 C_2_X ........................................................ 59, 60
angle............................................................56, 69 C_2_Y ........................................................ 59, 60
ARC_LENGTH_CONSTANT .........................44 C_2_Z......................................................... 59, 60
ARC_LENGTH_PREVIOUS_STEP_LENGTH C159, 60, 121, 123, 124, 136, 141, 142, 147, 149
......................................................................42
C2 ..................... 60, 121, 123, 124, 137, 141, 142
ARC_LENGTH_RESET_STEP_LENGTH....42,
43 C3 ................................................................... 142
FACTOR33, 34, 37, 38, 56, 58, 83, 85, 86, 89, GF83, 84, 86, 87, 90, 91, 94, 95, 98, 106, 113,
90, 93, 94, 97, 98, 105, 106, 108, 112, 113, 114, 117, 118, 121, 122, 123, 124
116, 117, 120, 121, 125, 150, 151, 152, 158,
GIBBS-POOLE................................................ 47
159, 160, 162, 163, 164, 288, 289
GLOBAL.. 55, 189, 190, 203, 209, 291, 292, 299
FATIGUE_BASE_STRESS ...........................113
GROUP22, 24, 25, 26, 51, 62, 66, 67, 76, 188,
FATIGUE_COD_LOAD_COEFF..............48, 49
190, 198, 203, 207, 228, 238, 244, 245, 257,
FATIGUE_CYCLES ........................................48 258, 263, 265, 266, 268, 269, 270, 271, 272,
273, 290, 292, 298, 299, 304
FATIGUE_CYCLES_TO_FAILURE ............113
GROUPS .................................................. 51, 184
FATIGUE_MAX_FRACT_STRAIN.......48, 113
FATIGUE_MAX_FRACT_STRAIN_MULT.48, H
49
FATIGUE_PARAMS ...........................29, 30, 48 HARDENING ........................................ 126, 127
FATIGUE_TASK .............................................48 HISTORY13, 15, 16, 151, 153, 155, 156, 158,
159, 162, 189, 190, 218, 290
FC83, 86, 87, 88, 90, 91, 94, 95, 98, 99, 101,
106, 107, 110, 113, 114, 115, 117, 118, 120, HUMIDITY147, 150, 151, 152, 153, 155, 156,
121, 122, 124, 143, 145 158, 159, 160, 161, 162, 163, 164, 165, 218,
219
FC086, 88, 90, 91, 94, 95, 98, 106, 113, 115,
117, 118, 120, 121
I
FCYL28150, 151, 152, 153, 154, 155, 156, 157,
158, 159, 160, 161, 162, 163, 164 ID 51, 53, 54, 67, 68, 77, 78, 181, 189, 190, 198,
FILE........................................203, 206, 245, 246 200, 201, 202, 203, 209, 220, 225, 228, 236,
239, 240, 244, 245, 252, 261, 263, 264, 265,
FIXED56, 57, 85, 89, 92, 96, 105, 112, 115, 119, 271, 272, 273, 274, 276, 277, 286, 287, 298,
125, 188, 190, 201, 202, 232, 233 299, 303, 304, 305, 306
FRACTURE....................................211, 212, 215 IDS ................. 29, 30, 51, 54, 61, 63, 64, 65, 209
FRICTION ........................................56, 132, 133 IMPORT13, 15, 16, 189, 190, 203, 206, 218, 290
FROM22, 23, 25, 26, 28, 45, 46, 64, 65, 203, INCIDENCES66, 67, 75, 76, 210, 217, 228, 238,
206, 218, 231, 247, 270, 271, 272 291, 298, 304
FT83, 85, 86, 88, 89, 90, 91, 92, 94, 95, 96, 98, INCREMENENTAL LOAD.................. 291, 293
105, 106, 107, 112, 113, 114, 115, 117, 118,
119, 120, 121, 122, 123, 124, 132, 133 INCREMENT32, 77, 78, 201, 202, 221, 222,
224, 231, 232, 233, 288, 302
FULL_NR...................................................39, 40
INCREMENTAL LOAD ............................... 181
FUNCTION15, 56, 57, 98, 99, 101, 102, 103,
104, 106, 108, 110, 111, 132, 133, 135, 136, INERTIA_Y ..................................................... 59
137, 138, 139, 140, 141, 143, 146, 147, 166, INERTIA_Z ..................................................... 59
167, 168, 169, 170, 180, 181, 183, 187, 188,
198, 240, 241, 244, 245, 277, 291, 292, 293, INITIAL68, 69, 180, 188, 189, 190, 220, 252,
294, 298, 299 253, 286, 287
INPUT .............................................. 15, 245, 246
G INTERFACE...................... 79, 80, 130, 131, 132
GAMMA_COEFF ......................................68, 69 INTERNAL.... 203, 207, 214, 231, 299, 304, 306
GAMMA_REF ...........................................68, 69 INTERVALS.................................................. 218
GEOMETRY13, 15, 54, 55, 56, 59, 61, 64, 65, IP 13, 22, 25, 40, 41, 189, 190, 203, 208, 210,
67, 189, 203, 216, 218, 225, 236, 244, 276, 211, 291, 293
277, 279, 290, 298, 299, 304 IPS .............. 13, 22, 25, 64, 65, 66, 203, 207, 211
ITEM .............................................. 203, 207, 231
310
ITERATION33, 34, 37, 39, 40, 41, 46, 203, 204, LOADING_DISPLACEMENT_SCALE_CONS
205, 208 TANT ...........................................................45
LOCAL .....15, 53, 54, 55, 59, 189, 190, 291, 292
J
LOCATION41, 45, 46, 56, 57, 203, 207, 231,
JOINT15, 22, 26, 53, 76, 226, 236, 244, 245, 299, 300, 304, 306
253, 256, 257, 261, 262, 263, 264, 266, 267, LOSS...............................150, 151, 153, 156, 208
268, 269, 270, 272, 297, 299, 302
M
K
M 9, 10, 51, 52, 54, 59, 82, 84, 88, 92, 96, 105,
K 51, 52, 128, 129, 132, 134, 139, 140, 141, 111, 115, 119, 126, 128, 130, 131, 135, 137,
143, 146, 147, 166, 167, 168, 169, 170, 213, 138, 139, 142, 145, 211, 213, 215, 219, 273
218, 244, 274, 275, 277, 290, 300, 302, 303
MASTER180, 181, 182, 184, 186, 187, 198, 199,
K1 ................................... 126, 127, 141, 143, 144 210, 297, 299
K2 ................................... 126, 128, 141, 143, 144 MATERIAL13, 15, 61, 62, 64, 65, 66, 67, 78,
K3 ................................................... 141, 143, 144 79, 80, 81, 98, 107, 130, 131, 132, 149, 166,
167, 168, 169, 170, 171, 174, 175, 178, 179,
K4 ................................... 141, 142, 143, 144, 147 197, 203, 210, 212, 216, 217, 225, 226, 235,
KSI_FATIGUE....................................... 113, 116 236, 244, 245, 274, 275, 287, 298, 299, 301,
303, 304
L MATERIALS .....................................................62
MAXIMUM_ETA ............................................47
LIMIT_ETA ............................................... 46, 47
MESSAGE................................................15, 246
LINE .......................................................... 39, 40
MICROPLANE...........................79, 80, 140, 141
LINE, ........................................................ 39, 40
MINIMUM_ETA..............................................47
LINE_SEARCH_ITERATION_LIMIT ........... 46
MODIFIED_NR..........................................39, 40
LINE_SEARCH_WITH_ITERATIONS.......... 46
MODULUS.............................106, 108, 126, 127
LINE_SEARCH_WITHOUT_ITERATIONS . 46
MOISTURE ............150, 151, 153, 158, 159, 160
LINEAR29, 33, 39, 40, 68, 78, 79, 81, 82, 149,
150, 151, 154, 156, 158, 159, 160, 161, 162, MOMENT.........................................................59
163, 164, 165, 264, 265 MONITOR......203, 204, 205, 206, 208, 209, 231
LINE-SEARCH ................................................ 40 MU81, 82, 83, 86, 90, 94, 98, 106, 107, 113,
LIST22, 23, 26, 28, 64, 65, 203, 207, 231, 304, 121, 122, 126, 127, 128, 129, 131, 141, 144
306
LOAD13, 14, 15, 41, 45, 51, 76, 150, 151, 153,
N
155, 156, 158, 159, 160, 161, 162, 163, 164,
NAME15, 16, 17, 54, 67, 68, 78, 181, 198, 200,
165, 180, 181, 182, 183, 184, 185, 186, 187,
201, 202, 206, 209, 231, 240, 263, 265, 266,
188, 189, 190, 195, 198, 199, 200, 201, 202,
269, 270, 271, 272, 273, 298, 299, 303, 304,
203, 208, 210, 220, 222, 229, 230, 232, 233,
305, 306
239, 240, 244, 245, 252, 254, 286, 290, 291,
292, 293, 294, 299, 300, 301, 302, 304, 305, NEWMARK ...........32, 33, 50, 51, 221, 222, 231
306
NEWTON-RAPHSON ...............................39, 40
LOAD_DISPLACEMENT_RATIO........... 41, 45
NODAL15, 47, 48, 49, 50, 76, 77, 78, 184, 186,
LOADING_DISPLACEMENT_BERGAN_CON 198, 203, 214, 217, 220, 221, 230, 252, 263,
STANT......................................................... 45 265, 286, 287, 299, 300, 301, 303
LOADING_DISPLACEMENT_RATIO_CONS
TANT........................................................... 45
ATENA Input File Format 311
NODE13, 14, 22, 24, 26, 45, 46, 77, 78, 187, Q
189, 190, 197, 198, 203, 207, 220, 221, 252,
269, 270, 286, 287, 301, 303, 305, 306 Q1 ................................................................... 147
NODES13, 22, 24, 25, 26, 51, 180, 181, 182,
Q2 ................................................................... 148
183, 184, 185, 186, 187, 199, 203, 207, 210,
212, 214, 217, 231, 257, 258, 259, 260, 263, Q3 ................................................................... 148
264, 265, 266, 267, 269, 270, 271, 272, 273,
Q4 ................................................................... 148
291, 293, 294, 295, 298, 300, 304, 306
NOMINAL_HC ................................................291 R
NONE .......................47, 133, 203, 206, 257, 258
R_C83, 86, 88, 90, 91, 94, 95, 98, 99, 106, 107,
NONLINEAR ...................................................68
113, 114, 115, 117, 118, 121, 122
NORMAL_UPDATE .......................................41
R_C086, 88, 90, 91, 94, 95, 98, 106, 113, 115,
NP ...................................................................141 117, 118
NUMBER .......44, 64, 65, 66, 147, 218, 223, 240 R_T83, 86, 90, 94, 98, 106, 107, 113, 114, 117,
121, 122, 132
O RADIUS ........................................................... 56
RATIO 41, 45, 137, 138, 139, 166, 274, 276, 303
OFF .....................................29, 30, 203, 207, 208
RC42, 44, 83, 86, 90, 91, 94, 98, 99, 106, 107,
ON.............................30, 108, 180, 203, 207, 208
113, 114, 117, 118, 121, 122
OPTIMIZE............................................29, 30, 47
RC086, 88, 90, 91, 94, 95, 98, 106, 113, 115,
OUTPUT13, 15, 16, 189, 190, 195, 202, 203, 117, 118
207, 208, 209, 210, 211, 212, 214, 216, 217,
REFERENCE_DLAMBDA....................... 42, 43
231, 256, 291, 292, 294, 295, 299, 300, 304,
306 REFERENCE_ETA ......................... 46, 288, 302
REFERENCE_NUMBER_OF_ITERATIONS 44
P
REGION ........................................................... 13
PASTERNAK...................................................59 REINFORCEMENT15, 61, 64, 65, 79, 80, 135,
PATCH ............................................................254 136, 137, 138, 209, 217, 262, 270, 276
PLANE_STRAIN81, 82, 83, 85, 86, 89, 90, 93, RELAX .......................................................... 181
94, 98, 106, 112, 113, 116, 117, 120, 126, REMOVE......... 22, 24, 26, 28, 45, 203, 206, 255
128, 130, 131, 141, 142, 143, 146
RESIDUAL ............ 34, 37, 38, 39, 214, 288, 289
PLANE_STRESS81, 82, 83, 85, 86, 89, 90, 93,
RESTORE .......................................... 15, 16, 247
94, 98, 106, 112, 113, 116, 117, 120, 126,
128, 130, 131, 141, 142, 143, 146 RETARD_TIMES_PER_DECADE .. 49, 50, 218
PLASTIC ................................211, 212, 213, 215 RETENTION ......................................... 121, 125
POINTS ..............................49, 50, 203, 207, 299 RHO81, 82, 83, 84, 86, 88, 90, 92, 94, 96, 98,
106, 113, 115, 121, 126, 128, 130, 131, 135,
POISSON81, 82, 83, 86, 90, 94, 98, 106, 107,
137, 138, 139, 142, 145
113, 121, 122, 126, 127, 128, 129, 131, 141,
144 RT83, 86, 90, 94, 98, 106, 107, 113, 114, 117,
121, 122, 132, 218, 290
POLAR .............................................................59
PRAGER.............................79, 80, 126, 128, 129 S
PRESTRESSING....................................188, 190
SAMPLE_TIMES_PER_DECADE........... 49, 50
PSI.............................77, 220, 252, 253, 286, 287
SBETA ........................................... 211, 212, 214
312
SECANT_PREDICTOR .................................. 40 128, 130, 131, 141, 142, 143, 146, 180, 188,
189, 190, 198, 211, 212, 213, 214, 215
SERIALIZE ........................ 29, 30, 47, 49, 50, 51
SUPPORT181, 210, 229, 239, 299, 301, 302,
SET15, 16, 22, 23, 28, 29, 30, 33, 38, 41, 209,
304, 305, 306
221, 223, 224, 231, 240, 288, 289, 299, 302
SURFACE........................................................254
SHAPE150, 151, 152, 158, 159, 160, 162, 163,
164, 184, 264, 265 switches
SHEAR59, 83, 85, 86, 89, 90, 93, 94, 97, 98, /batch_execute ..............................................11
103, 104, 105, 106, 108, 112, 113, 116, 117,
/execute.....................................................9, 11
120, 121, 125
/silent ............................................................11
SHEAR_Y ........................................................ 59
SHEAR_Z......................................................... 59 T
SHELL .............................................................. 81
T 98, 107
SHRINKAGE150, 151, 153, 155, 156, 158, 159,
160, 161, 162, 164, 165, 209 TASK15, 17, 48, 53, 55, 209, 224, 235, 261,
297, 302
SIMPLE180, 181, 182, 183, 229, 230, 239, 299,
301, 302, 304, 305, 306 TEMPERATURE51, 147, 151, 153, 155, 156,
158, 159, 162, 180, 188, 189, 190, 198, 211,
SLAVE180, 181, 182, 183, 184, 186, 187, 198,
213, 215, 218, 219, 274, 275, 291, 293, 303
199, 210, 256, 297, 299
TENSILE ........................................211, 213, 215
SLOAN............................................................. 47
THICKNESS54, 55, 61, 62, 64, 65, 150, 151,
SMEARED ......................... 79, 80, 135, 137, 138
153, 154, 156, 158, 160, 161, 162, 163, 164,
SOLVER_KEYS ............................ 29, 30, 49, 50 226, 236, 276, 284, 285, 286
SPRING54, 55, 79, 80, 139, 140, 166, 167, 169, TIME32, 49, 50, 51, 52, 150, 151, 153, 155, 156,
170, 174, 178, 180, 181, 197, 216 158, 159, 160, 161, 162, 163, 164, 165, 189,
190, 209, 218, 221, 222, 224, 231, 288, 291,
STANDARD .................. 47, 48, 49, 50, 203, 206
292, 293, 294, 302
STATIC32, 200, 201, 202, 224, 240, 299, 304,
TIME_INTEGRATION....32, 222, 223, 288, 302
306
TIMES...............................................................13
STEAM150, 151, 153, 158, 159, 160, 162, 163,
164, 165 TITLE .......................................17, 209, 297, 302
STEP13, 15, 16, 37, 38, 40, 41, 42, 43, 49, 200, TO22, 23, 26, 28, 45, 46, 64, 65, 113, 188, 190,
201, 202, 203, 204, 205, 206, 208, 209, 210, 200, 203, 206, 207, 218, 231, 247, 270, 272,
217, 222, 224, 231, 232, 233, 240, 244, 245, 290, 292, 306
288, 289, 299, 304, 306
TORGUE ..........................................................59
STEP_LENGTH......................................... 42, 43
TOTAL LOAD ...............................181, 291, 293
STEPS............................................................... 51
TOTAL_LOSS................................150, 151, 153
STOP_TIME....................... 49, 50, 221, 222, 231
TRACE ...........................................203, 207, 208
STORE ....................................... 15, 16, 247, 248
TRANSIENT32, 33, 221, 223, 224, 231, 288,
STRAIN48, 49, 68, 69, 81, 82, 83, 85, 86, 89, 302
90, 93, 94, 97, 98, 106, 112, 113, 116, 117,
TYPE13, 23, 28, 29, 32, 33, 40, 41, 54, 66, 67,
120, 126, 128, 130, 131, 141, 142, 143, 146,
68, 77, 78, 81, 83, 86, 90, 94, 98, 106, 107,
180, 188, 189, 190, 198, 211, 212, 213, 214,
113, 117, 121, 126, 128, 131, 132, 135, 136,
215, 295
137, 138, 139, 141, 143, 146, 147, 149, 166,
STRENGTH98, 104, 106, 111, 126, 127, 211, 167, 169, 170, 174, 178, 181, 183, 200, 201,
213, 215 202, 203, 216, 220, 221, 222, 223, 225, 228,
232, 233, 236, 239, 240, 244, 252, 263, 265,
STRESS68, 69, 81, 82, 83, 85, 86, 89, 90, 93, 94,
97, 98, 106, 112, 113, 116, 117, 120, 126,
ATENA Input File Format 313
V Y
VALUE13, 77, 181, 182, 183, 187, 189, 190, Y 15, 54, 59, 60, 69, 71, 77, 78, 81, 82, 83, 86,
198, 254, 291, 292, 295, 299, 305, 306 90, 94, 98, 106, 113, 121, 122, 126, 127, 128,
129, 131, 133, 141, 144, 180, 187, 189, 190,
VARIABLE44, 79, 81, 125, 166, 167, 169, 170,
198, 220, 221, 252, 253, 287, 291, 292, 302
174, 178
YIELD............................ 126, 127, 212, 213, 215
VARIATIONAL ..............................................208
YVALUES ............................................. 241, 298
W YX .......................................................... 189, 198
YZ .......................................................... 189, 198
WATER150, 151, 153, 158, 159, 160, 162, 163,
164, 165, 208
Z
WC147, 148, 150, 151, 152, 158, 159, 160, 162,
163, 164, 274, 276, 303
Z 59, 60, 77, 78, 187, 189, 190, 198, 220, 221,
WD83, 84, 86, 88, 90, 92, 94, 95, 98, 106, 113, 252, 253, 287, 291, 292
115, 121, 125, 128, 129
ZX .......................................................... 189, 198
WIDTH .............................................................47
ZY .......................................................... 189, 198
WINKLER ........................................................59