Code_Aster ®
Version
7.4
Titrate:
Macro-command CALC_PRECONT


Date:
01/02/05
Author (S):
S. MICHEL-PONNELLE, A. ASSIRE Key
:
U4.42.05-A Page
: 1/10

Organization (S): EDF-R & D/AMA

Handbook of Utilization
U4.4- booklet: Modeling
Document: U4.42.05

Macro-command CALC_PRECONT

1 Goal

This macro-command makes it possible to put in tension the cables in a structure (cf [R7.01.02]), of
such kind that, at the end of this calculation, the structure either balances some and the tension or equal to the tension
data by rules BPEL91, calculated by command DEFI_CABLE_BP.
It also allows:

·
to apply prestressing in a progressive way, in order to be able to treat the case where concrete
will plasticize or damage themselves according to the model of behavior selected,
·
to practice the phasage of setting in prestressing, i.e. the setting in tension of the cables
in a sequential way.
Handbook of Utilization
U4.4- booklet: Modeling
HT-66/05/004/A

Code_Aster ®
Version
7.4
Titrate:
Macro-command CALC_PRECONT


Date:
01/02/05
Author (S):
S. MICHEL-PONNELLE, A. ASSIRE Key
:
U4.42.05-A Page
: 2/10

Count

matters

1 Goal ......................................................................................................................................................... 1
2 Syntax .................................................................................................................................................. 3
3 Operands ............................................................................................................................................. 4
3.1 Operands MODELS/CHAM_MATER/CARA_ELEM .......................................................................... 4
3.2 Key word EXCIT ................................................................................................................................. 4
3.2.1 Operands CHARGES ................................................................................................................ 4
3.3 Key word COMP_INCR ......................................................................................................................... 4
3.4 Key word INCREMENT ......................................................................................................................... 4
3.4.1 Operands LIST_INST ......................................................................................................... 4
3.4.2 Operands INST_INIT/INST_FIN ................................................................................. 5
3.4.3 Operand PRECISION ........................................................................................................... 5
3.4.4 Operand SUBD_PAS/SUBD_PAS_MINI/COEF_SUBD_PAS_1 ................................... 5
3.5 Key word CABLE_BP ........................................................................................................................... 5
3.6 Key word CABLE_BP_INACTIF .......................................................................................................... 5
3.7 Key word NEWTON ............................................................................................................................... 6
3.8 Key word ETAT_INIT ......................................................................................................................... 6
3.9 Key word RECH_LINEAIRE ................................................................................................................ 6
3.10
Operand PARM_THETA ......................................................................................................... 6
3.11
Key word SOLVEUR .................................................................................................................... 6
3.12
Key word CONVERGENCE ........................................................................................................... 6
3.13
Operand SOLV_NON_LOCAL ................................................................................................ 6
3.14
Operand LAGR_NON_LOCAL ................................................................................................ 6
3.15
Operand INFORMATION ...................................................................................................................... 6
3.16
Operand TITRATES .................................................................................................................... 6
4 Example ................................................................................................................................................. 7
4.1 Before using the macro-command ............................................................................................... 7
4.2 Use of the macro-command: setting in successive tension of the 5 cables ............................. 9
4.3 Continuation of the loading after setting in tension of the cables ........................................................ 10

Handbook of Utilization
U4.4- booklet: Modeling
HT-66/05/004/A

Code_Aster ®
Version
7.4
Titrate:
Macro-command CALC_PRECONT


Date:
01/02/05
Author (S):
S. MICHEL-PONNELLE, A. ASSIRE Key
:
U4.42.05-A Page
: 3/10

2 Syntax

statnl [evol_noli] = CALC_PRECONT

