Code: Select all
wipe
if {[file exists output2]==0} {
file mkdir output2;
}
set g 9.81
set gravX 0.0
set gravY 0.0
set gravZ [expr -1*$g]
# --------------------------------------------------------------------------------------------------------------
# unit
# --------------------------------------------------------------------------------------------------------------
# Length : m
# Force : kN
# Moment : kNm
# Stress : kPa
# Mass : ton
model BasicBuilder -ndm 3 -ndf 4
node 1 0 1 1.5
node 2 0.333333 1 1.5
node 3 0 0.666667 1.5
node 4 0.25 0.75 1.5
node 5 0 1 1
node 6 0.25 0.482947 1.5
node 7 0.517053 0.75 1.5
node 8 0.333333 1 1
node 9 0 0.666667 1
node 10 0.25 0.75 1
node 11 0.666667 1 1.5
node 12 0 0.333333 1.5
node 13 0.25 0.482947 1
node 14 0.517053 0.75 1
node 15 0.765269 0.75 1.5
node 16 0.25 0.234731 1.5
node 17 0.666667 1 1
node 18 0 0.333333 1
node 19 0.646447 0.353553 1.5
node 20 0.566987 0.25 1.5
node 21 0.75 0.433013 1.5
node 22 0.765269 0.75 1
node 23 0.25 0.234731 1
node 24 1 1 1.5
node 25 0 0 1.5
node 26 0 1 0.5
node 27 0.87059 0.482963 1.5
node 28 0.517037 0.12941 1.5
node 29 0.717157 0.282843 1.5
node 30 0.25 0 1.5
node 31 1 0.75 1.5
node 32 0.65359 0.2 1.5
node 33 0.8 0.34641 1.5
node 34 0.646447 0.353553 1
node 35 0 0.666667 0.5
node 36 0.333333 1 0.5
node 37 0.25 0.75 0.5
node 38 0.566987 0.25 1
node 39 0.75 0.433013 1
node 40 0.896472 0.38637 1.5
node 41 0.61363 0.103528 1.5
node 42 0 0 1
node 43 1 1 1
node 44 1 0.5 1.5
node 45 0.5 0 1.5
node 46 0.87059 0.482963 1
node 47 0.517037 0.12941 1
node 48 0.717157 0.282843 1
node 49 0.25 0 1
node 50 1 0.75 1
node 51 0.8 0.34641 1
node 52 0.65359 0.2 1
node 53 0.25 0.482947 0.5
node 54 0.517053 0.75 0.5
node 55 0.6 -1.9916e-17 1.5
node 56 1 0.4 1.5
node 57 0.61363 0.103528 1
node 58 0.896472 0.38637 1
node 59 0.666667 1 0.5
node 60 0 0.333333 0.5
node 61 1 0.5 1
node 62 0.5 0 1
node 63 1.12941 0.482963 1.5
node 64 1.23473 0.75 1.5
node 65 1.10353 0.38637 1.5
node 66 0.6 -1.9916e-17 1
node 67 1 0.4 1
node 68 0.765269 0.75 0.5
node 69 0.25 0.234731 0.5
node 70 1.33333 1 1.5
node 71 1.12941 0.482963 1
node 72 0.646447 0.353553 0.5
node 73 1.23473 0.75 1
node 74 1.10353 0.38637 1
node 75 1.2 0.34641 1.5
node 76 1.25 0.433013 1.5
node 77 0.75 0.433013 0.5
node 78 0.566987 0.25 0.5
node 79 0 0 0.5
node 80 1 1 0.5
node 81 1 0 1.5
node 82 0.87059 0.482963 0.5
node 83 0.517037 0.12941 0.5
node 84 1.33333 1 1
node 85 0.717157 0.282843 0.5
node 86 1 0.75 0.5
node 87 0.25 0 0.5
node 88 0.65359 0.2 0.5
node 89 0.8 0.34641 0.5
node 90 1.2 0.34641 1
node 91 1.25 0.433013 1
node 92 1.28284 0.282843 1.5
node 93 0.61363 0.103528 0.5
node 94 0.896472 0.38637 0.5
node 95 1.35355 0.353553 1.5
node 96 1 0 1
node 97 1 0.5 0.5
node 98 0 1 0
node 99 0.5 0 0.5
node 100 1.48295 0.75 1.5
node 101 1 0.4 0.5
node 102 0.6 -1.9916e-17 0.5
node 103 0.333333 1 0
node 104 0 0.666667 0
node 105 0.25 0.75 0
node 106 1.28284 0.282843 1
node 107 1.34641 0.2 1.5
node 108 1.35355 0.353553 1
node 109 1.48295 0.75 1
node 110 1.12941 0.482963 0.5
node 111 0.25 0.482947 0
node 112 0.517053 0.75 0
node 113 1.23473 0.75 0.5
node 114 1.10353 0.38637 0.5
node 115 1.43301 0.25 1.5
node 116 0.666667 1 0
node 117 0 0.333333 0
node 118 1.34641 0.2 1
node 119 1.38637 0.103528 1.5
node 120 1.66667 1 1.5
node 121 1.33333 1 0.5
node 122 1.43301 0.25 1
node 123 1.2 0.34641 0.5
node 124 1.25 0.433013 0.5
node 125 0.25 0.234731 0
node 126 0.765269 0.75 0
node 127 1.48296 0.12941 1.5
node 128 1.4 7.66952e-18 1.5
node 129 1.38637 0.103528 1
node 130 1 0 0.5
node 131 1.66667 1 1
node 132 0.646447 0.353553 0
node 133 1.75 0.75 1.5
node 134 0.75 0.433013 0
node 135 0.566987 0.25 0
node 136 1.28284 0.282843 0.5
node 137 1.48296 0.12941 1
node 138 1.4 7.66952e-18 1
node 139 1.35355 0.353553 0.5
node 140 1 1 0
node 141 1.5 0 1.5
node 142 0 0 0
node 143 1.48295 0.75 0.5
node 144 0.517037 0.12941 0
node 145 0.87059 0.482963 0
node 146 0.717157 0.282843 0
node 147 0.25 0 0
node 148 1 0.75 0
node 149 0.65359 0.2 0
node 150 0.8 0.34641 0
node 151 1.75 0.482947 1.5
node 152 1.75 0.75 1
node 153 0.61363 0.103528 0
node 154 0.896472 0.38637 0
node 155 1.34641 0.2 0.5
node 156 0.5 0 0
node 157 1 0.5 0
node 158 1.5 0 1
node 159 1.75 0.482947 1
node 160 1 0.4 0
node 161 0.6 -1.9916e-17 0
node 162 1.43301 0.25 0.5
node 163 1.75 0.234731 1.5
node 164 1.38637 0.103528 0.5
node 165 1.66667 1 0.5
node 166 1.12941 0.482963 0
node 167 1.23473 0.75 0
node 168 1.10353 0.38637 0
node 169 1.75 0.234731 1
node 170 1.48296 0.12941 0.5
node 171 1.4 7.66952e-18 0.5
node 172 2 1 1.5
node 173 1.33333 1 0
node 174 1.75 0 1.5
node 175 2 0.666667 1.5
node 176 1.2 0.34641 0
node 177 1.75 0.75 0.5
node 178 1.25 0.433013 0
node 179 1 0 0
node 180 2 1 1
node 181 1.5 0 0.5
node 182 1.75 0 1
node 183 1.75 0.482947 0.5
node 184 2 0.666667 1
node 185 1.28284 0.282843 0
node 186 2 0.333333 1.5
node 187 1.35355 0.353553 0
node 188 1.48295 0.75 0
node 189 1.75 0.234731 0.5
node 190 2 0.333333 1
node 191 1.34641 0.2 0
node 192 1.43301 0.25 0
node 193 1.38637 0.103528 0
node 194 2 0 1.5
node 195 2 1 0.5
node 196 1.66667 1 0
node 197 1.75 0 0.5
node 198 2 0.666667 0.5
node 199 1.48296 0.12941 0
node 200 1.4 7.66952e-18 0
node 201 2 0 1
node 202 1.75 0.75 0
node 203 2 0.333333 0.5
node 204 1.5 0 0
node 205 1.75 0.482947 0
node 206 1.75 0.234731 0
node 207 2 0 0.5
node 208 2 1 0
node 209 1.75 0 0
node 210 2 0.666667 0
node 211 2 0.333333 0
node 212 2 0 0
set numnp 212
set SnodeList {}
for {set i 1} {$i <= $numnp} {incr i 1} {
lappend SnodeList [expr $i]
}
#-----------------------------------------------------------
#Boundary conditions
#-----------------------------------------------------------
fixX [expr 0.0] 1 0 0 0 -tol 1e-10
fixX [expr 2.0] 1 0 0 0 -tol 1e-10
fixY [expr 0.0] 0 1 0 0 -tol 1e-10
fixY [expr 1.0] 0 1 0 0 -tol 1e-10
fixZ [expr 1.5] 0 0 1 0 -tol 1e-10
#------------------------------------------------------------
#Hydraulic boundary conditions -no excess pore pressure at mudline
#-------------------------------------------------------------
fix 1 0 0 0 1
fix 2 0 0 0 1
fix 3 0 0 0 1
fix 4 0 0 0 1
fix 6 0 0 0 1
fix 7 0 0 0 1
fix 11 0 0 0 1
fix 12 0 0 0 1
fix 15 0 0 0 1
fix 16 0 0 0 1
fix 19 0 0 0 1
fix 20 0 0 0 1
fix 21 0 0 0 1
fix 24 0 0 0 1
fix 25 0 0 0 1
fix 27 0 0 0 1
fix 28 0 0 0 1
fix 29 0 0 0 1
fix 30 0 0 0 1
fix 31 0 0 0 1
fix 32 0 0 0 1
fix 33 0 0 0 1
fix 40 0 0 0 1
fix 41 0 0 0 1
fix 44 0 0 0 1
fix 45 0 0 0 1
fix 55 0 0 0 1
fix 56 0 0 0 1
fix 63 0 0 0 1
fix 64 0 0 0 1
fix 65 0 0 0 1
fix 70 0 0 0 1
fix 75 0 0 0 1
fix 76 0 0 0 1
fix 81 0 0 0 1
fix 92 0 0 0 1
fix 95 0 0 0 1
fix 100 0 0 0 1
fix 107 0 0 0 1
fix 115 0 0 0 1
fix 119 0 0 0 1
fix 120 0 0 0 1
fix 127 0 0 0 1
fix 128 0 0 0 1
fix 133 0 0 0 1
fix 141 0 0 0 1
fix 151 0 0 0 1
fix 163 0 0 0 1
fix 172 0 0 0 1
fix 174 0 0 0 1
fix 175 0 0 0 1
fix 186 0 0 0 1
fix 194 0 0 0 1
#Material Parameters
set G0 95
set nu 0.05
set e_init 0.740
set Mc 1.35
set c 0.81
set lambda_c 0.055
set e0 1.035
set ksi 0.36
set P_atm 101.3
set m 0.01
set h0 7.60
set ch 0.97
set nb 1.2
set A0 0.74
set nd 1.79
set z_max 5
set cz 800
set Den 1.71
set fBulk 2.2e6; #kN/m3
set n [expr $e_init/(1+$e_init)]
set Bc [expr $fBulk/$n]
set fmass 1
set perm [expr 1.0e-6/$g/$fmass]
set alpha 1.0e-5
#Material Definition
nDMaterial ManzariDafalias 1 $G0 $nu $e_init $Mc $c $lambda_c $e0 $ksi $P_atm $m $h0 $ch $nb $A0 $nd $z_max $cz $Den
#-------------------------------------------------------------------------------------------------------------------------------------------------
#Soil elements definition
#-------------------------------------------------------------------------------------------------------------------------------------------------
element SSPbrickUP 1 68 59 116 126 54 36 103 112 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 2 86 80 140 148 68 59 116 126 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 3 113 121 173 167 86 80 140 148 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 4 143 165 196 188 113 121 173 167 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 5 189 197 209 206 203 207 212 211 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 6 183 189 206 205 198 203 211 210 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 7 177 183 205 202 195 198 210 208 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 8 69 60 117 125 87 79 142 147 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 9 53 35 104 111 69 60 117 125 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 10 37 26 98 105 53 35 104 111 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 11 130 89 94 101 179 150 154 160 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 12 130 101 114 123 179 160 168 176 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 13 99 83 93 102 156 144 153 161 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 14 88 93 153 149 130 102 161 179 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 15 85 88 78 72 146 149 135 132 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 16 72 77 89 85 132 134 150 146 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 17 130 155 164 171 179 191 193 200 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 18 171 164 170 181 200 193 199 204 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 19 136 123 124 139 185 176 178 187 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 20 139 162 155 136 187 192 191 185 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 21 101 94 82 97 160 154 145 157 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 22 97 110 114 101 157 166 168 160 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 23 69 87 147 125 83 99 156 144 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 24 83 78 88 93 144 135 149 153 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 25 69 53 78 83 125 111 135 144 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 26 53 37 72 78 111 105 132 135 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 27 37 54 77 72 105 112 134 132 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 28 77 82 94 89 134 145 154 150 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 29 68 54 112 126 82 77 134 145 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 30 86 68 126 148 97 82 145 157 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 31 86 113 110 97 148 167 166 157 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 32 110 124 123 114 166 178 176 168 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 33 143 113 167 188 124 110 166 178 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 34 177 143 188 202 139 124 178 187 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 35 177 183 162 139 202 205 192 187 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 36 162 170 164 155 192 199 193 191 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 37 189 183 205 206 170 162 192 199 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 38 197 189 206 209 181 170 199 204 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 39 89 85 146 150 130 88 149 179 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 40 155 136 185 191 130 123 176 179 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 41 54 36 103 112 37 26 98 105 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 42 165 143 188 196 195 177 202 208 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 43 22 17 59 68 14 8 36 54 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 44 50 43 80 86 22 17 59 68 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 45 73 84 121 113 50 43 80 86 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 46 109 131 165 143 73 84 121 113 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 47 169 182 197 189 190 201 207 203 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 48 159 169 189 183 184 190 203 198 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 49 152 159 183 177 180 184 198 195 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 50 23 18 60 69 49 42 79 87 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 51 13 9 35 53 23 18 60 69 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 52 10 5 26 37 13 9 35 53 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 53 96 51 58 67 130 89 94 101 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 54 96 67 74 90 130 101 114 123 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 55 62 47 57 66 99 83 93 102 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 56 52 57 93 88 96 66 102 130 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 57 48 52 38 34 85 88 78 72 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 58 34 39 51 48 72 77 89 85 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 59 96 118 129 138 130 155 164 171 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 60 138 129 137 158 171 164 170 181 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 61 106 90 91 108 136 123 124 139 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 62 108 122 118 106 139 162 155 136 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 63 67 58 46 61 101 94 82 97 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 64 61 71 74 67 97 110 114 101 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 65 23 49 87 69 47 62 99 83 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 66 47 38 52 57 83 78 88 93 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 67 23 13 38 47 69 53 78 83 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 68 13 10 34 38 53 37 72 78 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 69 10 14 39 34 37 54 77 72 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 70 39 46 58 51 77 82 94 89 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 71 22 14 54 68 46 39 77 82 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 72 50 22 68 86 61 46 82 97 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 73 50 73 71 61 86 113 110 97 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 74 71 91 90 74 110 124 123 114 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 75 109 73 113 143 91 71 110 124 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 76 152 109 143 177 108 91 124 139 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 77 152 159 122 108 177 183 162 139 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 78 122 137 129 118 162 170 164 155 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 79 169 159 183 189 137 122 162 170 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 80 182 169 189 197 158 137 170 181 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 81 51 48 85 89 96 52 88 130 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 82 118 106 136 155 96 90 123 130 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 83 14 8 36 54 10 5 26 37 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 84 131 109 143 165 180 152 177 195 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 85 15 11 17 22 7 2 8 14 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 86 31 24 43 50 15 11 17 22 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 87 64 70 84 73 31 24 43 50 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 88 100 120 131 109 64 70 84 73 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 89 163 174 182 169 186 194 201 190 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 90 151 163 169 159 175 186 190 184 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 91 133 151 159 152 172 175 184 180 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 92 16 12 18 23 30 25 42 49 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 93 6 3 9 13 16 12 18 23 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 94 4 1 5 10 6 3 9 13 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 95 81 33 40 56 96 51 58 67 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 96 81 56 65 75 96 67 74 90 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 97 45 28 41 55 62 47 57 66 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 98 32 41 57 52 81 55 66 96 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 99 29 32 20 19 48 52 38 34 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 100 19 21 33 29 34 39 51 48 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 101 81 107 119 128 96 118 129 138 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 102 128 119 127 141 138 129 137 158 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 103 92 75 76 95 106 90 91 108 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 104 95 115 107 92 108 122 118 106 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 105 56 40 27 44 67 58 46 61 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 106 44 63 65 56 61 71 74 67 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 107 16 30 49 23 28 45 62 47 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 108 28 20 32 41 47 38 52 57 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 109 16 6 20 28 23 13 38 47 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 110 6 4 19 20 13 10 34 38 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 111 4 7 21 19 10 14 39 34 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 112 21 27 40 33 39 46 58 51 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 113 15 7 14 22 27 21 39 46 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 114 31 15 22 50 44 27 46 61 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 115 31 64 63 44 50 73 71 61 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 116 63 76 75 65 71 91 90 74 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 117 100 64 73 109 76 63 71 91 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 118 133 100 109 152 95 76 91 108 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 119 133 151 115 95 152 159 122 108 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 120 115 127 119 107 122 137 129 118 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 121 163 151 159 169 127 115 122 137 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 122 174 163 169 182 141 127 137 158 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 123 33 29 48 51 81 32 52 96 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 124 107 92 106 118 81 75 90 96 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 125 7 2 8 14 4 1 5 10 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
element SSPbrickUP 126 120 100 109 131 172 133 152 180 1 $fBulk $fmass $perm $perm $perm $e_init $alpha $gravX $gravY $gravZ; # soil
set numel 3720
set SoilEle {}
for {set i 1} {$i <= $numel} {incr i 1} {
lappend SoilEle [expr $i]
}
puts " ## Soil is built!! ## "
#---------------------------------------------------------------------------------------------------------------------------------------------------
# Stage 1 - Soil self-weight
#---------------------------------------------------------------------------------------------------------------------------------------------------
numberer RCM
system SparseSPD
test NormDispIncr 1.00e-000 100 1
algorithm KrylovNewton
constraints Penalty 1.0e18 1.0e18
integrator Newmark 1.5 1.0
analysis Transient
#output2 definition
#---------------------------------------------------------------------------
eval "recorder Node -file output2/name.dg0E -time -node $SnodeList -dof 1 2 3 disp"
eval "recorder Node -file output2/name.pg0E -time -node $SnodeList -dof 4 vel"
eval "recorder Element -eleRange 1 $numel -time -file output2/name.s1g0E stress"
eval "recorder Element -eleRange 1 $numel -time -file output2/name.str1g0E strain"
#----------------------------------------------------------------------------
#Analysis running
#----------------------------------------------------------------------------
updateMaterialStage -material 1 -stage 0
updateMaterialStage -material 2 -stage 0
set NumIncr0 5
for {set numincr 1} {$numincr <= $NumIncr0} {incr numincr 1} {
puts "**Elastic Geostatic Analysis step : $numincr **"
analyze 1 10
}
puts " Elastic Geostatic Analysis step success ";
# switch material stage from elastic (gravity) to plastic ------- 2nd run
updateMaterialStage -material 1 -stage 1
updateMaterialStage -material 2 -stage 1
#output2 definition
#---------------------------------------------------------------------------
eval "recorder Node -file output2/name.dg0P -time -node $SnodeList -dof 1 2 3 disp"
eval "recorder Node -file output2/name.pg0P -time -node $SnodeList -dof 4 vel"
eval "recorder Element -eleRange 1 $numel -time -file output2/name.s1g0P stress"
eval "recorder Element -eleRange 1 $numel -time -file output2/name.str1g0P strain"
#----------------------------------------------------------------------------
set NumIncr1 5
for {set numincr 1} {$numincr <= $NumIncr1} {incr numincr 1} {
puts "**Plastic Geostatic Analysis step : $numincr **"
analyze 1 5
}
puts " Plastic Geostatic Analysis step success ";
puts "First run done. Gravity Applied"
puts "####################################"
# remove element
remove element 97
remove element 99
remove element 100
remove element 102
remove element 103
remove element 104
remove element 105
remove element 106
remove element 108
remove element 112
remove element 116
remove element 120
#-------------------------------------------------------------------------------
# Pile nodes creation
#-------------------------------------------------------------------------------
model BasicBuilder -ndm 3 -ndf 3
node 213 0.646447 0.353553 1.5
node 214 0.566987 0.25 1.5
node 215 0.75 0.433013 1.5
node 216 0.87059 0.482963 1.5
node 217 0.517037 0.12941 1.5
node 218 0.717157 0.282843 1.5
node 219 0.65359 0.2 1.5
node 220 0.8 0.34641 1.5
node 221 0.646447 0.353553 1
node 222 0.566987 0.25 1
node 223 0.75 0.433013 1
node 224 0.896472 0.38637 1.5
node 225 0.61363 0.103528 1.5
node 226 1 0.5 1.5
node 227 0.5 0 1.5
node 228 0.87059 0.482963 1
node 229 0.517037 0.12941 1
node 230 0.717157 0.282843 1
node 231 0.8 0.34641 1
node 232 0.65359 0.2 1
node 233 0.6 -1.9916E-17 1.5
node 234 1 0.4 1.5
node 235 0.61363 0.103528 1
node 236 0.896472 0.38637 1
node 237 1 0.5 1
node 238 0.5 0 1
node 239 1.12941 0.482963 1.5
node 240 1.10353 0.38637 1.5
node 241 0.6 -1.9916E-17 1
node 242 1 0.4 1
node 243 1.12941 0.482963 1
node 244 1.10353 0.38637 1
node 245 1.2 0.34641 1.5
node 246 1.25 0.433013 1.5
node 247 1.2 0.34641 1
node 248 1.25 0.433013 1
node 249 1.28284 0.282843 1.5
node 250 1.35355 0.353553 1.5
node 251 1.28284 0.282843 1
node 252 1.34641 0.2 1.5
node 253 1.35355 0.353553 1
node 254 1.43301 0.25 1.5
node 255 1.34641 0.2 1
node 256 1.38637 0.103528 1.5
node 257 1.43301 0.25 1
node 258 1.48296 0.12941 1.5
node 259 1.4 7.66952E-18 1.5
node 260 1.38637 0.103528 1
node 261 1.48296 0.12941 1
node 262 1.4 7.66952E-18 1
node 263 1.5 0 1.5
node 264 1.5 0 1
#------------------------------------------------------------------------------
#Pile material definition
nDMaterial ElasticIsotropic 200 2e+08 0.3 7.84
#Pile element creation
#-------------------------------------------------------------------------------
element SSPbrick 127 227 217 225 233 238 229 235 241 200 0 0 -77.0085
element SSPbrick 128 218 219 214 213 230 232 222 221 200 0 0 -77.0085
element SSPbrick 129 213 215 220 218 221 223 231 230 200 0 0 -77.0085
element SSPbrick 130 259 256 258 263 262 260 261 264 200 0 0 -77.0085
element SSPbrick 131 249 245 246 250 251 247 248 253 200 0 0 -77.0085
element SSPbrick 132 250 254 252 249 253 257 255 251 200 0 0 -77.0085
element SSPbrick 133 234 224 216 226 242 236 228 237 200 0 0 -77.0085
element SSPbrick 134 226 239 240 234 237 243 244 242 200 0 0 -77.0085
element SSPbrick 135 217 214 219 225 229 222 232 235 200 0 0 -77.0085
element SSPbrick 136 215 216 224 220 223 228 236 231 200 0 0 -77.0085
element SSPbrick 137 239 246 245 240 243 248 247 244 200 0 0 -77.0085
element SSPbrick 138 254 258 256 252 257 261 260 255 200 0 0 -77.0085
#---------------------------------------------------------------------------------
#Soil-monopile nodes equalDOF connection
#----------------------------------------------------------------------------------
equalDOF 19 213 1 2 3
equalDOF 20 214 1 2 3
equalDOF 21 215 1 2 3
equalDOF 27 216 1 2 3
equalDOF 28 217 1 2 3
equalDOF 29 218 1 2 3
equalDOF 32 219 1 2 3
equalDOF 33 220 1 2 3
equalDOF 34 221 1 2 3
equalDOF 38 222 1 2 3
equalDOF 39 223 1 2 3
equalDOF 40 224 1 2 3
equalDOF 41 225 1 2 3
equalDOF 44 226 1 2 3
equalDOF 45 227 1 2 3
equalDOF 46 228 1 2 3
equalDOF 47 229 1 2 3
equalDOF 48 230 1 2 3
equalDOF 51 231 1 2 3
equalDOF 52 232 1 2 3
equalDOF 55 233 1 2 3
equalDOF 56 234 1 2 3
equalDOF 57 235 1 2 3
equalDOF 58 236 1 2 3
equalDOF 61 237 1 2 3
equalDOF 62 238 1 2 3
equalDOF 63 239 1 2 3
equalDOF 65 240 1 2 3
equalDOF 66 241 1 2 3
equalDOF 67 242 1 2 3
equalDOF 71 243 1 2 3
equalDOF 74 244 1 2 3
equalDOF 75 245 1 2 3
equalDOF 76 246 1 2 3
equalDOF 90 247 1 2 3
equalDOF 91 248 1 2 3
equalDOF 92 249 1 2 3
equalDOF 95 250 1 2 3
equalDOF 106 251 1 2 3
equalDOF 107 252 1 2 3
equalDOF 108 253 1 2 3
equalDOF 115 254 1 2 3
equalDOF 118 255 1 2 3
equalDOF 119 256 1 2 3
equalDOF 122 257 1 2 3
equalDOF 127 258 1 2 3
equalDOF 128 259 1 2 3
equalDOF 129 260 1 2 3
equalDOF 137 261 1 2 3
equalDOF 138 262 1 2 3
equalDOF 141 263 1 2 3
equalDOF 158 264 1 2 3
loadConst -time 0.0
set numSteps2ndRun 5
#Analysis running
#---------------------------------------------------------------------------------
for {set numincr 1} {$numincr <= $numSteps2ndRun} {incr numincr 1} {
puts "**Structure Installed step : $numincr **"
analyze 1 1
}
puts "Second run done. Structure Installed"
puts "##########################################"