bug fixes

Project.toml

@@ -1,7 +1,7 @@
 name = "LowLevelFEM"
 uuid = "6171b9fb-adbf-4751-adb9-5faded75de07"
 authors = ["Balázs Pere"]
-version = "1.3.0"
+version = "1.3.1"
 
 [deps]
 Arpack = "7d9fca2a-8960-54d3-9f78-7d1dccf2cb97"

docs/src/index.md

@@ -41,7 +41,7 @@ The above described steps can be easily performed using the LowLevelFEM package.
 ## Planned features
 
 - [ ] 3D (and  2D) truss structures
-- [ ] 2D axisymmetric problem
+- [x] 2D axisymmetric problem
 - [ ] 3D (and 2D) beam structures
 - [ ] Shells
 - [ ] Giving loads and prescribed displacements with functions
@@ -65,7 +65,7 @@ The above described steps can be easily performed using the LowLevelFEM package.
 
   - [ ] solving conductivity matrix,
   - [ ] solving heat capacity matrix.
-- [ ] Dynamic transient problems with HHT-α (or Newmark).
+- [x] Dynamic transient problems with HHT-α (or Newmark).
 - [ ] Linear buckling.
 - [ ] Modal analysis (eigenfrequencies, modal shapes).
 

examples/axi2D.ipynb

