difference between T(Opensees) & T(Sap) in 3D model!

[aminamini]

Hi everybody!

I've modeled a 3-D, 10-story steel moment frame using nonlinearBeamColumn elements. the model has 2 bays in X and Y directions. I,ve taken into account the torsional stiffness by section aggregator. the first three periods of structure are: 1.208, 1.168, 1.137. But, they are 1.700, 1.455, and 1.113 in Sap model. I don't now what's the problem. can anyone help me?

the shear stiffness has been removed in Sap model.

Here is my model:

wipe;
model basic -ndm 3 -ndf 6
file mkdir data

# Create nodes
# ------------
for {set x 0} {$x<=10} {incr x 1} {
node [expr 9*$x+1] 0. 0. [expr ($x)*3.2]
node [expr 9*$x+2] 5. 0. [expr ($x)*3.2]
node [expr 9*$x+3] 10. 0. [expr ($x)*3.2]
node [expr 9*$x+4] 0. 5. [expr ($x)*3.2]
node [expr 9*$x+5] 5. 5. [expr ($x)*3.2]
node [expr 9*$x+6] 10. 5. [expr ($x)*3.2]
node [expr 9*$x+7] 0. 10. [expr ($x)*3.2]
node [expr 9*$x+8] 5. 10. [expr ($x)*3.2]
node [expr 9*$x+9] 10. 10. [expr ($x)*3.2]
}

mass 10 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 11 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 12 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 13 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 14 7135.58 7135.58 1.e-9 1.e-9 1.e-9 [expr (1./12)*(70000./9.81)*(pow(10,2)+pow(10,2))]
mass 15 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 16 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 17 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 18 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 19 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 20 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 21 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 22 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 23 7135.58 7135.58 1.e-9 1.e-9 1.e-9 [expr (1./12)*(70000./9.81)*(pow(10,2)+pow(10,2))]
mass 24 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 25 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 26 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 27 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 28 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 29 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 30 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 31 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 32 7135.58 7135.58 1.e-9 1.e-9 1.e-9 [expr (1./12)*(70000./9.81)*(pow(10,2)+pow(10,2))]
mass 33 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 34 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 35 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 36 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 37 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 38 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 39 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 40 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 41 7135.58 7135.58 1.e-9 1.e-9 1.e-9 [expr (1./12)*(70000./9.81)*(pow(10,2)+pow(10,2))]
mass 42 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 43 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 44 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 45 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 46 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 47 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 48 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 49 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 50 7135.58 7135.58 1.e-9 1.e-9 1.e-9 [expr (1./12)*(70000./9.81)*(pow(10,2)+pow(10,2))]
mass 51 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 52 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 53 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 54 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 55 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 56 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 57 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 58 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 59 7135.58 7135.58 1.e-9 1.e-9 1.e-9 [expr (1./12)*(70000./9.81)*(pow(10,2)+pow(10,2))]
mass 60 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 61 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 62 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 63 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 64 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 65 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 66 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 67 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 68 7135.58 7135.58 1.e-9 1.e-9 1.e-9 [expr (1./12)*(70000./9.81)*(pow(10,2)+pow(10,2))]
mass 69 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 70 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 71 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 72 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 73 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 74 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 75 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 76 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 77 7135.58 7135.58 1.e-9 1.e-9 1.e-9 [expr (1./12)*(70000./9.81)*(pow(10,2)+pow(10,2))]
mass 78 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 79 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 80 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 81 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 82 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 83 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 84 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 85 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 86 7135.58 7135.58 1.e-9 1.e-9 1.e-9 [expr (1./12)*(70000./9.81)*(pow(10,2)+pow(10,2))]
mass 87 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 88 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 89 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 90 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 91 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 92 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 93 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 94 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 95 7135.58 7135.58 1.e-9 1.e-9 1.e-9 [expr (1./12)*(70000./9.81)*(pow(10,2)+pow(10,2))]
mass 96 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 97 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 98 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9
mass 99 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9 1.e-9

