Question on "Recorder element " for section analys

[dechoe]

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
}

[silvia]

please post the main program that would call this proc.

[silvia]

You should download and execute the latest version of OpenSees and will not have a crash problem, but do read below.

[fmk]

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

[silvia]

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???