[수치해석]FEM을 이용한 열전달 해석

소개글

수치해석- FEM을 이용한 열전달 해석

포트란 코딩 소스 포함

목차

1. 문제에 대한 지배방정식과 경계조건을 설정하고 weakform 및 matrix eqn.을 유
2. program 에 대한 flowchart, 각 routine에 사용된 algorithm 에 대해 설명
4. 해석에 사용한 요소 수 및 time step의 변화에 따른 solution의 변화를 나타내고, 해석에 적합한 값을 제시5. 더 정확한 해석을 위한 방안

본문내용

2. program 에 대한 flowchart, 각 routine에 사용된 algorithm 에 대해 설명

flow chart

main routine에서는 heatflux가 일어나는 부분과 convection이 일어나는 부분의 element 개수와 size를 입력받고 그것으로 mesh를 만들고, array에 node-element connectivity를 저장하였다. subroutine(name:selMK)에서 Jacobian과 적분에 필요한 값들을 계산하고 각 element마다의 volume integral과 surface integral를 하였다. 그리고 다시 main routine에서 각 element에 해당하는 global node에 대응하는 위치에 값들을 입력한다. 즉, assemble한다. 이 부분에서 B.C가 적용되는 부분의 assemble도 한다. 그 후 main routine에서 assemble된 , , 과 초기 온도 값을 입력받은 후 와 를 계산한다. 구해진 와 를 이용하여 time increment를 증가시키고, 를 계산한다. 구해진 를 이용하여 를 구한다. 위의 과정을 iteration하면서, 지정한 accuracy 를 만족하는 steady한 값을 계산한다. accuracy는 지정된 3개의 점에서의 온도와 전체평균온도에서 이전 step과 온도차이가 0.1K이하로 되었을 때, steady 상태로 간주했다.

참고 자료

없음