{ "cells": [ { "cell_type": "markdown", "id": "03bf5cde-0076-4fcc-a97b-05e5c2afaf5b", "metadata": {}, "source": [ "# Using LayeredUniaxial" ] }, { "cell_type": "code", "execution_count": 1, "id": "7b5ad552-a0fc-4f1c-8a58-e17a65e7582d", "metadata": {}, "outputs": [], "source": [ "import pymecht as pmt\n", "from matplotlib import pyplot as plt\n", "import numpy as np" ] }, { "cell_type": "markdown", "id": "8a7da9e9-eb77-4add-903d-be078622e4b8", "metadata": {}, "source": [ "We start by creating two uniaxial samples with different material properties and different **reference length**, one with `L0=20` and second with `L0=22`. We also set the fiber direction to the first axis for both." ] }, { "cell_type": "code", "execution_count": 2, "id": "53debe42-10bd-4f97-ab71-f6be4b9b0fcc", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Fiber directions set to 0 degrees ( 0.0 radians)\n", "Fiber directions set to 0 degrees ( 0.0 radians)\n" ] } ], "source": [ "mat1 = pmt.MatModel('nh','goh')\n", "strip_1 = pmt.UniaxialExtension(mat1,disp_measure='length', force_measure='force')\n", "pmt.specify_single_fiber(strip_1)\n", "params = strip_1.parameters\n", "params.set('mu_0',22.57)\n", "params.set('k1_1',276.45)\n", "params.set('k2_1',42.85)\n", "params.set('k3_1',0.246)\n", "params.set('L0',20)\n", "strip_1.parameters = params\n", "\n", "mat2 = pmt.MatModel('nh','goh')\n", "strip_2 = pmt.UniaxialExtension(mat2,disp_measure='length', force_measure='force')\n", "pmt.specify_single_fiber(strip_2)\n", "params = strip_2.parameters\n", "params.set('mu_0',14.30)\n", "params.set('k1_1',290.22)\n", "params.set('k2_1',4.87)\n", "params.set('k3_1',0.224)\n", "params.set('L0',22)\n", "strip_2.parameters = params" ] }, { "cell_type": "markdown", "id": "c17cda0e-f46c-44a8-b355-30a6f80ae3f5", "metadata": {}, "source": [ "Next, we combine them using `LayeredUniaxial` and check its parameters. We notice that all the parameters for the first \"layer\" have `_layer0` appended to their names, and those for second layer have `_layer1` appended, etc." ] }, { "cell_type": "code", "execution_count": 3, "id": "0f67cd5d-9614-4432-92be-c72b859004a3", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "------------------------------------------------------------------\n", "Keys Value Fixed? Lower bound Upper bound \n", "------------------------------------------------------------------\n", "L0_layer0 20.00 No 1.00e-04 1.00e+03 \n", "A0_layer0 1.00 No 1.00e-04 1.00e+03 \n", "mu_0_layer0 22.57 No 1.00e-04 1.00e+02 \n", "k1_1_layer0 2.76e+02 No 0.10 30.00 \n", "k2_1_layer0 42.85 No 0.10 30.00 \n", "k3_1_layer0 0.25 No 0.00 0.33 \n", "L0_layer1 22.00 No 1.00e-04 1.00e+03 \n", "A0_layer1 1.00 No 1.00e-04 1.00e+03 \n", "mu_0_layer1 14.30 No 1.00e-04 1.00e+02 \n", "k1_1_layer1 2.90e+02 No 0.10 30.00 \n", "k2_1_layer1 4.87 No 0.10 30.00 \n", "k3_1_layer1 0.22 No 0.00 0.33 \n", "------------------------------------------------------------------\n", "\n" ] } ], "source": [ "combined = pmt.LayeredUniaxial(strip_1,strip_2)\n", "print(combined.parameters)" ] }, { "cell_type": "markdown", "id": "b9c6b8a2-557f-46ca-9d6d-d30415a56630", "metadata": {}, "source": [ "We impose a deformed length varying between 20 and 25 and calculate and plot the resulting force." ] }, { "cell_type": "code", "execution_count": 4, "id": "d4f69603-d0d5-4816-97d9-4d6d004678dc", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAANIAAACfCAYAAACMX4fwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8g+/7EAAAACXBIWXMAAAsTAAALEwEAmpwYAAATQ0lEQVR4nO2de7iVVZ3HP1/OOcIRUEQQFUG8IKZdREkzLe0mTE+jPF5Sx2myetR6stEehzSnJ5uaypHRsnGyUce0SXN0NEZLxUvepkIFvCAKQkjZQYEkEOEA5/KbP961Y7N99z7v2ft99/X3eZ7znP2u9/ZbsL9nrfVbv7V+MjMcx6mMIbU2wHGaAReS46SAC8lxUsCF5Dgp4EJynBRwITlOCrTX2oBKGDNmjE2aNKnWZjgtwoIFC/5kZmPjzjW0kCZNmsT8+fNrbYbTIkj6fbFzDS0kx0mLOc90MXvuUlat72bvUZ3Mmj6FmVPHJ77fheS0PHOe6eKrdy2iu6cPgK713Xz1rkUAicXkzgan5Zk9d+lfRJSju6eP2XOXJn6GC8lpeVat7x5UeRwuJKfl2XtU56DK43AhOS3PrOlT6GjTDmWdHW3Mmj4l8TNcSE7LM3PqeCbtvjPtQ4SA8aM6+e7J73KvneMMhj9v2saKP23mvA/uz1dmHFzWM7xFclqeXy1ZQ1+/Mf3QPct+hgvJaXnmLn6dPXcZxrv32bXsZ7iQnJame1sfjy9bywmHjkPSwDcUwYXktDSPvbyWLT39FXXrwIXktDgPLH6dXTs7OHK/0RU9x4XktCw9ff089NJqPvKOPehoq0wKLiSnZXlyxTre3NJbcbcOfB7JaUFySya6Qizdhs09FT/TheS0FIVLJgAuu3sxO7UPGVQkQyHetXNaijSWTMThQnJaijSWTMThQnJaijSWTMThQnJailnTp7BThUsm4nAhOS3FzKnj2X/sCIaIspdMxOFeO6elWLW+m6WrN3L+hw7kohMqa4XySdQiSdpX0kfD505JI1OzwHGqyJ0L/ogZnHbEhFSfO6CQJJ0D/A/wH6FoH2BOqlY4ThXo7zduX/AqR++/OxN33znVZydpkb4IHAO8CWBmy4A9UrXCcarAvBVv8Oq6bk5/b7qtESQT0lYz25Y7kNQODJgvU9IESY9IelHSYkkXhPLRkh6UtCz83i2US9IPJC2X9Lykw8utlOPEcfv8Vxk5rJ0Z76w8tq6QJEJ6TNKlQKekjwF3APckuK8XuMjMDgHeB3xR0iHAJcDDZjYZeDgcA/wVMDn8nAtcO6iaOE4JNnT3cN8Lr3PSYXszrKMt9ecn8dpdAnwOWAScB9wL3DDQTWb2GvBa+LxR0kvAeOAk4Phw2c3Ao8DFofwnFmWHnidplKS9wnMcpywKA1THjRyWyXuSCKkTuNHMrgeQ1BbKNid9iaRJwFTgSWBcnjheB8aFz+OBV/Nu+2MocyE5ZREXoPrDR5czYfTOFc8bFZKka/cwkXBydAIPJX2BpBHAncCFZvZm/rnQ+gw43ip43rmS5kuav3bt2sHc6rQY8QGq/RUHqMaRREjDzOyt3EH4nMh3KKmDSES3mNldoXi1pL3C+b2ANaG8C8h3p+wTynbAzK4zs2lmNm3s2NicT44DZBegGkcSIW3K96BJOgIY0BJFW7L8J/CSmV2Vd+pu4NPh86eB/80r/7vgvXsfsMHHR04lZBWgGkeSMdIFwB2SVhGFJ+0JnJ7gvmOATwGLJD0byi4FLgdul/Q54PfAJ8O5e4GPA8uJxl+fSVgHx4ll1vQpXHTHc/T1bx89pBGgGkdJIQXHwgeAg4Hc25ea2YBrc83s/4iEF8dHYq43oslfx0mFDx40FmF0drSxpaevrEx8SSkpJDPrk3SmmX0PeCH1tztOhtz8m5X09sN9FxzD5HHZhocm6dr9WtI1wH8Dm3KFZrYwM6scp0I2b+vl5t+u5KPvGJe5iCCZkA4Lv7+ZV2bAh1O3xnFS4ranXmX95h6+cPwBVXnfgEIysw9VwxDHSYttvf3c8MQKjtxvNEfsu1tV3plkGcWukq7KTYJKulJS+dv2O06GzHmmiyO/8xCrNmzhd2veYs4zb5uKzIQk80g3AhuJ3NSfJFpO8eMsjXKccohCgp5nfdjw8Y1N2/jqXYuqIqYkQjrAzC4zsxXh55+A/bM2zHEGSxQS1L9DWRp71iUhiZC6JR2bO5B0DAkiGxyn2lQzJKiQJF67zwM/yRsX/ZntIT6OUzeMHNbOm1t631aeRUhQIUWFJOkCM7saGGFm75G0C0BhBLfj1ANvvLWVrT19DBHkRQRlFhJUSKmuXS7W7d8gEpCLyKlXvv/QMnoNLp4xhfGjOlPdsy4Jpbp2L0laBuwt6fm8chGFxr07W9McJxnL12zk1qf+wFlHTeS84w7kvOMOrLoNRYVkZmdK2hOYC5xYPZMcJxn5y8gFHDRuRM1sGSho9XXgPVWyxXESU7iM3IBv/3IJI4Z2VKUrV4jv/e00JFnlOSoXF5LTkNRyzigOF5LTkIwdOTS2vBpzRnGUmke6hxI7/JiZOyCcmtDb18/Qjre3AdWaM4qjVIv0r8CVwCtEIUHXh5+3gN9lb5rjxHP9E6/w6rpuPnX0vjWZM4qjlPv7MQBJV5rZtLxT90ian7lljhPDstUb+d6DLzPj0D355omH8q2T3llrk4BksXbDJe1vZisAJO0HDM/WLMfZkTnPdHHF3CWsWr+FIYL3H7g70Y5v9UESIX0ZeFTSCqKohn2J9gB3nKpQOGfUb/Dde5ewy7DazBnFkWSp+f2SJhNtyQWwxMy2ZmuW42yn1JxRvQgpyVLznYFZwPlm9hwwUdInMrfMcQL1NmcUR5J5pB8D24Cjw3EX8M+ZWeQ4efT29dPRHv81rdWcURxJl5pfAfQAmNlmiu+g6jipcvl9S9jW209H245fuVrOGcWRxNmwTVInYXJW0gGAj5GczChMDvaByWM45fB9mD13KavWd2e69XC5JBHSZcD9wARJtxBtjn92lkY5rUtccrCnV67jlMP34deX1O+epAN27czsQeBkIvH8DJhmZo9ma5bTqsR56LZklBwsTYoKSdLB4ffhRHNHrwGriLx2A2Ycl3SjpDWSXsgr84zmTkkawUMXR6mu3UXAOUTxdoUk2fv7JuAa4Cd5ZbmM5pdLuiQcX8yOGc2PIspoflQC+50montbHx1tQ9jW1/+2c/XkoYujVKzdOeF3WXt/m9njIQlzPp7R3NmBnGNh1fpudmqPRNTRJnr6sk8OlialllGcXOrGvJywg6HijOaSzgXOBZg4cWIZJjj1QqFjYWtwc5/+3gk8smRt3Xro4ijVtfvrEucMKEdI2x9gZpIGldE83HcdcB3AtGnTBn2/Uz/EORZ6+oxHlqytaw9dHKW6dlnkcF2d67KVk9HcaS4a1bEQR5JYu92DR22hpAWSrpa0e5nv84zmDgBbevoY2gChP0lJMiF7G/A4cEo4PosoDeZHS90k6WdEjoUxkv5INLHrGc1bmHzHQs4714iOhTgUOcpKXCC9YGbvLChbZGbvytSyBEybNs3mz/fFuo1AXMRCozkWJC0oWC3+F5K0SA9IOgO4PRyfSrT7quMkppkcC3Ekif4+B7iVaCnFNqKu3nmSNkryTfWdRHQ1kWMhjiQrZLPPre40HfnjoRFDi3/NGtGxEEeSrl1ucvZYovmjJ8xsTpZGOY1N4Xho49ZeBLQ3iWMhjiTu7x8SZe1bBLwAfF7Sv2dtmNO4xI2HDBi+U3vd7EOXNklapA8D7whxcEi6GVicqVVOQ5Hfjdtr1DBWrd8Se92G7h6eveyEKltXHZIIaTkwkWjeB6IIhOWZWeQ0FIXduGIiguYZD8WRREgjibL3PUXUQh8JzJd0N/ge4K1OXDcOQlrHvONmGg/FkURIX8/cCqdhyO/G7T2qs6hb24jGQY0w0ZoGSdzfj+UfSzoWONPMvpiZVU5dUtiNKyYiiETUDBOtSUnq/p4K/A1wGlF2ijuzNMqpT4p14wpp9m5cHKUW9h0EnBl+/kQUqKpyV8w6jUfSbhy0VjcujlIt0hLgCeATZrYcQNKXq2KVU3O8Gzc4SgnpZOAM4BFJ9xPF2PkOq01Mfgs0RKIvZmVAq3njklI0ssHM5pjZGURZKB4BLgT2kHStpOacVWthci1Q1/puDGJFBNu9cc0YnVAJSbx2m4iiv28N+9CdRrTzzwMZ2+ZkSOH4Z/O23kSOBO/GxZPIa5fDzP5MtPHIddmY41SDwYx/8vFuXHEGJSSnMSm39QFok+g3a1lvXFJcSE1Oua0PRC2Qj4GS4UJqMippfUZ1djB8aHtLzweViwupiai09fnGiYe6cMrEhdTAFLY+m7z1qRkupAahUDQfOngsdy7o8tanTnAh1SFJRPPTeX9I/DxvfbLHhVRj0hZNId76VAcXUhXJWjTgrU+tcCGlRJxI8rfizUI0cQGk3vrUhgH3/q5narX390AtSxbEieaUI8Y3zL7ZzUCle3+3NEm6Y7fM+wNp/jly0TQedSUkSTOAq4E24AYzuzyrdxUKJBeMWc4YphIRuWiag7rp2klqA14GPkaUQ/Zpok1WXix2T9KuXZKuWMcQgdhhS920cdE0No3StTsSWG5mKwAk3UaU7byokOIo1zPW01+ZgApFUoiLprmpJyHFZTY/qvCiUlnN42LNKvWMxZGkZSn02rlompt6ElIiSmU1T7pd1GDw7piThHoSUsWZzQeTtKpQIHFjJBeNk5R6EtLTwGRJ+xEJ6AyiTSkTU2zvtaStCvA2T56LxklC3QjJzHolnU+Un7YNuNHMBpU+Ztb0KW9L+DvYVsWF45RD3QgJwMzuBe4t9/6cCLxVcapNXQkpDWZOHe/CcapO3UzIloOktWxPgFbIGKI9y5sZr2N12dfMxsadaGghlULS/GKz0M2C17F+GDAZs+M4A+NCcpwUaGYhtcK2yl7HOqFpx0iOU02auUVynKrRFEKSNEHSI5JelLRY0gWhfLSkByUtC793q7Wt5VCifrMlLZH0vKSfSxpVY1PLplgd885fJMkkjamVjaVoiq6dpL2AvcxsoaSRwAJgJnA2sM7MLpd0CbCbmV1cO0vLo0T99gF+FcKr/gWgEesHxetoZi9KmgDcQJT07ggzq5d5pb/QFC2Smb1mZgvD543AS0Trm04Cbg6X3Uz05Ws4itXPzB4ws95w2TwiYTUkJf4PAb4HfIXKVvVnSlMIKR9Jk4CpwJPAODN7LZx6HRhXK7vSoqB++XwWuK/qBmVAfh0lnQR0mdlztbWqNE0VaydpBHAncKGZvSltzx1tZiapbv+iJaGwfnnl/wj0ArfUyra0yK8jUZ0uBeo+Z3HTtEiSOoj+A24xs7tC8erQ9871wdfUyr5KKVI/JJ0NfAI4yxp8wBtTxwOA/YDnJK0k6roulLRn7ayMp1mcDSIaA60zswvzymcDb+Q5G0ab2VdqZGbZlKjfDOAq4DgzW1sj81KhWB0LrlkJTKtHZ0OzCOlY4AlgEdAfii8lGkfcDkwkihL/pJmtq4mRFVCifj8AhgJvhLJ5Zvb56ltYOcXqGNao5a5ZiQvJcZqXphkjOU4tcSE5Tgq4kBwnBVxIjpMCLiTHSQEXUspI6pP0bIhgfi5ELQ/47xwiuReHua+qI+l4Sb9IWp7C+2ZKOiTv+FFJdb83QzGaKkSoTug2s8MAJO0B3ArsAlw2wH3nEk0YJ9q8XFJ7XsBqIzIT+AWDzDZSr3iLlCFmtoZIIOcroi20PE+HNUTnAUi6GxgBLJB0uqRJkn4VrnlY0sRw3U2SfiTpSeCKcHytpHmSVoTW40ZJL0m6KWeHpBMk/VbSQkl3hHg2JM0I65kWAicPVB9Jw8Pzn5L0TAgoRdLZku6SdH9Y+3VF3j2fk/RyuOd6SddIej9wIjA7tN4HhMtPC9e9LOkDlf8PVBEz858Uf4C3YsrWE0Wenwt8LZQNBeYD+xXeB9wDfDp8/iwwJ3y+ieiveFve8W1E25ufBLwJvIvoD+QC4DCifeEeB4aHey4Gvg4MI0qjMzncfzvwixjbj8+VA98B/jZ8HkWUGG440bqvFcCu4bm/J0qIsDewEhgNdBBFLlyTZ/upee95FLgyfP448FCt/y8H8+Ndu+pyAvBuSaeG412JvsivFFx3NNtbiP8Crsg7d4ft2P27x8xM0iJgtZktApC0GJhEFOh5CPDrEA2/E/BbokVyr5jZsnD9Twl5pwaw/0RJ/xCOhxGFXwE8bGYbwrNeBPYlEvFjFsKyJN0BHFTi+blg3AXB9obBhZQxkvYH+ogizwV8yczmVvDITQXHW8Pv/rzPueP28O4HzezMArsOK+PdAk4xs6UFzzqq4N19lPfdyj2j3Ptrho+RMkTSWOBHRN0ZI8q08YWwXABJB0kaHnPrb4jS2gCcRdQlKpd5wDGSDgzvHC7pIGAJMClvfHJmsQfkMRf4UojURtLUAa5/GjhO0m6S2oFT8s5tBEYOoh51TUOpvkHolPQs0Zigl6hrdlU4dwNRl2Vh+DKuJX75+5eAH0uaFa75TLnGmNnasGbpZ5KGhuKvmdnLitKI/lLSZiKxDvTF/hbwfeD54NJ/hWgtVLF3d0n6DvAUsI5IvBvC6duA6yX9PXBqkUc0DB797WSKpBFm9lZokX5OlPfq57W2K228a+dkzTdCC/0CUQs2p6bWZIS3SI6TAt4iOU4KuJAcJwVcSI6TAi4kx0kBF5LjpIALyXFS4P8BKhSIByxbJ20AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "def_length = np.linspace(20,25,30)\n", "force = combined.disp_controlled(def_length)\n", "\n", "fig,ax = plt.subplots(1,1,figsize=(4*0.7,3*0.7))\n", "ax.plot(def_length,force,'-o')\n", "ax.set_xlabel('Deformed length')\n", "ax.set_ylabel('Applied force')\n", "plt.show()" ] }, { "cell_type": "markdown", "id": "1dd6ec6d-06b8-4e45-8f08-65e2990e86cc", "metadata": {}, "source": [ "Moreover, we calculate the \"free\" length of the layered structure by applied a zero force and find that it is somewhere between the lengths of the two layers. " ] }, { "cell_type": "code", "execution_count": 5, "id": "fd019bad-7ffb-4562-ab5f-0ef36995bca5", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Length without applied force is 20.43300896426766\n" ] } ], "source": [ "free_length = combined.force_controlled(0,x0=20)\n", "print('Length without applied force is', free_length)" ] }, { "cell_type": "markdown", "id": "e643c92c-af67-4028-b2ef-d1ea63f54ed6", "metadata": {}, "source": [ "## Required `disp_measure` and `force_measure`" ] }, { "cell_type": "markdown", "id": "2e4b5d0b-7301-468a-8aff-f405cf390607", "metadata": {}, "source": [ "For the `LayeredUniaxial` to work correctly, we need the component layers to have lengths as their `disp_measure` and force as their `force_measure`. Lengths are required to remove ambiguity about what is the stretch/strain etc. is with respect to (i.e., is it the length of the first layer or the second layer?). The force is required, since the forces generated by the component layers will add up, while stresses do not.\n", "\n", "If we do not follow these rules, errors will be raised as shown below." ] }, { "cell_type": "code", "execution_count": 6, "id": "4d567e78-7d89-4c06-9e34-b89e45c2b988", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "ValueError raised: The input of all layers in LayeredUniaxial should be length to remove ambiguity about the reference length\n", "ValueError raised: The output of the LayeredUniaxial should be force, as stresses are not additive. The results may be spurious\n" ] } ], "source": [ "try:\n", " strip_1 = pmt.UniaxialExtension(mat1, disp_measure='stretch', force_measure='force')\n", " strip_2 = pmt.UniaxialExtension(mat2, disp_measure='stretch', force_measure='force')\n", " combined = pmt.LayeredUniaxial(strip_1,strip_2)\n", "except ValueError as e:\n", " print(\"ValueError raised:\", e)\n", "\n", "try:\n", " strip_1 = pmt.UniaxialExtension(mat1, disp_measure='length', force_measure='cauchy')\n", " strip_2 = pmt.UniaxialExtension(mat2, disp_measure='length', force_measure='cauchy')\n", " combined = pmt.LayeredUniaxial(strip_1,strip_2)\n", "except ValueError as e:\n", " print(\"ValueError raised:\", e)" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.9.16" } }, "nbformat": 4, "nbformat_minor": 5 }