7.2.4. Energy in Seismic Damper

A hands off example showing how to get cumulative energy in a seismic damper.The following data is used in the example:

Input xy data here.
Input load protocol here.
Reported energy in hysteresis here.

# -*- coding: utf-8 -*-
"""
An example where data from a siesmic damper is processed.

@author: Christian
"""

import numpy as np
import hysteresis as hys
import matplotlib.pyplot as plt

# load the .csv Data
InputName    = 'Hys1.csv'
LPName       = 'LoadProtocol.csv'
EnergyName   = 'Hys1_Energy.csv'
xyData       = np.loadtxt(InputName,         delimiter=',')
LoadProtocol = np.loadtxt(LPName,  delimiter=',',skiprows = 1)
Energy       = np.loadtxt(EnergyName,  delimiter=',')

# some small adjustments to account for the elastic portion of the curve that was skipped.
xyData          = xyData[:-21,:]
xyData[-1,0]    = xyData[-1,0] + 0.6

# Create the trace Hysteresis
hysTrace = hys.Hysteresis(xyData)
fig, ax  = hysTrace.initFig()
hysTrace.plotLoadProtocol()
fig, ax  = hysTrace.initFig()
hysTrace.plot(True)


# Expand the hysteresis Trace
FullHys = hys.exandHysTrace(hysTrace, np.array([3,3,3,3,3,3,3]), skipStart = 1, skipEnd = 0)
fig, ax = hysTrace.initFig()
FullHys.plot()
 
# Plot the load Protocol
fig, ax = hysTrace.initFig()
FullHys.plotLoadProtocol()
hysProt = FullHys.loadProtocol

# Get cumulative displacement and area
cumDisp = FullHys.getCumDisp() +10
cumArea = FullHys.getCumArea()/1000

# Make a nice plot
fig, ax = plt.subplots()
ax.plot(cumDisp, cumArea, label = 'Hysteresis Analysis')
ax.plot(Energy[:,0],Energy[:,1], label = 'Chan and Albermani')
ax.set_xlabel('Cumulative Deformation (mm)')
ax.set_ylabel('Cumulative Energy (kJ)')
plt.minorticks_on()
ax.grid(which='major', color='grey', linewidth=0.5, alpha = 0.8)
ax.grid(which='minor', linewidth=0.5, alpha = 0.4)
ax.legend()