샌드위치 복합소재 대형 안테나구조물의 유한요소 모델링 및 구조 특성 검증

FEA Modeling and Verification of Structural Property of Large Sandwich Composite Antenna Structure

Article information

J. KIMS Technol. 2025;28(1):20-27
Publication date (electronic) : 2025 February 05
doi : https://doi.org/10.9766/KIMST.2025.28.1.020
1)Department of mechanical engineering, Hanwha system, Korea
2)Department of mechanical engineering, UNIST, Korea
3)Department of space & defense engineering, Hankuk Spacepro, Korea
장성욱1),, 소재욱1), 이명규1), 배종인1), 지우석2), 김혜규2), 김동영3)
1)한화시스템㈜ 기계시스템1팀
2)울산과학기술원 기계공학과
3)한국화이바 우주방산사업팀
*Corresponding author, E-mail: sungwook84.jang@hanwha.com
Received 2024 September 4; Revised 2024 November 21; Accepted 2024 November 24.

Trans Abstract

The composite sandwich structure applied to the subject of this study is an effectively designed structure that enhances the specific stiffness of the panel. It features CFRP face sheets, which are strong under tensile and compressive loads, on both sides of an aluminum honeycomb core, a porous structure that is weak in tension and compression but strong in bending. The subject of this study weighs over 4 tons and possesses the largest weight and size among composite antenna structures in the country. Due to the lack of many references for large composite structures, accurate material properties are obtained through specimen tests, and the finite element analysis model is experimentally validated to ensure it accurately represents the actual structure. Finally, the structural robustness is verified through the structural analysis results of the validated model.

1. 서 론

우크라이나 전쟁 이후 미국의 리더십은 약해지고 자국 안보의 필요성이 강하게 대두되었다. 한국의 안보 환경을 고려하여 첨단 기술이 적용된 대공 무기의 수요가 늘고 있다. 다기능, 고성능이 레이다 시스템에 요구되며 레이다의 크기와 무게의 증가는 피할 수 없는 변화이다. 이를 최소화하기 위하여 비강성이 높은 복합소재의 적용을 검토하게 되었다. 복합소재의 높은 비강성으로 인하여 특히 경량화를 요구하는 항공 분야 및 스포츠, 건축 등 다양한 분야에서 사용되고 있다. 이러한 장점에도 복합소재는 고도의 신뢰도를 요구하는 산업 분야에서 2차 구조물로만 사용되었다. 최근에는 지속적인 기술의 발전으로 2차 구조물뿐만 아니라 주 구조물로 사용이 확대되고 있다[1]. 연구 대상에 적용된 복합소재 샌드위치 구조물은 인장, 압축에 약하지만 가벼운 다공성 구조체인 AL honeycomb core 의 양면에 인장, 압축 하중에 강한 CFRP 소재를 사용하여 패널의 비강도 및 비강성을 증가시킨 효과적으로 설계된 구조이다[2].

복합재료는 탄소섬유 구성 방법, 수지의 종류, 섬유의 방향, 제조 방법 등의 설계 요소에 영향을 받아 물성이 변화하기 때문에 일반적 등방성 재료의 설계 기준을 적용할 수 없다[3]. 게다가 복합소재는 작업자 간의 편차, 작업환경, 제품의 로트에 따라 물성의 편차가 존재한다. 물성의 편차가 존재하기 때문에 복합소재는 시험적인 방법을 통하여 대표 물성치를 확인해야 한다. 시험으로 확인된 물성치를 적용한 유한요소 해석 모델을 통하여 해석적인 방법으로 구조적 안전성 평가가 이루어져야 한다[4].

대표 물성치를 확인하기 위하여 시편 시험으로 획득한 데이터를 통계적 계산 방법을 적용하여 보수적인 설계 허용치를 산출한다[5].

H. T. Kim은 연구 대상과 유사한 목적의 복합소재 구조물을 대상으로 외부 진동 환경에 대한 구조물의 강건성을 유한요소 해석 모델을 통하여 검증하는 연구를 수행 하였다[6].

