ModifiedIMK Class Reference

Class that stores funcions and material properties of a steel double symmetric I-shape profile with modified Ibarra-Medina-Krawinkler as the material model for the nonlinear springs and the OpenSeesPy command type used to model it is Bilin. More...

Inheritance diagram for ModifiedIMK:

Public Member Functions
def	__init__ (self, int ID, str Type, d, bf, tf, tw, h_1, Iy_mod, iz, E, Fy, Npl, My, L, N_G=0, K_factor=3, L_0=-1, L_b=-1, Mc=-1, K=-1, theta_u=-1, safety_factors=False)
	Constructor of the class. More...
def	Bilin (self)
	Generate the material model Bilin (Modified IMK) using the computed parameters. More...
def	CheckApplicability (self)
	Implementation of the homonym abstract method. More...
def	Computea (self)
	Method that computes the strain hardening ratio with the n modification. More...
def	Computea_s (self)
	Method that computes the modified strain hardening ratio for the spring. More...
def	ComputeK (self)
	Method that computes the residual strength ratio. More...
def	ComputeKe (self)
	Method that computes the elastic stiffness. More...
def	ComputeMc (self)
	Method that computes the capping moment. More...
def	ComputeMyStar (self)
	Method that computes the effective yield moment. More...
def	ComputeRefEnergyDissipationCap (self)
	Method that computes the reference energy dissipation capacity. More...
def	ComputeTheta_p (self)
	Method that computes the plastic rotation. More...
def	ComputeTheta_pc (self)
	Method that computes the post capping rotation. More...
def	ComputeTheta_u (self)
	Method that computes the ultimate rotation. More...
def	ComputeTheta_y (self)
	Method that computes the yield rotation. More...
def	ReInit (self, Mc=-1, K=-1, theta_u=-1)
	Implementation of the homonym abstract method. More...
def	ShowInfo (self, plot=False, block=False)
	Implementation of the homonym abstract method. More...
def	UpdateStoredData (self)
	Implementation of the homonym abstract method. More...
def	CheckApplicability (self)
	Abstract function used to check the applicability of the material model. More...

Public Attributes
	a
	a_s
	bf
	d
	data
	E
	Fy
	gamma_rm
	h_1
	ID
	Initialized
	Iy_mod
	iz
	K
	K_factor
	Ke
	L
	L_0
	L_b
	Mc
	McMy
	My
	My_star
	N_G
	Npl
	prob_factor
	rate_det
	section_name_tag
	tf
	theta_p
	theta_pc
	theta_u
	theta_y
	tw
	Type

Static Public Attributes
float	n = 10.0

Detailed Description

Class that stores funcions and material properties of a steel double symmetric I-shape profile with modified Ibarra-Medina-Krawinkler as the material model for the nonlinear springs and the OpenSeesPy command type used to model it is Bilin.

The default values are valid for a simple cantelever. For more information about the empirical model for the computation of the parameters, see Lignos Krawinkler 2011. The parameter 'n' is used as global throughout the SteelIShape sections to optimise the program (given the fact that is constant everytime).

Parameters

MaterialModels

Parent abstract class.

Definition at line 34 of file MaterialModels.py.

Constructor & Destructor Documentation

__init__()

def __init__	(		self,
		int	ID,
		str	Type,
			d,
			bf,
			tf,
			tw,
			h_1,
			Iy_mod,
			iz,
			E,
			Fy,
			Npl,
			My,
			L,
			N_G = 0,
			K_factor = 3,
			L_0 = -1,
			L_b = -1,
			Mc = -1,
			K = -1,
			theta_u = -1,
			safety_factors = False
	)

Constructor of the class.

Every argument that is optional and is initialised as -1, will be computed in this class.

Parameters