@@ -2,7 +2,7 @@
  "cells": [
   {
    "cell_type": "code",
-   "execution_count": 53,
+   "execution_count": 19,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -14,7 +14,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 54,
+   "execution_count": 20,
    "metadata": {},
    "outputs": [
     {
@@ -30,11 +30,11 @@
       "Info    : [ 60%] Meshing curve 5 (Line)\n",
       "Info    : [ 80%] Meshing curve 6 (Line)\n",
       "Info    : [ 90%] Meshing curve 7 (Line)\n",
-      "Info    : Done meshing 1D (Wall 0.000360992s, CPU 0.000347s)\n",
+      "Info    : Done meshing 1D (Wall 0.000397551s, CPU 0.000398s)\n",
       "Info    : Meshing 2D...\n",
       "Info    : [  0%] Meshing surface 1 (Transfinite)\n",
       "Info    : [ 50%] Meshing surface 2 (Transfinite)\n",
-      "Info    : Done meshing 2D (Wall 5.264e-05s, CPU 5e-05s)\n",
+      "Info    : Done meshing 2D (Wall 5.3636e-05s, CPU 5.4e-05s)\n",
       "Info    : 15 nodes 28 elements\n",
       "Info    : Done reading 'axi2D-2.geo'\n"
      ]
@@ -46,7 +46,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 55,
+   "execution_count": 21,
    "metadata": {},
    "outputs": [
     {
@@ -75,7 +75,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 56,
+   "execution_count": 22,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -84,7 +84,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 57,
+   "execution_count": 23,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -102,7 +102,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 58,
+   "execution_count": 24,
    "metadata": {},
    "outputs": [
     {
@@ -129,7 +129,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 59,
+   "execution_count": 25,
    "metadata": {},
    "outputs": [
     {
@@ -156,13 +156,13 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 60,
+   "execution_count": 26,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "text/plain": [
-       "(\"load\", 0.001, 1, 0)"
+       "(\"load\", 0.0, 1, 0)"
       ]
      },
      "metadata": {},
@@ -172,12 +172,12 @@
    "source": [
     "supp = FEM.displacementConstraint(\"supp\", uy=0)\n",
     "supp2 = FEM.displacementConstraint(\"load\", uy=0.0001)\n",
-    "load = FEM.load(\"load\", fx=0.001, fy=1)"
+    "load = FEM.load(\"load\", fx=0.0000, fy=1)"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 61,
+   "execution_count": 27,
    "metadata": {},
    "outputs": [
     {
@@ -188,9 +188,9 @@
        "   0.0\n",
        "   0.0\n",
        "   0.0\n",
-       "   0.2879793265790645\n",
+       "   0.0\n",
        " 287.9793265790645\n",
-       "   0.1832595714594045\n",
+       "   0.0\n",
        " 183.25957145940447\n",
        "   0.0\n",
        "   0.0\n",
@@ -216,7 +216,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 62,
+   "execution_count": 28,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -225,33 +225,33 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 63,
+   "execution_count": 29,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "text/plain": [
        "30-element Vector{Float64}:\n",
-       " -1.5030922415197964e-5\n",
+       " -1.5000000000000072e-5\n",
        "  0.0\n",
-       " -3.002550724844994e-5\n",
+       " -3.0000000000000082e-5\n",
        "  0.0\n",
-       " -2.978597623883016e-5\n",
-       "  9.98873058636508e-5\n",
-       " -1.4753679882901774e-5\n",
-       "  0.00010008841052820427\n",
-       " -1.4978169334560818e-5\n",
-       "  5.00701699199387e-5\n",
+       " -2.9999999999999916e-5\n",
+       "  0.00010000000000000002\n",
+       " -1.4999999999999882e-5\n",
+       "  0.00010000000000000021\n",
+       " -1.5000000000000019e-5\n",
+       "  5.0000000000000104e-5\n",
        "  ⋮\n",
-       "  7.509798430097896e-5\n",
-       " -1.5019232232159304e-5\n",
-       "  2.5034385772055184e-5\n",
-       " -2.2486569287257476e-5\n",
-       "  5.0007176650063396e-5\n",
-       " -2.2522883208703935e-5\n",
-       "  2.500467877836452e-5\n",
-       " -2.2408238932206295e-5\n",
-       "  7.500392469137584e-5"
+       "  7.500000000000017e-5\n",
+       " -1.5000000000000044e-5\n",
+       "  2.5000000000000035e-5\n",
+       " -2.2500000000000045e-5\n",
+       "  5.000000000000009e-5\n",
+       " -2.2500000000000062e-5\n",
+       "  2.4999999999981028e-5\n",
+       " -2.2499999999999968e-5\n",
+       "  7.500000000000012e-5"
       ]
      },
      "metadata": {},
@@ -264,13 +264,13 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 64,
+   "execution_count": 30,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "text/plain": [
-       "LowLevelFEM.StressField([[0.00032219739849290163; 0.00017984896982550848; … ; 0.0; 0.00019147364595539784;;], [0.000407286612302854; 0.0001747062142590867; … ; 0.0; 0.001219691650751309;;], [0.00025985582385595896; 0.0001564010484018923; … ; 0.0; -0.0001303804410069942;;], [0.000283819864447421; 0.00014622369230693044; … ; 0.0; 0.00044494756753321596;;], [0.0004414574800232496; 0.0003227490874791514; … ; 0.0; 0.003166366611555671;;], [0.0002437824094488339; 0.0007147366383763147; … ; 0.0; 0.005898702898659847;;], [0.0003694953412559654; 0.0003190110333445031; … ; 0.0; 0.0015567927901792289;;], [0.0002988470420785261; 0.000671515322579748; … ; 0.0; 0.003417976337485973;;]], [5, 6, 7, 8, 9, 10, 11, 12], 1)"
+       "LowLevelFEM.StressField([[-1.1102230246251565e-16; 4.2351647362715017e-16; … ; 0.0; -2.220446049250313e-16;;], [-1.6653345369377348e-16; 9.773457083703465e-17; … ; 0.0; 1.6653345369377348e-16;;], [5.551115123125783e-17; 3.257819027901155e-16; … ; 0.0; -7.771561172376096e-16;;], [-3.608224830031759e-16; -1.303127611160462e-16; … ; 0.0; 2.220446049250313e-16;;], [4.996003610813204e-16; 7.167201861382541e-16; … ; 0.0; 2.831068712794149e-15;;], [3.0531133177191805e-16; 9.773457083703465e-17; … ; 0.0; 3.1086244689504383e-15;;], [-3.3306690738754696e-16; -1.303127611160462e-16; … ; 0.0; 9.43689570931383e-16;;], [5.551115123125783e-17; -5.864074250222079e-16; … ; 0.0; 1.609823385706477e-15;;]], [5, 6, 7, 8, 9, 10, 11, 12], 1)"
       ]
      },
      "metadata": {},
@@ -283,7 +283,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 65,
+   "execution_count": 31,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -293,7 +293,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 66,
+   "execution_count": 32,
    "metadata": {},
    "outputs": [
     {
@@ -343,7 +343,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 67,
+   "execution_count": 33,
    "metadata": {},
    "outputs": [
     {
@@ -366,9 +366,9 @@
       " [2] getEntitiesForPhysicalName(name::String)\n",
       "   @ LowLevelFEM.gmsh.model ~/.julia/artifacts/48f681ec04b6ed4aad46a56b1e6dc4bf716585f1/lib/gmsh.jl:640\n",
       " [3] plotOnPath(problem::LowLevelFEM.Problem, pathName::String, field::Int32, points::Int64; step::Complex{Int64}, plot::Bool, name::String, visible::Bool)\n",
-      "   @ LowLevelFEM ~/Dokumentumok/GitHub/LowLevelFEM.jl/src/LowLevelFEM.jl:1826\n",
+      "   @ LowLevelFEM ~/Dokumentumok/GitHub/LowLevelFEM.jl/src/LowLevelFEM.jl:1828\n",
       " [4] top-level scope\n",
-      "   @ ~/Dokumentumok/GitHub/LowLevelFEM.jl/examples/jl_notebook_cell_df34fa98e69747e1a8f8a730347b8e2f_X16sZmlsZQ==.jl:1"
+      "   @ ~/Dokumentumok/GitHub/LowLevelFEM.jl/examples/jl_notebook_cell_df34fa98e69747e1a8f8a730347b8e2f_X20sZmlsZQ==.jl:1"
      ]
     }
    ],
@@ -380,7 +380,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 68,
+   "execution_count": 34,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -394,7 +394,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 69,
+   "execution_count": 35,
    "metadata": {},
    "outputs": [
     {
@@ -430,7 +430,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 70,
+   "execution_count": 36,
    "metadata": {
     "slideshow": {
      "slide_type": ""

examples/test/2D-1.geo

@@ -0,0 +1,22 @@
+//+
+SetFactory("OpenCASCADE");
+Rectangle(1) = {1, 0, 0, 10, 10, 0};
+
+
+//+
+Transfinite Curve {4, 2} = 100 Using Progression 1;
+//+
+Transfinite Curve {3, 1} = 100 Using Progression 1;
+//+
+Transfinite Surface {1};
+//+
+Recombine Surface {1};
+
+Mesh.ElementOrder = 1;
+Mesh 2;
+//+
+Physical Surface("body", 5) = {1};
+//+
+Physical Curve("supp", 6) = {1};
+//+
+Physical Curve("load", 7) = {3};

examples/test/2D-1.geo~

@@ -0,0 +1,17 @@
+//+
+SetFactory("OpenCASCADE");
+Rectangle(1) = {1, 0, 0, 10, 10, 0};
+
+
+//+
+Transfinite Curve {4, 3, 2, 1} = 10 Using Progression 1;
+//+
+Transfinite Surface {1};
+//+
+Recombine Surface {1};
+//+
+Physical Surface("body", 5) = {1};
+//+
+Physical Curve("supp", 6) = {1};
+//+
+Physical Curve("load", 7) = {3};