can,t run

Forum for OpenSees users to post questions, comments, etc. on the use of the OpenSees interpreter, OpenSees.exe

Moderators: silvia, selimgunay, Moderators

Post Reply
shafeli
Posts: 6
Joined: Mon Sep 07, 2020 10:47 pm

can,t run

Post by shafeli »

hi
I would be happy if anyone can help me in this case.
I finish modelling my structure however, tried to run the model and get the outcome but there is no sign of it and i cannot get any output...
what am i have done wrong ??

my model is relating to shear walls using line element (displacement element) and each wall includes of 6 DBE,s along the wall
here I send you one story model to ease examination process


model basic -ndm 2 -ndf 3
###nodes for beams
node 101 0.0 180
node 102 144 180
node 103 161 180
node 104 162.7 180
node 105 171.35 180
node 106 180 180
node 107 188.65 180
node 108 197.29 180
node 109 198.99 180
node 110 215.99 180
node 111 360 180
node 1011 0.0 180
node 1021 144 180
node 1031 161 180
node 1041 162.7 180
node 1051 171.35 180
node 1061 180 180
node 1071 188.65 180
node 1081 197.29 180
node 1091 198.99 180
node 1101 215.99 180
node 1111 360 180


###EqualDof for Beams
equalDOF 101 1011 1 3 5
equalDOF 102 1021 1 3 5
equalDOF 103 1031 1 3 5
equalDOF 104 1041 1 3 5
equalDOF 105 1051 1 3 5
equalDOF 106 1061 1 3 5
equalDOF 107 1071 1 3 5
equalDOF 108 1081 1 3 5
equalDOF 109 1091 1 3 5
equalDOF 110 1101 1 3 5
equalDOF 111 1111 1 3 5

node 1 0 0
node 10 0 30
node 11 0 60
node 12 0 90
node 13 0 120
node 14 0 150
node 2 360 0.0
node 15 360 30
node 16 360 60
node 17 360 90
node 18 360 120
node 19 360 150




# Fix supports at base of columns
# tag DX DY RZ
fix 1 1 1 1
fix 2 1 1 1

# Define materials for nonlinear columns# ------------------------------------------

# CONCRETE matTag fpc epsc0 fpcu epsU lambda ft Ets
# Core concrete (confined)
uniaxialMaterial Concrete02 1 -7.0 -0.01046 -1.4 -0.0263 0.1 0.01206 220

# Cover concrete (unconfined)
uniaxialMaterial Concrete02 2 -5.0 -0.00247 -0.05 -0.0200 0.1 0.01019 185

# STEEL
# Reinforcing steel

# matTag Fy E b R0 cR1 cR2 <$a1 $a2 $a3 $a4 $sigInit>
uniaxialMaterial Steel02 3 60 29000 0.006 20 0.925 0.15
#for bar bucking
uniaxialMaterial MinMax 4 3
set wallLength 288
set wallThickness 24

set cover 1.5
#area of no. 8 no. 11 bars area of no. 9 bars area of no. 4 bars

set Asi 0.79
set Ase 1.56
set Asn 1
set Asf 0.2
# some variables derived from the parameters
set y1 [expr $wallLength/2.0]
set z1 [expr $wallThickness/2.0]