ID	(int): ID of the material model.
Type	(str): Type of the section. It can be 'Col' for column or 'Beam' for beams.
d	(float): Depth of the section.
bf	(float): Flange's width of the section
tf	(float): Flange's thickness of the section
tw	(float): Web's thickness of the section
h_1	(float): Depth excluding the flange's thicknesses and the weld fillets.
Iy_mod	(float): n modified moment of inertia (strong axis)
iz	(float): Radius of gyration (weak axis).
E	(float): Young modulus.
Fy	(float): Yield strength.
Npl	(float): Maximal vertical axial load.
My	(float): Yielding moment.
L	(float): Effective length of the element associated with this section. If the panel zone is present, exclude its dimension.
N_G	(float, optional): Gravity axial load. Defaults to 0.
K_factor	(float, optional): Rigidity factor. Defaults to 3 (assuming cantilever).
L_0	(float, optional): Position of the inflection point. Defaults to -1, e.g. computed as the total length, assuming cantilever.
L_b	(float, optional): Maximal unbraced lateral torsional buckling length. Defaults to -1, e.g. computed as the total length, assuming cantilever with no bracing support.
Mc	(float, optional): Capping moment. Defaults to -1, e.g. computed in ComputeMc.
K	(float, optional): Residual strength ratio. Defaults to -1, e.g. computed in ComputeK.
theta_u	(float, optional): Ultimate rotation. Defaults to -1, e.g. computed in ComputeTheta_u.
safety_factors	(bool, optional): Safety factors used if standard mechanical parameters are used (not test results). Defaults to False.

Exceptions

NegativeValue	ID needs to be a positive integer.
WrongArgument	Type needs to be 'Col' or 'Beam'.
NegativeValue	d needs to be positive.
NegativeValue	bf needs to be positive.
NegativeValue	tf needs to be positive.
NegativeValue	tw needs to be positive.
NegativeValue	h_1 needs to be positive.
NegativeValue	Iy_mod needs to be positive.
NegativeValue	iz needs to be positive.
NegativeValue	E needs to be positive.
NegativeValue	Fy needs to be positive.
NegativeValue	Npl needs to be positive.
NegativeValue	My needs to be positive.
NegativeValue	L needs to be positive.
NegativeValue	N_G needs to be positive.
NegativeValue	L_0 needs to be positive if different from -1.
NegativeValue	L_b needs to be positive if different from -1.
NegativeValue	Mc needs to be positive if different from -1.
NegativeValue	K needs to be positive if different from -1.
NegativeValue	theta_u needs to be positive if different from -1.
InconsistentGeometry	h_1 can't be bigger than d
MemberFailure	N_G can't be bigger than Npl (section failure).
InconsistentGeometry	L_0 can't be bigger than L

Reimplemented in ModifiedIMKSteelIShape.