(
reuse
=
statnl,

MODELE
=
Mo
,
[model]
CHAM_MATER

=
chmat
,
[cham_mater]

CARA_ELEM


=
carac
,
[cara_elem]

EXCIT

= (_F (

CHARGE
=
chi
/
[char_meca]
















/[char_cine_meca]






),),
COMP_INCR = (see the document [U4.51.03]
)

INCREMENT =_F (
LIST_INST
=
litps
,


[listr8]

INST_INIT
=
instini,
[R]
INST_FIN
=
instfin,
[R]






PRECISION
=/1.0E-3
, [DEFAUT]
/
prec,
[R]


SUBD_PAS
=
/
1,
[DEFAUT]

/subpas
, [I]

SUBD_PAS_MINI
=
submini,
[R]

COEF_SUBD_PAS_1
=/1.
, [DEFAUT]
/coefsub, [R]
),

CABLE_BP



=
l_cabl_pr
,


[cable_precont]

CABLE_BP_INACTIF
=
l_cabl_pr
,


[cable_precont]

ETAT_INIT
= (see the document [U4.51.03]
),

NEWTON
=
(
to see the document [U4.51.03]
),
RECH_LINEAIRE
= (see the document [U4.51.03]
),


PARM_THETA
= (see the document [U4.51.03]
),



SOLVEUR

=
(
to see the document [U4.50.01]),


CONVERGENCE
= (see the document [U4.51.03]
),

LAGR_NON_LOCAL = (see the document [U4.51.03]
),

SOLV_NON_LOCAL = (see the document [U4.50.01]
),


INFO =
/1,
[DEFAUT]





/2,


TITER
=
tx [kN]

)
Handbook of Utilization
U4.4- booklet: Modeling
HT-66/05/004/A

Code_Aster ®
Version
7.4
Titrate:
Macro-command CALC_PRECONT


Date:
01/02/05
Author (S):
S. MICHEL-PONNELLE, A. ASSIRE Key
:
U4.42.05-A Page
: 4/10

3 Operands

3.1 Operands
MODEL/CHAM_MATER/CARA_ELEM

MODELE = Mo

Name of the model whose elements are the subject of mechanical calculation.

CHAM_MATER = chmat

Name of the affected material field on the grid. Attention, all meshs of the model
must be associated a material.

CARA_ELEM = carac

Name of the characteristics of the elements of hull, beam, pipe, bars, discrete cable and elements
affected on the model Mo.

3.2 Word
key
EXCIT

EXCIT =

This key word factor makes it possible to describe a load with each occurrence. It is necessary to provide them
boundary conditions for the structure, possibly of the instantaneous loadings like
gravity as well as the connections kinematics related to the cables having already been put in tension by
a preceding call to CALC_PRECONT. To in no case, one should not include the loading of the cable
that one wants to put in tension by the call to this macro-command.

3.2.1 Operands
CHARGE

CHARGE = chi

CH is the mechanical loading specified with the ième occurrence of
I
EXCIT.

3.3 Word
key
COMP_INCR

The syntax of this key word common to command STAT_NON_LINE is described in the document
[U4.51.03].

3.4 Word
key
INCREMENT

INCREMENT =

Defines the intervals of time taken in the incremental method.

3.4.1 Operands
LIST_INST

LIST_INST = litps

The moments of calculations are those defined in the concept litps by the operator
DEFI_LIST_REEL [U4.34.01]. This list must be ordered in a chronological way (increasing).

Note:

Even if calculation is carried out with several steps of time, only the last step of time is
stored in the concept result.
During the process of setting in tension of the cables, moments are generated
automatically in addition to those provided by the user. It is thus completely normal to see
to appear in the file of message of calculations with STAT_NON_LINE at moments that
the user did not specify. These moments are not stored in the concept result.
Handbook of Utilization
U4.4- booklet: Modeling
HT-66/05/004/A

Code_Aster ®
Version
7.4
Titrate:
Macro-command CALC_PRECONT


Date:
01/02/05
Author (S):
S. MICHEL-PONNELLE, A. ASSIRE Key
:
U4.42.05-A Page
: 5/10

3.4.2 Operands
INST_INIT/INST_FIN



INST_INIT

= instini

The initial moment of the calculation (which thus (Re) is not calculated) must be indicated by its value.

If this initial moment is not recomputed, the defect is calculated in the following way:

·
if an initial state is specified (operand ETAT_INIT) and if it definite one moment
corresponding (by EVOL_NOLI or INST_ETAT_INIT) then the initial moment is that
defined by the initial state,
·
if there is no initial state (operand ETAT_INIT) or that it does not define a moment
corresponding (the fields are given in ETAT_INIT without specifying
INST_ETAT_INIT), then one takes the first moment of the list of moments litps
(NUME_INST_INIT =0).