본 연구의 연구 대상은 4톤 이상으로, 복합소재 안테나 구조물 중 국내 최대의 중량과 크기를 가지고 있다. 복합소재 대형구조물 중 참고할 예시가 많지 않아 시편 시험을 통하여 정확한 물성치를 획득하는 과정을 서술하고, 시험을 통해 획득한 물성치를 반영한 유한요소 해석 모델이 실제 제품의 구조적 특징을 반영하는지 확인한다. 확인된 모델의 구조해석을 통하여 구조 강건성을 검증한다.

2. 본 론

2.1 연구 목적

본 연구는 경량화를 위해 설계 시 적용된 샌드위치 복합재 안테나 구조물의 유한요소 해석 모델을 제시하고 시험적 방법을 통하여 해석 결과와 비교 검토하여 해석 모델의 신뢰도를 확인하는 것이 목적이다.

Fig. 1.

Form of the sandwich composite housing of antenna

2.2 연구 대상의 구조

연구 대상은 4톤 이상의 고중량 안테나의 외부 하우징 구조물이다. 고중량의 전자 부품이 수백 개가 설치되어 높은 강성을 요구하고 수송 시 무게의 제약 조건이 있어서 일반적으로 안테나 외부 하우징에 많이 사용되는 소재인 금속 소재는 적용할 수 없었다. 해당 조건을 만족하기 위하여 샌드위치 복합재를 적용하였다. 해당 하우징 구조물은 주요 하중 경로가 되는 외벽과 각종 구성품의 설치부와 내부 구조 분할을 위한 외부 격벽으로 구성된다.

샌드위치 구조물은 굽힘과 비틀림 강성을 증가시키기 위하여 Al 허니콤 코어가 중심에 삽입되어 있고 Carbon/Epoxy 스킨(HPW193/RS1222)이 양면을 덮는 형상으로 구성된다. 샌드위치 구조물의 단면은 Fig. 2와 같다.

Fig. 2.

The formation of sandwich structure

초기 설계 시 복합소재 레이다 하우징은 일체 성형을 하려 하였으나 복잡한 내부 형상으로 인하여 불가능하였다. 해당 하우징은 외부 구조물과 내부 구조물이 별도로 적층, 성형되며 성형된 내·외부 구조물을 고온, 고압 본딩을 통하여 접착하여 완성한다.

2.3 연구 대상의 물성치 획득시험 규격과 결과

레이다 하우징을 구성하는 라미네이트의 물성을 획득하기 위하여 인장, 압축, 굽힘, 전단 시험을 수행한다.

HPW193/RS1222는 이방성 소재이므로 인장시험은 섬유의 직조 방향에 따라 0°, 90°의 두 가지 시편을 제작하여 진행한다. 전단 시험의 경우 섬유의 적층 방향에 의한 영향성이 적다고 판단하여 0°에 해당하는 시편만을 제작하여 수행하였다.

2.3.1 시편 규격

먼저 각 시험 규격에 해당하는 시편을 제작하기 위해서 Lamina를 적층하고, 오토클레이브로 압축 성형한다. 한 장의 Prepreg의 두께는 0.2 mm 정도이다. 압축 성형된 시편은 Table 1의 규격에서 제시하는 권장 크기에 맞게 절단한다. 예를 들어 인장시험의 경우 길이 25.00 mm, 두께 2.8 mm로 제작되었고 제작공차는 0.1 mm 이내로 제작하였다. 해당 규격에 따라 시편은 각각 5개씩 제작하였다. 인장, 압축, 굽힘, 전단 시험의 시편의 규격은 하기 표를 참고한다.

Standards for specimen property testing

2.3.2 복합소재 물성치 획득

구조물 강도의 적절성을 평가하기 위하여 시편의 파단 강도를 확인한다. Table 1의 규격에 따라 시험을 수행하였다. 해당 소재의 시편의 제작 및 시험의 세팅 과정에서 이상적인 조건을 충족하기 어렵기 때문에 금속 소재에 비하여 결과의 편차가 크다. 그러므로 인장, 압축, 굽힘, 전단 시험 후 통계적 방법을 사용하여 대푯값을 획득한다.

