Modified Ibarra-Medina-Krawinkler Deterioration Model with Bilinear Hysteretic Response (Bilin Material)

• Command_Manual
  • Tcl Commands
  • Modeling_Commands
    • model
    • uniaxialMaterial
    • ndMaterial
    • frictionModel
    • section
    • geometricTransf
    • element
    • node
    • sp commands
      • fix
      • fixX
      • fixY
      • fixZ
    • mp commands
      • equalDOF
      • rigidDiaphragm
      • rigidLink
    • timeSeries
    • pattern
    • mass
    • block commands
      • block2D
      • block3D
    • region
    • rayleigh
  • Analysis Commands
    • constraints Command
    • numberer Command
    • system Command
    • test Command
    • algorithm Command
    • integrator Command
    • analysis Command
    • eigen Command
    • analyze Command
  • Output Commands
    • recorder
    • record
    • print
    • printA
    • logFile
    • RealTime_Output_Commands
  • Misc Commands
    • eleResponse Command
    • getTime Command
    • getEleTags Command
    • getNodeTags
    • loadConst Command
    • nodeCoord Command
    • nodeDisp Command
    • nodeVel Command
    • nodeAccel Command
    • nodeEigenvector Command
    • nodeBounds
    • print Command
    • remove Command
    • reset Command
    • testIter
    • testNorms
    • wipe Command
    • exit Command
  • DataBase Commands
    • Template:DataBaseCommandsMenu

---

The original Bilin material has been revised and should no longer be used. As the input format for the new material is slightly different, the old material will remain in the executable so that we do not break old scripts. However we will no longer document the command to use the old code. The new and improved material documentation is here:

https://opensees.github.io/OpenSeesDocumentation/user/manual/material/uniaxialMaterials/IMKBilin.html

---

The original command and documentation are included below for old users who still wish to use the material.

