#DEFINE PROJECT = Circle-Benchmark
+PROG AQUA urs:1
HEAD Materials
PAGE UNII 0
ECHO FULL no
ECHO MAT FULL

MATE 1 E 10000[mpa]
MATE 2 E 54.722432[mpa]
MATE 3 E 137.278634[mpa]
MATE 4 E 250.696152[mpa]
MATE 5 E 487.131105[mpa]
MATE 6 E 610.285625[mpa]
MATE 7 E 838.110375[mpa]
MATE 8 E 1254.437394[mpa]
MATE 9 E 1647.924909[mpa]
MATE 10 E 2469.182279[mpa]
MATE 11 E 1766.077907[mpa]
MATE 12 E 1011.648779[mpa]
MATE 13 E 682.861546[mpa]
MATE 14 E 422.096887[mpa]
MATE 15 E 354.104102[mpa]
MATE 16 E 314.965581[mpa]
MATE 17 E 260.892805[mpa]
MATE 18 E 206.465406[mpa]
MATE 19 E 93.364139[mpa]
END
+PROG SOFIMSHC urs:3
HEAD
PAGE UNII 0
SYST 3D GDIR NEGZ GDIV -1000
CTRL INIT 7
CTRL TOLG 0.0001
CTRL MESH 1
CTRL HMIN 0.05
CTRL TOPO XFLG 1
CTRL NODE 1000
Spt No 101 -0.947076819681804 0.168173525797887 0
Spt No 102 -0.862158128481804 0.168173525797887 0
Spt No 103 -0.831989646081804 0.268173525797887 0
Spt No 104 -0.947076819681804 0.268173525797887 0
Spt No 105 -0.862158128481804 0.168173525797887 0
Spt No 106 -0.747549819681804 0.168173525797887 0
Spt No 107 -0.747549819681804 0.268173525797887 0
Spt No 108 -0.831989646081804 0.268173525797887 0
Spt No 109 -0.662870560881804 0.168173525797887 0
Spt No 110 -0.632702078481804 0.268173525797887 0
Spt No 111 -0.662870560881804 0.168173525797887 0
Spt No 112 -0.548022819681804 0.168173525797887 0
Spt No 113 -0.548022819681804 0.268173525797887 0
Spt No 114 -0.632702078481804 0.268173525797887 0
Spt No 115 -0.463343560881804 0.168173525797887 0
Spt No 116 -0.433175078481804 0.268173525797887 0
Spt No 117 -0.463343560881804 0.168173525797887 0
Spt No 118 -0.348495819681804 0.168173525797887 0
Spt No 119 -0.348495819681804 0.268173525797887 0
Spt No 120 -0.433175078481804 0.268173525797887 0
Spt No 121 -0.264175709481804 0.168173525797887 0
Spt No 122 -0.233648078481804 0.268173525797887 0
Spt No 123 -0.264175709481804 0.168173525797887 0
Spt No 124 -0.148968819681804 0.168173525797887 0
Spt No 125 -0.148968819681804 0.268173525797887 0
Spt No 126 -0.233648078481804 0.268173525797887 0
Spt No 127 -0.0642895608818039 0.168173525797887 0
Spt No 128 -0.0341210784818041 0.268173525797887 0
Spt No 129 -0.0642895608818039 0.168173525797887 0
Spt No 130 0.0505581803181961 0.168173525797887 0
Spt No 131 0.0505581803181961 0.268173525797887 0
Spt No 132 -0.0341210784818041 0.268173525797887 0
Spt No 133 0.135237439118196 0.168173525797887 0
Spt No 134 0.165405921518196 0.268173525797887 0
Spt No 135 0.135237439118196 0.168173525797887 0
Spt No 136 0.250085180318196 0.168173525797887 0
Spt No 137 0.250085180318196 0.268173525797887 0
Spt No 138 0.165405921518196 0.268173525797887 0
Spt No 139 0.334764439118196 0.168173525797887 0
Spt No 140 0.364932921518196 0.268173525797887 0
Spt No 141 0.334764439118196 0.168173525797887 0
Spt No 142 0.449612180318196 0.168173525797887 0
Spt No 143 0.449612180318196 0.268173525797887 0
Spt No 144 0.364932921518196 0.268173525797887 0
Spt No 145 0.534291439118196 0.168173525797887 0
Spt No 146 0.564100772918196 0.268173525797887 0
Spt No 147 0.534291439118196 0.168173525797887 0
Spt No 148 0.649139180318196 0.168173525797887 0
Spt No 149 0.649139180318196 0.268173525797887 0
Spt No 150 0.564100772918196 0.268173525797887 0
Spt No 151 0.733459290518196 0.168173525797887 0
Spt No 152 0.763627772918196 0.268173525797887 0
Spt No 153 0.733459290518196 0.168173525797887 0
Spt No 154 0.848666180318196 0.168173525797887 0
Spt No 155 0.848666180318196 0.268173525797887 0
Spt No 156 0.763627772918196 0.268173525797887 0
Spt No 157 0.933704587718196 0.168173525797887 0
Spt No 158 0.963873070118197 0.268173525797887 0
Spt No 159 0.933704587718196 0.168173525797887 0
Spt No 160 1.0481931803182 0.168173525797887 0
Spt No 161 1.0481931803182 0.268173525797887 0
Spt No 162 0.963873070118197 0.268173525797887 0
Spt No 163 1.1332315877182 0.168173525797887 0
Spt No 164 1.1630409215182 0.268173525797887 0
Spt No 165 1.1332315877182 0.168173525797887 0
Spt No 166 1.2477201803182 0.168173525797887 0
Spt No 167 1.2477201803182 0.268173525797887 0
Spt No 168 1.1630409215182 0.268173525797887 0
Spt No 169 1.3327585877182 0.168173525797887 0
Spt No 170 1.3625679215182 0.268173525797887 0
Spt No 171 1.3327585877182 0.168173525797887 0
Spt No 172 1.4472471803182 0.168173525797887 0
Spt No 173 1.4472471803182 0.268173525797887 0
Spt No 174 1.3625679215182 0.268173525797887 0
Spt No 175 1.5319264391182 0.168173525797887 0
Spt No 176 1.5620949215182 0.268173525797887 0
Spt No 177 1.5319264391182 0.168173525797887 0
Spt No 178 1.6467741803182 0.168173525797887 0
Spt No 179 1.6467741803182 0.268173525797887 0
Spt No 180 1.5620949215182 0.268173525797887 0
Spt No 181 1.7314534391182 0.168173525797887 0
Spt No 182 1.7612627729182 0.268173525797887 0
Spt No 183 1.7314534391182 0.168173525797887 0
Spt No 184 1.8463011803182 0.168173525797887 0
Spt No 185 1.8463011803182 0.268173525797887 0
Spt No 186 1.7612627729182 0.268173525797887 0
Spt No 187 1.9309804391182 0.168173525797887 0
Spt No 188 1.9611489215182 0.268173525797887 0
Spt No 189 1.9309804391182 0.168173525797887 0
Spt No 190 2.0458281803182 0.168173525797887 0
Spt No 191 2.0458281803182 0.268173525797887 0
Spt No 192 1.9611489215182 0.268173525797887 0
Spt No 193 2.1305074391182 0.168173525797887 0
Spt No 194 2.1603167729182 0.268173525797887 0
Spt No 195 2.1305074391182 0.168173525797887 0
Spt No 196 2.2453551803182 0.168173525797887 0
Spt No 197 2.2453551803182 0.268173525797887 0
Spt No 198 2.1603167729182 0.268173525797887 0
Spt No 199 2.3300344391182 0.168173525797887 0
Spt No 200 2.3598437729182 0.268173525797887 0
Spt No 201 2.3300344391182 0.168173525797887 0
Spt No 202 2.4448821803182 0.168173525797887 0
Spt No 203 2.4448821803182 0.268173525797887 0
Spt No 204 2.3598437729182 0.268173525797887 0
Spt No 205 2.5296013445182 0.168173525797887 0
Spt No 206 2.5594904891182 0.268173525797887 0
Spt No 207 2.5296013445182 0.168173525797887 0
Spt No 208 2.6444091803182 0.168173525797887 0
Spt No 209 2.6444091803182 0.268173525797887 0
Spt No 210 2.5594904891182 0.268173525797887 0
Sln No  101 101 102
Sln No  102 102 103
Sln No  103 103 104
Sln No  104 104 101
Sln No  105 105 106
Sln No  106 106 107
Sln No  107 107 108
Sln No  108 108 105
Sln No  109 106 109
Sln No  110 109 110
Sln No  111 110 107
Sln No  112 107 106
Sln No  113 111 112
Sln No  114 112 113
Sln No  115 113 114
Sln No  116 114 111
Sln No  117 112 115
Sln No  118 115 116
Sln No  119 116 113
Sln No  120 113 112
Sln No  121 117 118
Sln No  122 118 119
Sln No  123 119 120
Sln No  124 120 117
Sln No  125 118 121
Sln No  126 121 122
Sln No  127 122 119
Sln No  128 119 118
Sln No  129 123 124
Sln No  130 124 125
Sln No  131 125 126
Sln No  132 126 123
Sln No  133 124 127
Sln No  134 127 128
Sln No  135 128 125
Sln No  136 125 124
Sln No  137 129 130
Sln No  138 130 131
Sln No  139 131 132
Sln No  140 132 129
Sln No  141 130 133
Sln No  142 133 134
Sln No  143 134 131
Sln No  144 131 130
Sln No  145 135 136
Sln No  146 136 137
Sln No  147 137 138
Sln No  148 138 135
Sln No  149 136 139
Sln No  150 139 140
Sln No  151 140 137
Sln No  152 137 136
Sln No  153 141 142
Sln No  154 142 143
Sln No  155 143 144
Sln No  156 144 141
Sln No  157 142 145
Sln No  158 145 146
Sln No  159 146 143
Sln No  160 143 142
Sln No  161 147 148
Sln No  162 148 149
Sln No  163 149 150
Sln No  164 150 147
Sln No  165 148 151
Sln No  166 151 152
Sln No  167 152 149
Sln No  168 149 148
Sln No  169 153 154
Sln No  170 154 155
Sln No  171 155 156
Sln No  172 156 153
Sln No  173 154 157
Sln No  174 157 158
Sln No  175 158 155
Sln No  176 155 154
Sln No  177 159 160
Sln No  178 160 161
Sln No  179 161 162
Sln No  180 162 159
Sln No  181 160 163
Sln No  182 163 164
Sln No  183 164 161
Sln No  184 161 160
Sln No  185 165 166
Sln No  186 166 167
Sln No  187 167 168
Sln No  188 168 165
Sln No  189 166 169
Sln No  190 169 170
Sln No  191 170 167
Sln No  192 167 166
Sln No  193 171 172
Sln No  194 172 173
Sln No  195 173 174
Sln No  196 174 171
Sln No  197 172 175
Sln No  198 175 176
Sln No  199 176 173
Sln No  200 173 172
Sln No  201 177 178
Sln No  202 178 179
Sln No  203 179 180
Sln No  204 180 177
Sln No  205 178 181
Sln No  206 181 182
Sln No  207 182 179
Sln No  208 179 178
Sln No  209 183 184
Sln No  210 184 185
Sln No  211 185 186
Sln No  212 186 183
Sln No  213 184 187
Sln No  214 187 188
Sln No  215 188 185
Sln No  216 185 184
Sln No  217 189 190
Sln No  218 190 191
Sln No  219 191 192
Sln No  220 192 189
Sln No  221 190 193
Sln No  222 193 194
Sln No  223 194 191
Sln No  224 191 190
Sln No  225 195 196
Sln No  226 196 197
Sln No  227 197 198
Sln No  228 198 195
Sln No  229 196 199
Sln No  230 199 200
Sln No  231 200 197
Sln No  232 197 196
Sln No  233 201 202
Sln No  234 202 203
Sln No  235 203 204
Sln No  236 204 201
Sln No  237 202 205
Sln No  238 205 206
Sln No  239 206 203
Sln No  240 203 202
Sln No  241 207 208
Sln No  242 208 209
Sln No  243 209 210
Sln No  244 210 207
Sar No 1001 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 101,102,103,104,
Sar No 1002 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 105,106,107,108,
Sar No 2003 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 109,110,111,112,
Sar No 2004 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 113,114,115,116,
Sar No 3005 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 117,118,119,120,
Sar No 3006 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 121,122,123,124,
Sar No 4007 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 125,126,127,128,
Sar No 4008 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 129,130,131,132,
Sar No 5009 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 133,134,135,136,
Sar No 5010 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 137,138,139,140,
Sar No 6011 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 141,142,143,144,
Sar No 6012 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 145,146,147,148,
Sar No 7013 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 149,150,151,152,
Sar No 7014 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 153,154,155,156,
Sar No 8015 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 157,158,159,160,
Sar No 8016 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 161,162,163,164,
Sar No 9017 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 165,166,167,168,
Sar No 9018 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 169,170,171,172,
Sar No 10019 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 173,174,175,176,
Sar No 10020 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 177,178,179,180,
Sar No 11021 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 181,182,183,184,
Sar No 11022 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 185,186,187,188,
Sar No 12023 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 189,190,191,192,
Sar No 12024 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 193,194,195,196,
Sar No 13025 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 197,198,199,200,
Sar No 13026 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 201,202,203,204,
Sar No 14027 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 205,206,207,208,
Sar No 14028 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 209,210,211,212,
Sar No 15029 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 213,214,215,216,
Sar No 15030 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 217,218,219,220,
Sar No 16031 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 221,222,223,224,
Sar No 16032 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 225,226,227,228,
Sar No 17033 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 229,230,231,232,
Sar No 17034 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 233,234,235,236,
Sar No 18035 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 237,238,239,240,
Sar No 18036 t 0.02 mno 1 grp 10 MCTL REGM
SarB Out NL 241,242,243,244,
Sln No -102
SlnS fix 'PXPYPZ' RefT >Sln Ref 108 CD 0.183847
Sln No -110
SlnS fix 'PXPYPZ' RefT >Sln Ref 116 CD 0.465065
Sln No -118
SlnS fix 'PXPYPZ' RefT >Sln Ref 124 CD 0.859174
Sln No -126
SlnS fix 'PXPYPZ' RefT >Sln Ref 132 CD 1.710966
Sln No -134
SlnS fix 'PXPYPZ' RefT >Sln Ref 140 CD 2.17164
Sln No -142
SlnS fix 'PXPYPZ' RefT >Sln Ref 148 CD 3.056492
Sln No -150
SlnS fix 'PXPYPZ' RefT >Sln Ref 156 CD 4.792568
Sln No -158
SlnS fix 'PXPYPZ' RefT >Sln Ref 164 CD 6.592499
Sln No -166
SlnS fix 'PXPYPZ' RefT >Sln Ref 172 CD 10.955144
Sln No -174
SlnS fix 'PXPYPZ' RefT >Sln Ref 180 CD 7.166551
Sln No -182
SlnS fix 'PXPYPZ' RefT >Sln Ref 188 CD 3.760597
Sln No -190
SlnS fix 'PXPYPZ' RefT >Sln Ref 196 CD 2.448822
Sln No -198
SlnS fix 'PXPYPZ' RefT >Sln Ref 204 CD 1.472478
Sln No -206
SlnS fix 'PXPYPZ' RefT >Sln Ref 212 CD 1.226579
Sln No -214
SlnS fix 'PXPYPZ' RefT >Sln Ref 220 CD 1.086598
Sln No -222
SlnS fix 'PXPYPZ' RefT >Sln Ref 228 CD 0.895056
Sln No -230
SlnS fix 'PXPYPZ' RefT >Sln Ref 236 CD 0.704393
Sln No -238
SlnS fix 'PXPYPZ' RefT >Sln Ref 244 CD 0.314892
Spt No -101 Fix PXPYPZ
Spt No -104 Fix PXPYPZ
Spt No -208 Fix PXPYPZ
Spt No -209 Fix PXPYPZ
END


