Anybody can help me? wipe
#-----------------------------------------------------------------------------------------
# 1. CREATE SOIL NODES
#-----------------------------------------------------------------------------------------
model BasicBuilder -ndm 2 -ndf 2
# define soil nodes
# Node tag xCrd yCrd
node 1 +0.000000E+000 +0.000000E+000
node 2 +1.100000E+000 +0.000000E+000
node 3 +2.200000E+000 +0.000000E+000
node 4 +3.300000E+000 +0.000000E+000
....
node 1116 +5.690000E+001 +1.894000E+001
node 1117 +5.800000E+001 +1.894000E+001
node 1118 +5.910000E+001 +1.894000E+001
node 1119 +6.020000E+001 +1.894000E+001
node 1120 +6.130000E+001 +1.894000E+001
puts "Finshed creating all soil nodes..."
#-----------------------------------------------------------------------------------------
# 2. Define SOIL materials
#-----------------------------------------------------------------------------------------
# Material "Nevada dry sand (Dr=74%)": matTag nd rho Gr Br phi gammamax pr d phipt contrac dilat1 dilat2 liqf1 liqf2 liqf3
# noYS e cs1 cs2 cs3 pa
nDMaterial PressureDependMultiYield 1 2 +1.692 +5.300000E+004 +1.150000E+005 +3.500000E+001 +1.000000E-001 +5.400000E+001 +5.000000E-001 +2.700000E+001 +5.000000E-002 +6.000000E-001 +3.000000E+000 +0.000000E+000 +0.000000E+000 +0.000000E+000 \
11 +5.660000E-001 +9.000000E-001 +2.000000E-002 +7.000000E-001 +1.010000E+002
updateMaterialStage -material 1 -stage 0
# body force in x-direction
set xWgt1 0.00
# body force in y-direction
set yWgt1 [expr -9.81*1.692]
puts "Finished creating all soil materials..."
#-----------------------------------------------------------------------------------------
# 3. define quadElement for soil
#-----------------------------------------------------------------------------------------
element quad 1 1 2 12 11 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 2 2 3 13 12 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 3 3 4 14 13 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
....
element quad 1017 1106 1107 1117 1116 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 1018 1107 1108 1118 1117 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 1019 1108 1109 1119 1118 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
element quad 1020 1109 1110 1120 1119 1.00E+00 PlaneStrain 1 0 1.692 $xWgt1 $yWgt1
puts "Finished creating all soil elements..."
#-----------------------------------------------------------------------------------------
# 4. define fixities for Bottom boundary
#-----------------------------------------------------------------------------------------
# SPC tag Dx Dy Rz
fix 1 1 1 1
fix 2 1 1 1
fix 3 1 1 1
fix 4 1 1 1
fix 5 1 1 1
fix 6 1 1 1
fix 7 1 1 1
fix 8 1 1 1
fix 9 1 1 1
fix 10 1 1 1
fix 281 1 1 1
fix 282 1 1 1
fix 283 1 1 1
fix 284 1 1 1
fix 285 1 1 1
fix 286 1 1 1
fix 365 1 1 1
fix 366 1 1 1
fix 367 1 1 1
fix 368 1 1 1
fix 369 1 1 1
fix 370 1 1 1
fix 371 1 1 1
fix 372 1 1 1
fix 373 1 1 1
fix 374 1 1 1
fix 375 1 1 1
fix 376 1 1 1
fix 377 1 1 1
fix 378 1 1 1
fix 771 1 1 1
fix 772 1 1 1
fix 773 1 1 1
fix 774 1 1 1
fix 775 1 1 1
fix 776 1 1 1
fix 855 1 1 1
fix 856 1 1 1
fix 857 1 1 1
fix 858 1 1 1
fix 859 1 1 1
fix 860 1 1 1
fix 861 1 1 1
fix 862 1 1 1
fix 863 1 1 1
puts "Finished creating all fixities..."
#-----------------------------------------------------------------------------------------
# 5. define equal DOF (soil to soil)
#-----------------------------------------------------------------------------------------
# Equal DOF: boundary: mNodeTag sNodeTag dof
equalDOF 1 863 1 2
equalDOF 11 872 1 2
equalDOF 21 881 1 2
equalDOF 31 890 1 2
equalDOF 41 899 1 2
equalDOF 51 908 1 2
equalDOF 61 917 1 2
equalDOF 71 926 1 2
equalDOF 81 935 1 2
equalDOF 91 944 1 2
equalDOF 101 953 1 2
equalDOF 111 962 1 2
equalDOF 121 971 1 2
equalDOF 131 980 1 2
equalDOF 141 990 1 2
equalDOF 151 1000 1 2
equalDOF 161 1010 1 2
equalDOF 171 1020 1 2
equalDOF 181 1030 1 2
equalDOF 191 1040 1 2
equalDOF 201 1050 1 2
equalDOF 211 1060 1 2
equalDOF 221 1070 1 2
equalDOF 231 1080 1 2
equalDOF 241 1090 1 2
equalDOF 251 1100 1 2
equalDOF 261 1110 1 2
equalDOF 271 1120 1 2
puts "Finished creating equal DOF (soil to soil)..."
#-----------------------------------------------------------------------------------------
# 6. creat zerolength elements nodes
#-----------------------------------------------------------------------------------------
node 1121 +9.900000E+000 +1.334000E+001
node 1122 +9.900000E+000 +1.374000E+001
node 1123 +9.900000E+000 +1.414000E+001
node 1124 +9.900000E+000 +1.454000E+001
node 1125 +9.900000E+000 +1.494000E+001
node 1126 +9.900000E+000 +1.534000E+001
node 1127 +9.900000E+000 +1.574000E+001
node 1128 +9.900000E+000 +1.614000E+001
node 1129 +9.900000E+000 +1.654000E+001
node 1130 +9.900000E+000 +1.694000E+001
node 1131 +9.900000E+000 +1.734000E+001
node 1132 +9.900000E+000 +1.774000E+001
node 1133 +9.900000E+000 +1.814000E+001
node 1134 +9.900000E+000 +1.854000E+001
node 1135 +9.900000E+000 +1.894000E+001
node 1136 +2.190000E+001 +1.334000E+001
node 1137 +2.190000E+001 +1.374000E+001
node 1138 +2.190000E+001 +1.414000E+001
node 1139 +2.190000E+001 +1.454000E+001
node 1140 +2.190000E+001 +1.494000E+001
node 1141 +2.190000E+001 +1.534000E+001
node 1142 +2.190000E+001 +1.574000E+001
node 1143 +2.190000E+001 +1.614000E+001
node 1144 +2.190000E+001 +1.654000E+001
node 1145 +2.190000E+001 +1.694000E+001
node 1146 +2.190000E+001 +1.734000E+001
node 1147 +2.190000E+001 +1.774000E+001
node 1148 +2.190000E+001 +1.814000E+001
node 1149 +2.190000E+001 +1.854000E+001
node 1150 +2.190000E+001 +1.894000E+001
node 1151 +3.940000E+001 +1.334000E+001
node 1152 +3.940000E+001 +1.374000E+001
node 1153 +3.940000E+001 +1.414000E+001
node 1154 +3.940000E+001 +1.454000E+001
node 1155 +3.940000E+001 +1.494000E+001
node 1156 +3.940000E+001 +1.534000E+001
node 1157 +3.940000E+001 +1.574000E+001
node 1158 +3.940000E+001 +1.614000E+001
node 1159 +3.940000E+001 +1.654000E+001
node 1160 +3.940000E+001 +1.694000E+001
node 1161 +3.940000E+001 +1.734000E+001
node 1162 +3.940000E+001 +1.774000E+001
node 1163 +3.940000E+001 +1.814000E+001
node 1164 +3.940000E+001 +1.854000E+001
node 1165 +3.940000E+001 +1.894000E+001
node 1166 +5.140000E+001 +1.334000E+001
node 1167 +5.140000E+001 +1.374000E+001
node 1168 +5.140000E+001 +1.414000E+001
node 1169 +5.140000E+001 +1.454000E+001
node 1170 +5.140000E+001 +1.494000E+001
node 1171 +5.140000E+001 +1.534000E+001
node 1172 +5.140000E+001 +1.574000E+001
node 1173 +5.140000E+001 +1.614000E+001
node 1174 +5.140000E+001 +1.654000E+001
node 1175 +5.140000E+001 +1.694000E+001
node 1176 +5.140000E+001 +1.734000E+001
node 1177 +5.140000E+001 +1.774000E+001
node 1178 +5.140000E+001 +1.814000E+001
node 1179 +5.140000E+001 +1.854000E+001
node 1180 +5.140000E+001 +1.894000E+001
node 1181 +1.190000E+001 +1.334000E+001
node 1182 +1.390000E+001 +1.334000E+001
node 1183 +1.590000E+001 +1.334000E+001
node 1184 +1.790000E+001 +1.334000E+001
node 1185 +1.990000E+001 +1.334000E+001
node 1186 +4.140000E+001 +1.334000E+001
node 1187 +4.340000E+001 +1.334000E+001
node 1188 +4.540000E+001 +1.334000E+001
node 1189 +4.740000E+001 +1.334000E+001
node 1190 +4.940000E+001 +1.334000E+001
puts "Finished creating zerolength elements nodes..."
#-----------------------------------------------------------------------------------------
# 7. Define zerolength elements materials
#-----------------------------------------------------------------------------------------
# Material: matTag E
uniaxialMaterial ENT 2 +4.000000E+004
# Material: matTag C alpha
uniaxialMaterial Elastic 3 10
# Material "interaction": matTag tags
uniaxialMaterial Parallel 4 2 3
puts "Finished creating all Zerolength elements materials..."
#-----------------------------------------------------------------------------------------
# 8. Define zerolength elements
#-----------------------------------------------------------------------------------------
element zeroLength 1021 140 1121 -mat 4 -dir 1
element zeroLength 1022 150 1122 -mat 4 -dir 1
element zeroLength 1023 160 1123 -mat 4 -dir 1
element zeroLength 1024 170 1124 -mat 4 -dir 1
element zeroLength 1025 180 1125 -mat 4 -dir 1
element zeroLength 1026 190 1126 -mat 4 -dir 1
element zeroLength 1027 200 1127 -mat 4 -dir 1
element zeroLength 1028 210 1128 -mat 4 -dir 1
element zeroLength 1029 220 1129 -mat 4 -dir 1
element zeroLength 1030 230 1130 -mat 4 -dir 1
element zeroLength 1031 240 1131 -mat 4 -dir 1
element zeroLength 1032 250 1132 -mat 4 -dir 1
element zeroLength 1033 260 1133 -mat 4 -dir 1
element zeroLength 1034 270 1134 -mat 4 -dir 1
element zeroLength 1035 280 1135 -mat 4 -dir 1
element zeroLength 1036 359 1181 -mat 4 -dir 2
element zeroLength 1037 360 1182 -mat 4 -dir 2
element zeroLength 1038 361 1183 -mat 4 -dir 2
element zeroLength 1039 362 1184 -mat 4 -dir 2
element zeroLength 1040 363 1185 -mat 4 -dir 2
element zeroLength 1041 756 1150 -mat 4 -dir 1
element zeroLength 1042 741 1149 -mat 4 -dir 1
element zeroLength 1043 726 1148 -mat 4 -dir 1
element zeroLength 1044 711 1147 -mat 4 -dir 1
element zeroLength 1045 696 1146 -mat 4 -dir 1
element zeroLength 1046 681 1145 -mat 4 -dir 1
element zeroLength 1047 666 1144 -mat 4 -dir 1
element zeroLength 1048 651 1143 -mat 4 -dir 1
element zeroLength 1049 636 1142 -mat 4 -dir 1
element zeroLength 1050 621 1141 -mat 4 -dir 1
element zeroLength 1051 606 1140 -mat 4 -dir 1
element zeroLength 1052 591 1139 -mat 4 -dir 1
element zeroLength 1053 576 1138 -mat 4 -dir 1
element zeroLength 1054 561 1137 -mat 4 -dir 1
element zeroLength 1055 364 1136 -mat 4 -dir 1
element zeroLength 1056 560 1151 -mat 4 -dir 1
element zeroLength 1057 575 1152 -mat 4 -dir 1
element zeroLength 1058 590 1153 -mat 4 -dir 1
element zeroLength 1059 605 1154 -mat 4 -dir 1
element zeroLength 1060 620 1155 -mat 4 -dir 1
element zeroLength 1061 635 1156 -mat 4 -dir 1
element zeroLength 1062 650 1157 -mat 4 -dir 1
element zeroLength 1063 665 1158 -mat 4 -dir 1
element zeroLength 1064 680 1159 -mat 4 -dir 1
element zeroLength 1065 695 1160 -mat 4 -dir 1
element zeroLength 1066 710 1161 -mat 4 -dir 1
element zeroLength 1067 725 1162 -mat 4 -dir 1
element zeroLength 1068 740 1163 -mat 4 -dir 1
element zeroLength 1069 755 1164 -mat 4 -dir 1
element zeroLength 1070 770 1165 -mat 4 -dir 1
element zeroLength 1071 849 1186 -mat 4 -dir 2
element zeroLength 1072 850 1187 -mat 4 -dir 2
element zeroLength 1073 851 1188 -mat 4 -dir 2
element zeroLength 1074 852 1189 -mat 4 -dir 2
element zeroLength 1075 853 1190 -mat 4 -dir 2
element zeroLength 1076 854 1166 -mat 4 -dir 1
element zeroLength 1077 981 1167 -mat 4 -dir 1
element zeroLength 1078 991 1168 -mat 4 -dir 1
element zeroLength 1079 1001 1169 -mat 4 -dir 1
element zeroLength 1080 1011 1170 -mat 4 -dir 1
element zeroLength 1081 1021 1171 -mat 4 -dir 1
element zeroLength 1082 1031 1172 -mat 4 -dir 1
element zeroLength 1083 1041 1173 -mat 4 -dir 1
element zeroLength 1084 1051 1174 -mat 4 -dir 1
element zeroLength 1085 1061 1175 -mat 4 -dir 1
element zeroLength 1086 1071 1176 -mat 4 -dir 1
element zeroLength 1087 1081 1177 -mat 4 -dir 1
element zeroLength 1088 1091 1178 -mat 4 -dir 1
element zeroLength 1089 1101 1179 -mat 4 -dir 1
element zeroLength 1090 1111 1180 -mat 4 -dir 1
puts "Finished creating zerolength elements ..."
#-----------------------------------------------------------------------------------------
# 9. create the 1st analysis method & analyze
#-----------------------------------------------------------------------------------------
set gamma 1.6
# create the SOE, ConstraintHandler, Integrator, Algorithm and Numberer
integrator Newmark $gamma [expr pow($gamma+0.5, 2)/4]
test EnergyIncr 1.0e-8 400 1;
constraints Transformation
algorithm Newton
numberer RCM
system ProfileSPD
analysis Transient
analyze 2000 5.0e-5
puts "finished 1st analysis (static)"
#-----------------------------------------------------------------------------------------
# 10. update Material Stage
#-----------------------------------------------------------------------------------------
updateMaterialStage -material 1 -stage 1
#-----------------------------------------------------------------------------------------
# 11. Creat Wall and base nodes
#-----------------------------------------------------------------------------------------
model BasicBuilder -ndm 2 -ndf 3
node 1191 9.90E+00 1.33E+01
node 1192 9.90E+00 1.37E+01
node 1193 9.90E+00 1.41E+01
node 1194 9.90E+00 1.45E+01
node 1195 9.90E+00 1.49E+01
node 1196 9.90E+00 1.53E+01
node 1197 9.90E+00 1.57E+01
node 1198 9.90E+00 1.61E+01
node 1199 9.90E+00 1.65E+01
node 1200 9.90E+00 1.69E+01
node 1201 9.90E+00 1.73E+01
node 1202 9.90E+00 1.77E+01
node 1203 9.90E+00 1.81E+01
node 1204 9.90E+00 1.85E+01
node 1205 9.90E+00 1.89E+01
node 1206 2.19E+01 1.33E+01
node 1207 2.19E+01 1.37E+01
node 1208 2.19E+01 1.41E+01
node 1209 2.19E+01 1.45E+01
node 1210 2.19E+01 1.49E+01
node 1211 2.19E+01 1.53E+01
node 1212 2.19E+01 1.57E+01
node 1213 2.19E+01 1.61E+01
node 1214 2.19E+01 1.65E+01
node 1215 2.19E+01 1.69E+01
node 1216 2.19E+01 1.73E+01
node 1217 2.19E+01 1.77E+01
node 1218 2.19E+01 1.81E+01
node 1219 2.19E+01 1.85E+01
node 1220 2.19E+01 1.89E+01
node 1221 3.94E+01 1.33E+01
node 1222 3.94E+01 1.37E+01
node 1223 3.94E+01 1.41E+01
node 1224 3.94E+01 1.45E+01
node 1225 3.94E+01 1.49E+01
node 1226 3.94E+01 1.53E+01
node 1227 3.94E+01 1.57E+01
node 1228 3.94E+01 1.61E+01
node 1229 3.94E+01 1.65E+01
node 1230 3.94E+01 1.69E+01
node 1231 3.94E+01 1.73E+01
node 1232 3.94E+01 1.77E+01
node 1233 3.94E+01 1.81E+01
node 1234 3.94E+01 1.85E+01
node 1235 3.94E+01 1.89E+01
node 1236 5.14E+01 1.33E+01
node 1237 5.14E+01 1.37E+01
node 1238 5.14E+01 1.41E+01
node 1239 5.14E+01 1.45E+01
node 1240 5.14E+01 1.49E+01
node 1241 5.14E+01 1.53E+01
node 1242 5.14E+01 1.57E+01
node 1243 5.14E+01 1.61E+01
node 1244 5.14E+01 1.65E+01
node 1245 5.14E+01 1.69E+01
node 1246 5.14E+01 1.73E+01
node 1247 5.14E+01 1.77E+01
node 1248 5.14E+01 1.81E+01
node 1249 5.14E+01 1.85E+01
node 1250 5.14E+01 1.89E+01
node 1251 1.19E+01 1.33E+01
node 1252 1.39E+01 1.33E+01
node 1253 1.59E+01 1.33E+01
node 1254 1.79E+01 1.33E+01
node 1255 1.99E+01 1.33E+01
node 1256 4.14E+01 1.33E+01
node 1257 4.34E+01 1.33E+01
node 1258 4.54E+01 1.33E+01
node 1259 4.74E+01 1.33E+01
node 1260 4.94E+01 1.33E+01
puts "Finished creating wall and base nodes ..."
#-----------------------------------------------------------------------------------------
# 12. define equal DOF (soil to wall)
#-----------------------------------------------------------------------------------------
equalDOF 1191 1121 1 2
equalDOF 1192 1122 1 2
equalDOF 1193 1123 1 2
equalDOF 1194 1124 1 2
equalDOF 1195 1125 1 2
equalDOF 1196 1126 1 2
equalDOF 1197 1127 1 2
equalDOF 1198 1128 1 2
equalDOF 1199 1129 1 2
equalDOF 1200 1130 1 2
equalDOF 1201 1131 1 2
equalDOF 1202 1132 1 2
equalDOF 1203 1133 1 2
equalDOF 1204 1134 1 2
equalDOF 1205 1135 1 2
equalDOF 1206 1181 1 2
equalDOF 1207 1182 1 2
equalDOF 1208 1183 1 2
equalDOF 1209 1184 1 2
equalDOF 1210 1185 1 2
equalDOF 1211 1150 1 2
equalDOF 1212 1149 1 2
equalDOF 1213 1148 1 2
equalDOF 1214 1147 1 2
equalDOF 1215 1146 1 2
equalDOF 1216 1145 1 2
equalDOF 1217 1144 1 2
equalDOF 1218 1143 1 2
equalDOF 1219 1142 1 2
equalDOF 1220 1141 1 2
equalDOF 1221 1140 1 2
equalDOF 1222 1139 1 2
equalDOF 1223 1138 1 2
equalDOF 1224 1137 1 2
equalDOF 1225 1136 1 2
equalDOF 1226 1151 1 2
equalDOF 1227 1152 1 2
equalDOF 1228 1153 1 2
equalDOF 1229 1154 1 2
equalDOF 1230 1155 1 2
equalDOF 1231 1156 1 2
equalDOF 1232 1157 1 2
equalDOF 1233 1158 1 2
equalDOF 1234 1159 1 2
equalDOF 1235 1160 1 2
equalDOF 1236 1161 1 2
equalDOF 1237 1162 1 2
equalDOF 1238 1163 1 2
equalDOF 1239 1164 1 2
equalDOF 1240 1165 1 2
equalDOF 1241 1186 1 2
equalDOF 1242 1187 1 2
equalDOF 1243 1188 1 2
equalDOF 1244 1189 1 2
equalDOF 1245 1190 1 2
equalDOF 1246 1166 1 2
equalDOF 1247 1167 1 2
equalDOF 1248 1168 1 2
equalDOF 1249 1169 1 2
equalDOF 1250 1170 1 2
equalDOF 1251 1171 1 2
equalDOF 1252 1172 1 2
equalDOF 1253 1173 1 2
equalDOF 1254 1174 1 2
equalDOF 1255 1175 1 2
equalDOF 1256 1176 1 2
equalDOF 1257 1177 1 2
equalDOF 1258 1178 1 2
equalDOF 1259 1179 1 2
equalDOF 1260 1180 1 2
puts "Finished defining equal DOF (soil to wall) ..."
#-----------------------------------------------------------------------------------------
# 13. define wall and base elements
#-----------------------------------------------------------------------------------------
geomTransf Linear 1
# stiff wall elements
element elasticBeamColumn 1091 1191 1192 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1092 1192 1193 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1093 1193 1194 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1094 1194 1195 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1095 1195 1196 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1096 1196 1197 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1097 1197 1198 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1098 1198 1199 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1099 1199 1200 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1100 1200 1201 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1101 1201 1202 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1102 1202 1203 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1103 1203 1204 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1104 1204 1205 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1105 1205 1206 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1106 1206 1207 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1107 1207 1208 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1108 1208 1209 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1109 1209 1210 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1110 1210 1211 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1111 1211 1212 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1112 1212 1213 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1113 1213 1214 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1114 1214 1215 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1115 1215 1216 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1116 1216 1217 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1117 1217 1218 1.40E-01 70000000 2.43E-03 1
element elasticBeamColumn 1118 1218 1219 1.40E-01 70000000 2.43E-03 1
# base elements for stiff wall
element elasticBeamColumn 1119 1191 1251 2.50E-01 70000000 1.42E-02 1
element elasticBeamColumn 1120 1251 1252 2.50E-01 70000000 1.42E-02 1
element elasticBeamColumn 1121 1252 1253 2.50E-01 70000000 1.42E-02 1
element elasticBeamColumn 1122 1253 1254 2.50E-01 70000000 1.42E-02 1
element elasticBeamColumn 1123 1254 1255 2.50E-01 70000000 1.42E-02 1
element elasticBeamColumn 1124 1255 1205 2.50E-01 70000000 1.42E-02 1
# flexible wall elements
element elasticBeamColumn 1125 1220 1221 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1126 1221 1222 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1127 1222 1223 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1128 1223 1224 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1129 1224 1225 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1130 1225 1226 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1131 1226 1227 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1132 1227 1228 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1133 1228 1229 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1134 1229 1230 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1135 1230 1231 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1136 1231 1232 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1137 1232 1233 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1138 1233 1234 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1139 1234 1235 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1140 1235 1236 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1141 1236 1237 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1142 1237 1238 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1143 1238 1239 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1144 1239 1240 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1145 1240 1241 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1146 1241 1242 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1147 1242 1243 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1148 1243 1244 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1149 1244 1245 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1150 1245 1246 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1151 1246 1247 8.00E-02 70000000 4.26E-04 1
element elasticBeamColumn 1152 1247 1248 8.00E-02 70000000 4.26E-04 1
# base elements for flexible wall
element elasticBeamColumn 1153 1220 1255 2.50E-01 70000000 1.42E-02 1
element elasticBeamColumn 1154 1255 1256 2.50E-01 70000000 1.42E-02 1
element elasticBeamColumn 1155 1256 1257 2.50E-01 70000000 1.42E-02 1
element elasticBeamColumn 1156 1257 1258 2.50E-01 70000000 1.42E-02 1
element elasticBeamColumn 1157 1258 1259 2.50E-01 70000000 1.42E-02 1
element elasticBeamColumn 1158 1259 1235 2.50E-01 70000000 1.42E-02 1
puts "Finished creating wall and base elements ..."
#-----------------------------------------------------------------------------------------
# 14. settime and wipeAnalysis
#-----------------------------------------------------------------------------------------
wipeAnalysis
loadConst -time 0.0
#-----------------------------------------------------------------------------------------
# 15. apply the weight of structure
#-----------------------------------------------------------------------------------------
set wstiffwall -2.38
set wflexiblewall -2.03
set wbase -19.62
set wbs -22
set wbf -21.65
set wstiffwallh -1.19
set wflexiblewallh -1.015
pattern Plain 5 Linear {
load 1191 0.00E+00 $wbs 0.00E+00
load 1192 0.00E+00 $wstiffwall 0.00E+00
load 1193 0.00E+00 $wstiffwall 0.00E+00
load 1194 0.00E+00 $wstiffwall 0.00E+00
load 1195 0.00E+00 $wstiffwall 0.00E+00
load 1196 0.00E+00 $wstiffwall 0.00E+00
load 1197 0.00E+00 $wstiffwall 0.00E+00
load 1198 0.00E+00 $wstiffwall 0.00E+00
load 1199 0.00E+00 $wstiffwall 0.00E+00
load 1200 0.00E+00 $wstiffwall 0.00E+00
load 1201 0.00E+00 $wstiffwall 0.00E+00
load 1202 0.00E+00 $wstiffwall 0.00E+00
load 1203 0.00E+00 $wstiffwall 0.00E+00
load 1204 0.00E+00 $wstiffwall 0.00E+00
load 1205 0.00E+00 $wstiffwallh 0.00E+00
load 1206 0.00E+00 $wbase 0.00E+00
load 1207 0.00E+00 $wbase 0.00E+00
load 1208 0.00E+00 $wbase 0.00E+00
load 1209 0.00E+00 $wbase 0.00E+00
load 1210 0.00E+00 $wbase 0.00E+00
load 1211 0.00E+00 $wstiffwallh 0.00E+00
load 1212 0.00E+00 $wstiffwall 0.00E+00
load 1213 0.00E+00 $wstiffwall 0.00E+00
load 1214 0.00E+00 $wstiffwall 0.00E+00
load 1215 0.00E+00 $wstiffwall 0.00E+00
load 1216 0.00E+00 $wstiffwall 0.00E+00
load 1217 0.00E+00 $wstiffwall 0.00E+00
load 1218 0.00E+00 $wstiffwall 0.00E+00
load 1219 0.00E+00 $wstiffwall 0.00E+00
load 1220 0.00E+00 $wstiffwall 0.00E+00
load 1221 0.00E+00 $wstiffwall 0.00E+00
load 1222 0.00E+00 $wstiffwall 0.00E+00
load 1223 0.00E+00 $wstiffwall 0.00E+00
load 1224 0.00E+00 $wstiffwall 0.00E+00
load 1225 0.00E+00 $wbs 0.00E+00
load 1226 0.00E+00 $wbf 0.00E+00
load 1227 0.00E+00 $wflexiblewall 0.00E+00
load 1228 0.00E+00 $wflexiblewall 0.00E+00
load 1229 0.00E+00 $wflexiblewall 0.00E+00
load 1230 0.00E+00 $wflexiblewall 0.00E+00
load 1231 0.00E+00 $wflexiblewall 0.00E+00
load 1232 0.00E+00 $wflexiblewall 0.00E+00
load 1233 0.00E+00 $wflexiblewall 0.00E+00
load 1234 0.00E+00 $wflexiblewall 0.00E+00
load 1235 0.00E+00 $wflexiblewall 0.00E+00
load 1236 0.00E+00 $wflexiblewall 0.00E+00
load 1237 0.00E+00 $wflexiblewall 0.00E+00
load 1238 0.00E+00 $wflexiblewall 0.00E+00
load 1239 0.00E+00 $wflexiblewall 0.00E+00
load 1240 0.00E+00 $wflexiblewallh 0.00E+00
load 1241 0.00E+00 $wbase 0.00E+00
load 1242 0.00E+00 $wbase 0.00E+00
load 1243 0.00E+00 $wbase 0.00E+00
load 1244 0.00E+00 $wbase 0.00E+00
load 1245 0.00E+00 $wbase 0.00E+00
load 1246 0.00E+00 $wbf 0.00E+00
load 1247 0.00E+00 $wflexiblewall 0.00E+00
load 1248 0.00E+00 $wflexiblewall 0.00E+00
load 1249 0.00E+00 $wflexiblewall 0.00E+00
load 1250 0.00E+00 $wflexiblewall 0.00E+00
load 1251 0.00E+00 $wflexiblewall 0.00E+00
load 1252 0.00E+00 $wflexiblewall 0.00E+00
load 1253 0.00E+00 $wflexiblewall 0.00E+00
load 1254 0.00E+00 $wflexiblewall 0.00E+00
load 1255 0.00E+00 $wflexiblewall 0.00E+00
load 1256 0.00E+00 $wflexiblewall 0.00E+00
load 1257 0.00E+00 $wflexiblewall 0.00E+00
load 1258 0.00E+00 $wflexiblewall 0.00E+00
load 1259 0.00E+00 $wflexiblewall 0.00E+00
load 1260 0.00E+00 $wflexiblewallh 0.00E+00
}
puts "Finished applying the structures weight ..."
#-----------------------------------------------------------------------------------------
# 16. Recorder
#-----------------------------------------------------------------------------------------
source recorders1.txt
#-----------------------------------------------------------------------------------------
# 17. create the 2nd analysis method & analyze
#-----------------------------------------------------------------------------------------
constraints Penalty +1.000000E+012 +1.000000E+012
test NormDispIncr +1.000000E-006 10 0
integrator LoadControl +1.000000E+000 1 +1.000000E+000 +1.000000E+000
algorithm Newton
numberer Plain
system ProfileSPD
analysis Static
analyze 100
#-----------------------------------------------------------------------------------------
# 18. earthquake excitation
#-----------------------------------------------------------------------------------------
wipeAnalysis
loadConst -time 0.0
timeSeries Path 3 -dt +5.000000E-2 -filePath tabasFP.txt -factor +9.810000E+000
# LoadPattern "tabas": patternTag dir tsTag
pattern UniformExcitation 3 1 -accel 3
#-----------------------------------------------------------------------------------------
# 19. Define recorder
#-----------------------------------------------------------------------------------------
source recorders2.txt
#-----------------------------------------------------------------------------------------
# 20. create the 3rd analysis method & analyze
#-----------------------------------------------------------------------------------------
constraints Penalty +1.000000E+012 +1.000000E+012
test NormDispIncr +1.000000E-004 10 0
integrator Newmark 0.6 0.3025 0 0 0.001 0
algorithm Newton
numberer Plain
system ProfileSPD
analysis Transient
# Define damping parameters
# -------------------------
# parameter set "DampingParam01": alphaM betaK betaKinit betaKcomm
rayleigh +0.000000E+000 +4.400000E-003 +0.000000E+000 +0.000000E+000
analyze 2500 +6.000000E-003