section Fiber 1 {
patch rect 2 10 1 -102 -$z1 102 $z1
patch rect 1 10 1 -$y1 -$z1 [expr -$y1+54] $z1
patch rect 1 10 1 [expr $y1-30] -$z1 $y1 $z1
# Create the reinforcing fibers (left, middle, right) first one is right

#right column longitudinal bar
###layer straight $matTag $numFiber $areaFiber $yStart $zStart $yEnd $zEnd
layer straight 3 4 $Asi [expr $y1-$cover] [expr $z1-$cover] [expr $y1-$cover] [expr $cover-$z1]
layer straight 3 2 $Asi [expr $y1-$cover-9] [expr $z1-$cover] [expr $y1-$cover-9] [expr $cover-$z1]
layer straight 3 2 $Asi [expr $y1-$cover-18] [expr $z1-$cover] [expr $y1-$cover-18] [expr $cover-$z1]
layer straight 3 4 $Asi [expr $y1-$cover-27] [expr $z1-$cover] [expr $y1-$cover-27] [expr $cover-$z1]

#right wall longitudinal bar

layer straight 3 4 $Asf [expr $y1-$cover-34.5] [expr $z1-$cover] [expr $y1-$cover-34.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $y1-$cover-40.5] [expr $z1-$cover] [expr $y1-$cover-40.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $y1-$cover-46.5] [expr $z1-$cover] [expr $y1-$cover-46.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $y1-$cover-52.5] [expr $z1-$cover] [expr $y1-$cover-52.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $y1-$cover-58.5] [expr $z1-$cover] [expr $y1-$cover-58.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $y1-$cover-64.5] [expr $z1-$cover] [expr $y1-$cover-64.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $y1-$cover-70.5] [expr $z1-$cover] [expr $y1-$cover-70.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $y1-$cover-76.5] [expr $z1-$cover] [expr $y1-$cover-76.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $y1-$cover-82.5] [expr $z1-$cover] [expr $y1-$cover-82.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $y1-$cover-88.5] [expr $z1-$cover] [expr $y1-$cover-88.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $y1-$cover-94.5] [expr $z1-$cover] [expr $y1-$cover-94.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $y1-$cover-100.5] [expr $z1-$cover] [expr $y1-$cover-100.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $y1-$cover-106.5] [expr $z1-$cover] [expr $y1-$cover-106.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $y1-$cover-112.5] [expr $z1-$cover] [expr $y1-$cover-112.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $y1-$cover-118.5] [expr $z1-$cover] [expr $y1-$cover-118.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $y1-$cover-124.5] [expr $z1-$cover] [expr $y1-$cover-124.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $y1-$cover-130.5] [expr $z1-$cover] [expr $y1-$cover-130.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $y1-$cover-136.5] [expr $z1-$cover] [expr $y1-$cover-136.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $y1-$cover-142.5] [expr $z1-$cover] [expr $y1-$cover-142.5] [expr $cover-$z1]
#left column longitudinal bar
layer straight 3 4 $Ase [expr $cover-$y1] [expr $z1-$cover] [expr $cover-$y1] [expr $cover-$z1]
layer straight 3 2 $Ase [expr $cover-$y1+6.375] [expr $z1-$cover] [expr $cover-$y1+6.375] [expr $cover-$z1]
layer straight 3 2 $Ase [expr $cover-$y1+12.75] [expr $z1-$cover] [expr $cover-$y1+12.75] [expr $cover-$z1]
layer straight 3 2 $Ase [expr $cover-$y1+19.125] [expr $z1-$cover] [expr $cover-$y1+19.125] [expr $cover-$z1]
layer straight 3 2 $Ase [expr $cover-$y1+25.5] [expr $z1-$cover] [expr $cover-$y1+25.5] [expr $cover-$z1]
layer straight 3 2 $Ase [expr $cover-$y1+31.875] [expr $z1-$cover] [expr $cover-$y1+31.875] [expr $cover-$z1]
layer straight 3 2 $Ase [expr $cover-$y1+38.25] [expr $z1-$cover] [expr $cover-$y1+38.25] [expr $cover-$z1]
layer straight 3 2 $Ase [expr $cover-$y1+44.625] [expr $z1-$cover] [expr $cover-$y1+44.625] [expr $cover-$z1]
layer straight 3 4 $Ase [expr $cover-$y1+51] [expr $z1-$cover] [expr $cover-$y1+51] [expr $cover-$z1]


#left wall longitudinal bar
layer straight 3 2 $Asf [expr $cover-$y1+58.5] [expr $z1-$cover] [expr $cover-$y1+58.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $cover-$y1+64.5] [expr $z1-$cover] [expr $cover-$y1+64.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $cover-$y1+70.5] [expr $z1-$cover] [expr $cover-$y1+70.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $cover-$y1+76.5] [expr $z1-$cover] [expr $cover-$y1+76.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $cover-$y1+82.5] [expr $z1-$cover] [expr $cover-$y1+82.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $cover-$y1+88.5] [expr $z1-$cover] [expr $cover-$y1+88.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $cover-$y1+94.5] [expr $z1-$cover] [expr $cover-$y1+94.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $cover-$y1+100.5] [expr $z1-$cover] [expr $cover-$y1+100.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $cover-$y1+106.5] [expr $z1-$cover] [expr $cover-$y1+106.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $cover-$y1+112.5] [expr $z1-$cover] [expr $cover-$y1+112.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $cover-$y1+118.5] [expr $z1-$cover] [expr $cover-$y1+118.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $cover-$y1+124.5] [expr $z1-$cover] [expr $cover-$y1+124.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $cover-$y1+130.5] [expr $z1-$cover] [expr $cover-$y1+130.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $cover-$y1+136.5] [expr $z1-$cover] [expr $cover-$y1+136.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr $cover-$y1+142.5] [expr $z1-$cover] [expr $cover-$y1+142.5] [expr $cover-$z1]
}



####another side of wall (right wing)


