Analysis Class Reference

Class dedicated to the analysis of the OpenSeesPy model. More...

Public Member Functions
def	__init__ (self, str data_dir, str name_ODB, algo="KrylovNewton", test_type="NormDispIncr", test_opt=0, max_iter=MAX_ITER, tol=TOL, allow_smaller_step=False)
	The constructor of the class. More...
def	DeformedShape (self, scale=1, animate=False, dt=0.01)
	Method that shows the final deformed shape of the model. More...
def	FiberResponse (self, int ele_fiber_ID_analysed, fiber_section=1, animate_stress=False, animate_strain=False, fps=25)
	Method that shows the final stress response of the fiber section chosen. More...
def	Gravity (self, list loaded_nodes, list Fy, int timeSeries_ID, int pattern_ID, n_step=10, timeSeries_type="Linear", pattern_type="Plain", constraints_type="Plain", numberer_type="RCM", system_type="BandGeneral", analysis_type="Static", show_plot=False)
	Method to perform the gravity analyisis with vertical loadings (load-control). More...
def	LateralForce (self, list loaded_nodes, list Fx, int timeSeries_ID, int pattern_ID, n_step=1000, fiber_ID_analysed=-1, fiber_section=1, timeSeries_type="Linear", pattern_type="Plain", constraints_type="Plain", numberer_type="RCM", system_type="BandGeneral", analysis_type="Static", show_plot=True, block=False)
	Method to perform the lateral force analyisis with lateral loading (load-control). More...
def	LoadingProtocol (self, int CtrlNode, np.ndarray discr_LP, int timeSeries_ID, int pattern_ID, Fx=1 *kN_unit, ele_fiber_ID_analysed=-1, fiber_section=1, timeSeries_type="Linear", pattern_type="Plain", constraints_type="Plain", numberer_type="RCM", system_type="UmfPack", analysis_type="Static", show_plot=True, block=False)
	Method to perform a loading protocol analysis (displacement-control). More...
def	Pushover (self, int CtrlNode, Dmax, Dincr, int timeSeries_ID, int pattern_ID, Fx=1 *kN_unit, ele_fiber_ID_analysed=-1, fiber_section=1, timeSeries_type="Linear", pattern_type="Plain", constraints_type="Plain", numberer_type="RCM", system_type="UmfPack", analysis_type="Static", show_plot=True, block=False)
	Method to perform a pushover analysis (displacement-control). More...

Public Attributes
	algo
	allow_smaller_step
	data_dir
	load_case
	max_iter
	name_ODB
	test_opt
	test_type
	tol

Detailed Description

Class dedicated to the analysis of the OpenSeesPy model.

The Gravity method should be run first to perform the Load-control analysis (apply the vertical load). If no vertical load, this method can be omitted.
Then only one of the Displacement-control (Pushover or LoadingProtocol) or Load-control (LateralForce) analysis can ran.
After the analysis reach convergence in the last step, for the postprocessing, the DeformedShape method can be used to see the final deformed shape and the animation of the entire loading protocol; the FiberResponse method can be used to see the animation of the same fiber section recorded during the analysis (strain and/or stress).

Definition at line 16 of file AnalysisAndPostProcessing.py.

Constructor & Destructor Documentation

__init__()

def __init__	(		self,
		str	data_dir,
		str	name_ODB,
			algo = "KrylovNewton",
			test_type = "NormDispIncr",
			test_opt = 0,
			max_iter = MAX_ITER,
			tol = TOL,
			allow_smaller_step = False
	)

The constructor of the class.

Parameters

data_dir	(str): Directory in which the results from the analysis will be stored. Use the recorders (from OpenSeesPy) or the Record method from MemberModel.
name_ODB	(str): Name for the folder in which the data for the animations and the fibers are stored.
algo	(str, optional): Type of alghoritm chosen for the analysis. It detemines how to construct a SolutionAlgorithm object, which determines the sequence of steps taken to solve the non-linear equation. For more information on the available types, see the OpenSeesPy documentation. Defaults to "KrylovNewton".
test_type	(str, optional): Type of test chosen for the analysis. It determines how to construct a ConvergenceTest object. Certain SolutionAlgorithm objects require a ConvergenceTest object to determine if convergence has been achieved at the end of an iteration step. For more information on the available types, see the OpenSeesPy documentation. Defaults to "NormDispIncr".
test_opt	(int, optional): Print-flag from 0 to 5 used to receive more info during the iteration (for example: 0 print nothing and 2 print information on norms and number of iterations at end of successful test). For more information, see the OpenSeesPy documentation. Defaults to 0.
max_iter	(float, optional): Maximal number of iterations to check. Defaults to MAX_ITER (from Constants Module).
tol	(float, optional): Tolerance criteria used to check for convergence. Defaults to TOL (from Constants Module).
allow_smaller_step	(bool, optional): Allow smaller steps in the displacement-control analysis. Defaults to False.

Exceptions

NegativeValue	The argument max_iter should be positive.
NegativeValue	The argument tol should be positive.

