Bridge example

gpexamples/gpbridge/bridge.pro

@@ -5,16 +5,16 @@ gpinit =
   mesh =
   {
     type = gmsh;
-    file = meshes/beam_fine2.msh;
+    file = meshes/bridge-q4-r2.msh;
   };
 
   coarseMesh =
   {
     type = gmsh;
-    file = meshes/beam_coarse.msh;
+    file = meshes/bridge-q4-r0.msh;
   };
 
-  nodeGroups = [ lb, rb, tm ];
+  nodeGroups = [ lb, rb, m1, m2 ];
 
   lb =
   {
@@ -28,17 +28,23 @@ gpinit =
     ytype = min;
   };
 
-  tm =
+  m1 =
   {
-    ytype = max;
-    xtype = mid;
+    xtype = 3;
+    ytype = min;
+  };
+
+  m2 =
+  {
+    xtype = 7;
+    ytype = min;
   };
 };
 
 gpsolver =
 {
   type = GPSolver;
-  nsample = 3;
+  nsample = 100;
   seed = 0;
   tables = [ stress, strain ];
 };
@@ -65,12 +71,10 @@ model =
       nu = 0.2;
     };
 
-    thickness = 0.2;
-
     shape =
     {
-      type = Triangle3;
-      intScheme = Gauss1;
+      type = Quad4;
+      intScheme = Gauss4;
     };
   };
 
@@ -81,12 +85,12 @@ model =
     elements = all;
 
     dofs   = [ dy ];
-    values = [ -0.2 ];
+    values = [ -10. ];
 
     shape =
     {
-      type = Triangle3;
-      intScheme = Gauss1;
+      type = Quad4;
+      intScheme = Gauss4;
     };
   };
 
@@ -108,8 +112,8 @@ model =
 
     shape =
     {
-      type = Triangle3;
-      intScheme = Gauss1;
+      type = Quad4;
+      intScheme = Gauss4;
     };
   };
 
@@ -117,15 +121,8 @@ model =
   {
     type = Dirichlet;
 
-    groups = [ lb, lb, rb ];
-    dofs   = [ dx, dy, dy ];
-    values = [ 0., 0., 0. ];
+    groups = [ lb, lb, m1, m2, rb ];
+    dofs   = [ dx, dy, dy, dy, dy ];
+    values = [ 0., 0., 0., 0., 0. ];
   };
 };
-
-vtkout =
-{
-  file = results;
-  type = VTKOut;
-  tables = [ stress, strain ];
-};