Ask for 3Dversion MPM

zhulinshenchu · 2020 年1 月 17 日 16:10

Can you provide me with a simplified 3D version of MPM method code, just like taichi88 lines. I really need it, Thank you so much.

yuanming · 2020 年1 月 17 日 17:17

Please checkout this file: https://github.com/taichi-dev/taichi_elements/blob/master/demo_3d.py

zhulinshenchu · 2020 年1 月 17 日 20:10

Thank you very much for your reply. I ran the above code with the following questions:

Can you help me to figure out the cause of the problem? Thank you very much again

yuanming · 2020 年1 月 17 日 21:45

Which Taichi version are you on? v0.3.21 should work.

zhulinshenchu · 2020 年1 月 18 日 07:54

It works now after I changing it to version 3.2.1. Thank you so much

archibate · 2020 年12 月 21 日 05:47

Now available in examples/mpm3d.py, a 3D version of mpm88. Run ti example mpm3d to have fun

binYang320 · 2022 年10 月 20 日 02:24

Now where is the file

Proton · 2022 年10 月 20 日 03:37

github.com

taichi-dev/taichi/blob/master/python/taichi/examples/simulation/mpm3d.py

export_file = ''  # use '/tmp/mpm3d.ply' for exporting result to disk

import numpy as np

import taichi as ti

ti.init(arch=ti.gpu)

#dim, n_grid, steps, dt = 2, 128, 20, 2e-4
#dim, n_grid, steps, dt = 2, 256, 32, 1e-4
dim, n_grid, steps, dt = 3, 32, 25, 4e-4
#dim, n_grid, steps, dt = 3, 64, 25, 2e-4
#dim, n_grid, steps, dt = 3, 128, 25, 8e-5

n_particles = n_grid**dim // 2**(dim - 1)
dx = 1 / n_grid

p_rho = 1
p_vol = (dx * 0.5)**2
p_mass = p_vol * p_rho

This file has been truncated. show original

You can also run it directly via ti example (at index #36) in command line.

Proton · 2022 年10 月 21 日 06:44

Could you please open up a new thread for your question? This seems off-topic here.