Definition at line 22 of file AnalysisAndPostProcessing.py.

    def __init__(self, data_dir: str,  name_ODB: str, algo = "KrylovNewton", test_type = "NormDispIncr", test_opt = 0, max_iter = MAX_ITER, tol = TOL, allow_smaller_step = False):
        """
        The constructor of the class.
 
        @param data_dir (str): Directory in which the results from the analysis will be stored. Use the recorders (from OpenSeesPy) or the Record method from MemberModel.
        @param name_ODB (str): Name for the folder in which the data for the animations and the fibers are stored.
        @param algo (str, optional): Type of alghoritm chosen for the analysis. It detemines how to construct a SolutionAlgorithm object, which determines the sequence of steps taken to solve the non-linear equation.
            For more information on the available types, see the OpenSeesPy documentation. Defaults to "KrylovNewton".
        @param test_type (str, optional): Type of test chosen for the analysis. It determines how to construct a ConvergenceTest object.
            Certain SolutionAlgorithm objects require a ConvergenceTest object to determine if convergence has been achieved at the end of an iteration step.
            For more information on the available types, see the OpenSeesPy documentation. Defaults to "NormDispIncr".
        @param test_opt (int, optional): Print-flag from 0 to 5 used to receive more info during the iteration
            (for example: 0 print nothing and 2 print information on norms and number of iterations at end of successful test).
            For more information, see the OpenSeesPy documentation. Defaults to 0.
        @param max_iter (float, optional): Maximal number of iterations to check. Defaults to MAX_ITER (from Constants Module).
        @param tol (float, optional): Tolerance criteria used to check for convergence. Defaults to TOL (from Constants Module).
        @param allow_smaller_step (bool, optional): Allow smaller steps in the displacement-control analysis. Defaults to False.
 
        @exception NegativeValue: The argument max_iter should be positive.
        @exception NegativeValue: The argument tol should be positive.
        """
        if max_iter < 0: raise NegativeValue()
        if tol < 0: raise NegativeValue()
        if not os.path.exists(data_dir):
            print("Folder {} not found in this directory; creating one".format(data_dir))
            os.makedirs(data_dir)
 
        self.data_dir = data_dir
        self.name_ODB = name_ODB
        self.algo = algo
        self.test_type = test_type
        self.tol = tol
        self.test_opt = test_opt
        self.max_iter = max_iter
        self.allow_smaller_step = allow_smaller_step
        self.load_case = "None"
 
 

Member Function Documentation

DeformedShape()

def DeformedShape	(		self,
			scale = 1,
			animate = False,
			dt = 0.01
	)

Method that shows the final deformed shape of the model.

It can also show the animation that shows how the model behaved during the analysis.

Parameters

scale	(int, optional): The scaling factor to magnify the deformation. The value should be adjusted for each model. Defaults to 1.
animate	(bool, optional): Option to show the animation of the model during the analysis. Defaults to False.
dt	(float, optional): The time step between every iteration. Defaults to 0.01.

Exceptions

NameError

The methods for the analysis were not called.

Definition at line 652 of file AnalysisAndPostProcessing.py.

    def DeformedShape(self, scale = 1, animate = False, dt = 0.01):
        """
        Method that shows the final deformed shape of the model. It can also show the animation that shows how the model behaved during the analysis.
 
        @param scale (int, optional): The scaling factor to magnify the deformation. The value should be adjusted for each model. Defaults to 1.
        @param animate (bool, optional): Option to show the animation of the model during the analysis. Defaults to False.
        @param dt (float, optional): The time step between every iteration. Defaults to 0.01.
        
        @exception NameError: The methods for the analysis were not called.
        """
        if self.load_case == "None": raise NameError("The analysis is not complete.")
        
        # Display deformed shape, the scaling factor needs to be adjusted for each model
        opsplt.plot_deformedshape(Model = self.name_ODB, LoadCase=self.load_case, scale = scale)
        if animate:
            opsplt.animate_deformedshape(Model = self.name_ODB, LoadCase=self.load_case, dt = dt, scale = scale)
 
 

FiberResponse()

def FiberResponse	(		self,
		int	ele_fiber_ID_analysed,
			fiber_section = 1,
			animate_stress = False,
			animate_strain = False,
			fps = 25
	)

Method that shows the final stress response of the fiber section chosen.

It can also show the animation that shows how the fiber section behaved during the analysis. The fiber ID and section needs to be recorded during the analysis, thus if the method LateralForce, Pushover or LoadingProtocol was used, the same fiber ID and section need to be used.

Parameters

ele_fiber_ID_analysed	(int): The ID of the analysed fiber. If fibers are present in the model and the user wants to save ODB data (to use in the post-processing with for example FiberResponse), assign to this argument the ID of the fiber chosen. -1 will ignore the storage of data for fibers.
fiber_section	(int, optional): The section number, i.e. the Gauss integratio number. If the fiber_ID_analysed is equal to -1, this argument is not used. Defaults to 1.
animate_stress	(bool, optional): Option to show the animation of the fiber stress during the analysis. Defaults to False.
animate_strain	(bool, optional): Option to show the animation of the fiber strain during the analysis. Defaults to False.
fps	(int, optional): Number of frame per seconds for the animations. Defaults to 25.

Exceptions

NameError

The methods for the analysis were not called.

Definition at line 670 of file AnalysisAndPostProcessing.py.

    def FiberResponse(self, ele_fiber_ID_analysed: int, fiber_section = 1, animate_stress = False, animate_strain = False, fps = 25):
        """
        Method that shows the final stress response of the fiber section chosen.
        It can also show the animation that shows how the fiber section behaved during the analysis. The fiber ID and section needs to be recorded during the analysis,
        thus if the method LateralForce, Pushover or LoadingProtocol was used, the same fiber ID and section need to be used. 
 
        @param ele_fiber_ID_analysed (int): The ID of the analysed fiber. If fibers are present in the model and the user wants to save ODB data
            (to use in the post-processing with for example FiberResponse), assign to this argument the ID of the fiber chosen.
            -1 will ignore the storage of data for fibers.
        @param fiber_section (int, optional): The section number, i.e. the Gauss integratio number.
            If the fiber_ID_analysed is equal to -1, this argument is not used. Defaults to 1.
        @param animate_stress (bool, optional): Option to show the animation of the fiber stress during the analysis. Defaults to False.
        @param animate_strain (bool, optional): Option to show the animation of the fiber strain during the analysis. Defaults to False.
        @param fps (int, optional): Number of frame per seconds for the animations. Defaults to 25.
        @exception NameError: The methods for the analysis were not called.
        """
        if self.load_case == "None": raise NameError("The analysis is not complete.")
 
        opsplt.plot_fiberResponse2D(self.name_ODB, self.load_case, ele_fiber_ID_analysed, fiber_section, InputType = 'stress')
        if animate_stress:
            ani1 = opsplt.animate_fiberResponse2D(self.name_ODB, self.load_case, ele_fiber_ID_analysed, fiber_section, InputType = 'stress', fps = fps)
        if animate_strain:
            ani1 = opsplt.animate_fiberResponse2D(self.name_ODB, self.load_case, ele_fiber_ID_analysed, fiber_section, InputType = 'strain', fps = fps)
 
 

