SSPbrickUP instead of brickUP can't run correctly

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

Moderators: silvia, selimgunay, Moderators

Post Reply
xiaoweimeng
Posts: 4
Joined: Wed Sep 30, 2020 12:27 am

SSPbrickUP instead of brickUP can't run correctly

Post by xiaoweimeng »

I can't run the following code ,but when I change the SSPbrickUP into brickUP,I can run correctly.who can tell me the reason?

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 "##########################################"
Post Reply