Moment curvature with uncoupled fiber elements

[giorgioirotina]

I need to calculate moment curvature of a section where concrete and steel behavior is uncoupled in two different elements.
It may look nonsense because I could just add concrete and steel fibers to one section, but is part of a more complicated model that I cannot explain here for reasons of time.
Let m know if you need more detail on the problem I want to solve.


This is the script


# Remove existing model
wipe

# Create ModelBuilder

model BasicBuilder -ndm 2 -ndf 3

########################################## MATERIALS ################################################
#
#
########################################## MATERIALS ################################################
# concrete
#uniaxialMaterial Concrete01 $matTag $fpc $epsc0 $fpcu $epsU
uniaxialMaterial Concrete01 1 -5.0 -0.002 0.0 -0.006

# steel
#uniaxialMaterial Steel01 $matTag $Fy $E0 $b
uniaxialMaterial Steel01 2 60 30000 0.01

########################################## SECTIONS ################################################
#
#
########################################## SECTIONS ################################################

# section Elastic $secTag $E $A $Iz
section Elastic 3 10000 100 100000
section Elastic 4 10000 100 .000001

# concrete section
section Fiber 1 {

#patch rect $matTag $numSubdivY $numSubdivZ $yI $zI $yJ $zJ
patch rect 1 10 2 -5 -2 5 2

}

section Fiber 2 {

# fiber $yLoc $zLoc $A $matTag
fiber 0 0 1 2

}


########################################## NODES ################################################
#
#
########################################## NODES ################################################


# Define concrete nodes
node 1 0 0
node 2 1 0

# Define steel nodes
node 3 0 -4
node 4 1 -4


########################################## RESTRAINS/CONSTRAINS ################################################
#
#
########################################## RESTRAINS/CONSTRAINS ################################################


# Fix degrees of point 1
fix 1 1 1 1

equalDOF 4 2 1 2 3
equalDOF 1 3 2 3
# geometric transformation
geomTransf Linear 1


########################################## ELEMENTS ################################################
#
#
########################################## ELEMENTS ################################################


# Define element

#element forceBeamColumn $eleTag $iNode $jNode $transfTag "IntegrationType arg1 arg2 ..."


element forceBeamColumn 1 1 2 1 HingeMidpoint 1 1 1 1 3
element forceBeamColumn 2 3 4 1 HingeMidpoint 2 1 2 1 3



########################################## RECORDER ################################################
#
#
########################################## RECORDER ################################################

# Create recorder
recorder Node -file M.out -node 2 -dof 3 disp
recorder Node -file R.out -node 1 -dof 3 reaction



########################################## ANALYSIS PARAMETERS ################################################
#
#
########################################## ANALYSIS PARAMETERS ################################################

# Define analysis parameters
integrator LoadControl 0.0
system SparseGeneral -piv;
test NormUnbalance 1.0e-6 400
numberer Plain
constraints Plain
algorithm Newton
analysis Static
analyze 1

########################################## LOAD CASE ################################################
#
#
########################################## LOAD CASE ################################################


# Define reference moment
pattern Plain 1 "Linear" {

load 2 0 0 1
}


# Use displacement control at node 2
#integrator DisplacementControl $node $dof $incr
integrator DisplacementControl 2 3 -0.001

# Do the analysis
analyze 10