section Fiber 9 {
patch rect 2 10 1 -102 -$z1 102 $z1
patch rect 1 10 1 -$y1 -$z1 [expr 30-$y1] $z1
patch rect 1 10 1 [expr $y1-54] -$z1 $y1 $z1
# Create the reinforcing fibers (left, middle, right) first one is right

#left column longitudinal bar boundry B
###layer straight $matTag $numFiber $areaFiber $yStart $zStart $yEnd $zEnd
layer straight 3 4 $Asi [expr -$y1+$cover] [expr $z1-$cover] [expr -$y1+$cover] [expr $cover-$z1]
layer straight 3 2 $Asi [expr -$y1+$cover-9] [expr $z1-$cover] [expr -$y1+$cover-9] [expr $cover-$z1]
layer straight 3 2 $Asi [expr -$y1+$cover-18] [expr $z1-$cover] [expr -$y1+$cover-18] [expr $cover-$z1]
layer straight 3 4 $Asi [expr -$y1+$cover-27] [expr $z1-$cover] [expr -$y1+$cover-27] [expr $cover-$z1]

#left wall longitudinal bar

layer straight 3 4 $Asf [expr -$y1+$cover+34.5] [expr $z1-$cover] [expr -$y1+$cover+34.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$y1+$cover+40.5] [expr $z1-$cover] [expr -$y1+$cover+40.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$y1+$cover+46.5] [expr $z1-$cover] [expr -$y1+$cover+46.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$y1+$cover+52.5] [expr $z1-$cover] [expr -$y1+$cover+52.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$y1+$cover+58.5] [expr $z1-$cover] [expr -$y1+$cover+58.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$y1+$cover+64.5] [expr $z1-$cover] [expr -$y1+$cover+64.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$y1+$cover+70.5] [expr $z1-$cover] [expr -$y1+$cover+70.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$y1+$cover+76.5] [expr $z1-$cover] [expr -$y1+$cover+76.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$y1+$cover+82.5] [expr $z1-$cover] [expr -$y1+$cover+82.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$y1+$cover+88.5] [expr $z1-$cover] [expr -$y1+$cover+88.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$y1+$cover+94.5] [expr $z1-$cover] [expr -$y1+$cover+94.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$y1+$cover+100.5] [expr $z1-$cover] [expr -$y1+$cover+100.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$y1+$cover+106.5] [expr $z1-$cover] [expr -$y1+$cover+106.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$y1+$cover+112.5] [expr $z1-$cover] [expr -$y1+$cover+112.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$y1+$cover+118.5] [expr $z1-$cover] [expr -$y1+$cover+118.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$y1+$cover+124.5] [expr $z1-$cover] [expr -$y1+$cover+124.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$y1+$cover+130.5] [expr $z1-$cover] [expr -$y1+$cover+130.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$y1+$cover+136.5] [expr $z1-$cover] [expr -$y1+$cover+136.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$y1+$cover+142.5] [expr $z1-$cover] [expr -$y1+$cover+142.5] [expr $cover-$z1]
#right column longitudinal bar Boundry A
layer straight 3 4 $Ase [expr -$cover+$y1] [expr $z1-$cover] [expr -$cover+$y1] [expr $cover-$z1]
layer straight 3 2 $Ase [expr -$cover+$y1-6.375] [expr $z1-$cover] [expr -$cover+$y1-6.375] [expr $cover-$z1]
layer straight 3 2 $Ase [expr -$cover+$y1-12.75] [expr $z1-$cover] [expr -$cover+$y1-12.75] [expr $cover-$z1]
layer straight 3 2 $Ase [expr -$cover+$y1-19.125] [expr $z1-$cover] [expr -$cover+$y1-19.125] [expr $cover-$z1]
layer straight 3 2 $Ase [expr -$cover+$y1-25.5] [expr $z1-$cover] [expr -$cover+$y1-25.5] [expr $cover-$z1]
layer straight 3 2 $Ase [expr -$cover+$y1-31.875] [expr $z1-$cover] [expr -$cover+$y1-31.875] [expr $cover-$z1]
layer straight 3 2 $Ase [expr -$cover+$y1-38.25] [expr $z1-$cover] [expr -$cover+$y1-38.25] [expr $cover-$z1]
layer straight 3 2 $Ase [expr -$cover+$y1-44.625] [expr $z1-$cover] [expr -$cover+$y1-44.625] [expr $cover-$z1]
layer straight 3 4 $Ase [expr -$cover+$y1-51] [expr $z1-$cover] [expr -$cover+$y1-51] [expr $cover-$z1]


#right wall longitudinal bar
layer straight 3 2 $Asf [expr -$cover+$y1-58.5] [expr $z1-$cover] [expr -$cover+$y1-58.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$cover+$y1-64.5] [expr $z1-$cover] [expr -$cover+$y1-64.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$cover+$y1-70.5] [expr $z1-$cover] [expr -$cover+$y1-70.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$cover+$y1-76.5] [expr $z1-$cover] [expr -$cover+$y1-76.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$cover+$y1-82.5] [expr $z1-$cover] [expr -$cover+$y1-82.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$cover+$y1-88.5] [expr $z1-$cover] [expr -$cover+$y1-88.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$cover+$y1-94.5] [expr $z1-$cover] [expr -$cover+$y1-94.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$cover+$y1-100.5] [expr $z1-$cover] [expr -$cover+$y1-100.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$cover+$y1-106.5] [expr $z1-$cover] [expr -$cover+$y1-106.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$cover+$y1-112.5] [expr $z1-$cover] [expr -$cover+$y1-112.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$cover+$y1-118.5] [expr $z1-$cover] [expr -$cover+$y1-118.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$cover+$y1-124.5] [expr $z1-$cover] [expr -$cover+$y1-124.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$cover+$y1-130.5] [expr $z1-$cover] [expr -$cover+$y1-130.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$cover+$y1-136.5] [expr $z1-$cover] [expr -$cover+$y1-136.5] [expr $cover-$z1]
layer straight 3 2 $Asf [expr -$cover+$y1-142.5] [expr $z1-$cover] [expr -$cover+$y1-142.5] [expr $cover-$z1]
}

