Force Chain 그리기

먼저 필요한 데이터

dump dump_forceChainFabricTensor all local/gran/vtk ${screen_nstep} post/forcechain*.vtp cpgl

 

Atom 먼저 Clip으로 잘라서 단면을 보여주고,

 

Glyph Scaling도 여기서 한다.

 

 

그 다음, cpgl로 dump한 vtp 파일을 PARAVIEW에서 load하면, 다음의 cell data를 포함하는 것을 알 수 있다.

 

우리에게 필요한 것은 force_normal의 합력

 

다음으로 Calculator 필터를 사용하여 force_normal의 합력을 계산한다. 계산 결과는 magF라는 Cell Data로 저장된다.

 

Force Chain을 Tube형태로 그리기 위해서는 magF가 Cell Data에서 Point Data로 변환되어야 한다.

그러기 위해서 CellDatatoPointData 라는 필터를 사용한다.

 

 

 

너무 많은 Tube가 생기면 지저분해서 보기 힘드니까, 주된 Force Chain 만 보기위해서는 또 필터를 써야 된다.

이때 사용되는 필터는 Threshold 이다.

 

 

이후 Tube 만들기 위한 Extract Surface 필터를 한번 써 주고,

 

Tube 필터 적용

여기서 빡치는게, Paraview는 매 타임스텝마다 max, min값을 가지고 max radius, min radius로 그린다. 절대적인 값으로 그리지 않는 다는 이야기다.

radius의 min은 radius, max는 radius X radius Factor이다.

paraview의 병신같은 기능으로 몇일을 고생하지 갑자기 빡치네.

 

로그스케일이 좋다니까,

 

 

 

잘라서 보려면 Filter중 Clip, Box로

 

이렇게 나온다

 

 

필터트리를 정리하면

 

 

 

세상에 쉬운일이 없다.