공대생 도전 일지

2024.11.16 - [FEM/HOB] - HOB 기록 저장소 - HOB Shell Node 연결, 벽 연결까지 + Modal HOB 기록 저장소 - HOB Shell Node 연결, 벽 연결까지 + Modal이제 드디어 HOB에 끝이 보이네요작년에 "C++로 돈 벌어보자" 라고 시작한 CAE 프로그램 제작이 슬슬 마지막이 보입니다.Beam, Shell을 벽면에 접합하는 것으로 시작하여 Node 연결도 진행하고, Modal에yoonschallenge.tistory.com기존에 제작한 버전에서 좌표계 때문에 문제가 생긴 것이 있었습니다.제가 사용한 좌표계와 기존 제작자가 사용한 좌표계가 살짝 달라서 오류가 생겼고, 빔의 방향이 다를 때 떨어지는 버그가 있었습니다.제가 사용한 좌표계가 절대적으로 맞지 않다..

이제 드디어 HOB에 끝이 보이네요작년에 "C++로 돈 벌어보자" 라고 시작한 CAE 프로그램 제작이 슬슬 마지막이 보입니다.Beam, Shell을 벽면에 접합하는 것으로 시작하여 Node 연결도 진행하고, Modal에 대한 시각화까지 진행하였네요. 교수님이 제가 궁금할 때 마다 궁금증을 해결해 주셨으니 만들 수 있었지만 임금 체불은 많이 섭섭했어요.........그래도 일단 거의 완성했으니 막바지 버그만 조금 잡으면 될 것 같네요 일단 사진으로 좀 확인하겠습니다.이 사진은 차량 Modal 해석입니다.전체적인 부분은 못 보여주니 그냥 해석이 돌아갔다, 내가 차량 해석도 돌려 봤다 정도 기록으로 남겨놓네요이제 인공지능만 하고 차량, 자율주행이랑은 거의 끝이니... 서브 프레임 Modal 해석입니다.전문 지..

어제 코드를 정신없이 작성했고...오늘은 그에 대한 기록을 해놔야 기억을 할 수 있을 것 같네요일단 Beam 부터 시작하겠습니다. Beam은 가운데 있는 노드마다 n개의 자유도를 가집니다. n은 적을 땐 20개가 될 수도 있고, 100개가 될 수 도 있습니다.Beam의 단면 형상따라 모습이 바뀌는 것으로 알고 있습니다.이렇게 Node마다 n개의 자유도와 미분 자유도로 구성되어 있습니다.그럼 이제 빔의 특점 지점 즉 Nodal section의 3차원 변위는 이러한 방식으로 구할 수 있습니다.지정된 Psi함수에 edge정보와 scoord 좌표를 넣어주면 특정 위치의 변위를 구할 수 있습니다. 그럼 특정 지점 변위 두 개를 잡아서 같다고 두면 두 지점의 변위는 동시에 움직이게 됩니다이렇게 각각 변위를 뽑아서 ..

빔에는 3가지 좌표계가 있다.Global, Beam, sectional edge 좌표계가 존재한다.이 Beam 노드는 각각 n개의 자유도를 가지게 된다. 그럼 이 6개 노드의 빔은 아래와 같은 총 12 * n 자유도를 가지게 된다. 자유도를 미분한 것도 가지게 된다. 이제 여기서 쉘 좌표계를 정리해보자!!HOB input file에서 마지막 숫자 3개는 Beam 좌표계의 x 축 방향을 나타내는 단위백터이다.ex)       1,        1,        1,        3,      0.0,    500.0,      0.0        2,        1,        3,        4,      0.0,    500.0,      0.0        3,        1,        4..

HOB 이론에 기반한 C++ CAE 프로그램 제작 기록용입니다.modal 모드에 따른 형상 변화처음에 노드를 잘 못 맞춰서 오류가 있었는데 그 부분 수정하니 정상적으로 돌아갑니다.여긴 원통형 modal 해석인데 쉘이 사각 모양입니다.이건 원통형인데 쉘이 삼각형 모양입니다.사각 모양 쉘과 차이가 있어서 존재하는 matlab 코드로 비교해보았습니다.matlab 코드도 동일하게 변형되고, freq도 동일한 것을 볼 수 있습니다.이건 교수님께 보고드리면서 말씀드려보고, 수정할 것이 있으면 수정 진행하겠습니다.삼각쉘과 사각쉘이 섞여도 정상적으로 돌아가는 모습을 보입니다.마지막으로 삼각쉘과 사각쉘, 재료까지 섞인 원통형 modal 해석입니다.제대로 돌아가고, 오류가 없는 것을 볼 수 있습니다.여기서도 삼각쉘이 뭔가..

// matlabfunction [Ke, Me, TLG] = MITC4(Young, Poisson, thick, density, node, Vn) nNPE = 4; % # of nodes per element nDPN = 5; % # of DOFs per node nDIM = 3; % Problem dimension nGP = 2; % # of Gauss points % Vn = Make_Vn(node); [V1,V2] = Director(Vn(1,:)); % norn vec 구하기 TLG = [V1' V2' Vn(1,:)']; Ke = Shell_Stiffness(Young, Poisson, thick, node, Vn); ..

node - 각 노드들 위치 정보etable - 쉘을 구성하는 노드 번호들 fnode- 힘이 들어간 노드 번호 fix_node - boundary condition 노드 뽑기 물성치 뽑기 shell 자유도(Ndof)가 5이다 -> 총 자유도(tdof)는 5 * 노드 수(Nn) Vn - 각 노드 들의 norm 벡터 뽑기 근데 왜 쉘이 아니라 노드들의 벡터를 뽑지 ?K matrix - sparse 메트릭스로 총 자유도의 열과, 행으로 이루어져있다.반복문 마다 그 etable에 맞는 노드들과 노말 벡터들 뽑아서 MITC4에 넣기shell_stiffness 물성치 구하기! i, j, k 가 전부 1,2로 for 문 반복합니다.가우시안22 포인트 ?? - 1이면 -0.577- 2이면 +0.577stain_Disp..

20200323노드 당 자유도 * 노드 수 = K matrix 행과 열의 길이1. K element 만들기2. assemble 해서 k matrix 만들기3. boundary condition 4. kd = F 풀기 5. ub, Fa 풀기!!문제 푸는 Matlab 코드입니다.% step 1 : element stiffness matrices generationk = 200ek1 = [k , -k; -k, k];ek2 = ek1ek3 = ek2ek4 = ek3% step 2 : assembletdof = 5K = zeros(tdof);K([1,2],[1,2]) = K([1,2],[1,2]) + ek1 ;K([2,3],[2,3]) = K([2,3],[2,3]) + ek2 ;K([3,4],[3,4]) = K(..

이전에 받았던 일을 다시 진행하게 되어 조금만 정리 하겠습니다.. 20200316FEM - 컴퓨터를 이용해 구조 해석(고체, 동)을 진행하여 디자인을 하는 과목 PDE  - 수학적으로는 전부 미분 방정식을 푸는 문제이다   구조에 힘을 주면 변형(deformed shape)이 일어난다- 탄성 영역에서만 공부한다.변형체는 스프링으로 볼 수 있다. stiffness Method == Spring elementk = spring constant = spring stiffness = 강성 u = displacementF = Force F = ku 구조 하나하나를 쪼갠 것 = element절점들 = nodes 이렇게 하나하나 나누게 되면 백터와 메트릭스로 변하게 된다. 스프링!!!K == stiffnessku =..