###########################beam of 1st story#############################################################################

set beamWidth 60
set beamDepth 24



set y2 [expr $beamDepth/2.0]
set z2 [expr $beamWidth/2.0]


section Fiber 2 {

# Create the concrete core fibers confined sma start
patch rect 1 10 1 -$y2 -$z2 $y2 $z2

# Create the reinforcing fibers (left, middle, right)
layer straight 3 8 $Asf [expr $y2-$cover] [expr $z2-$cover-6] [expr $y2-$cover] [expr $cover-$z2+6]
layer straight 3 2 $Asf 3 [expr $z2-$cover] 3 [expr $cover-$z2]
layer straight 3 2 $Asf -3 [expr $z2-$cover] -3 [expr $cover-$z2]
layer straight 3 8 $Asf [expr $cover-$y2] [expr $z2-$cover-6] [expr $cover-$y2] [expr $cover-$z2+6]
layer straight 3 2 $Asf [expr $y2-$cover] [expr $z2-$cover] [expr $y2-$cover] [expr $cover-$z2]
layer straight 3 2 $Asf [expr $cover-$y2] [expr $z2-$cover] [expr $cover-$y2] [expr $cover-$z2]
layer straight 3 2 $Ase [expr $y2-$cover-6] [expr -$z2+$cover+6] [expr $y2-$cover-6] [expr $cover-$z2+12]
layer straight 3 2 $Ase [expr $y2-$cover-12] [expr -$z2+$cover+6] [expr $y2-$cover-12] [expr $cover-$z2+12]
layer straight 3 2 $Ase [expr -$y2+$cover+6] [expr $z2-$cover-6] [expr -$y2+$cover+6] [expr -$cover+$z2-12]
layer straight 3 2 $Ase [expr -$y2+$cover+12] [expr $z2-$cover-6] [expr -$y2+$cover+12] [expr -$cover+$z2-12]
}

section Fiber 3 {

# Create the concrete core fibers confined sma finish
patch rect 1 10 1 -$y2 -$z2 $y2 $z2

# Create the reinforcing fibers (left, middle, right)
layer straight 3 8 $Asf [expr $y2-$cover] [expr $z2-$cover-6] [expr $y2-$cover] [expr $cover-$z2+6]
layer straight 3 2 $Asf 3 [expr $z2-$cover] 3 [expr $cover-$z2]
layer straight 3 2 $Asf -3 [expr $z2-$cover] -3 [expr $cover-$z2]
layer straight 3 8 $Asf [expr $cover-$y2] [expr $z2-$cover-6] [expr $cover-$y2] [expr $cover-$z2+6]
layer straight 3 2 $Asf [expr $y2-$cover] [expr $z2-$cover] [expr $y2-$cover] [expr $cover-$z2]
layer straight 3 2 $Asf [expr $cover-$y2] [expr $z2-$cover] [expr $cover-$y2] [expr $cover-$z2]
layer straight 3 2 $Ase [expr $y2-$cover-6] [expr -$z2+$cover+20] [expr $y2-$cover-6] [expr $cover-$z2+26]
layer straight 3 2 $Ase [expr $y2-$cover-12] [expr -$z2+$cover+20] [expr $y2-$cover-12] [expr $cover-$z2+26]
layer straight 3 2 $Ase [expr -$y2+$cover+6] [expr $z2-$cover-20] [expr -$y2+$cover+6] [expr -$cover+$z2-26]
layer straight 3 2 $Ase [expr -$y2+$cover+12] [expr $z2-$cover-20] [expr -$y2+$cover+12] [expr -$cover+$z2-26]
}

section Fiber 4 {

# Create the concrete core fibers confined coupler start
patch rect 1 10 1 -$y2 -$z2 $y2 $z2

# Create the reinforcing fibers (left, middle, right)
layer straight 3 8 $Asf [expr $y2-$cover] [expr $z2-$cover-6] [expr $y2-$cover] [expr $cover-$z2+6]
layer straight 3 2 $Asf 3 [expr $z2-$cover] 3 [expr $cover-$z2]
layer straight 3 2 $Asf -3 [expr $z2-$cover] -3 [expr $cover-$z2]
layer straight 3 8 $Asf [expr $cover-$y2] [expr $z2-$cover-6] [expr $cover-$y2] [expr $cover-$z2+6]
layer straight 3 2 $Asf [expr $y2-$cover] [expr $z2-$cover] [expr $y2-$cover] [expr $cover-$z2]
layer straight 3 2 $Asf [expr $cover-$y2] [expr $z2-$cover] [expr $cover-$y2] [expr $cover-$z2]
layer straight 3 2 $Ase [expr $y2-$cover-6] [expr -$z2+$cover+20] [expr $y2-$cover-6] [expr $cover-$z2+26]
layer straight 3 2 $Ase [expr $y2-$cover-12] [expr -$z2+$cover+20] [expr $y2-$cover-12] [expr $cover-$z2+26]
layer straight 3 2 $Ase [expr -$y2+$cover+6] [expr $z2-$cover-20] [expr -$y2+$cover+6] [expr -$cover+$z2-26]
layer straight 3 2 $Ase [expr -$y2+$cover+12] [expr $z2-$cover-20] [expr -$y2+$cover+12] [expr -$cover+$z2-26]
}

