Cylinder modeling using shell elements

[hamidreza]

Dear Moderators
I modeled a cylinder with 32 m diameter and 46 m height with shell elements and applied time history loading on it. But the model does not converge with the real values for E and ro(density). When I use E=2.486E10 N/m and ro= 2400 (Kg/M^3) the model diverges in 4 second. While if I change the (Density) ro=2.4 or increase the E the model would converge.
I would be thankful if you could take a look at model and tell me the problem.
Best Regards,

wipe all ;
# Definition of model------------------------------------------
model basic -ndm 3 -ndf 6
set H 46
set h [expr $H/10]
set dataDir edc-0; # set up name of data directory
file mkdir $dataDir;
# Definition of node-------------------------------------------
#BOTTOM DIAPHRAGM
#**********************
###OUTER CIRCLE
node 56 16 0 0
node 57 15.035 5.47 0
node 58 12.26 10.28 0
node 59 8 13.86 0
node 60 2.78 15.76 0
node 61 -2.78 15.76 0
node 62 -8 13.86 0
node 63 -12.26 10.28 0
node 64 -15.035 5.47 0
node 65 -16 0 0
node 66 -15.035 -5.47 0
node 67 -12.26 -10.28 0
node 68 -8 -13.86 0
node 69 -2.78 -15.76 0
node 70 2.78 -15.76 0
node 71 8 -13.86 0
node 72 12.26 -10.28 0
node 73 15.035 -5.47 0
################################
##############################################################
node 1056 16 0 $h
node 1057 15.035 5.47 $h
node 1058 12.26 10.28 $h
node 1059 8 13.86 $h
node 1060 2.78 15.76 $h
node 1061 -2.78 15.76 $h
node 1062 -8 13.86 $h
node 1063 -12.26 10.28 $h
node 1064 -15.035 5.47 $h
node 1065 -16 0 $h
node 1066 -15.035 -5.47 $h
node 1067 -12.26 -10.28 $h
node 1068 -8 -13.86 $h
node 1069 -2.78 -15.76 $h
node 1070 2.78 -15.76 $h
node 1071 8 -13.86 $h
node 1072 12.26 -10.28 $h
node 1073 15.035 -5.47 $h
##############################################################
node 2056 16 0 [expr 2*$h]
node 2057 15.035 5.47 [expr 2*$h]
node 2058 12.26 10.28 [expr 2*$h]
node 2059 8 13.86 [expr 2*$h]
node 2060 2.78 15.76 [expr 2*$h]
node 2061 -2.78 15.76 [expr 2*$h]
node 2062 -8 13.86 [expr 2*$h]
node 2063 -12.26 10.28 [expr 2*$h]
node 2064 -15.035 5.47 [expr 2*$h]
node 2065 -16 0 [expr 2*$h]
node 2066 -15.035 -5.47 [expr 2*$h]
node 2067 -12.26 -10.28 [expr 2*$h]
node 2068 -8 -13.86 [expr 2*$h]
node 2069 -2.78 -15.76 [expr 2*$h]
node 2070 2.78 -15.76 [expr 2*$h]
node 2071 8 -13.86 [expr 2*$h]
node 2072 12.26 -10.28 [expr 2*$h]
node 2073 15.035 -5.47 [expr 2*$h]
##############################################################
node 3056 16 0 [expr 3*$h]
node 3057 15.035 5.47 [expr 3*$h]
node 3058 12.26 10.28 [expr 3*$h]
node 3059 8 13.86 [expr 3*$h]
node 3060 2.78 15.76 [expr 3*$h]
node 3061 -2.78 15.76 [expr 3*$h]
node 3062 -8 13.86 [expr 3*$h]
node 3063 -12.26 10.28 [expr 3*$h]
node 3064 -15.035 5.47 [expr 3*$h]
node 3065 -16 0 [expr 3*$h]
node 3066 -15.035 -5.47 [expr 3*$h]
node 3067 -12.26 -10.28 [expr 3*$h]
node 3068 -8 -13.86 [expr 3*$h]
node 3069 -2.78 -15.76 [expr 3*$h]
node 3070 2.78 -15.76 [expr 3*$h]
node 3071 8 -13.86 [expr 3*$h]
node 3072 12.26 -10.28 [expr 3*$h]
node 3073 15.035 -5.47 [expr 3*$h]
###############################################################
node 4056 16 0 [expr 4*$h]
node 4057 15.035 5.47 [expr 4*$h]
node 4058 12.26 10.28 [expr 4*$h]
node 4059 8 13.86 [expr 4*$h]
node 4060 2.78 15.76 [expr 4*$h]
node 4061 -2.78 15.76 [expr 4*$h]
node 4062 -8 13.86 [expr 4*$h]
node 4063 -12.26 10.28 [expr 4*$h]
node 4064 -15.035 5.47 [expr 4*$h]
node 4065 -16 0 [expr 4*$h]
node 4066 -15.035 -5.47 [expr 4*$h]
node 4067 -12.26 -10.28 [expr 4*$h]
node 4068 -8 -13.86 [expr 4*$h]
node 4069 -2.78 -15.76 [expr 4*$h]
node 4070 2.78 -15.76 [expr 4*$h]
node 4071 8 -13.86 [expr 4*$h]
node 4072 12.26 -10.28 [expr 4*$h]
node 4073 15.035 -5.47 [expr 4*$h]
##############################################################
node 5056 16 0 [expr 5*$h]
node 5057 15.035 5.47 [expr 5*$h]
node 5058 12.26 10.28 [expr 5*$h]
node 5059 8 13.86 [expr 5*$h]
node 5060 2.78 15.76 [expr 5*$h]
node 5061 -2.78 15.76 [expr 5*$h]
node 5062 -8 13.86 [expr 5*$h]
node 5063 -12.26 10.28 [expr 5*$h]
node 5064 -15.035 5.47 [expr 5*$h]
node 5065 -16 0 [expr 5*$h]
node 5066 -15.035 -5.47 [expr 5*$h]
node 5067 -12.26 -10.28 [expr 5*$h]
node 5068 -8 -13.86 [expr 5*$h]
node 5069 -2.78 -15.76 [expr 5*$h]
node 5070 2.78 -15.76 [expr 5*$h]
node 5071 8 -13.86 [expr 5*$h]
node 5072 12.26 -10.28 [expr 5*$h]
node 5073 15.035 -5.47 [expr 5*$h]
#############################################################
node 6056 16 0 [expr 6*$h]
node 6057 15.035 5.47 [expr 6*$h]
node 6058 12.26 10.28 [expr 6*$h]
node 6059 8 13.86 [expr 6*$h]
node 6060 2.78 15.76 [expr 6*$h]
node 6061 -2.78 15.76 [expr 6*$h]
node 6062 -8 13.86 [expr 6*$h]
node 6063 -12.26 10.28 [expr 6*$h]
node 6064 -15.035 5.47 [expr 6*$h]
node 6065 -16 0 [expr 6*$h]
node 6066 -15.035 -5.47 [expr 6*$h]
node 6067 -12.26 -10.28 [expr 6*$h]
node 6068 -8 -13.86 [expr 6*$h]
node 6069 -2.78 -15.76 [expr 6*$h]
node 6070 2.78 -15.76 [expr 6*$h]
node 6071 8 -13.86 [expr 6*$h]
node 6072 12.26 -10.28 [expr 6*$h]
node 6073 15.035 -5.47 [expr 6*$h]
###############################################################
node 7056 16 0 [expr 7*$h]
node 7057 15.035 5.47 [expr 7*$h]
node 7058 12.26 10.28 [expr 7*$h]
node 7059 8 13.86 [expr 7*$h]
node 7060 2.78 15.76 [expr 7*$h]
node 7061 -2.78 15.76 [expr 7*$h]
node 7062 -8 13.86 [expr 7*$h]
node 7063 -12.26 10.28 [expr 7*$h]
node 7064 -15.035 5.47 [expr 7*$h]
node 7065 -16 0 [expr 7*$h]
node 7066 -15.035 -5.47 [expr 7*$h]
node 7067 -12.26 -10.28 [expr 7*$h]
node 7068 -8 -13.86 [expr 7*$h]
node 7069 -2.78 -15.76 [expr 7*$h]
node 7070 2.78 -15.76 [expr 7*$h]
node 7071 8 -13.86 [expr 7*$h]
node 7072 12.26 -10.28 [expr 7*$h]
node 7073 15.035 -5.47 [expr 7*$h]
##############################################################
node 8056 16 0 [expr 8*$h]
node 8057 15.035 5.47 [expr 8*$h]
node 8058 12.26 10.28 [expr 8*$h]
node 8059 8 13.86 [expr 8*$h]
node 8060 2.78 15.76 [expr 8*$h]
node 8061 -2.78 15.76 [expr 8*$h]
node 8062 -8 13.86 [expr 8*$h]
node 8063 -12.26 10.28 [expr 8*$h]
node 8064 -15.035 5.47 [expr 8*$h]
node 8065 -16 0 [expr 8*$h]
node 8066 -15.035 -5.47 [expr 8*$h]
node 8067 -12.26 -10.28 [expr 8*$h]
node 8068 -8 -13.86 [expr 8*$h]
node 8069 -2.78 -15.76 [expr 8*$h]
node 8070 2.78 -15.76 [expr 8*$h]
node 8071 8 -13.86 [expr 8*$h]
node 8072 12.26 -10.28 [expr 8*$h]
node 8073 15.035 -5.47 [expr 8*$h]
###############################################################
node 9056 16 0 [expr 9*$h]
node 9057 15.035 5.47 [expr 9*$h]
node 9058 12.26 10.28 [expr 9*$h]
node 9059 8 13.86 [expr 9*$h]
node 9060 2.78 15.76 [expr 9*$h]
node 9061 -2.78 15.76 [expr 9*$h]
node 9062 -8 13.86 [expr 9*$h]
node 9063 -12.26 10.28 [expr 9*$h]
node 9064 -15.035 5.47 [expr 9*$h]
node 9065 -16 0 [expr 9*$h]
node 9066 -15.035 -5.47 [expr 9*$h]
node 9067 -12.26 -10.28 [expr 9*$h]
node 9068 -8 -13.86 [expr 9*$h]
node 9069 -2.78 -15.76 [expr 9*$h]
node 9070 2.78 -15.76 [expr 9*$h]
node 9071 8 -13.86 [expr 9*$h]
node 9072 12.26 -10.28 [expr 9*$h]
node 9073 15.035 -5.47 [expr 9*$h]


