background image
Code_Aster
®
Version
7.4
Titrate:
Macro-control
CALC_PRECONT
Date:
01/02/05
Author (S):
S. MICHEL-PONNELLE, A. ASSIRE
Key
:
U4.42.05-A
Page
:
1/10
Instruction manual
U4.4- booklet: Modeling
HT-66/05/004/A
Organization (S):
EDF-R & D/AMA














Instruction manual
U4.4- booklet: Modeling
Document: U4.42.05



Macro-control
CALC_PRECONT








1 Goal
This macro-control makes it possible to put in voltage 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 voltage or equal to the voltage
data by rules BPEL91, calculated by control 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 voltage of the cables
in a sequential way.
background image
Code_Aster
®
Version
7.4
Titrate:
Macro-control
CALC_PRECONT
Date:
01/02/05
Author (S):
S. MICHEL-PONNELLE, A. ASSIRE
Key
:
U4.42.05-A
Page
:
2/10
Instruction manual
U4.4- booklet: Modeling
HT-66/05/004/A
Count
matters
1
......................................................................................................................................................... 1 drank
2
Syntax .................................................................................................................................................. 3
3
Operands ............................................................................................................................................. 4
3.1
Operands
MODEL
/
CHAM_MATER
/
CARA_ELEM
.......................................................................... 4
3.2
Key word
EXCIT
................................................................................................................................. 4
3.2.1
Operands
CHARGE
................................................................................................................ 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
TITRATE
.................................................................................................................... 6
4
Example ................................................................................................................................................. 7
4.1
Before using the macro-control ............................................................................................... 7
4.2
Use of the macro-control: setting in successive voltage of the 5 cables ............................. 9
4.3
Continuation of the loading after setting in voltage of the cables ........................................................ 10
background image
Code_Aster
®
Version
7.4
Titrate:
Macro-control
CALC_PRECONT
Date:
01/02/05
Author (S):
S. MICHEL-PONNELLE, A. ASSIRE
Key
:
U4.42.05-A
Page
:
3/10
Instruction manual
U4.4- booklet: Modeling
HT-66/05/004/A
2 Syntax
statnl [evol_noli] = CALC_PRECONT
(
reuse
=
statnl,
MODEL
=
Mo
,
[model]
CHAM_MATER
=
chmat
,
[cham_mater]
CARA_ELEM
=
carac
,
[cara_elem]
EXCIT
= (_F (
CHARGE
=
chi
/
[char_meca]
/[char_cine_meca]
),),
COMP_INCR = (
to see the document [U4.51.03]
)
INCREMENT =_F (
LIST_INST
=
litps
,
[listr8]
INST_INIT
=
instini,
[R]
INST_FIN
=
instfin,
[R]
PRECISION
=/1.0E-3
, [DEFECT]
/
prec,
[R]
SUBD_PAS
=
/
1,
[DEFECT]
/subpas
, [I]
SUBD_PAS_MINI
=
submini,
[R]
COEF_SUBD_PAS_1
=/1.
, [DEFECT]
/coefsub, [R]
),
CABLE_BP
=
l_cabl_pr
,
[cable_precont]

CABLE_BP_INACTIF
=
l_cabl_pr
,
[cable_precont]
ETAT_INIT
= (
to see the document [U4.51.03]
),
NEWTON
=
(
to see the document [U4.51.03]
),
RECH_LINEAIRE
= (
to see the document [U4.51.03]
),
PARM_THETA
= (
to see the document [U4.51.03]
),
SOLVEUR
=
(
to see the document [U4.50.01]
),
CONVERGENCE
= (
to see the document [U4.51.03]
),
LAGR_NON_LOCAL = (
to see the document [U4.51.03]
),
SOLV_NON_LOCAL = (
to see the document [U4.50.01]
),
INFORMATION =
/1,
[DEFECT]
/2,
TITRATE
=
tx [kN]
)
background image
Code_Aster
®
Version
7.4
Titrate:
Macro-control
CALC_PRECONT
Date:
01/02/05
Author (S):
S. MICHEL-PONNELLE, A. ASSIRE
Key
:
U4.42.05-A
Page
:
4/10
Instruction manual
U4.4- booklet: Modeling
HT-66/05/004/A
3 Operands
3.1 Operands
MODEL
/
CHAM_MATER
/
CARA_ELEM
MODEL = Mo
Name of the model whose elements are the subject of mechanical calculation.
CHAM_MATER = chmat
Name of the affected material field on the mesh. 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 voltage 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 voltage by the call to this macro-control.
3.2.1 Operands
CHARGE
CHARGE = CH
I
CH
I
is the mechanical loading specified with I
ème
occurrence of
EXCIT
.
3.3 Word
key
COMP_INCR
The syntax of this key word common to the control
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 pitches of time, only the last pitch of time is
stored in the concept result.
During the process of setting in voltage 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.
background image
Code_Aster
®
Version
7.4
Titrate:
Macro-control
CALC_PRECONT
Date:
01/02/05
Author (S):
S. MICHEL-PONNELLE, A. ASSIRE
Key
:
U4.42.05-A
Page
:
5/10
Instruction manual
U4.4- booklet: Modeling
HT-66/05/004/A
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 pitch) 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 pitch of time when the algorithm of Newton
do not converge.
The pitch of time is redécoupé in
subpas
under pitch. By defect there is no recutting
(SUBD_PAS = 1). The automatic subdivision stops when the new pitches created are more
small that
SUBD_PAS_MINI
. The new pitches 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
STAMP = “ELASTIC”
under the 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-control 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.
background image
Code_Aster
®
Version
7.4
Titrate:
Macro-control
CALC_PRECONT
Date:
01/02/05
Author (S):
S. MICHEL-PONNELLE, A. ASSIRE
Key
:
U4.42.05-A
Page
:
6/10
Instruction manual
U4.4- booklet: Modeling
HT-66/05/004/A
3.7 Word
key
NEWTON
The syntax of this key word common to the control
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 the control
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 the control
STAT_NON_LINE
is described in the document
[U4.51.03].
3.10 Operand
PARM_THETA
The syntax of this key word common to the control
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 controls is described in the document [U4.50.01].
3.12 Word
key
CONVERGENCE
The syntax of this key word common to the control
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 the control
STAT_NON_LINE
is described in the document
[U4.51.03].
3.15 Operand
INFORMATION
The syntax of this key word common to the control
STAT_NON_LINE
is described in the document
[U4.51.03].
3.16 Operand
TITRATE
TITRATE = tx
tx
is the title of calculation. It will be printed at the head results. See [U4.03.01].
background image
Code_Aster
®
Version
7.4
Titrate:
Macro-control
CALC_PRECONT
Date:
01/02/05
Author (S):
S. MICHEL-PONNELLE, A. ASSIRE
Key
:
U4.42.05-A
Page
:
7/10
Instruction manual
U4.4- booklet: Modeling
HT-66/05/004/A
4 Example
One details here the main stages of a structural analysis containing of the cables of prestressing
(case-test SSNV164)
4.1
Before using the macro-control
·
Definition of the cables
In the event of setting in not-simultaneous voltage 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,
)
background image
Code_Aster
®
Version
7.4
Titrate:
Macro-control
CALC_PRECONT
Date:
01/02/05
Author (S):
S. MICHEL-PONNELLE, A. ASSIRE
Key
:
U4.42.05-A
Page
:
8/10
Instruction manual
U4.4- booklet: Modeling
HT-66/05/004/A
·
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 voltage 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 voltage of the cables (optional)
It is a question here of taking into account the loadings which apply to the structure before the setting in
voltage 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 WITHOUT, which imposes that the stress 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 (LOAD = CMCAB1),
_F (LOAD = CMCAB3),
_F (LOAD = CMCAB5),),
INCREMENT=_F (LIST_INST = LINST,
INST_FIN = 150.),)
background image
Code_Aster
®
Version
7.4
Titrate:
Macro-control
CALC_PRECONT
Date:
01/02/05
Author (S):
S. MICHEL-PONNELLE, A. ASSIRE
Key
:
U4.42.05-A
Page
:
9/10
Instruction manual
U4.4- booklet: Modeling
HT-66/05/004/A
4.2
Use of the macro-control: setting in successive voltage of
5 cables
One calls upon the macro-control 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 (LOAD = 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 (LOAD = CMCAB1,),
_F (LOAD = CMCAB3,),
CABLE_BP = (CAB_BP5,),
INCREMENT=_F (LIST_INST = LINST,
INST_FIN = 600.,
SUBD_PAS = 4,
SUBD_PAS_MINI = 0.01,),)
background image
Code_Aster
®
Version
7.4
Titrate:
Macro-control
CALC_PRECONT
Date:
01/02/05
Author (S):
S. MICHEL-PONNELLE, A. ASSIRE
Key
:
U4.42.05-A
Page
:
10/10
Instruction manual
U4.4- booklet: Modeling
HT-66/05/004/A
4.3
Continuation of the loading after setting in voltage 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 (LOAD = CMCAB1),
_F (LOAD = CMCAB3),
_F (LOAD = CMCAB5),
_F (CHARGE=CHMECA,
FONC_MULT = FCT,),
INCREMENT=_F (LIST_INST = LINST,),)