Push over of 3-story R.C. shear wall

[abdallahmontaser]

Hello,
I'm a beginner in OpenSees and I'm trying to conducted a Push over on a 3-story R.C. shear wall using FBE. Each story has a mass of 12t affecting on lateral direction only. When I run the model, it closes so fast after applying the vertical load. Can any one revise my commands and tell me about my mistakes! units in mm and N
Any help will be much appreciated!!

model basic -ndm 2 -ndf 3
# Nodes Formation and Restrains
node 1 0 0
node 2 0 1360
node 3 0 2720
node 4 0 4080
fix 1 1 1 1
puts "Nodes Formation and Restrains Done"

#Material Definition

# Concrete
set fc -36.5; set ec -0.0021; set Ec 32700;
set ft 1.87; set et 0.00008; set xp 23;
set xn 2; set rc 10;
uniaxialMaterial Concrete07 1 $fc $ec $Ec $ft $et $xp $xn $rc

#Steel
set Ev 215700; set R0 12
set cR1 0.925; set cR2 0.15
set Fyp3 474.2; set Fyn3 474.2;
set bp3 0.01; set bn3 0.01
uniaxialMaterial SteelMPF 12 $Fyp3 $Fyn3 $Ev $bp3 $bn3 $R0 $cR1 $cR2
uniaxialMaterial MinMax 2 12 -min -0.062 -max 0.062
puts "Material Definition Done"

#Section Definition
geomTransf PDelta 1
section Fiber 1 {
patch quad 1 90 10 -450 -50 450 -50 450 50 -450 50

layer straight 2 10 17 -433 -33 433 -33

layer straight 2 10 17 -433 33 433 33
}
puts "Section Definition Done"

#MASS
mass 2 12 0 0
mass 3 12 0 0
mass 4 12 0 0

set integration1 "Lobatto 1 5"
element forceBeamColumn 1 1 2 1 $integration1
element forceBeamColumn 2 2 3 1 $integration1
element forceBeamColumn 3 3 4 1 $integration1

# Apply Axial Load

set tolAx 1.0e-7

set iterAx 10

pattern Plain 1 Linear { load 4 0.0 -100000 0.0 }

integrator LoadControl 0.1;

system SparseGeneral -piv; test NormUnbalance $tolAx $iterAx 1;

numberer RCM; constraints Transformation;

algorithm KrylovNewton; analysis Static;

puts "Gravity Load Application Done"

# perform the gravity load analysis,

analyze 10

# Set the gravity loads to be constant & reset the time in the domain

loadConst -time 0.0

puts "Gravity load analysis completed"

# Set Lateral Loads

recorder Node -file TopDisp4.out -time -node 4 -dof 1 disp
recorder Node -file TopDisp3.out -time -node 3 -dof 1 disp
recorder Node -file TopDisp2.out -time -node 2 -dof 1 disp
recorder Node -file eigen.out -node 1 2 3 4 -dof 1 "eigen 0 1 2 3 "
recorder Node -file Vbase.out -node 1 -dof 1 reaction


set tolLat 1.0e-12; set iterLat 100

recorder Node -file TopDisp2.out -time -node 4 -dof 1 disp

pattern Plain 2 Linear { load 2 111.1 0.0 0.0
load 3 222.2 0.0 0.0
load 4 333.3 0.0 0.0} # triangular load distribution factors= w*hi/sigma w*htotal= 1/9 , 2/9 and 1/3 Then they are multiplied by 1000 to get base shear in KN. Is that right?

test NormDispIncr $tolLat $iterLat 1

puts "Lateral Load Application Done"

#Analyze lateral Loads
#U=160 mm
integrator DisplacementControl 4 1 0.1
analyze 1600

[selimgunay]

The quickest way of debugging your tcl file is to add puts commands in several lines to see until which point it proceeds fine. Also it is needed to run opensees from the command prompt to avoid sudden closing of the window

[abdallahmontaser]

Thank you so much for your reply !!
I have added (puts) command and there is error in reading the load pattern information of lateral pushover. I think I'm defining the pushover in wrong manner. Could you please check it and tell me what's wrong in it and how can I fix it?
Kind Regards

[ahmadbsr]

The quickest way of debugging your tcl file is to add puts commands in several lines to see until which point it proceeds fine. Also it is needed to run opensees from the command prompt to avoid sudden closing of the window

Can I model the Shear wall in the bay of frame with SFI element ?