|
| 1 | +use slosh_testbed2d::{RapierData, slosh}; |
| 2 | + |
| 3 | +use nalgebra::{point, vector}; |
| 4 | +use rapier2d::prelude::{ColliderBuilder, RigidBodyBuilder}; |
| 5 | +use slang_hal::backend::WebGpu; |
| 6 | +use slosh::{ |
| 7 | + pipeline::MpmData, |
| 8 | + solver::{Particle, SimulationParams}, |
| 9 | +}; |
| 10 | +use slosh_testbed2d::{AppState, PhysicsContext}; |
| 11 | +use slosh2d::solver::{GpuBoundaryCondition, ParticleModel}; |
| 12 | + |
| 13 | +#[allow(dead_code)] |
| 14 | +fn main() { |
| 15 | + panic!("Run the `testbed2` example instead."); |
| 16 | +} |
| 17 | + |
| 18 | +pub fn beam_demo(backend: &WebGpu, app_state: &mut AppState) -> PhysicsContext { |
| 19 | + let mut rapier_data = RapierData::default(); |
| 20 | + |
| 21 | + let width = 10.0; |
| 22 | + let height = 2.0; |
| 23 | + let fixed_part = 1.0; |
| 24 | + let cell_width = 0.2; |
| 25 | + let particle_per_cell_dim = 2; |
| 26 | + let young_modulus = 1.0e8; |
| 27 | + let poisson_ratio = 0.3; |
| 28 | + |
| 29 | + let diameter = cell_width / particle_per_cell_dim as f32; |
| 30 | + let ni = ((width + fixed_part) / diameter).ceil() as usize; |
| 31 | + let nj = (height / diameter).ceil() as usize; |
| 32 | + |
| 33 | + let mut particles = vec![]; |
| 34 | + for i in 0..ni { |
| 35 | + for j in 0..nj { |
| 36 | + let position = point![i as f32, j as f32] * diameter; |
| 37 | + let density = 1000.0; |
| 38 | + let radius = diameter / 2.0; |
| 39 | + let model = ParticleModel::elastic_neo_hookean(young_modulus, poisson_ratio); |
| 40 | + particles.push(Particle::new(position, radius, density, model)); |
| 41 | + } |
| 42 | + } |
| 43 | + |
| 44 | + if !app_state.restarting { |
| 45 | + app_state.min_num_substeps = 150; |
| 46 | + app_state.max_num_substeps = 150; |
| 47 | + app_state.gravity_factor = 1.0; |
| 48 | + }; |
| 49 | + |
| 50 | + let params = SimulationParams { |
| 51 | + gravity: vector![0.0, -9.81] * app_state.gravity_factor, |
| 52 | + dt: 1.0 / 60.0, |
| 53 | + padding: 0.0, |
| 54 | + }; |
| 55 | + |
| 56 | + let rb = RigidBodyBuilder::fixed() |
| 57 | + .translation(vector![0.0, height / 2.0]) |
| 58 | + .build(); |
| 59 | + let rb_handle = rapier_data.bodies.insert(rb); |
| 60 | + let co = ColliderBuilder::cuboid(fixed_part, height); |
| 61 | + let ground = rapier_data |
| 62 | + .colliders |
| 63 | + .insert_with_parent(co, rb_handle, &mut rapier_data.bodies); |
| 64 | + |
| 65 | + let data = MpmData::new( |
| 66 | + backend, |
| 67 | + params, |
| 68 | + &particles, |
| 69 | + &rapier_data.bodies, |
| 70 | + &rapier_data.colliders, |
| 71 | + &[(ground, GpuBoundaryCondition::stick())], |
| 72 | + cell_width, |
| 73 | + 30_000, |
| 74 | + ) |
| 75 | + .unwrap(); |
| 76 | + PhysicsContext { |
| 77 | + data, |
| 78 | + rapier_data, |
| 79 | + callbacks: vec![], |
| 80 | + hooks_state: None, |
| 81 | + } |
| 82 | +} |
0 commit comments