INST_FIN
=

instfin

The final moment (last calculated step) is indicated same manner as the initial moment, except
that it is not possible to refer to the moment of the initial state.

3.4.3 Operand
PRECISION

PRECISION
=
prec (Cf. [U4.71.00])

3.4.4 Operand
SUBD_PAS/SUBD_PAS_MINI/COEF_SUBD_PAS_1

SUBD_PAS = subpas,
SUBD_PAS_MINI = submini,
COEF_SUBD_PAS_1 = coefsub,

Allows to carry out an automatic recutting of the step of time when the algorithm of Newton
do not converge.
The step of time is redécoupé in subpas under step. By defect there is no recutting
(SUBD_PAS = 1). The automatic subdivision stops when the new steps created are more
small that SUBD_PAS_MINI. The new steps created are of identical size, except the first
who is equal to this size multiplied by COEF_SUBD_PAS_1 (by defect 1). This allows best
to take into account the problems of discharge of the structure (change of tangent matrix)
without using the elastic matrix (PREDICTION = ' ELASTIQUE' or MATRICE = “ELASTIQUE”
under operand NEWTON).

3.5 Word
key
CABLE_BP

CABLE_BP = l_cabl_pr

It is a question here of providing a list of the concepts of the cabl_precont type produced by the operator
DEFI_CABLE_BP [U4.42.04]. All the cables concerned will be tended at the end of this calculation.

3.6 Word
key
CABLE_BP_INACTIF


CABLE_BP_INACTIF = l_cabl_pr

It is a question here of providing a list of the concepts of the cabl_precont type produced by the operator
DEFI_CABLE_BP [U4.42.04]. The macro-command is given the responsability to generate the connections
kinematics related to these inactive cables, and does not take into account the rigidity of these cables.
Handbook of Utilization
U4.4- booklet: Modeling
HT-66/05/004/A

Code_Aster ®
Version
7.4
Titrate:
Macro-command CALC_PRECONT


Date:
01/02/05
Author (S):
S. MICHEL-PONNELLE, A. ASSIRE Key
:
U4.42.05-A Page
: 6/10

3.7 Word
key
NEWTON

The syntax of this key word common to command STAT_NON_LINE is described in the document
[U4.51.03].

3.8 Word
key
ETAT_INIT

The syntax of this key word common to command STAT_NON_LINE is described in the document
[U4.51.03].

3.9 Word
key
RECH_LINEAIRE

The syntax of this key word common to command STAT_NON_LINE is described in the document
[U4.51.03].

3.10 Operand
PARM_THETA

The syntax of this key word common to command STAT_NON_LINE is described in the document
[U4.51.03].

3.11 Word
key
SOLVEUR

The syntax of this key word common to several commands is described in the document [U4.50.01].

3.12 Word
key
CONVERGENCE

The syntax of this key word common to command STAT_NON_LINE is described in the document
[U4.51.03].

3.13 Operand
SOLV_NON_LOCAL

The syntax of this key word is identical to key word SOLVEUR describes in the document [U4.50.01]. With
to use for a nonlocal model.

3.14 Operand
LAGR_NON_LOCAL

The syntax of this key word common to command STAT_NON_LINE is described in the document
[U4.51.03].

3.15 Operand
INFO

The syntax of this key word common to command STAT_NON_LINE is described in the document
[U4.51.03].

3.16 Operand
TITER

TITER = tx

tx is the title of calculation. It will be printed at the head results. See [U4.03.01].
Handbook of Utilization
U4.4- booklet: Modeling
HT-66/05/004/A

Code_Aster ®
Version
7.4
Titrate:
Macro-command CALC_PRECONT


Date:
01/02/05
Author (S):
S. MICHEL-PONNELLE, A. ASSIRE Key
:
U4.42.05-A Page
: 7/10

4 Example

One details here the principal stages of a structural analysis containing of the cables of prestressing
(case-test SSNV164)

4.1
Before using the macro-command

·
Definition of the cables