Gravity()

def Gravity	(		self,
		list	loaded_nodes,
		list	Fy,
		int	timeSeries_ID,
		int	pattern_ID,
			n_step = 10,
			timeSeries_type = "Linear",
			pattern_type = "Plain",
			constraints_type = "Plain",
			numberer_type = "RCM",
			system_type = "BandGeneral",
			analysis_type = "Static",
			show_plot = False
	)

Method to perform the gravity analyisis with vertical loadings (load-control).

It can be used before calling the Pushover or LoadingProtocol methods that perform the actual anlysis. If no vertical loadings present, this method can be avoided.

Parameters

loaded_nodes	(list): List of nodes that are loaded by the the forces in Fy. The first node will be recorded (thus usually should be in the roof).
Fy	(list): List of vertical loadings (negative is toward the ground, thus compression; see global coordinate system).
timeSeries_ID	(int): ID of the timeseries.
pattern_ID	(int): ID of the pattern.
n_step	(int, optional): Number of steps used to during the analysis to reach the objective state (with 100% vertical loadings imposed). Defaults to 10.
timeSeries_type	(str, optional): Type of timeseries chosen. For more information, see the OpenSeesPy documentation. Defaults to "Linear".
pattern_type	(str, optional): Type of pattern chosen. For more information, see the OpenSeesPy documentation. Defaults to "Plain".
constraints_type	(str, optional): Type of contraints chosen. It detemines how the constraint equations are enforced in the analysis. For more information, see the OpenSeesPy documentation. Defaults to "Plain".
numberer_type	(str, optional): Type of numberer chosen. It determines the mapping between equation numbers and degrees-of-freedom. For more information, see the OpenSeesPy documentation. Defaults to "RCM".
system_type	(str, optional): Type of system of equations chosen. It determines how to construct the LinearSOE and LinearSolver objects to store and solve the system of equations in the analysis. For more information, see the OpenSeesPy documentation. Defaults to "BandGeneral".
analysis_type	(str, optional): Type of analysis chosen. It determines how to construct the Analysis object, which defines what type of analysis is to be performed. For more information, see the OpenSeesPy documentation. Defaults to "Static".
show_plot	(bool, optional): Option to show the 'vertical displacement vs. vertical loading' curve after the analysis. Defaults to False.

Exceptions

WrongDimension	The dimension of the loaded_nodes and Fy arguments needs to be the same.
NegativeValue	The ID of timeSeries_ID needs to be a positive integer.
NegativeValue	The ID of pattern_ID needs to be a positive integer.

Definition at line 60 of file AnalysisAndPostProcessing.py.

        constraints_type = "Plain", numberer_type = "RCM", system_type = "BandGeneral", analysis_type = "Static", show_plot = False):
        """
        Method to perform the gravity analyisis with vertical loadings (load-control).
        It can be used before calling the Pushover or LoadingProtocol methods that perform the actual anlysis. If no vertical loadings present, this method can be avoided.
 
        @param loaded_nodes (list): List of nodes that are loaded by the the forces in Fy. The first node will be recorded (thus usually should be in the roof).
        @param Fy (list): List of vertical loadings (negative is toward the ground, thus compression; see global coordinate system).
        @param timeSeries_ID (int): ID of the timeseries.
        @param pattern_ID (int): ID of the pattern.
        @param n_step (int, optional): Number of steps used to during the analysis to reach the objective state (with 100% vertical loadings imposed). Defaults to 10.
        @param timeSeries_type (str, optional): Type of timeseries chosen.
            For more information, see the OpenSeesPy documentation. Defaults to "Linear".
        @param pattern_type (str, optional): Type of pattern chosen.
            For more information, see the OpenSeesPy documentation. Defaults to "Plain".
        @param constraints_type (str, optional): Type of contraints chosen. It detemines how the constraint equations are enforced in the analysis.
            For more information, see the OpenSeesPy documentation. Defaults to "Plain".
        @param numberer_type (str, optional): Type of numberer chosen. It determines the mapping between equation numbers and degrees-of-freedom.
            For more information, see the OpenSeesPy documentation. Defaults to "RCM".
        @param system_type (str, optional): Type of system of equations chosen. It determines how to construct the LinearSOE and LinearSolver objects to store and solve the system of equations in the analysis.
            For more information, see the OpenSeesPy documentation. Defaults to "BandGeneral".
        @param analysis_type (str, optional): Type of analysis chosen. It determines how to construct the Analysis object, which defines what type of analysis is to be performed.
            For more information, see the OpenSeesPy documentation. Defaults to "Static".
        @param show_plot (bool, optional): Option to show the 'vertical displacement vs. vertical loading' curve after the analysis. Defaults to False.
 
        @exception WrongDimension: The dimension of the loaded_nodes and Fy arguments needs to be the same.
        @exception NegativeValue: The ID of timeSeries_ID needs to be a positive integer.
        @exception NegativeValue: The ID of pattern_ID needs to be a positive integer.
        """
        if len(loaded_nodes) != len(Fy): raise WrongDimension()
        if timeSeries_ID < 1: raise NegativeValue()
        if pattern_ID < 1: raise NegativeValue()
 
        # for mass defined: opsplt.createODB(self.name_ODB, "Gravity", Nmodes = nEigen); 
        # for tracking gravity with ODB: opsplt.createODB(self.name_ODB, "Gravity");
 
        # Create load pattern
        timeSeries(timeSeries_type, timeSeries_ID)
        pattern(pattern_type, timeSeries_ID, pattern_ID)
        for ii, node_ID in enumerate(loaded_nodes):
            load(node_ID, 0.0, Fy[ii], 0.0)     # load(IDNode, Fx, Fy, Mz)
        DGravity = 1.0/n_step                   # load increment                  
 
        # Set up analysis options
        constraints(constraints_type)   # how it handles boundary conditions
        numberer(numberer_type)         # renumber dof's to minimize band-width (optimization)
        system(system_type)             # how to store and solve the system of equations in the analysis
                                        # For static model, BandGeneral, for transient and/or big model, UmfPack
        integrator("LoadControl", DGravity) # LoadControl and DisplacementControl only with static model, linear TimeSeries w/ factor of 1
                                            # Newmark used for transient model  
        algorithm("Newton")             # placeholder
        analysis(analysis_type)         # define type of analysis: static for pushover
 
        # Analysis
        dataG = np.zeros((n_step+1,2))
        print("")
        print("Gravity analysis starts")
        for iteration in range(n_step):
            convergence = self.__LoadCtrlLoop(DGravity, iteration,
                self.algo, self.test_type, self.tol, self.test_opt, self.max_iter)
            if convergence != 0: break
            dataG[iteration+1,0] = nodeDisp(loaded_nodes[0], 2)/mm_unit
            dataG[iteration+1,1] = getLoadFactor(pattern_ID)*Fy[0]/kN_unit
 
        if show_plot:
            plt.plot(dataG[:,0], dataG[:,1])
            plt.xlabel('Vertical Displacement [mm]')
            plt.ylabel('Vertical Load [kN]')
            plt.title('Gravity curve')
            plt.show()
 
        loadConst("-time", 0.0)
 
        print("")
        print("Gravity complete")
 
 

