How to tackle the convergence problems in OpenSees?

[Mohsenkalatehahani]

I modeled a four-story frame using the Ibarra-Krawinkler deterioration model (including concrete slab contribution / panel zone) and I have faced convergence problems in dynamic analysis! I adopted the adaptive strategy which is introduced in OpenSees Message Board, However, it is not successful. The tcl file is here. I would be grateful if anybody could help me.

wipe all;
model BasicBuilder -ndm 2 -ndf 3;
set dataDir DataDynamic;
file mkdir $dataDir;
set GMdir "D:/ResidualCapacity/ModelA/YekMordad/";
source rotSpring2DModIKModel.tcl;
source rotLeaningCol.tcl;
source rotPanelZone2D.tcl;
source elemPanelZone2D.tcl;
source DisplayModel2D.tcl;
source DisplayPlane.tcl;
node 11 0 0
node 12 9.14 0
node 13 18.28 0
node 14 27.42 0
node 24 27.42 4.57
node 34 27.42 8.23
node 44 27.42 11.89
node 54 27.42 15.55
node 117 0 0
node 127 9.14 0
node 137 18.28 0
node 215 0 4.226
node 216 0 4.226
node 217 0 4.914
node 218 0 4.914
node 211 0.70127 4.57
node 212 0.70127 4.57
node 225 9.14 4.226
node 226 9.14 4.226
node 227 9.14 4.914
node 228 9.14 4.914
node 221 9.84127 4.57
node 222 9.84127 4.57
node 223 8.43873 4.57
node 224 8.43873 4.57
node 235 18.28 4.226
node 236 18.28 4.226
node 237 18.28 4.914
node 238 18.28 4.914
node 233 17.5787 4.57
node 234 17.5787 4.57
node 315 0 7.886
node 316 0 7.886
node 317 0 8.574
node 318 0 8.574
node 311 0.70127 8.23
node 312 0.70127 8.23
node 325 9.14 7.886
node 326 9.14 7.886
node 327 9.14 8.574
node 328 9.14 8.574
node 321 9.84127 8.23
node 322 9.84127 8.23
node 323 8.43873 8.23
node 324 8.43873 8.23
node 335 18.28 7.886
node 336 18.28 7.886
node 337 18.28 8.574
node 338 18.28 8.574
node 333 17.5787 8.23
node 334 17.5787 8.23
node 415 0 11.546
node 416 0 11.546
node 417 0 12.234
node 418 0 12.234
node 411 0.70127 11.89
node 412 0.70127 11.89
node 425 9.14 11.546
node 426 9.14 11.546
node 427 9.14 12.234
node 428 9.14 12.234
node 421 9.84127 11.89
node 422 9.84127 11.89
node 423 8.43873 11.89
node 424 8.43873 11.89
node 435 18.28 11.546
node 436 18.28 11.546
node 437 18.28 12.234
node 438 18.28 12.234
node 433 17.5787 11.89
node 434 17.5787 11.89
node 515 0 15.206
node 516 0 15.206
node 517 0 15.894
node 511 0.70127 15.55
node 512 0.70127 15.55
node 525 9.14 15.206
node 526 9.14 15.206
node 527 9.14 15.894
node 521 9.84127 15.55
node 522 9.84127 15.55
node 523 8.43873 15.55
node 524 8.43873 15.55
node 535 18.28 15.206
node 536 18.28 15.206
node 537 18.28 15.894
node 533 17.5787 15.55
node 534 17.5787 15.55
node 246 27.42 4.57
node 247 27.42 4.57
node 346 27.42 8.23
node 347 27.42 8.23
node 446 27.42 11.89
node 447 27.42 11.89
node 546 27.42 15.55
node 2101 -0.311 4.914
node 2102 -0.311 4.914
node 2103 0.311 4.914
node 2104 0.311 4.914
node 2105 0.311 4.57
node 2106 0.311 4.226
node 2107 0.311 4.226
node 2108 -0.311 4.226
node 2109 -0.311 4.226
node 2110 -0.311 4.57
node 2201 8.829 4.914
node 2202 8.829 4.914
node 2203 9.451 4.914
node 2204 9.451 4.914
node 2205 9.451 4.57
node 2206 9.451 4.226
node 2207 9.451 4.226
node 2208 8.829 4.226
node 2209 8.829 4.226
node 2210 8.829 4.57
node 2301 17.969 4.914
node 2302 17.969 4.914
node 2303 18.591 4.914
node 2304 18.591 4.914
node 2305 18.591 4.57
node 2306 18.591 4.226
node 2307 18.591 4.226
node 2308 17.969 4.226
node 2309 17.969 4.226
node 2310 17.969 4.57
node 3101 -0.311 8.574
node 3102 -0.311 8.574
node 3103 0.311 8.574
node 3104 0.311 8.574
node 3105 0.311 8.23
node 3106 0.311 7.886
node 3107 0.311 7.886
node 3108 -0.311 7.886
node 3109 -0.311 7.886
node 3110 -0.311 8.23
node 3201 8.829 8.574
node 3202 8.829 8.574
node 3203 9.451 8.574
node 3204 9.451 8.574
node 3205 9.451 8.23
node 3206 9.451 7.886
node 3207 9.451 7.886
node 3208 8.829 7.886
node 3209 8.829 7.886
node 3210 8.829 8.23
node 3301 17.969 8.574
node 3302 17.969 8.574
node 3303 18.591 8.574
node 3304 18.591 8.574
node 3305 18.591 8.23
node 3306 18.591 7.886
node 3307 18.591 7.886
node 3308 17.969 7.886
node 3309 17.969 7.886
node 3310 17.969 8.23
node 4101 -0.311 12.234
node 4102 -0.311 12.234
node 4103 0.311 12.234
node 4104 0.311 12.234
node 4105 0.311 11.89
node 4106 0.311 11.546
node 4107 0.311 11.546
node 4108 -0.311 11.546
node 4109 -0.311 11.546
node 4110 -0.311 11.89
node 4201 8.829 12.234
node 4202 8.829 12.234
node 4203 9.451 12.234
node 4204 9.451 12.234
node 4205 9.451 11.89
node 4206 9.451 11.546
node 4207 9.451 11.546
node 4208 8.829 11.546
node 4209 8.829 11.546
node 4210 8.829 11.89
node 4301 17.969 12.234
node 4302 17.969 12.234
node 4303 18.591 12.234
node 4304 18.591 12.234
node 4305 18.591 11.89
node 4306 18.591 11.546
node 4307 18.591 11.546
node 4308 17.969 11.546
node 4309 17.969 11.546
node 4310 17.969 11.89
node 5101 -0.311 15.894
node 5102 -0.311 15.894
node 5103 0.311 15.894
node 5104 0.311 15.894
node 5105 0.311 15.55
node 5106 0.311 15.206
node 5107 0.311 15.206
node 5108 -0.311 15.206
node 5109 -0.311 15.206
node 5110 -0.311 15.55
node 5201 8.829 15.894
node 5202 8.829 15.894
node 5203 9.451 15.894
node 5204 9.451 15.894
node 5205 9.451 15.55
node 5206 9.451 15.206
node 5207 9.451 15.206
node 5208 8.829 15.206
node 5209 8.829 15.206
node 5210 8.829 15.55
node 5301 17.969 15.894
node 5302 17.969 15.894
node 5303 18.591 15.894
node 5304 18.591 15.894
node 5305 18.591 15.55
node 5306 18.591 15.206
node 5307 18.591 15.206
node 5308 17.969 15.206
node 5309 17.969 15.206
node 5310 17.969 15.55
equalDOF 2105 2205 1;
equalDOF 2105 2305 1;
equalDOF 3105 3205 1;
equalDOF 3105 3305 1;
equalDOF 4105 4205 1;
equalDOF 4105 4305 1;
equalDOF 5105 5205 1;
equalDOF 5105 5305 1;
fix 11 1 1 1;
fix 12 1 1 1;
fix 13 1 1 1;
fix 14 1 1 0;
set PDeltaTransf 1;
geomTransf PDelta $PDeltaTransf ;
set Es 2.039432e+10;
set n 10;
element elasticBeamColumn 20211 2105 211 0.0194 $Es 0.00151 $PDeltaTransf;
element elasticBeamColumn 2021 212 223 0.0194 $Es 0.001661 $PDeltaTransf;
element elasticBeamColumn 20212 224 2210 0.0194 $Es 0.00151 $PDeltaTransf;
element elasticBeamColumn 20221 2205 221 0.0194 $Es 0.00151 $PDeltaTransf;
element elasticBeamColumn 2022 222 233 0.0194 $Es 0.001661 $PDeltaTransf;
element elasticBeamColumn 20222 234 2310 0.0194 $Es 0.00151 $PDeltaTransf;
element elasticBeamColumn 20311 3105 311 0.0194 $Es 0.00151 $PDeltaTransf;
element elasticBeamColumn 2031 312 323 0.0194 $Es 0.001661 $PDeltaTransf;
element elasticBeamColumn 20312 324 3210 0.0194 $Es 0.00151 $PDeltaTransf;
element elasticBeamColumn 20321 3205 321 0.0194 $Es 0.00151 $PDeltaTransf;
element elasticBeamColumn 2032 322 333 0.0194 $Es 0.001661 $PDeltaTransf;
element elasticBeamColumn 20322 334 3310 0.0194 $Es 0.00151 $PDeltaTransf;
element elasticBeamColumn 20411 4105 411 0.0176 $Es 0.000862 $PDeltaTransf;
element elasticBeamColumn 2041 412 423 0.0176 $Es 0.0009482 $PDeltaTransf;
element elasticBeamColumn 20412 424 4210 0.0176 $Es 0.000862 $PDeltaTransf;
element elasticBeamColumn 20421 4205 421 0.0176 $Es 0.000862 $PDeltaTransf;
element elasticBeamColumn 2042 422 433 0.0176 $Es 0.0009482 $PDeltaTransf;
element elasticBeamColumn 20422 434 4310 0.0176 $Es 0.000862 $PDeltaTransf;
element elasticBeamColumn 20511 5105 511 0.0176 $Es 0.000862 $PDeltaTransf;
element elasticBeamColumn 2051 512 523 0.0176 $Es 0.0009482 $PDeltaTransf;
element elasticBeamColumn 20512 524 5210 0.0176 $Es 0.000862 $PDeltaTransf;
element elasticBeamColumn 20521 5205 521 0.0176 $Es 0.000862 $PDeltaTransf;
element elasticBeamColumn 2052 522 533 0.0176 $Es 0.0009482 $PDeltaTransf;
element elasticBeamColumn 20522 534 5310 0.0176 $Es 0.000862 $PDeltaTransf;
element elasticBeamColumn 1011 117 215 0.0249 $Es 0.001837 $PDeltaTransf;
element elasticBeamColumn 1012 127 225 0.0249 $Es 0.001837 $PDeltaTransf;
element elasticBeamColumn 1013 137 235 0.0249 $Es 0.001837 $PDeltaTransf;
element elasticBeamColumn 1021 218 315 0.0249 $Es 0.001837 $PDeltaTransf;
element elasticBeamColumn 1022 228 325 0.0249 $Es 0.001837 $PDeltaTransf;
element elasticBeamColumn 1023 238 335 0.0249 $Es 0.001837 $PDeltaTransf;
element elasticBeamColumn 1031 318 415 0.0145 $Es 0.0009614 $PDeltaTransf;
element elasticBeamColumn 1032 328 425 0.0145 $Es 0.0009614 $PDeltaTransf;
element elasticBeamColumn 1033 338 435 0.0145 $Es 0.0009614 $PDeltaTransf;
element elasticBeamColumn 1041 418 515 0.0145 $Es 0.0009614 $PDeltaTransf;
element elasticBeamColumn 1042 428 525 0.0145 $Es 0.0009614 $PDeltaTransf;
element elasticBeamColumn 1043 438 535 0.0145 $Es 0.0009614 $PDeltaTransf;
set TrussMatID 600;
set Arigid 5;
set Irigid 1;
uniaxialMaterial Elastic $TrussMatID $Es;
element truss 3023 2305 24 $Arigid $TrussMatID;
element truss 3033 3305 34 $Arigid $TrussMatID;
element truss 3043 4305 44 $Arigid $TrussMatID;
element truss 3053 5305 54 $Arigid $TrussMatID;
element elasticBeamColumn 4014 14 246 $Arigid $Es $Irigid $PDeltaTransf;
element elasticBeamColumn 4024 247 346 $Arigid $Es $Irigid $PDeltaTransf;
element elasticBeamColumn 4034 347 446 $Arigid $Es $Irigid $PDeltaTransf;
element elasticBeamColumn 4044 447 546 $Arigid $Es $Irigid $PDeltaTransf;
set Apz 5;
set Ipz 1;
elemPanelZone2D 500211 2101 $Es $Apz $Ipz $PDeltaTransf;
elemPanelZone2D 500221 2201 $Es $Apz $Ipz $PDeltaTransf;
elemPanelZone2D 500231 2301 $Es $Apz $Ipz $PDeltaTransf;
elemPanelZone2D 500311 3101 $Es $Apz $Ipz $PDeltaTransf;
elemPanelZone2D 500321 3201 $Es $Apz $Ipz $PDeltaTransf;
elemPanelZone2D 500331 3301 $Es $Apz $Ipz $PDeltaTransf;
elemPanelZone2D 500411 4101 $Es $Apz $Ipz $PDeltaTransf;
elemPanelZone2D 500421 4201 $Es $Apz $Ipz $PDeltaTransf;
elemPanelZone2D 500431 4301 $Es $Apz $Ipz $PDeltaTransf;
elemPanelZone2D 500511 5101 $Es $Apz $Ipz $PDeltaTransf;
elemPanelZone2D 500521 5201 $Es $Apz $Ipz $PDeltaTransf;
elemPanelZone2D 500531 5301 $Es $Apz $Ipz $PDeltaTransf;
mass 2105 3798.75 1e-9 1e-9
mass 2205 7525.93 1e-9 1e-9
mass 2305 3798.75 1e-9 1e-9
mass 3105 3773.85 1e-9 1e-9
mass 3205 7501.03 1e-9 1e-9
mass 3305 3773.85 1e-9 1e-9
mass 4105 3755.9 1e-9 1e-9
mass 4205 7477.56 1e-9 1e-9
mass 4305 3755.9 1e-9 1e-9
mass 5105 3738.78 1e-9 1e-9
mass 5205 7460.43 1e-9 1e-9
mass 5305 3738.78 1e-9 1e-9
set Ks_col 5.31913e+08;
set MycolPos 257766;
set MycolNeg 257766;
set McMyPos 1.1;
set McMyNeg 1.1;
set LS 1000.0;
set LK 1000.0;
set LA 1000.0;
set LD 1000.0;
set cS 1.0;
set cK 1.0;
set cA 1.0;
set cD 1.0;
set th_pP 0.232571;
set th_pN 0.232571;
set th_pcP 1.33699;
set th_pcN 1.33699;
set ResP 0.4;
set ResN 0.4;
set th_uP 0.55;
set th_uN 0.55;
set DP 1.0;
set DN 1.0;
set a_memPos [expr ($n+1.0)*($MycolPos*($McMyPos-1.0)) / ($Ks_col*$th_pP)];
set asPos [expr ($a_memPos)/(1.0+$n*(1.0-$a_memPos))];
set a_memNeg [expr ($n+1.0)*($MycolNeg*($McMyNeg-1.0)) / ($Ks_col*$th_pN)];
set asNeg [expr ($a_memNeg)/(1.0+$n*(1.0-$a_memNeg))];
rotSpring2DModIKModel 50111 11 117 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
rotSpring2DModIKModel 50112 215 216 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
set Ks_col 5.31913e+08;
set MycolPos 257766;
set MycolNeg 257766;
set McMyPos 1.1;
set McMyNeg 1.1;
set LS 1000.0;
set LK 1000.0;
set LA 1000.0;
set LD 1000.0;
set cS 1.0;
set cK 1.0;
set cA 1.0;
set cD 1.0;
set th_pP 0.232571;
set th_pN 0.232571;
set th_pcP 1.33699;
set th_pcN 1.33699;
set ResP 0.4;
set ResN 0.4;
set th_uP 0.55;
set th_uN 0.55;
set DP 1.0;
set DN 1.0;
set a_memPos [expr ($n+1.0)*($MycolPos*($McMyPos-1.0)) / ($Ks_col*$th_pP)];
set asPos [expr ($a_memPos)/(1.0+$n*(1.0-$a_memPos))];
set a_memNeg [expr ($n+1.0)*($MycolNeg*($McMyNeg-1.0)) / ($Ks_col*$th_pN)];
set asNeg [expr ($a_memNeg)/(1.0+$n*(1.0-$a_memNeg))];
rotSpring2DModIKModel 50121 12 127 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
rotSpring2DModIKModel 50122 225 226 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
set Ks_col 5.31913e+08;
set MycolPos 257766;
set MycolNeg 257766;
set McMyPos 1.1;
set McMyNeg 1.1;
set LS 1000.0;
set LK 1000.0;
set LA 1000.0;
set LD 1000.0;
set cS 1.0;
set cK 1.0;
set cA 1.0;
set cD 1.0;
set th_pP 0.232571;
set th_pN 0.232571;
set th_pcP 1.33699;
set th_pcN 1.33699;
set ResP 0.4;
set ResN 0.4;
set th_uP 0.55;
set th_uN 0.55;
set DP 1.0;
set DN 1.0;
set a_memPos [expr ($n+1.0)*($MycolPos*($McMyPos-1.0)) / ($Ks_col*$th_pP)];
set asPos [expr ($a_memPos)/(1.0+$n*(1.0-$a_memPos))];
set a_memNeg [expr ($n+1.0)*($MycolNeg*($McMyNeg-1.0)) / ($Ks_col*$th_pN)];
set asNeg [expr ($a_memNeg)/(1.0+$n*(1.0-$a_memNeg))];
rotSpring2DModIKModel 50131 13 137 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
rotSpring2DModIKModel 50132 235 236 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
set Ks_col 7.56347e+08;
set MycolPos 257766;
set MycolNeg 257766;
set McMyPos 1.1;
set McMyNeg 1.1;
set LS 1000.0;
set LK 1000.0;
set LA 1000.0;
set LD 1000.0;
set cS 1.0;
set cK 1.0;
set cA 1.0;
set cD 1.0;
set th_pP 0.232571;
set th_pN 0.232571;
set th_pcP 1.33699;
set th_pcN 1.33699;
set ResP 0.4;
set ResN 0.4;
set th_uP 0.55;
set th_uN 0.55;
set DP 1.0;
set DN 1.0;
set a_memPos [expr ($n+1.0)*($MycolPos*($McMyPos-1.0)) / ($Ks_col*$th_pP)];
set asPos [expr ($a_memPos)/(1.0+$n*(1.0-$a_memPos))];
set a_memNeg [expr ($n+1.0)*($MycolNeg*($McMyNeg-1.0)) / ($Ks_col*$th_pN)];
set asNeg [expr ($a_memNeg)/(1.0+$n*(1.0-$a_memNeg))];
rotSpring2DModIKModel 50211 217 218 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
rotSpring2DModIKModel 50212 315 316 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
set Ks_col 7.56347e+08;
set MycolPos 257766;
set MycolNeg 257766;
set McMyPos 1.1;
set McMyNeg 1.1;
set LS 1000.0;
set LK 1000.0;
set LA 1000.0;
set LD 1000.0;
set cS 1.0;
set cK 1.0;
set cA 1.0;
set cD 1.0;
set th_pP 0.232571;
set th_pN 0.232571;
set th_pcP 1.33699;
set th_pcN 1.33699;
set ResP 0.4;
set ResN 0.4;
set th_uP 0.55;
set th_uN 0.55;
set DP 1.0;
set DN 1.0;
set a_memPos [expr ($n+1.0)*($MycolPos*($McMyPos-1.0)) / ($Ks_col*$th_pP)];
set asPos [expr ($a_memPos)/(1.0+$n*(1.0-$a_memPos))];
set a_memNeg [expr ($n+1.0)*($MycolNeg*($McMyNeg-1.0)) / ($Ks_col*$th_pN)];
set asNeg [expr ($a_memNeg)/(1.0+$n*(1.0-$a_memNeg))];
rotSpring2DModIKModel 50221 227 228 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
rotSpring2DModIKModel 50222 325 326 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
set Ks_col 7.56347e+08;
set MycolPos 257766;
set MycolNeg 257766;
set McMyPos 1.1;
set McMyNeg 1.1;
set LS 1000.0;
set LK 1000.0;
set LA 1000.0;
set LD 1000.0;
set cS 1.0;
set cK 1.0;
set cA 1.0;
set cD 1.0;
set th_pP 0.232571;
set th_pN 0.232571;
set th_pcP 1.33699;
set th_pcN 1.33699;
set ResP 0.4;
set ResN 0.4;
set th_uP 0.55;
set th_uN 0.55;
set DP 1.0;
set DN 1.0;
set a_memPos [expr ($n+1.0)*($MycolPos*($McMyPos-1.0)) / ($Ks_col*$th_pP)];
set asPos [expr ($a_memPos)/(1.0+$n*(1.0-$a_memPos))];
set a_memNeg [expr ($n+1.0)*($MycolNeg*($McMyNeg-1.0)) / ($Ks_col*$th_pN)];
set asNeg [expr ($a_memNeg)/(1.0+$n*(1.0-$a_memNeg))];
rotSpring2DModIKModel 50231 237 238 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
rotSpring2DModIKModel 50232 335 336 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
set Ks_col 3.95837e+08;
set MycolPos 139517;
set MycolNeg 139517;
set McMyPos 1.1;
set McMyNeg 1.1;
set LS 1000.0;
set LK 1000.0;
set LA 1000.0;
set LD 1000.0;
set cS 1.0;
set cK 1.0;
set cA 1.0;
set cD 1.0;
set th_pP 0.197059;
set th_pN 0.197059;
set th_pcP 1.08941;
set th_pcN 1.08941;
set ResP 0.4;
set ResN 0.4;
set th_uP 0.55;
set th_uN 0.55;
set DP 1.0;
set DN 1.0;
set a_memPos [expr ($n+1.0)*($MycolPos*($McMyPos-1.0)) / ($Ks_col*$th_pP)];
set asPos [expr ($a_memPos)/(1.0+$n*(1.0-$a_memPos))];
set a_memNeg [expr ($n+1.0)*($MycolNeg*($McMyNeg-1.0)) / ($Ks_col*$th_pN)];
set asNeg [expr ($a_memNeg)/(1.0+$n*(1.0-$a_memNeg))];
rotSpring2DModIKModel 50311 317 318 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
rotSpring2DModIKModel 50312 415 416 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
set Ks_col 3.95837e+08;
set MycolPos 139517;
set MycolNeg 139517;
set McMyPos 1.1;
set McMyNeg 1.1;
set LS 1000.0;
set LK 1000.0;
set LA 1000.0;
set LD 1000.0;
set cS 1.0;
set cK 1.0;
set cA 1.0;
set cD 1.0;
set th_pP 0.197059;
set th_pN 0.197059;
set th_pcP 1.08941;
set th_pcN 1.08941;
set ResP 0.4;
set ResN 0.4;
set th_uP 0.55;
set th_uN 0.55;
set DP 1.0;
set DN 1.0;
set a_memPos [expr ($n+1.0)*($MycolPos*($McMyPos-1.0)) / ($Ks_col*$th_pP)];
set asPos [expr ($a_memPos)/(1.0+$n*(1.0-$a_memPos))];
set a_memNeg [expr ($n+1.0)*($MycolNeg*($McMyNeg-1.0)) / ($Ks_col*$th_pN)];
set asNeg [expr ($a_memNeg)/(1.0+$n*(1.0-$a_memNeg))];
rotSpring2DModIKModel 50321 327 328 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
rotSpring2DModIKModel 50322 425 426 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
set Ks_col 3.95837e+08;
set MycolPos 139517;
set MycolNeg 139517;
set McMyPos 1.1;
set McMyNeg 1.1;
set LS 1000.0;
set LK 1000.0;
set LA 1000.0;
set LD 1000.0;
set cS 1.0;
set cK 1.0;
set cA 1.0;
set cD 1.0;
set th_pP 0.197059;
set th_pN 0.197059;
set th_pcP 1.08941;
set th_pcN 1.08941;
set ResP 0.4;
set ResN 0.4;
set th_uP 0.55;
set th_uN 0.55;
set DP 1.0;
set DN 1.0;
set a_memPos [expr ($n+1.0)*($MycolPos*($McMyPos-1.0)) / ($Ks_col*$th_pP)];
set asPos [expr ($a_memPos)/(1.0+$n*(1.0-$a_memPos))];
set a_memNeg [expr ($n+1.0)*($MycolNeg*($McMyNeg-1.0)) / ($Ks_col*$th_pN)];
set asNeg [expr ($a_memNeg)/(1.0+$n*(1.0-$a_memNeg))];
rotSpring2DModIKModel 50331 337 338 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
rotSpring2DModIKModel 50332 435 436 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
set Ks_col 3.95837e+08;
set MycolPos 139517;
set MycolNeg 139517;
set McMyPos 1.1;
set McMyNeg 1.1;
set LS 1000.0;
set LK 1000.0;
set LA 1000.0;
set LD 1000.0;
set cS 1.0;
set cK 1.0;
set cA 1.0;
set cD 1.0;
set th_pP 0.197059;
set th_pN 0.197059;
set th_pcP 1.08941;
set th_pcN 1.08941;
set ResP 0.4;
set ResN 0.4;
set th_uP 0.55;
set th_uN 0.55;
set DP 1.0;
set DN 1.0;
set a_memPos [expr ($n+1.0)*($MycolPos*($McMyPos-1.0)) / ($Ks_col*$th_pP)];
set asPos [expr ($a_memPos)/(1.0+$n*(1.0-$a_memPos))];
set a_memNeg [expr ($n+1.0)*($MycolNeg*($McMyNeg-1.0)) / ($Ks_col*$th_pN)];
set asNeg [expr ($a_memNeg)/(1.0+$n*(1.0-$a_memNeg))];
rotSpring2DModIKModel 50411 417 418 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
rotSpring2DModIKModel 50412 515 516 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
set Ks_col 3.95837e+08;
set MycolPos 139517;
set MycolNeg 139517;
set McMyPos 1.1;
set McMyNeg 1.1;
set LS 1000.0;
set LK 1000.0;
set LA 1000.0;
set LD 1000.0;
set cS 1.0;
set cK 1.0;
set cA 1.0;
set cD 1.0;
set th_pP 0.197059;
set th_pN 0.197059;
set th_pcP 1.08941;
set th_pcN 1.08941;
set ResP 0.4;
set ResN 0.4;
set th_uP 0.55;
set th_uN 0.55;
set DP 1.0;
set DN 1.0;
set a_memPos [expr ($n+1.0)*($MycolPos*($McMyPos-1.0)) / ($Ks_col*$th_pP)];
set asPos [expr ($a_memPos)/(1.0+$n*(1.0-$a_memPos))];
set a_memNeg [expr ($n+1.0)*($MycolNeg*($McMyNeg-1.0)) / ($Ks_col*$th_pN)];
set asNeg [expr ($a_memNeg)/(1.0+$n*(1.0-$a_memNeg))];
rotSpring2DModIKModel 50421 427 428 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
rotSpring2DModIKModel 50422 525 526 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
set Ks_col 3.95837e+08;
set MycolPos 139517;
set MycolNeg 139517;
set McMyPos 1.1;
set McMyNeg 1.1;
set LS 1000.0;
set LK 1000.0;
set LA 1000.0;
set LD 1000.0;
set cS 1.0;
set cK 1.0;
set cA 1.0;
set cD 1.0;
set th_pP 0.197059;
set th_pN 0.197059;
set th_pcP 1.08941;
set th_pcN 1.08941;
set ResP 0.4;
set ResN 0.4;
set th_uP 0.55;
set th_uN 0.55;
set DP 1.0;
set DN 1.0;
set a_memPos [expr ($n+1.0)*($MycolPos*($McMyPos-1.0)) / ($Ks_col*$th_pP)];
set asPos [expr ($a_memPos)/(1.0+$n*(1.0-$a_memPos))];
set a_memNeg [expr ($n+1.0)*($MycolNeg*($McMyNeg-1.0)) / ($Ks_col*$th_pN)];
set asNeg [expr ($a_memNeg)/(1.0+$n*(1.0-$a_memNeg))];
rotSpring2DModIKModel 50431 437 438 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
rotSpring2DModIKModel 50432 535 536 $Ks_col $asPos $asNeg $MycolPos [expr -$MycolNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
region 1 -ele 50111 50112 50121 50122 50131 50132 50211 50212 50221 50222 50231 50232 50311 50312 50321 50322 50331 50332 50411 50412 50421 50422 50431 50432;
set Ks_beam 2.62683e+08;
set MybeamPos 287116;
set MybeamNeg 265848;
set McMyPos 1.3;
set McMyNeg 1.05;
set LS 1000.0;
set LK 1000.0;
set LA 1000.0;
set LD 1000.0;
set cS 1.0;
set cK 1.0;
set cA 1.0;
set cD 1.0;
set th_pP 0.444857;
set th_pN 0.234785;
set th_pcP 3.18636;
set th_pcN 1.68169;
set ResP 0.3;
set ResN 0.2;
set th_uP 0.65;
set th_uN 0.65;
set DP 1.0;
set DN 1.0;
set a_memPos [expr ($n+1.0)*($MybeamPos*($McMyPos-1.0)) / ($Ks_beam*$th_pP)];
set asPos [expr ($a_memPos)/(1.0+$n*(1.0-$a_memPos))];
set a_memNeg [expr ($n+1.0)*($MybeamNeg*($McMyNeg-1.0)) / ($Ks_beam*$th_pN)];
set asNeg [expr ($a_memNeg)/(1.0+$n*(1.0-$a_memNeg))];
rotSpring2DModIKModel 60211 211 212 $Ks_beam $asPos $asNeg $MybeamPos [expr -$MybeamNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
rotSpring2DModIKModel 60212 223 224 $Ks_beam $asPos $asNeg $MybeamPos [expr -$MybeamNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
set Ks_beam 2.62683e+08;
set MybeamPos 287116;
set MybeamNeg 265848;
set McMyPos 1.3;
set McMyNeg 1.05;
set LS 1000.0;
set LK 1000.0;
set LA 1000.0;
set LD 1000.0;
set cS 1.0;
set cK 1.0;
set cA 1.0;
set cD 1.0;
set th_pP 0.444857;
set th_pN 0.234785;
set th_pcP 3.18636;
set th_pcN 1.68169;
set ResP 0.3;
set ResN 0.2;
set th_uP 0.65;
set th_uN 0.65;
set DP 1.0;
set DN 1.0;
set a_memPos [expr ($n+1.0)*($MybeamPos*($McMyPos-1.0)) / ($Ks_beam*$th_pP)];
set asPos [expr ($a_memPos)/(1.0+$n*(1.0-$a_memPos))];
set a_memNeg [expr ($n+1.0)*($MybeamNeg*($McMyNeg-1.0)) / ($Ks_beam*$th_pN)];
set asNeg [expr ($a_memNeg)/(1.0+$n*(1.0-$a_memNeg))];
rotSpring2DModIKModel 60221 221 222 $Ks_beam $asPos $asNeg $MybeamPos [expr -$MybeamNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
rotSpring2DModIKModel 60222 233 234 $Ks_beam $asPos $asNeg $MybeamPos [expr -$MybeamNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
set Ks_beam 2.62683e+08;
set MybeamPos 287116;
set MybeamNeg 265848;
set McMyPos 1.3;
set McMyNeg 1.05;
set LS 1000.0;
set LK 1000.0;
set LA 1000.0;
set LD 1000.0;
set cS 1.0;
set cK 1.0;
set cA 1.0;
set cD 1.0;
set th_pP 0.444857;
set th_pN 0.234785;
set th_pcP 3.18636;
set th_pcN 1.68169;
set ResP 0.3;
set ResN 0.2;
set th_uP 0.65;
set th_uN 0.65;
set DP 1.0;
set DN 1.0;
set a_memPos [expr ($n+1.0)*($MybeamPos*($McMyPos-1.0)) / ($Ks_beam*$th_pP)];
set asPos [expr ($a_memPos)/(1.0+$n*(1.0-$a_memPos))];
set a_memNeg [expr ($n+1.0)*($MybeamNeg*($McMyNeg-1.0)) / ($Ks_beam*$th_pN)];
set asNeg [expr ($a_memNeg)/(1.0+$n*(1.0-$a_memNeg))];
rotSpring2DModIKModel 60311 311 312 $Ks_beam $asPos $asNeg $MybeamPos [expr -$MybeamNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
rotSpring2DModIKModel 60312 323 324 $Ks_beam $asPos $asNeg $MybeamPos [expr -$MybeamNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
set Ks_beam 2.62683e+08;
set MybeamPos 287116;
set MybeamNeg 265848;
set McMyPos 1.3;
set McMyNeg 1.05;
set LS 1000.0;
set LK 1000.0;
set LA 1000.0;
set LD 1000.0;
set cS 1.0;
set cK 1.0;
set cA 1.0;
set cD 1.0;
set th_pP 0.444857;
set th_pN 0.234785;
set th_pcP 3.18636;
set th_pcN 1.68169;
set ResP 0.3;
set ResN 0.2;
set th_uP 0.65;
set th_uN 0.65;
set DP 1.0;
set DN 1.0;
set a_memPos [expr ($n+1.0)*($MybeamPos*($McMyPos-1.0)) / ($Ks_beam*$th_pP)];
set asPos [expr ($a_memPos)/(1.0+$n*(1.0-$a_memPos))];
set a_memNeg [expr ($n+1.0)*($MybeamNeg*($McMyNeg-1.0)) / ($Ks_beam*$th_pN)];
set asNeg [expr ($a_memNeg)/(1.0+$n*(1.0-$a_memNeg))];
rotSpring2DModIKModel 60321 321 322 $Ks_beam $asPos $asNeg $MybeamPos [expr -$MybeamNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
rotSpring2DModIKModel 60322 333 334 $Ks_beam $asPos $asNeg $MybeamPos [expr -$MybeamNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
set Ks_beam 1.49955e+08;
set MybeamPos 207872;
set MybeamNeg 192474;
set McMyPos 1.3;
set McMyNeg 1.05;
set LS 1000.0;
set LK 1000.0;
set LA 1000.0;
set LD 1000.0;
set cS 1.0;
set cK 1.0;
set cA 1.0;
set cD 1.0;
set th_pP 0.609666;
set th_pN 0.321768;
set th_pcP 4.86636;
set th_pcN 2.56836;
set ResP 0.3;
set ResN 0.2;
set th_uP 0.65;
set th_uN 0.65;
set DP 1.0;
set DN 1.0;
set a_memPos [expr ($n+1.0)*($MybeamPos*($McMyPos-1.0)) / ($Ks_beam*$th_pP)];
set asPos [expr ($a_memPos)/(1.0+$n*(1.0-$a_memPos))];
set a_memNeg [expr ($n+1.0)*($MybeamNeg*($McMyNeg-1.0)) / ($Ks_beam*$th_pN)];
set asNeg [expr ($a_memNeg)/(1.0+$n*(1.0-$a_memNeg))];
rotSpring2DModIKModel 60411 411 412 $Ks_beam $asPos $asNeg $MybeamPos [expr -$MybeamNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
rotSpring2DModIKModel 60412 423 424 $Ks_beam $asPos $asNeg $MybeamPos [expr -$MybeamNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
set Ks_beam 1.49955e+08;
set MybeamPos 207872;
set MybeamNeg 192474;
set McMyPos 1.3;
set McMyNeg 1.05;
set LS 1000.0;
set LK 1000.0;
set LA 1000.0;
set LD 1000.0;
set cS 1.0;
set cK 1.0;
set cA 1.0;
set cD 1.0;
set th_pP 0.609666;
set th_pN 0.321768;
set th_pcP 4.86636;
set th_pcN 2.56836;
set ResP 0.3;
set ResN 0.2;
set th_uP 0.65;
set th_uN 0.65;
set DP 1.0;
set DN 1.0;
set a_memPos [expr ($n+1.0)*($MybeamPos*($McMyPos-1.0)) / ($Ks_beam*$th_pP)];
set asPos [expr ($a_memPos)/(1.0+$n*(1.0-$a_memPos))];
set a_memNeg [expr ($n+1.0)*($MybeamNeg*($McMyNeg-1.0)) / ($Ks_beam*$th_pN)];
set asNeg [expr ($a_memNeg)/(1.0+$n*(1.0-$a_memNeg))];
rotSpring2DModIKModel 60421 421 422 $Ks_beam $asPos $asNeg $MybeamPos [expr -$MybeamNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
rotSpring2DModIKModel 60422 433 434 $Ks_beam $asPos $asNeg $MybeamPos [expr -$MybeamNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
set Ks_beam 1.49955e+08;
set MybeamPos 207872;
set MybeamNeg 192474;
set McMyPos 1.3;
set McMyNeg 1.05;
set LS 1000.0;
set LK 1000.0;
set LA 1000.0;
set LD 1000.0;
set cS 1.0;
set cK 1.0;
set cA 1.0;
set cD 1.0;
set th_pP 0.609666;
set th_pN 0.321768;
set th_pcP 4.86636;
set th_pcN 2.56836;
set ResP 0.3;
set ResN 0.2;
set th_uP 0.65;
set th_uN 0.65;
set DP 1.0;
set DN 1.0;
set a_memPos [expr ($n+1.0)*($MybeamPos*($McMyPos-1.0)) / ($Ks_beam*$th_pP)];
set asPos [expr ($a_memPos)/(1.0+$n*(1.0-$a_memPos))];
set a_memNeg [expr ($n+1.0)*($MybeamNeg*($McMyNeg-1.0)) / ($Ks_beam*$th_pN)];
set asNeg [expr ($a_memNeg)/(1.0+$n*(1.0-$a_memNeg))];
rotSpring2DModIKModel 60511 511 512 $Ks_beam $asPos $asNeg $MybeamPos [expr -$MybeamNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
rotSpring2DModIKModel 60512 523 524 $Ks_beam $asPos $asNeg $MybeamPos [expr -$MybeamNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
set Ks_beam 1.49955e+08;
set MybeamPos 207872;
set MybeamNeg 192474;
set McMyPos 1.3;
set McMyNeg 1.05;
set LS 1000.0;
set LK 1000.0;
set LA 1000.0;
set LD 1000.0;
set cS 1.0;
set cK 1.0;
set cA 1.0;
set cD 1.0;
set th_pP 0.609666;
set th_pN 0.321768;
set th_pcP 4.86636;
set th_pcN 2.56836;
set ResP 0.3;
set ResN 0.2;
set th_uP 0.65;
set th_uN 0.65;
set DP 1.0;
set DN 1.0;
set a_memPos [expr ($n+1.0)*($MybeamPos*($McMyPos-1.0)) / ($Ks_beam*$th_pP)];
set asPos [expr ($a_memPos)/(1.0+$n*(1.0-$a_memPos))];
set a_memNeg [expr ($n+1.0)*($MybeamNeg*($McMyNeg-1.0)) / ($Ks_beam*$th_pN)];
set asNeg [expr ($a_memNeg)/(1.0+$n*(1.0-$a_memNeg))];
rotSpring2DModIKModel 60521 521 522 $Ks_beam $asPos $asNeg $MybeamPos [expr -$MybeamNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
rotSpring2DModIKModel 60522 533 534 $Ks_beam $asPos $asNeg $MybeamPos [expr -$MybeamNeg] $LS $LK $LA $LD $cS $cK $cA $cD $th_pP $th_pN $th_pcP $th_pcN $ResP $ResN $th_uP $th_uN $DP $DN;
region 2 -ele 60211 60212 60221 60222 60311 60312 60321 60322 60411 60412 60421 60422 60511 60512 60521 60522;
source rotPanelZone2D.tcl;
set Fy 3.515348e+07;
set Ry 1.100000e+00;
set as_pz 3.000000e-02;
set dc 622;
set bf_c 328;
set tf_c 24.4;
set tp 15.4;
set db 688;
rotPanelZone2D 7021 2103 2104 $Es $Fy $dc $bf_c $tf_c $tp $db $Ry $as_pz;
set dc 622;
set bf_c 328;
set tf_c 24.4;
set tp 15.4;
set db 688;
rotPanelZone2D 7022 2203 2204 $Es $Fy $dc $bf_c $tf_c $tp $db $Ry $as_pz;
set dc 622;
set bf_c 328;
set tf_c 24.4;
set tp 15.4;
set db 688;
rotPanelZone2D 7023 2303 2304 $Es $Fy $dc $bf_c $tf_c $tp $db $Ry $as_pz;
set dc 622;
set bf_c 328;
set tf_c 24.4;
set tp 15.4;
set db 688;
rotPanelZone2D 7031 3103 3104 $Es $Fy $dc $bf_c $tf_c $tp $db $Ry $as_pz;
set dc 622;
set bf_c 328;
set tf_c 24.4;
set tp 15.4;
set db 688;
rotPanelZone2D 7032 3203 3204 $Es $Fy $dc $bf_c $tf_c $tp $db $Ry $as_pz;
set dc 622;
set bf_c 328;
set tf_c 24.4;
set tp 15.4;
set db 688;
rotPanelZone2D 7033 3303 3304 $Es $Fy $dc $bf_c $tf_c $tp $db $Ry $as_pz;
set dc 607;
set bf_c 228;
set tf_c 17.3;
set tp 11.2;
set db 549;
rotPanelZone2D 7041 4103 4104 $Es $Fy $dc $bf_c $tf_c $tp $db $Ry $as_pz;
set dc 607;
set bf_c 228;
set tf_c 17.3;
set tp 11.2;
set db 549;
rotPanelZone2D 7042 4203 4204 $Es $Fy $dc $bf_c $tf_c $tp $db $Ry $as_pz;
set dc 607;
set bf_c 228;
set tf_c 17.3;
set tp 11.2;
set db 549;
rotPanelZone2D 7043 4303 4304 $Es $Fy $dc $bf_c $tf_c $tp $db $Ry $as_pz;
set dc 607;
set bf_c 228;
set tf_c 17.3;
set tp 11.2;
set db 549;
rotPanelZone2D 7051 5103 5104 $Es $Fy $dc $bf_c $tf_c $tp $db $Ry $as_pz;
set dc 607;
set bf_c 228;
set tf_c 17.3;
set tp 11.2;
set db 549;
rotPanelZone2D 7052 5203 5204 $Es $Fy $dc $bf_c $tf_c $tp $db $Ry $as_pz;
set dc 607;
set bf_c 228;
set tf_c 17.3;
set tp 11.2;
set db 549;
rotPanelZone2D 7053 5303 5304 $Es $Fy $dc $bf_c $tf_c $tp $db $Ry $as_pz;
rotLeaningCol 80142 246 24;
rotLeaningCol 80241 24 247;
rotLeaningCol 80242 346 34;
rotLeaningCol 80341 34 347;
rotLeaningCol 80342 446 44;
rotLeaningCol 80441 44 447;
rotLeaningCol 80442 546 54;
region 3 -ele 80142 80241 80242 80341 80342 80441 80442;
set eigenValues "EigenValuesNA.txt"
set EigenValues [open $eigenValues "w"]
puts $EigenValues "[eigen 10]"
close $EigenValues
pattern Plain 101 Constant {
load 216 0.0 -75182.4 0.0;
load 226 0.0 -75181.4 0.0;
load 236 0.0 -75180.4 0.0;
load 316 0.0 -74943.6 0.0;
load 326 0.0 -74942.6 0.0;
load 336 0.0 -74941.6 0.0;
load 416 0.0 -74708.3 0.0;
load 426 0.0 -74707.3 0.0;
load 436 0.0 -74706.3 0.0;
load 516 0.0 -74501.5 0.0;
load 526 0.0 -74500.5 0.0;
load 536 0.0 -74499.5 0.0;
load 24 0.0 -97939.7 0.0;
load 34 0.0 -97939.7 0.0;
load 44 0.0 -97939.7 0.0;
load 54 0.0 -97939.7 0.0;
}
recorder Drift -file $dataDir/StoryDrift1.txt -iNode 11 -jNode 2105 -dof 1 -perpDirn 2;
recorder Drift -file $dataDir/StoryDrift2.txt -iNode 2105 -jNode 3105 -dof 1 -perpDirn 2;
recorder Drift -file $dataDir/StoryDrift3.txt -iNode 3105 -jNode 4105 -dof 1 -perpDirn 2;
recorder Drift -file $dataDir/StoryDrift4.txt -iNode 4105 -jNode 5105 -dof 1 -perpDirn 2;
recorder Drift -file $dataDir/Drift.txt -time -iNode 11 2105 3105 4105 -jNode 2105 3105 4105 5105 -dof 1 -perpDirn 2;
recorder Drift -file $dataDir/RoofDisp.txt -iNode 12 -jNode 5205 -dof 1 -perpDirn 2;
recorder Node -file $dataDir/Vbase.txt -node 117 127 137 14 -dof 1 reaction;
set Tol 1.0e-6;
constraints Plain;
numberer RCM;
system BandGeneral ;
test NormDispIncr $Tol 6 ;
algorithm Newton;
set NstepGravity 10;
set DGravity [expr 1./$NstepGravity];
integrator LoadControl $DGravity;
analysis Static;
analyze $NstepGravity;
loadConst -time 0.0
puts "Model Built"
puts "Running dynamic analysis..."
set ViewScale 5;
set pi [expr 2.0*asin(1.0)];
set nEigenI 1;
set nEigenJ 2;
set lambdaN [eigen [expr $nEigenJ]];
set lambdaI [lindex $lambdaN [expr 0]];
set lambdaJ [lindex $lambdaN [expr $nEigenJ-1]];
set w1 [expr pow($lambdaI,0.5)];
set w2 [expr pow($lambdaJ,0.5)];
set T1 [expr 2.0*$pi/$w1];
set T2 [expr 2.0*$pi/$w2];
set zeta 0.05;
set a0 [expr $zeta*2.0*$w1*$w2/($w1 + $w2)];
set a1 [expr $zeta*2.0/($w1 + $w2)];
set a1_mod [expr $a1*(1.0+$n)/$n];
region 4 -ele 2021 2022 2031 2032 2041 2042 2051 2052 1011 1012 1013 1021 1022 1023 1031 1032 1033 1041 1042 1043 rayleigh 0.0 0.0 $a1_mod 0.0;
region 5 -ele 20211 20212 20221 20222 20311 20312 20321 20322 20411 20412 20421 20422 20511 20512 20521 20522 rayleigh 0.0 0.0 $a1 0.0;
rayleigh $a0 0.0 0.0 0.0;
source ReadSMDfile.tcl;
set GMdirection 1;
set GMfile "SylmarX" ;
set GMfact "29.418";
set dt 0.005;
set TmaxAnalysis 40;
set IDloadTag 400;
set inFile $GMdir$GMfile.at2
set outFile $GMdir$GMfile.g3;
ReadSMDFile $inFile $outFile dt;
set AccelSeries "Series -dt $dt -filePath $outFile -factor $GMfact";
pattern UniformExcitation $IDloadTag $GMdirection -accel $AccelSeries ;
wipeAnalysis;
set dt_analysis 0.005
set TolDynamic 1.0e-8;
constraints Plain;
numberer RCM;
system BandGeneral;
set NewmarkGamma 0.5;
set NewmarkBeta 0.25;
set maxNumIter 1000;
set moreIter [expr $maxNumIter*10];
set printFlag 0;
set TestType NormDispIncr;
set algorithmType Newton;
test $TestType $TolDynamic $maxNumIter;
algorithm $algorithmType -initial ;
integrator Newmark $NewmarkGamma $NewmarkBeta;
analysis Transient;
set Nsteps [expr int($TmaxAnalysis/$dt_analysis)];
for {set ik 1} {$ik <= $Nsteps} {incr ik 1} {
puts "$ik/$Nsteps"
set ok [analyze 1 $dt_analysis]
if {$ok != 0} {
puts "Trying InitialInterpolated ...";
algorithm NewtonLineSearch <-type InitialInterpolated>;
set ok [analyze 1 $dt_analysis]
algorithm $algorithmType -initial;
};
if {$ok != 0} {
puts "Trying Bisection ...";
algorithm NewtonLineSearch <-type Bisection>;
set ok [analyze 1 $dt_analysis]
algorithm $algorithmType -initial;
};
if {$ok != 0} {
puts "Trying Secant ...";
algorithm NewtonLineSearch <-type Secant>;
set ok [analyze 1 $dt_analysis]
algorithm $algorithmType -initial;
};
if {$ok != 0} {
puts "Trying RegulaFalsi ...";
algorithm NewtonLineSearch <-type RegulaFalsi>;
set ok [analyze 1 $dt_analysis]
algorithm $algorithmType -initial;
};
if {$ok != 0} {
puts "Trying KrylovNewton ...";
algorithm KrylovNewton;
set ok [analyze 1 $dt_analysis]
algorithm $algorithmType -initial;
};
if {$ok != 0} {
puts "Trying BFGS ...";
algorithm BFGS;
set ok [analyze 1 $dt_analysis]
algorithm $algorithmType -initial;
};
if {$ok != 0} {
puts "Trying Broyden ...";
algorithm Broyden;
set ok [analyze 1 $dt_analysis]
algorithm $algorithmType -initial;
};
if {$ok != 0} {
puts "Trying BackwardEuler ...";
integrator BackwardEuler;
set ok [analyze 1 $dt_analysis]
integrator Newmark $NewmarkGamma $NewmarkBeta;
};
if {$ok != 0} {
puts "Trying HHT 0.9 ...";
integrator HHT 0.9;
set ok [analyze 1 $dt_analysis]
integrator Newmark $NewmarkGamma $NewmarkBeta;
};
if {$ok != 0} {
puts "Trying OS ...";
integrator AlphaOS 1.00;
algorithm Linear;
set ok [analyze 1 $dt_analysis]
integrator Newmark $NewmarkGamma $NewmarkBeta;
algorithm $algorithmType -initial;
};
if {$ok != 0} {
puts "Trying OSG ...";
integrator AlphaOSGeneralized 1.00;
algorithm Linear;
set ok [analyze 1 $dt_analysis]
integrator Newmark $NewmarkGamma $NewmarkBeta;
algorithm $algorithmType -initial;
};
if {$ok != 0} {
puts "Trying more iterations...";
test $TestType $Tol $moreIter ;
set ok [analyze 1 $dt_analysis]
test $TestType $Tol $maxNumIter ;
};
if {$ok != 0} {
puts "Trying tolerance 1.0e-7 ...";
test $TestType 1.0e-7 $maxNumIter 0;
set ok [analyze 1 $dt_analysis]
test $TestType $Tol $maxNumIter ;
};
if {$ok != 0} {
puts "Trying tolerance 1.0e-6 ...";
test $TestType 1.0e-6 $maxNumIter 0;
set ok [analyze 1 $dt_analysis]
test $TestType $Tol $maxNumIter ;
};
if {$ok != 0} {
puts "Trying tolerance 1.0e-5 ...";
test $TestType 1.0e-5 $maxNumIter 0;
set ok [analyze 1 $dt_analysis]
test $TestType $Tol $maxNumIter ;
};
if {$ok != 0} {
puts "Trying tolerance 1.0e-3 ...";
test $TestType 1.0e-3 $maxNumIter 0;
set ok [analyze 1 $dt_analysis]
test $TestType $Tol $maxNumIter ;
};
if {$ok != 0} {
puts "Trying tolerance 1.0e-3 & dt_analysis/2 ...";
test $TestType 1.0e-3 $maxNumIter 0;
set ok [analyze 1 [expr $dt_analysis/2.0]]
test $TestType $Tol $maxNumIter ;
};
if {$ok != 0} {
puts "Trying tolerance 1.0e-3 & dt_analysis/4 ...";
test $TestType 1.0e-3 $maxNumIter 0;
set ok [analyze 1 [expr $dt_analysis/4.0]]
test $TestType $Tol $maxNumIter ;
};
if {$ok != 0} {
puts "Trying tolerance 1.0e-3 & dt_analysis/8...";
test $TestType 1.0e-3 $maxNumIter 0;
set ok [analyze 1 [expr $dt_analysis/8.0]]
test $TestType $Tol $maxNumIter ;
};
if {$ok != 0} {
puts "Trying tolerance 1.0e-3 & dt_analysis/16 ...";
test $TestType 1.0e-3 $maxNumIter 0;
set ok [analyze 1 [expr $dt_analysis/16.0]]
test $TestType $Tol $maxNumIter ;
};
if {$ok != 0} {
puts "Trying tolerance 1.0e-3 & dt_analysis/32 ...";
test $TestType 1.0e-3 $maxNumIter 0;
set ok [analyze 1 [expr $dt_analysis/32.0]]
test $TestType $Tol $maxNumIter ;
};
if {$ok != 0} {
puts "Trying tolerance 1.0e-3 & dt_analysis/64 ...";
test $TestType 1.0e-3 $maxNumIter 0;
set ok [analyze 1 [expr $dt_analysis/64.0]]
test $TestType $Tol $maxNumIter ;
};
if {$ok != 0} {
set Nstepsmax [expr $ik-1]
break;
}
};
if {$ok == 0} {
puts "Dynamic analysis complete";
} else {
puts "Dynamic analysis did not converge";
}
wipe all;

[selimgunay]

Several questions about your model/analysis:

1. Did you make an eigenvalue analysis? If yes, are the periods reasonable?

2. At what part of the ground motion does the convergence problem happen? Is it there from the beginning?

[Mohsenkalatehahani]

Thanks for your quick reply.

1- Yes. The periods are close to those reported in the previous works; T1= 0.92 sec and T2= 0.30 sec.
2- No, the problem occurs when the response (Drift/Displacement) is too large. Actually, I want to calculate the residual deformation under different scales of a single earthquake record. When the scale factor is small, there is no problem. However, for large ones, I face with convergence problems.

How can we identify when a collapse has occurred? Maybe the frame will experience a sidesway collapse under significant earthquakes!

[apapaso]

The easy way is to check the inter-story drift for unexpectedly large values.
The hard way is to calculate and compare the input energy and the dissipated energy: They are almost equal to each other, until gravity takes over...

[selimgunay]

You are right, if the adaptive strategies cannot provide a solution and if the convergence happens at really large drifts, it may be indicative of the sidesway collapse.

You can conduct a pushover analysis to have some idea about the displacements corresponding to collapse.

[Mohsenkalatehahani]

Thank you Apapaso and Selim.

I performed a pushover analysis in order to find when the base shear becomes zero. The corresponding maximum interstory drift is the collapse drift. The nonlinear dynamic analysis continues beyond the collapse drift and then convergence problems occur.