금속 소재는 일정 수 이상의 시편 시험을 진행하고 해당 데이터 집단의 제 1 백분위수 혹은 제 5 백분위수를 소재의 대표 물성치로 사용한다. 제 5 백분위 수를 구하기 위해서는 최소 30개 이상의 시편이, 제 1 백분위 수를 구하기 위해서는 최소 90개 이상의 시편이 필요하다. 복합소재는 금속보다 시험을 구성하는 규격이 복잡하고 시편을 제작하는 비용 또한 고가이기 때문에 시험마다 30개 이상의 시편을 구성하는 데에 어려움이 있을 수 있다. MIL-HDBK-17F Chapter 8에서는 이를 염두하여 적은 수의 시편 데이터로도 제 1 백분위수 혹은 제 5 백분위수를 구하는 통계적 처리 방법을 제시하고 있다. 다음은 통계적 처리 과정을 통해 복합소재의 대표 물성치를 산출하는 과정이다.

Fig. 3.

Tensile strength & interlaminar shear strength test & flexure test

2.3.2.1 통계적 분포 모델 선택

측정한 데이터와 일치하는 통계적 분포 모델로 정규분포도를 선택한다. 정규분포도의 주요 인자는 다음과 같다.

  • 1. 모집단의 평균

    (1) x¯=1ni=1nxi

  • 2. 표준 편차

    (2) s=1n1i=1n(xix¯)2

  • 3. 정규 분포 함수

    (3) F(x)=12πs2e(xx¯)22s2

그려진 정규분포도에서 제 1 백분위 수를 대표 물성치로 결정한다. 더 보수적인 설계, 현실적인 시험의 제약사항으로 인하여 제 1 백분위수를 파단 강도로 결정하였다. 그 결과는 하기 Table 3과 같다. 복합소재의 경우 취성이 높아 별도의 항복강도가 존재하지 않으므로 획득한 파단 강도를 기준으로 구조물의 강도 설계가 적절한지 평가한다.

Breaking strength table of CFRP

파단 강도와 달리 강성의 경우 ASTM D 3039을 통하여 young's modulus를 확인한다. young's modulus의 경우 실물의 현상을 반영해야 하므로 보수적인 결과를 획득하는 통계적 방법 대신 시험 결과의 평균값을 적용한다.

시험을 통하여 0도와 90도 적층 방향의 물성치 만을 획득했지만, 강성행렬의 회전변환을 통해 45도 방향 등 임의의 각도에서의 물성치를 획득할 수 있다. 0도의 young's modulus 평균값을 대입하여 45도 방향의 물성치를 확보하였다. 면내 방향으로 각도만큼 회전하였을 때 강성행렬은 회전변환행렬 T를 이용하여 다음과 같이 계산된다. 45도 방향의 young's modulus는 Table 2를 참고한다.

The directional young's modulus of CFRP

(4) C¯=TCTTT=[cos2θsin2θ0002sinθcosθsin2θcos2θ0002sinθcosθ001000000100000010sinθcosθsinθcosθ000cos2θsin2θ]

2.4 유한요소 해석 모델링 구성

시험을 통하여 획득한 young's modulus를 유한요소해석 모델에 반영한다. 복합재 스킨은 외벽의 경우 [0/45/0/45/0/45/0]의 7층 적층 내벽의 경우 [0/45/0/45/0]의 5층 적층을 사용하였다. 초기 해석 모델은 전체 Shell mesh로 구성하였다. 각각의 파트는 노드를 공유하여 연결하였다.

2.5 모드 시험과 modal 해석 결과

2.5.1 임팩트 해머를 통한 모드 시험

시험 대상의 구조적 특성을 실험적으로 확인하기 위하여 임팩트 해머 충격시험을 수행한다. 임팩트 해머 충격시험은 시험 대상에 모드 형상을 표현할 수 있도록 모든 절점에 가속도 센서를 설치(66개)하고 해머로 대상에 충격을 가하여 센서의 응답 데이터를 수집한다. 시험의 신뢰도를 높이기 위하여 동일한 시험을 10회 반복하고, 유의미한 값을 추출하여 고유진동수와 진동의 형상을 확인한다. 시험 대상은 free-free 조건의 유한요소 해석 모델에 유사하도록 매달아두었다(Fig. 4 참고).

