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
- mp commands
- timeSeries
- pattern
- mass
- block commands
- region
- rayleigh
- Analysis Commands
- Output Commands
- Misc Commands
- DataBase Commands
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:
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