Definition at line 47 of file MaterialModels.py.

        N_G = 0, K_factor = 3, L_0 = -1, L_b = -1, Mc = -1, K = -1, theta_u = -1, safety_factors = False):
        """
        Constructor of the class. Every argument that is optional and is initialised as -1, will be computed in this class.
 
        @param ID (int): ID of the material model.
        @param Type (str): Type of the section. It can be 'Col' for column or 'Beam' for beams.
        @param d (float): Depth of the section.
        @param bf (float): Flange's width of the section
        @param tf (float): Flange's thickness of the section
        @param tw (float): Web's thickness of the section
        @param h_1 (float): Depth excluding the flange's thicknesses and the weld fillets. 
        @param Iy_mod (float): n modified moment of inertia (strong axis)
        @param iz (float): Radius of gyration (weak axis).
        @param E (float): Young modulus.
        @param Fy (float): Yield strength.
        @param Npl (float): Maximal vertical axial load.
        @param My (float): Yielding moment.
        @param L (float): Effective length of the element associated with this section.
            If the panel zone is present, exclude its dimension.
        @param N_G (float, optional): Gravity axial load. Defaults to 0.
        @param K_factor (float, optional): Rigidity factor. Defaults to 3 (assuming cantilever).
        @param L_0 (float, optional): Position of the inflection point.
            Defaults to -1, e.g. computed as the total length, assuming cantilever.
        @param L_b (float, optional): Maximal unbraced lateral torsional buckling length.
            Defaults to -1, e.g. computed as the total length, assuming cantilever with no bracing support.
        @param Mc (float, optional): Capping moment. Defaults to -1, e.g. computed in ComputeMc.
        @param K (float, optional): Residual strength ratio. Defaults to -1, e.g. computed in ComputeK.
        @param theta_u (float, optional): Ultimate rotation. Defaults to -1, e.g. computed in ComputeTheta_u.
        @param safety_factors (bool, optional): Safety factors used if standard mechanical parameters are used (not test results). Defaults to False.
 
        @exception NegativeValue: ID needs to be a positive integer.
        @exception WrongArgument: Type needs to be 'Col' or 'Beam'.
        @exception NegativeValue: d needs to be positive.
        @exception NegativeValue: bf needs to be positive.
        @exception NegativeValue: tf needs to be positive.
        @exception NegativeValue: tw needs to be positive.
        @exception NegativeValue: h_1 needs to be positive.
        @exception NegativeValue: Iy_mod needs to be positive.
        @exception NegativeValue: iz needs to be positive.
        @exception NegativeValue: E needs to be positive.
        @exception NegativeValue: Fy needs to be positive.
        @exception NegativeValue: Npl needs to be positive.
        @exception NegativeValue: My needs to be positive.
        @exception NegativeValue: L needs to be positive.
        @exception NegativeValue: N_G needs to be positive.
        @exception NegativeValue: L_0 needs to be positive if different from -1.
        @exception NegativeValue: L_b needs to be positive if different from -1.
        @exception NegativeValue: Mc needs to be positive if different from -1.
        @exception NegativeValue: K needs to be positive if different from -1.
        @exception NegativeValue: theta_u needs to be positive if different from -1.
        @exception InconsistentGeometry: h_1 can't be bigger than d
        @exception MemberFailure: N_G can't be bigger than Npl (section failure).
        @exception InconsistentGeometry: L_0 can't be bigger than L
        """
        # Check
        if ID < 0: raise NegativeValue()
        if Type != "Beam" and Type != "Col": raise WrongArgument()
        if d < 0: raise NegativeValue()
        if bf < 0: raise NegativeValue()
        if tf < 0: raise NegativeValue()
        if tw < 0: raise NegativeValue()
        if h_1 < 0: raise NegativeValue()
        if Iy_mod < 0: raise NegativeValue()
        if iz < 0: raise NegativeValue()
        if E < 0: raise NegativeValue()
        if Fy < 0: raise NegativeValue()
        if Npl < 0: raise NegativeValue()
        if My < 0: raise NegativeValue()
        if L < 0: raise NegativeValue()
        if N_G < 0: raise NegativeValue()
        if L_0 != -1 and L_0 < 0: raise NegativeValue()
        if L_b != -1 and L_b < 0: raise NegativeValue()
        if Mc != -1 and Mc < 0: raise NegativeValue()
        if K != -1 and K < 0: raise NegativeValue()
        if theta_u != -1 and theta_u < 0: raise NegativeValue()
        if h_1 > d: raise InconsistentGeometry()
        if N_G > Npl: raise MemberFailure()
        if L_0 > L: raise InconsistentGeometry()
 
        # Arguments
        self.Type = Type
        self.ID = ID
        self.d = d
        self.bf = bf
        self.tf = tf
        self.tw = tw
        self.h_1 = h_1
        self.Iy_mod = Iy_mod
        self.iz = iz
        self.E = E
        self.Fy = Fy
        self.Npl = Npl
        self.My = My
        self.L = L
        self.N_G = N_G
        self.K_factor = K_factor
        self.L_0 = L if L_0 == -1 else L_0
        self.L_b = L if L_b == -1 else L_b
 
        # Initialized the parameters that are dependent from others
        self.section_name_tag = "None"
        self.Initialized = False
        if safety_factors:
            self.gamma_rm = 1.25
            self.prob_factor = 1.15
        else:
            self.gamma_rm = 1.0
            self.prob_factor = 1.0
        self.ReInit(Mc, K, theta_u)
 
 

Member Function Documentation

Bilin()

def Bilin

(

self

)

Generate the material model Bilin (Modified IMK) using the computed parameters.

See _Bilin function for more information.

Definition at line 478 of file MaterialModels.py.

    def Bilin(self):
        """
        Generate the material model Bilin (Modified IMK) using the computed parameters.
        See _Bilin function for more information.
        """
        _Bilin(self.ID, self.Ke, self.a_s, self.My_star, self.theta_p, self.theta_pc, self.K, self.theta_u, self.rate_det)
        self.Initialized = True
        self.UpdateStoredData()
 
 

CheckApplicability()

def CheckApplicability

(

self

)

