#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 25.539392[mpa]
MATE 3 E 64.376245[mpa]
MATE 4 E 121.124128[mpa]
MATE 5 E 225.471054[mpa]
MATE 6 E 285.954754[mpa]
MATE 7 E 413.360933[mpa]
MATE 8 E 597.838091[mpa]
MATE 9 E 837.870407[mpa]
MATE 10 E 1147.009484[mpa]
MATE 11 E 537.198688[mpa]
MATE 12 E 351.770186[mpa]
MATE 13 E 203.919001[mpa]
MATE 14 E 151.181845[mpa]
MATE 15 E 137.303463[mpa]
MATE 16 E 112.171076[mpa]
MATE 17 E 88.680602[mpa]
MATE 18 E 40.221783[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.851287833011216 0.168173525797887 0
Spt No 103 -0.838901326114157 0.268173525797887 0
Spt No 104 -0.947076819681804 0.268173525797887 0
Spt No 105 -0.851287833011216 0.168173525797887 0
Spt No 106 -0.743183702034745 0.168173525797887 0
Spt No 107 -0.743183702034745 0.268173525797887 0
Spt No 108 -0.838901326114157 0.268173525797887 0
Spt No 109 -0.647476272611216 0.168173525797887 0
Spt No 110 -0.634998013811216 0.268173525797887 0
Spt No 111 -0.647476272611216 0.168173525797887 0
Spt No 112 -0.539290584387686 0.168173525797887 0
Spt No 113 -0.539290584387686 0.268173525797887 0
Spt No 114 -0.634998013811216 0.268173525797887 0
Spt No 115 -0.443583154964157 0.168173525797887 0
Spt No 116 -0.431104896164157 0.268173525797887 0
Spt No 117 -0.443583154964157 0.168173525797887 0
Spt No 118 -0.335397466740627 0.168173525797887 0
Spt No 119 -0.335397466740627 0.268173525797887 0
Spt No 120 -0.431104896164157 0.268173525797887 0
Spt No 121 -0.239690037317098 0.168173525797887 0
Spt No 122 -0.227211778517098 0.268173525797887 0
Spt No 123 -0.239690037317098 0.168173525797887 0
Spt No 124 -0.131504349093569 0.168173525797887 0
Spt No 125 -0.131504349093569 0.268173525797887 0
Spt No 126 -0.227211778517098 0.268173525797887 0
Spt No 127 -0.0357969196700395 0.168173525797887 0
Spt No 128 -0.0233186608700394 0.268173525797887 0
Spt No 129 -0.0357969196700395 0.168173525797887 0
Spt No 130 0.0723887685534901 0.168173525797887 0
Spt No 131 0.0723887685534901 0.268173525797887 0
Spt No 132 -0.0233186608700394 0.268173525797887 0
Spt No 133 0.16809619797702 0.168173525797887 0
Spt No 134 0.180574456777019 0.268173525797887 0
Spt No 135 0.16809619797702 0.168173525797887 0
Spt No 136 0.276281886200549 0.168173525797887 0
Spt No 137 0.276281886200549 0.268173525797887 0
Spt No 138 0.180574456777019 0.268173525797887 0
Spt No 139 0.371989315624078 0.168173525797887 0
Spt No 140 0.384467574424078 0.268173525797887 0
Spt No 141 0.371989315624078 0.168173525797887 0
Spt No 142 0.480175003847607 0.168173525797887 0
Spt No 143 0.480175003847607 0.268173525797887 0
Spt No 144 0.384467574424078 0.268173525797887 0
Spt No 145 0.575882433271137 0.168173525797887 0
Spt No 146 0.588360692071137 0.268173525797887 0
Spt No 147 0.575882433271137 0.168173525797887 0
Spt No 148 0.684068121494667 0.168173525797887 0
Spt No 149 0.684068121494667 0.268173525797887 0
Spt No 150 0.588360692071137 0.268173525797887 0
Spt No 151 0.779775550918196 0.168173525797887 0
Spt No 152 0.792253809718196 0.268173525797887 0
Spt No 153 0.779775550918196 0.168173525797887 0
Spt No 154 0.887961239141725 0.168173525797887 0
Spt No 155 0.887961239141725 0.268173525797887 0
Spt No 156 0.792253809718196 0.268173525797887 0
Spt No 157 0.983668668565255 0.168173525797887 0
Spt No 158 0.996146927365255 0.268173525797887 0
Spt No 159 0.983668668565255 0.168173525797887 0
Spt No 160 1.09185435678878 0.168173525797887 0
Spt No 161 1.09185435678878 0.268173525797887 0
Spt No 162 0.996146927365255 0.268173525797887 0
Spt No 163 1.18756178621231 0.168173525797887 0
Spt No 164 1.20004004501231 0.268173525797887 0
Spt No 165 1.18756178621231 0.168173525797887 0
Spt No 166 1.29574747443584 0.168173525797887 0
Spt No 167 1.29574747443584 0.268173525797887 0
Spt No 168 1.20004004501231 0.268173525797887 0
Spt No 169 1.39145490385937 0.168173525797887 0
Spt No 170 1.40393316265937 0.268173525797887 0
Spt No 171 1.39145490385937 0.168173525797887 0
Spt No 172 1.4996405920829 0.168173525797887 0
Spt No 173 1.4996405920829 0.268173525797887 0
Spt No 174 1.40393316265937 0.268173525797887 0
Spt No 175 1.59534802150643 0.168173525797887 0
Spt No 176 1.60782628030643 0.268173525797887 0
Spt No 177 1.59534802150643 0.168173525797887 0
Spt No 178 1.70353370972996 0.168173525797887 0
Spt No 179 1.70353370972996 0.268173525797887 0
Spt No 180 1.60782628030643 0.268173525797887 0
Spt No 181 1.79924113915349 0.168173525797887 0
Spt No 182 1.81171939795349 0.268173525797887 0
Spt No 183 1.79924113915349 0.168173525797887 0
Spt No 184 1.90742682737702 0.168173525797887 0
Spt No 185 1.90742682737702 0.268173525797887 0
Spt No 186 1.81171939795349 0.268173525797887 0
Spt No 187 2.00313425680055 0.168173525797887 0
Spt No 188 2.01561251560055 0.268173525797887 0
Spt No 189 2.00313425680055 0.168173525797887 0
Spt No 190 2.11131994502408 0.168173525797887 0
Spt No 191 2.11131994502408 0.268173525797887 0
Spt No 192 2.01561251560055 0.268173525797887 0
Spt No 193 2.20702737444761 0.168173525797887 0
Spt No 194 2.21950563324761 0.268173525797887 0
Spt No 195 2.20702737444761 0.168173525797887 0
Spt No 196 2.31521306267114 0.168173525797887 0
Spt No 197 2.31521306267114 0.268173525797887 0
Spt No 198 2.21950563324761 0.268173525797887 0
Spt No 199 2.41093068675055 0.168173525797887 0
Spt No 200 2.42331719364761 0.268173525797887 0
Spt No 201 2.41093068675055 0.168173525797887 0
Spt No 202 2.5191061803182 0.168173525797887 0
Spt No 203 2.5191061803182 0.268173525797887 0
Spt No 204 2.42331719364761 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
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,
Sln No -102
SlnS fix 'PXPYPZ' RefT >Sln Ref 108 CD 0.08372
Sln No -110
SlnS fix 'PXPYPZ' RefT >Sln Ref 116 CD 0.211854
Sln No -118
SlnS fix 'PXPYPZ' RefT >Sln Ref 124 CD 0.400894
Sln No -126
SlnS fix 'PXPYPZ' RefT >Sln Ref 132 CD 0.754225
Sln No -134
SlnS fix 'PXPYPZ' RefT >Sln Ref 140 CD 0.962506
Sln No -142
SlnS fix 'PXPYPZ' RefT >Sln Ref 148 CD 1.409838
Sln No -150
SlnS fix 'PXPYPZ' RefT >Sln Ref 156 CD 2.079036
Sln No -158
SlnS fix 'PXPYPZ' RefT >Sln Ref 164 CD 2.990106
Sln No -166
SlnS fix 'PXPYPZ' RefT >Sln Ref 172 CD 4.236266
Sln No -174
SlnS fix 'PXPYPZ' RefT >Sln Ref 180 CD 1.856186
Sln No -182
SlnS fix 'PXPYPZ' RefT >Sln Ref 188 CD 1.192113
Sln No -190
SlnS fix 'PXPYPZ' RefT >Sln Ref 196 CD 0.680631
Sln No -198
SlnS fix 'PXPYPZ' RefT >Sln Ref 204 CD 0.501905
Sln No -206
SlnS fix 'PXPYPZ' RefT >Sln Ref 212 CD 0.455189
Sln No -214
SlnS fix 'PXPYPZ' RefT >Sln Ref 220 CD 0.370925
Sln No -222
SlnS fix 'PXPYPZ' RefT >Sln Ref 228 CD 0.292552
Sln No -230
SlnS fix 'PXPYPZ' RefT >Sln Ref 236 CD 0.132044
Spt No -101 Fix PXPYPZ
Spt No -104 Fix PXPYPZ
Spt No -202 Fix PXPYPZ
Spt No -203 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 PZZ 20
AREA SAR 1002 TYPE PZZ 20
AREA SAR 2003 TYPE PZZ 20
AREA SAR 2004 TYPE PZZ 20
AREA SAR 3005 TYPE PZZ 20
AREA SAR 3006 TYPE PZZ 20
AREA SAR 4007 TYPE PZZ 20
AREA SAR 4008 TYPE PZZ 20
AREA SAR 5009 TYPE PZZ 20
AREA SAR 5010 TYPE PZZ 20
AREA SAR 6011 TYPE PZZ 20
AREA SAR 6012 TYPE PZZ 20
AREA SAR 7013 TYPE PZZ 20
AREA SAR 7014 TYPE PZZ 20
AREA SAR 8015 TYPE PZZ 20
AREA SAR 8016 TYPE PZZ 20
AREA SAR 9017 TYPE PZZ 20
AREA SAR 9018 TYPE PZZ 20
AREA SAR 10019 TYPE PZZ 20
AREA SAR 10020 TYPE PZZ 20
AREA SAR 11021 TYPE PZZ 20
AREA SAR 11022 TYPE PZZ 20
AREA SAR 12023 TYPE PZZ 20
AREA SAR 12024 TYPE PZZ 20
AREA SAR 13025 TYPE PZZ 20
AREA SAR 13026 TYPE PZZ 20
AREA SAR 14027 TYPE PZZ 20
AREA SAR 14028 TYPE PZZ 20
AREA SAR 15029 TYPE PZZ 20
AREA SAR 15030 TYPE PZZ 20
AREA SAR 16031 TYPE PZZ 20
AREA SAR 16032 TYPE PZZ 20
AREA SAR 17033 TYPE PZZ 20
AREA SAR 17034 TYPE PZZ 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 PZZ 10
AREA SAR 1002 TYPE PZZ 10
AREA SAR 2003 TYPE PZZ 10
AREA SAR 2004 TYPE PZZ 10
AREA SAR 3005 TYPE PZZ 10
AREA SAR 3006 TYPE PZZ 10
AREA SAR 4007 TYPE PZZ 10
AREA SAR 4008 TYPE PZZ 10
AREA SAR 5009 TYPE PZZ 10
AREA SAR 5010 TYPE PZZ 10
AREA SAR 6011 TYPE PZZ 10
AREA SAR 6012 TYPE PZZ 10
AREA SAR 7013 TYPE PZZ 10
AREA SAR 7014 TYPE PZZ 10
AREA SAR 8015 TYPE PZZ 10
AREA SAR 8016 TYPE PZZ 10
AREA SAR 9017 TYPE PZZ 10
AREA SAR 9018 TYPE PZZ 10
AREA SAR 10019 TYPE PZZ 10
AREA SAR 10020 TYPE PZZ 10
AREA SAR 11021 TYPE PZZ 10
AREA SAR 11022 TYPE PZZ 10
AREA SAR 12023 TYPE PZZ 10
AREA SAR 12024 TYPE PZZ 10
AREA SAR 13025 TYPE PZZ 10
AREA SAR 13026 TYPE PZZ 10
AREA SAR 14027 TYPE PZZ 10
AREA SAR 14028 TYPE PZZ 10
AREA SAR 15029 TYPE PZZ 10
AREA SAR 15030 TYPE PZZ 10
AREA SAR 16031 TYPE PZZ 10
AREA SAR 16032 TYPE PZZ 10
AREA SAR 17033 TYPE PZZ 10
AREA SAR 17034 TYPE PZZ 10
NODE no 104 TYPE wxx p1 1.935
NODE no 104 TYPE wyy p1 -0.382
NODE no 104 TYPE wzz p1 -5.817
NODE no 101 TYPE wxx p1 -1.935
NODE no 101 TYPE wyy p1 0.382
NODE no 101 TYPE wzz p1 5.817
NODE no 202 TYPE wxx p1 -1711.722
NODE no 202 TYPE wyy p1 421.348
NODE no 202 TYPE wzz p1 -5.04
NODE no 203 TYPE wxx p1 -1705.03
NODE no 203 TYPE wyy p1 420.612
NODE no 203 TYPE wzz p1 5.04
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
'Iteration / Convergence / Lambda / Lambda increment'
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
#lc 'yes' #(#lambda,.3) #(#lambdaInc,.3)
else
let#n (#n-0.5*#s)
let#s 0.5*#s
let#lambdaInc #s/(1.0-#n0)
#lc 'no' #(#lambda,.3) #(#lambdaInc,.3)
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
sto#fflc #lc-1
END
+PROG ASE urs:8
HEAD Dynamic relaxation
let#lf 201 ; let#plf 29
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.0001 TITL 'Relaxation with dynamic damping'
END
let#plf #lf ; let#lf #lf+1
endloop
END