LateralForce()

def LateralForce	(		self,
		list	loaded_nodes,
		list	Fx,
		int	timeSeries_ID,
		int	pattern_ID,
			n_step = 1000,
			fiber_ID_analysed = -1,
			fiber_section = 1,
			timeSeries_type = "Linear",
			pattern_type = "Plain",
			constraints_type = "Plain",
			numberer_type = "RCM",
			system_type = "BandGeneral",
			analysis_type = "Static",
			show_plot = True,
			block = False
	)

Method to perform the lateral force analyisis with lateral loading (load-control).

If this method is called, the LoadingProtocol and Pushover methods should be avoided.

Parameters

loaded_nodes	(list): List of nodes that are loaded by the the forces in Fx. The first node will be recorded (thus usually should be in the roof).
Fx	(list): List of horizontal loadings (negative is toward left; see global coordinate system).
timeSeries_ID	(int): ID of the timeseries.
pattern_ID	(int): ID of the pattern.
n_step	(int, optional): Number of steps used to during the analysis to reach the objective state (with 100% horizontal loadings imposed). Defaults to 1000.
fiber_ID_analysed	(int, optional): The ID of the analysed fiber. If fibers are present in the model and the user wants to save ODB data (to use in the post-processing with for example FiberResponse), assign to this argument the ID of the fiber chosen. -1 will ignore the storage of data for fibers. Defaults to -1.
fiber_section	(int, optional): The section number, i.e. the Gauss integratio number. If the fiber_ID_analysed is equal to -1, this argument is not used. Defaults to 1.
timeSeries_type	(str, optional): Type of timeseries chosen. For more information, see the OpenSeesPy documentation. Defaults to "Linear".
pattern_type	(str, optional): Type of pattern chosen. For more information, see the OpenSeesPy documentation. Defaults to "Plain".
constraints_type	(str, optional): Type of contraints chosen. It detemines how the constraint equations are enforced in the analysis. For more information, see the OpenSeesPy documentation. Defaults to "Plain".
numberer_type	(str, optional): Type of numberer chosen. It determines the mapping between equation numbers and degrees-of-freedom. For more information, see the OpenSeesPy documentation. Defaults to "RCM".
system_type	(str, optional): Type of system of equations chosen. It determines how to construct the LinearSOE and LinearSolver objects to store and solve the system of equations in the analysis. For more information, see the OpenSeesPy documentation. Defaults to "BandGeneral".
analysis_type	(str, optional): Type of analysis chosen. It determines how to construct the Analysis object, which defines what type of analysis is to be performed. For more information, see the OpenSeesPy documentation. Defaults to "Static".
show_plot	(bool, optional): Option to show the 'Horizontal displacement vs. Horizontal loading' curve after the analysis. Defaults to True.
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.

Exceptions

WrongDimension	The dimension of the loaded_nodes and Fx arguments needs to be the same.
NegativeValue	The ID of timeSeries_ID needs to be a positive integer.
NegativeValue	The ID of pattern_ID needs to be a positive integer.
NegativeValue	The ID of fiber_ID_analysed needs to be a positive integer.