Implementation of the homonym abstract method.

See parent class MaterialModels for detailed information.

Reimplemented from MaterialModels.

Definition at line 289 of file MaterialModels.py.

    def CheckApplicability(self):
        """
        Implementation of the homonym abstract method.
        See parent class MaterialModels for detailed information.
        """
        Check = True
        if self.Type == "Beam":
            if self.d/self.tw < 20 or self.d/self.tw > 55:
                Check = False
                print("The d/tw check was not fullfilled")
            if self.L_b/self.iz < 20 or self.L_b/self.iz > 80:
                Check = False
                print("The Lb/iz check was not fullfilled")
            if self.bf/2/self.tf < 4 or self.bf/2/self.tf  > 8:
                Check = False
                print("The bf/2/tf check was not fullfilled")
            if self.L/self.d < 2.5 or self.L/self.d > 7:
                Check = False
                print("The check L/d was not fullfilled")
            if self.d < 102*mm_unit or self.d > 914*mm_unit:
                Check = False
                print("The d check was not fullfilled")
            if self.Fy < 240*MPa_unit or self.Fy > 450*MPa_unit:
                Check = False
                print("The Fy check was not fullfilled")
        else:
            if self.h_1/self.tw < 3.71 or self.d/self.tw > 57.5:
                Check = False
                print("The h1/tw check was not fullfilled")
            if self.L_b/self.iz < 38.4 or self.L_b/self.iz > 120:
                Check = False
                print("The Lb/iz check was not fullfilled")
            if self.N_G/self.Npl < 0 or self.N_G/self.Npl > 0.75:
                Check = False
                print("The NG/Npl check was not fullfilled")
        if not Check:
            print("The validity of the equations is not fullfilled.")
            print("!!!!!!! WARNING !!!!!!! Check material model of Modified IMK, ID=", self.ID)
            print("")
 

Computea()

def Computea

(

self

)

Method that computes the strain hardening ratio with the n modification.

Returns

float: Strain hardening ratio.

Definition at line 337 of file MaterialModels.py.

    def Computea(self):
        """
        Method that computes the strain hardening ratio with the n modification.
 
        @returns float: Strain hardening ratio.
        """
        # strain hardening ratio of spring
        return (n+1.0)*self.My_star*(self.McMy-1.0)/(self.Ke*self.theta_p)
 

Computea_s()

def Computea_s

(

self

)

Method that computes the modified strain hardening ratio for the spring.

For more info see Ibarra & Krawinkler 2005.

Returns

float: Strain hardening ratio.

Definition at line 346 of file MaterialModels.py.

    def Computea_s(self):
        """
        Method that computes the modified strain hardening ratio for the spring.
        For more info see Ibarra & Krawinkler 2005.
 
        @returns float: Strain hardening ratio.
        """
        return self.a/(1.0+n*(1.0-self.a))
 

ComputeK()

def ComputeK

(

self

)

Method that computes the residual strength ratio.

For more info see Lignos & Krawinkler 2011 and Lignos et Al. 2019.

Returns

float: Residual strength ratio.

Definition at line 384 of file MaterialModels.py.

    def ComputeK(self):
        """
        Method that computes the residual strength ratio.
        For more info see Lignos & Krawinkler 2011 and Lignos et Al. 2019.
 
        @returns float: Residual strength ratio.
        """
        if self.Type == "Beam":
            return 0.4
        else:
            tmp = 0.5-0.4*self.N_G/self.Npl
            return max(tmp, 0)
    

ComputeKe()

def ComputeKe

(

self

)

Method that computes the elastic stiffness.

Returns

float: The stiffness

Definition at line 329 of file MaterialModels.py.

    def ComputeKe(self):
        """
        Method that computes the elastic stiffness.
 
        @returns float: The stiffness
        """
        return self.K_factor*n*self.E*self.Iy_mod/self.L
 

ComputeMc()

def ComputeMc

(

self

)

Method that computes the capping moment.

For more info see Lignos & Krawinkler 2011 and Lignos et Al. 2019.

Returns

float: Capping moment.

