{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Maximum Dynamic Pressure" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "1. https://en.wikipedia.org/wiki/Max_Q\n", "2. https://en.wikipedia.org/wiki/Density_of_air" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 75, "metadata": {}, "outputs": [], "source": [ "def get_ro(h):\n", " p0 = 101.325e3 # kPa\n", " T0 = 288.15 # K\n", " g = 9.806 # m/s\n", " L = 0.0065 # K/m\n", " R = 8.31447 # J/(mol K)\n", " M = 0.0289644 # kg/mol\n", " \n", " T = T0 - L * h \n", " p = p0 * (1 - L * h / T0) ** (g * M / (R * L))\n", " return p * M / (R * T)" ] }, { "cell_type": "code", "execution_count": 76, "metadata": {}, "outputs": [], "source": [ "# vertical acceleration\n", "ay = 3 # m/s2" ] }, { "cell_type": "code", "execution_count": 79, "metadata": {}, "outputs": [], "source": [ "t = np.arange(0, 90, 1)\n", "vy = ay * t\n", "h = vy * t\n", "ro = get_ro(h)\n", "q = 1 / 2 * (ro * vy) ** 2" ] }, { "cell_type": "code", "execution_count": 80, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEWCAYAAABxMXBSAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd8VFX6x/HPk0ZCEnpCCxBK6J1IEVREULB3UVFEXazr\nou666hZdy09d26qruKgodlkrKgsiAoLUIL0mkAAJJQmhJZD+/P6Yi0YFQiCTO+V5v17zytwzd2a+\nCUOenHvPPUdUFWOMMeZ4hbgdwBhjjH+xwmGMMaZKrHAYY4ypEiscxhhjqsQKhzHGmCqxwmGMMaZK\nrHAYY4ypEiscxhhjqsQKhzHHSUQyRKRYRBr9qn2ZiKiIJHrhPWuJyBMislVEDolIqoj8UUSkut/L\nmONlhcOYqkkHrj68ISLdgNpefL//AmcB5wKxwHXALcCzXnxPY47JCocxVfMOcH2F7dHA2xV3EJHz\nnF7IfhHZJiIPV3jsKhFJF5E6zvYIEdkpInG/fiMROQs4G7hMVVeraqmqLgRGAX8QkTbV/+0ZUzkr\nHMZUzUKgjoh0EpFQYCTw7q/2KcBTXOoB5wG3icjFAKr6ETAfeFFEGgJvADeras4R3msYsEhVt1Vs\nVNVFQCaenogxNc4KhzFVd7jXMQxYB2RVfFBVZ6vqKlUtV9WVwAfAGRV2uQMYAswGvlTVr47yPo2A\nHUd5bAfwm16KMTUhzO0Axvihd4Dvgdb86jAVgIj0A54EugIRQC085yoAUNW9IvJf4B7gsmO8Ty6Q\ndJTHmjqPG1PjrMdhTBWp6hY8J8nPBT49wi7vA1OAFqpaF3gV+GkUlIj0BG7E0xN58Rhv9S3QT0Ra\nVGx0ClNLYM5JfBvGnDArHMacmJuAIapacITHYoE8VS0Ukb7ANYcfEJFIPOdEHgTGAM1F5PYjvYGq\nfgvMBD4RkS4iEioi/Z3nv62qG6r3WzLm+FjhMOYEqOomVU05ysO3A4+IyAHg78DkCo89AWxT1fGq\nWoRnhNRjInK0Q1KXAbOAaUAhsMC5P7Yavg1jTojYCoDG+A8RmQQ0A85T1WK385jgZD0OY/zLzXjO\nffR2O4gJXtbjMMYYUyXW4zDGGFMlAXkdR6NGjTQxMdHtGMYY41eWLl2aq6qVXlgakIUjMTGRlJSj\nDXgxxhhzJCKy5Xj2s0NVxhhjqsQKhzHGmCqxwmGMMaZKrHAYY4ypEiscxhhjqsQKhzHGmCqxwmGM\nMaZKAvI6DmNOVlm5sn3vITbl5LM5p4C9h0oIFSEsVKgVFkLbuBg6No2lSZ1IRKTyFzQmgFjhMMaR\nX1TKrPXZTFuzk9nrsykoLqv0OXWjwunXugHn92jGWR3jia5l/6VM4PP6p1xEQoEUIEtVzxeR1sCH\nQENgKXCdqhaLSC08y3D2AXYDV6lqhvMaD+BZOKcMuEtVp3s7twkem3PyeXnWJr5cuZ3i0nIaxURw\nYc/m9EioS5u4GNrERdMwOoKycqVMlUPFZaRm57N+x37W7tjPzHXZfLN2F5HhIQzr3IRbz2hDl2Z1\n3f62jPGamvjz6A/AOqCOs/0U8Lyqfigir+IpCOOdr3tUtZ2IjHT2u0pEOgMjgS541iH4VkTaq2rl\nfw4acwxp2Qd46bs0vlyxnYiwEK5KbsEFPZrRp1V9QkN+e/gpLFQIA2qFhXJKYgNOSWwAeA5rLcnI\n4+uVO/h8WRZfrtjO2Z0bc9dZSXRtbgXEBB6vTqsuIgnAJOBx4B7gAiAHaKKqpSIyAHhYVc8RkenO\n/QUiEgbsBOKA+wFU9QnnNX/a72jvm5ycrDZXlTmaotIyXpqZxqtzNhERFsJ1A1rxu9Pa0Cim1km/\n9r6DJUz8IZ2JP6RzoLCUS3s35+/nd6Ze7YhqSG6Md4nIUlVNrmw/b/c4/gXch2cNZvAcntqrqqXO\ndibQ3LnfHNgG4BSVfc7+zYGFFV6z4nN+IiJjcZbTbNmyZfV+FyZgLNu6h/s+Xklqdj6X9U7gL+d1\nokF09f1Sr1s7nLuHteem01rz6uxNTPh+M99vzOGRi7oyomsTO5FuAoLXhuOKyPlAtqou9dZ7VKSq\nE1Q1WVWT4+IqnRXYBBlV5ZXZaVw2fj75RaW8OeYUnr2yR7UWjYrqRIZz3/COTLlzEE3rRnH7ez9y\nx/s/cqCwxCvvZ0xN8maPYyBwoYicC0TiOcfxAlBPRMKcXkcCkOXsnwW0ADKdQ1V18ZwkP9x+WMXn\nGFOpwpIy7vt4JVNWbOeCHs14/JKu1IkMr5H37tysDp/dfiqvzU3nmW82sH7nD0y4Lpl28TE18v7G\neIPXehyq+oCqJqhqIp6T29+p6rXALOByZ7fRwBfO/SnONs7j36nnBMwUYKSI1HJGZCUBi72V2wSW\nnfsKueLVBXy5cjt/OqcDL47sWWNF47Cw0BBuG9yW927ux/5DJVz88g9MX7OzRjMYU53cuHL8z8A9\nIpKG5xzGG077G0BDp/0efj4pvgaYDKwFpgF32Igqczwycgu45JUf2JyTz4TrkrnjzHaunmPo36Yh\nX/5+EG3jY7jlnaVM+H6Ta1mMORleHVXlFhtVZTbn5HP1awspLi3n3Zv7+dR1FYUlZdz73xV8vXIH\ntw9uy5/O6WAnzY1P8JVRVcbUuLRsT9EoL1c+GNufjk3qVP6kGhQZHsqLI3tRJzKMV2ZvYn9hCY9c\n2JWQI1w7YowvssJhAkp6bgEjJ3hGb384tj9JjWMreYY7QkOE/7ukG3WiwvnPnM0cLCrjmSt6WPEw\nfsEKhwkYu/OLuOHNxZSrMvmWAT4/cklEeGBEJ6IjwnhuxkaiIkJ57OKudtjK+DwrHCYgHCou46ZJ\nKezcV8gHY/v7fNGo6PdD2lFQXMp/5mwmNjKc+0d0dDuSMcdkhcP4vbJyZdxHy1iRuZfx1/ahd8v6\nbkeqEhHh/uEdyS8s5dU5m4iNDOOOM9u5HcuYo7LCYfzek/9bx/Q1u3j4gs4M79rE7TgnRER49KKu\n5BeV8vT0DTSuE8nlfRLcjmXMEVnhMH7t65U7eG1uOqMHtOKGga3djnNSQkKEZ67oQc6BIh74dCUt\nG9Smb+sGbscy5jds6Vjjt9Ky87nv4xX0blmPv5zX2e041SI8NITx1/ahRf3a3PJOClt2F7gdyZjf\nsMJh/FJBUSm3vbuUWuGhvHxtbyLCAuejXLd2OG/ccArlCjdNSmHfIZsY0fiWwPnfZoKGqvLAp6vY\nlJPPS1f3omndKLcjVbvWjaJ5dVQfMnILGPfhMsrLA2+GB+O/rHAYv/PJj1lMWbGde4a1Z2C7Rm7H\n8ZoBbRvy9ws6M2tDDq/MTnM7jjE/scJh/Mq2vIM8PGUNfVs34LbBgT9k9br+rbioZzOem7GRH9Jy\n3Y5jDGCFw/iRsnLl3skrAHj2ih5HXBc80IgIT1zajbZxMdz1wTJ27it0O5IxVjiM/3h97mYWZ+Tx\n8IVdaNGgtttxakztiDDGj+pDYUkZd7z/IyVl5W5HMkHOCofxC+t27OfZbzYyvEsTLuv9myXnA167\n+BievKw7S7fs4aWZqW7HMUHOm2uOR4rIYhFZISJrROQfTvtbIpIuIsudW0+nXUTkRRFJE5GVItK7\nwmuNFpFU5zb6aO9pAlNpWTn3fbySOlFh/N+l3YJ2EsALejTj8j4J/HtWGovT89yOY4KYN3scRcAQ\nVe0B9ASGi0h/57E/qWpP57bcaRuBZ1nYJGAsMB5ARBoADwH9gL7AQyLiX5MRmZPy1vwMVmXt46EL\nutAgOsLtOK46fJju7o+W2/UdxjXeXHNcVTXf2Qx3bscajH4R8LbzvIVAPRFpCpwDzFDVPFXdA8wA\nhnsrt/Et2/IO8uw3GxnSMZ7zuzd1O47rYmqF8a+rerJzfyF/+3w1gbiCp/F9Xj3HISKhIrIcyMbz\ny3+R89DjzuGo50WkltPWHNhW4emZTtvR2n/9XmNFJEVEUnJycqr9ezE1T1X56+erEYFHbZ2Kn/Rq\nWZ+7hyYxZcV2vli+3e04Jgh5tXCoapmq9gQSgL4i0hV4AOgInAI0AP5cTe81QVWTVTU5Li6uOl7S\nuGzKiu3M2ZjDH8/uQPN6gXd1+Mm4bXA7+rSqz0NT1pC934bomppVI6OqVHUvMAsYrqo7nMNRRcCb\neM5bAGQBLSo8LcFpO1q7CWD7C0t49Ku1dE+oy+hTE92O43NCQ4SnL+9OYUkZD362yg5ZmRrlzVFV\ncSJSz7kfBQwD1jvnLRDPcYeLgdXOU6YA1zujq/oD+1R1BzAdOFtE6jsnxc922kwAe+HbVHYXFPP4\nxd2C4kK/E9EmLoY/ndOBb9dl89ky+1vK1BxvrsfRFJgkIqF4CtRkVf1KRL4TkThAgOXArc7+U4Fz\ngTTgIDAGQFXzRORRYImz3yOqamMRA1jqrgNMmp/ByFNa0C2hrttxfNqYga2ZtnonD09Zw8B2jWhc\nJ9LtSCYISCB2cZOTkzUlJcXtGOYEqCrXvbGYlZl7mfXHwTSMqVX5k4Lc5px8Rrwwl0HtGvH66GQb\nRGBOmIgsVdXkyvazK8eNT5m+Zhfz0nK5Z1h7KxrHqU1cDH88uwMz12czddVOt+OYIGCFw/iMwpIy\nHvt6Le0bxzCqfyu34/iVMQMT6da8Lg9NWcO+g3ZhoPEuKxzGZ7wxL53MPYd4+IIuhIXaR7MqwkJD\neOLSbuw5WMyT09a5HccEOPvfaXzC7vwixs/exNBOjTk1gBdn8qauzety06DWfLB4G4s273Y7jglg\nVjiMT3jpuzQOlZRx/4iObkfxa+OGJpFQP4oHPltFYUmZ23FMgLLCYVyXnlvAuwu3cNUpLWgXH+N2\nHL9WOyKMxy7uyuacAiZ8v9ntOCZAWeEwrnt6+noiwkIYNzTJ7SgBYXCHeM7r1pSXZ6WxLe+g23FM\nALLCYVy1dMsepq7aydjT2xAfaxevVZe/nt+J0BDh4Slr3I5iApAVDuMaVeWp/60nLrYWvzutjdtx\nAkrTulGMG5rEzPXZzFi7y+04JsBY4TCu+T41l8UZedw1pB3Rtbw5+01wGjOwNe0bx/DwlDUcKrYT\n5ab6WOEwrlBVnv1mAwn1o7jqlJZuxwlI4aEhPHJRV7L2HuLlWWluxzEBxAqHccX0NbtYmbmPP5yV\nRESYfQy9pX+bhlzcsxkT5m5m6247UW6qh/2PNTWurFx5bsYG2sRFc0mv3yzmaKrZA+d2IixEeOzr\ntW5HMQHCCoepcV+t3M7GXfncPbS9TS1SAxrXieSOM9vxzdpdzE21ZZXNybP/taZGlZaV8/yMjXRs\nEst53Zq6HSdo3DSoNa0a1uYfX66lpKzc7TjGz3lzBcBIEVksIitEZI2I/MNpby0ii0QkTUQ+EpEI\np72Ws53mPJ5Y4bUecNo3iMg53spsvO/z5dvJ2H2Qe4a1J8RW9qsxkeGh/PW8zqRl5/POgi1uxzF+\nzps9jiJgiKr2AHoCw50lYZ8CnlfVdsAe4CZn/5uAPU77885+iEhnYCTQBRgOvOKsKmj8TGlZOf/+\nLpXOTeswrHNjt+MEnaGd4jktqRHPf7uR3flFbscxfsxrhUM98p3NcOemwBDgY6d9Ep51xwEucrZx\nHj/LWZf8IuBDVS1S1XQ8S8v29VZu4z1frvT0Nu46K8lWqXOBiPDQBZ05WFzG899udDuO8WNePcch\nIqEishzIBmYAm4C9qlrq7JIJHB5W0xzYBuA8vg9oWLH9CM+p+F5jRSRFRFJycuwEoK8pK1de+i6N\njk1iOdt6G65pFx/LqH4teX/RVjbuOuB2HOOnvFo4VLVMVXsCCXh6CV6bM1tVJ6hqsqomx8XFeett\nzAn6auV2NucUcNdZSXZuw2XjhrYnplYYj361FlV1O47xQzUyqkpV9wKzgAFAPRE5PL9EApDl3M8C\nWgA4j9cFdldsP8JzjB8od3ob7RvHMLxLE7fjBL360RHcdVYSc1Nzmb3Beuem6rw5qipOROo596OA\nYcA6PAXkcme30cAXzv0pzjbO49+p58+hKcBIZ9RVayAJWOyt3Kb6TV29g7TsfH4/xHobvuL6AYm0\nbhTNY1/b8FxTdd7scTQFZonISmAJMENVvwL+DNwjIml4zmG84ez/BtDQab8HuB9AVdcAk4G1wDTg\nDlW1Gdv8hKry7+/SaBMXzbl23YbPiAgL4cFzO7Epp4D3F211O47xM16bklRVVwK9jtC+mSOMilLV\nQuCKo7zW48Dj1Z3ReN+sDdms33mAZ67oQaj1NnzK0E7xDGjTkH99u5FLejenTmS425GMn7Arx41X\nvTJrE83rRXFRz2ZuRzG/IiL85bxO7DlYwvjZm9yOY/yIFQ7jNYvT80jZsoexp7ch3Oak8kldm9fl\n0l7NeWNeOll7D7kdx/gJ+99svOblWWk0jI7gyuQWle9sXHPvOR0AeGb6BpeTGH9hhcN4xeqsfczZ\nmMONg1oTFWEzxPiy5vWiuGlQaz5blsWqzH1uxzF+wAqH8YrxszcRWyuM6wa0cjuKOQ63DW5Lg+gI\nHp9qFwWaylnhMNUuI7eAqat3MGpAKxup4yfqRIbzh7OSWLg5zy4KNJWywmGq3WtzNxMeEsKYgYlu\nRzFVcHXfliQ2rM0T/1tHWbn1OszRWeEw1So3v4iPl2Zyae/mxMdGuh3HVEFEWAh/OqcjG3fl88nS\nTLfjGB9mhcNUq7fnZ1BcVs7vTm/jdhRzAs7t1oQeLerx3IyNHCq2CRrMkVnhMNXmYHEpby/cwtBO\njWkbF+N2HHMCRIQHR3Rk5/5CJv6Q7nYc46OscJhqM3nJNvYeLOEW6234tX5tGjK0Uzyvzt5EXkGx\n23GMD7LCYapFaVk5r89Lp0+r+iQnNnA7jjlJfx7ekYLiUv79XZrbUYwPssJhqsXU1TvJ3HOIsdbb\nCAhJjWO5ok8L3lmYwba8g27HMT7GCoc5aarK63M307pRNMM62bKwgWLcsCRCRHhuhq1Pbn7JCoc5\naUsy9rAycx83DmptCzUFkKZ1o7hxUGs+X57Fmu02FYn5mTdXAGwhIrNEZK2IrBGRPzjtD4tIlogs\nd27nVnjOAyKSJiIbROScCu3DnbY0EbnfW5nNiXl97mbq1Q7n8t4Jbkcx1ezWM9pSJzKcp6bZBIjm\nZ97scZQC96pqZ6A/cIeIdHYee15Vezq3qQDOYyOBLsBw4BURCRWRUOBlYATQGbi6wusYl6XnFjBj\n3S5G9WtlkxkGoLpR4dx5Zju+35jDD2m5bscxPsJrhUNVd6jqj879A3jWG29+jKdcBHyoqkWqmg6k\n4VkpsC+QpqqbVbUY+NDZ1/iAN39IJzwkhOtPtckMA9V1A1rRrG4kT01bbxMgGqCGznGISCKeZWQX\nOU13ishKEZkoIvWdtubAtgpPy3TajtZuXLb3YDH/Tcnkwp7NbHqRABYZHsrdw9qzMnMfU1ftdDuO\n8QFeLxwiEgN8AoxT1f3AeKAt0BPYATxbTe8zVkRSRCQlJ8dm96wJ7y3ayqGSMm4+rbXbUYyXXdo7\ngfaNY3h6+npKysrdjmNc5tXCISLheIrGe6r6KYCq7lLVMlUtB17DcygKIAuouFRcgtN2tPZfUNUJ\nqpqsqslxcXHV/82YXyguLWfS/AxOS2pExyZ13I5jvCw0RLjvnI5k7D7IR0u2Vf4EE9C8OapKgDeA\ndar6XIX2phV2uwRY7dyfAowUkVoi0hpIAhYDS4AkEWktIhF4TqBP8VZuc3ymrtpB9oEibhxkvY1g\ncVaneE5JrM8LM1M5WFzqdhzjIm/2OAYC1wFDfjX09p8iskpEVgJnAncDqOoaYDKwFpgG3OH0TEqB\nO4HpeE6wT3b2NS5RVd6Yl07buGjOSLLeXbAQEe4f0ZGcA0VMnGcTIAazMG+9sKrOA450NdjUYzzn\nceDxI7RPPdbzTM1K2bKHVVn7eOzirnbBX5Dp06oBwzo35j9zNnNNv1Y0iI5wO5JxgV05bqps4rx0\n6kaFc2lvG9wWjO47pwMFxaW8MssmQAxWVjhMlWzLO8j0NTu5pl9Lakd4rcNqfFhS41gu653A2wu2\nkLnHJkAMRlY4TJVMmp+BiHD9ALvgL5jdPaw9CDw/I9XtKMYFxywcItLgOG71aiqscVd+USkfLdnG\nud2a0rRulNtxjIua1YvihlMT+XRZJut37nc7jqlhlfU4tgMpwNJj3FZ6M6DxHR+nbONAUSk3Dkx0\nO4rxAbcPbktMrTCetgkQg05lB6nXqWqvY+0gIsuqMY/xUeXlyqQFW+jVsh69Wtav/Akm4NWrHcFt\ng9vyz2kbWJKRxym28mPQqKzHMeA4XuN49jF+bvbGbNJzCxgz0C74Mz8bc2prGtepxRNT19kEiEHk\nmIVDVQsrbotIvIi0PHw70j4mML35QwZN6kQyomsTt6MYHxIVEcq4oe35cetevlm7y+04poYc16gq\nEblQRFKBdGAOkAH8z4u5jA9J3XWAuam5XDegFeGhNhDP/NIVfRJoGxfNP6etp9QmQAwKx/tb4FE8\nizFtVNXWwFnAD15LZXzKm/MzqBUWwtV9W7odxfigsNAQ7hvekU05Bfx3aabbcUwNON7CUaKqu4EQ\nEQlR1Vl4pkU3AW7vwWI+/TGTi3s2t+klzFGd3bkxvVvW4/kZGzlUXOZ2HONlx1s49jrranwPvCci\nL+BZGtYEuA+XbKOwpJwxgxLdjmJ8mIjwwLmdyD5QxMQfbALEQHe8heMi4CCemWynAZuAC7wVyviG\n0rJy3p6fQf82DWzNDVOpUxIbMLRTY16dvYm8gmK34xgvqrRwiMjFwG3AMFUtVdVJqvqic+jKBLAZ\na3exfV+hDcE1x+3Pwz0TIP77O5sAMZBVNuXIK3h6GQ2BR0XkbzWSyviEN+dnkFA/iqGdGrsdxfiJ\npMaxXJncgncWZrB1t02AGKgq63GcDgxR1QeAwcDFXk9kfMKa7ftYnJ7H6AGJhNqaG6YK7h7WntAQ\n4elvbCqSQFVZ4ShW1TIAVT3IkRdmOiIRaSEis0RkrYisEZE/OO0NRGSGiKQ6X+s77SIiL4pImois\nFJHeFV5rtLN/qoiMrvq3aarqrR8yiAoP5cpTWlS+szEVNK4Tyc2D2vDliu2szNzrdhzjBZUVjo7O\nL/GVIrKqwvbhpV+PpRS4V1U747kG5A4R6QzcD8xU1SRgprMNMALPOuNJwFhgPHgKDfAQ0A/oCzx0\nuNgY79idX8QXK7ZzWZ/m1I0KdzuO8UO3nNGGBtERPDF1vU1FEoAqm+Sw04m+sKruAHY49w+IyDqg\nOZ4RWoOd3SYBs4E/O+1vq+dTtlBE6olIU2ffGaqaByAiM4DhwAcnms0c2weLt1JcWs4Npya6HcX4\nqdjIcO4a0o6Hv1zL7A05nNkx3u1IphpVNlfVlmPdjvdNRCQR6AUsAho7RQVgJ3D4zGtzYFuFp2U6\nbUdr//V7jBWRFBFJycnJOd5o5ldKysp5Z+EWTktqRLv4WLfjGD92Tb9WtGpYmyf+t86mIgkwlY2q\n+qqyF6hsH+fCwU+Acar6ixVfnN5FtfRjVXWCqiaranJcXFx1vGRQ+t/qnezaX8QYW3PDnKSIsBDu\nH96Rjbvy+dimIgkolR2qGiQiU47xuACdj/qgSDieovGeqn7qNO8SkaaqusM5FJXttGcBFc/EJjht\nWfx8aOtw++xKcpsT9NYP6SQ2rM3g9nZowZy84V2b0KdVfZ6dsZELejQjupatUx8IKvtXvOg4XuOI\nl4iKiABv4FkM6rkKD00BRgNPOl+/qNB+p4h8iOdE+D6nuEwH/q/CCfGzgQeOI5epohXb9vLj1r08\ndEFnQmwIrqkGIsKD53bisvHzeW3uZsYNbe92JFMNjlk4VHXOSbz2QOA6YJWILHfaHsRTMCaLyE3A\nFuBK57GpwLlAGp7pTcY4GfJE5FFgibPfI4dPlJvq9eYP6cTUCuPyPgluRzEBpE+r+pzXrSn/mbOZ\na/q2JL5OpNuRzEnyWr9RVedx9Os+zjrC/grccZTXmghMrL505tey9xfy9aodXNuvFbGRNgTXVK/7\nhnfgm7U7eW7GRp68rLvbccxJslV5DADvLtpKabnaEFzjFa0aRnNd/0Qmp2xj3Y79lT/B+LTjXQGw\ntoh0d261vB3K1Kyi0jLeX7SFMzvEk9go2u04JkDddVY7YiPDefxrW5/c31U2HDdcRP6F59qJN4G3\ngM0icr/zuC3mFAC+WrGD3Pxi620Yr6pXO4JxQ5OYl5bLrA3ZlT/B+KzKehzPAjFAK1Xto6q98VxN\n3kZExgOfeTug8S5V5c356bSNi+a0pEZuxzEBblT/VrRpFM1jX6+jxC4K9FuVFY5zgd+p6oHDDc5F\nfLcBI4GrvZjN1ICULXtYnbWfMQNb4xlBbYz3hIeG8OC5ndicU8D7i7a6HcecoMoKR7ke4WCkM2Nu\njqou9E4sU1MmzkunblQ4l/b+zSwuxnjFWZ3iGdiuIc9/u5F9B0vcjmNOQGWFY62IXP/rRhEZBazz\nTiRTUzL3HGT6mp2M7NuC2hF2Ra+pGSLCX8/rzP5DJbwwM9XtOOYEVPbb4g7gUxG5EVjqtCUDUcAl\n3gxmvO+dBVsQEa4fkOh2FBNkOjWtw8i+LXl7QQbX9GthE2r6mcpmx81S1X7AI0CGc3tEVfuqapb3\n4xlvOVhcygeLtzK8SxOa14tyO44JQvcOa09URCiPfGXDc/3NcV3HoarfqepLzm2mt0MZ7/vkxyz2\nF5baLLjGNQ1jajFuaHu+35hjw3P9jF05HoTKy5W3fkine0Jd+rSyxRSNe64f0Iq2cdE8+tU6iktt\neK6/sMIRhOak5rApp4AxAxNtCK5xVXhoCH87vzPpuQW8NT/d7TjmOFnhCEIT56UTH1uL87o1czuK\nMQzuEM+QjvG8ODON7P2Fbscxx8EKR5DZsPMAc1NzGX1qIhFh9s9vfMPfz+9McWk5T/xvvdtRzHGw\n3xxBZuK8dCLDQ7imb0u3oxjzk8RG0Yw9vQ2fLcticbott+PrvFY4RGSiiGSLyOoKbQ+LSJaILHdu\n51Z47AERSRORDSJyToX24U5b2uHJFc2Jyc0v4rPlWVzaO4H60RFuxzHmF24/sy3N6kby9y9WU2rz\nWPk0b/ajFgc8AAAWbElEQVQ43gKGH6H9eVXt6dymAohIZzxzX3VxnvOKiISKSCjwMjACz9rmVzv7\nmhPw3sKtFJeWc+PA1m5HMeY3akeE8bfzO7N+5wHes3msfJrXCoeqfg8cb5/zIuBDVS1S1XQ8y8f2\ndW5pqrpZVYuBDzm+ddDNrxSVlvHOwi0M7hBHu/gYt+MYc0TDuzZhULtGPPPNBnIOFLkdxxyFG+c4\n7hSRlc6hrMMXETQHtlXYJ9NpO1r7b4jIWBFJEZGUnJwcb+T2a1OWbyc3v4ibBllvw/guEeHhC7tQ\nWFLGE1NtOjxfVdOFYzzQFugJ7MCz3ke1UNUJqpqsqslxcXHV9bIBQVV5fW46HZvEMqidrblhfFu7\n+BhuOb0tny7LYsGm3W7HMUdQo4VDVXepapmqlgOv4TkUBZAFtKiwa4LTdrR2UwXfp+ayYdcBbj6t\njV3wZ/zCnUPa0aJBFH/9fJVdUe6DarRwiEjTCpuXAIdHXE0BRopILRFpDSQBi4ElQJKItBaRCDwn\n0KfUZOZA8PrczcTH1uLCHnbBn/EPkeGhPHJhVzblFPDa3M1uxzG/4rVFGETkA2Aw0EhEMoGHgMHO\nOuWKZ6bdWwBUdY2ITAbWAqXAHc5iUYjIncB0IBSYqKprvJU5EK3bsZ+5qbn86ZwOdsGf8Stndoxn\nRNcmvDgzlQt7NKNFg9puRzIOCcTpjJOTkzUlJcXtGD7hnsnLmbZ6J/PvH0K92nbthvEvO/YdYuiz\nc0hObMBbY06xQ61eJiJLVTW5sv3sT9AAtmt/IV+u2M6VyS2saBi/1LRuFH88pwNzNuYwZcV2t+MY\nhxWOAPbW/AzKytUu+DN+7foBifRsUY9/fLmWvIJit+MYrHAErAOFJby7cAvDuzahZUM7Nmz8V2iI\n8ORl3dh/qITHvl7rdhyDFY6A9eHibRwoLOWW09u6HcWYk9axSR1uPaMtn/6YxdxUu8DXbVY4AlBx\naTlvzEtnQJuG9GhRz+04xlSLO4e0o02jaB78bBUFRaVuxwlqVjgC0OfLs9i5v5BbB1tvwwSOyPBQ\nnrysO5l7DvH09A1uxwlqVjgCTHm58p85m+jUtA6nJ9n0Iiaw9G3dgNEDEnlrfgYLN9t0JG6xwhFg\nZq7PZlNOAbeeYdOLmMB03/AOtGpYm/s+XsnBYjtk5QYrHAHm1TmbSKgfxXndmla+szF+qHZEGP+8\nrDtb8w7yz2l2yMoNVjgCyOL0PJZu2cPvTmtDWKj905rA1a9NQ2441Q5ZucV+uwSQf89Ko2F0BFcm\nt6h8Z2P83H3DO5DYsDb3Tl7BgcISt+MEFSscAWJl5l6+35jDTae1Jioi1O04xnhd7YgwnruqJzv2\nHeIfX9qFgTXJCkeAeHlWGnUiw7iufyu3oxhTY3q3rM+dZ7bj46WZTFu9w+04QcMKRwDYuOsA09fs\n4oaBrYmNDHc7jjE16vdnJdE9oS4PfLqK7P2FbscJClY4AsDLs9KoHRHKmFMT3Y5iTI0LDw3h+at6\ncqikjD99vJLy8sBbKsLXeK1wiMhEEckWkdUV2hqIyAwRSXW+1nfaRUReFJE0EVkpIr0rPGe0s3+q\niIz2Vl5/lZFbwJcrtjOqfyvqR9vU6SY4tY2L4S/ndWbOxhwm/pDudpyA580ex1vA8F+13Q/MVNUk\nYKazDTACz3KxScBYYDx4Cg2elQP74Vmf/KHDxcZ4jJ+9ibDQEG4eZFOnm+A2ql9LzunSmKemrWdl\n5l634wQ0rxUOVf0eyPtV80XAJOf+JODiCu1vq8dCoJ6zPvk5wAxVzVPVPcAMfluMgta2vIN88mMm\n1/RtSXydSLfjGOMqEeGpy7oTF1OL33+wzIboelFNn+NorKqHhz7sBBo795sD2yrsl+m0Ha39N0Rk\nrIikiEhKTk5wTLv80nephIQIt9lkhsYAUK92BC9c3YtteQf56+erCcSlsX2BayfH1fMvWm3/qqo6\nQVWTVTU5Li6uul7WZ23ZXcAnP2ZxTd+WNLbehjE/OSWxAXcPbc8Xy7fz4ZJtlT/BVFlNF45dziEo\nnK/ZTnsWUPFy5wSn7WjtQe+l79IICxFut96GMb9x+5ntOC2pEQ99scbOd3hBTReOKcDhkVGjgS8q\ntF/vjK7qD+xzDmlNB84WkfrOSfGznbaglpFbwGfLsri2Xys7t2HMEYSGCC+M7EWjmAhue/dH9h60\ntcqrkzeH434ALAA6iEimiNwEPAkME5FUYKizDTAV2AykAa8BtwOoah7wKLDEuT3itAW1w72NWwe3\ncTuKMT6rQXQEr4zqQ/aBQsZ9tNyu76hGYd56YVW9+igPnXWEfRW44yivMxGYWI3R/Fpadj6fLctk\nzMDWxMdab8OYY+nZoh5/v6ALf/t8NS9+l8q4oe3djhQQ7MpxP/PcjA1EhYfaSCpjjtOofi25tHdz\n/vVtKtNW73Q7TkCwwuFHVmzby9RVO7n5tDY0iqnldhxj/IKI8H+XdKNHi3rcM3k563bsdzuS37PC\n4Ueenr6BBtER3HyaXSVuTFVEhocy4bo+xEaGcfOkFHbnF7kdya9Z4fATP6TlMi8tl9sHt7UZcI05\nAY3rRDLhumRy8ou47b0fKS4tdzuS37LC4QdUlX9OW0+zupGMsvU2jDlhPVrU4+nLu7M4PY8/f7LS\nriw/QV4bVWWqz7TVO1mRuY9/Xt6dyHBb3c+Yk3FRz+ZsyzvIM99sJKF+FPee3cHtSH7HCoePKyot\n44n/rScpPoZLex1xmi5jTBXdcWY7Mvcc4qXv0mheL4qRfVu6HcmvWOHwcZPmZ7A17yCTbuxLWKgd\nWTSmOogIj17clR37CvnL56tpXCeSMzvGux3Lb9hvIh+2O7+Il2amMbhDHGe0D/yJG42pSeGhIbx8\nbW86NY3ltveWsjg96CelOG5WOHzY899u5GBJGX89r5PbUYwJSDG1wpg0pi/N60Vx01tLWJ21z+1I\nfsEKh4/auOsA7y/ayrX9WtIuPtbtOMYErIYxtXjnpn7UiQrn+omLScvOdzuSz7PC4YNUlUe/WktM\nrTCbW8eYGtCsXhTv3tyPEBFGvb6I9NwCtyP5NCscPmja6p3MTc1l3ND2NIiOcDuOMUGhdaNo3r25\nL8Vl5YycsMCKxzFY4fAx+UWl/OPLtXRqWofrB9jFfsbUpI5N6vDB7/pTWqZc9Z8FbMqxw1ZHYoXD\nx7w4M5Wd+wt57OKuNvzWGBd0aBLL+7/rT1m5MnLCQlJ3HXA7ks+x30w+ZMPOA7wxL52Rp7SgT6v6\nbscxJmh1aBLLB2P7A3DFfxawbOselxP5FlcKh4hkiMgqEVkuIilOWwMRmSEiqc7X+k67iMiLIpIm\nIitFpLcbmb2tvFz56+erqBMZxp+Hd3Q7jjFBr33jWD659VTqRIZz7euL+H5jjtuRfIabPY4zVbWn\nqiY72/cDM1U1CZjpbAOMAJKc21hgfI0nrQEfpWxjScYeHhjRifp2QtwYn9CyYW0+vm0ArRpGc9Ok\nJUxZsd3tSD7Blw5VXQRMcu5PAi6u0P62eiwE6olIUzcCekvW3kM8/vU6Tm3bkMv7JLgdxxhTQXxs\nJB+O7U+vlvW564NlvDQzNehn1XWrcCjwjYgsFZGxTltjVd3h3N8JNHbuNwe2VXhuptP2CyIyVkRS\nRCQlJ8d/upSqyv3O9M5PXdadkBBxO5Ix5lfqRoXzzk19uaRXc56dsZF7Jq+gqLTM7ViucWuSw0Gq\nmiUi8cAMEVlf8UFVVRGpUklX1QnABIDk5GS/+XPgg8XbmJuay2MXd6VFg9puxzHGHEWtsFCeu7IH\nbeOieeabjWzLO8gro3oTHxvpdrQa50qPQ1WznK/ZwGdAX2DX4UNQztdsZ/csoEWFpyc4bX4vc89B\nHv96LQPbNeTafjatszG+TkS4c0gSL1/Tm9Xb93H+i/NYkhF8kyPWeOEQkWgRiT18HzgbWA1MAUY7\nu40GvnDuTwGud0ZX9Qf2VTik5bfKypU//ncFAE9d1h0RO0RljL84r3tTPr9jILUjQhk5YSGvz90c\nVOc93OhxNAbmicgKYDHwtapOA54EholIKjDU2QaYCmwG0oDXgNtrPnL1e+m7VBZuzuOhC7uQUN8O\nURnjbzo2qcOU3w/irI7xPPb1Om55Zyl5BcVux6oREohVMjk5WVNSUtyOcVTzN+Vy7euLuKRnc569\nsof1NozxY6rKG/PS+ee0DdSrHc4zV/TgdD9dP0dElla4ROKofGk4blDIzS9i3IfLad0omkcv7mpF\nwxg/JyLcfFobPr9jIHWdqdkfnrKGg8WlbkfzGiscNai8XLln8gr2HSrh5Wt6E13LVu41JlB0blaH\nL38/iBtOTeSt+Rmc86/vmZvqP5cGVIUVjhr01PT1fL8xh4cu6EKnpnXcjmOMqWaR4aE8fGEXJt8y\ngPDQEK57YzH3TF7O7vwit6NVKyscNWTykm38Z85mruvfiqv7tqj8CcYYv9W3dQOm3nUavx/SjinL\ntzP4mdm8MS+dkrJyt6NVCyscNWDBpt08+NkqTktqxEMXdLbzGsYEgcjwUO49uwPTxp1Gr5b1efSr\ntQz/1/d8t36X3w/dtcLhZem5Bdz23lJaNazNv6/pbWtsGBNk2sXHMmnMKUy8IZlyhRvfSuGKVxew\nYNNut6OdMPst5kXb8g4y6vVFCDDxhlOoGxXudiRjjAtEhCEdGzN93Ok8dnFXtu05yNWvLWTU64tY\nsGm33/VA7DoOL8nae4ir/rOA/YdKeP93/enavK6reYwxvqOwpIx3F27h1TmbyM0vpkdCXW45oy3n\ndGlCqIsTnR7vdRxWOLxgx75DjJywkLyCYt67uR/dE+q5lsUY47sKS8r45MdMXvt+Mxm7D5JQP4qr\n+7bkqlNa0CimVo3nscLhUuFIzy1gzJuLyc0v5p2b+tKrpS0Ba4w5trJy5Zs1O3l7wRYWbN5NeKhw\nTpcmXNY7gdOSGtXYudHjLRx2BVo1mr8pl9ve/ZHQEOFtKxrGmOMUGiKM6NaUEd2akpadz3uLtvDp\nj1l8tXIHjWIiOL97My7o0ZReLer7xJo91uOoJh8u3spfP19NYqNoJo4+hZYNbeJCY8yJKyotY/aG\nHD5flsXMddkUl5XTKKYWwzo35uzOjenXpgG1I6r3b387VFVDhWPfwRIe/XotHy/N5LSkRrx8bW/q\nRNroKWNM9dlfWMKs9dl8s2YXszdkU1BcRkRoCMmJ9RmU1Ij+bRrStVldIsJO7pCWFY4aKBzfrt3F\ng5+tYndBMbee0Ya7h7a36zSMMV5VWFLGkow85qbm8v3GHNbvPABAZHgIPVvU44z28dw2uO0JvbYV\nDi8WjvU79/PSzDS+XrWDjk1iefryHnRLsOG2xpial32gkKUZe1ickUdKxh4axUTw5pi+J/RaAXdy\nXESGAy8AocDrqvpkJU+pVqrKj1v3MH72Jr5dl010RCjjhiZx++B2J909NMaYExUfG/nTiXXwzMLt\nbX5ROEQkFHgZGAZkAktEZIqqrvXWexaXlpObX8SqrH3M3pDDnA3ZbN9XSP3a4dwzrD2jByRSt7ad\nyzDG+JaaGHXlF4UD6AukqepmABH5ELgIqNbCkXOgiFGvLyL7QCF7Dpb81B5TK4xB7Rrxh6FxnN+9\nma2jYYwJav7yG7A5sK3CdibQr+IOIjIWGAvQsmXLE3qT2MgwWjWsTXJifeJjI4mLrUWbuGj6tKpP\nuJ30NsYYwH8KR6VUdQIwATwnx0/kNSLDQ5lwfaXnhYwxJqj5y5/RWUDF1Y8SnDZjjDE1zF8KxxIg\nSURai0gEMBKY4nImY4wJSn5xqEpVS0XkTmA6nuG4E1V1jcuxjDEmKPlF4QBQ1anAVLdzGGNMsPOX\nQ1XGGGN8hBUOY4wxVWKFwxhjTJVY4TDGGFMlATk7rojkAFtO4iUaAbnVFCcQ2M/jl+zn8Vv2M/kl\nf/15tFLVuMp2CsjCcbJEJOV4phYOFvbz+CX7efyW/Ux+KdB/HnaoyhhjTJVY4TDGGFMlVjiObILb\nAXyM/Tx+yX4ev2U/k18K6J+HneMwxhhTJdbjMMYYUyVWOIwxxlSJFY4KRGS4iGwQkTQRud/tPG4Q\nkRYiMktE1orIGhH5g9PeQERmiEiq87W+21lrkoiEisgyEfnK2W4tIoucz8pHznT/QUFE6onIxyKy\nXkTWicgA+3zI3c7/l9Ui8oGIRAbyZ8QKh0NEQoGXgRFAZ+BqEensbipXlAL3qmpnoD9wh/NzuB+Y\nqapJwExnO5j8AVhXYfsp4HlVbQfsAW5yJZU7XgCmqWpHoAeen0vQfj5EpDlwF5Csql3xLP0wkgD+\njFjh+FlfIE1VN6tqMfAhcJHLmWqcqu5Q1R+d+wfw/FJojudnMcnZbRJwsTsJa56IJADnAa872wIM\nAT52dgman4eI1AVOB94AUNViVd1LEH8+HGFAlIiEAbWBHQTwZ8QKx8+aA9sqbGc6bUFLRBKBXsAi\noLGq7nAe2gk0dimWG/4F3AeUO9sNgb2qWupsB9NnpTWQA7zpHLp7XUSiCeLPh6pmAc8AW/EUjH3A\nUgL4M2KFwxyRiMQAnwDjVHV/xcfUM4Y7KMZxi8j5QLaqLnU7i48IA3oD41W1F1DArw5LBdPnA8A5\nn3MRnqLaDIgGhrsaysuscPwsC2hRYTvBaQs6IhKOp2i8p6qfOs27RKSp83hTINutfDVsIHChiGTg\nOXw5BM8x/nrOYQkIrs9KJpCpqouc7Y/xFJJg/XwADAXSVTVHVUuAT/F8bgL2M2KF42dLgCRnJEQE\nnpNbU1zOVOOc4/dvAOtU9bkKD00BRjv3RwNf1HQ2N6jqA6qaoKqJeD4T36nqtcAs4HJnt2D6eewE\ntolIB6fpLGAtQfr5cGwF+otIbef/z+GfScB+RuzK8QpE5Fw8x7NDgYmq+rjLkWqciAwC5gKr+PmY\n/oN4znNMBlrimbL+SlXNcyWkS0RkMPBHVT1fRNrg6YE0AJYBo1S1yM18NUVEeuIZKBABbAbG4Pkj\nNGg/HyLyD+AqPKMSlwE34zmnEZCfESscxhhjqsQOVRljjKkSKxzGGGOqxAqHMcaYKrHCYYwxpkqs\ncBhjjKkSKxzGGGOqxAqHMccgIg1FZLlz2ykiWRW253vh/W4QkRwRef0Y+0Q5718sIo2qO4MxlQmr\nfBdjgpeq7gZ6AojIw0C+qj7j5bf9SFXvPEamQ0BPZxoUY2qc9TiMOUEiku98HSwic0RksohsFJEn\nReRaEVksIqtEpK2zX5yIfCIiS5zbwON4jy7O6ywXkZUikuTt78uYyliPw5jq0QPoBOThmYbjdVXt\n66yg+HtgHJ7JEZ9X1Xki0hKY7jznWG4FXlDV95w51EK99h0Yc5yscBhTPZYcXo9CRDYB3zjtq4Az\nnftDgc6eefAAqCMiMaqaf4zXXQD8xVlM6lNVTa3+6MZUjR2qMqZ6VJy8rrzCdjk//4EWAvRX1Z7O\nrXklRQNVfR+4EDgETBeRIdWc25gqs8JhTM35Bs9hK+CnWWaPyZmFd7Oqvohn6vLu3otnzPGxwmFM\nzbkLSHZOcq/Fc/6iMlcCq0VkOdAReNubAY05HjatujE+RERuAJKPNRy3wr4Zzr653s5lTEXW4zDG\ntxwCRhzPBYBAOD8vtmVMjbEehzHGmCqxHocxxpgqscJhjDGmSqxwGGOMqRIrHMYYY6rk/wFdR018\nTxLAzwAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(t, q)\n", "plt.title('Max Q')\n", "plt.xlabel('Time [s]')\n", "plt.ylabel('Q [Pa]')\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.6" } }, "nbformat": 4, "nbformat_minor": 4 }