#OUTER CIRCLE
node 10056 16 0 $H
node 10057 15.035 5.47 $H
node 10058 12.26 10.28 $H
node 10059 8 13.86 $H
node 10060 2.78 15.76 $H
node 10061 -2.78 15.76 $H
node 10062 -8 13.86 $H
node 10063 -12.26 10.28 $H
node 10064 -15.035 5.47 $H
node 10065 -16 0 $H
node 10066 -15.035 -5.47 $H
node 10067 -12.26 -10.28 $H
node 10068 -8 -13.86 $H
node 10069 -2.78 -15.76 $H
node 10070 2.78 -15.76 $H
node 10071 8 -13.86 $H
node 10072 12.26 -10.28 $H
node 10073 15.035 -5.47 $H
##############################################################
rigidDiaphragm 3 10056 10057 10058 10059 10060 10061 10062 10063 10064 10065 10066 10067 10068 10069 10070 10071 10072 10073
######################################################

##############################
#OUTER CIRCLE IS FIXED***
fix 56 1 1 1 1 1 1
fix 57 1 1 1 1 1 1
fix 58 1 1 1 1 1 1
fix 59 1 1 1 1 1 1
fix 60 1 1 1 1 1 1
fix 61 1 1 1 1 1 1
fix 62 1 1 1 1 1 1
fix 63 1 1 1 1 1 1
fix 64 1 1 1 1 1 1
fix 65 1 1 1 1 1 1
fix 66 1 1 1 1 1 1
fix 67 1 1 1 1 1 1
fix 68 1 1 1 1 1 1
fix 69 1 1 1 1 1 1
fix 70 1 1 1 1 1 1
fix 71 1 1 1 1 1 1
fix 72 1 1 1 1 1 1
fix 73 1 1 1 1 1 1

