[논문 리뷰] Free Vibration analysis of Curvilinearly Stiffened Composite plates with an arbitrarily shaped cutout using Isogeometric Analysis

Balakrishnan Devarajan|arXiv (Cornell University)|2021. 04. 26.

Advanced Numerical Analysis Techniques인용 수 10

한 줄 요약

이 논문은 NURBS 기반 단일 패치 모델링을 사용한 고등 이소지오메트릭 해석(IsoGeometric Analysis, IGA) 프레임워크를 제시하여 곡선형 보강재가 부착된 복합재판의 자유 진동 해석을 수행한다. 이 방법은 판과 보강재를 별도로 모델링하며, 새로운 매핑 기법을 통해 변위 호환성을 확보하고, 정확도와 계산 효율성이 뛰어나 문헌 자료와의 검증을 통해 입증되었으며, 플라이 방향, 컷아웃 기하학, 경계 조건에 대한 파arametric 연구에 응용되었다.

ABSTRACT

This paper focuses on the isogeometric vibration analysis of curvilinearly stiffened composite panels. The stiffness matrices and the mass matrices are derived using the first-order shear deformation theory (FSDT). The present method models the plate and the stiffener separately, which allows the stiffener element nodes to not coincide with the plate shell-element nodes. The stiffness and mass matrices of a stiffener are transformed to those of the plate through the displacement compatibility conditions at the plate/stiffener interface by interpolation using NURBS basis functions. Cutouts are modeled using a single NURBS patch generated by creating a ruled surface between two curves. The proposed formulation is first validated by comparing it with available literature. The effects of width-to-thickness ratio, fiber orientation, ply layups, shape and size of the cutouts and the boundary conditions on the response of stiffened composite plates are then analyzed and the numerical results are used to derive useful conclusions.

연구 동기 및 목표

• 복잡한 컷아웃을 가진 곡선형 보강 복합재판의 자유 진동 해석을 위한 계산 효율성과 정확도가 높은 이소지오메트릭 해석 방법을 개발한다.
• 기존의 유한요소법(FEM)의 한계를 극복하기 위해 NURBS 기반 이소지오메트릭 해석을 통해 정확한 기하 표현과 더 높은 연속성을 확보한다.
• 판과 보강재를 별도로 모델링하면서도 새로운 매핑 기법을 통해 인터페이스에서 변위 호환성을 유지한다.
• 컷아웃 형상, 크기, 위치, 섬유 방향, 플라이 레이업이 보강 복합재 패널의 동적 거동에 미치는 영향을 조사한다.
• 기존 문헌과의 비교를 통해 제안된 IGA 프레임워크의 타당성을 검증하고, 수치 기준시험을 통해 복잡한 기하 구조에 대한 적용 가능성을 입증한다.

제안 방법

• 전체 판, 보강재, 컷아웃을 단일 패치로 표현하기 위해 비균일 유리 B-스플라인(NURBS)을 사용하여 정확한 기하 표현을 확보한다.
• 판 및 보강재 요소의 강성 및 질량 행렬을 도출하기 위해 1차 전단변형 이론(FSDT)을 적용한다.
• 판과 보강재 노드의 메esh가 일치하지 않을 경우에도 변위 호환성을 확보하기 위해 새로운 매핑 기법을 구현한다.
• 두 개의 경계 곡선 사이에 생성된 룰드 서피스를 사용하여 컷아웃을 모델링함으로써 단일 NURBS 패치 내에서 임의의 컷아웃 형상을 구현한다.
• 기하 및 장 근사에 모두 NURBS 기반 함수를 사용하는 등가요소 형식을 적용하여 높은 연속성을 확보하고 메쉬 생성을 회피한다.
• 일반화된 고유값 문제를 해결하여 고유 진동수와 진동형상을 추출하며, 결과는 정규화되어 기준 데이터와 비교된다.

실험 결과

연구 질문

• RQ1 단일 NURBS 패치를 사용한 제안된 IGA 프레임워크는 복잡한 컷아웃을 가진 곡선형 보강 복합재판의 고유 진동수를 예측하는 데 얼마나 효과적인가?
• RQ2 컷아웃 형상, 크기, 위치는 보강 복합재 패널의 자유 진동 거동에 어떤 영향을 미치는가?
• RQ3 다양한 플라이 레이업(대칭, 반대칭, 각도 플라이, 크로스 플라이)과 섬유 방향은 동적 거동에 어떤 영향을 미치는가?
• RQ4 경계 조건과 폭 두께 비율(γ)은 고유 진동수와 진동형상에 어떤 영향을 미치는가?
• RQ5 제안된 매핑 기법이 별도로 모델링된 보강재와 판 사이의 정확한 변위 호환성을 얼마나 잘 보장하는가?

주요 결과

• 제안된 IGA 방법은 고유 진동수 예측에서 높은 정확도를 달성하였으며, 문헌 자료 및 기타 고급 수치 해석 방법과의 결과가 매우 유사하게 일치하였다.
• γ = 5 및 Δε = 0일 때, 반대칭 크로스 플라이 레이업의 첫 번째 모드 진동수는 23.986(정규화된 값)이며, 대칭 각도 플라이 레이업의 경우 23.572로, 플라이 방향에 민감함을 보였다.
• 폭 두께 비율(γ)을 5에서 10으로 증가시키면 고유 진동수가 크게 상승하였으며, 반대칭 크로스 플라이 형상에서 첫 번째 모드 진동수는 약 24에서 약 25로 증가하였다.
• 컷아웃 곡률(Δε = 0.25)을 포함시킬 경우 Δε = 0일 때 대비 첫 번째 모드 진동수가 최대 약 10% 증가하여, 컷아웃 형상이 동적 스태틱 강성에 상당한 영향을 미친다는 것을 시사하였다.
• 모드 형상도(예: 다섯 번째 모드)는 플라이 레이업 및 기하 구조에 따라 다릅니다. 반대칭 및 각도 플라이 형상에서는 대칭 형상보다 더 복잡한 절점 패턴을 보였다.
• 이 방법은 복잡한 컷아웃 형상과 다양한 보강재 곡률을 포함한 다양한 구성 조건에서도 뛰어난 강건성과 효율성을 보였으며, 이를 통해 구조 최적화 워크플로우에서의 활용 가능성을 뒷받침한다.