Dear all
i have modeled SAC9 steel model. the modal periods are the same as experimental results but during the static analysis, the compatibility error happens and the model does not converge. please send me some advice about the error.
#Units: N, m
wipe;
model BasicBuilder -ndm 2 -ndf 3;
# define GEOMETRY ---------------------------------------------------
set Y1 3.65;
set Y2 5.49;
set Y3 3.96;
set X 9.15;
set i1m 96484.96408
set e1m 48330.03645
set i2m 102088.4466
set e2m 50956.66887
set im 98060.94353
set em 49030.47176
set irm 108042.1467
set erm 53933.51894
# Base
#Definition of Nodes Id and cordinate with mass at the connection between column to the foundation
# ID x, y massx massy massz
node 01 0 0 -mass 0 0 0
node 02 [expr 1*$X] 0 -mass 0 0 0
node 03 [expr 2*$X] 0 -mass 0 0 0
node 04 [expr 3*$X] 0 -mass 0 0 0
node 05 [expr 4*$X] 0 -mass 0 0 0
node 06 [expr 5*$X] 0 -mass 0 0 0
# Story 1
node 11 0 $Y1 -mass $e1m 0 0
node 12 [expr 1*$X] $Y1 -mass $i1m 0 0
node 13 [expr 2*$X] $Y1 -mass $i1m 0 0
node 14 [expr 3*$X] $Y1 -mass $i1m 0 0
node 15 [expr 4*$X] $Y1 -mass $i1m 0 0
node 16 [expr 5*$X] $Y1 -mass $e1m 0 0
node 17 [expr 5*$X] $Y1 -mass 0 0 0
#Using Node 17 to model hinge restraint
# Story 2
node 21 0 [expr $Y1+$Y2] -mass $e2m 0 0
node 22 [expr 1*$X] [expr $Y1+$Y2] -mass $i2m 0 0
node 23 [expr 2*$X] [expr $Y1+$Y2] -mass $i2m 0 0
node 24 [expr 3*$X] [expr $Y1+$Y2] -mass $i2m 0 0
node 25 [expr 4*$X] [expr $Y1+$Y2] -mass $i2m 0 0
node 26 [expr 5*$X] [expr $Y1+$Y2] -mass $e2m 0 0
node 27 [expr 5*$X] [expr $Y1+$Y2] -mass 0 0 0
#Splice column between story 2 nd 3
node 931 0 [expr $Y1+$Y2+0.46*$Y3] -mass 0 0 0
node 932 [expr 1*$X] [expr $Y1+$Y2+0.46*$Y3] -mass 0 0 0
node 933 [expr 2*$X] [expr $Y1+$Y2+0.46*$Y3] -mass 0 0 0
node 934 [expr 3*$X] [expr $Y1+$Y2+0.46*$Y3] -mass 0 0 0
node 935 [expr 4*$X] [expr $Y1+$Y2+0.46*$Y3] -mass 0 0 0
node 936 [expr 5*$X] [expr $Y1+$Y2+0.46*$Y3] -mass 0 0 0
# Story 3
node 31 0 [expr $Y1+$Y2+$Y3] -mass $em 0 0
node 32 [expr 1*$X] [expr $Y1+$Y2+$Y3] -mass $im 0 0
node 33 [expr 2*$X] [expr $Y1+$Y2+$Y3] -mass $im 0 0
node 34 [expr 3*$X] [expr $Y1+$Y2+$Y3] -mass $im 0 0
node 35 [expr 4*$X] [expr $Y1+$Y2+$Y3] -mass $im 0 0
node 36 [expr 5*$X] [expr $Y1+$Y2+$Y3] -mass $em 0 0
node 37 [expr 5*$X] [expr $Y1+$Y2+$Y3] -mass 0 0 0
# Story 4
node 41 0 [expr $Y1+$Y2+2*$Y3] -mass $em 0 0
node 42 [expr 1*$X] [expr $Y1+$Y2+2*$Y3] -mass $im 0 0
node 43 [expr 2*$X] [expr $Y1+$Y2+2*$Y3] -mass $im 0 0
node 44 [expr 3*$X] [expr $Y1+$Y2+2*$Y3] -mass $im 0 0
node 45 [expr 4*$X] [expr $Y1+$Y2+2*$Y3] -mass $im 0 0
node 46 [expr 5*$X] [expr $Y1+$Y2+2*$Y3] -mass $em 0 0
node 47 [expr 5*$X] [expr $Y1+$Y2+2*$Y3] -mass 0 0 0
#Splice column between story 4 nd 5
node 951 0 [expr $Y1+$Y2+2.46*$Y3] -mass 0 0 0
node 952 [expr 1*$X] [expr $Y1+$Y2+2.46*$Y3] -mass 0 0 0
node 953 [expr 2*$X] [expr $Y1+$Y2+2.46*$Y3] -mass 0 0 0
node 954 [expr 3*$X] [expr $Y1+$Y2+2.46*$Y3] -mass 0 0 0
node 955 [expr 4*$X] [expr $Y1+$Y2+2.46*$Y3] -mass 0 0 0
node 956 [expr 5*$X] [expr $Y1+$Y2+2.46*$Y3] -mass 0 0 0
# Story 5
node 51 0 [expr $Y1+$Y2+3*$Y3] -mass $em 0 0
node 52 [expr 1*$X] [expr $Y1+$Y2+3*$Y3] -mass $im 0 0
node 53 [expr 2*$X] [expr $Y1+$Y2+3*$Y3] -mass $im 0 0
node 54 [expr 3*$X] [expr $Y1+$Y2+3*$Y3] -mass $im 0 0
node 55 [expr 4*$X] [expr $Y1+$Y2+3*$Y3] -mass $im 0 0
node 56 [expr 5*$X] [expr $Y1+$Y2+3*$Y3] -mass $em 0 0
node 57 [expr 5*$X] [expr $Y1+$Y2+3*$Y3] -mass 0 0 0
# Story 6
node 61 0 [expr $Y1+$Y2+4*$Y3] -mass $em 0 0
node 62 [expr 1*$X] [expr $Y1+$Y2+4*$Y3] -mass $im 0 0
node 63 [expr 2*$X] [expr $Y1+$Y2+4*$Y3] -mass $im 0 0
node 64 [expr 3*$X] [expr $Y1+$Y2+4*$Y3] -mass $im 0 0
node 65 [expr 4*$X] [expr $Y1+$Y2+4*$Y3] -mass $im 0 0
node 66 [expr 5*$X] [expr $Y1+$Y2+4*$Y3] -mass $em 0 0
node 67 [expr 5*$X] [expr $Y1+$Y2+4*$Y3] -mass 0 0 0
#Splice column between story 7 nd 8
node 971 0 [expr $Y1+$Y2+4.46*$Y3] -mass 0 0 0
node 972 [expr 1*$X] [expr $Y1+$Y2+4.46*$Y3] -mass 0 0 0
node 973 [expr 2*$X] [expr $Y1+$Y2+4.46*$Y3] -mass 0 0 0
node 974 [expr 3*$X] [expr $Y1+$Y2+4.46*$Y3] -mass 0 0 0
node 975 [expr 4*$X] [expr $Y1+$Y2+4.46*$Y3] -mass 0 0 0
node 976 [expr 5*$X] [expr $Y1+$Y2+4.46*$Y3] -mass 0 0 0
# Story 7
node 71 0 [expr $Y1+$Y2+5*$Y3] -mass $em 0 0
node 72 [expr 1*$X] [expr $Y1+$Y2+5*$Y3] -mass $im 0 0
node 73 [expr 2*$X] [expr $Y1+$Y2+5*$Y3] -mass $im 0 0
node 74 [expr 3*$X] [expr $Y1+$Y2+5*$Y3] -mass $im 0 0
node 75 [expr 4*$X] [expr $Y1+$Y2+5*$Y3] -mass $im 0 0
node 76 [expr 5*$X] [expr $Y1+$Y2+5*$Y3] -mass $em 0 0
node 77 [expr 5*$X] [expr $Y1+$Y2+5*$Y3] -mass 0 0 0
# Story 8
node 81 0 [expr $Y1+$Y2+6*$Y3] -mass $em 0 0
node 82 [expr 1*$X] [expr $Y1+$Y2+6*$Y3] -mass $im 0 0
node 83 [expr 2*$X] [expr $Y1+$Y2+6*$Y3] -mass $im 0 0
node 84 [expr 3*$X] [expr $Y1+$Y2+6*$Y3] -mass $im 0 0
node 85 [expr 4*$X] [expr $Y1+$Y2+6*$Y3] -mass $im 0 0
node 86 [expr 5*$X] [expr $Y1+$Y2+6*$Y3] -mass $em 0 0
node 87 [expr 5*$X] [expr $Y1+$Y2+6*$Y3] -mass 0 0 0
#Splice column between story 8 nd 9
node 991 0 [expr $Y1+$Y2+6.46*$Y3] -mass 0 0 0
node 992 [expr 1*$X] [expr $Y1+$Y2+6.46*$Y3] -mass 0 0 0
node 993 [expr 2*$X] [expr $Y1+$Y2+6.46*$Y3] -mass 0 0 0
node 994 [expr 3*$X] [expr $Y1+$Y2+6.46*$Y3] -mass 0 0 0
node 995 [expr 4*$X] [expr $Y1+$Y2+6.46*$Y3] -mass 0 0 0
node 996 [expr 5*$X] [expr $Y1+$Y2+6.46*$Y3] -mass 0 0 0
# Story 9
node 91 0 [expr $Y1+$Y2+7*$Y3] -mass $em 0 0
node 92 [expr 1*$X] [expr $Y1+$Y2+7*$Y3] -mass $im 0 0
node 93 [expr 2*$X] [expr $Y1+$Y2+7*$Y3] -mass $im 0 0
node 94 [expr 3*$X] [expr $Y1+$Y2+7*$Y3] -mass $im 0 0
node 95 [expr 4*$X] [expr $Y1+$Y2+7*$Y3] -mass $im 0 0
node 96 [expr 5*$X] [expr $Y1+$Y2+7*$Y3] -mass $em 0 0
node 97 [expr 5*$X] [expr $Y1+$Y2+7*$Y3] -mass 0 0 0
# Story 10
node 101 0 [expr $Y1+$Y2+8*$Y3] -mass $erm 0 0
node 102 [expr 1*$X] [expr $Y1+$Y2+8*$Y3] -mass $irm 0 0
node 103 [expr 2*$X] [expr $Y1+$Y2+8*$Y3] -mass $irm 0 0
node 104 [expr 3*$X] [expr $Y1+$Y2+8*$Y3] -mass $irm 0 0
node 105 [expr 4*$X] [expr $Y1+$Y2+8*$Y3] -mass $irm 0 0
node 106 [expr 5*$X] [expr $Y1+$Y2+8*$Y3] -mass $erm 0 0
node 107 [expr 5*$X] [expr $Y1+$Y2+8*$Y3] -mass 0 0 0
# M1, M2, Mf and Mr are the total mass of each storey.
set M1 [expr 2*$e1m+4*$i1m];
set M2 [expr 2*$e2m+4*$i2m];
set Mf [expr 2*$em+4*$im];
set Mr [expr 2*$erm+4*$irm];
#Display Mass of Different Stories
puts " M1: $M1 "
puts " M2: $M2 "
puts " Mf: $Mf "
puts " Mr: $Mr "
# Single point constraints -- Boundary Conditions
# All connecton are pinned at the base floor
fix 01 1 1 0
fix 02 1 1 0
fix 03 1 1 0
fix 04 1 1 0
fix 05 1 1 0
fix 06 1 1 0
# Boundary condition of node 11 and 16 at floor 1 (both sides of building)
fix 11 1 0 0
fix 16 1 0 0
#impose the traslational displacements in x and y directions node 3 to be the same as those of
#node2
equalDOF 16 17 1 2
equalDOF 26 27 1 2
equalDOF 36 37 1 2
equalDOF 46 47 1 2
equalDOF 56 57 1 2
equalDOF 66 67 1 2
equalDOF 76 77 1 2
equalDOF 86 87 1 2
equalDOF 96 97 1 2
equalDOF 106 107 1 2
# tag E fy Hiso Hkin
#cloumn
uniaxialMaterial Hardening 1 1.998e11 3.45e8 0 1000
#beam
uniaxialMaterial Hardening 2 1.998e11 2.48e8 0 1000
# define SECTION ---------------------------------------------------
source wsection.tcl
#Beam Sections
# Section "w36x160": secTag
#Wsection $GirdSecTagFiber $matIDhard $d $bf $tf $tw $nfdw $nftw $nfbf $nftf
Wsection 1 2 0.9144 0.3048 0.025908 0.01651 16 4 16 4
# Section "w36x135": secTag
#Wsection $GirdSecTagFiber $matIDhard $d $bf $tf $tw $nfdw $nftw $nfbf $nftf
Wsection 2 2 0.90424 0.3048 0.020066 0.01524 16 4 16 4
# Section "w30x99": secTag
#Wsection $GirdSecTagFiber $matIDhard $d $bf $tf $tw $nfdw $nftw $nfbf $nftf
Wsection 3 2 0.75438 0.2667 0.017018 0.013208 16 4 16 4
# Section "w27x84": secTag
#Wsection $GirdSecTagFiber $matIDhard $d $bf $tf $tw $nfdw $nftw $nfbf $nftf
Wsection 4 2 0.67818 0.254 0.016256 0.011684 16 4 16 4
# Section "w24x68": secTag
#Wsection $GirdSecTagFiber $matIDhard $d $bf $tf $tw $nfdw $nftw $nfbf $nftf
Wsection 5 2 0.60198 0.227838 0.014859 0.010541 16 4 16 4
#Column Sections
# Section "w14x500": secTag
#Wsection $GirdSecTagFiber $matIDhard $d $bf $tf $tw $nfdw $nftw $nfbf $nftf
Wsection 6 1 0.49784 0.4318 0.0889 0.055626 16 4 16 4
# Section "w14x455": secTag
#Wsection $GirdSecTagFiber $matIDhard $d $bf $tf $tw $nfdw $nftw $nfbf $nftf
Wsection 7 1 0.4826 0.42672 0.081534 0.051308 16 4 16 4
# Section "w14x370": secTag
#Wsection $GirdSecTagFiber $matIDhard $d $bf $tf $tw $nfdw $nftw $nfbf $nftf
Wsection 8 1 0.45466 0.4191 0.067564 0.042164 16 4 16 4
# Section "w14x283": secTag
#Wsection $GirdSecTagFiber $matIDhard $d $bf $tf $tw $nfdw $nftw $nfbf $nftf
Wsection 9 1 0.42418 0.40894 0.052578 0.032766 16 4 16 4
# Section "w14x257": secTag
#Wsection $GirdSecTagFiber $matIDhard $d $bf $tf $tw $nfdw $nftw $nfbf $nftf
Wsection 10 1 0.41656 0.4064 0.048006 0.029972 16 4 16 4
# define ELEMENTS ---------------------------------------------------
set IDcolTrans 1;
set IDbeamTrans 2;
# options, Linear PDelta Corotational
geomTransf Corotational $IDcolTrans ;
geomTransf Corotational $IDbeamTrans ;
#Column Elements
# Eleme nt "cloumns 1": eleTag NodeI NodeJ NIP secTag geoTranTag
element nonlinearBeamColumn 110 1 11 5 6 1
element nonlinearBeamColumn 120 2 12 5 6 1
element nonlinearBeamColumn 130 3 13 5 6 1
element nonlinearBeamColumn 144 4 14 5 6 1
element nonlinearBeamColumn 150 5 15 5 6 1
element nonlinearBeamColumn 160 6 16 5 6 1
# Eleme nt "cloumns 2": eleTag NodeI NodeJ NIP secTag geoTranTag
element nonlinearBeamColumn 210 11 21 5 6 1
element nonlinearBeamColumn 220 12 22 5 6 1
element nonlinearBeamColumn 230 13 23 5 6 1
element nonlinearBeamColumn 240 14 24 5 6 1
element nonlinearBeamColumn 250 15 25 5 6 1
element nonlinearBeamColumn 260 16 26 5 6 1
# Eleme nt "cloumns 3": eleTag NodeI NodeJ NIP secTag geoTranTag
element nonlinearBeamColumn 311 21 931 5 6 1
element nonlinearBeamColumn 321 22 932 5 6 1
element nonlinearBeamColumn 331 23 933 5 6 1
element nonlinearBeamColumn 341 24 934 5 6 1
element nonlinearBeamColumn 351 25 935 5 6 1
element nonlinearBeamColumn 361 26 936 5 6 1
element nonlinearBeamColumn 312 931 31 5 7 1
element nonlinearBeamColumn 322 932 32 5 7 1
element nonlinearBeamColumn 332 933 33 5 7 1
element nonlinearBeamColumn 342 934 34 5 7 1
element nonlinearBeamColumn 352 935 35 5 7 1
element nonlinearBeamColumn 362 936 36 5 7 1
# Eleme nt "cloumns 4": eleTag NodeI NodeJ NIP secTag geoTranTag
element nonlinearBeamColumn 410 31 41 5 7 1
element nonlinearBeamColumn 420 32 42 5 7 1
element nonlinearBeamColumn 430 33 43 5 7 1
element nonlinearBeamColumn 440 34 44 5 7 1
element nonlinearBeamColumn 450 35 45 5 7 1
element nonlinearBeamColumn 460 36 46 5 7 1
# Eleme nt "cloumns 5": eleTag NodeI NodeJ NIP secTag geoTranTag
element nonlinearBeamColumn 511 41 951 5 7 1
element nonlinearBeamColumn 521 42 952 5 7 1
element nonlinearBeamColumn 531 43 953 5 7 1
element nonlinearBeamColumn 541 44 954 5 7 1
element nonlinearBeamColumn 551 45 955 5 7 1
element nonlinearBeamColumn 561 46 956 5 7 1
element nonlinearBeamColumn 512 951 51 5 8 1
element nonlinearBeamColumn 522 952 52 5 8 1
element nonlinearBeamColumn 532 953 53 5 8 1
element nonlinearBeamColumn 542 954 54 5 8 1
element nonlinearBeamColumn 552 955 55 5 8 1
element nonlinearBeamColumn 562 956 56 5 8 1
# Eleme nt "cloumns 6": eleTag NodeI NodeJ NIP secTag geoTranTag
element nonlinearBeamColumn 610 51 61 5 8 1
element nonlinearBeamColumn 620 52 62 5 8 1
element nonlinearBeamColumn 630 53 63 5 8 1
element nonlinearBeamColumn 640 54 64 5 8 1
element nonlinearBeamColumn 650 55 65 5 8 1
element nonlinearBeamColumn 660 56 66 5 8 1
# Eleme nt "cloumns 7": eleTag NodeI NodeJ NIP secTag geoTranTag
element nonlinearBeamColumn 711 61 971 5 8 1
element nonlinearBeamColumn 721 62 972 5 8 1
element nonlinearBeamColumn 731 63 973 5 8 1
element nonlinearBeamColumn 741 64 974 5 8 1
element nonlinearBeamColumn 751 65 975 5 8 1
element nonlinearBeamColumn 761 66 976 5 8 1
element nonlinearBeamColumn 712 971 71 5 9 1
element nonlinearBeamColumn 722 972 72 5 9 1
element nonlinearBeamColumn 732 973 73 5 9 1
element nonlinearBeamColumn 742 974 74 5 9 1
element nonlinearBeamColumn 752 975 75 5 9 1
element nonlinearBeamColumn 762 976 76 5 9 1
# Eleme nt "cloumns 8": eleTag NodeI NodeJ NIP secTag geoTranTag
element nonlinearBeamColumn 810 71 81 5 9 1
element nonlinearBeamColumn 820 72 82 5 9 1
element nonlinearBeamColumn 830 73 83 5 9 1
element nonlinearBeamColumn 840 74 84 5 9 1
element nonlinearBeamColumn 850 75 85 5 9 1
element nonlinearBeamColumn 860 76 86 5 9 1
# Eleme nt "cloumns 9": eleTag NodeI NodeJ NIP secTag geoTranTag
element nonlinearBeamColumn 911 81 991 5 9 1
element nonlinearBeamColumn 921 82 992 5 9 1
element nonlinearBeamColumn 931 83 993 5 9 1
element nonlinearBeamColumn 941 84 994 5 9 1
element nonlinearBeamColumn 951 85 995 5 9 1
element nonlinearBeamColumn 961 86 996 5 9 1
element nonlinearBeamColumn 912 991 91 5 10 1
element nonlinearBeamColumn 922 992 92 5 10 1
element nonlinearBeamColumn 932 993 93 5 10 1
element nonlinearBeamColumn 942 994 94 5 10 1
element nonlinearBeamColumn 952 995 95 5 10 1
element nonlinearBeamColumn 962 996 96 5 10 1
# Eleme nt "cloumns 10": eleTag NodeI NodeJ NIP secTag geoTranTag
element nonlinearBeamColumn 1010 91 101 5 10 1
element nonlinearBeamColumn 1020 92 102 5 10 1
element nonlinearBeamColumn 1030 93 103 5 10 1
element nonlinearBeamColumn 1040 94 104 5 10 1
element nonlinearBeamColumn 1050 95 105 5 10 1
element nonlinearBeamColumn 1060 96 106 5 10 1
#Beams Elements
#Element beams 1: eleTag NodeI NodeJ NIP secTag geoTranTag
element nonlinearBeamColumn 11 11 12 5 1 2
element nonlinearBeamColumn 12 12 13 5 1 2
element nonlinearBeamColumn 13 13 14 5 1 2
element nonlinearBeamColumn 14 14 15 5 1 2
element nonlinearBeamColumn 15 15 17 5 1 2
#Element beams 2: eleTag NodeI NodeJ NIP secTag geoTranTag
element nonlinearBeamColumn 21 21 22 5 1 2
element nonlinearBeamColumn 22 22 23 5 1 2
element nonlinearBeamColumn 23 23 24 5 1 2
element nonlinearBeamColumn 24 24 25 5 1 2
element nonlinearBeamColumn 25 25 27 5 1 2
#Element beams 3: eleTag NodeI NodeJ NIP secTag geoTranTag
element nonlinearBeamColumn 31 31 32 5 1 2
element nonlinearBeamColumn 32 32 33 5 1 2
element nonlinearBeamColumn 33 33 34 5 1 2
element nonlinearBeamColumn 34 34 35 5 1 2
element nonlinearBeamColumn 35 35 37 5 1 2
#Element beams 4: eleTag NodeI NodeJ NIP secTag geoTranTag
element nonlinearBeamColumn 41 41 42 5 2 2
element nonlinearBeamColumn 42 42 43 5 2 2
element nonlinearBeamColumn 43 43 44 5 2 2
element nonlinearBeamColumn 44 44 45 5 2 2
element nonlinearBeamColumn 45 45 47 5 2 2
#Element beams 5: eleTag NodeI NodeJ NIP secTag geoTranTag
element nonlinearBeamColumn 51 51 52 5 2 2
element nonlinearBeamColumn 52 52 53 5 2 2
element nonlinearBeamColumn 53 53 54 5 2 2
element nonlinearBeamColumn 54 54 55 5 2 2
element nonlinearBeamColumn 55 55 57 5 2 2
#Element beams 6: eleTag NodeI NodeJ NIP secTag geoTranTag
element nonlinearBeamColumn 61 61 62 5 2 2
element nonlinearBeamColumn 62 62 63 5 2 2
element nonlinearBeamColumn 63 63 64 5 2 2
element nonlinearBeamColumn 64 64 65 5 2 2
element nonlinearBeamColumn 65 65 67 5 2 2
#Element beams 7: eleTag NodeI NodeJ NIP secTag geoTranTag
element nonlinearBeamColumn 71 71 72 5 2 2
element nonlinearBeamColumn 72 72 73 5 2 2
element nonlinearBeamColumn 73 73 74 5 2 2
element nonlinearBeamColumn 74 74 75 5 2 2
element nonlinearBeamColumn 75 75 77 5 2 2
#Element beams 8: eleTag NodeI NodeJ NIP secTag geoTranTag
element nonlinearBeamColumn 81 81 82 5 3 2
element nonlinearBeamColumn 82 82 83 5 3 2
element nonlinearBeamColumn 83 83 84 5 3 2
element nonlinearBeamColumn 84 84 85 5 3 2
element nonlinearBeamColumn 85 85 87 5 3 2
#Element beams 9: eleTag NodeI NodeJ NIP secTag geoTranTag
element nonlinearBeamColumn 91 91 92 5 4 2
element nonlinearBeamColumn 92 92 93 5 4 2
element nonlinearBeamColumn 93 93 94 5 4 2
element nonlinearBeamColumn 94 94 95 5 4 2
element nonlinearBeamColumn 95 95 97 5 4 2
#Element beams 10: eleTag NodeI NodeJ NIP secTag geoTranTag
element nonlinearBeamColumn 101 101 102 5 5 2
element nonlinearBeamColumn 102 102 103 5 5 2
element nonlinearBeamColumn 103 103 104 5 5 2
element nonlinearBeamColumn 104 104 105 5 5 2
element nonlinearBeamColumn 105 105 107 5 5 2
# ---------------------------------------------------------------------
#Calculation of Final Periods and Frequencies:
# ---------------------------------------------------------------------
source Eigen.txt
# ---------------------------------------------------------------------
# Constant gravity load
#----------------------------------------------------------------------
set W1 103482
set W2 109414
set W3 105134
set W4 115797
pattern Plain 1 Constant {eleLoad -ele 11 -type -beamUniform -$W1
eleLoad -ele 12 -type -beamUniform -$W1
eleLoad -ele 13 -type -beamUniform -$W1
eleLoad -ele 14 -type -beamUniform -$W1
eleLoad -ele 15 -type -beamUniform -$W1
eleLoad -ele 21 -type -beamUniform -$W2
eleLoad -ele 22 -type -beamUniform -$W2
eleLoad -ele 23 -type -beamUniform -$W2
eleLoad -ele 24 -type -beamUniform -$W2
eleLoad -ele 25 -type -beamUniform -$W2
eleLoad -ele 31 -type -beamUniform -$W3
eleLoad -ele 32 -type -beamUniform -$W3
eleLoad -ele 33 -type -beamUniform -$W3
eleLoad -ele 34 -type -beamUniform -$W3
eleLoad -ele 35 -type -beamUniform -$W3
eleLoad -ele 41 -type -beamUniform -$W3
eleLoad -ele 42 -type -beamUniform -$W3
eleLoad -ele 43 -type -beamUniform -$W3
eleLoad -ele 44 -type -beamUniform -$W3
eleLoad -ele 45 -type -beamUniform -$W3
eleLoad -ele 51 -type -beamUniform -$W3
eleLoad -ele 52 -type -beamUniform -$W3
eleLoad -ele 53 -type -beamUniform -$W3
eleLoad -ele 54 -type -beamUniform -$W3
eleLoad -ele 55 -type -beamUniform -$W3
eleLoad -ele 61 -type -beamUniform -$W3
eleLoad -ele 62 -type -beamUniform -$W3
eleLoad -ele 63 -type -beamUniform -$W3
eleLoad -ele 64 -type -beamUniform -$W3
eleLoad -ele 65 -type -beamUniform -$W3
eleLoad -ele 71 -type -beamUniform -$W3
eleLoad -ele 72 -type -beamUniform -$W3
eleLoad -ele 73 -type -beamUniform -$W3
eleLoad -ele 74 -type -beamUniform -$W3
eleLoad -ele 75 -type -beamUniform -$W3
eleLoad -ele 81 -type -beamUniform -$W3
eleLoad -ele 82 -type -beamUniform -$W3
eleLoad -ele 83 -type -beamUniform -$W3
eleLoad -ele 84 -type -beamUniform -$W3
eleLoad -ele 85 -type -beamUniform -$W3
eleLoad -ele 91 -type -beamUniform -$W3
eleLoad -ele 92 -type -beamUniform -$W3
eleLoad -ele 93 -type -beamUniform -$W3
eleLoad -ele 94 -type -beamUniform -$W3
eleLoad -ele 95 -type -beamUniform -$W3
eleLoad -ele 101 -type -beamUniform -$W4
eleLoad -ele 102 -type -beamUniform -$W4
eleLoad -ele 103 -type -beamUniform -$W4
eleLoad -ele 104 -type -beamUniform -$W4
eleLoad -ele 105 -type -beamUniform -$W4
}
# ---------------------------------------------------------------------
# Analysis Generation for Gravity Analysis:
# ---------------------------------------------------------------------
test NormDispIncr 1.0e-8 10 0; #convergence test: 1e-8 and a max number of iterations of 10
algorithm Newton
integrator LoadControl 0.1 1 0.1 0.1; #integration scheme: the LoadControl scheme using steps of 0.1
system BandGeneral; #system of equation: a SPD using a profile storage scheme
numberer RCM; #DOF numberer: the reverse Cuthill-McKee algorithm
constraints Plain; #constraint handler: the transformation method
analysis Static; #the analysis object
#
analyze 10; # perform the gravity load analysis,by 10 steps
loadConst -time 0.0
puts "END of Gravity Load Analyze"
Problem in SAC 9 Steel model- SI unit
Moderators: silvia, selimgunay, Moderators
-
n_siahpolo
- Posts: 12
- Joined: Tue Nov 09, 2010 10:53 am
- Location: ACECR Institute for higher education-Ahwaz