HOB - MITC4 Shell 계산 과정

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.577

stain_Displacement - 가우시안, thick, node, norm vec 들고간다. -> 쉘 요소의 공변 변형률 계산 
Assumed_Convariant_Strain - 가정된 공 변형률 계산 ? 
Covariant_Strain - 변형률 연산...? (r,s,t, rs, st, tr)
Derivative_Displacement - r,s,t에 대한 변위를 각각 설정해준다.
Shape_Function - H를 구해준다. == 각 노드에 대응하는 모양 함수 값 
dHrs_Matrix - H에서 r,s로 각각 미분한 dHrs을 보내준다. == 기하학적 특성 계산 

Jacobian - H, dHrs를 가지고 r,s,t 국소 좌표계를 실제 물리적 좌표계로의 변환 

이 것들 가지고 covariant_strain 구하기

 

Contravariant_Basis - 반 공변 기저 백터(contravariant basis vectors)를 계산 

Local_Basis - 기저 벡터 반환 

 

어 음

이게 뭐지

일단 MITC4로 시작 

Shell_stiffness
Material_law
for i j k 2
Gauss22
Strain_Displacement
	Assumed_Convariant_Strain
		Covariant_Strain*4
			Derivative_Displacement	H, dHrs를 통해 r,s,t 변위 각각 설정
			Jacobian	r,s,t좌표계를 실 좌표계로
			Covariant_Basis	r,s,t 방향 각각 뽑기
			변형률 연산하여 추출
		같은 Cov인데 r,s,t
	Jacobian
	Contravariant_Basis	반 공변 기저 벡터
	Local_Basis	기저 벡터
	공변 변형률 성분 계산
Jacovian
K-matrix 연산