Fig. 4.

Location of Check point for accelerometers

2.5.2 유한요소 해석/모드 시험 결과 비교

임팩트 해머를 통한 모드 시험과 유한요소 해석의 결과를 비교한다. 5차 모드까지의 고유진동수의 오차율을 비교하고 추출된 mode shape의 일치성을 확인한다. 고유진동수는 유한요소 해석과 시험의 결과의 오차율이 최대 11.45 %임을 확인하였다(Table 4 참조).

Comparison of initial modal analysis results with experimental results

2.5.3 FEM 개선

유한요소 해석을 통한 시험 결과와 측정된 고유진동수의 오차가 크다고 판단하여 FEM 모델링을 개선하였다. 초기 해석 모델은 shell mesh로 구성되며, 내·외부를 나누어 모델링을 구성하였다. 내·외부 연결부는 노드를 공유하는 형태로 하나의 파트처럼 인식되는 형태이다.

이를 개선하기 위해서 하우징의 외벽은 Al 코어에 내부 스킨과 외부 스킨을 나누어 모델링하고 Tie constraints로 결합한다. 모델링 된 외벽에 수직하게 연결되는 구조도 Tie constraints으로 고정한다(Fig. 5 참고).

Fig. 5.

Improvement of finite element modeling, a) initial model, b) Improved model

Tie constraint는 shell 요소와 solid 요소 간 접착하는 노드 간 병진 및 회전 자유도를 결합한다. 이 과정에서 solid 요소와 결합되는 shell 요소의 가장자리 노드는 solid 요소 표면에 서로 가까운 곳에 붙어 두께의 중심이 된다. 표면을 따라 두께를 형상하므로 실제 형상에 가깝게 요소가 구성되며 복합소재와 같이 shell 요소로 모델링 되어야 하는 부품과 형상으로 인하여 solid 요소로 모델링 되어야 하는 부품의 결합에 적합하다.

하우징 외벽은 국소적으로 두께와 곡률이 변화하는 복잡한 형상을 가졌기 때문에 중심 면을 정의하기 어렵다. 따라서 외곽 면과 내곽 면을 분리하여 각각 쉘 요소로 구성된 곡면으로 정의한 뒤 두 곡면을 결합하는 방법을 사용하였다(Fig. 5, b 참고). 시험 대상은 하우징의 외부 파트와 내부 격벽을 별도로 제작하여 본딩으로 결합하였다. 본딩 면은 별도로 모델링 하여 정의하였다. 접착 면에 가해지는 응력을 분석하기 위한 목적도 있다. 복합소재로 내부 파티션과 본드는 서로 수직 한 면이므로 Tie constraints를 적용해 결합하였다(Fig. 6 참고).

Fig. 6.

Improvement of finite element modeling(case 2, Bond modeling)

상기 외에 형상으로 인하여 쉘 요소로 구현할 수 없는 부분이 존재한다. 해당 부위는 3차원 솔리드 요소로 모델링 되었다. 하우징 정면의 블록은 외벽의 각 위치에 Tie constraints로 결합되었다. 해당 솔리드 요소의 형상은 Fig. 7을 참고한다.

Fig. 7.

Section of FEM, a) shell element, b) solid element

기존 모델에서 상기 내용 반영 시 모델링의 무게는 약 920 kg에서 약 950 kg으로 30 kg 정도의 증량이 있었다. 실물 중량은 948 kg으로 변경된 중량과의 중량 차이는 0.2 %이다.

Material properties table of antenna housing

