공대생 도전 일지

// 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 =..