section ElasticMembranePlateSection 1 2.486E10 0.2 0.5 2.4

#########################################################_____________
## Mass lumped at master node(N/m/s^2)
set pi 3.14159
set g 9.806

# Definition of elements--------------------------------------------------
# Definition of shell

element ShellMITC4 1 56 1056 1057 57 1
element ShellMITC4 2 57 1057 1058 58 1
element ShellMITC4 3 58 1058 1059 59 1
element ShellMITC4 4 59 1059 1060 60 1
element ShellMITC4 5 60 1060 1061 61 1
element ShellMITC4 6 61 1061 1062 62 1
element ShellMITC4 7 62 1062 1063 63 1
element ShellMITC4 8 63 1063 1064 64 1
element ShellMITC4 9 64 1064 1065 65 1
element ShellMITC4 10 65 1065 1066 66 1
element ShellMITC4 11 66 1066 1067 67 1
element ShellMITC4 12 67 1067 1068 68 1
element ShellMITC4 13 68 1068 1069 69 1
element ShellMITC4 14 69 1069 1070 70 1
element ShellMITC4 15 70 1070 1071 71 1
element ShellMITC4 16 71 1071 1072 72 1
element ShellMITC4 17 72 1072 1073 73 1
element ShellMITC4 18 73 1073 1056 56 1

