opseestools 0.60

A collection of OpenSeesPy routines for performing several types of analyses and other tools

opseestools

This is a library to be used with OpenSeesPy.

It includes a collection of analysis functions that are useful if you are analyzing frames or wall buildings both in two and three dimensions, storing the results in Python variables. It also includes another useful tools for processing results and creating scripts

Installation

pip install opseestools

Description and use

ops tools comprises a set of functions in four categories:

1. 2D Analysis functions. This library is called analysis. You can import it as:

import opstools.analisis as an

2. 3D Analysis functions. This library is called analysis. You can import it as:

import opstools.analisis3D as an3D (or any name you want)

3. A set of utilities to support model building. You can import it as:

import opstools.utilidades as ut (or any other name)

4. A library to calculate and plot fragility functions. You can import it as:

import opstools.Lib_frag as lf (or any other name)

The full documentation is available on the GitHub website: https://github.com/odarroyo/opseestools Tutorials of how to use the functions are available in the examples folder in: https://github.com/odarroyo/opseestools

2D Analysis functions

Here's a brief description of some of the functions:

gravedad()

To perform gravity analysis of a model you created. In a more general sense, it runs a static analysis, but returns no output. The idea of this script is to perform the gravity analysis before a pushover or other nonlinear analysis.

pushover2(Dmax,Dincr,IDctrlNode,IDctrlDOF)

This function performs a pushover analysis. Dmax = maximum displacement Dincr = increment of displacement IDctrlNode = control node for the pushover IDctrlDOF = DOF for pushover

The function returns two variables as numpy arrays: roof displacement and base shear. In that order. An example call of this function is:

roof,shear = an.pushover2(1.5,0.01,9,1), which performs a pushover analysis using node 9 as a control node in the DOF = 1 until a displacement of 1.5 in 0.01 increments

dinamicoIDA4P(recordName,dtrec,nPts,dtan,fact,damp,IDctrlNode,IDctrlDOF,elements,nodes_control,modes = [0,2],Kswitch = 1,Tol=1e-4)

This is a function to be used with any 2D model where each element has six degrees of freedom, e.g., any 2D frame (steel or concrete), with or without infills (as long as they're modeled as 6DOF elements) should work.

Inputs are:

• recordName: name of the txt file with the values of acceleration. e.g., 'GM01.txt'. Only one value per line
• dtrec: dt of the record, i.e, time interval in seconds
• nPts: number of points of the file
• dtan: dt to be used by OpenSeesPy in the analysis
• fact: scale factor to be applied to the record
• damp: damping ratio to be applied to the model
• IDctrlNode = control node for the pushover
• IDctrlDOF = DOF for pushover
• elements = list with elements to record forces. e.g, [3,5,6] records forces for elements 3, 5 and 6
• nodes_control = nodes to record displacement in the DOF direction. This is meant to calculate interstory drifts, thereforce you should input one node per building floor
• modes = list to include the number of modes for computing the damping in a Rayleight scheme. Default is [0,2], which uses the first and third mode. Adjust according to your needs.
• Kswitch = 1 to use initial stiffness, 2 for current stiffness. Default is 1
• Tol = tolerance for equilibrium. The algorithm uses the NormUnbalance test, which is appropriate for most circumstances and a 1e-4 value. Adjust the tolerance according to your problem size.

Outputs are:

• time: time of the analysis
• roof: roof displacement
• element_forces: numpy 3D array with the element forces. Think of this as a stack of 2D matrices, where each 2D matrix contains the response of one element: in the rows the analysis time time and in the columns the element DOF.
• node_disp: 2D matrix where each row is time and columns have the node displacement in the specified DOF
• node_vel: same as node_disp but for velocity
• node_acel: same as node_disp but for acceleration. Returns the floor relative acceleration.
• drift: same as node_disp but for interstory drifts.

The outputs should all be specified and in that order. Example:

time,roof,element_forces,node_disp,node_vel,node_acel,drift = an.dinamicoIDA4P(....)