fix 1 1 1 1 1 1 1
fix 2 1 1 1 1 1 1
fix 3 1 1 1 1 1 1
fix 4 1 1 1 1 1 1
fix 5 1 1 1 1 1 1
fix 6 1 1 1 1 1 1
fix 7 1 1 1 1 1 1
fix 8 1 1 1 1 1 1
fix 9 1 1 1 1 1 1

fix 14 0 0 1 1 1 0
fix 23 0 0 1 1 1 0
fix 32 0 0 1 1 1 0
fix 41 0 0 1 1 1 0
fix 50 0 0 1 1 1 0
fix 59 0 0 1 1 1 0
fix 68 0 0 1 1 1 0
fix 77 0 0 1 1 1 0
fix 86 0 0 1 1 1 0
fix 95 0 0 1 1 1 0

rigidDiaphragm 3 14 10 11 12 13 15 16 17 18
rigidDiaphragm 3 23 19 20 21 22 24 25 26 27
rigidDiaphragm 3 32 28 29 30 31 33 34 35 36
rigidDiaphragm 3 41 37 38 39 40 42 43 44 45
rigidDiaphragm 3 50 46 47 48 49 51 52 53 54
rigidDiaphragm 3 59 55 56 57 58 60 61 62 63
rigidDiaphragm 3 68 64 65 66 67 69 70 71 72
rigidDiaphragm 3 77 73 74 75 76 78 79 80 81
rigidDiaphragm 3 86 82 83 84 85 87 88 89 90
rigidDiaphragm 3 95 91 92 93 94 96 97 98 99

set ColumnTransfTag 1
geomTransf PDelta $ColumnTransfTag 1 0 0

set BeamTransfTag 2
geomTransf Linear $BeamTransfTag 0 -1 0

set GirderTransfTag 3
geomTransf Linear $GirderTransfTag 1 0 0

uniaxialMaterial Steel01 1 24000000.0 21000000000.0 0.005

# SC4(column 40X2.5)
# set some paramaters
set h 0.4
set t .025
set y [expr $h/2.0]
section Fiber 1 {
# create top,bot,left,right
patch quad 1 30 4 [expr -$y] [expr $y] [expr -$y] [expr $y-$t] [expr $y] [expr $y-$t] [expr $y] [expr $y]
patch quad 1 30 4 [expr -$y] [expr -$y+$t] [expr -$y] [expr -$y] [expr $y] [expr -$y] [expr $y] [expr -$y+$t]
patch quad 1 4 25 [expr -$y] [expr $y-$t] [expr -$y] [expr -$y+$t] [expr -$y+$t] [expr -$y+$t] [expr -$y+$t] [expr $y-$t]
patch quad 1 4 25 [expr $y-$t] [expr $y-$t] [expr $y-$t] [expr -$y+$t] [expr $y] [expr -$y+$t] [expr $y] [expr $y-$t]
}

# SC3(column 35X2.5)
# set some paramaters
set h 0.35
set t .025
set y [expr $h/2.0]
section Fiber 2 {
# create top,bot,left,right
patch quad 1 30 4 [expr -$y] [expr $y] [expr -$y] [expr $y-$t] [expr $y] [expr $y-$t] [expr $y] [expr $y]
patch quad 1 30 4 [expr -$y] [expr -$y+$t] [expr -$y] [expr -$y] [expr $y] [expr -$y] [expr $y] [expr -$y+$t]
patch quad 1 4 25 [expr -$y] [expr $y-$t] [expr -$y] [expr -$y+$t] [expr -$y+$t] [expr -$y+$t] [expr -$y+$t] [expr $y-$t]
patch quad 1 4 25 [expr $y-$t] [expr $y-$t] [expr $y-$t] [expr -$y+$t] [expr $y] [expr -$y+$t] [expr $y] [expr $y-$t]
}