element ShellMITC4 19 1056 2056 2057 1057 1
element ShellMITC4 20 1057 2057 2058 1058 1
element ShellMITC4 21 1058 2058 2059 1059 1
element ShellMITC4 22 1059 2059 2060 1060 1
element ShellMITC4 23 1060 2060 2061 1061 1
element ShellMITC4 24 1061 2061 2062 1062 1
element ShellMITC4 25 1062 2062 2063 1063 1
element ShellMITC4 26 1063 2063 2064 1064 1
element ShellMITC4 27 1064 2064 2065 1065 1
element ShellMITC4 28 1065 2065 2066 1066 1
element ShellMITC4 29 1066 2066 2067 1067 1
element ShellMITC4 30 1067 2067 2068 1068 1
element ShellMITC4 31 1068 2068 2069 1069 1
element ShellMITC4 32 1069 2069 2070 1070 1
element ShellMITC4 33 1070 2070 2071 1071 1
element ShellMITC4 34 1071 2071 2072 1072 1
element ShellMITC4 35 1072 2072 2073 1073 1
element ShellMITC4 36 1073 2073 2056 1056 1
#
element ShellMITC4 37 2056 3056 3057 2057 1
element ShellMITC4 38 2057 3057 3058 2058 1
element ShellMITC4 39 2058 3058 3059 2059 1
element ShellMITC4 40 2059 3059 3060 2060 1
element ShellMITC4 41 2060 3060 3061 2061 1
element ShellMITC4 42 2061 3061 3062 2062 1
element ShellMITC4 43 2062 3062 3063 2063 1
element ShellMITC4 44 2063 3063 3064 2064 1
element ShellMITC4 45 2064 3064 3065 2065 1
element ShellMITC4 46 2065 3065 3066 2066 1
element ShellMITC4 47 2066 3066 3067 2067 1
element ShellMITC4 48 2067 3067 3068 2068 1
element ShellMITC4 49 2068 3068 3069 2069 1
element ShellMITC4 50 2069 3069 3070 2070 1
element ShellMITC4 51 2070 3070 3071 2071 1
element ShellMITC4 52 2071 3071 3072 2072 1
element ShellMITC4 53 2072 3072 3073 2073 1
element ShellMITC4 54 2073 3073 3056 2056 1
#
element ShellMITC4 55 3056 4056 4057 3057 1
element ShellMITC4 56 3057 4057 4058 3058 1
element ShellMITC4 57 3058 4058 4059 3059 1
element ShellMITC4 58 3059 4059 4060 3060 1
element ShellMITC4 59 3060 4060 4061 3061 1
element ShellMITC4 60 3061 4061 4062 3062 1
element ShellMITC4 61 3062 4062 4063 3063 1
element ShellMITC4 62 3063 4063 4064 3064 1
element ShellMITC4 63 3064 4064 4065 3065 1
element ShellMITC4 64 3065 4065 4066 3066 1
element ShellMITC4 65 3066 4066 4067 3067 1
element ShellMITC4 66 3067 4067 4068 3068 1
element ShellMITC4 67 3068 4068 4069 3069 1
element ShellMITC4 68 3069 4069 4070 3070 1
element ShellMITC4 69 3070 4070 4071 3071 1
element ShellMITC4 70 3071 4071 4072 3072 1
element ShellMITC4 71 3072 4072 4073 3073 1
element ShellMITC4 72 3073 4073 4056 3056 1
#
element ShellMITC4 73 4056 5056 5057 4057 1
element ShellMITC4 74 4057 5057 5058 4058 1
element ShellMITC4 75 4058 5058 5059 4059 1
element ShellMITC4 76 4059 5059 5060 4060 1
element ShellMITC4 77 4060 5060 5061 4061 1
element ShellMITC4 78 4061 5061 5062 4062 1
element ShellMITC4 79 4062 5062 5063 4063 1
element ShellMITC4 80 4063 5063 5064 4064 1
element ShellMITC4 81 4064 5064 5065 4065 1
element ShellMITC4 82 4065 5065 5066 4066 1
element ShellMITC4 83 4066 5066 5067 4067 1
element ShellMITC4 84 4067 5067 5068 4068 1
element ShellMITC4 85 4068 5068 5069 4069 1
element ShellMITC4 86 4069 5069 5070 4070 1
element ShellMITC4 87 4070 5070 5071 4071 1
element ShellMITC4 88 4071 5071 5072 4072 1
element ShellMITC4 89 4072 5072 5073 4073 1
element ShellMITC4 90 4073 5073 5056 4056 1
#
element ShellMITC4 91 5056 6056 6057 5057 1
element ShellMITC4 92 5057 6057 6058 5058 1
element ShellMITC4 93 5058 6058 6059 5059 1
element ShellMITC4 94 5059 6059 6060 5060 1
element ShellMITC4 95 5060 6060 6061 5061 1
element ShellMITC4 96 5061 6061 6062 5062 1
element ShellMITC4 97 5062 6062 6063 5063 1
element ShellMITC4 98 5063 6063 6064 5064 1
element ShellMITC4 99 5064 6064 6065 5065 1
element ShellMITC4 100 5065 6065 6066 5066 1
element ShellMITC4 101 5066 6066 6067 5067 1
element ShellMITC4 102 5067 6067 6068 5068 1
element ShellMITC4 103 5068 6068 6069 5069 1
element ShellMITC4 104 5069 6069 6070 5070 1
element ShellMITC4 105 5070 6070 6071 5071 1
element ShellMITC4 106 5071 6071 6072 5072 1
element ShellMITC4 107 5072 6072 6073 5073 1
element ShellMITC4 108 5073 6073 6056 5056 1
#
element ShellMITC4 109 6056 7056 7057 6057 1
element ShellMITC4 110 6057 7057 7058 6058 1
element ShellMITC4 111 6058 7058 7059 6059 1
element ShellMITC4 112 6059 7059 7060 6060 1
element ShellMITC4 113 6060 7060 7061 6061 1
element ShellMITC4 114 6061 7061 7062 6062 1
element ShellMITC4 115 6062 7062 7063 6063 1
element ShellMITC4 116 6063 7063 7064 6064 1
element ShellMITC4 117 6064 7064 7065 6065 1
element ShellMITC4 118 6065 7065 7066 6066 1
element ShellMITC4 119 6066 7066 7067 6067 1
element ShellMITC4 120 6067 7067 7068 6068 1
element ShellMITC4 121 6068 7068 7069 6069 1
element ShellMITC4 122 6069 7069 7070 6070 1
element ShellMITC4 123 6070 7070 7071 6071 1
element ShellMITC4 124 6071 7071 7072 6072 1
element ShellMITC4 125 6072 7072 7073 6073 1
element ShellMITC4 126 6073 7073 7056 6056 1
#
element ShellMITC4 127 7056 8056 8057 7057 1
element ShellMITC4 128 7057 8057 8058 7058 1
element ShellMITC4 129 7058 8058 8059 7059 1
element ShellMITC4 130 7059 8059 8060 7060 1
element ShellMITC4 131 7060 8060 8061 7061 1
element ShellMITC4 132 7061 8061 8062 7062 1
element ShellMITC4 133 7062 8062 8063 7063 1
element ShellMITC4 134 7063 8063 8064 7064 1
element ShellMITC4 135 7064 8064 8065 7065 1
element ShellMITC4 136 7065 8065 8066 7066 1
element ShellMITC4 137 7066 8066 8067 7067 1
element ShellMITC4 138 7067 8067 8068 7068 1
element ShellMITC4 139 7068 8068 8069 7069 1
element ShellMITC4 140 7069 8069 8070 7070 1
element ShellMITC4 141 7070 8070 8071 7071 1
element ShellMITC4 142 7071 8071 8072 7072 1
element ShellMITC4 143 7072 8072 8073 7073 1
element ShellMITC4 144 7073 8073 8056 7056 1
#
element ShellMITC4 145 8056 9056 9057 8057 1
element ShellMITC4 146 8057 9057 9058 8058 1
element ShellMITC4 147 8058 9058 9059 8059 1
element ShellMITC4 148 8059 9059 9060 8060 1
element ShellMITC4 149 8060 9060 9061 8061 1
element ShellMITC4 150 8061 9061 9062 8062 1
element ShellMITC4 151 8062 9062 9063 8063 1
element ShellMITC4 152 8063 9063 9064 8064 1
element ShellMITC4 153 8064 9064 9065 8065 1
element ShellMITC4 154 8065 9065 9066 8066 1
element ShellMITC4 155 8066 9066 9067 8067 1
element ShellMITC4 156 8067 9067 9068 8068 1
element ShellMITC4 157 8068 9068 9069 8069 1
element ShellMITC4 158 8069 9069 9070 8070 1
element ShellMITC4 159 8070 9070 9071 8071 1
element ShellMITC4 160 8071 9071 9072 8072 1
element ShellMITC4 161 8072 9072 9073 8073 1
element ShellMITC4 162 8073 9073 9056 8056 1
#
element ShellMITC4 163 9056 10056 10057 9057 1
element ShellMITC4 164 9057 10057 10058 9058 1
element ShellMITC4 165 9058 10058 10059 9059 1
element ShellMITC4 166 9059 10059 10060 9060 1
element ShellMITC4 167 9060 10060 10061 9061 1
element ShellMITC4 168 9061 10061 10062 9062 1
element ShellMITC4 169 9062 10062 10063 9063 1
element ShellMITC4 170 9063 10063 10064 9064 1
element ShellMITC4 171 9064 10064 10065 9065 1
element ShellMITC4 172 9065 10065 10066 9066 1
element ShellMITC4 173 9066 10066 10067 9067 1
element ShellMITC4 174 9067 10067 10068 9068 1
element ShellMITC4 175 9068 10068 10069 9069 1
element ShellMITC4 176 9069 10069 10070 9070 1
element ShellMITC4 177 9070 10070 10071 9071 1
element ShellMITC4 178 9071 10071 10072 9072 1
element ShellMITC4 179 9072 10072 10073 9073 1
element ShellMITC4 180 9073 10073 10056 9056 1