section Fiber 5 {

# Create the concrete core fibers confined coupler finish
patch rect 1 10 1 -$y2 -$z2 $y2 $z2

# Create the reinforcing fibers (left, middle, right)
layer straight 3 8 $Asf [expr $y2-$cover] [expr $z2-$cover-6] [expr $y2-$cover] [expr $cover-$z2+6]
layer straight 3 2 $Asf 3 [expr $z2-$cover] 3 [expr $cover-$z2]
layer straight 3 2 $Asf -3 [expr $z2-$cover] -3 [expr $cover-$z2]
layer straight 3 8 $Asf [expr $cover-$y2] [expr $z2-$cover-6] [expr $cover-$y2] [expr $cover-$z2+6]
layer straight 3 2 $Asf [expr $y2-$cover] [expr $z2-$cover] [expr $y2-$cover] [expr $cover-$z2]
layer straight 3 2 $Asf [expr $cover-$y2] [expr $z2-$cover] [expr $cover-$y2] [expr $cover-$z2]
layer straight 3 2 $Ase [expr $y2-$cover-6] [expr -$z2+$cover+21.9] [expr $y2-$cover-6] [expr $cover-$z2+27.9]
layer straight 3 2 $Ase [expr $y2-$cover-12] [expr -$z2+$cover+21.9] [expr $y2-$cover-12] [expr $cover-$z2+27.9]
layer straight 3 2 $Ase [expr -$y2+$cover+6] [expr $z2-$cover-21.9] [expr -$y2+$cover+6] [expr -$cover+$z2-27.9]
layer straight 3 2 $Ase [expr -$y2+$cover+12] [expr $z2-$cover-21.9] [expr -$y2+$cover+12] [expr -$cover+$z2-27.9]
}
section Fiber 6 {

# Create the concrete core fibers only bar start
patch rect 1 10 1 -$y2 -$z2 $y2 $z2

# Create the reinforcing fibers (left, middle, right)
layer straight 3 10 $Asf [expr $y2-$cover] [expr $z2-$cover-6] [expr $y2-$cover] [expr $cover-$z2+6]
layer straight 3 2 $Asf 3 [expr $z2-$cover] 3 [expr $cover-$z2]
layer straight 3 2 $Asf -3 [expr $z2-$cover] -3 [expr $cover-$z2]
layer straight 3 10 $Asf [expr $cover-$y2] [expr $z2-$cover-6] [expr $cover-$y2] [expr $cover-$z2+6]
layer straight 3 2 $Ase [expr $y2-$cover-6] [expr -$z2+$cover+21.9] [expr $y2-$cover-6] [expr $cover-$z2+27.9]
layer straight 3 2 $Ase [expr $y2-$cover-12] [expr -$z2+$cover+21.9] [expr $y2-$cover-12] [expr $cover-$z2+27.9]
layer straight 3 2 $Ase [expr -$y2+$cover+6] [expr $z2-$cover-21.9] [expr -$y2+$cover+6] [expr -$cover+$z2-27.9]
layer straight 3 2 $Ase [expr -$y2+$cover+12] [expr $z2-$cover-21.9] [expr -$y2+$cover+12] [expr -$cover+$z2-27.9]

}

section Fiber 7 {

# Create the concrete core fibers only bar finish
patch rect 1 10 1 -$y2 -$z2 $y2 $z2

# Create the reinforcing fibers (left, middle, right)
layer straight 3 10 $Asf [expr $y2-$cover] [expr $z2-$cover-6] [expr $y2-$cover] [expr $cover-$z2+6]
layer straight 3 2 $Asf 3 [expr $z2-$cover] 3 [expr $cover-$z2]
layer straight 3 2 $Asf -3 [expr $z2-$cover] -3 [expr $cover-$z2]
layer straight 3 10 $Asf [expr $cover-$y2] [expr $z2-$cover-6] [expr $cover-$y2] [expr $cover-$z2+6]
layer straight 3 2 $Ase [expr $y2-$cover-6] [expr -$z2+$cover+28.75] [expr $y2-$cover-6] [expr $cover-$z2+34.75]
layer straight 3 2 $Ase [expr $y2-$cover-12] [expr -$z2+$cover+28.75] [expr $y2-$cover-12] [expr $cover-$z2+34.75]
layer straight 3 2 $Ase [expr -$y2+$cover+6] [expr $z2-$cover-28.75] [expr -$y2+$cover+6] [expr -$cover+$z2-34.75]
layer straight 3 2 $Ase [expr -$y2+$cover+12] [expr $z2-$cover-28.75] [expr -$y2+$cover+12] [expr -$cover+$z2-34.75]

}