모델링 개선 과정에서 강성이 낮은 본딩 부위를 모델링에 반영하고 개선 전·후 증량의 차이가 있었으나 1차 모드가 발생하는 후방 부에 강성이 높은 솔리드 요소를 적용하며 오히려 1차 고유진동수는 증가하였다. 이로 인하여 임팩트 해머 모드 시험 결과와의 일치성이 상승하였다. 오차율은 최대 4.11 % 임을 확인할 수 있다. 이는 모델링의 구속 조건은 6방향 구속 조건이 풀려있는 조건이나, 실제 시험 대상은 매달려 있음을 고려하였을 때 수용할 만한 오차라고 판단하였다.

5차까지의 모드 형상은 시험과 해석 결과가 일치한다. 1차 모드는 강성이 약한 후면에서, 2차 모드는 구조물 전체의 뒤틀림, 3차와 4차 모드는 후면에서, 5차 모드는 상단에서 발생한다. 모드 형상의 일치함을 통하여 유한요소 해석 모델이 충분히 안테나 하우징 구조물의 구조적 특성과 유사하게 구현하였다고 판단할 수 있다. 이를 통하여 유한요소 해석의 신뢰성을 확보하였다.

Comparison of modal analysis results with experimental results

2.6 정하중 해석을 통한 구조안전성 확인

앞서 신뢰성이 확인된 유한요소 모델을 통하여 정하중 해석을 수행한다. 경계조건으로 외부 힌지 설치부의 모든 자유도를 구속하였다. 운용 시 발생할 수 있는 하중을 고려하여 -X 방향으로 1 g, -Z 방향으로 3 g의 가속도가 적용되었다. 이때 발생하는 최대응력을 복합소재 스킨의 적층 방향에 따라 1, 2, 3방향으로 나누어 파단 강도와 비교하였다. 그리고 방향에 무관하게 내외부 구조 간 접합을 위하여 적용된 본드의 최대응력을 데이터 시트에서 제공하는 강도와 비교하였다.

이때 최대 응력은 내부 알루미늄 구조물이 장착되는 내벽과 외부 힌지 설치부에서 발생하였다(Fig. 8). 정하중 해석을 통하여 확인한 최대응력에 대한 안전율은 파단 강도 대비 3.0 이상으로 충분히 강건한 구조 설계가 되었다는 것을 확인하였다. 하기 Table 7을 통하여 방향별 최대 응력을 확인한다.

Fig. 8.

0° direction stress distribution of composite structure under static load

Analysis results and safety factor of composite structure under static load (X: 1 g, Z: 3 g)

Fig. 9.

90° direction stress distribution of composite structure under static load

Fig. 10.

Shear stress distribution of composite structure under static load

3. 결 론

이 연구에서 재료 물성치 시험을 통하여 복합소재의 물성치와 파단 강도를 확인하고 이를 반영한 유한요소모델링을 통한 modal 해석 결과와 임팩트해머 시험 결과를 비교하였다. 비교를 통하여 신뢰도가 확보된 유한요소 해석 모델링에 운용 시 발생할 하중을 적용하여 안테나 유니트 하우징의 구조 강건성을 확인하였다.

4톤 이상, 약 4 m × 3 m × 1 m의 대형, 복합소재 구조물의 경우 크기, 무게, 인가주파수의 영역으로 인하여 진동, 충격시험에 대한 환경시험을 수행할 수 있는 장비가 드물다. 그래서 환경시험을 수행하는 대신 유한요소 해석을 통하여 구조적 강건함을 검증하고자 유한요소 해석 모델링의 신뢰도를 확인하였다.

modal 해석을 통하여 확인한 5차까지의 고유진동수와 시험 결과의 최대 오차는 4.11 %이고 모드 형상은 일치한다. 충분히 신뢰할 만한 모델링을 구성하였다고 판단하고 이 모델링을 통하여 연구 대상의 구조안전성을 검증하였다. 발생한 오차는 구조안전성을 확인하는 과정에서 안전율을 높여 보수적으로 판단하는 요소로 반영할 예정이다.

References

