\n",
"\n",
"You can find the original Examples:
\n",
"https://opensees.berkeley.edu/wiki/index.php/Examples_Manual
\n",
"Original Examples by By Silvia Mazzoni & Frank McKenna, 2006, in Tcl
\n",
"Converted to OpenSeesPy by SilviaMazzoni, 2020
\n",
"
\n", "\n", "
\n", " \n", "Introductory Examples\n", "The objective of Example 1a and Example 1b is to give an overview of input-file format in OpenSees using simple scripts.\n", "These scripts do not take advantage of the Tcl scripting capabilities shown in the later examples. However, they do provide starting a place where the input file is similar to that of more familiar Finite-Element Analysis software. Subsequent examples should be used as the basis for user input files.\n", " " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
![](\"https://opensees.berkeley.edu/wiki/images/d/d5/Example1b_EQ.GIF\")
\n",
"\n",
" "
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Done!\n"
]
}
],
"source": [
"############################################################\n",
"# EXAMPLE: \n",
"# pyEx1b.Portal2D.EQ.tcl.py\n",
"# for OpenSeesPy\n",
"# --------------------------------------------------------#\n",
"# by: Silvia Mazzoni, 2020\n",
"# silviamazzoni@yahoo.com\n",
"############################################################\n",
"\n",
"# configure Python workspace\n",
"import openseespy.opensees as ops\n",
"import eSEESminiPy\n",
"import os\n",
"import math\n",
"import numpy as numpy\n",
"import matplotlib.pyplot as plt\n",
"ops.wipe()\n",
"# --------------------------------------------------------------------------------------------------\n",
"# Example 1. portal frame in 2D\n",
"# dynamic earthquake analysis of Portal Frame, with gravity.\n",
"# all units are in kip, inch, second\n",
"# elasticBeamColumn ELEMENT\n",
"# Silvia Mazzoni and Frank McKenna, 2006\n",
"#\n",
"# ^Y\n",
"# or\n",
"# 3_________(3)________4 __\n",
"# or | |\n",
"# or | |\n",
"# or | |\n",
"# (1) (2) LCol\n",
"# or | |\n",
"# or | |\n",
"# or | |\n",
"# =1= =2= _or_ -------->X\n",
"# or----------LBeam------------|\n",
"#\n",
"\n",
"# SET UP ----------------------------------------------------------------------------\n",
"ops.wipe() # clear opensees model\n",
"ops.model('basic','-ndm',2,'-ndf',3) # 2 dimensions, 3 dof per node\n",
"if not os.path.exists('Data'):\n",
" os.mkdir('Data')\n",
"\n",
"\n",
"# define GEOMETRY -------------------------------------------------------------\n",
"# nodal coordinates:\n",
"ops.node(1,0,0) # node , X Y\n",
"ops.node(2,504,0)\n",
"ops.node(3,0,432)\n",
"ops.node(4,504,432)\n",
"\n",
"# Single point constraints -- Boundary Conditions\n",
"ops.fix(1,1,1,1) # node DX DY RZ\n",
"ops.fix(2,1,1,1) # node DX DY RZ\n",
"ops.fix(3,0,0,0)\n",
"ops.fix(4,0,0,0)\n",
"\n",
"# nodal masses:\n",
"ops.mass(3,5.18,0.,0.) # node , Mx My Mz, Mass=Weight/g.\n",
"ops.mass(4,5.18,0.,0.)\n",
"\n",
"# Define ELEMENTS -------------------------------------------------------------\n",
"# define geometric transformation: performs a linear geometric transformation of beam stiffness and resisting force from the basic system to the global-coordinate system\n",
"ops.geomTransf('Linear',1) # associate a tag to transformation\n",
"\n",
"# connectivity: (make A very large, 10e6 times its actual value)\n",
"ops.element('elasticBeamColumn',1,1,3,3600000000,4227,1080000,1) # element elasticBeamColumn eleTag iNode jNode A E Iz transfTag\n",
"ops.element('elasticBeamColumn',2,2,4,3600000000,4227,1080000,1)\n",
"ops.element('elasticBeamColumn',3,3,4,5760000000,4227,4423680,1)\n",
"\n",
"# Define RECORDERS -------------------------------------------------------------\n",
"ops.recorder('Node','-file','Data/DFreeEx1bEQ.out','-time','-node',3,4,'-dof',1,2,3,'disp') # displacements of free nodes\n",
"ops.recorder('Node','-file','Data/DBaseEx1bEQ.out','-time','-node',1,2,'-dof',1,2,3,'disp') # displacements of support nodes\n",
"ops.recorder('Node','-file','Data/RBaseEx1bEQ.out','-time','-node',1,2,'-dof',1,2,3,'reaction') # support reaction\n",
"ops.recorder('Element','-file','Data/FColEx1bEQ.out','-time','-ele',1,2,'globalForce') # element forces -- column\n",
"ops.recorder('Element','-file','Data/FBeamEx1bEQ.out','-time','-ele',3,'globalForce') # element forces -- beam\n",
"\n",
"# define GRAVITY -------------------------------------------------------------\n",
"ops.timeSeries('Linear',1) # timeSeries Linear 1;\n",
"# define Load Pattern\n",
"ops.pattern('Plain',1,1) # \n",
"ops.eleLoad('-ele',3,'-type','-beamUniform',-7.94) # distributed superstructure-weight on beam\n",
"\n",
"ops.wipeAnalysis() # adding this to clear Analysis module \n",
"ops.constraints('Plain') # how it handles boundary conditions\n",
"ops.numberer('Plain') # renumber dofs to minimize band-width (optimization), if you want to\n",
"ops.system('BandGeneral') # how to store and solve the system of equations in the analysis\n",
"ops.test('NormDispIncr',1.0e-8,6) # determine if convergence has been achieved at the end of an iteration step\n",
"ops.algorithm('Newton') # use Newtons solution algorithm: updates tangent stiffness at every iteration\n",
"ops.integrator('LoadControl',0.1) # determine the next time step for an analysis, apply gravity in 10 steps\n",
"ops.analysis('Static') # define type of analysis static or transient\n",
"ops.analyze(10) # perform gravity analysis\n",
"ops.loadConst('-time',0.0) # hold gravity constant and restart time\n",
"\n",
"# DYNAMIC ground-motion analysis -------------------------------------------------------------\n",
"# create load pattern\n",
"accelSeries = 900\n",
"ops.timeSeries('Path',accelSeries,'-dt',0.01,'-filePath','BM68elc.acc','-factor',1) # define acceleration vector from file (dt=0.01 is associated with the input file gm)\n",
"ops.pattern('UniformExcitation',2,1,'-accel',accelSeries) # define where and how (pattern tag, dof) acceleration is applied\n",
"ops.rayleigh(0.,0.,0.,2*0.02/math.pow(ops.eigen('-fullGenLapack',1)[0],0.5)) # set damping based on first eigen mode\n",
"\n",
"# create the analysis\n",
"ops.wipeAnalysis() # clear previously-define analysis parameters\n",
"ops.wipeAnalysis() # adding this to clear Analysis module \n",
"ops.constraints('Plain') # how it handles boundary conditions\n",
"ops.numberer('Plain') # renumber dofs to minimize band-width (optimization), if you want to\n",
"ops.system('BandGeneral') # how to store and solve the system of equations in the analysis\n",
"ops.test('NormDispIncr',1.0e-8,10) # determine if convergence has been achieved at the end of an iteration step\n",
"ops.algorithm('Newton') # use Newtons solution algorithm: updates tangent stiffness at every iteration\n",
"ops.integrator('Newmark',0.5,0.25) # determine the next time step for an analysis\n",
"ops.analysis('Transient') # define type of analysis: time-dependent\n",
"ops.analyze(1000,0.02) # apply 1000 0.02-sec time steps in analysis\n",
"\n",
"\n",
"print('Done!')\n",
"\n",
"\n",
"\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEICAYAAACktLTqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAATaUlEQVR4nO3dfYxd9X3n8feHMSFhgQYXG7kYbLK1VB7aEjICqiBCSFMeggrK05p1WtPQ9R9GuyGsVMGy2lWktcpGSxe0AiRviNYNNJa1pcKiiJZCSQRawppgWrBhcYLBBhc7TaOEpHVt890/7jG5mBnPjOfOjOfn90u6uuf8ztP3O7r+zPG5Z+5NVSFJastRM12AJGnwDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7tIEJFmcpJLMGce61yZ5Yjrqkg5kuKtpSbYm+eckJx0wvrEL6cUzU5k0tQx3HQleAa7ZP5PkV4EPzFw50tQz3HUk+Abwu33zy4E/3j+T5BeS/HGSXUleTfIfkxzVLRtK8t+S/CDJ94FP9e+42/aeJDuSvJ7kvyQZmo6mpIMx3HUkeAo4IckZXfD+K+DevuX/A/gF4EPAx+j9Ivi9btm/Aa4EPgwMA589YN9rgL3AL3fr/Bbw+1PThjR+hruOFPvP3j8JvAi83o3vD/ubq+onVbUVuA34nW7554Hbq2pbVf0Q+MP9O0xyMnA5cENV/bSqdgL/HVg6Df1IBzXmO/5SI74BfBs4nb5LMsBJwPuAV/vGXgVO6aZ/Cdh2wLL9FgFHAzuS7B876oD1pRlhuOuIUFWvJnkFuAK4rm/RD4A99IJ6Uzd2Gj8/s98BnNq3/ml909uA3cBJVbV3KuqWDpWXZXQkuQ64pKp+2je2D1gHrEpyfJJFwI38/Jr8OuDfJVmY5ETgpv0bVtUO4C+B25KckOSoJP8yycempRvpIAx3HTGq6ntVtWGERf8W+CnwfeAJ4E+Ar3fL/ifwF8BzwHeB+w/Y9nfpXdbZBPwD8L+BBQMvXpqg+GUdktQez9wlqUGGuyQ1yHCXpAYZ7pLUoMPiPveTTjqpFi9ePNNlSNKs8swzz/ygquaNtOywCPfFixezYcNId6iN0z/9E1x0EezeDXv3wmc/C1/5yuAK1BHNl5emxb59MDwMp5wCDz44rk2SvDrassMi3CftmGPgscfguONgzx648EK4/HK44IKZrkwN8OWlaXHHHXDGGfDjHw9kd21cc096//Kg969vz57emDQAvrw05bZvhz//c/j9wX2gaBvhDr3/0pxzDsyfD5/8JJx//kxXpIb48tKUuuEG+OpX4ajBRXI74T40BBs39n4DPv00PP/8TFekhvjy0pR58MHeWcNHPjLQ3bYT7vt98INw8cXw8MMzXYka5MtLA/fkk7B+PSxeDEuX9t7g+cIXJr3bWR3uK1euZM6cOcxL+MWhIVauXAn/+I/wV38Fv/IrM12eGrFrF/zoR71pX14alP35lVtvZc7f/R0rr7gC1q6FSy6Be+8dewdjmLV3y6xcuZK7774b6H0E35q332bo7rt5Y+1afumGG+DKK2e2QDVjxw5Yvrx33f3tt+Hzn/flpcnpzy+Affv2cffdd7Pk9df58oCOcVh8KuTw8HBN9D73OXPmsG/fvhGXfexjfpy2pMPXt771rRHHh4aG2Lt3/N/7kuSZqhoeadmsvSwzWrBL0mw1yFybtZdlhoaGRvxBDA0N8fjjj09/QZI0TqNdeRgaGhrYMWbtmfuKFSsmNC5Jh4vpyK9Ze+Z+1113AbzzpsTQ0BArVqx4Z1ySDlfTkV+z9g3V/S6++GIAL8VImnUmm19NvqEqSRqd4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaNO5wTzKU5NkkD3bzc5M8kuTl7vnEvnVvTrIlyUtJLp2KwiVJo5vImfuXgM198zcBj1bVEuDRbp4kZwJLgbOAy4C7kgzuAxMkSWMaV7gnWQh8Cvha3/BVwJpueg1wdd/42qraXVWvAFuA8wZTriRpPMZ75n478AfA231jJ1fVDoDueX43fgqwrW+97d3YuyRZkWRDkg27du2acOGSpNGNGe5JrgR2VtUz49xnRhh7zwfYVNXqqhququF58+aNc9eSpPEYz6dCfhT47SRXAO8HTkhyL/BmkgVVtSPJAmBnt/524NS+7RcCbwyyaEnSwY155l5VN1fVwqpaTO+N0seq6gvAemB5t9py4IFuej2wNMkxSU4HlgBPD7xySdKoJvN57rcC65JcB7wGfA6gql5Isg7YBOwFrq8qvxNPkqbRhMK9qh4HHu+m/x74xCjrrQJWTbI2SdIh8i9UJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktSgMcM9yfuTPJ3kuSQvJPlKNz43ySNJXu6eT+zb5uYkW5K8lOTSqWxAkvRe4zlz3w1cUlW/DpwDXJbkAuAm4NGqWgI82s2T5ExgKXAWcBlwV5KhqShekjSyMcO9et7qZo/uHgVcBazpxtcAV3fTVwFrq2p3Vb0CbAHOG2jVkqSDGtc19yRDSTYCO4FHquo7wMlVtQOge57frX4KsK1v8+3dmCRpmowr3KtqX1WdAywEzkty9kFWz0i7eM9KyYokG5Js2LVr1/iqlSSNy4TulqmqHwGP07uW/maSBQDd885ute3AqX2bLQTeGGFfq6tquKqG582bdwilS5JGM567ZeYl+WA3/QHgN4EXgfXA8m615cAD3fR6YGmSY5KcDiwBnh504ZKk0c0ZxzoLgDXdHS9HAeuq6sEk/wdYl+Q64DXgcwBV9UKSdcAmYC9wfVXtm5ryJUkjGTPcq+pvgA+PMP73wCdG2WYVsGrS1UmSDol/oSpJDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWpQM+H+xS/C/Plw9tkzXYkkTcy2bfDxj8MZZ8BZZ8Edd0x+n82E+7XXwsMPz3QVkjRxc+bAbbfB5s3w1FNw552wadPk9tlMuF90EcydO9NVSNLELVgA557bmz7++N4Z/OuvT26fzYS7JLVg61Z49lk4//zJ7cdwl6TDxFtvwWc+A7ffDiecMLl9Ge6SdBjYs6cX7MuWwac/Pfn9Ge6SNMOq4Lrretfab7xxMPtsJtyvuQZ+4zfgpZdg4UK4556ZrkiSxufJJ+Eb34DHHoNzzuk9HnpocvucM5jSZt43vznTFUjSobnwwt7Z+yA1c+YuSfo5w12SGmS4S1KDDHdJatCY4Z7k1CR/nWRzkheSfKkbn5vkkSQvd88n9m1zc5ItSV5KculUNiBJeq/xnLnvBf59VZ0BXABcn+RM4Cbg0apaAjzazdMtWwqcBVwG3JVkaCqKlySNbMxwr6odVfXdbvonwGbgFOAqYE232hrg6m76KmBtVe2uqleALcB5gy5ckjS6CV1zT7IY+DDwHeDkqtoBvV8AwPxutVOAbX2bbe/GDtzXiiQbkmzYtWvXxCuXJI1q3OGe5DjgT4EbqurHB1t1hLH33J5fVaurariqhufNmzfeMiRJ4zCucE9yNL1gv6+q7u+G30yyoFu+ANjZjW8HTu3bfCHwxmDKlSSNx3julglwD7C5qv6ob9F6YHk3vRx4oG98aZJjkpwOLAGeHlzJkqSxjOezZT4K/A7wt0k2dmP/AbgVWJfkOuA14HMAVfVCknXAJnp32lxfVfsGXrkkaVRjhntVPcHI19EBPjHKNquAVZOoS5I0Cf6FqiQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoPGDPckX0+yM8nzfWNzkzyS5OXu+cS+ZTcn2ZLkpSSXTlXhkqTRjefM/X8Blx0wdhPwaFUtAR7t5klyJrAUOKvb5q4kQwOrVpI0LmOGe1V9G/jhAcNXAWu66TXA1X3ja6tqd1W9AmwBzhtQrZKkcTrUa+4nV9UOgO55fjd+CrCtb73t3dh7JFmRZEOSDbt27TrEMiRJIxn0G6oZYaxGWrGqVlfVcFUNz5s3b8BlSNKR7VDD/c0kCwC6553d+Hbg1L71FgJvHHp5kqRDcajhvh5Y3k0vBx7oG1+a5JgkpwNLgKcnV6IkaaLmjLVCkm8CFwMnJdkO/GfgVmBdkuuA14DPAVTVC0nWAZuAvcD1VbVvimqXJI1izHCvqmtGWfSJUdZfBayaTFGSpMnxL1QlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWpQO+H+xS/C/Plw9tkzXYkkTcy2bfDxj8MZZ8BZZ8Edd0x6l+2E+7XXwsMPz3QVkjRxc+bAbbfB5s3w1FNw552wadOkdtlOuF90EcydO9NVSNLELVgA557bmz7++N4Z/OuvT2qX7YS7JLVg61Z49lk4//xJ7cZwl6TDxVtvwWc+A7ffDiecMKldGe6SdDjYs6cX7MuWwac/PendzRlASSNKchlwBzAEfK2qbh30Me677z6eeuopdu/ezeLFi7n9hhu4etAHkaQp8K78WrSIx047jQ995CNw440D2f+UnLknGQLuBC4HzgSuSXLmII9x3333sWLFCnbv3g3AH776Khd8+cu8/eKLsHAh3HPPIA8nSQNzYH4tfO01PvTEE/zD/ffDOef0Hg89NKljTNWZ+3nAlqr6PkCStcBVwOTu7elzyy238LOf/eyd+X/dPS9auJCtW7cO6jCSNHAH5teTQIBFRx3F1o0bB3KMqbrmfgqwrW9+ezf2jiQrkmxIsmHXrl0TPsBrr702oXFJOlxMR35NVbhnhLF610zV6qoarqrhefPmTfgAp5122oTGJelwMR35NVXhvh04tW9+IfDGIA+watUqjj322HeNHXvssaxatWqQh5GkgZuO/JqqcP+/wJIkpyd5H7AUWD/IAyxbtozVq1ezaNEikrBo0SJWr17NsmXLBnkYSRq46civVNXYax3KjpMrgNvp3Qr59aoa9VfS8PBwbdiwYUrqkKRWJXmmqoZHWjZl97lX1UPA5O7lkSQdEv9CVZIaZLhLUoMMd0lqkOEuSQ2asrtlJlREsgt4dRK7OAn4wYDKmQ3st31HWs/2e2gWVdWIfwV6WIT7ZCXZMNrtQC2y3/YdaT3b7+B5WUaSGmS4S1KDWgn31TNdwDSz3/YdaT3b74A1cc1dkvRurZy5S5L6GO6S1KBZHe5JLkvyUpItSW6a6XoGIcnXk+xM8nzf2NwkjyR5uXs+sW/ZzV3/LyW5dGaqnpwkpyb56ySbk7yQ5EvdeJN9J3l/kqeTPNf1+5VuvMl+90sylOTZJA928633uzXJ3ybZmGRDNzZ9PVfVrHzQ+yjh7wEfAt4HPAecOdN1DaCvi4Bzgef7xr4K3NRN3wT81276zK7vY4DTu5/H0Ez3cAg9LwDO7aaPB/5f11uTfdP7prLjuumjge8AF7Tab1/fNwJ/AjzYzbfe71bgpAPGpq3n2Xzm/s6XcFfVPwP7v4R7VquqbwM/PGD4KmBNN70GuLpvfG1V7a6qV4At9H4us0pV7aiq73bTPwE20/vO3Sb7rp63utmju0fRaL8ASRYCnwK+1jfcbL8HMW09z+ZwH/NLuBtyclXtgF4QAvO78eZ+BkkWAx+mdzbbbN/dJYqNwE7gkapqul96X9zzB8DbfWMt9wu9X9h/meSZJCu6sWnrecq+rGMajPkl3EeApn4GSY4D/hS4oap+nIzUXm/VEcZmVd9VtQ84J8kHgT9LcvZBVp/V/Sa5EthZVc8kuXg8m4wwNmv67fPRqnojyXzgkSQvHmTdgfc8m8/cp/xLuA8jbyZZANA97+zGm/kZJDmaXrDfV1X3d8PN911VPwIeBy6j3X4/Cvx2kq30Lp9ekuRe2u0XgKp6o3veCfwZvcss09bzbA73Kf8S7sPIemB5N70ceKBvfGmSY5KcDiwBnp6B+iYlvVP0e4DNVfVHfYua7DvJvO6MnSQfAH4TeJFG+62qm6tqYVUtpvfv9LGq+gKN9guQ5F8kOX7/NPBbwPNMZ88z/Y7yJN+NvoLenRXfA26Z6XoG1NM3gR3AHnq/za8DfhF4FHi5e57bt/4tXf8vAZfPdP2H2POF9P4L+jfAxu5xRat9A78GPNv1+zzwn7rxJvs9oPeL+fndMs32S+8uvue6xwv782k6e/bjBySpQbP5sowkaRSGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWrQ/weCTboaHWIhiAAAAABJRU5ErkJggg==\n",
"text/plain": [
"