Definition at line 137 of file AnalysisAndPostProcessing.py.

        show_plot = True, block = False):
        """
        Method to perform the lateral force analyisis with lateral loading (load-control).
        If this method is called, the LoadingProtocol and Pushover methods should be avoided.
 
        @param loaded_nodes (list): List of nodes that are loaded by the the forces in Fx. The first node will be recorded (thus usually should be in the roof).
        @param Fx (list): List of horizontal loadings (negative is toward left; see global coordinate system).
        @param timeSeries_ID (int): ID of the timeseries.
        @param pattern_ID (int): ID of the pattern.
        @param n_step (int, optional): Number of steps used to during the analysis to reach the objective state (with 100% horizontal loadings imposed). Defaults to 1000.
        @param fiber_ID_analysed (int, optional): The ID of the analysed fiber. If fibers are present in the model and the user wants to save ODB data
            (to use in the post-processing with for example FiberResponse), assign to this argument the ID of the fiber chosen.
            -1 will ignore the storage of data for fibers. Defaults to -1.
        @param fiber_section (int, optional): The section number, i.e. the Gauss integratio number.
            If the fiber_ID_analysed is equal to -1, this argument is not used. Defaults to 1.
        @param timeSeries_type (str, optional): Type of timeseries chosen.
            For more information, see the OpenSeesPy documentation. Defaults to "Linear".
        @param pattern_type (str, optional): Type of pattern chosen.
            For more information, see the OpenSeesPy documentation. Defaults to "Plain".
        @param constraints_type (str, optional): Type of contraints chosen. It detemines how the constraint equations are enforced in the analysis.
            For more information, see the OpenSeesPy documentation. Defaults to "Plain".
        @param numberer_type (str, optional): Type of numberer chosen. It determines the mapping between equation numbers and degrees-of-freedom.
            For more information, see the OpenSeesPy documentation. Defaults to "RCM".
        @param system_type (str, optional): Type of system of equations chosen. It determines how to construct the LinearSOE and LinearSolver objects to store and solve the system of equations in the analysis.
            For more information, see the OpenSeesPy documentation. Defaults to "BandGeneral".
        @param analysis_type (str, optional): Type of analysis chosen. It determines how to construct the Analysis object, which defines what type of analysis is to be performed.
            For more information, see the OpenSeesPy documentation. Defaults to "Static".
        @param show_plot (bool, optional): Option to show the 'Horizontal displacement vs. Horizontal loading' curve after the analysis. Defaults to True.
        @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.
 
        @exception WrongDimension: The dimension of the loaded_nodes and Fx arguments needs to be the same.
        @exception NegativeValue: The ID of timeSeries_ID needs to be a positive integer.
        @exception NegativeValue: The ID of pattern_ID needs to be a positive integer.
        @exception NegativeValue: The ID of fiber_ID_analysed needs to be a positive integer.
        """
        if len(loaded_nodes) != len(Fx): raise WrongDimension()
        if timeSeries_ID < 1: raise NegativeValue()
        if pattern_ID < 1: raise NegativeValue()
        if fiber_ID_analysed != -1 and fiber_ID_analysed < 1: raise NegativeValue()
 
        # for mass defined: opsplt.createODB(self.name_ODB, "LateralForce", Nmodes = nEigen); 
        opsplt.createODB(self.name_ODB, "LateralForce");
        if fiber_ID_analysed != -1: opsplt.saveFiberData2D(self.name_ODB, "LateralForce", fiber_ID_analysed, fiber_section)
 
        # Create load pattern
        timeSeries(timeSeries_type, timeSeries_ID)
        pattern(pattern_type, timeSeries_ID, pattern_ID)
        for ii, node_ID in enumerate(loaded_nodes):
            load(node_ID, Fx[ii], 0.0, 0.0)     # load(IDNode, Fx, Fy, Mz)
        force = 1.0/n_step                   # load increment                  
 
        # Set up analysis options
        constraints(constraints_type)   # how it handles boundary conditions
        numberer(numberer_type)         # renumber dof's to minimize band-width (optimization)
        system(system_type)             # how to store and solve the system of equations in the analysis
                                        # For static model, BandGeneral, for transient and/or big model, UmfPack
        integrator("LoadControl", force)# LoadControl and DisplacementControl only with static model, linear TimeSeries w/ factor of 1
                                        # Newmark used for transient model  
        algorithm("Newton")             # placeholder
        analysis(analysis_type)         # define type of analysis: static for pushover
 
        # Analysis
        dataLF = np.zeros((n_step+1,2))
        print("")
        print("Lateral Force analysis starts")
        for iteration in range(n_step):
            convergence = self.__LoadCtrlLoop(force, iteration,
                self.algo, self.test_type, self.tol, self.test_opt, self.max_iter)
            if convergence != 0: break
            dataLF[iteration+1,0] = nodeDisp(loaded_nodes[0], 1)/mm_unit
            dataLF[iteration+1,1] = getLoadFactor(pattern_ID)*Fx[0]/kN_unit
 
        if show_plot:
            plt.plot(dataLF[:,0], dataLF[:,1])
            plt.xlabel('Lateral Displacement [mm]')
            plt.ylabel('Lateral Load [kN]')
            plt.title('Lateral force curve')
            if block:
                plt.show()
 
        loadConst("-time", 0.0)
 
        print("")
        print("Lateral force complete")
        self.load_case = "LateralForce"
 
        wipe()
 
 

LoadingProtocol()

def LoadingProtocol	(		self,
		int	CtrlNode,
		np.ndarray	discr_LP,
		int	timeSeries_ID,
		int	pattern_ID,
			Fx = 1*kN_unit,
			ele_fiber_ID_analysed = -1,
			fiber_section = 1,
			timeSeries_type = "Linear",
			pattern_type = "Plain",
			constraints_type = "Plain",
			numberer_type = "RCM",
			system_type = "UmfPack",
			analysis_type = "Static",
			show_plot = True,
			block = False
	)

Method to perform a loading protocol analysis (displacement-control).

If this method is called, the Pushover and LateralForce methods should be avoided.

Parameters

