CST요소를 이용한 캔틸레버보 FEM 해석
공학공부 / 2010. 11. 17. 17:57
어젯밤 꼬박 새서 숙제를 완성했다.
이 "유한요소해석입문"이란 FEM수업은 공대 이해성 교수의 수업이었다고 하는데 지금은 박원석이라는 강사가 가르치고 있다. 이해성 교수 수업을 들어본 것은 아니지만 난 오히려 지금 강사님 수업도 매우 마음에 든다. 명료하고 재미있게 가르친다. 물론 그렇다고 내 수업태도가 좋은건 아니다. 한시간에 다섯번정도씩 꾸벅꾸벅 존다. 미안해 죽겠다.
그리고 중간고사는 거의 꼴지한 것 같다. 반타작했다.
뭐 그래도 수업 자체에서 얻고 있는건 굉장히 많고 만족한다.
어제 낮 네시쯤이면 코딩 다하고 레포트를 여섯시쯤까지 쓸 줄 알았더니 왠걸.... 여기저기 숨어있는 오류들이 어찌나 많은지..... 컴파일 하는 과정에서 코딩 잘못했다고 에러가 떠주면 좋으련만, 내가 만든 프로그램은 아무런 얘기도 없이 결과까지 띡 내놓는데 그 결과가 맞는 것 같으면서도 안 맞더라는....
시간에는 쫓기고 잠은 자고 싶고... 점점 코드는 걸레가 되어가고 코딩하기 싫은건 인터넷 여기저기서 긁어오고 난리도 아니었다.
그래도 나름의 성과는 있었다고 생각하는게 조금은 포트란에 친숙해졌다는 느낌이 든다는 것이다.
"hello world"도 찍을 줄 모르던 나에게 프로그램을 맨 처음으로 가르쳐준 xeraph선생께서는 OOP(객체지향프로그래밍)를 강조하며, "클래스는 붕어빵 틀, 인스턴스는 붕어빵", 그리고 "한 함수에서는 한가지 기능만 하게 만들라"는 명언을 남기셨다. 물론 내가 잘못 이해해서 왜곡해서 기억하고 있는지도 모른다.
맨처음 배운게 C++이었고, 그 이후에는 JAVA로 갈아탔다. 나의 두번째 스승님인 연구실의 gonny가 java를 하고 있었기에... 이사람도 OOP 매니아. 시뮬레이션 연구하는 사람이다보니까 그렇겠지...
어쨌든 하고 싶은 이야기는 난 순차적인 프로그래밍을 해본적이 없다.
어떤 게임프로그래밍 책에서는 어설프게 OOP따라하려고 애쓰지 말고 그냥 촌스러워도 순차적 프로그래밍하라고 그게 더 빨리 만든다고 조언되어 있는것까지 봤지만, 어디까지나 그건 순차적 프로그래밍으로 몇십년을 먹고 살았던 사람들 이야기. 난 오히려 그렇게 체계적으로 처음부터 쫘악 파악하고 쓰는 프로그래밍하기엔 머리가 너무 후달린단 말이다. 그러한 연유로 나한테 포트란은 정말 짜증나는 언어였다. 포트란 책을 봐도 갑갑하긴 마찬가지. 포트란이 손에 익은 사람들이라면 어쨌든 OOP이전 세대들일 테니까..
근데 구글링을 열심히 하다보니 포트란으로 클래스 만들고 어쩌고 하는 방법을 알게 되었다. 하여간 포트란에 대한 불평불만은 다음에 또 쓰기로 하자. 어차피 엄청 긴 포스팅이 될테니까.
숙제는 다음과 같다.
세가지 하중 조건에서 저 구조물이 어떻게 변형될까 풀어내는거다. 초록색일 때, 빨간색일때, 주황색일때.
그리고 이걸 CST요소로 풀어내는거였다.
문제에 대해 자세히 알고 싶다면 이해성 교수님의 연구실 홈페이지에 가서 찾는게 좋겠지.
몹시 창피한 레포트이지만 어쨌든 이렇게 정리 해놓지 않으면 굴러다니다가 결국 어디갔는지 모르게 될테니 이렇게 올린다. 포트란으로 작성한 코드 역시. 물론 내가 작성한 코드는 여기저기서 긁어온 코드들도 섞여있는 걸레같은 부끄러운 코드다. 게다가 시간에 쫓겨서 여기저기 테스트를 위한 코드들도 남아있다. 정리할 시간이 없었다. 프로그램을 완성한게 새벽 다섯시쯤이었고, 보고서 쓰고 어쩌고 하니 일곱시 반이었다.
포트란으로 작성한 코드는 다음과 같다.
Beam.f90
결과는 이렇다. 이건 데이터 받아서 엑셀에 넣어서 그린거다. 어쩌다보니 이 그래프는 레포트에 첨부 못했네.
CASE 2그림은 변위가 너무 작아서 30배 뻥튀기 해서 과장해서 그린거다.
여기 안올려놓으면 보나마나 잃어버릴테니 PDF 파일도 여기에 걸어놓자.
2010fem_hw6.pdf'공학공부' 카테고리의 다른 글
| 공대 FEM 숙제 fortran을 이용한 근사식 구하기 (0) | 2010.09.29 |
|---|