본문 바로가기

LIGGGHTS

WALL의 Information 획득하기

https://nmge.tistory.com/66

 

예전에 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을 설정하거나.

각 어레이가 나타내는 값은 아래를 참조

https://nmge.tistory.com/49

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