CtrlNode	(int): The node that will be used to impose the displacement from the discr_LP to perform the analysis.
discr_LP	(np.ndarray): The loading protocol array (1 dimension) discretised. It needs to be filled with imposed displacement, not SDR. Use the functions DiscretizeLoadProtocol and DiscretizeLinearly in FunctionalFeatures module to help create and/or discretise one.
timeSeries_ID	(int): ID of the timeseries.
pattern_ID	(int): ID of the pattern.
Fx	(float, optional): The force imposed at the control node CtrlNode. It is used for convergence reasons and it can be arbitrarly small. Defaults to 1*kN_unit.
ele_fiber_ID_analysed	(int, optional): The ID of the analysed element with fibers. If fibers are present in the model and the user wants to save ODB data (to use in the post-processing with for example FiberResponse), assign to this argument the ID of the element with fibers chosen. -1 will ignore the storage of data for fibers. Defaults to -1.
fiber_section	(int, optional): The section number, i.e. the Gauss integratio number. If the fiber_ID_analysed is equal to -1, this argument is not used. Defaults to 1.
timeSeries_type	(str, optional): Type of timeseries chosen. For more information, see the OpenSeesPy documentation. Defaults to "Linear".
pattern_type	(str, optional): Type of pattern chosen. For more information, see the OpenSeesPy documentation. Defaults to "Plain".
constraints_type	(str, optional): Type of contraints chosen. It detemines how the constraint equations are enforced in the analysis. For more information, see the OpenSeesPy documentation. Defaults to "Plain".
numberer_type	(str, optional): Type of numberer chosen. It determines the mapping between equation numbers and degrees-of-freedom. For more information, see the OpenSeesPy documentation. Defaults to "RCM".
system_type	(str, optional): Type of system of equations chosen. It determines how to construct the LinearSOE and LinearSolver objects to store and solve the system of equations in the analysis. For more information, see the OpenSeesPy documentation. Defaults to "UmfPack".
analysis_type	(str, optional): Type of analysis chosen. It determines how to construct the Analysis object, which defines what type of analysis is to be performed. For more information, see the OpenSeesPy documentation. Defaults to "Static".
show_plot	(bool, optional): Option to show the 'lateral displacement vs. lateral loading' curve after the analysis. Defaults to True.
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.

Exceptions

NegativeValue	The ID of CtrlNode needs to be a positive integer.
NegativeValue	The ID of timeSeries_ID needs to be a positive integer.
NegativeValue	The ID of pattern_ID needs to be a positive integer.
NegativeValue	The ID of fiber_ID_analysed needs to be a positive integer if is different from -1.

Definition at line 319 of file AnalysisAndPostProcessing.py.

        show_plot = True, block = False):
        """
        Method to perform a loading protocol analysis (displacement-control). If this method is called, the Pushover and LateralForce methods should be avoided.
 
        @param CtrlNode (int): The node that will be used to impose the displacement from the discr_LP to perform the analysis.
        @param discr_LP (np.ndarray): The loading protocol array (1 dimension) discretised. It needs to be filled with imposed displacement, not SDR.
            Use the functions DiscretizeLoadProtocol and DiscretizeLinearly in FunctionalFeatures module to help create and/or discretise one.
        @param timeSeries_ID (int): ID of the timeseries.
        @param pattern_ID (int): ID of the pattern.
        @param Fx (float, optional): The force imposed at the control node CtrlNode. It is used for convergence reasons and it can be arbitrarly small.
            Defaults to 1*kN_unit.
        @param ele_fiber_ID_analysed (int, optional): The ID of the analysed element with fibers. If fibers are present in the model and the user wants to save ODB data
            (to use in the post-processing with for example FiberResponse), assign to this argument the ID of the element with fibers chosen.
            -1 will ignore the storage of data for fibers. Defaults to -1.
        @param fiber_section (int, optional): The section number, i.e. the Gauss integratio number.
            If the fiber_ID_analysed is equal to -1, this argument is not used. Defaults to 1.
        @param timeSeries_type (str, optional): Type of timeseries chosen.
            For more information, see the OpenSeesPy documentation. Defaults to "Linear".
        @param pattern_type (str, optional): Type of pattern chosen.
            For more information, see the OpenSeesPy documentation. Defaults to "Plain".
        @param constraints_type (str, optional): Type of contraints chosen. It detemines how the constraint equations are enforced in the analysis.
            For more information, see the OpenSeesPy documentation. Defaults to "Plain".
        @param numberer_type (str, optional): Type of numberer chosen. It determines the mapping between equation numbers and degrees-of-freedom.
            For more information, see the OpenSeesPy documentation. Defaults to "RCM".
        @param system_type (str, optional): Type of system of equations chosen. It determines how to construct the LinearSOE and LinearSolver objects to store and solve the system of equations in the analysis.
            For more information, see the OpenSeesPy documentation. Defaults to "UmfPack".
        @param analysis_type (str, optional): Type of analysis chosen. It determines how to construct the Analysis object, which defines what type of analysis is to be performed.
            For more information, see the OpenSeesPy documentation. Defaults to "Static".
        @param show_plot (bool, optional): Option to show the 'lateral displacement vs. lateral loading' curve after the analysis. Defaults to True.
        @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.
 
        @exception NegativeValue: The ID of CtrlNode needs to be a positive integer.
        @exception NegativeValue: The ID of timeSeries_ID needs to be a positive integer.
        @exception NegativeValue: The ID of pattern_ID needs to be a positive integer.
        @exception NegativeValue: The ID of fiber_ID_analysed needs to be a positive integer if is different from -1.
        """
        if CtrlNode < 1: raise NegativeValue()
        if timeSeries_ID < 1: raise NegativeValue()
        if pattern_ID < 1: raise NegativeValue()
        if ele_fiber_ID_analysed != -1 and ele_fiber_ID_analysed < 1: raise NegativeValue()
 
        # for mass defined: opsplt.createODB(self.name_ODB, "LoadingProtocol", Nmodes = nEigen); 
        opsplt.createODB(self.name_ODB, "LoadingProtocol");
        if ele_fiber_ID_analysed != -1: opsplt.saveFiberData2D(self.name_ODB, "LoadingProtocol", ele_fiber_ID_analysed, fiber_section)
 
        # Create load pattern
        timeSeries(timeSeries_type, timeSeries_ID)
        pattern(pattern_type, timeSeries_ID, pattern_ID)
        load(CtrlNode, Fx, 0.0, 0.0)    # load(IDNode, Fx, Fy, Mz)
        dU_prev = 0
        Nsteps = np.size(discr_LP)      # number of pushover analysis steps
 
        # Set up analysis options
        constraints(constraints_type)   # how it handles boundary conditions
        numberer(numberer_type)         # renumber dof's to minimize band-width (optimization)
        system(system_type)             # how to store and solve the system of equations in the analysis
                                        # For static model, BandGeneral, for transient and/or big model, UmfPack
        integrator("LoadControl", 1)    # placeholder
        algorithm("Newton")             # placeholder
        analysis(analysis_type)         # define type of analysis: static for LoadingProtocol
 
        # Analysis
        dataLP = np.zeros((Nsteps+1,2))
        next_step = 0
        print("")
        print("Loading Protocol analysis starts")
        for iteration in range(Nsteps):
            # Compute displacement usinf the given loading protocol (discretized)
            dU_next = discr_LP[iteration]
            dU = dU_next - dU_prev
            dU_prev = dU_next
            
            next_step = ProgressingPercentage(Nsteps, iteration, next_step)
            convergence = self.__LatDispCtrlLoop(CtrlNode, dU, iteration, 
                self.algo, self.test_type, self.tol, self.test_opt, self.max_iter, self.allow_smaller_step)
            if convergence != 0: break
            dataLP[iteration+1,0] = nodeDisp(CtrlNode, 1)/mm_unit
            dataLP[iteration+1,1] = getLoadFactor(pattern_ID)*Fx/kN_unit
        
        if show_plot:
            plt.plot(dataLP[:,0], dataLP[:,1])
            plt.xlabel('Horizontal Displacement [mm]')
            plt.ylabel('Horizontal Load [kN]')
            plt.title('Loading Protocol curve')
            if block:
                plt.show()
        
        print("")
        print("Loading Protocol complete")
        self.load_case = "LoadingProtocol"
        
        wipe()
 
 