Definition at line 370 of file MaterialModels.py.

    def ComputeMc(self):
        """
        Method that computes the capping moment.
        For more info see Lignos & Krawinkler 2011 and Lignos et Al. 2019.
 
        @returns float: Capping moment.
        """
        if self.Type == "Beam":
            return self.My_star*1.11
            # For RBS: My_star*1.09
        else:
            tmp = 12.5*(self.h_1/self.tw)**(-0.2)*(self.L_b/self.iz)**(-0.4)*(1-self.N_G/self.Npl)**0.4
            return max(min(1.3, tmp), 1.0)*self.My_star
 

ComputeMyStar()

def ComputeMyStar

(

self

)

Method that computes the effective yield moment.

For more info see Lignos & Krawinkler 2011 and Lignos et Al. 2019.

Returns

float: Effective yield moment.

Definition at line 355 of file MaterialModels.py.

    def ComputeMyStar(self):
        """
        Method that computes the effective yield moment.
        For more info see Lignos & Krawinkler 2011 and Lignos et Al. 2019.
 
        @returns float: Effective yield moment.
        """
        if self.Type == "Beam":
            return self.prob_factor*self.My*self.gamma_rm*1.1
        else:
            if self.N_G/self.Npl > 0.2:
                return 1.15*self.prob_factor*self.My*self.gamma_rm*(1-self.N_G/self.Npl)*9.0/8.0
            else:
                return 1.15*self.prob_factor*self.My*self.gamma_rm*(1-self.N_G/2.0/self.Npl)
 

ComputeRefEnergyDissipationCap()

def ComputeRefEnergyDissipationCap

(

self

)

Method that computes the reference energy dissipation capacity.

For more info see Lignos & Krawinkler 2011 and Lignos et Al. 2019.

Returns

float: Reference energy dissipation capacity.

Definition at line 456 of file MaterialModels.py.

    def ComputeRefEnergyDissipationCap(self):
        """
        Method that computes the reference energy dissipation capacity.
        For more info see Lignos & Krawinkler 2011 and Lignos et Al. 2019.
 
        @returns float: Reference energy dissipation capacity.
        """
        if self.Type == "Beam":
            if self.d < 533.0*mm_unit:
                return 495.0*(self.h_1/self.tw)**(-1.34)*(self.bf/2.0/self.tf)**(-0.595)*(self.Fy/(355.0*MPa_unit))**(-0.360)
            else:
                return 536.0*(self.h_1/self.tw)**(-1.26)*(self.bf/2.0/self.tf)**(-0.525)*(self.L_b/self.iz)**(-0.130)*(self.Fy/(355.0*MPa_unit))**(-0.291)
                # With RBS: ...
        else:
            if self.N_G/self.Npl > 0.35:
                tmp = 268000.0*(self.h_1/self.tw)**(-2.30)*(self.L_b/self.iz)**(-1.130)*(1.0-self.N_G/self.Npl)**(1.19)
                return min(tmp, 3.0)
            else:
                tmp = 25000.0*(self.h_1/self.tw)**(-2.14)*(self.L_b/self.iz)**(-0.53)*(1.0-self.N_G/self.Npl)**(4.92)
                return min(tmp, 3.0)
 
 

ComputeTheta_p()

def ComputeTheta_p

(

self

)

Method that computes the plastic rotation.

For more info see Lignos & Krawinkler 2011 and Lignos et Al. 2019.

Returns

float: Plastic rotation.

Definition at line 406 of file MaterialModels.py.

    def ComputeTheta_p(self):
        """
        Method that computes the plastic rotation.
        For more info see Lignos & Krawinkler 2011 and Lignos et Al. 2019.
 
        @returns float: Plastic rotation.
        """
        if self.Type == "Beam":
            if self.d < 533.0*mm_unit:
                return 0.0865*(self.h_1/self.tw)**(-0.365)*(self.bf/2.0/self.tf)**(-0.14)*(self.L_0/self.d)**(0.34)*(self.d/(533.0*mm_unit))**(-0.721)*(self.Fy/(355.0*MPa_unit))**(-0.23)
            else:
                return 0.318*(self.h_1/self.tw)**(-0.550)*(self.bf/2.0/self.tf)**(-0.345)*(self.L_0/self.d)**(0.090)*(self.L_b/self.iz)**(-0.023)*(self.d/(533.0*mm_unit))**(-0.330)*(self.Fy/(355.0*MPa_unit))**(-0.130)
                # With RBS: ...
        else:
            tmp = 294.0*(self.h_1/self.tw)**(-1.7)*(self.L_b/self.iz)**(-0.7)*(1.0-self.N_G/self.Npl)**(1.6) # *(self.E/self.Fy/gamma_rm)**(0.2) # EC8
            if tmp > 0.2:
                tmp = 0.2
            # if tmp > self.theta_u-self.theta_y:
            #     tmp = (self.theta_u-self.theta_y)*0.799 # convergence issue
            return tmp
 

