Fiber section element recorder
Moderators: silvia, selimgunay, Moderators
Fiber section element recorder
I've got a problem with my element recorders for fiber section elements, output file only contains timesteps and nothing else! Any hints and suggestions will be appreciated. here is one of the Recorders:
recorder Element -file ele348Force.out -time -ele 348 section 6 fiber 0.33 0.33 22 force
and here is the fiber section:
set HSec 1.2192; # Column Depth
set BSec 1.0668; # Column Width
set coverH 0.0762; # Column cover to reinforcing steel NA, parallel to H
set coverB 0.0762; # Column cover to reinforcing steel NA, parallel to B
set numBarsTop 6; # number of longitudinal-reinforcement bars in steel layer. -- top
set numBarsBot 9; # number of longitudinal-reinforcement bars in steel layer. -- bot
set numBarsIntTot 4; # number of longitudinal-reinforcement bars in steel layer. -- total intermediate skin reinforcement, symm about y-axis
set barAreaTop [expr 0.0286*0.0286*3.1415*0.25]; # area of longitudinal-reinforcement bars -- top
set barAreaBot [expr 0.0286*0.0286*3.1415*0.25]; # area of longitudinal-reinforcement bars -- bot
set barAreaInt [expr 0.0159*0.0159*3.1415*0.25]; # area of longitudinal-reinforcement bars -- intermediate skin reinf
set CapBeamSecTag 4; # set tag for symmetric section
set coverY [expr $HSec/2.0]; # The distance from the section z-axis to the edge of the cover concrete -- outer edge of cover concrete
set coverZ [expr $BSec/2.0]; # The distance from the section y-axis to the edge of the cover concrete -- outer edge of cover concrete
set coreY [expr $coverY-$coverH]; # The distance from the section z-axis to the edge of the core concrete -- edge of the core concrete/inner edge of cover concrete
set coreZ [expr $coverZ-$coverB]; # The distance from the section y-axis to the edge of the core concrete -- edge of the core concrete/inner edge of cover concreteset nfY 16; # number of fibers for concrete in y-direction
set nfY 14; # number of fibers for concrete in y-direction
set nfZ 14; # number of fibers for concrete in z-direction
set numBarsInt [expr $numBarsIntTot/2]; # number of intermediate bars per side
section Fiber $CapBeamSecTag {; # Define the fiber section
patch quadr 22 $nfZ $nfY -$coreY $coreZ -$coreY -$coreZ $coreY -$coreZ $coreY $coreZ; # Define the core patch
patch quadr 21 1 $nfY -$coverY $coverZ -$coreY $coreZ $coreY $coreZ $coverY $coverZ; # Define the four cover patches
patch quadr 21 1 $nfY -$coreY -$coreZ -$coverY -$coverZ $coverY -$coverZ $coreY -$coreZ
patch quadr 21 $nfZ 1 -$coverY $coverZ -$coverY -$coverZ -$coreY -$coreZ -$coreY $coreZ
patch quadr 21 $nfZ 1 $coreY $coreZ $coreY -$coreZ $coverY -$coverZ $coverY $coverZ
layer straight 23 $numBarsInt $barAreaInt -$coreY $coreZ $coreY $coreZ; # intermediate skin reinf. +z
layer straight 23 $numBarsInt $barAreaInt -$coreY -$coreZ $coreY -$coreZ; # intermediate skin reinf. -z
layer straight 23 $numBarsTop $barAreaTop $coreY $coreZ $coreY -$coreZ; # top layer reinfocement
layer straight 23 $numBarsBot $barAreaBot -$coreY $coreZ -$coreY -$coreZ; # bottom layer reinforcement
}; # end of fibersection definition
######### assign torsional Stiffness for 3D Model
set CapBeamSecTagTorsion 30; # ID tag for torsional section behavior
set CapBeamSecTag3D 6; # ID tag for combined behavior for 3D model
uniaxialMaterial Elastic $CapBeamSecTagTorsion 96775483.; # define elastic torsional stiffness
section Aggregator $CapBeamSecTag3D $CapBeamSecTagTorsion T -section $CapBeamSecTag; # combine section properties
recorder Element -file ele348Force.out -time -ele 348 section 6 fiber 0.33 0.33 22 force
and here is the fiber section:
set HSec 1.2192; # Column Depth
set BSec 1.0668; # Column Width
set coverH 0.0762; # Column cover to reinforcing steel NA, parallel to H
set coverB 0.0762; # Column cover to reinforcing steel NA, parallel to B
set numBarsTop 6; # number of longitudinal-reinforcement bars in steel layer. -- top
set numBarsBot 9; # number of longitudinal-reinforcement bars in steel layer. -- bot
set numBarsIntTot 4; # number of longitudinal-reinforcement bars in steel layer. -- total intermediate skin reinforcement, symm about y-axis
set barAreaTop [expr 0.0286*0.0286*3.1415*0.25]; # area of longitudinal-reinforcement bars -- top
set barAreaBot [expr 0.0286*0.0286*3.1415*0.25]; # area of longitudinal-reinforcement bars -- bot
set barAreaInt [expr 0.0159*0.0159*3.1415*0.25]; # area of longitudinal-reinforcement bars -- intermediate skin reinf
set CapBeamSecTag 4; # set tag for symmetric section
set coverY [expr $HSec/2.0]; # The distance from the section z-axis to the edge of the cover concrete -- outer edge of cover concrete
set coverZ [expr $BSec/2.0]; # The distance from the section y-axis to the edge of the cover concrete -- outer edge of cover concrete
set coreY [expr $coverY-$coverH]; # The distance from the section z-axis to the edge of the core concrete -- edge of the core concrete/inner edge of cover concrete
set coreZ [expr $coverZ-$coverB]; # The distance from the section y-axis to the edge of the core concrete -- edge of the core concrete/inner edge of cover concreteset nfY 16; # number of fibers for concrete in y-direction
set nfY 14; # number of fibers for concrete in y-direction
set nfZ 14; # number of fibers for concrete in z-direction
set numBarsInt [expr $numBarsIntTot/2]; # number of intermediate bars per side
section Fiber $CapBeamSecTag {; # Define the fiber section
patch quadr 22 $nfZ $nfY -$coreY $coreZ -$coreY -$coreZ $coreY -$coreZ $coreY $coreZ; # Define the core patch
patch quadr 21 1 $nfY -$coverY $coverZ -$coreY $coreZ $coreY $coreZ $coverY $coverZ; # Define the four cover patches
patch quadr 21 1 $nfY -$coreY -$coreZ -$coverY -$coverZ $coverY -$coverZ $coreY -$coreZ
patch quadr 21 $nfZ 1 -$coverY $coverZ -$coverY -$coverZ -$coreY -$coreZ -$coreY $coreZ
patch quadr 21 $nfZ 1 $coreY $coreZ $coreY -$coreZ $coverY -$coverZ $coverY $coverZ
layer straight 23 $numBarsInt $barAreaInt -$coreY $coreZ $coreY $coreZ; # intermediate skin reinf. +z
layer straight 23 $numBarsInt $barAreaInt -$coreY -$coreZ $coreY -$coreZ; # intermediate skin reinf. -z
layer straight 23 $numBarsTop $barAreaTop $coreY $coreZ $coreY -$coreZ; # top layer reinfocement
layer straight 23 $numBarsBot $barAreaBot -$coreY $coreZ -$coreY -$coreZ; # bottom layer reinforcement
}; # end of fibersection definition
######### assign torsional Stiffness for 3D Model
set CapBeamSecTagTorsion 30; # ID tag for torsional section behavior
set CapBeamSecTag3D 6; # ID tag for combined behavior for 3D model
uniaxialMaterial Elastic $CapBeamSecTagTorsion 96775483.; # define elastic torsional stiffness
section Aggregator $CapBeamSecTag3D $CapBeamSecTagTorsion T -section $CapBeamSecTag; # combine section properties
Re: Fiber section element recorder
Change force to stress in your recorder
Re: Fiber section element recorder
Thanks for your quick reply. I tried all below combinations , none of them worked:
recorder Element -file ele348Force.out -time -ele 348 section 6 fiber 0.33 0.33 22 stress
recorder Element -file ele348Force.out -time -ele 348 section 6 fiber 0.33 0.33 22 Stress
recorder Element -file ele348Force.out -time -ele 348 section 6 fiber 0.33 0.33 22 stressStrain
recorder Element -file ele348Force.out -time -ele 348 section 6 fiber 0.33 0.33 22 Strain
recorder Element -file ele348Force.out -time -ele 348 section 6 fiber 0.33 0.33 22 strain
recorder Element -file ele348Force.out -time -ele 348 section 6 fiber 0.33 0.33 22 deformation
but this one gives forces in output file, however I don't for what fibers are they:
recorder Element -file ele348Force.out -time -ele 348 force
and the output is like:
0.01 -2092.81 -87058.7 -127010 -72478.5 1448.69 -23598.2 2092.81 141286 127010 -121614 2109.09 23598.2
0.02 -2092.1 -87011.3 -126959 -72445.7 1448.19 -23592.4 2092.1 141238 126959 -121567 2108.37 23592.4
0.03 -2091.81 -86941.6 -126885 -72398.6 1447.96 -23585.5 2091.81 141169 126885 -121495 2108.11 23585.5
0.04 -2091.64 -86790.5 -126731 -72296.8 1447.77 -23579.5 2091.64 141018 126731 -121340 2108.01 23579.5
0.05 -2090.8 -86586.5 -126533 -72161.9 1447.05 -23570.5 2090.8 140814 126533 -121128 2107.31 23570.5
0.06 -2089 -86254.1 -126225 -71943.3 1445.55 -23554.1 2089 140481 126225 -120782 2105.75 23554.1
0.07 -2087.11 -85761.9 -125782 -71621.7 1443.81 -23531.9 2087.11 139989 125782 -120266 2104.27 23531.9
recorder Element -file ele348Force.out -time -ele 348 section 6 fiber 0.33 0.33 22 stress
recorder Element -file ele348Force.out -time -ele 348 section 6 fiber 0.33 0.33 22 Stress
recorder Element -file ele348Force.out -time -ele 348 section 6 fiber 0.33 0.33 22 stressStrain
recorder Element -file ele348Force.out -time -ele 348 section 6 fiber 0.33 0.33 22 Strain
recorder Element -file ele348Force.out -time -ele 348 section 6 fiber 0.33 0.33 22 strain
recorder Element -file ele348Force.out -time -ele 348 section 6 fiber 0.33 0.33 22 deformation
but this one gives forces in output file, however I don't for what fibers are they:
recorder Element -file ele348Force.out -time -ele 348 force
and the output is like:
0.01 -2092.81 -87058.7 -127010 -72478.5 1448.69 -23598.2 2092.81 141286 127010 -121614 2109.09 23598.2
0.02 -2092.1 -87011.3 -126959 -72445.7 1448.19 -23592.4 2092.1 141238 126959 -121567 2108.37 23592.4
0.03 -2091.81 -86941.6 -126885 -72398.6 1447.96 -23585.5 2091.81 141169 126885 -121495 2108.11 23585.5
0.04 -2091.64 -86790.5 -126731 -72296.8 1447.77 -23579.5 2091.64 141018 126731 -121340 2108.01 23579.5
0.05 -2090.8 -86586.5 -126533 -72161.9 1447.05 -23570.5 2090.8 140814 126533 -121128 2107.31 23570.5
0.06 -2089 -86254.1 -126225 -71943.3 1445.55 -23554.1 2089 140481 126225 -120782 2105.75 23554.1
0.07 -2087.11 -85761.9 -125782 -71621.7 1443.81 -23531.9 2087.11 139989 125782 -120266 2104.27 23531.9
Re: Fiber section element recorder
Do you intend section 6 to refer to the section tag or to the integration point number (1,...,Np)?
Re: Fiber section element recorder
it is the section number.
Re: Fiber section element recorder
Thanks for your tip! it should be the integration point number not the section number.
Re: Fiber section element recorder
just one more question, for circular sections, should I mention Y-Z location of the fiber or R-Theta??