예전에 GranWall에 작용하는 Force를 뽑아내는 방법에 대해 포스팅 한 적이 있다.
몇번 사용하다 보니, GranWall로 한정되는 것이 아니라 mesh/surface/stress 로 정의된 Wall이 가능하다.
주의! : mesh/surface 로 정의된 Wall에서는 information을 뽑을 수 없다 (단정하지는 못한다. 아직 발견하지 못했을 뿐이다)
방법은 mesh/surface/stress로 wall을 정의하고 stress on 옵션을 활성화 한다.
fix box_sieve all mesh/surface/stress file sieve_1.stl type 2 scale 1.0 stress on
참고로 본 예에서 위의 box_sieve mesh는 위로 움직이며, atom이 통과할 수 있는 sieve로 설정하였다.
fix sieve_D10 all massflow/mesh/sieve mesh box_sieve point_at_outlet 0 0 -1000 count once delete_atoms no sieveSize 4e-3 sieveSpacing 4.5e-3 sieveStiffness 1e5 sieveDamping 0.01
fix move_sieve all move/mesh mesh box_sieve linear 0 0 0.008
이후, z-position을 뽑아내려면 아래와 같이 하면 된다. 또는 print, region을 설정하거나.
각 어레이가 나타내는 값은 아래를 참조
variable sieve_position_z equal f_box_sieve[9] # Z-position of box sieve
print "${sieve_position_z}"
region above_sieve block -0.033 0.033 -0.033 0.033 ${sieve_position_z} 0.06
'LIGGGHTS' 카테고리의 다른 글
unfix move/mesh (0) | 2024.01.02 |
---|---|
mesh/surface/stress 기준 coordinate (0) | 2023.11.11 |
servo controls (0) | 2023.10.19 |
Contact Model Optimization (0) | 2023.08.30 |
Log 파일의 이해 (1) | 2023.08.27 |