# SC2(column 30X2)
# set some paramaters
set h 0.3
set t .02
set y [expr $h/2.0]
section Fiber 3 {
# create top,bot,left,right
patch quad 1 30 4 [expr -$y] [expr $y] [expr -$y] [expr $y-$t] [expr $y] [expr $y-$t] [expr $y] [expr $y]
patch quad 1 30 4 [expr -$y] [expr -$y+$t] [expr -$y] [expr -$y] [expr $y] [expr -$y] [expr $y] [expr -$y+$t]
patch quad 1 4 25 [expr -$y] [expr $y-$t] [expr -$y] [expr -$y+$t] [expr -$y+$t] [expr -$y+$t] [expr -$y+$t] [expr $y-$t]
patch quad 1 4 25 [expr $y-$t] [expr $y-$t] [expr $y-$t] [expr -$y+$t] [expr $y] [expr -$y+$t] [expr $y] [expr $y-$t]
}

# SB5(beam 40 1 22.5 2)
# set some paramaters
set ht 0.4
set tw .01
set bf .225
set tf .02
set y1 [expr $ht/2.0]
set y2 [expr $bf/2.0]
set y3 [expr $y1-$tf]
set y4 [expr $tw/2.0]
section Fiber 4 {
# create bot flange,web, top flange
patch quad 1 4 20 [expr -$y1] [expr $y2] [expr -$y1] [expr -$y2] [expr -$y3] [expr -$y2] [expr -$y3] [expr $y2]
patch quad 1 30 2 [expr -$y3] [expr $y4] [expr -$y3] [expr -$y4] [expr $y3] [expr -$y4] [expr $y3] [expr $y4]
patch quad 1 4 20 [expr $y3] [expr $y2] [expr $y3] [expr -$y2] [expr $y1] [expr -$y2] [expr $y1] [expr $y2]
}

# SB4(beam 35 .8 22.5 2)
# set some paramaters
set ht 0.35
set tw .008
set bf .225
set tf .02
set y1 [expr $ht/2.0]
set y2 [expr $bf/2.0]
set y3 [expr $y1-$tf]
set y4 [expr $tw/2.0]
section Fiber 5 {
# create bot flange,web, top flange
patch quad 1 4 20 [expr -$y1] [expr $y2] [expr -$y1] [expr -$y2] [expr -$y3] [expr -$y2] [expr -$y3] [expr $y2]
patch quad 1 30 2 [expr -$y3] [expr $y4] [expr -$y3] [expr -$y4] [expr $y3] [expr -$y4] [expr $y3] [expr $y4]
patch quad 1 4 20 [expr $y3] [expr $y2] [expr $y3] [expr -$y2] [expr $y1] [expr -$y2] [expr $y1] [expr $y2]
}

# SB3(beam 30 0.8 20 1.5)
# set some paramaters
set ht 0.3
set tw .008
set bf .2
set tf .015
set y1 [expr $ht/2.0]
set y2 [expr $bf/2.0]
set y3 [expr $y1-$tf]
set y4 [expr $tw/2.0]
section Fiber 6
# create bot flange,web, top flange
patch quad 1 4 20 [expr -$y1] [expr $y2] [expr -$y1] [expr -$y2] [expr -$y3] [expr -$y2] [expr -$y3] [expr $y2]
patch quad 1 30 2 [expr -$y3] [expr $y4] [expr -$y3] [expr -$y4] [expr $y3] [expr -$y4] [expr $y3] [expr $y4]
patch quad 1 4 20 [expr $y3] [expr $y2] [expr $y3] [expr -$y2] [expr $y1] [expr -$y2] [expr $y1] [expr $y2]
}

