공대생 도전 일지

어제 코드를 정신없이 작성했고...오늘은 그에 대한 기록을 해놔야 기억을 할 수 있을 것 같네요일단 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..