\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/e/ec/Example1a_Push.GIF\")
\n"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Done!\n"
]
}
],
"source": [
"############################################################\n",
"# EXAMPLE: \n",
"# pyEx1a.Canti2D.Push.tcl.py\n",
"# for OpenSeesPy\n",
"# --------------------------------------------------------#\n",
"# by: Silvia Mazzoni, 2020\n",
"# silviamazzoni@yahoo.com\n",
"############################################################\n",
"# This file was obtained from a conversion of the updated Tcl script\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. cantilever 2D\n",
"# static pushover analysis 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",
"# 2 __\n",
"# or |\n",
"# or |\n",
"# or |\n",
"# (1) 36'\n",
"# or |\n",
"# or |\n",
"# or |\n",
"# =1= ---- -------->X\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",
"# define GEOMETRY -------------------------------------------------------------\n",
"# nodal coordinates:\n",
"ops.node(1,0,0) # node , X Y\n",
"ops.node(2,0,432)\n",
"\n",
"# Single point constraints -- Boundary Conditions\n",
"ops.fix(1,1,1,1) # node DX DY RZ\n",
"\n",
"# nodal masses:\n",
"ops.mass(2,5.18,0.,0.) # node , Mx My Mz, Mass=Weight/g.\n",
"\n",
"# Define ELEMENTS -------------------------------------------------------------\n",
"# define geometric transformation: performs a linear geometric transformation of beam stiffness\n",
"# 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",
"# element elasticBeamColumn eleTag iNode jNode A E Iz transfTag\n",
"ops.element('elasticBeamColumn',1,1,2,3600000000,4227,1080000,1) # element elasticBeamColumn 1 1 2 3600000000 4227 1080000 1;\n",
"\n",
"# Define RECORDERS -------------------------------------------------------------\n",
"ops.recorder('Node','-file','Data/DFreeEx1aPush.out','-time','-node',2,'-dof',1,2,3,'disp') # displacements of free nodes\n",
"ops.recorder('Node','-file','Data/DBaseEx1aPush.out','-time','-node',1,'-dof',1,2,3,'disp') # displacements of support nodes\n",
"ops.recorder('Node','-file','Data/RBaseEx1aPush.out','-time','-node',1,'-dof',1,2,3,'reaction') # support reaction\n",
"ops.recorder('Element','-file','Data/FColEx1aPush.out','-time','-ele',1,'globalForce') # element forces -- column\n",
"ops.recorder('Element','-file','Data/DColEx1aPush.out','-time','-ele',1,'deformation') # element deformations -- column\n",
"\n",
"# define GRAVITY -------------------------------------------------------------\n",
"ops.timeSeries('Linear',1) # timeSeries Linear 1;\n",
"# define Load Pattern\n",
"ops.pattern('Plain',1,1) # \n",
"ops.load(2,0.,-2000.,0.) # node , FX FY MZ -- superstructure-weight\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",
"# define LATERAL load -------------------------------------------------------------\n",
"# Lateral load pattern\n",
"ops.timeSeries('Linear',2) # timeSeries Linear 2;\n",
"# define Load Pattern\n",
"ops.pattern('Plain',2,2) # \n",
"ops.load(2,2000.,0.0,0.0) # node , FX FY MZ -- representative lateral load at top node\n",
"\n",
"\n",
"# pushover: diplacement controlled static analysis\n",
"ops.integrator('DisplacementControl',2,1,0.1) # switch to displacement control, for node 11, dof 1, 0.1 increment\n",
"ops.analyze(1000) # apply 100 steps of pushover analysis to a displacement of 10\n",
"\n",
"print('Done!')\n"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEICAYAAACktLTqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAARh0lEQVR4nO3df4xdZ33n8fcnNoRaJSWux1k3djwp9R84qQrUSul2RVFTIA2IRFBap9PiNNm1KkW7oXRVJeuVqko7Km2XbbSrompo0JriElkLhQillNRtqKg2ZCeQ0ObXxhD/Skw8lK4aEsnF4bt/zDG9ce7M3PlxfceP3y/p6pzzPM859/voSh8fn3vunFQVkqS2XDDqAiRJK89wl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuLUKS8SSVZO0AY29M8sWzUZd0JsNdTUtyKMk/J9lwRvtDXUiPj6YyabgMd50PngJuOL2R5EeB7xtdOdLwGe46H/wJ8L6e7V3Ax05vJPmBJB9LMpPkcJL/nOSCrm9Nkv+a5JtJvg68o/fA3b53Jjme5Okk/yXJmrMxKWk+hrvOB/cDFyV5XRe8vwh8vKf/fwA/APww8NPM/kPwq13fvwPeCbwB2AH8/BnH3gucAn6kG/M24N8OZxrS4Ax3nS9On72/FXgceLprPx32t1fVc1V1CPgQ8Ctd/y8Ad1TV0ar6FvA7pw+Y5BLg54D3V9XzVXUC+ANg51mYjzSvBb/xlxrxJ8DfAJfTc0kG2AC8Ejjc03YYuLRb/yHg6Bl9p20FXgEcT3K67YIzxksjYbjrvFBVh5M8BVwL3NzT9U3gO8wG9aNd22X8y5n9cWBLz/jLetaPAieBDVV1ahh1S0vlZRmdT24Gfqaqnu9pexHYD0wmeXWSrcAH+Jdr8vuB/5Bkc5KLgdtO71hVx4HPAx9KclGSC5K8NslPn5XZSPMw3HXeqKqvVdV0n65/DzwPfB34IvCnwEe7vo8AfwE8DHwZ+NQZ+76P2cs6jwL/CPwvYNOKFy8tUnxYhyS1xzN3SWqQ4S5JDTLcJalBhrskNWhV3Oe+YcOGGh8fH3UZknROefDBB79ZVWP9+lZFuI+PjzM93e8ONWkVOHoU3vc++MY34IILYPduuPXWUVclkeTwXH2rItylVW3tWvjQh+CNb4TnnoMf/3F461th+/ZRVybNyWvu0kI2bZoNdoBXvxpe9zp4+un595FGzHCXFuPQIfjKV+AnfmLUlUjzMtylQX372/Ce98Add8BFF426Gmlehrs0iO98ZzbYJybg3e8edTXSggx3aQ779u1jfHycCxI+efHFPJ7ABz4w6rKkgRjuUh/79u1j9+7dHD58mH8NvOf55zn1+c/zra1b4fWvh3vuGXWJ0ry8FVLqY8+ePbzwwgsA/C0QgCq2Jhx66KFRliYNxDN3qY8jR44sql1abQx3qY/LLrtsUe3SamO4S31MTk6ybt26l7StW7eOycnJEVUkLY7hLvUxMTHB1NQUF154IQBbt25lamqKiYmJEVcmDcYvVKU5TExM8JGPfASA++67b7TFSIvkmbskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0aONyTrEnylSSf7bbXJ7k3yZPd8uKesbcnOZjkiSRvH0bhkqS5LebM/VbgsZ7t24ADVbUNONBtk2Q7sBO4ArgG+HCSNStTriRpEAOFe5LNwDuAP+5pvg7Y263vBa7vab+rqk5W1VPAQeCqlSlXkjSIQc/c7wB+E/huT9slVXUcoFtu7NovBY72jDvWtb1Ekt1JppNMz8zMLLpwSdLcFgz3JO8ETlTVgwMeM33a6mUNVVNVtaOqdoyNjQ14aEnSIAb5w2E/BbwrybXAq4CLknwceDbJpqo6nmQTcKIbfwzY0rP/ZuCZlSxakjS/Bc/cq+r2qtpcVePMflH6V1X1y8DdwK5u2C7gM9363cDOJBcmuRzYBjyw4pVLkua0nD/5+0Fgf5KbgSPAewGq6pEk+4FHgVPALVX14rIrlSQNbFHhXlX3Afd16/8AXD3HuEnAR9ZI0oj4C1VJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNWjBcE/yqiQPJHk4ySNJfrtrX5/k3iRPdsuLe/a5PcnBJE8kefswJyBJerlBztxPAj9TVT8GvB64JsmbgNuAA1W1DTjQbZNkO7ATuAK4BvhwkjXDKF6S1N+C4V6zvt1tvqJ7FXAdsLdr3wtc361fB9xVVSer6ingIHDVilYtSZrXQNfck6xJ8hBwAri3qr4EXFJVxwG65cZu+KXA0Z7dj3VtkqSzZKBwr6oXq+r1wGbgqiRXzjM8/Q7xskHJ7iTTSaZnZmYGq1aSNJBF3S1TVf8PuI/Za+nPJtkE0C1PdMOOAVt6dtsMPNPnWFNVtaOqdoyNjS2hdEnSXAa5W2YsyWu69e8DfhZ4HLgb2NUN2wV8plu/G9iZ5MIklwPbgAdWunBJ0tzWDjBmE7C3u+PlAmB/VX02yf8G9ie5GTgCvBegqh5Jsh94FDgF3FJVLw6nfElSPwuGe1V9FXhDn/Z/AK6eY59JYHLZ1UmSlsRfqEpSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEsDuOkm2LgRrrxy1JVIgzHcpQHceCN87nOjrkIanOEuDeDNb4b160ddhTQ4w12SGmS4S1KDDHdJapDhLkkNMtylAdxwA/zkT8ITT8DmzXDnnaOuSJrf2lEXIJ0LPvGJUVcgLY5n7pLUIMNdkhpkuEtSgwx3SWrQguGeZEuSv07yWJJHktzata9Pcm+SJ7vlxT373J7kYJInkrx9mBOQJL3cIGfup4DfqKrXAW8CbkmyHbgNOFBV24AD3TZd307gCuAa4MNJ1gyjeElSfwuGe1Udr6ovd+vPAY8BlwLXAXu7YXuB67v164C7qupkVT0FHASuWunCJUlzW9Q19yTjwBuALwGXVNVxmP0HANjYDbsUONqz27Gu7cxj7U4ynWR6ZmZm8ZVLkuY0cLgn+X7gk8D7q+qf5hvap61e1lA1VVU7qmrH2NjYoGVIkgYwULgneQWzwb6vqj7VNT+bZFPXvwk40bUfA7b07L4ZeGZlypUkDWKQu2UC3Ak8VlX/rafrbmBXt74L+ExP+84kFya5HNgGPLByJUuSFjLI35b5KeBXgL9L8lDX9p+ADwL7k9wMHAHeC1BVjyTZDzzK7J02t1TViyteuSRpTguGe1V9kf7X0QGunmOfSWByGXVJkpbBX6hKUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1aMFwT/LRJCeS/H1P2/ok9yZ5slte3NN3e5KDSZ5I8vZhFS5JmtsgZ+7/E7jmjLbbgANVtQ040G2TZDuwE7ii2+fDSdasWLWSpIEsGO5V9TfAt85ovg7Y263vBa7vab+rqk5W1VPAQeCqFapVkjSgpV5zv6SqjgN0y41d+6XA0Z5xx7q2l0myO8l0kumZmZklliFJ6melv1BNn7bqN7CqpqpqR1XtGBsbW+EyJOn8ttRwfzbJJoBueaJrPwZs6Rm3GXhm6eVJkpZiqeF+N7CrW98FfKanfWeSC5NcDmwDHlheiZKkxVq70IAknwDeAmxIcgz4LeCDwP4kNwNHgPcCVNUjSfYDjwKngFuq6sUh1S5JmsOC4V5VN8zRdfUc4yeByeUUJUlaHn+hKkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3aRA33QQbN8KVV466Emkghrs0iBtvhM99btRVSAMz3KVBvPnNsH79qKuQBma4S1KDDHdJapDhLkkNGlq4J7kmyRNJDia5bVjvIw3Lvn37uP/++/nCF77A+Pg4n/70p0ddkjSwoYR7kjXAHwI/B2wHbkiyfRjvJQ3Dvn372L17NydPngTgdw4f5k2//ut89/HHYfNmuPPOEVcozW/tkI57FXCwqr4OkOQu4Drg0SG9n7Si9uzZwwsvvPC97V/qlls3b+bQoUMjqUlajGFdlrkUONqzfaxr+54ku5NMJ5memZkZUhnS0hw5cmRR7dJqM6xwT5+2eslG1VRV7aiqHWNjY0MqQ1qayy67bFHt0mozrHA/Bmzp2d4MPDOk95JW3OTkJOvWrXtJ27p165icnBxRRdLiDCvc/w+wLcnlSV4J7ATuHtJ7SStuYmKCqakptm7dShK2bt3K1NQUExMToy5NGkiqauFRSzlwci1wB7AG+GhVzXnKs2PHjpqenh5KHZLUqiQPVtWOfn3DuluGqroHuGdYx5ckzc1fqEpSgwx3SWqQ4S5JDTLcJalBQ7tbZlFFJDPA4VHXsQQbgG+OuoizzDmfH863OZ+r891aVX1/Bboqwv1clWR6rtuQWuWczw/n25xbnK+XZSSpQYa7JDXIcF+eqVEXMALO+fxwvs25ufl6zV2SGuSZuyQ1yHCXpAYZ7gtIsj7JvUme7JYXzzFu3geCJ/mPSSrJhuFXvTzLnXOS30/yeJKvJvmzJK85e9UPboDPLEn+e9f/1SRvHHTf1Wqpc06yJclfJ3ksySNJbj371S/Ncj7nrn9Nkq8k+ezZq3oFVJWveV7A7wG3deu3Ab/bZ8wa4GvADwOvBB4Gtvf0bwH+gtkfam0Y9ZyGPWfgbcDabv13++0/6tdCn1k35lrgz5l9stibgC8Nuu9qfC1zzpuAN3brrwb+b+tz7un/APCnwGdHPZ/FvDxzX9h1wN5ufS9wfZ8x33sgeFX9M3D6geCn/QHwm5zxqMFVbFlzrqrPV9Wpbtz9zD6Ja7VZ6DOj2/5YzbofeE2STQPuuxotec5VdbyqvgxQVc8Bj3HGc5FXqeV8ziTZDLwD+OOzWfRKMNwXdklVHQfolhv7jJnzgeBJ3gU8XVUPD7vQFbSsOZ/hJmbPilabQeqfa8ygc19tljPn70kyDrwB+NKKV7jyljvnO5g9MfvusAoclqE9rONckuQvgX/Vp2vPoIfo01ZJ1nXHeNtSaxuWYc35jPfYA5wC9i2uurNiwfrnGTPIvqvRcuY825l8P/BJ4P1V9U8rWNuwLHnOSd4JnKiqB5O8ZcUrGzLDHaiqn52rL8mzp/9b2v1X7USfYXM9EPy1wOXAw0lOt385yVVV9Y0Vm8ASDHHOp4+xC3gncHV1Fy5XmUEe4j7XmFcOsO9qtJw5k+QVzAb7vqr61BDrXEnLmfPPA+/qHhn6KuCiJB+vql8eYr0rZ9QX/Vf7C/h9Xvrl4u/1GbMW+DqzQX76S5sr+ow7xLnxheqy5gxcAzwKjI16LvPMccHPjNlrrb1ftD2wmM97tb2WOecAHwPuGPU8ztaczxjzFs6xL1RHXsBqfwE/CBwAnuyW67v2HwLu6Rl3LbN3EHwN2DPHsc6VcF/WnIGDzF7DfKh7/dGo5zTHPF9WP/BrwK916wH+sOv/O2DHYj7v1fha6pyBf8Ps5Yyv9nyu1456PsP+nHuOcc6Fu39+QJIa5N0yktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ16P8DfIo8+L/ri9UAAAAASUVORK5CYII=\n",
"text/plain": [
"