#Eigenvalue analysis----------------------------------------------------------------------------------
set numModes 5
set lambda [eigen $numModes];
# calculate frequencies and periods of the structure
set omega {}
set f {}
set T {}
set pi 3.141593

foreach lam $lambda {
lappend omega [expr sqrt($lam)]
lappend f [expr sqrt($lam)/(2*$pi)]
lappend T [expr (2*$pi)/sqrt($lam)]
}
puts "periods are $T"

#Create the recorder objects-------------------------------------------------------------------------------
recorder Node -file "top_dsp.txt " -time -node 10056 -dof 1 disp

# Display of Model3D --------------------------------------------------------------------------------
recorder display animation 10 10 600 600 -file 1
prp 20 90 60
vup 0 0 2
fill 4
display 10 4 1000

# Definition of analysis-----------------------------------------------------------------------------------
set north "Path -filePath Elcentro.txt -dt 0.02 -factor $g"

rayleigh 0.95 0.00224 0 0 ;


pattern UniformExcitation 2 1 -accel $north
test EnergyIncr 1.0e-4 200 3
algorithm Newton
system SparseGeneral
constraints Transformation
integrator Newmark 0.5 0.25
numberer RCM
analysis Transient

analyze 2000 .01

puts "Analysis is done"
exit

[vesna]

