WheelDPSimplified-Error

[sumaya]

Hello,

I am looking for some suggestions about how to fix my error for the script wheelDPsimplified. I have changed the script based on my research requirement. My terrain bed is made of glass beads of 3 mm radius, which is very small. The beads are a clump of two and purposely, I have overlapped them slightly. Now the volume is very small compare to the values provided in the demo script. My IXX, IYY, IZZ are also very small  and I am not currently scaling them. There was issue with the bin size but I have temporarily fixed it with the following lines. However, I am currently getting unexpected velocity that's crashing the entire system. I have played with the terrain properties, step size and HPC sampler spacing but nothing seems to fix this problem. Can you please suggest what could be possibly happening?

The error message :

DEME_ERROR(
            "System max velocity is %.7g, exceeded max allowance (%.7g).\nIf this velocity is not abnormal and you "
            "want to increase this allowance, use SetErrorOutVelocity before initializing simulation.\nOtherwise, the "
            "simulation may have diverged and relaxing the physics may help, such as decreasing the step size and "
            "modifying material properties.\nIf this happens at the start of simulation, check if there are initial "
            "penetrations, a.k.a. elements initialized inside walls.",

 Parameters Output  : 

Current number of clumps: 4216

Number of total active devices: 1
User-specified X-dimension range: [-0.45, 0.45]
User-specified Y-dimension range: [-0.102, 0.102]
User-specified Z-dimension range: [-0.155, 0.155]
User-specified dimensions should NOT be larger than the following simulation world.
The dimension of the simulation world: 1.9584000110626221, 0.24480000138282776, 0.48960000276565552
Simulation world X range: [-0.54, 1.4184]
Simulation world Y range: [-0.1224, 0.1224]
Simulation world Z range: [-0.186, 0.3036]
The length unit in this simulation is: 3.5623088680269976e-12
The edge length of a voxel: 2.3345947397501732e-07
The initial time step size: 5e-06
The initial edge length of a bin: 0.012287999153137209
The initial number of bins: 128000
The total number of clumps: 4216
The combined number of component spheres: 8432
The total number of analytical objects: 2
The total number of meshes: 1
Grand total number of owners: 4219
The number of material types: 2
History-based Hertzian contact model is in use.
The solver to set to adaptively change the contact margin size. 

The script crashes shortly after.

Terrain Properties and Co-ordinates :  

x            y             z      r
0.0015 0.0012   0    0.003   

0.0045   0.0028 0    0.003

 Glass beads : auto mat_type_terrain = DEMSim.LoadMaterial({{"E", 6.5e10}, {"nu", 0.245}, {"CoR", 0.926}, {"mu", 0.16}, {"Crr", 0.01}});

// Define the terrain particle templates
   // Calculate its mass and MOI
    float terrain_density = 2500;
    float volume1 = 1.99e-07 ;
    float mass1 = terrain_density * volume1;
    float3 MOI1 = make_float3(1.431e-12, 9.058e-13, 1.282e-12) * terrain_density;

HCPSampler sampler(0.025f);
    float sample_halfheight = 0.124;
    float sample_halfwidth_x = 0.42;
    float sample_halfwidth_y = 0.098;
    float offset_z = bottom + sample_halfheight + 0.01f;

[Ruochun Zhang]

Sumaya,

The test scene you are targeting appears to be perfectly doable, but in order for us to help you, you need to present two things:

1. The error message you gave seems to be the corresponding part of the code, rather than what got output to the screen... Please provide the actual error message.

2. Please provide a visual rendering of the simulation scene, before the crash. Even if it crashes right at the beginning of the simulation, you can still visualize the initial configuration. This is important and without visuals, anything we can say is speculative.

Thank you,

Ruochun

[sumaya]

Hi Ruochun,

Here is the error message(attached). Secondly, no vtk files were generated. I am running the script in debug mode in VS. The script crashed with the last message outputting frame=0(image attached). Thank you in advance.

Sumaiya

[sumaya]

This is the final exit message.

[sumaya]

Hi Ruochun, 

Sorry for the previous email. Although the error stays the same. However, I was able to get the one vtk file that was generated. From the visuals, it looks like the wheel is slightly bigger than the soil. I am not sure if this is the reason that the unexpected velocity error is happening. Also, I have changed the wheel dimensions specific to my research. However, I have used already available viper wheel .obj  file created by chrono. Do you suggest I replace that .obj file based on my research specifc wheel?  Below is the screenshot and csv file attached. 

[Ruochun Zhang]

Sumaya,

Somehow the message you paste is still not the final one, and it appears to have stopped before outputting the final metrics due to being stopped by the debug console or something. But it matters not, usually these kinds of problems are caused by initial penetrations anyway, as you would've seen if the error message were output to screen. The rendering you show is insufficient to show anything, as the spheres are not scaled by their real size. Please render it with the Scale Array being the radius, and the Scale Factor being 2. After that we may be able to see which part of the configuration is problematic.

Whichever obj file you use, it has to make physical sense, a.k.a. no initial penetration and has good mesh quality.

Thank you,

Ruochun

[sumaya]

Hi Ruochun, 

Below is the message that shows up in my debug console(Details are posted in the github link). I have been trying to fix this for a while thinking the dimensions that I have entered is incorrect. However, I have done the reverse calculations multiple times to verify the soil bed dimensions and the halfwidth for my own personalized soil bed and I get the same error message. I have also uploaded the WheelDPSimplified.cpp file as well as the csv template with the read me file, in case you want to generate the same error message in your desktop. Here is the github repo link : https://github.com/abedin17/Chrono-DEM-Engine_WheelDPSimplified. Thank you for your help.

'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Program Files\NVIDIA GPU Computing Toolkit\CUDA\v12.8\bin\nvrtc64_120_0.dll'. Module was built without symbols.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\user32.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\win32u.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\gdi32.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\gdi32full.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\msvcp_win.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\ucrtbase.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\ole32.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\combase.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\bcryptprimitives.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\shell32.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\cfgmgr32.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\SHCore.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\windows.storage.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\profapi.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\powrprof.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\umpdc.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\shlwapi.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\kernel.appcore.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\cryptsp.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\nvcuda.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Users\umroot\chrono_DEM-Engine\build\bin\Debug\DEMERuntimeDataHelper.dll'. Symbols loaded.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\msvcp140d.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\vcruntime140d.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\vcruntime140_1d.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\ucrtbased.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\cryptbase.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\imm32.dll'.
The thread 0x4490 has exited with code 0 (0x0).
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\DXCore.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\DriverStore\FileRepository\nvddi.inf_amd64_b26dcb0ba9fe7785\nvcuda64.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\version.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\msasn1.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\cryptnet.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\crypt32.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\wldp.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\wintrust.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\drvstore.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\devobj.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\nvapi64.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\setupapi.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\bcrypt.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\oleaut32.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\DriverStore\FileRepository\nvddi.inf_amd64_b26dcb0ba9fe7785\nvdxgdmal64.dll'.
The thread 0x4f10 has exited with code 0 (0x0).
The thread 0x4f18 has exited with code 0 (0x0).
The thread 0x2b4c has exited with code 0 (0x0).
The thread 0x457c has exited with code 0 (0x0).
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Program Files\NVIDIA GPU Computing Toolkit\CUDA\v12.8\bin\nvrtc-builtins64_128.dll'. Module was built without symbols.
The thread 0x4540 has exited with code 0 (0x0).
The thread 0x53d8 has exited with code 0 (0x0).
The thread 0x4774 has exited with code 0 (0x0).
The thread 0x384c has exited with code 0 (0x0).
The thread 0x4ed0 has exited with code 0 (0x0).
Exception thrown at 0x00007FF9343A3FA9 in DEMdemo_WheelDPSimplified.exe: Microsoft C++ exception: std::runtime_error at memory location 0x000000FE6FAFF030.

Debug Error!

Program: ...rono_DEM-Engine\build\bin\Debug\DEMdemo_WheelDPSimplified.exe

abort() has been called

(Press Retry to debug the application)
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\uxtheme.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\msctf.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\TextInputFramework.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\CoreUIComponents.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\CoreMessaging.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\ntmarta.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\WinTypes.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\WinTypes.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Unloaded 'C:\Windows\System32\WinTypes.dll'
'DEMdemo_WheelDPSimplified.exe' (Win32): Loaded 'C:\Windows\System32\WinTypes.dll'.
'DEMdemo_WheelDPSimplified.exe' (Win32): Unloaded 'C:\Windows\System32\WinTypes.dll'
The thread 0x575c has exited with code 0 (0x0).
The thread 0x5578 has exited with code 0 (0x0).
The thread 0x44d4 has exited with code 0 (0x0).
The thread 0x54b4 has exited with code 0 (0x0).
The thread 0x55c8 has exited with code 0 (0x0).
The thread 0x48d8 has exited with code 0 (0x0).

Sumaiya

[Ruochun Zhang]

Hi Sumaiya,

Instead of reverse calculations, I suggest making a rendering first in debugging. Attached is what I got by rendering your initial configuration. Obviously, there's a huge initial overlap between the mesh and particles (the mesh is also largely out-of-the-simulation-world, by the way), and the physics cannot be right. For a start, I suggest removing the meshed wheel from your simulation, and testing how your particle bed settles. Then add back the mesh with a proper size and location.

Thank you,

Ruochun