Hi,
I have a 3D bridge model with concrete columns models with fiber section.
I use eleResponse command to save the curvature of the elements in tcl lists for further analysis using the following format:
eleResponse $element_number section $integration_point deformation
I used this command in the last version (2.4.4) and it worked. However, in the latest version (2.4.6) I get zeros for the curvature.
Is it possible that the new version has a bug?
I would appreciate it if you can help me with this.
Zero curvature recorded byeleResponse in version 2.4.6
Moderators: silvia, selimgunay, Moderators
-
- Posts: 24
- Joined: Mon Jun 18, 2012 9:29 am
- Location: Lehigh University
Re: Zero curvature recorded byeleResponse in version 2.4.6
it should not have changed .. can you post a 1 element example
-
- Posts: 24
- Joined: Mon Jun 18, 2012 9:29 am
- Location: Lehigh University
Re: Zero curvature recorded byeleResponse in version 2.4.6
Thank you for your prompt response.
I double checked and it seems that the number of digits stored for the output in version 2.4.6 is less than version 2.4.4
If I want the outputs in a text file the "Recorder" command works just fine. However, "eleResponse" not, since the values of curvature is very small and the numbers are not captured as the number of digits are not enough.
The following is a very simple example I just put together. This is a one element circular steel column with fiber section. A concentrated load is applied at the top and the curvature of the 5th integration point along the column is recorded.
I would appreciate it if you run it with both versions and compare the results and let me know if I am right or not, and what I should do in such cases.
wipe
model BasicBuilder -ndm 2 -ndf 3
set fy 0.414; # fy steel yield stress in KN/mm^2
set Es 200.0; # Es steel modulus of elasticity in KN/mm^2
set hr 0.018; # hr hardening ratio
# uniaxialMaterial Elastic 1 $Es
# uniaxialMaterial Steel01 $matTag $Fy $E0 $b <$a1 $a2 $a3 $a4>
uniaxialMaterial Steel01 1 $fy $Es $hr
section Fiber 1 {
# patch circ $matTag $numSubdivCirc $numSubdivRad $yCenter $zCenter $intRad $extRad <$startAng $endAng>
patch circ 1 12 8 0.0 0.0 0.0 50.0 0.0 360.0
}
node 1 0.0 0.0;
node 2 0.0 2000.0;
fix 1 1 1 1;
fix 2 0 0 0;
geomTransf Linear 1;
element dispBeamColumn 1 1 2 10 1 1;
recorder Element -file def.out -ele 1 section 5 deformation
# recorder Node -file nodesD.out -node 2 -dof 1 disp;
pattern Plain 1 Linear {
load 2 5.0 0.0 0.0
}
constraints Transformation
numberer RCM
system BandGeneral
test NormDispIncr 1.0e-10 1000
algorithm Newton
integrator LoadControl 0.01
# integrator DisplacementControl 2 1 0.01
analysis Static
# analyze 100
set def {};
set node {};
set analyze_count 1;
set n_analyze 100;
while {$analyze_count<=$n_analyze} {
set ok [analyze 1];
set defout [eleResponse 1 section 5 deformation];
puts "[lindex $defout 1]"
incr analyze_count;
}
wipe
I double checked and it seems that the number of digits stored for the output in version 2.4.6 is less than version 2.4.4
If I want the outputs in a text file the "Recorder" command works just fine. However, "eleResponse" not, since the values of curvature is very small and the numbers are not captured as the number of digits are not enough.
The following is a very simple example I just put together. This is a one element circular steel column with fiber section. A concentrated load is applied at the top and the curvature of the 5th integration point along the column is recorded.
I would appreciate it if you run it with both versions and compare the results and let me know if I am right or not, and what I should do in such cases.
wipe
model BasicBuilder -ndm 2 -ndf 3
set fy 0.414; # fy steel yield stress in KN/mm^2
set Es 200.0; # Es steel modulus of elasticity in KN/mm^2
set hr 0.018; # hr hardening ratio
# uniaxialMaterial Elastic 1 $Es
# uniaxialMaterial Steel01 $matTag $Fy $E0 $b <$a1 $a2 $a3 $a4>
uniaxialMaterial Steel01 1 $fy $Es $hr
section Fiber 1 {
# patch circ $matTag $numSubdivCirc $numSubdivRad $yCenter $zCenter $intRad $extRad <$startAng $endAng>
patch circ 1 12 8 0.0 0.0 0.0 50.0 0.0 360.0
}
node 1 0.0 0.0;
node 2 0.0 2000.0;
fix 1 1 1 1;
fix 2 0 0 0;
geomTransf Linear 1;
element dispBeamColumn 1 1 2 10 1 1;
recorder Element -file def.out -ele 1 section 5 deformation
# recorder Node -file nodesD.out -node 2 -dof 1 disp;
pattern Plain 1 Linear {
load 2 5.0 0.0 0.0
}
constraints Transformation
numberer RCM
system BandGeneral
test NormDispIncr 1.0e-10 1000
algorithm Newton
integrator LoadControl 0.01
# integrator DisplacementControl 2 1 0.01
analysis Static
# analyze 100
set def {};
set node {};
set analyze_count 1;
set n_analyze 100;
while {$analyze_count<=$n_analyze} {
set ok [analyze 1];
set defout [eleResponse 1 section 5 deformation];
puts "[lindex $defout 1]"
incr analyze_count;
}
wipe