ComputeTheta_pc()

def ComputeTheta_pc

(

self

)

Method that computes the post capping rotation.

For more info see Lignos & Krawinkler 2011 and Lignos et Al. 2019.

Returns

float: Post capping rotation.

Definition at line 427 of file MaterialModels.py.

    def ComputeTheta_pc(self):
        """
        Method that computes the post capping rotation.
        For more info see Lignos & Krawinkler 2011 and Lignos et Al. 2019.
 
        @returns float: Post capping rotation.
        """
        if self.Type == "Beam":
            if self.d < 533.0*mm_unit:
                return 5.63*(self.h_1/self.tw)**(-0.565)*(self.bf/2.0/self.tf)**(-0.800)*(self.d/(533.0*mm_unit))**(-0.280)*(self.Fy/(355.0*MPa_unit))**(-0.430)
            else:
                return 7.50*(self.h_1/self.tw)**(-0.610)*(self.bf/2.0/self.tf)**(-0.710)*(self.L_b/self.iz)**(-0.110)*(self.d/(533.0*mm_unit))**(-0.161)*(self.Fy/(355.0*MPa_unit))**(-0.320)
                # With RBS: ...
        else:
            tmp =  90.0*(self.h_1/self.tw)**(-0.8)*(self.L_b/self.iz)**(-0.8)*(1.0-self.N_G/self.Npl)**(2.5) # *(self.E/self.Fy/gamma_rm)**(0.07) # EC8
            return min(tmp, 0.3)
 

ComputeTheta_u()

def ComputeTheta_u

(

self

)

Method that computes the ultimate rotation.

For more info see Lignos & Krawinkler 2011 and Lignos et Al. 2019.

Returns

float: Ultimate rotation.

Definition at line 444 of file MaterialModels.py.

    def ComputeTheta_u(self):
        """
        Method that computes the ultimate rotation.
        For more info see Lignos & Krawinkler 2011 and Lignos et Al. 2019.
 
        @returns float: Ultimate rotation.
        """
        if self.Type == "Beam":
            return 0.2
        else:
            return 0.15
 

ComputeTheta_y()

def ComputeTheta_y

(

self

)

Method that computes the yield rotation.

For more info see Lignos & Krawinkler 2011 and Lignos et Al. 2019.

Returns

float: Yield rotation.

Definition at line 397 of file MaterialModels.py.

    def ComputeTheta_y(self):
        """
        Method that computes the yield rotation.
        For more info see Lignos & Krawinkler 2011 and Lignos et Al. 2019.
 
        @returns float: Yield rotation.
        """
        return self.My_star/self.Ke*(n+1)
 

ReInit()

def ReInit	(		self,
			Mc = -1,
			K = -1,
			theta_u = -1
	)

Implementation of the homonym abstract method.

See parent class DataManagement for detailed information.

Parameters

Mc	(float, optional): Capping moment. Defaults to -1, e.g. computed in ComputeMc.
K	(float, optional): Residual strength ratio. Defaults to -1, e.g. computed in ComputeK.
theta_u	(float, optional): Ultimate rotation. Defaults to -1, e.g. computed in ComputeTheta_u.