section Fiber 8 {

# Create the concrete core fibers only bar middle
patch rect 1 10 1 -$y2 -$z2 $y2 $z2

# Create the reinforcing fibers (left, middle, right)
layer straight 3 10 $Asf [expr $y2-$cover] [expr $z2-$cover-6] [expr $y2-$cover] [expr $cover-$z2+6]
layer straight 3 2 $Asf 3 [expr $z2-$cover] 3 [expr $cover-$z2]
layer straight 3 2 $Asf -3 [expr $z2-$cover] -3 [expr $cover-$z2]
layer straight 3 10 $Asf [expr $cover-$y2] [expr $z2-$cover-6] [expr $cover-$y2] [expr $cover-$z2+6]
layer straight 3 2 $Ase [expr $y2-$cover-6] 0 [expr $y2-$cover-6] 0
layer straight 3 2 $Ase [expr $y2-$cover-12] 0 [expr $y2-$cover-12] 0
layer straight 3 2 $Ase [expr -$y2+$cover+6] 0 [expr -$y2+$cover+6] 0
layer straight 3 2 $Ase [expr -$y2+$cover+12] 0 [expr -$y2+$cover+12] 0
}


section Fiber 26 {

# Create the concrete core fibers confined sma start
patch rect 1 10 1 -$y2 -$z2 $y2 $z2

# Create the reinforcing fibers (left, middle, right)
layer straight 3 8 $Asf [expr $y2-$cover] [expr $z2-$cover-6] [expr $y2-$cover] [expr $cover-$z2+6]
layer straight 3 2 $Asf 3 [expr $z2-$cover] 3 [expr $cover-$z2]
layer straight 3 2 $Asf -3 [expr $z2-$cover] -3 [expr $cover-$z2]
layer straight 3 8 $Asf [expr $cover-$y2] [expr $z2-$cover-6] [expr $cover-$y2] [expr $cover-$z2+6]
layer straight 3 2 $Asf [expr $y2-$cover] [expr $z2-$cover] [expr $y2-$cover] [expr $cover-$z2]
layer straight 3 2 $Asf [expr $cover-$y2] [expr $z2-$cover] [expr $cover-$y2] [expr $cover-$z2]
layer straight 3 2 $Ase [expr -$y2+$cover+6] [expr -$z2+$cover+6] [expr -$y2+$cover+6] [expr $cover-$z2+12]
layer straight 3 2 $Ase [expr -$y2+$cover+12] [expr -$z2+$cover+6] [expr -$y2+$cover+12] [expr $cover-$z2+12]
layer straight 3 2 $Ase [expr $y2-$cover-6] [expr $z2-$cover-6] [expr $y2-$cover-6] [expr -$cover+$z2-12]
layer straight 3 2 $Ase [expr $y2-$cover-12] [expr $z2-$cover-6] [expr $y2-$cover-12] [expr -$cover+$z2-12]
}

section Fiber 27 {

# Create the concrete core fibers confined sma finish
patch rect 1 10 1 -$y2 -$z2 $y2 $z2

# Create the reinforcing fibers (left, middle, right)
layer straight 3 8 $Asf [expr $y2-$cover] [expr $z2-$cover-6] [expr $y2-$cover] [expr $cover-$z2+6]
layer straight 3 2 $Asf 3 [expr $z2-$cover] 3 [expr $cover-$z2]
layer straight 3 2 $Asf -3 [expr $z2-$cover] -3 [expr $cover-$z2]
layer straight 3 8 $Asf [expr $cover-$y2] [expr $z2-$cover-6] [expr $cover-$y2] [expr $cover-$z2+6]
layer straight 3 2 $Asf [expr $y2-$cover] [expr $z2-$cover] [expr $y2-$cover] [expr $cover-$z2]
layer straight 3 2 $Asf [expr $cover-$y2] [expr $z2-$cover] [expr $cover-$y2] [expr $cover-$z2]
layer straight 3 2 $Ase [expr -$y2+$cover+6] [expr -$z2+$cover+20] [expr -$y2+$cover+6] [expr $cover-$z2+26]
layer straight 3 2 $Ase [expr -$y2+$cover+12] [expr -$z2+$cover+20] [expr -$y2+$cover+12] [expr $cover-$z2+26]
layer straight 3 2 $Ase [expr $y2-$cover-6] [expr $z2-$cover-20] [expr $y2-$cover-6] [expr -$cover+$z2-26]
layer straight 3 2 $Ase [expr $y2-$cover-12] [expr $z2-$cover-20] [expr $y2-$cover-12] [expr -$cover+$z2-26]
}