# SB2(beam 30 0.8 15 1.5)
# set some paramaters
set ht 0.3
set tw .008
set bf .15
set tf .015
set y1 [expr $ht/2.0]
set y2 [expr $bf/2.0]
set y3 [expr $y1-$tf]
set y4 [expr $tw/2.0]
section Fiber 7 {
# create bot flange,web, top flange
patch quad 1 4 20 [expr -$y1] [expr $y2] [expr -$y1] [expr -$y2] [expr -$y3] [expr -$y2] [expr -$y3] [expr $y2]
patch quad 1 30 2 [expr -$y3] [expr $y4] [expr -$y3] [expr -$y4] [expr $y3] [expr -$y4] [expr $y3] [expr $y4]
patch quad 1 4 20 [expr $y3] [expr $y2] [expr $y3] [expr -$y2] [expr $y1] [expr -$y2] [expr $y1] [expr $y2]
}

set e 21000000000.0
set nu 0.3
set G [expr $e/(2.*(1+$nu))]

set j1 1.318e-3
set j2 5.582e-4
set j3 4.39e-4
set j4 1.251e-6
set j5 1.185e-6
set j6 4.74e-7
set j7 3.615e-7

uniaxialMaterial Elastic 21 [expr $G*$j1]
uniaxialMaterial Elastic 22 [expr $G*$j2]
uniaxialMaterial Elastic 23 [expr $G*$j3]
uniaxialMaterial Elastic 24 [expr $G*$j4]
uniaxialMaterial Elastic 25 [expr $G*$j5]
uniaxialMaterial Elastic 26 [expr $G*$j6]
uniaxialMaterial Elastic 27 [expr $G*$j7]

section Aggregator 31 21 T -section 1
section Aggregator 32 22 T -section 2
section Aggregator 33 23 T -section 3
section Aggregator 34 24 T -section 4
section Aggregator 35 25 T -section 5
section Aggregator 36 26 T -section 6
section Aggregator 37 27 T -section 7


# Define column elements
# ---------------------------------

# Geometry of column elements
set n 10
#story 1 to 4

for {set x 0} {$x<=3} {incr x 1} {

# tag ndI ndJ nsecs secID transfTag
element nonlinearBeamColumn [expr 1+9*$x] [expr 1+9*$x] [expr 1+9*($x+1)] $n 32 $ColumnTransfTag
element nonlinearBeamColumn [expr 2+9*$x] [expr 2+9*$x] [expr 2+9*($x+1)] $n 31 $ColumnTransfTag
element nonlinearBeamColumn [expr 3+9*$x] [expr 3+9*$x] [expr 3+9*($x+1)] $n 32 $ColumnTransfTag
element nonlinearBeamColumn [expr 4+9*$x] [expr 4+9*$x] [expr 4+9*($x+1)] $n 31 $ColumnTransfTag
element nonlinearBeamColumn [expr 5+9*$x] [expr 5+9*$x] [expr 5+9*($x+1)] $n 31 $ColumnTransfTag
element nonlinearBeamColumn [expr 6+9*$x] [expr 6+9*$x] [expr 6+9*($x+1)] $n 31 $ColumnTransfTag
element nonlinearBeamColumn [expr 7+9*$x] [expr 7+9*$x] [expr 7+9*($x+1)] $n 32 $ColumnTransfTag
element nonlinearBeamColumn [expr 8+9*$x] [expr 8+9*$x] [expr 8+9*($x+1)] $n 31 $ColumnTransfTag
element nonlinearBeamColumn [expr 9+9*$x] [expr 9+9*$x] [expr 9+9*($x+1)] $n 32 $ColumnTransfTag

}

#story 5 & 6

