Zero length section element

[jiachenghan]

Hello everyone
I am trying to simulate the bond slip of reinforced concrete using ZeroLengthSection. But the hysteresis curve of the model has not changed. It seems that the ZeroLengthSection is not working.I don't know why. I would be grateful if you could help me. This is my script:
wipe;
model BasicBuilder -ndm 3 -ndf 6;
set dataDir Data_Hyster;
file mkdir $dataDir;
node 1 0. 0. 0.;
node 2 0. 0. 1600.;
node 3 0. 0. 0.;
set fpc -29.05;
set epsc0 -0.002;
set fpcu -23.24;
set epsU -0.004;
set u 0.0086;
set Fubar 315;
set f [expr $u*0.75*$Fubar]
set fp2c [expr $fpc*(-1.254+2.254*sqrt(1+7.94*$f/(-$fpc))-2*$f/(-$fpc))];
set epsc20 [expr 0.002*(1+5*($f/(-$fpc)-1))];
set eps2U [expr 0.004+0.9*$u*$Fubar/300];
set fpc2u [expr 2*$fp2c*($eps2U/$epsc20)/(1+($eps2U/$epsc20)*($eps2U/$epsc20))];
set R0 20.;
set cR1 0.925;
set cR2 0.15;
set IDconcCore1 1;
set IDconcCore2 2;
uniaxialMaterial Concrete01 $IDconcCore1 $fpc $epsc0 $fpcu $epsU;
uniaxialMaterial Concrete01 $IDconcCore2 $fp2c $epsc20 $fpc2u $eps2U;
set Fybar 407;
set Esbar 200000;
set Bsbar 0.01;
set IDSteelBar1 3;
uniaxialMaterial Steel02 $IDSteelBar1 $Fybar $Esbar $Bsbar $R0 $cR1 $cR2;
uniaxialMaterial Bond_SP01 6 407 0.49 573.25 17.14 0.4 0.35;
set bardimeth 20.;
set barArea [expr 3.1415926*$bardimeth*$bardimeth/4.];
set WidthColum 400;
set depthColum 300;
set coverdepth 20.;
set coverY [expr $WidthColum/2.0];
set coverZ [expr $depthColum/2.0];
set coreY [expr $coverY-$coverdepth];
set coreZ [expr $coverZ-$coverdepth];
set normSec 1;
set zeroLengthSection 3;
set coreID 2;
set coverID 1;
set coreIDy 5;
section fiberSec $normSec {
patch quadr $coreID 36 26 -$coreY -$coreZ $coreY -$coreZ $coreY $coreZ -$coreY -$coreZ;
patch quadr $coverID 40 1 -$coverY -$coverZ $coverY -$coverZ $coreY -$coreZ -$coreY -$coreZ;
patch quadr $coverID 1 30 $coreY -$coreZ $coverY -$coverZ $coverY $coverZ $coreY $coreZ;
patch quadr $coverID 40 1 -$coreY $coreZ $coreY $coreZ $coverY $coverZ -$coverY $coverZ;
patch quadr $coverID 1 30 -$coverY -$coverZ -$coreY -$coreZ -$coreY $coreZ -$coverY $coverZ;
layer straight 3 3 $barArea -$coreY -$coreZ $coreY -$coreZ;
layer straight 3 3 $barArea -$coreY $coreZ $coreY $coreZ;
layer straight 3 2 $barArea -$coreY 0 $coreY 0;
};
section fiberSec $zeroLengthSection {
patch quadr $coreIDy 36 26 -$coreY -$coreZ $coreY -$coreZ $coreY $coreZ -$coreY -$coreZ;
patch quadr $coverID 40 1 -$coverY -$coverZ $coverY -$coverZ $coreY -$coreZ -$coreY -$coreZ;
patch quadr $coverID 1 30 $coreY -$coreZ $coverY -$coverZ $coverY $coverZ $coreY $coreZ;
patch quadr $coverID 40 1 -$coreY $coreZ $coreY $coreZ $coverY $coverZ -$coverY $coverZ;
patch quadr $coverID 1 30 -$coverY -$coverZ -$coreY -$coreZ -$coreY $coreZ -$coverY $coverZ;
layer straight 6 3 $barArea -$coreY -$coreZ $coreY -$coreZ;
layer straight 6 3 $barArea -$coreY $coreZ $coreY $coreZ;
layer straight 6 2 $barArea -$coreY 0 $coreY 0;
};
set MatTorsion 99;
set Ubig 1.189E8;
set MatTorsion1 97;
set Ubig1 1.32E9;
set MatTorsion2 98;
set torsnormSec 101;
set torsstrengthenSec 102;
set torszeroLengthSection 103;
uniaxialMaterial Elastic $MatTorsion $Ubig;
uniaxialMaterial Elastic $MatTorsion1 $Ubig1;
uniaxialMaterial Elastic $MatTorsion2 $Ubig1;
section Aggregator $torsnormSec $MatTorsion1 Vy $MatTorsion2 Vz $MatTorsion T -section $normSec;
section Aggregator $torszeroLengthSection $MatTorsion1 Vy $MatTorsion2 Vz $MatTorsion T -section $zeroLengthSection;
set IDColTransf 1;
set ColTransfType PDelta ;
geomTransf $ColTransfType $IDColTransf 0 1 0;
element dispBeamColumn 1 1 2 5 $normSec $IDColTransf;
element zeroLengthSection 2 3 1 $torszeroLengthSection -orient 0 0 1 1 0 0;
recorder Node -file $dataDir/Displace.txt -time -node 2 -dof 1 2 3 disp;
fix 1 1 1 1 1 1 1;
fix 2 0 0 0 1 0 1;
fix 3 0 1 0 1 0 1;
mass 2 1. 1. 1. 0. 0. 0.;
pattern Plain 1 Linear {
load 2 0. 0. -558000. 0. 0. 0.;
}
set Tol 1.0e-8;
constraints Plain;
numberer RCM;
system BandGeneral;
test EnergyIncr $Tol 60;
algorithm Newton;
set NstepGravity 10;
set DGravity [expr 1./$NstepGravity];
integrator LoadControl $DGravity;
analysis Static;
analyze $NstepGravity;
loadConst -time 0.0
puts "Model Built\n"
constraints Plain;
numberer RCM;
system BandGeneral;
variable TolStatic 1E-5;
variable maxNumIterStatic 30;
variable printFlagStatic 1;
test NormDispIncr $TolStatic $maxNumIterStatic $printFlagStatic;
set algorithmType Newton;
algorithm $algorithmType;
analysis Static;
set IDctrlNode 2;
pattern Plain 100 Linear {
load $IDctrlNode 1. 0. 0. 0. 0. 0.;
}
integrator DisplacementControl $IDctrlNode 1 0.25;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 -0.5;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 0.75;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 -1;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 1.5;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 -2;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 2.5;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 -3;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 3.5;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 -4;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 4.5;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 -5;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 5.5;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 -6;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 6.5;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 -7;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 7.5;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 -8;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 8.5;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 -9;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 9.5;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 -10;
analyze 20;
integrator DisplacementControl $IDctrlNode 1 5;
analyze 20;

[luochuizi]

hi!i am a new boy of opensees and i fond your topic by accident.i want to help you so i checked your model.and i fond a mistake in your model.but i don't know if it will solve your problem.anyway i will tell you the issue as follow:
the define of the zero length section commends:
section fiberSec $zeroLengthSection {
patch quadr $coreIDy 36 26 -$coreY -$coreZ $coreY -$coreZ $coreY $coreZ -$coreY -$coreZ;
in the model,i fond the $coreIDy is 5,but i didn't fond the material 5 that defined the core concrete.i think you should check it anyway.
hope it will help you.

[Eric]

the fibersec should between node 2 and node 3 ,not node 2 and 1