Hello,
I'm using Bilin material model as structural element, and I ask your advice for my script.
The problem is... values of inter-story drift ratio from my modeling are too small. It is close to zero!
I guess the modeling is too stiff and input data for various parameters of the material model has some problems, but I couldn't find a proper solution.
You can also check the information for this modeling: A.B. Liel (2008), chapter 5 (p. 133~),
https://static1.squarespace.com/static/ ... Thesis.pdf
I attached my script as below:
model BasicBuilder -ndm 2 -ndf 3
# Set parameters for overall model geomerty
# Unit : in, ksi
# Create nodes
# tag X Y
node 15 0.000000 0.000000
node 14 0.000000 180.000000
node 13 0.000000 336.000000
node 12 0.000000 492.000000
node 11 0.000000 648.000000
node 25 300.000000 0.000000
node 24 300.000000 180.000000
node 23 300.000000 336.000000
node 22 300.000000 492.000000
node 21 300.000000 648.000000
node 35 600.000000 0.000000
node 34 600.000000 180.000000
node 33 600.000000 336.000000
node 32 600.000000 492.000000
node 31 600.000000 648.000000
node 45 900.000000 0.000000
node 44 900.000000 180.000000
node 43 900.000000 336.000000
node 42 900.000000 492.000000
node 41 900.000000 648.000000
# Create zerolength nodes for beam and column
# tag X Y
node 11014 0.000000 0.000000
node 11013 0.000000 180.000000
node 11012 0.000000 336.000000
node 11011 0.000000 492.000000
node 12014 0.000000 180.000000
node 12013 0.000000 336.000000
node 12012 0.000000 492.000000
node 12011 0.000000 648.000000
node 11024 300.000000 0.000000
node 11023 300.000000 180.000000
node 11022 300.000000 336.000000
node 11021 300.000000 492.000000
node 12024 300.000000 180.000000
node 12023 300.000000 336.000000
node 12022 300.000000 492.000000
node 12021 300.000000 648.000000
node 11034 600.000000 0.000000
node 11033 600.000000 180.000000
node 11032 600.000000 336.000000
node 11031 600.000000 492.000000
node 12034 600.000000 180.000000
node 12033 600.000000 336.000000
node 12032 600.000000 492.000000
node 12031 600.000000 648.000000
node 11044 900.000000 0.000000
node 11043 900.000000 180.000000
node 11042 900.000000 336.000000
node 11041 900.000000 492.000000
node 12044 900.000000 180.000000
node 12043 900.000000 336.000000
node 12042 900.000000 492.000000
node 12041 900.000000 648.000000
node 21014 0.000000 180.000000
node 22014 300.000000 180.000000
node 21013 0.000000 336.000000
node 22013 300.000000 336.000000
node 21012 0.000000 492.000000
node 22012 300.000000 492.000000
node 21011 0.000000 648.000000
node 22011 300.000000 648.000000
node 21024 300.000000 180.000000
node 22024 600.000000 180.000000
node 21023 300.000000 336.000000
node 22023 600.000000 336.000000
node 21022 300.000000 492.000000
node 22022 600.000000 492.000000
node 21021 300.000000 648.000000
node 22021 600.000000 648.000000
node 21034 600.000000 180.000000
node 22034 900.000000 180.000000
node 21033 600.000000 336.000000
node 22033 900.000000 336.000000
node 21032 600.000000 492.000000
node 22032 900.000000 492.000000
node 21031 600.000000 648.000000
node 22031 900.000000 648.000000
# Fix supports at base of columns
# tag DX DY MZ
fix 15 1 1 1
fix 25 1 1 1
fix 35 1 1 1
fix 45 1 1 1
# Bilin $matTag $K0 $as_Plus $as_Neg $My_Plus $My_Neg $Lamda_S $Lamda_C $Lamda_A $Lamda_K $c_S $c_C $c_A $c_K $theta_p_Plus $theta_p_Neg $theta_pc_Plus $theta_pc_Neg $Res_Pos $Res_Neg $theta_u_Plus $theta_u_Neg $D_Plus $D_Neg
uniaxialMaterial Bilin 9011 1.354620e+07 0.050000 0.050000 5502000000 -5502000000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.032000 0.032000 0.069000 0.069000 1 1 406.266500 406.266500 1 1
uniaxialMaterial Bilin 9012 1.393324e+07 0.050000 0.050000 5935000000 -5935000000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.030000 0.030000 0.061000 0.061000 1 1 406.266500 406.266500 1 1
uniaxialMaterial Bilin 9013 1.470731e+07 0.050000 0.050000 6361000000 -6361000000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.028000 0.028000 0.053000 0.053000 1 1 406.266500 406.266500 1 1
uniaxialMaterial Bilin 9014 1.625544e+07 0.050000 0.050000 8587000000 -8587000000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.035000 0.035000 0.076000 0.076000 1 1 406.266500 406.266500 1 1
uniaxialMaterial Bilin 9021 1.393324e+07 0.050000 0.050000 3113000000 -3113000000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.023000 0.023000 0.044000 0.044000 1 1 406.266500 406.266500 1 1
uniaxialMaterial Bilin 9022 1.664248e+07 0.050000 0.050000 5190000000 -5190000000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.021000 0.021000 0.034000 0.034000 1 1 406.266500 406.266500 1 1
uniaxialMaterial Bilin 9023 1.896468e+07 0.050000 0.050000 7604000000 -7604000000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.025000 0.025000 0.043000 0.043000 1 1 406.266500 406.266500 1 1
uniaxialMaterial Bilin 9024 2.206096e+07 0.050000 0.050000 8167000000 -8167000000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.021000 0.021000 0.033000 0.033000 1 1 406.266500 406.266500 1 1
uniaxialMaterial Bilin 9031 1.393324e+07 0.050000 0.050000 3113000000 -3113000000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.023000 0.023000 0.044000 0.044000 1 1 406.266500 406.266500 1 1
uniaxialMaterial Bilin 9032 1.664248e+07 0.050000 0.050000 5190000000 -5190000000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.021000 0.021000 0.034000 0.034000 1 1 406.266500 406.266500 1 1
uniaxialMaterial Bilin 9033 1.896468e+07 0.050000 0.050000 7604000000 -7604000000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.025000 0.025000 0.043000 0.043000 1 1 406.266500 406.266500 1 1
uniaxialMaterial Bilin 9034 2.206096e+07 0.050000 0.050000 8167000000 -8167000000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.021000 0.021000 0.033000 0.033000 1 1 406.266500 406.266500 1 1
uniaxialMaterial Bilin 9041 1.393324e+07 0.050000 0.050000 5502000000 -5502000000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.032000 0.032000 0.069000 0.069000 1 1 406.266500 406.266500 1 1
uniaxialMaterial Bilin 9042 1.393324e+07 0.050000 0.050000 5935000000 -5935000000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.030000 0.030000 0.061000 0.061000 1 1 406.266500 406.266500 1 1
uniaxialMaterial Bilin 9043 1.470731e+07 0.050000 0.050000 6361000000 -6361000000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.028000 0.028000 0.053000 0.053000 1 1 406.266500 406.266500 1 1
uniaxialMaterial Bilin 9044 1.625544e+07 0.050000 0.050000 8544000000 -8544000000 41.000000 41.000000 41.000000 41.000000 1 1 1 1 0.035000 0.035000 0.076000 0.076000 1 1 406.266500 406.266500 1 1
uniaxialMaterial Bilin 9111 2.976101e+07 0.050000 0.050000 4530000000 -4530000000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.084400 0.084400 0.100000 0.100000 1 1 152.396995 152.396995 1 1
uniaxialMaterial Bilin 9112 2.976101e+07 0.050000 0.050000 4177000000 -4177000000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.072000 0.072000 0.100000 0.100000 1 1 152.396995 152.396995 1 1
uniaxialMaterial Bilin 9113 1.354620e+07 0.050000 0.050000 4392000000 -4392000000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.046000 0.046000 0.095000 0.095000 1 1 152.396995 152.396995 1 1
uniaxialMaterial Bilin 9114 1.354620e+07 0.050000 0.050000 4886000000 -4886000000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.047000 0.047000 0.095000 0.095000 1 1 152.396995 152.396995 1 1
uniaxialMaterial Bilin 9121 2.976101e+07 0.050000 0.050000 4530000000 -4530000000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.074300 0.074300 0.100000 0.100000 1 1 152.396995 152.396995 1 1
uniaxialMaterial Bilin 9122 2.976101e+07 0.050000 0.050000 4177000000 -4177000000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.072000 0.072000 0.100000 0.100000 1 1 152.396995 152.396995 1 1
uniaxialMaterial Bilin 9123 1.354620e+07 0.050000 0.050000 4392000000 -4392000000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.046000 0.046000 0.095000 0.095000 1 1 152.396995 152.396995 1 1
uniaxialMaterial Bilin 9124 1.354620e+07 0.050000 0.050000 4886000000 -4886000000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.047000 0.047000 0.095000 0.095000 1 1 152.396995 152.396995 1 1
uniaxialMaterial Bilin 9131 2.976101e+07 0.050000 0.050000 4530000000 -4530000000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.084400 0.084400 0.100000 0.100000 1 1 152.396995 152.396995 1 1
uniaxialMaterial Bilin 9132 2.976101e+07 0.050000 0.050000 4177000000 -4177000000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.072000 0.072000 0.100000 0.100000 1 1 152.396995 152.396995 1 1
uniaxialMaterial Bilin 9133 1.354620e+07 0.050000 0.050000 4392000000 -4392000000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.046000 0.046000 0.095000 0.095000 1 1 152.396995 152.396995 1 1
uniaxialMaterial Bilin 9134 1.354620e+07 0.050000 0.050000 4886000000 -4886000000 131.000000 131.000000 131.000000 131.000000 1 1 1 1 0.047000 0.047000 0.095000 0.095000 1 1 152.396995 152.396995 1 1
uniaxialMaterial Elastic 9911 1161103.047896
uniaxialMaterial Elastic 9811 504827.412128
section Aggregator 11 9911 P 9011 Mz 9811 Vy
uniaxialMaterial Elastic 9912 1161103.047896
uniaxialMaterial Elastic 9812 504827.412128
section Aggregator 12 9912 P 9012 Mz 9812 Vy
uniaxialMaterial Elastic 9913 1161103.047896
uniaxialMaterial Elastic 9813 504827.412128
section Aggregator 13 9913 P 9013 Mz 9813 Vy
uniaxialMaterial Elastic 9914 1161103.047896
uniaxialMaterial Elastic 9814 504827.412128
section Aggregator 14 9914 P 9014 Mz 9814 Vy
uniaxialMaterial Elastic 9921 1161103.047896
uniaxialMaterial Elastic 9821 504827.412128
section Aggregator 21 9921 P 9021 Mz 9821 Vy
uniaxialMaterial Elastic 9922 1161103.047896
uniaxialMaterial Elastic 9822 504827.412128
section Aggregator 22 9922 P 9022 Mz 9822 Vy
uniaxialMaterial Elastic 9923 1161103.047896
uniaxialMaterial Elastic 9823 504827.412128
section Aggregator 23 9923 P 9023 Mz 9823 Vy
uniaxialMaterial Elastic 9924 1161103.047896
uniaxialMaterial Elastic 9824 504827.412128
section Aggregator 24 9924 P 9024 Mz 9824 Vy
uniaxialMaterial Elastic 9931 1161103.047896
uniaxialMaterial Elastic 9831 504827.412128
section Aggregator 31 9931 P 9031 Mz 9831 Vy
uniaxialMaterial Elastic 9932 1161103.047896
uniaxialMaterial Elastic 9832 504827.412128
section Aggregator 32 9932 P 9032 Mz 9832 Vy
uniaxialMaterial Elastic 9933 1161103.047896
uniaxialMaterial Elastic 9833 504827.412128
section Aggregator 33 9933 P 9033 Mz 9833 Vy
uniaxialMaterial Elastic 9934 1161103.047896
uniaxialMaterial Elastic 9834 504827.412128
section Aggregator 34 9934 P 9034 Mz 9834 Vy
uniaxialMaterial Elastic 9941 1161103.047896
uniaxialMaterial Elastic 9841 504827.412128
section Aggregator 41 9941 P 9041 Mz 9841 Vy
uniaxialMaterial Elastic 9942 1161103.047896
uniaxialMaterial Elastic 9842 504827.412128
section Aggregator 42 9942 P 9042 Mz 9842 Vy
uniaxialMaterial Elastic 9943 1161103.047896
uniaxialMaterial Elastic 9843 504827.412128
section Aggregator 43 9943 P 9043 Mz 9843 Vy
uniaxialMaterial Elastic 9944 1161103.047896
uniaxialMaterial Elastic 9844 504827.412128
section Aggregator 44 9944 P 9044 Mz 9844 Vy
uniaxialMaterial Elastic 9711 1161103.047896
uniaxialMaterial Elastic 9611 504827.412128
section Aggregator 211 9711 P 9111 Mz 9611 Vy
uniaxialMaterial Elastic 9712 1161103.047896
uniaxialMaterial Elastic 9612 504827.412128
section Aggregator 212 9712 P 9112 Mz 9612 Vy
uniaxialMaterial Elastic 9713 1509433.962264
uniaxialMaterial Elastic 9613 656275.635767
section Aggregator 213 9713 P 9113 Mz 9613 Vy
uniaxialMaterial Elastic 9714 1509433.962264
uniaxialMaterial Elastic 9614 656275.635767
section Aggregator 214 9714 P 9114 Mz 9614 Vy
uniaxialMaterial Elastic 9721 1161103.047896
uniaxialMaterial Elastic 9621 504827.412128
section Aggregator 221 9721 P 9121 Mz 9621 Vy
uniaxialMaterial Elastic 9722 1161103.047896
uniaxialMaterial Elastic 9622 504827.412128
section Aggregator 222 9722 P 9122 Mz 9622 Vy
uniaxialMaterial Elastic 9723 1509433.962264
uniaxialMaterial Elastic 9623 656275.635767
section Aggregator 223 9723 P 9123 Mz 9623 Vy
uniaxialMaterial Elastic 9724 1509433.962264
uniaxialMaterial Elastic 9624 656275.635767
section Aggregator 224 9724 P 9124 Mz 9624 Vy
uniaxialMaterial Elastic 9731 1161103.047896
uniaxialMaterial Elastic 9631 504827.412128
section Aggregator 231 9731 P 9131 Mz 9631 Vy
uniaxialMaterial Elastic 9732 1161103.047896
uniaxialMaterial Elastic 9632 504827.412128
section Aggregator 232 9732 P 9132 Mz 9632 Vy
uniaxialMaterial Elastic 9733 1509433.962264
uniaxialMaterial Elastic 9633 656275.635767
section Aggregator 233 9733 P 9133 Mz 9633 Vy
uniaxialMaterial Elastic 9734 1509433.962264
uniaxialMaterial Elastic 9634 656275.635767
section Aggregator 234 9734 P 9134 Mz 9634 Vy
geomTransf PDelta 1
# Area of element
set Modul 2.902758e+03
# Create the column element
# tag i_node j_node A(mm^2) E(MPa) Iz(mm^4) transfTag
element elasticBeamColumn 11 11011 12011 400 $Modul 1.333333e+04 1
element elasticBeamColumn 12 11012 12012 400 $Modul 1.333333e+04 1
element elasticBeamColumn 13 11013 12013 400 $Modul 1.333333e+04 1
element elasticBeamColumn 14 11014 12014 400 $Modul 1.333333e+04 1
element elasticBeamColumn 21 11021 12021 400 $Modul 1.333333e+04 1
element elasticBeamColumn 22 11022 12022 400 $Modul 1.333333e+04 1
element elasticBeamColumn 23 11023 12023 400 $Modul 1.333333e+04 1
element elasticBeamColumn 24 11024 12024 400 $Modul 1.333333e+04 1
element elasticBeamColumn 31 11031 12031 400 $Modul 1.333333e+04 1
element elasticBeamColumn 32 11032 12032 400 $Modul 1.333333e+04 1
element elasticBeamColumn 33 11033 12033 400 $Modul 1.333333e+04 1
element elasticBeamColumn 34 11034 12034 400 $Modul 1.333333e+04 1
element elasticBeamColumn 41 11041 12041 400 $Modul 1.333333e+04 1
element elasticBeamColumn 42 11042 12042 400 $Modul 1.333333e+04 1
element elasticBeamColumn 43 11043 12043 400 $Modul 1.333333e+04 1
element elasticBeamColumn 44 11044 12044 400 $Modul 1.333333e+04 1
element elasticBeamColumn 2011 21011 22011 400 $Modul 1.333333e+04 1
element elasticBeamColumn 2012 21012 22012 400 $Modul 1.333333e+04 1
element elasticBeamColumn 2013 21013 22013 520 $Modul 2.929333e+04 1
element elasticBeamColumn 2014 21014 22014 520 $Modul 2.929333e+04 1
element elasticBeamColumn 2021 21021 22021 400 $Modul 1.333333e+04 1
element elasticBeamColumn 2022 21022 22022 400 $Modul 1.333333e+04 1
element elasticBeamColumn 2023 21023 22023 520 $Modul 2.929333e+04 1
element elasticBeamColumn 2024 21024 22024 520 $Modul 2.929333e+04 1
element elasticBeamColumn 2031 21031 22031 400 $Modul 1.333333e+04 1
element elasticBeamColumn 2032 21032 22032 400 $Modul 1.333333e+04 1
element elasticBeamColumn 2033 21033 22033 520 $Modul 2.929333e+04 1
element elasticBeamColumn 2034 21034 22034 520 $Modul 2.929333e+04 1
element zeroLengthSection 3011 12 11011 11 -orient 0 1 0 -1 0 0
element zeroLengthSection 4011 12011 11 11 -orient 0 1 0 -1 0 0
element zeroLengthSection 3012 13 11012 12 -orient 0 1 0 -1 0 0
element zeroLengthSection 4012 12012 12 12 -orient 0 1 0 -1 0 0
element zeroLengthSection 3013 14 11013 13 -orient 0 1 0 -1 0 0
element zeroLengthSection 4013 12013 13 13 -orient 0 1 0 -1 0 0
element zeroLengthSection 3014 15 11014 14 -orient 0 1 0 -1 0 0
element zeroLengthSection 4014 12014 14 14 -orient 0 1 0 -1 0 0
element zeroLengthSection 3021 22 11021 21 -orient 0 1 0 -1 0 0
element zeroLengthSection 4021 12021 21 21 -orient 0 1 0 -1 0 0
element zeroLengthSection 3022 23 11022 22 -orient 0 1 0 -1 0 0
element zeroLengthSection 4022 12022 22 22 -orient 0 1 0 -1 0 0
element zeroLengthSection 3023 24 11023 23 -orient 0 1 0 -1 0 0
element zeroLengthSection 4023 12023 23 23 -orient 0 1 0 -1 0 0
element zeroLengthSection 3024 25 11024 24 -orient 0 1 0 -1 0 0
element zeroLengthSection 4024 12024 24 24 -orient 0 1 0 -1 0 0
element zeroLengthSection 3031 32 11031 31 -orient 0 1 0 -1 0 0
element zeroLengthSection 4031 12031 31 31 -orient 0 1 0 -1 0 0
element zeroLengthSection 3032 33 11032 32 -orient 0 1 0 -1 0 0
element zeroLengthSection 4032 12032 32 32 -orient 0 1 0 -1 0 0
element zeroLengthSection 3033 34 11033 33 -orient 0 1 0 -1 0 0
element zeroLengthSection 4033 12033 33 33 -orient 0 1 0 -1 0 0
element zeroLengthSection 3034 35 11034 34 -orient 0 1 0 -1 0 0
element zeroLengthSection 4034 12034 34 34 -orient 0 1 0 -1 0 0
element zeroLengthSection 3041 42 11041 41 -orient 0 1 0 -1 0 0
element zeroLengthSection 4041 12041 41 41 -orient 0 1 0 -1 0 0
element zeroLengthSection 3042 43 11042 42 -orient 0 1 0 -1 0 0
element zeroLengthSection 4042 12042 42 42 -orient 0 1 0 -1 0 0
element zeroLengthSection 3043 44 11043 43 -orient 0 1 0 -1 0 0
element zeroLengthSection 4043 12043 43 43 -orient 0 1 0 -1 0 0
element zeroLengthSection 3044 45 11044 44 -orient 0 1 0 -1 0 0
element zeroLengthSection 4044 12044 44 44 -orient 0 1 0 -1 0 0
element zeroLengthSection 5011 11 21011 211 -orient 0 1 0 -1 0 0
element zeroLengthSection 6011 22011 21 211 -orient 0 1 0 -1 0 0
element zeroLengthSection 5012 12 21012 212 -orient 0 1 0 -1 0 0
element zeroLengthSection 6012 22012 22 212 -orient 0 1 0 -1 0 0
element zeroLengthSection 5013 13 21013 213 -orient 0 1 0 -1 0 0
element zeroLengthSection 6013 22013 23 213 -orient 0 1 0 -1 0 0
element zeroLengthSection 5014 14 21014 214 -orient 0 1 0 -1 0 0
element zeroLengthSection 6014 22014 24 214 -orient 0 1 0 -1 0 0
element zeroLengthSection 5021 21 21021 221 -orient 0 1 0 -1 0 0
element zeroLengthSection 6021 22021 31 221 -orient 0 1 0 -1 0 0
element zeroLengthSection 5022 22 21022 222 -orient 0 1 0 -1 0 0
element zeroLengthSection 6022 22022 32 222 -orient 0 1 0 -1 0 0
element zeroLengthSection 5023 23 21023 223 -orient 0 1 0 -1 0 0
element zeroLengthSection 6023 22023 33 223 -orient 0 1 0 -1 0 0
element zeroLengthSection 5024 24 21024 224 -orient 0 1 0 -1 0 0
element zeroLengthSection 6024 22024 34 224 -orient 0 1 0 -1 0 0
element zeroLengthSection 5031 31 21031 231 -orient 0 1 0 -1 0 0
element zeroLengthSection 6031 22031 41 231 -orient 0 1 0 -1 0 0
element zeroLengthSection 5032 32 21032 232 -orient 0 1 0 -1 0 0
element zeroLengthSection 6032 22032 42 232 -orient 0 1 0 -1 0 0
element zeroLengthSection 5033 33 21033 233 -orient 0 1 0 -1 0 0
element zeroLengthSection 6033 22033 43 233 -orient 0 1 0 -1 0 0
element zeroLengthSection 5034 34 21034 234 -orient 0 1 0 -1 0 0
element zeroLengthSection 6034 22034 44 234 -orient 0 1 0 -1 0 0
equalDOF 11 21 1
equalDOF 21 31 1
equalDOF 31 41 1
equalDOF 12 22 1
equalDOF 22 32 1
equalDOF 32 42 1
equalDOF 13 23 1
equalDOF 23 33 1
equalDOF 33 43 1
equalDOF 14 24 1
equalDOF 24 34 1
equalDOF 34 44 1
mass 11 2.671905e+01 0 0
mass 12 2.897039e+01 0 0
mass 13 3.032120e+01 0 0
mass 14 3.160273e+01 0 0
mass 21 5.108284e+01 0 0
mass 22 5.333418e+01 0 0
mass 23 5.613971e+01 0 0
mass 24 5.742124e+01 0 0
mass 31 5.108284e+01 0 0
mass 32 5.333418e+01 0 0
mass 33 5.613971e+01 0 0
mass 34 5.742124e+01 0 0
mass 41 2.671905e+01 0 0
mass 42 2.897039e+01 0 0
mass 43 3.032120e+01 0 0
mass 44 3.160273e+01 0 0
set GMfile "Dataset_1"
set d 5
set analysisType "pushover"
set numModes 3
recorder Node -file [format "mode1.out" ] -node 15 14 13 12 11 -dof 1 "eigen 1"
recorder Node -file [format "mode2.out" ] -node 15 14 13 12 11 -dof 1 "eigen 2"
recorder Node -file [format "mode3.out" ] -node 15 14 13 12 11 -dof 1 "eigen 3"
# perform eigen analysis
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)]
}
set PI 3.14159265
set lambda [eigen 3]
set outFileID [open lambda.out w]
puts $outFileID $lambda
close $outFileID
# $nodeTag
pattern Plain 101 Constant {
load 11 0 -10315903.000393 0
load 12 0 -11185120.186858 0
load 13 0 -11706650.498737 0
load 14 0 -12201435.666417 0
load 21 0 -19722471.098023 0
load 22 0 -20591688.284488 0
load 23 0 -21674866.624544 0
load 24 0 -22169651.792224 0
load 31 0 -19722471.098023 0
load 32 0 -20591688.284488 0
load 33 0 -21674866.624544 0
load 34 0 -22169651.792224 0
load 41 0 -10315903.000393 0
load 42 0 -11185120.186858 0
load 43 0 -11706650.498737 0
load 44 0 -12201435.666417 0
}
# Gravity-analysis: load-controlled static analysis
set Tol 1.0e-8; # convergence tolerance for testconstraints $constraintsTypeGravity ; # how it handles boundary conditions
numberer RCM; # renumber dofs to minimize band-width (optimization), if you want to
system BandGeneral ; # how to store and solve the system of equations in the analysis (large model: try UmfPack)
test NormDispIncr $Tol 6 ; # determine if convergence has been achieved at the end of an iteration step
algorithm Newton; # use Newtons solution algorithm: updates tangent stiffness at every iteration
set NstepGravity 10; # apply gravity in 10 steps
set DGravity [expr 1./$NstepGravity]; # first load increment;
integrator LoadControl $DGravity; # determine the next time step for an analysis
analysis Static; # define type of analysis static or transient
analyze $NstepGravity; # apply gravity
loadConst -time 0.0
puts "Model Built"
recorder Drift -file EQDrift.out -time -iNode 12 13 14 15 -jNode 11 12 13 14 -dof 1 -perDirn 2
recorder Node -file node.out -time -node 11 12 13 14 15 -dof 1 disp
recorder Node -file Vbase.out -node 15 25 35 45 14 24 34 44 13 23 33 43 12 22 32 42 11 21 31 41 -dof 1 2 3 reaction
recorder Node -file Disp.out -node 11 12 13 14 15 -dof 1 2 3 disp
set xDamp 0.05;
set MpropSwitch 1.0;
set KcurrSwitch 0.0;
set KcommSwitch 1.0;
set KinitSwitch 0.0;
set nEigenI 1;
set nEigenJ 2;
set lambdaN [eigen [expr $nEigenJ]];
set lambdaI [lindex $lambdaN [expr $nEigenI-1]];
set lambdaJ [lindex $lambdaN [expr $nEigenJ-1]];
set omegaI [expr pow($lambdaI,0.5)];
set omegaJ [expr pow($lambdaJ,0.5)];
set alphaM [expr $MpropSwitch*$xDamp*(2*$omegaI*$omegaJ)/($omegaI+$omegaJ)];
set betaKcurr [expr $KcurrSwitch*2.*$xDamp/($omegaI+$omegaJ)];
set betaKcomm [expr $KcommSwitch*2.*$xDamp/($omegaI+$omegaJ)];
set betaKinit [expr $KinitSwitch*2.*$xDamp/($omegaI+$omegaJ)];
rayleigh $alphaM $betaKcurr $betaKinit $betaKcomm;
set g 386
# define ground motion parameters
set patternID 1; # load pattern ID
set GMdirection 1; # ground motion direction (1 = x)
set dt 0.005000; # timestep of input GM file
set TotalNumberOfSteps 2999; # number of steps in ground motion
set GMtime [expr $dt*$TotalNumberOfSteps ]; # total time of ground motion + 10 sec of free vibration
# define the acceleration series for the ground motion
# syntax: "Series -dt $timestep_of_record -filePath $filename_with_acc_history -factor $scale_record_by_this_amount
set accelSeries "Series -dt $dt -filePath $GMfile.txt"
# create load pattern: apply acceleration to all fixed nodes with UniformExcitation
# command: pattern UniformExcitation $patternID $GMdir -accel $timeSeriesID
pattern UniformExcitation $patternID $GMdirection -accel $accelSeries;
# define dynamic analysis parameters
set dt_analysis 0.005; # timestep of analysis
wipeAnalysis; # destroy all components of the Analysis object, i.e. any objects created with system, numberer, constraints, integrator, algorithm, and analysis commands
constraints Plain; # how it handles boundary conditions
numberer RCM; # renumber dofs to minimize band-width (optimization)
system UmfPack; # how to store and solve the system of equations in the analysis
test NormDispIncr 1.0e-8 50; # type of convergence criteria with tolerance, max iterations
algorithm NewtonLineSearch; # use NewtonLineSearch solution algorithm: updates tangent stiffness at every iteration and introduces line search to the Newton-Raphson algorithm to solve the nonlinear residual equation. Line search increases the effectiveness of the Newton method
integrator Newmark 0.5 0.25; # uses Newmarks average acceleration method to compute the time history
analysis Transient; # type of analysis: transient or static
set NumSteps [expr round(($GMtime + 0.0)/$dt_analysis)]; # number of steps in analysis
# perform the dynamic analysis and display whether analysis was successful
set ok [analyze $NumSteps $dt_analysis]; # ok = 0 if analysis was completed
if {$ok == 0} {
puts "Dynamic analysis complete";
} else {
puts "Dynamic analysis did not converge";
}
# output time at end of analysis
set currentTime [getTime]; # get current analysis time (after dynamic analysis)
puts "The current time is: $currentTime";
---------------------------------------------------------------------------------------------------------------------------------------------------------------------
Regards,
Sangwook
Bilin material
Moderators: silvia, selimgunay, Moderators
Re: Bilin material
i suggest you stick with a portal frame before you launch into something larger .. start with elastic elements, check your drift, then add the zerolength elements, check drift again.