[1]. . Olsson R.. Engineering Method for Prediction of Imoact Response and Damage in Sandwich Panels. Jounal of Sandwich Structures&Materials 4(1):3–29. 2002;
[2]. . Vinson J. R.. The behavior of sandwich structures of isotropic and composite materials. Routledge, Inc. England 1999.
[3]. . Park S. B., Park J. H., Ahn K. C.. Strength Prediction of Bolted Woven Composite Joint Using Characteristic Length. Jounal of the Korean Society of Safety 18(3):8–15. 2003;
[4]. . Ko H. Y., Shin K. B., Jeong J. C.. A Study on the Comparison of Structural Performance Test and Analysis for Design Verification of Bimodal Tram Vehicle made of Sandwich Composite. Journal of the Korean Society for Railway 12(4):518–525. 2009;
[5]. . “Department of Defense Handbook: Composite Materials Handbook - Polymer matrix composites guidelines for characterization, MIL-HDBK-17-1F Volume 1 of 5 Departement of Defense. US: 2002.
[6]. . Kim H. T.. Characteristic Evaluation Method of Massive Sandwich Composite-strucure under Transfer Vibration Condition. Transactions of the Korean Society for Noise and Vibration Engineering 27(3):297–305. 2017;

Article information Continued

Fig. 1.

Form of the sandwich composite housing of antenna

Fig. 2.

The formation of sandwich structure

Table 1.

Standards for specimen property testing

Testing Method No. of specimens
Tensile Test ASTM D 3039 5(2종)
Compression Test ASTM D 695 5
Flexure Test ASTM D 790 5
Interlaminar Shear Test ASTM D 2344 5

Fig. 3.

Tensile strength & interlaminar shear strength test & flexure test

Table 2.

The directional young's modulus of CFRP

Young's modulus (GPa)
Young's modulus test of CFRP Average
45.5 45.8 45.2 45.7 47.5 46.0
90° 43.0 43.2 44.4 43.9 43.3 43.6
45° 24.2

Table 3.

Breaking strength table of CFRP

Test Direction Parameter Breaking Strength (MPa)
Tensile Test Average 887.2
Design Allowable value 776.5
90° Average 790.4
Design Allowable value 749.2
Compression Test Average 520
Design Allowable value 315.6
Flexure test Average 757
Design Allowable value 622.1
Interlaminar Shear Test Average 66.5
Design Allowable value 53.0

Fig. 4.

Location of Check point for accelerometers

Table 4.

Comparison of initial modal analysis results with experimental results

Mode Natural frequency(Hz) Error(%)
Test Analysis
1 15.2 14.229 7.83
2 19.3 18.762 2.86
3 27.4 26.472 3.51
4 33.4 32.067 4.16
5 35.7 40.316 11.45

Fig. 5.

Improvement of finite element modeling, a) initial model, b) Improved model

Fig. 6.

Improvement of finite element modeling(case 2, Bond modeling)

Fig. 7.

Section of FEM, a) shell element, b) solid element

Table 5.

Material properties table of antenna housing

Bond Al6061 Resin
Density ρ (g/cm3) 1.5 2.7 1.6
Young's Modulus E (GPa) 4.237 68.9 46.0
Poisson Ratio ν12 0.4 0.27 0.31
Shear Modulus G12(GPa) 1.513 27.56 3.73
Breaking Strength (MPa) 32.0 110 550

Table 6.

Comparison of modal analysis results with experimental results

Mode Natural frequency (Hz) Error(%)(after-test)
Test Analysis (before) Analysis (after)
1 15.2 14.2 15.7 3.5
2 19.3 18.8 19.2 0.5
3 27.4 26.5 27.6 0.6
4 33.4 32.1 34.8 4.1
5 35.7 40.3 36.2 1.5

Fig. 8.

0° direction stress distribution of composite structure under static load

Fig. 9.

90° direction stress distribution of composite structure under static load

Fig. 10.

Shear stress distribution of composite structure under static load

Table 7.

Analysis results and safety factor of composite structure under static load (X: 1 g, Z: 3 g)

90° Shear direction Bond
Breaking Strength (MPa) 776.5 749.2 66.5 32.0
Max. Stress (MPa) 103.7 91.91 13.75 10.49
Factor of Safety (MPa) 7.49 8.15 4.84 3.05