joining 2dof with 3dof elements

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

Moderators: silvia, selimgunay, Moderators

Post Reply
partla
Posts: 35
Joined: Wed Feb 01, 2012 3:04 am
Location: Urmia university

joining 2dof with 3dof elements

Post by partla »

Hi,
I’m trying to model a retaining wall and I want to use zero length elements with a parallel material (an elastic no tension material parallel with a viscous material) as interface elements between wall and soil elements. I use quadrilateral elements for soil and elastic beam column elements for wall. I know that in this model there are two types of domains, the two-DOF domain where the quad elements and their respective nodes have been created, and the three-DOF domain where the wall elements and their respective nodes have been created.
In order for the two domains to be connected, I use equalDOF command to connect nodes on the common boundaries of the two domains. But I got a warning. The warning is that “ nodes … and … have differing dof at ends for Zerolength … “.
In your idea what’s the problem? How can I fix it?

Thanks in advance
Top
vesna
Posts: 3033
Joined: Tue May 23, 2006 11:23 am
Location: UC Berkeley

Re: joining 2dof with 3dof elements

Post by vesna »

Can you copy the full warning message?
Top
partla
Posts: 35
Joined: Wed Feb 01, 2012 3:04 am
Location: Urmia university

Re: joining 2dof with 3dof elements

Post by partla »

Hi Vesna,

I think I solved this problem. I defined 3 nodes with same coordinate and joined 2 of nodes with a zero length element and joined one side of zerolength element with the third node. After this changes it hasn't that warning but I've got other errors with analyze commands. Is defining 3 nodes right about that problem?
Top
vesna
Posts: 3033
Joined: Tue May 23, 2006 11:23 am
Location: UC Berkeley

Re: joining 2dof with 3dof elements

Post by vesna »

You can have three nodes with the same coordinates but you need to know how to set up the model and analysis properly. I think it is easier to have only two nodes but to use Series material (http://opensees.berkeley.edu/wiki/index ... s_Material) to model the same effect.
Top
partla
Posts: 35
Joined: Wed Feb 01, 2012 3:04 am
Location: Urmia university

Re: joining 2dof with 3dof elements

Post by partla »

I'm trying to model the project of "PEER 2008/104" for verifying my own project. I wrote the script of that but I couldn't setup the analysis format. Here is the script. 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
Top
partla
Posts: 35
Joined: Wed Feb 01, 2012 3:04 am
Location: Urmia university

Re: joining 2dof with 3dof elements

Post by partla »

No body to help me? How should I set up the analysis commands in the upper script?
Top
partla
Posts: 35
Joined: Wed Feb 01, 2012 3:04 am
Location: Urmia university

Re: joining 2dof with 3dof elements

Post by partla »

any idea about this problem?
Top
vesna
Posts: 3033
Joined: Tue May 23, 2006 11:23 am
Location: UC Berkeley

Re: joining 2dof with 3dof elements

Post by vesna »

You will hardly find someone to look at as big script as yours.
Top
partla
Posts: 35
Joined: Wed Feb 01, 2012 3:04 am
Location: Urmia university

Re: joining 2dof with 3dof elements

Post by partla »

Thank you vesna.
OK you are right. let me explain my problem and I hope you will reply. I really need it.
As I explained in my 1st post I used 3 nodes with same coordination for defining a wall (element elasticBeamColumn) in a soil mesh. and I used a zero length element with an spring and dashpot in parallel as interface elements (I'm trying to model the project of "PEER 2008/104" for verifying my own project). But I couldn't setup the analyses commands (as you said before) for static and dynamic analyses.
I will be happy hearing from you
Top
partla
Posts: 35
Joined: Wed Feb 01, 2012 3:04 am
Location: Urmia university

Re: joining 2dof with 3dof elements

Post by partla »

and here is the warnings that I got from analyses of above code.

ProfileSPDLinDirectSolver::solve() - aii < 0 (i, aii): (1, 0)
WARNING Linear::solveCurrentStep() -the LinearSOE failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 0.1
OpenSees > analyze failed, returned: -3 error flag
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 17
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 1
OpenSees > analyze failed, returned: -3 error flag
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 17
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 0.006
OpenSees > analyze failed, returned: -3 error flag
-3
Top
vesna
Posts: 3033
Joined: Tue May 23, 2006 11:23 am
Location: UC Berkeley

Re: joining 2dof with 3dof elements

Post by vesna »

Why don't you try using two nodes with the same coordinates and define your spring and a dashpot using "Parallel material".

If you have problem with your analysis, post the warning messages that you get from OpenSees.
Top
partla
Posts: 35
Joined: Wed Feb 01, 2012 3:04 am
Location: Urmia university

Re: joining 2dof with 3dof elements

Post by partla »

thank you vesna

I've tried 2 nodes with same material but I got error with Dof. I'll try it again and post it's problems.

Thank you again
Top
vesna
Posts: 3033
Joined: Tue May 23, 2006 11:23 am
Location: UC Berkeley

Re: joining 2dof with 3dof elements

Post by vesna »

the error message you posted indicates instability of the model. make sure to tie nodes with the same coordinates using equalDOF command.
Top
partla
Posts: 35
Joined: Wed Feb 01, 2012 3:04 am
Location: Urmia university

Re: joining 2dof with 3dof elements

Post by partla »

I tie the left side of soil mesh to right side of it with equalDOF command. Is it right?
Top
vesna
Posts: 3033
Joined: Tue May 23, 2006 11:23 am
Location: UC Berkeley

Re: joining 2dof with 3dof elements

Post by vesna »

yes, that is ok. but you also need to tie two nodes with the same coordinates in all directions in which you did not assign a material.

i skimmed through your model and noticed that you are fixing nodes of 2dof domain in 3 dofs. Check on this!
Top
Post Reply

Return to “OpenSees.exe Users”