DO NOT MERGE; test with KW

ORNL · Jul 23, 2024 · b29d902 · b29d902
1 parent c4b7c9c
commit b29d902
Showing 1 changed file with 14 additions and 8 deletions.
diff --git a/examples/mechanics/kalthoff_winkler.cpp b/examples/mechanics/kalthoff_winkler.cpp
@@ -76,23 +76,26 @@ void kalthoffWinklerExample( const std::string filename )
     CabanaPD::Prenotch<2> prenotch( v1, v2, notch_positions );
 
     // ====================================================
-    //                    Force model
+    //                    Force model options
     // ====================================================
-    using model_type = CabanaPD::ForceModel<CabanaPD::PMB, CabanaPD::Fracture>;
-    model_type force_model( delta, K, G0 );
-    // using model_type =
-    //     CabanaPD::ForceModel<CabanaPD::LPS, CabanaPD::Fracture>;
-    // model_type force_model( delta, K, G, G0 );
+    using model_type = CabanaPD::PMB;
+    using thermal_type = CabanaPD::TemperatureIndependent;
 
     // ====================================================
     //                 Particle generation
     // ====================================================
     // Does not set displacements, velocities, etc.
     auto particles = std::make_shared<
-        CabanaPD::Particles<memory_space, typename model_type::base_model,
-                            typename model_type::thermal_type>>(
+        CabanaPD::Particles<memory_space, model_type, thermal_type>>(
         exec_space(), low_corner, high_corner, num_cells, halo_width );
 
+    std::array<double, 2> same_d = { delta, delta };
+    std::array<double, 2> fake_K = { K, K * 1.2 };
+    std::array<double, 2> fake_G0 = { G0, G0 * 1.8 };
+    auto force_model =
+        createForceModel( model_type{}, CabanaPD::Fracture{}, thermal_type{},
+                          *particles, same_d, fake_K, fake_G0 );
+
     // ====================================================
     //                Boundary conditions
     // ====================================================
@@ -112,6 +115,7 @@ void kalthoffWinklerExample( const std::string filename )
     auto x = particles->sliceReferencePosition();
     auto v = particles->sliceVelocity();
     auto f = particles->sliceForce();
+    auto type = particles->sliceType();
 
     double v0 = inputs["impactor_velocity"];
     auto init_functor = KOKKOS_LAMBDA( const int pid )
@@ -122,6 +126,8 @@ void kalthoffWinklerExample( const std::string filename )
         if ( x( pid, 1 ) > y_prenotch1 && x( pid, 1 ) < y_prenotch2 &&
              x( pid, 0 ) < -0.5 * height + dx )
             v( pid, 0 ) = v0;
+        if ( x( pid, 1 ) > 0.0 )
+            type( pid ) = 1;
     };
     particles->updateParticles( exec_space{}, init_functor );