frenelli

[elisapaso]

Hi dear all,
I make a strengthening masonry arches with composite materials and now i shape the "frenelli" and the program gives me an error! I don't understand why!
This is the file tcl:

# ARCO VC_SRG

wipe;

# UNITA' DI MISURA: N, mm

model BasicBuilder -ndm 2 -ndf 3;

source DatiArco.tcl ;
source NodiArco.tcl ;

set PP 24.e-6;
set PI [expr 2*asin(1.0)] ; # pi.greco
set ss [expr .5*$s] ;
set dd [expr .5*$d] ;
set NNA [expr $NEA+1] ;
set PPA 24.e-6 ;
set P1E [expr $PPA*$L1E*$s*$d] ; # [N] Peso di un singolo elemento dell'arco
set PA [expr $P1E*$NEA] ; # [N] Peso totale dell'arco
set E0 1193. ; # Modulo elastico per il materiale PPGap [MPa]
set NC [expr int($NEA/2 + 1)] ; # Nodo di chiave
set sec 1 ;
set fmt 0.173 ;
set epsc0 -[expr 2*$fc/1193.] ; # Deformazione corrispondente al picco di resistenza per i materiali Concrete

#VINCOLI
fix 1 1 1 1;
fix $NNA 1 1 1 ;

#MATERIALI
uniaxialMaterial Concrete01 1 -$fc $epsc0 [expr $fc*.2] [expr 4*$epsc0] ;
uniaxialMaterial ElasticPPGap 2 $Ef $Fy 0. ;

#SEZIONE DI FIBRE
section Fiber 1 {patch quad 1 $NF_IJ $NF_JK -$ss -$dd $ss -$dd $ss $dd -$ss $dd ;
if {$SRG_extr==1} {
fiber $ss 0. $Af 2 }
if {$SRG_intr==1} {
fiber -$ss 0. $Af 2 } };


if {$ancor==0} {
section Fiber 2 {
patch quad $mat $NF_IJ $NF_JK -$ss -$dd $ss -$dd $ss $dd -$ss $dd
if {$SRG_intr==1} {
fiber -$ss 0. $Af 1 } } };

if {$ancor==1} {
for {set j 1} {$j<$NNA} {incr j 1} {
geomTransf Linear $j;
element nonlinearBeamColumn $j $j [expr $j+1] 2 1 $j }
} else {
for {set j 2} {$j<$NEA} {incr j 1} {
geomTransf Linear $j;
element nonlinearBeamColumn $j $j [expr $j+1] 2 2 $j };
geomTransf Linear 100 ;
element nonlinearBeamColumn 1 1 2 2 2 100 ;
geomTransf Linear 101 ;
element nonlinearBeamColumn $NEA $NEA $NNA 2 2 101 };

#FRENELLI
set HR 330; # [mm] Altezza del frenello

# Trova l'ultimo nodo del frenello
for {set n 1} {$n < [expr $NNA/2.]} {incr n 1} {

if {$Y($n) < $HR} {

set nR $n; # Ultimo nodo dei frenelli

} };


# Definisco i nodi dei frenelli

# Frenello di sinistra

for {set i 1001} {$i <= [expr 1000+$nR]} {incr i 1} {

set X(i) [expr -300.] ;

set Y(i) [expr $Y([expr $i-1000])] ;

node $i $X(i) $Y(i) 0. ;

fix $i 1 1 1 ;

};

# Frenello di destra

for {set i 2001} {$i <= [expr 2000+$nR]} {incr i 1} {

set X(i) [expr $X($NNA) + 300.] ;

set Y(i) [expr $Y([expr $i-(2000+$nR-$NNA)])] ;

node $i $X(i) $Y(i) 0. ;

fix $i 1 1 1 ;

};

# Definisco i frenelli

# Trasformazione Geometrica

geomTransf Corotational 900;

# Sezione dei frenelli

# Sezione di muratura mai rinforzata

# Frenello di sinistra

for {set i 2} {$i <= $nR} {incr i 1} {

element truss $i [expr $NNA+1-$i] [expr $NNA+$i] 2 ;

} };

puts ciao;
# Frenello di destra

for {set j 1001} {$j < [expr 1000+$nR]} {incr j 1} {

element truss $j $j [expr $j-1000] 2 ;

};

for {set j 2001} {$j < [expr 2000+$nR]} {incr j 1} {
element truss $j $j [expr $j-1000] 2 ;

};

[vitocasam]

Hi Elisa
I am trying on to do the same modellation. I will inform you of results.