for {set x 4} {$x<=5} {incr x 1} {

# tag ndI ndJ nsecs secID transfTag
element nonlinearBeamColumn [expr 1+9*$x] [expr 1+9*$x] [expr 1+9*($x+1)] $n 32 $ColumnTransfTag
element nonlinearBeamColumn [expr 2+9*$x] [expr 2+9*$x] [expr 2+9*($x+1)] $n 32 $ColumnTransfTag
element nonlinearBeamColumn [expr 3+9*$x] [expr 3+9*$x] [expr 3+9*($x+1)] $n 32 $ColumnTransfTag
element nonlinearBeamColumn [expr 4+9*$x] [expr 4+9*$x] [expr 4+9*($x+1)] $n 32 $ColumnTransfTag
element nonlinearBeamColumn [expr 5+9*$x] [expr 5+9*$x] [expr 5+9*($x+1)] $n 32 $ColumnTransfTag
element nonlinearBeamColumn [expr 6+9*$x] [expr 6+9*$x] [expr 6+9*($x+1)] $n 32 $ColumnTransfTag
element nonlinearBeamColumn [expr 7+9*$x] [expr 7+9*$x] [expr 7+9*($x+1)] $n 32 $ColumnTransfTag
element nonlinearBeamColumn [expr 8+9*$x] [expr 8+9*$x] [expr 8+9*($x+1)] $n 32 $ColumnTransfTag
element nonlinearBeamColumn [expr 9+9*$x] [expr 9+9*$x] [expr 9+9*($x+1)] $n 32 $ColumnTransfTag

}

#story 7 to 10

for {set x 6} {$x<=9} {incr x 1} {

# tag ndI ndJ nsecs secID transfTag
element nonlinearBeamColumn [expr 1+9*$x] [expr 1+9*$x] [expr 1+9*($x+1)] $n 33 $ColumnTransfTag
element nonlinearBeamColumn [expr 2+9*$x] [expr 2+9*$x] [expr 2+9*($x+1)] $n 33 $ColumnTransfTag
element nonlinearBeamColumn [expr 3+9*$x] [expr 3+9*$x] [expr 3+9*($x+1)] $n 33 $ColumnTransfTag
element nonlinearBeamColumn [expr 4+9*$x] [expr 4+9*$x] [expr 4+9*($x+1)] $n 33 $ColumnTransfTag
element nonlinearBeamColumn [expr 5+9*$x] [expr 5+9*$x] [expr 5+9*($x+1)] $n 33 $ColumnTransfTag
element nonlinearBeamColumn [expr 6+9*$x] [expr 6+9*$x] [expr 6+9*($x+1)] $n 33 $ColumnTransfTag
element nonlinearBeamColumn [expr 7+9*$x] [expr 7+9*$x] [expr 7+9*($x+1)] $n 33 $ColumnTransfTag
element nonlinearBeamColumn [expr 8+9*$x] [expr 8+9*$x] [expr 8+9*($x+1)] $n 33 $ColumnTransfTag
element nonlinearBeamColumn [expr 9+9*$x] [expr 9+9*$x] [expr 9+9*($x+1)] $n 33 $ColumnTransfTag

}

#Beams
#1-8

for {set x 0} {$x<=7} {incr x 1} {

# tag ndI ndJ nsecs secID transfTag
element nonlinearBeamColumn [expr 91+6*$x] [expr 10+9*$x] [expr 11+9*$x] $n 36 $BeamTransfTag
element nonlinearBeamColumn [expr 92+6*$x] [expr 11+9*$x] [expr 12+9*$x] $n 36 $BeamTransfTag
element nonlinearBeamColumn [expr 93+6*$x] [expr 13+9*$x] [expr 14+9*$x] $n 36 $BeamTransfTag
element nonlinearBeamColumn [expr 94+6*$x] [expr 14+9*$x] [expr 15+9*$x] $n 36 $BeamTransfTag
element nonlinearBeamColumn [expr 95+6*$x] [expr 16+9*$x] [expr 17+9*$x] $n 36 $BeamTransfTag
element nonlinearBeamColumn [expr 96+6*$x] [expr 17+9*$x] [expr 18+9*$x] $n 36 $BeamTransfTag

}

#9,10