This command is used to construct a bilin material. The bilin material simulates the modified Ibarra-Krawinkler deterioration model with bilinear hysteretic response. Note that the hysteretic response of this material has been calibrated with respect to more than 350 experimental data of steel beam-to-column connections and multivariate regression formulas are provided to estimate the deterioration parameters of the model for different connection types. These relationships were developed by Lignos and Krawinkler (2009, 2011) and have been adopted by PEER/ATC (2010). The input parameters for this component model can be computed interactively from this [link: http://resslabtools.epfl.ch/]. Use the module Component Model.

uniaxialMaterial 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 <$nFactor>

---

$matTag	integer tag identifying material
$K0	elastic stiffness
$as_Plus	strain hardening ratio for positive loading direction
$as_Neg	strain hardening ratio for negative loading direction
$My_Plus	effective yield strength for positive loading direction
$My_Neg	effective yield strength for negative loading direction (negative value)
$Lamda_S	Cyclic deterioration parameter for strength deterioration [E_t=Lamda_S*M_y; set Lamda_S = 0 to disable this mode of deterioration]
$Lamda_C	Cyclic deterioration parameter for post-capping strength deterioration [E_t=Lamda_C*M_y; set Lamda_C = 0 to disable this mode of deterioration]
$Lamda_A	Cyclic deterioration parameter for acceleration reloading stiffness deterioration (is not a deterioration mode for a component with Bilinear hysteretic response) [Input value is required, but not used; set Lamda_A = 0].
$Lamda_K	Cyclic deterioration parameter for unloading stiffness deterioration [E_t=Lamda_K*M_y; set Lamda_k = 0 to disable this mode of deterioration]
$c_S	rate of strength deterioration. The default value is 1.0.
$c_C	rate of post-capping strength deterioration. The default value is 1.0.
$c_A	rate of accelerated reloading deterioration. The default value is 1.0.
$c_K	rate of unloading stiffness deterioration. The default value is 1.0.
$theta_p_Plus	pre-capping rotation for positive loading direction (often noted as plastic rotation capacity)
$theta_p_Neg	pre-capping rotation for negative loading direction (often noted as plastic rotation capacity) (positive value)
$theta_pc_Plus	post-capping rotation for positive loading direction
$theta_pc_Neg	post-capping rotation for negative loading direction (positive value)
$Res_Pos	residual strength ratio for positive loading direction
$Res_Neg	residual strength ratio for negative loading direction (positive value)
$theta_u_Plus	ultimate rotation capacity for positive loading direction
$theta_u_Neg	ultimate rotation capacity for negative loading direction (positive value)
$D_Plus	rate of cyclic deterioration in the positive loading direction (this parameter is used to create assymetric hysteretic behavior for the case of a composite beam). For symmetric hysteretic response use 1.0.
$D_Neg	rate of cyclic deterioration in the negative loading direction (this parameter is used to create assymetric hysteretic behavior for the case of a composite beam). For symmetric hysteretic response use 1.0.
$nFactor	elastic stiffness amplification factor, mainly for use with concentrated plastic hinge elements (optional, default = 0).

Image from: Lignos et al. (2011)

The deterioration model parameters can be calibrated based on actual experimental data of steel components. Instead, a user can determine the same parameters with predictive equations that have been developed by Lignos and Krawinkler (2009, 2011). These equations relate deterioration modelling parameters with geometric and material properties of steel components and have been developed with multivariate regression analysis after extensive calibration of a recently developed steel component database of more than 300 steel specimens (Lignos and Krawinkler 2007, 2009, 2011).

---

References:

[1]	Ibarra L.F., and Krawinkler, H. (2005). “Global collapse of frame structures under seismic excitations”, Rep. No. TB 152, The John A. Blume Earthquake Engineering Center, Stanford University, Stanford, CA.
[2]	Ibarra L.F., Medina R. A., and Krawinkler H. (2005). “Hysteretic models that incorporate strength and stiffness deterioration”, Earthquake Engineering and Structural Dynamics, 34(12), 1489-1511.
[3]	Karavasilis T.L., Ricles J.M., Sause R. (2009). "Implementation of deterioration elements in OpenSEES for collapse simulations" ATLSS Engineering Research Center, Rep. No. 09-11.
[4]	Lignos, D.G., and Krawinkler, H. (2007), “A database in support of modeling of component deterioration for collapse prediction of steel frame structures”, Proc. ASCE Structures Congress, Long Beach CA, SEI institute.
[5]	Lignos, D.G. and Krawinkler, H. (2012). “Sidesway collapse of deteriorating structural systems under seismic excitations,” Rep.No.TB 177, The John A. Blume Earthquake Engineering Research Center, Stanford University, Stanford, CA. [electronic version: https://blume.stanford.edu/tech_reports]
[6]	Lignos, D. G., Krawinkler, H. and Whittaker, A.S. (2011). "Prediction and validation of sidesway collapse of two scale models of a 4-story steel moment frame", Earthquake Engineering and Structural Dynamics, EESD, Vol. 40 (7), 807-825.
[7]	Lignos, D.G., and Krawinkler, H. (2010). “A steel database for component deterioration of tubular hollow square steel columns under varying axial load for collapse assessment of steel structures under earthquakes”, Proc. 7th International Conference on Urban Earthquake Engineering (7CUEE), March 3-5, Tokyo, Japan.
[8]	Lignos, D.G., and Krawinkler, H. (2011). “Deterioration modeling of steel components in support of collapse prediction of steel moment frames under earthquake loading”, Journal of Structural Engineering, ASCE, Vol. 137 (11), 1291-1302. File:BilinMaterialOpenSees.pdf
[9]	PEER/ATC (2010). “Modeling and acceptance criteria for seismic design and analysis of tall buildings”, PEER/ATC-72-1, prepared by the Applied Technology Council in cooperation with the Pacific Earthquake Engineering Research Center, Redwood City, CA..
[10]	Lignos, D.G., Hartloper, A.R., Elkady, A.M.A., Deirlein, G.G., Hamburger, R. (2019). “Proposed Updates to the ASCE 41 Nonlinear Modeling Parameters for Wide-Flange Steel Columns in Support of Performance-Based Seismic Engineering”, Journal of Structural Engineering, 04019083, https://doi.org/10.1061/(ASCE)ST.1943-541X.0002353.

Code Developed by : by Dr. Dimitrios G. Lignos, McGill University and was implemented by Dr. Theodore Karavasilis, Oxford University

nFactor Implemented by: by Dr. Andre R. Barbosa, Oregon State University and Mr. Filipe L.A. Ribeiro, Universidade Nova de Lisboa