WeiyiVirtonomy
WeiyiVirtonomy
The previous PR was closed due to the incorrect merge. ## TODO: - Add test cases back - Check sliding under different resolutions of solid and shell
A 2d test case of multi-time resolution contact problem ## Description The left brick represents a catheter with higher stiffness, while the right brick represents the softer arterial wall. The...
The split cell linked list is replaced by a function called `particle_for_split`, where the cells are accessed by a stride of 3 in the indices. TODO: Add sequential version Add...
@Xiangyu-Hu Hello, Prof. Hu, I still have some questions about the update of CLL we discussed last time. I remember the last time you said that to get rid of...
**Describe the bug** The inflow region has velocity and pressure singularities if the fluid particles are relaxed. E.g. T-shaped pipe flow: This is within the inflow region, where the velocity...
I added a function `Real ContactStiffness(size_t index_i)` in `Solid`. For composite material, it will return `composite_materials_[(*material_id_)[index_i]]->ContactStiffness();`. In `RepulsionForce`, `solid_material->ContactStiffness()` is replaced by `solid_material->ContactStiffness(index_i)`. I think the contact impedance of self-contact...
Aim: test the behavior of solid multi-resolution simulation - [x] Use minimum smoothing length in time step size - [ ] Use smoothing_length_ / h_ratio in `Integration1stHalf` (maybe not necessary)...
A study of stiff-soft contact. The stiffer body has stronger compressibility, while the softer one is nearly incompressible. Aim: study if it is possible to replace Young's modulus of the...
The function `CalculateAllLambda()` defines the compliance matrix as: ``` Matd Compliance = Matd::Zero(); Compliance.col(0) = Vecd(1 / E_[0], -poisson_[0] / E_[1], -poisson_[1] / E_[2]); Compliance.col(1) = Vecd(-poisson_[0] / E_[0], 1...
The relaxation using `RelaxationStepLevelSetCorrectionInner` produces weird particle distribution Lattice:  After relaxation:  Gaps:  `RelaxationStepInner` doesn't have this problem: There might be a...