In the event of setting in not-simultaneous tension of all the cables, it is necessary to make (at least)
as many DEFI_CABLE_BP as of stages of loading.
CAB_BP1=DEFI_CABLE_BP (MODELE=MO,
CHAM_MATER=CMAT,
CARA_ELEM=CE,
GROUP_MA_BETON=' VOLTOT',
TYPE_ANCRAGE= (“ACTIVE”, “PASSIVE”,),
DEFI_CABLE= (_F (GROUP_MA=' CAB1',
GROUP_NO_ANCRAGE= (“PC1D', “PC1F”,),),
_F (GROUP_MA=' CAB2',
GROUP_NO_ANCRAGE= (“PC 2D”, “PC2F”,),),),
TENSION_INIT=3.750000E6,
RECUL_ANCRAGE=0.001,)

CAB_BP3=DEFI_CABLE_BP (MODELE=MO,
CHAM_MATER=CMAT,
CARA_ELEM=CE,
GROUP_MA_BETON=' VOLTOT',
TYPE_ANCRAGE= (“ACTIVE”, “PASSIVE”,),
DEFI_CABLE= (_F (GROUP_MA=' CAB3',
GROUP_NO_ANCRAGE= (“PC 3D”, “PC3F”),),
_F (GROUP_MA=' CAB4',
GROUP_NO_ANCRAGE= (“PC4D', “PC4F”),),),
TENSION_INIT=3.750000E6,
RECUL_ANCRAGE=0.001,)

CAB_BP5=DEFI_CABLE_BP (MODELE=MO,
CHAM_MATER=CMAT,
CARA_ELEM=CE,
GROUP_MA_BETON=' VOLTOT',
TYPE_ANCRAGE= (“ACTIVE”, “ACTIVE”,),
DEFI_CABLE=_F (GROUP_MA=' CAB5',
GROUP_NO_ANCRAGE= (“PC5D', “PC5F”,),),
TENSION_INIT=3.750000E6,
RECUL_ANCRAGE=0.001,
)
Handbook of Utilization
U4.4- booklet: Modeling
HT-66/05/004/A

Code_Aster ®
Version
7.4
Titrate:
Macro-command CALC_PRECONT


Date:
01/02/05
Author (S):
S. MICHEL-PONNELLE, A. ASSIRE Key
:
U4.42.05-A Page
: 8/10

·
Definition of the loadings

It is necessary to distinguish:
-
loadings related to the boundary conditions plus the possible instantaneous loadings,
-
the loading related to the cables containing only the connections kinematics,
-
posterior not-instantaneous loadings with the setting in tension of the cables.

CLIM =AFFE_CHAR_MECA (MODELE=MO,
DDL_IMPO= (
_F (GROUP_NO=' PP',
DX=0.0, DY=0.0,),
_F (GROUP_NO=' PX',
DY=0.0,),
_F (GROUP_NO=' PY',
DX=0.0,),
_F (GROUP_NO=' SU3',
DZ=0.0,),),
PESANTEUR= (9.8100000000000005, 0.0, 0.0, - 1.0,),)

CMCAB1=AFFE_CHAR_MECA (MODELE=MO,
RELA_CINE_BP=_F (CABLE_BP=CAB_BP1,
SIGM_BPEL=' NON',
RELA_CINE=' OUI',),)
CMCAB3=AFFE_CHAR_MECA (MODELE=MO,
RELA_CINE_BP=_F (CABLE_BP=CAB_BP3,
SIGM_BPEL=' NON',
RELA_CINE=' OUI',),)
CMCAB5=AFFE_CHAR_MECA (MODELE=MO,
RELA_CINE_BP=_F (CABLE_BP=CAB_BP5,
SIGM_BPEL=' NON',
RELA_CINE=' OUI',),)

CHMECA =AFFE_CHAR_MECA (MODELE=MO,
DDL_IMPO=_F (GROUP_NO=' SU2',
DZ=1.0,),)

·
Preliminary calculation before setting in tension of the cables (optional)

It is a question here of taking into account the loadings which apply to the structure before the setting in
tension of the cables. Not to take into account the rigidity of the cables, and not to introduce one
loading on the cables, it is necessary to carry out calculation is on a model not containing the cables
maybe in their affecting a law of behavior SANS, which imposes that the constraint remains null in these
elements.

RES1 = STAT_NON_LINE (MODELE=MO,
CHAM_MATER=CMAT,
CARA_ELEM=CE,
COMP_INCR= (_F (RELATION = “ELAS”,),
_F (RELATION = “WITHOUT”,
GROUP_MA= (“CAB1”, “CAB3”, “CAB5”),),),
EXCIT = (_F (LOAD = CLIM,),
_F (CHARGE = CMCAB1),
_F (CHARGE = CMCAB3),
_F (CHARGE = CMCAB5),),
INCREMENT=_F (LIST_INST = LINST,
INST_FIN = 150.),)
Handbook of Utilization
U4.4- booklet: Modeling
HT-66/05/004/A

Code_Aster ®
Version
7.4
Titrate:
Macro-command CALC_PRECONT


Date:
01/02/05
Author (S):
S. MICHEL-PONNELLE, A. ASSIRE Key
:
U4.42.05-A Page
: 9/10

4.2
Use of the macro-command: setting in successive tension of
5 cables

One calls upon the macro-command as many once as necessary.
RES1 = CALC_PRECONT (reuse=RES1,
ETAT_INIT=_F (EVOL_NOLI=RES1),
MODELE=MO,
CHAM_MATER=CMAT,
CARA_ELEM=CE,
COMP_INCR= (_F (RELATION = “ELAS”,),),
EXCIT = (_F (LOAD = CLIM,),),
CABLE_BP = (CAB_BP1,),
CABLE_BP_INACTIF = (CAB_BP3, CAB_BP5,),
INCREMENT=_F (LIST_INST = LINST,
INST_FIN = 300.,
SUBD_PAS = 4,
SUBD_PAS_MINI = 0.01,),)
RES1 = CALC_PRECONT (reuse=RES1,
ETAT_INIT=_F (EVOL_NOLI=RES1),
MODELE=MO,
CHAM_MATER=CMAT,
CARA_ELEM=CE,
COMP_INCR= (_F (RELATION = “ELAS”,),),
EXCIT = (_F (LOAD = CLIM,),
_F (CHARGE = CMCAB1,),
CABLE_BP = (CAB_BP3,),
CABLE_BP_INACTIF = (CAB_BP5,),
INCREMENT=_F (LIST_INST = LINST,
INST_FIN = 450.,
SUBD_PAS = 4,
SUBD_PAS_MINI = 0.01,),)
RES1 = CALC_PRECONT (reuse=RES1,
ETAT_INIT=_F (EVOL_NOLI=RES1),
MODELE=MO,
CHAM_MATER=CMAT,
CARA_ELEM=CE,
COMP_INCR= (_F (RELATION = “ELAS”,),),
EXCIT = (_F (LOAD = CLIM,),
_F (CHARGE = CMCAB1,),
_F (CHARGE = CMCAB3,),
CABLE_BP = (CAB_BP5,),
INCREMENT=_F (LIST_INST = LINST,
INST_FIN = 600.,
SUBD_PAS = 4,
SUBD_PAS_MINI = 0.01,),)
Handbook of Utilization
U4.4- booklet: Modeling
HT-66/05/004/A

Code_Aster ®
Version
7.4
Titrate:
Macro-command CALC_PRECONT


Date:
01/02/05
Author (S):
S. MICHEL-PONNELLE, A. ASSIRE Key
:
U4.42.05-A Page
: 10/10

4.3
Continuation of the loading after setting in tension of the cables

The cables being tended, there is not any more but to continue calculation by always including the connections
kinematics binding the nodes of the cable to the concrete.

RES1 = STAT_NON_LINE (reuse = RES1,
ETAT_INIT=_F (EVOL_NOLI=RES1),
MODELE=MO,
CHAM_MATER=CMAT,
CARA_ELEM=CE,
COMP_INCR= _F (RELATION = “ELAS”,),
EXCIT = (_F (LOAD = CLIM,),
_F (CHARGE = CMCAB1),
_F (CHARGE = CMCAB3),
_F (CHARGE = CMCAB5),
_F (CHARGE=CHMECA,
FONC_MULT = FCT,),
INCREMENT=_F (LIST_INST = LINST,),)

Handbook of Utilization
U4.4- booklet: Modeling
HT-66/05/004/A

Outline document