Pushover()

def Pushover	(		self,
		int	CtrlNode,
			Dmax,
			Dincr,
		int	timeSeries_ID,
		int	pattern_ID,
			Fx = 1*kN_unit,
			ele_fiber_ID_analysed = -1,
			fiber_section = 1,
			timeSeries_type = "Linear",
			pattern_type = "Plain",
			constraints_type = "Plain",
			numberer_type = "RCM",
			system_type = "UmfPack",
			analysis_type = "Static",
			show_plot = True,
			block = False
	)

Method to perform a pushover analysis (displacement-control).

If this method is called, the LoadingProtocol and LateralForce methods should be avoided.

Parameters

CtrlNode	(int): The node that will be used to impose the displacement Dmax of the pushover analysis. If the show_plot option is True, the curve displayed follows this node.
Dmax	(float): The imposed displacement.
Dincr	(float): The incremental displacement to reach Dmax. To converge, it should be small enough (1000 times smaller of Dmax).
timeSeries_ID	(int): ID of the timeseries.
pattern_ID	(int): ID of the pattern.
Fx	(float, optional): The force imposed at the control node CtrlNode. It is used for convergence reasons and it can be arbitrarly small. Defaults to 1*kN_unit.
ele_fiber_ID_analysed	(int, optional): The ID of the analysed element with fibers. If fibers are present in the model and the user wants to save ODB data (to use in the post-processing with for example FiberResponse), assign to this argument the ID of the element with fibers chosen. -1 will ignore the storage of data for fibers. Defaults to -1.
fiber_section	(int, optional): The section number, i.e. the Gauss integratio number. If the fiber_ID_analysed is equal to -1, this argument is not used. Defaults to 1.
timeSeries_type	(str, optional): Type of timeseries chosen. For more information, see the OpenSeesPy documentation. Defaults to "Linear".
pattern_type	(str, optional): Type of pattern chosen. For more information, see the OpenSeesPy documentation. Defaults to "Plain".
constraints_type	(str, optional): Type of contraints chosen. It detemines how the constraint equations are enforced in the analysis. For more information, see the OpenSeesPy documentation. Defaults to "Plain".
numberer_type	(str, optional): Type of numberer chosen. It determines the mapping between equation numbers and degrees-of-freedom. For more information, see the OpenSeesPy documentation. Defaults to "RCM".
system_type	(str, optional): Type of system of equations chosen. It determines how to construct the LinearSOE and LinearSolver objects to store and solve the system of equations in the analysis. For more information, see the OpenSeesPy documentation. Defaults to "UmfPack".
analysis_type	(str, optional): Type of analysis chosen. It determines how to construct the Analysis object, which defines what type of analysis is to be performed. For more information, see the OpenSeesPy documentation. Defaults to "Static".
show_plot	(bool, optional): Option to show the 'lateral displacement vs. lateral loading' curve after the analysis. Defaults to True.
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.

Exceptions

NegativeValue	The ID of CtrlNode needs to be a positive integer.
NegativeValue	The ID of timeSeries_ID needs to be a positive integer.
NegativeValue	The ID of pattern_ID needs to be a positive integer.
NegativeValue	The ID of ele_fiber_ID_analysed needs to be a positive integer if is different from -1.