Definition at line 160 of file MaterialModels.py.

    def ReInit(self, Mc = -1, K = -1, theta_u = -1):
        """
        Implementation of the homonym abstract method.
        See parent class DataManagement for detailed information.
 
        @param Mc (float, optional): Capping moment. Defaults to -1, e.g. computed in ComputeMc.
        @param K (float, optional): Residual strength ratio. Defaults to -1, e.g. computed in ComputeK.
        @param theta_u (float, optional): Ultimate rotation. Defaults to -1, e.g. computed in ComputeTheta_u.
        """
        # Precompute some members
        self.My_star = self.ComputeMyStar()
 
        # Arguments
        self.Mc = self.ComputeMc() if Mc == -1 else Mc
        self.K = self.ComputeK() if K == -1 else K
        self.theta_u = self.ComputeTheta_u() if theta_u == -1 else theta_u
 
        # Check applicability
        self.CheckApplicability()
 
        # Members
        self.Ke = self.ComputeKe()
        self.theta_y = self.ComputeTheta_y()
        self.theta_p = self.ComputeTheta_p()
        self.theta_pc = self.ComputeTheta_pc()
        self.McMy = self.Mc/self.My_star
        self.rate_det = self.ComputeRefEnergyDissipationCap()
        self.a = self.Computea()
        self.a_s = self.Computea_s()
        if self.section_name_tag != "None": self.section_name_tag = self.section_name_tag + " (modified)"
 
        # Data storage for loading/saving
        self.UpdateStoredData()
        
 

ShowInfo()

def ShowInfo	(		self,
			plot = False,
			block = False
	)

Implementation of the homonym abstract method.

See parent class DataManagement for detailed information.

Parameters

plot	(bool, optional): Option to show the plot of the material model. Defaults to False.
block	(bool, optional): Option to wait the user command 'plt.show()' (avoiding the stop of the program everytime that a plot should pop up). Defaults to False.

Definition at line 237 of file MaterialModels.py.

    def ShowInfo(self, plot = False, block = False):
        """
        Implementation of the homonym abstract method.
        See parent class DataManagement for detailed information.
        
        @param plot (bool, optional): Option to show the plot of the material model. Defaults to False.
        @param block (bool, optional): Option to wait the user command 'plt.show()' (avoiding the stop of the program everytime that a plot should pop up). Defaults to False.
        """
        Mr = self.K*self.My_star
        theta_p_plot = self.theta_p
        if self.theta_p > self.theta_u-self.theta_y:
            theta_p_plot = self.theta_u-self.theta_y
        theta_r = self.theta_y + theta_p_plot + self.theta_pc*(1.0-Mr/self.Mc)
        if theta_r > self.theta_u:
            theta_r = self.theta_u
            Mr = self.Mc*(1.0-1.0/self.theta_pc*(self.theta_u-self.theta_y-theta_p_plot))
 
        print("")
        print("Requested info for Modified IMK (Ibarra-Medina-Krawinkler) material model Parameters, ID = {}".format(self.ID))
        print("Section associated: {}".format(self.section_name_tag))
        print('theta y = {}'.format(self.theta_y))
        print('theta p = {}'.format(self.theta_p))
        print('theta r = {}'.format(theta_r))
        print('theta pc = {}'.format(self.theta_pc))
        print('theta u = {}'.format(self.theta_u))
        print('My star = {} kNm'.format(self.My_star/kNm_unit))
        print('Mc = {} kNm'.format(self.Mc/kNm_unit))
        print('Mr = {} kNm'.format(Mr/kNm_unit))
        print('a = {} '.format(self.a))
        print('as = {} '.format(self.a_s))
        print('lambda (deterioration rate) = {} '.format(self.rate_det))
        print("")
        
        if plot:
            # Data for plotting
            x_axis = np.array([0.0, self.theta_y, self.theta_y + theta_p_plot, theta_r, self.theta_u, self.theta_u])
            x_axis2 = np.array([self.theta_y + theta_p_plot, self.theta_y + theta_p_plot + self.theta_pc])
            y_axis = np.array([0.0, self.My_star, self.Mc, Mr, Mr, 0.0])/kNm_unit
            y_axis2 = np.array([self.Mc, 0.0])/kNm_unit
 
            fig, ax = plt.subplots()
            ax.plot(x_axis, y_axis, 'k-')
            ax.plot(x_axis2, y_axis2, 'k--')
 
            ax.set(xlabel='Rotation [rad]', ylabel='Moment [kNm]', 
                title='Modified IMK deterioration model (ID={})'.format(self.ID))
            ax.grid()
 
            if block:
                plt.show()
 
 

UpdateStoredData()

def UpdateStoredData

(

self

)

Implementation of the homonym abstract method.

See parent class DataManagement for detailed information.

