Question on "Recorder element " for section analys
Moderators: silvia, selimgunay, Moderators
Question on "Recorder element " for section analys
We have a trouble with using
1. Objective of Analysis
Section analysis --> record "stress-strain"
2. We did run successfully,
Section analysis--> recorder node.
But, failed to run
Section analysis--> recorder element
3. Error massage
is in seperately poped message box.
4. The way we checked for error is
1. Move the place record command itself "mc_test.tcl" to"test.tcl" and to other lines
2. Check with all example command in Opensees Manual (different order of option )
3. Check with other coordinate (y,z) -->(Radius, angle)
4. Check with "fiber"--> 'layer'(steel)
5. Check with insert material tag (opensees manual page341) and if coordinate (y,z) exist
--> finally, all failed
I'm attaching the file we used. Please check and run the files.
The main file is "test.tcl"
The file which include the record command is "mc_test.tcl"
Thank you so much for your help.
And, I'm looking forward to your answer.
==========================================
#mc_test.tcl
# Arguments
# secTag -- tag identifying section to be analyzed
# axialLoad -- axial load applied to section (negative is compression)
# maxK -- maximum curvature reached during analysis
# numIncr -- number of increments used to reach maxK (default 100)
#
# Sets up a recorder which writes moment-curvature results to file
# section$secTag.out ... the moment is in column 1, and curvature in column 2
proc MomentCurvature {secTag axialLoad maxK numIncr } {
# Define two nodes at (0,0)
node 1 0.0 0.0
node 2 0.0 0.0
# Fix all degrees of freedom except axial and bending
fix 1 1 1 1
fix 2 0 1 0
# Define element
# tag ndI ndJ secTag
element zeroLengthSection 1 1 2 $secTag
# Create recorder
recorder Node -file section$secTag.out -time -node 2 -dof 3 disp
# recorder Element -file ele1sec1Force.out -ele 1 section 1fiber 0.0 0.0 1 stressStrain
# Define constant axial load
pattern Plain 1 "Constant" {
load 2 $axialLoad 0.0 0.0
}
# Define analysis parameters
integrator LoadControl 0.0
system SparseGeneral -piv; # Overkill, but may need the pivoting!
test NormUnbalance 1.0e-9 10
numberer Plain
constraints Plain
algorithm Newton
analysis Static
# Do one analysis for constant axial load
analyze 1
# Define reference moment
pattern Plain 2 "Linear" {
load 2 0.0 0.0 1.0
}
# Compute curvature increment
set dK [expr $maxK/$numIncr]
# Use displacement control at node 2 for section analysis
integrator DisplacementControl 2 3 $dK 1 $dK $dK
# Do the section analysis
analyze $numIncr
}
1. Objective of Analysis
Section analysis --> record "stress-strain"
2. We did run successfully,
Section analysis--> recorder node.
But, failed to run
Section analysis--> recorder element
3. Error massage
is in seperately poped message box.
4. The way we checked for error is
1. Move the place record command itself "mc_test.tcl" to"test.tcl" and to other lines
2. Check with all example command in Opensees Manual (different order of option )
3. Check with other coordinate (y,z) -->(Radius, angle)
4. Check with "fiber"--> 'layer'(steel)
5. Check with insert material tag (opensees manual page341) and if coordinate (y,z) exist
--> finally, all failed
I'm attaching the file we used. Please check and run the files.
The main file is "test.tcl"
The file which include the record command is "mc_test.tcl"
Thank you so much for your help.
And, I'm looking forward to your answer.
==========================================
#mc_test.tcl
# Arguments
# secTag -- tag identifying section to be analyzed
# axialLoad -- axial load applied to section (negative is compression)
# maxK -- maximum curvature reached during analysis
# numIncr -- number of increments used to reach maxK (default 100)
#
# Sets up a recorder which writes moment-curvature results to file
# section$secTag.out ... the moment is in column 1, and curvature in column 2
proc MomentCurvature {secTag axialLoad maxK numIncr } {
# Define two nodes at (0,0)
node 1 0.0 0.0
node 2 0.0 0.0
# Fix all degrees of freedom except axial and bending
fix 1 1 1 1
fix 2 0 1 0
# Define element
# tag ndI ndJ secTag
element zeroLengthSection 1 1 2 $secTag
# Create recorder
recorder Node -file section$secTag.out -time -node 2 -dof 3 disp
# recorder Element -file ele1sec1Force.out -ele 1 section 1fiber 0.0 0.0 1 stressStrain
# Define constant axial load
pattern Plain 1 "Constant" {
load 2 $axialLoad 0.0 0.0
}
# Define analysis parameters
integrator LoadControl 0.0
system SparseGeneral -piv; # Overkill, but may need the pivoting!
test NormUnbalance 1.0e-9 10
numberer Plain
constraints Plain
algorithm Newton
analysis Static
# Do one analysis for constant axial load
analyze 1
# Define reference moment
pattern Plain 2 "Linear" {
load 2 0.0 0.0 1.0
}
# Compute curvature increment
set dK [expr $maxK/$numIncr]
# Use displacement control at node 2 for section analysis
integrator DisplacementControl 2 3 $dK 1 $dK $dK
# Do the section analysis
analyze $numIncr
}
You should download and execute the latest version of OpenSees and will not have a crash problem, but do read below.
Last edited by silvia on Wed Jun 01, 2005 10:23 am, edited 1 time in total.
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
the zeroLength element only has 1 section .. the section identifier is therefore not needed
in the recorder command .. this is explained as a note in the manual .. also you either provide the fiber number or a coordinate, not both .. correct command is:
recorder Element -file ele1sec1Force.out -ele 1 section fiber 0.0 0.0 stressStrain
in the recorder command .. this is explained as a note in the manual .. also you either provide the fiber number or a coordinate, not both .. correct command is:
recorder Element -file ele1sec1Force.out -ele 1 section fiber 0.0 0.0 stressStrain
in general, though, you should always check your file, as you didn't have a space in the following line:
recorder Element -file ele1sec1Force.out -ele 1 section 1fiber 0.0 0.0 1 stressStrain
ALSO,
why not use the moment-curvature script I had posted a while ago???
recorder Element -file ele1sec1Force.out -ele 1 section 1fiber 0.0 0.0 1 stressStrain
ALSO,
why not use the moment-curvature script I had posted a while ago???
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104