section Fiber 28 {

# Create the concrete core fibers confined coupler start
patch rect 1 10 1 -$y2 -$z2 $y2 $z2

# Create the reinforcing fibers (left, middle, right)
layer straight 3 8 $Asf [expr $y2-$cover] [expr $z2-$cover-6] [expr $y2-$cover] [expr $cover-$z2+6]
layer straight 3 2 $Asf 3 [expr $z2-$cover] 3 [expr $cover-$z2]
layer straight 3 2 $Asf -3 [expr $z2-$cover] -3 [expr $cover-$z2]
layer straight 3 8 $Asf [expr $cover-$y2] [expr $z2-$cover-6] [expr $cover-$y2] [expr $cover-$z2+6]
layer straight 3 2 $Asf [expr $y2-$cover] [expr $z2-$cover] [expr $y2-$cover] [expr $cover-$z2]
layer straight 3 2 $Asf [expr $cover-$y2] [expr $z2-$cover] [expr $cover-$y2] [expr $cover-$z2]
layer straight 3 2 $Ase [expr -$y2+$cover+6] [expr -$z2+$cover+20] [expr -$y2+$cover+6] [expr $cover-$z2+26]
layer straight 3 2 $Ase [expr -$y2+$cover+12] [expr -$z2+$cover+20] [expr -$y2+$cover+12] [expr $cover-$z2+26]
layer straight 3 2 $Ase [expr $y2-$cover-6] [expr $z2-$cover-20] [expr $y2-$cover-6] [expr -$cover+$z2-26]
layer straight 3 2 $Ase [expr $y2-$cover-12] [expr $z2-$cover-20] [expr $y2-$cover-12] [expr -$cover+$z2-26]
}

section Fiber 29 {

# Create the concrete core fibers confined coupler finish
patch rect 1 10 1 -$y2 -$z2 $y2 $z2

# Create the reinforcing fibers (left, middle, right)
layer straight 3 8 $Asf [expr $y2-$cover] [expr $z2-$cover-6] [expr $y2-$cover] [expr $cover-$z2+6]
layer straight 3 2 $Asf 3 [expr $z2-$cover] 3 [expr $cover-$z2]
layer straight 3 2 $Asf -3 [expr $z2-$cover] -3 [expr $cover-$z2]
layer straight 3 8 $Asf [expr $cover-$y2] [expr $z2-$cover-6] [expr $cover-$y2] [expr $cover-$z2+6]
layer straight 3 2 $Asf [expr $y2-$cover] [expr $z2-$cover] [expr $y2-$cover] [expr $cover-$z2]
layer straight 3 2 $Asf [expr $cover-$y2] [expr $z2-$cover] [expr $cover-$y2] [expr $cover-$z2]
layer straight 3 2 $Ase [expr -$y2+$cover+6] [expr -$z2+$cover+21.9] [expr -$y2+$cover+6] [expr $cover-$z2+27.9]
layer straight 3 2 $Ase [expr -$y2+$cover+12] [expr -$z2+$cover+21.9] [expr -$y2+$cover+12] [expr $cover-$z2+27.9]
layer straight 3 2 $Ase [expr $y2-$cover-6] [expr $z2-$cover-21.9] [expr $y2-$cover-6] [expr -$cover+$z2-27.9]
layer straight 3 2 $Ase [expr $y2-$cover-12] [expr $z2-$cover-21.9] [expr $y2-$cover-12] [expr -$cover+$z2-27.9]
}
section Fiber 30 {

# Create the concrete core fibers only bar start
patch rect 1 10 1 -$y2 -$z2 $y2 $z2

# Create the reinforcing fibers (left, middle, right)
layer straight 3 10 $Asf [expr $y2-$cover] [expr $z2-$cover-6] [expr $y2-$cover] [expr $cover-$z2+6]
layer straight 3 2 $Asf 3 [expr $z2-$cover] 3 [expr $cover-$z2]
layer straight 3 2 $Asf -3 [expr $z2-$cover] -3 [expr $cover-$z2]
layer straight 3 10 $Asf [expr $cover-$y2] [expr $z2-$cover-6] [expr $cover-$y2] [expr $cover-$z2+6]
layer straight 3 2 $Ase [expr -$y2+$cover+6] [expr -$z2+$cover+21.9] [expr -$y2+$cover+6] [expr $cover-$z2+27.9]
layer straight 3 2 $Ase [expr -$y2+$cover+12] [expr -$z2+$cover+21.9] [expr -$y2+$cover+12] [expr $cover-$z2+27.9]
layer straight 3 2 $Ase [expr $y2-$cover-6] [expr $z2-$cover-21.9] [expr $y2-$cover-6] [expr -$cover+$z2-27.9]
layer straight 3 2 $Ase [expr $y2-$cover-12] [expr $z2-$cover-21.9] [expr $y2-$cover-12] [expr -$cover+$z2-27.9]

}