+PROG TEMPLATE urs:2
HEAD Variables definition
 #define ini = 1         $load/displacement first increment
 #define max = 1         $final value of load/displacement
 sto#l1 2              $Reference loadcase
 sto#lc 10                $Start loadcase
 sto#x 1000              $maximum iteration loops
 sto#I 50                $optimal number of inner iterations
END
+PROG SOFILOAD urs:4
HEAD Initial lifting
PAGE UNII 0

LC 1 TYPE NONE

AREA SAR 1001 TYPE PZ -20
AREA SAR 1002 TYPE PZ -20
AREA SAR 2003 TYPE PZ -20
AREA SAR 2004 TYPE PZ -20
AREA SAR 3005 TYPE PZ -20
AREA SAR 3006 TYPE PZ -20
AREA SAR 4007 TYPE PZ -20
AREA SAR 4008 TYPE PZ -20
AREA SAR 5009 TYPE PZ -20
AREA SAR 5010 TYPE PZ -20
AREA SAR 6011 TYPE PZ -20
AREA SAR 6012 TYPE PZ -20
AREA SAR 7013 TYPE PZ -20
AREA SAR 7014 TYPE PZ -20
AREA SAR 8015 TYPE PZ -20
AREA SAR 8016 TYPE PZ -20
AREA SAR 9017 TYPE PZ -20
AREA SAR 9018 TYPE PZ -20
AREA SAR 10019 TYPE PZ -20
AREA SAR 10020 TYPE PZ -20
AREA SAR 11021 TYPE PZ -20
AREA SAR 11022 TYPE PZ -20
AREA SAR 12023 TYPE PZ -20
AREA SAR 12024 TYPE PZ -20
AREA SAR 13025 TYPE PZ -20
AREA SAR 13026 TYPE PZ -20
AREA SAR 14027 TYPE PZ -20
AREA SAR 14028 TYPE PZ -20
AREA SAR 15029 TYPE PZ -20
AREA SAR 15030 TYPE PZ -20
AREA SAR 16031 TYPE PZ -20
AREA SAR 16032 TYPE PZ -20
AREA SAR 17033 TYPE PZ -20
AREA SAR 17034 TYPE PZ -20
AREA SAR 18035 TYPE PZ -20
AREA SAR 18036 TYPE PZ -20