Definition at line 195 of file MaterialModels.py.

    def UpdateStoredData(self):
        """
        Implementation of the homonym abstract method.
        See parent class DataManagement for detailed information.
        """
        self.data = [["INFO_TYPE", "ModifiedIMK"], # Tag for differentiating different data
            ["ID", self.ID],
            ["section_name_tag", self.section_name_tag], 
            ["Type", self.Type],
            ["d", self.d],
            ["bf", self.bf],
            ["tf", self.tf],
            ["tw", self.tw],
            ["h_1", self.h_1],
            ["Iy_mod", self.Iy_mod],
            ["iz", self.iz],
            ["E", self.E],
            ["Fy", self.Fy],
            ["L", self.L],
            ["N_G", self.N_G],
            ["K_factor", self.K_factor],
            ["Ke", self.Ke],
            ["L_0", self.L_0],
            ["L_b", self.L_b],
            ["gamma_rm", self.gamma_rm],
            ["prob_factor", self.prob_factor],
            ["Npl", self.Npl],
            ["My", self.My],
            ["My_star", self.My_star],
            ["Mc", self.Mc],
            ["McMy", self.McMy],
            ["K", self.K],
            ["theta_y", self.theta_y],
            ["theta_p", self.theta_p],
            ["theta_pc", self.theta_pc],
            ["theta_u", self.theta_u],
            ["rate_det", self.rate_det],
            ["a", self.a],
            ["a_s", self.a_s],
            ["Initialized", self.Initialized]]
 
 

Member Data Documentation

a

a

Definition at line 187 of file MaterialModels.py.

a_s

a_s

Definition at line 188 of file MaterialModels.py.

bf

bf

Definition at line 131 of file MaterialModels.py.

d

d

Definition at line 130 of file MaterialModels.py.

data

data

Definition at line 200 of file MaterialModels.py.

E

E

Definition at line 137 of file MaterialModels.py.

Fy

Fy

Definition at line 138 of file MaterialModels.py.

gamma_rm

gamma_rm

Definition at line 151 of file MaterialModels.py.

h_1

h_1

Definition at line 134 of file MaterialModels.py.

ID

ID

Definition at line 129 of file MaterialModels.py.

Initialized

Initialized

Definition at line 149 of file MaterialModels.py.

Iy_mod

Iy_mod

Definition at line 135 of file MaterialModels.py.

iz

iz

Definition at line 136 of file MaterialModels.py.

K

K

Definition at line 174 of file MaterialModels.py.

K_factor

K_factor

Definition at line 143 of file MaterialModels.py.

Ke

Ke

Definition at line 181 of file MaterialModels.py.

L

L

Definition at line 141 of file MaterialModels.py.

L_0

L_0

Definition at line 144 of file MaterialModels.py.

L_b

L_b

Definition at line 145 of file MaterialModels.py.

Mc

Mc

Definition at line 173 of file MaterialModels.py.

McMy

McMy

Definition at line 185 of file MaterialModels.py.

My

My

Definition at line 140 of file MaterialModels.py.

My_star

My_star

Definition at line 170 of file MaterialModels.py.

n

float n = 10.0

static

Definition at line 45 of file MaterialModels.py.

N_G

N_G

Definition at line 142 of file MaterialModels.py.

Npl

Npl

Definition at line 139 of file MaterialModels.py.

prob_factor

prob_factor

Definition at line 152 of file MaterialModels.py.

rate_det

rate_det

Definition at line 186 of file MaterialModels.py.

section_name_tag

section_name_tag

Definition at line 148 of file MaterialModels.py.

tf

tf

Definition at line 132 of file MaterialModels.py.

theta_p

theta_p

Definition at line 183 of file MaterialModels.py.

theta_pc

theta_pc

Definition at line 184 of file MaterialModels.py.

theta_u

theta_u

Definition at line 175 of file MaterialModels.py.

theta_y

theta_y

Definition at line 182 of file MaterialModels.py.

tw

tw

Definition at line 133 of file MaterialModels.py.

Type

Type

Definition at line 128 of file MaterialModels.py.

---

The documentation for this class was generated from the following file:

• /media/carmine/DATA/Programmi/OpenSeesPyAssistant/MaterialModels.py