section Fiber 31 {

# Create the concrete core fibers only bar finish
patch rect 1 10 1 -$y2 -$z2 $y2 $z2

# Create the reinforcing fibers (left, middle, right)
layer straight 3 10 $Asf [expr $y2-$cover] [expr $z2-$cover-6] [expr $y2-$cover] [expr $cover-$z2+6]
layer straight 3 2 $Asf 3 [expr $z2-$cover] 3 [expr $cover-$z2]
layer straight 3 2 $Asf -3 [expr $z2-$cover] -3 [expr $cover-$z2]
layer straight 3 10 $Asf [expr $cover-$y2] [expr $z2-$cover-6] [expr $cover-$y2] [expr $cover-$z2+6]
layer straight 3 2 $Ase [expr -$y2+$cover+6] [expr -$z2+$cover+28.75] [expr -$y2+$cover+6] [expr $cover-$z2+34.75]
layer straight 3 2 $Ase [expr -$y2+$cover+12] [expr -$z2+$cover+28.75] [expr -$y2+$cover+12] [expr $cover-$z2+34.75]
layer straight 3 2 $Ase [expr $y2-$cover-6] [expr $z2-$cover-28.75] [expr $y2-$cover-6] [expr -$cover+$z2-34.75]
layer straight 3 2 $Ase [expr $y2-$cover-12] [expr $z2-$cover-28.75] [expr $y2-$cover-12] [expr -$cover+$z2-34.75]

}

geomTransf Linear 1

# Number of integration points along length of element
set np 5

##Create the coulumns using Beam-column elements
##element dispBeamColumn $eleTag $iNode $jNode $numIntgrPts $secTag $transfTag
element dispBeamColumn 1 1 10 $np 1 1
element dispBeamColumn 2 10 11 $np 1 1
element dispBeamColumn 3 11 12 $np 1 1
element dispBeamColumn 4 12 13 $np 1 1
element dispBeamColumn 5 13 14 $np 1 1
element dispBeamColumn 6 14 101 $np 1 1

element dispBeamColumn 7 2 15 $np 9 1
element dispBeamColumn 8 15 16 $np 9 1
element dispBeamColumn 9 16 17 $np 9 1
element dispBeamColumn 10 17 18 $np 9 1
element dispBeamColumn 11 18 19 $np 9 1
element dispBeamColumn 12 19 111 $np 9 1

geomTransf Linear 2
#rigidLink $type $rNodeTag $cNodeTag
rigidLink beam 101 102

##element truss 301 101 102 1e10 1
#element dispBeamColumn $eleTag $iNode $jNode $numIntgrPts -sections $secTag1 $secTag2 ... $transfTag


element dispBeamColumn 13 1021 103 $np -sections 2 3 2
element dispBeamColumn 14 1031 104 $np -sections 4 5 2
element dispBeamColumn 15 1041 105 $np -sections 6 7 2
element dispBeamColumn 16 1051 106 $np -sections 7 8 2

element dispBeamColumn 17 1106 107 $np -sections 8 31 2
element dispBeamColumn 18 1107 1081 $np -sections 31 30 2
element dispBeamColumn 19 108 1091 $np -sections 29 28 2
element dispBeamColumn 20 109 1101 $np -sections 27 26 2


####element truss 302 110 111 1e10 1
rigidLink beam 110 111
#

set Pt 236.5;
set pr 171 # 10% of axial capacity of columns

# Create a Plain load pattern with a Linear TimeSeries
timeSeries Linear 1
pattern Plain 1 1 {

# Create nodal loads at nodes 3 & 4
# nd FX FY MZ
load 101 0.0 [expr -$Pt] 0.0
load 111 0.0 [expr -$Pt] 0.0

}
recorder Node -file data/[expr $i./10]pga/pga[expr $i./10]node4disp.out -time -node 111 -dof 1 3 5 disp
recorder Node -file data/[expr $i./10]pga/pga[expr $i./10]node1reac.out -time -node 1 -dof 1 2 3 reaction
recorder Node -file data/[expr $i./10]pga/pga[expr $i./10]node2reac.out -time -node 2 -dof 1 2 3 reaction

recorder Drift -file data/[expr $i./10]pga/pga[expr $i./10]drift.out -time -iNode 1 -jNode 3 -dof 1 -perpDirn 2

recorder display frame 200 0 500 500 -wipe
vup 0 1 0 # dirn defining up direction of view plane(ghab bordar movazi)
vpn 0 0 1 # direction of outward normal to view plane(amud bar safhe ghab)
prp 0 0 1 # eye location in local coord sys defined by viewing system
viewWindow -2 2 -2 2
display 1 5 20


set accel "Series -dt 0.02 -filePath 125-7277-005.txt -factor [expr $i.*9.81/10]"


pattern UniformExcitation 3 1 -accel $accel


rayleigh 1.1938628480692413 0.0 0.0 0.0020940428827675574;


puts "groundmotion start!.Time: [getTime]"

set dt 0.02

constraints Plain
numberer Plain
system BandGeneral
test NormDispIncr 1.0e-8 10
algorithm Newton
integrator Newmark 0.5 0.25
analysis Transient
analyze [expr int(32.82/$dt)] $dt

puts "groundmotion [expr $i./10]pga done!.End Time: [getTime]"
Borislav
Posts: 4
Joined: Tue May 07, 2019 10:34 am

Re: can,t run

Post by Borislav »

looks to me like you haven't defined any masses and applied gravity loading
Post Reply