END

+PROG ASE urs:6
HEAD Initial lifting
PAGE UNII 0
SYST prob TH3 tol 0.01

LC 1


END

+PROG SOFILOAD urs:7
HEAD
PAGE UNII 0

LC 2 TYPE NONE

AREA SAR 1001 TYPE PZ -10
AREA SAR 1002 TYPE PZ -10
AREA SAR 2003 TYPE PZ -10
AREA SAR 2004 TYPE PZ -10
AREA SAR 3005 TYPE PZ -10
AREA SAR 3006 TYPE PZ -10
AREA SAR 4007 TYPE PZ -10
AREA SAR 4008 TYPE PZ -10
AREA SAR 5009 TYPE PZ -10
AREA SAR 5010 TYPE PZ -10
AREA SAR 6011 TYPE PZ -10
AREA SAR 6012 TYPE PZ -10
AREA SAR 7013 TYPE PZ -10
AREA SAR 7014 TYPE PZ -10
AREA SAR 8015 TYPE PZ -10
AREA SAR 8016 TYPE PZ -10
AREA SAR 9017 TYPE PZ -10
AREA SAR 9018 TYPE PZ -10
AREA SAR 10019 TYPE PZ -10
AREA SAR 10020 TYPE PZ -10
AREA SAR 11021 TYPE PZ -10
AREA SAR 11022 TYPE PZ -10
AREA SAR 12023 TYPE PZ -10
AREA SAR 12024 TYPE PZ -10
AREA SAR 13025 TYPE PZ -10
AREA SAR 13026 TYPE PZ -10
AREA SAR 14027 TYPE PZ -10
AREA SAR 14028 TYPE PZ -10
AREA SAR 15029 TYPE PZ -10
AREA SAR 15030 TYPE PZ -10
AREA SAR 16031 TYPE PZ -10
AREA SAR 16032 TYPE PZ -10
AREA SAR 17033 TYPE PZ -10
AREA SAR 17034 TYPE PZ -10
AREA SAR 18035 TYPE PZ -10
AREA SAR 18036 TYPE PZ -10
NODE no 104 TYPE wxx p1 4.568
NODE no 104 TYPE wyy p1 -2.144
NODE no 104 TYPE wzz p1 -13.731
NODE no 101 TYPE wxx p1 -4.568
NODE no 101 TYPE wyy p1 2.144
NODE no 101 TYPE wzz p1 13.731
NODE no 208 TYPE wxx p1 -1841.585
NODE no 208 TYPE wyy p1 1032.407
NODE no 208 TYPE wzz p1 -11.907
NODE no 209 TYPE wxx p1 -1825.774
NODE no 209 TYPE wyy p1 1028.233
NODE no 209 TYPE wzz p1 11.907