for {set x 8} {$x<=9} {incr x 1} {

# tag ndI ndJ nsecs secID transfTag
element nonlinearBeamColumn [expr 91+6*$x] [expr 10+9*$x] [expr 11+9*$x] $n 37 $BeamTransfTag
element nonlinearBeamColumn [expr 92+6*$x] [expr 11+9*$x] [expr 12+9*$x] $n 37 $BeamTransfTag
element nonlinearBeamColumn [expr 93+6*$x] [expr 13+9*$x] [expr 14+9*$x] $n 37 $BeamTransfTag
element nonlinearBeamColumn [expr 94+6*$x] [expr 14+9*$x] [expr 15+9*$x] $n 37 $BeamTransfTag
element nonlinearBeamColumn [expr 95+6*$x] [expr 16+9*$x] [expr 17+9*$x] $n 37 $BeamTransfTag
element nonlinearBeamColumn [expr 96+6*$x] [expr 17+9*$x] [expr 18+9*$x] $n 37 $BeamTransfTag
}

#Girders
#1-8
for {set x 0} {$x<=7} {incr x 1} {

# tag ndI ndJ nsecs secID transfTag
element nonlinearBeamColumn [expr 151+6*$x] [expr 10+9*$x] [expr 13+9*$x] $n 36 $GirderTransfTag
element nonlinearBeamColumn [expr 152+6*$x] [expr 11+9*$x] [expr 14+9*$x] $n 34 $GirderTransfTag
element nonlinearBeamColumn [expr 153+6*$x] [expr 12+9*$x] [expr 15+9*$x] $n 36 $GirderTransfTag
element nonlinearBeamColumn [expr 154+6*$x] [expr 13+9*$x] [expr 16+9*$x] $n 36 $GirderTransfTag
element nonlinearBeamColumn [expr 155+6*$x] [expr 14+9*$x] [expr 17+9*$x] $n 34 $GirderTransfTag
element nonlinearBeamColumn [expr 156+6*$x] [expr 15+9*$x] [expr 18+9*$x] $n 36 $GirderTransfTag

}

# 9,10

for {set x 8} {$x<=9} {incr x 1} {

# tag ndI ndJ nsecs secID transfTag
element nonlinearBeamColumn [expr 151+6*$x] [expr 10+9*$x] [expr 13+9*$x] $n 37 $GirderTransfTag
element nonlinearBeamColumn [expr 152+6*$x] [expr 11+9*$x] [expr 14+9*$x] $n 35 $GirderTransfTag
element nonlinearBeamColumn [expr 153+6*$x] [expr 12+9*$x] [expr 15+9*$x] $n 37 $GirderTransfTag
element nonlinearBeamColumn [expr 154+6*$x] [expr 13+9*$x] [expr 16+9*$x] $n 37 $GirderTransfTag
element nonlinearBeamColumn [expr 155+6*$x] [expr 14+9*$x] [expr 17+9*$x] $n 35 $GirderTransfTag
element nonlinearBeamColumn [expr 156+6*$x] [expr 15+9*$x] [expr 18+9*$x] $n 37 $GirderTransfTag

}

#recorder Node -file eigen1.out -node 13 -dof 1 2 6 "eigen1"
#recorder Node -file eigen2.out -node 13 -dof 1 2 6 "eigen2"
#recorder Node -file eigen3.out -node 13 -dof 1 2 6 "eigen3"

set eig [eigen 3]
set pi [expr 4*atan(1)]
puts "First Mode Period is:[expr 2*$pi/pow([lindex $eig 0],0.5)]"
puts "Second Mode Period is:[expr 2*$pi/pow([lindex $eig 1],0.5)]"
puts "Third Mode Period is:[expr 2*$pi/pow([lindex $eig 2],0.5)]"

[fmk]

the models are different! how can one help if you only show the OpenSees model?? this is really your problem. Like everyone who runs into this issue, i suggest setting the materials elastic to see if that solves it, if it does it is the materials definitions; if not it is the model geometry or element definitions.