Try doing the analysis in the loop like in this example: http://opensees.berkeley.edu/wiki/index ... e_Analysis

[hamidreza]

Dear Vesna

I applied the changes you asked me, but unfortunately It did not run completely. Apparently it has other problem with shell element, or if its problem is related to analysis, it also has the problem by the changes you suggested me..
I would be thankful if you could go through the model once again and see the problem.

I applied following changes in analysis:

#set MotionInterval 0.005
##set MotionFile "Elcentro.txt"
##set MotionScale 0.01
#set MotionStep 4000
#set TmaxAnalysis 20
#set Tol 0.001
#constraints Transformation
#numberer RCM
#system SparseGeneral
#test EnergyIncr 1.0e-4 200
#algorithm Newton
#integrator Newmark 0.5 0.25
#analysis Transient
#set Nsteps $MotionStep ;
#set ok [analyze $Nsteps $MotionInterval];
#if {$ok != 0} {
#set ok 0;
#set controlTime [getTime];
#while {$controlTime < $TmaxAnalysis && $ok == 0} {
#set controlTime [getTime]
#set ok [analyze 1 $MotionInterval]
#if {$ok != 0} {
#puts "Trying Newton with Initial Tangent .."
#test NormDispIncr $Tol 1000 0
#algorithm Newton -initial
#set ok [analyze 1 $MotionInterval]
#test EnergyIncr 1.0e-3 500 3
##test $testTypeDynamic $TolDynamic $maxNumIterDynamic 0
#algorithm ModifiedNewton
#}
#if {$ok != 0} {
#puts "Trying Broyden .."
#algorithm Broyden 8
#set ok [analyze 1 $MotionInterval]
#algorithm NewtonLineSearch 0.6
#}
#if {$ok != 0} {
#puts "Trying NewtonWithLineSearch .."
#algorithm NewtonLineSearch .8
#set ok [analyze 1 $MotionInterval]
#algorithm KrylovNewton
#}
#}
#}; # end if ok !0