END

+PROG ASE urs:5
HEAD Solve
 ECHO val off

 let#n 0.001       $initial load factor
 let#iter 0        $iteration counter
 let#plc  1        $initial primary load case
 let#s #n

 let#lambda0 0.0        $initial factor for adaptive increments
 let#lambda  0.001      $initial load factor
 let#lambdaInc #lambda  $sum of the increments
 let#lambdaIncDiff #lambda  $sum of the increments

 let#tol 0.01

<TEXT,FILE=loadIncrement.txt>
'Iteration / Convergence / Lambda / Lambda increment'

</TEXT>

 loop #x
  let#break 1
   loop
    SYST  prob th3  iter 100  plc #plc  fmax 1.0  tol #tol  $tol4 5*#tol  tol8 10*#tol
    CTRL opt iter  val 3  v2 1          $ update tangential stiffness by each iteration
    $CTRL iter v4 10

    CTRL opt cabl 0
    CTRL opt beam 5

    GRP - yes

    $let#lambda ($(ini)*#n)
    LC  no #lc  dlz 0.0  titl 'iter #iter, lambda #(#lambda,.3)'
     $if #lambda<0.5
      LCC  no #l1  fact #lambda
     $else
      $LCC  no #l1  fact 1.001-#lambda
     $endif
     $LCC  no #l2  fact #lambdaInc
     $LCC  no #l3  fact #lambdaIncDiff
    end

    sto#finalLc #lc

    @key LC_CTRL #lc
    let#converge @cri1
     if (#converge > 0)
      let#break 0

<TEXT,FILE=+loadIncrement.txt>
#lc       'yes'           #(#lambda,.3)    #(#lambdaInc,.3)
</TEXT>
     else
      let#n (#n-0.5*#s)
      let#s 0.5*#s
      let#lambdaInc #s/(1.0-#n0)

<TEXT,FILE=+loadIncrement.txt>
#lc       'no'            #(#lambda,.3)    #(#lambdaInc,.3)
</TEXT>
     endif

    endloop #break

    $__________options of increasing increments__________
    if #lambda < 1.0
     $ let#n   #n+0.01                $ linear increment
     $let#n   #iter*#iter            $ quadratic increment
     $let#n   #n*1.1                 $ nonl increment

     @key LC_CTRL #lc                $ adaptive increments
     let#Ik @cri1
     let#s0 #lambda-#lambda0
     let#s (#s0*(#I/#Ik)^0.1)
     if (#lambda+#s) > 1.0
      let#s 1.0-#lambda
     endif
     $ lambda -> for external loads
     let#lambda0 #lambda
     let#lambda #lambda+#s
     $ Incremental lambda -> for prestress and temperature
     let#lambdaInc #lambda-#lambda0
     $ Incremental differential lambda -> for strain
     let#lambdaIncDiff (#lambda-#lambda0)/(1.0-#lambda0)

    $____________________________________________________

     let#plc #lc     $ use previous lc as primary

    elseif #lambda == 1.0
     end
    endif

    let#iter #iter+1
    let#lc #lc+1
    let#tol #tol*#iter^2

  endloop

END


+PROG ASE urs:8
HEAD Dynamic relaxation
let#lf 201 ; let#plf #finalLc
loop#i0 20
  let#i1 #i0+1         $ 1,2,3, 4, 5... 20
  let#dt 0.01*#i1      $ 0.01,0.02, ...  0.25 sec.
  STEP 1 DT #dt DIV 0 BET 0.4 0.55 1.0
  SYST PLC #plf PROB TH3 FMAX 0.5 TOL -0.1 TOL8 -1
  ECHO GRP FULL
  ECHO STAT FULL
  GRP - RADB 0.001 $ nicht RADA verwenden, sonst zu hohe Dämpfung(Festhaltung) bei t>10 sec !
      $ RADB 0.001 : Dämpfung ca. 100 bei 0.01sec, Dämpfung ca. 2 bei 1.0sec,
  LC #lf FACD 0.0 TITL 'Relaxation with dynamic damping'
  END
  let#plf #lf ; let#lf #lf+1
endloop
END

+prog ase urs:9
head Update coordiante system to form found geometry
syst PLC 220 STOR yes
ctrl diff
end