본문 바로가기

LIGGGHTS

Restart 이후 Atom property 변경 https://www.cfdem.com/media/DEM/docu/set.html 더보기
Pair granular requires ghost atoms store velocity SuperQuadric 또는 Multisphere를 사용한 restart 파일을 읽어서 다시 실행할 경우 종종 보게 되는 에러 메세지이다. 이를 해결하기 위해서는 아래의 명령어를 read_restart 이전에 실행하여야 한다. communicate single vel yes 그럼 20,000 더보기
Change box LIGGGHTS의 Change box기능은 Project Space의 크기를 변화시킬 수 있는 명령이다. Change box를 통해 CPU가 담당하는 영역 배분을 균등하게 하여 MPI 해석의 효율성을 늘릴 수 있다. 그런데 LIGGHTS의 Change box는 stl mesh를 사용하는 경우, 뜻하지 않는 에러를 발생시키기도 한다. 예를 들어, 아래와 같은 에러가 가끔 나오는데, 내가 불러들인 stl파일에 기반한 mesh가 region 안에 들어와도 저런 에러가 종종 발생한다. 왜 그럴까? ERROR: Mesh (id cad2): Mesh elements have been lost / left the domain. Please use 'boundary m m m' or scale/translate/rot.. 더보기
libboost libraries LIGGGHTS의 SuperQuadrics Simulation 최적화를 위해서는 C++ libboost library를 설치해서 컴파일 해야 한다고 한다. https://www.cfdem.com/media/DEM/docu/superquadric_simulations.html How-To for “SUPERQUADRIC” simulations in LIGGGHTS(R) — LIGGGHTS v3.X documentation Docs » How-To for “SUPERQUADRIC” simulations in LIGGGHTS(R) View page source Website Commands How-To for “SUPERQUADRIC” simulations in LIGGGHTS(R) Description: T.. 더보기
fix particletemplate/superquadric LIGGGHTS의 superquadric particle template 생성을 위한 공식 명령어 reference에 오류가 있다. 아래와 같이 Example이 적혀 있으나, fix pts1 all particletemplate/superquadric 123457 atom_type 1 density constant 2500 shape constant 0.0015 0.0025 0.0035 blockiness 4.0 4.0 실제는 blockiness 값도 아래와 같이 constant 지정 후 입력하여야 한다. fix pts1 all particletemplate/superquadric 123457 atom_type 1 density constant 2500 shape constant 0.0015 0.0025.. 더보기
unfix move/mesh LIGGGHTS에서 동일 atom (wall)에 다수의 move/mesh fix가 지정되는 경우, unfix를 시키기 위해서는 reverse order로 하라고 한다. 예를 들어, 이렇게 하면, # Define moving of funnel fix funnel_xmove_1 all move/mesh mesh funnel_d10 viblin axis 1 0 0 order 1 amplitude 0.005 phase -1.57 period 0.1 fix funnel_ymove_1 all move/mesh mesh funnel_d10 viblin axis 0 1 0 order 1 amplitude 0.005 phase -1.57 period 1.0 fix funnel_zmove all move/mesh mesh.. 더보기
mesh/surface/stress 기준 coordinate https://nmge.tistory.com/87 WALL의 Information 획득하기 https://nmge.tistory.com/66 예전에 GranWall에 작용하는 Force를 뽑아내는 방법에 대해 포스팅 한 적이 있다. 몇번 사용하다 보니, GranWall로 한정되는 것이 아니라 mesh/surface/stress 로 정의된 Wall이 가능하다. nmge.tistory.com https://nmge.tistory.com/49 SERVO VARIABLE This mesh module stores a global vector with 9 components for access by various output commands. The first 3 components are equal to the.. 더보기
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 m.. 더보기