Definition at line 228 of file AnalysisAndPostProcessing.py.

        show_plot = True, block = False):
        """
        Method to perform a pushover analysis (displacement-control). If this method is called, the LoadingProtocol and LateralForce methods should be avoided.
 
        @param CtrlNode (int): The node that will be used to impose the displacement Dmax of the pushover analysis.
            If the show_plot option is True, the curve displayed follows this node.
        @param Dmax (float): The imposed displacement.
        @param Dincr (float): The incremental displacement to reach Dmax. To converge, it should be small enough (1000 times smaller of Dmax).
        @param timeSeries_ID (int): ID of the timeseries.
        @param pattern_ID (int): ID of the pattern.
        @param Fx (float, optional): The force imposed at the control node CtrlNode. It is used for convergence reasons and it can be arbitrarly small.
            Defaults to 1*kN_unit.
        @param ele_fiber_ID_analysed (int, optional): The ID of the analysed element with fibers. If fibers are present in the model and the user wants to save ODB data
            (to use in the post-processing with for example FiberResponse), assign to this argument the ID of the element with fibers chosen.
            -1 will ignore the storage of data for fibers. Defaults to -1.
        @param fiber_section (int, optional): The section number, i.e. the Gauss integratio number.
            If the fiber_ID_analysed is equal to -1, this argument is not used. Defaults to 1.
        @param timeSeries_type (str, optional): Type of timeseries chosen.
            For more information, see the OpenSeesPy documentation. Defaults to "Linear".
        @param pattern_type (str, optional): Type of pattern chosen.
            For more information, see the OpenSeesPy documentation. Defaults to "Plain".
        @param constraints_type (str, optional): Type of contraints chosen. It detemines how the constraint equations are enforced in the analysis.
            For more information, see the OpenSeesPy documentation. Defaults to "Plain".
        @param numberer_type (str, optional): Type of numberer chosen. It determines the mapping between equation numbers and degrees-of-freedom.
            For more information, see the OpenSeesPy documentation. Defaults to "RCM".
        @param system_type (str, optional): Type of system of equations chosen. It determines how to construct the LinearSOE and LinearSolver objects to store and solve the system of equations in the analysis.
            For more information, see the OpenSeesPy documentation. Defaults to "UmfPack".
        @param analysis_type (str, optional): Type of analysis chosen. It determines how to construct the Analysis object, which defines what type of analysis is to be performed.
            For more information, see the OpenSeesPy documentation. Defaults to "Static".
        @param show_plot (bool, optional): Option to show the 'lateral displacement vs. lateral loading' curve after the analysis. Defaults to True.
        @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.
 
        @exception NegativeValue: The ID of CtrlNode needs to be a positive integer.
        @exception NegativeValue: The ID of timeSeries_ID needs to be a positive integer.
        @exception NegativeValue: The ID of pattern_ID needs to be a positive integer.
        @exception NegativeValue: The ID of ele_fiber_ID_analysed needs to be a positive integer if is different from -1.
        """
        if CtrlNode < 1: raise NegativeValue()
        if timeSeries_ID < 1: raise NegativeValue()
        if pattern_ID < 1: raise NegativeValue()
        if ele_fiber_ID_analysed != -1 and ele_fiber_ID_analysed < 1: raise NegativeValue()
 
        # for mass defined: opsplt.createODB(self.name_ODB, "Pushover", Nmodes = nEigen); 
        opsplt.createODB(self.name_ODB, "Pushover");
        if ele_fiber_ID_analysed != -1: opsplt.saveFiberData2D(self.name_ODB, "Pushover", ele_fiber_ID_analysed, fiber_section)
 
        # Create load pattern
        timeSeries(timeSeries_type, timeSeries_ID)
        pattern(pattern_type, timeSeries_ID, pattern_ID)
        load(CtrlNode, Fx, 0.0, 0.0)    # load(IDNode, Fx, Fy, Mz)
        Nsteps = int(abs(Dmax/Dincr))   # number of pushover analysis steps
        
        # Set up analysis options
        constraints(constraints_type)   # how it handles boundary conditions
        numberer(numberer_type)         # renumber dof's to minimize band-width (optimization)
        system(system_type)             # how to store and solve the system of equations in the analysis
                                        # For static model, BandGeneral, for transient and/or big model, UmfPack
        integrator("LoadControl", 1)    # placeholder
        algorithm("Newton")             # placeholder
        analysis(analysis_type)         # define type of analysis: static for pushover
 
        # Analysis
        dataPO = np.zeros((Nsteps+1,2))
        next_step = 0
        print("")
        print("Pushover analysis starts")
        for iteration in range(Nsteps):
            next_step = ProgressingPercentage(Nsteps, iteration, next_step)
            convergence = self.__LatDispCtrlLoop(CtrlNode, Dincr, iteration,
                self.algo, self.test_type, self.tol, self.test_opt, self.max_iter, self.allow_smaller_step)
            if convergence != 0: break
            dataPO[iteration+1,0] = nodeDisp(CtrlNode, 1)/mm_unit
            dataPO[iteration+1,1] = getLoadFactor(pattern_ID)*Fx/kN_unit
        
        if show_plot:
            plt.plot(dataPO[:,0], dataPO[:,1])
            plt.xlabel('Horizontal Displacement [mm]')
            plt.ylabel('Horizontal Load [kN]')
            plt.title('Pushover curve')
            if block:
                plt.show()
        
        print("")
        print("Pushover complete")
        self.load_case = "Pushover"
 
        wipe()
 
 

Member Data Documentation

algo

algo

Definition at line 51 of file AnalysisAndPostProcessing.py.

allow_smaller_step

allow_smaller_step

Definition at line 56 of file AnalysisAndPostProcessing.py.

data_dir

data_dir

Definition at line 49 of file AnalysisAndPostProcessing.py.

load_case

load_case

Definition at line 57 of file AnalysisAndPostProcessing.py.

max_iter

max_iter

Definition at line 55 of file AnalysisAndPostProcessing.py.

name_ODB

name_ODB

Definition at line 50 of file AnalysisAndPostProcessing.py.

test_opt

test_opt

Definition at line 54 of file AnalysisAndPostProcessing.py.

test_type

test_type

Definition at line 52 of file AnalysisAndPostProcessing.py.

tol

tol

Definition at line 53 of file AnalysisAndPostProcessing.py.

---

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

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