In the present study, the fiber and matrix failure of composite laminates under arbitrary biaxial stresses were evaluated based on separate mode criteria such as Hasnin and Puck theories. There is a limitation to predict the fiber-dominant and/or matrix-dominant failures under arbitrary stress states using limit criteria (maximum stress and maximum strain theories) and interactive criteria (Tsai-Hill and Tsai-Wu theories). There is little literature for failure analysis of ships and offshore composite structures considering advanced failure theories such as Hashin and Puck theories. Furthermore, there is not enough practical commercial finite element analysis (FEA) code which is basically adopted the separate mode criteria. Hence, in the present study, the user-defined subroutine of commercial FEA code ABAQUS for evaluation of fiber and matrix failures of composite structures was developed based on Hashin and Puck failure criteria. And then, the proposed subroutine was validated by comparing with a series of experimental results of carbon- and glass-implemented composite laminates to guarantee the reliability and usefulness of the developed method.
복합재료의 여러 장점, 이를테면, 경량성을 비롯한 우수한 비강성·비강도, 유연성, 내부식성, 진동·소음 저감성능, 전자파 흡수성능 등으로 인해 최근 조선해양분야에서의 복합재료 활용 사례가 꾸준히 증가하고 있다. 특히, 기존의 복합재료가 주로 경량성에 초점이 맞춰졌다면, 최근에는 경량성 뿐만 아니라 상기에 전술한 기타 우수한 물성에 그 관심이 집중되고 있다. 예를 들면, 초심해용 URF(Umbilical·Riser·Flowline)는 릴링(reeling) 성능이 그 무엇보다 중요한데, 금속재료보다 유연하지만 (낮은 탄성계수) 강한 (높은 인장강도) 복합재료의 적용을 통해 금속재료 URF에 비해 우수한 Reeling 성능을 가진 복합재료 URF의 제작이 가능하다 (Pavlou, 2013). 또한 복합재료 프로펠러는 금속재료 프로펠러에 비해 가벼울 뿐만 아니라 내부식성이 우수하며 수중방사소음 저감능력이 탁월하다 (Lee, et al., 2012).
조선해양분야에 채택되는 복합재료 구조물의 안전한 운용을 위해서는 해석단계에서 정밀한 구조파손평가를 수행하여야 한다. 즉, 복합재료 구조물에서 나타날 수 있는 파손유형인 섬유파손(fiber breakage), 횡기지파손(transverse matrix cracking), 전단기지파손(shear matrix cracking) 등을 정밀하게 평가하여야 한다. 2차원 문제에서의 각각의 파손유형을 살펴보면, 섬유파손은 종방향(longitudinal) 응력
현재까지 매우 다양한 복합재료 파손이론이 많은 과학자들에 의해 제시되었으며, 이들을 크게 세 가지로 분류하면 다음과 같다 (Nahas, 1986; Hinton, et al., 2004; Lee, et al., 2007; Barbero, 2011; 2013). 첫째로 한계 기준(limit criteria)은 각각의 적층(lamina) 응력성분
둘째로 상호작용 기준 (interactive criteria)은 모든 응력 (혹은 변형률) 성분을 포함한 단일 이차방정식(single quadratic equation)이나 텐서 다항방정식(tensor polynomial equation)을 활용하여 등식이 만족될 때 파손이 일어나는 것으로 판단하는 기준이다. 특히, 다항식을 구성하는 종방향 응력 성분항, 횡방향 응력 성분항 및 전단응력 성분항을 각각 도출하여 어떠한 성분이 파손에 크게 기여했는지를 분석함으로써 섬유파손 혹은 기지파손을 판단할 수 있다. 대표적인 이론으로는 Tsai-Hill, Tsai-Wu 등이 있다.
마지막으로 분리모드 기준(separate mode criteria)은 상기의 두 기준을 보강한 방법으로서, 섬유파손과 횡기지파손, 전단기지파손을 별개의 다항방정식으로 분리하여 평가할 수 있을 뿐만 아니라, 다항식 선택 기준식이 명확하다는 장점이 존재한다. 특히, Tsai-Hill 및 Tsai-Wu 이론 등과 같은 상호작용 기준이 가지는 문제점, 예를 들면,
수많은 파손이론 가운데 현재까지 가장 널리 활용되고 있는 기준은 최대응력이론, 최대변형률이론, Tsai-Hill이론 및 Tsai-Wu이론인데 (Sun, et al., 1996) 그 이유는 이들 기준이 비교적 단순하여 적용이 매우 간편하고 식에서 요구되는 재료물성치도 그리 많지 않기 때문이다.
그러나 전술한 바와 같은 문제점도 또한 존재하는데, 이러한 문제를 극복하기 위해 적층 복합재료의 섬유 및 기지 파손을 각각 분리하여 고려하는 새로운 시도가 1980년 Hashin에 의해 제안된 바 있으며 (Hashin, 1980) 1998년 Puck에 의해 더욱 구체적으로 체계화되었다 (Puck & Schürmann, 1998; 2002). 나아가 이들 파손 기준은 여러 복합재료 역학자들에 의해 실험적 및 수치적으로 검증된 바 있다 (Hinton & Soden, 1998; Soden, et al., 1998a; 1998b; Hinton, et al., 2002; Kaddour, et al., 2004).
이에 최근 항공우주분야와 기계분야에서 분리모드 기준 기반 파손 평가 사례가 꾸준히 늘고 있으며, 특히, 비교적 파손 방정식이 단순한 Hashin이론이 많이 채택되고 있다 (Kim & Kwon, 2013; Park, et al., 2014). 그러나 조선해양분야에서는 아직까지 Hashin이론의 적용 사례가 많지 않을 뿐만 아니라, Hashin이론보다 더욱 정밀한 Puck이론의 적용 사례는 극히 일부만 있을 뿐 (Son, et al., 2012), 그 예를 찾아보기가쉽지 않다. 또한, 이들 파손기준식이 기본적으로 탑재된 범용 유한요소해석코드는 극히 일부에 지나지 않아 복합재료 구조물의 파손평가가 쉽지 않은 실정이다.
그럼에도 불구하고 풍력발전기용 블레이드를 포함한 복합 재료 기반 산업구조물의 정적해석 시 상기 파손기준식의 채용을 의무화하고 있으며 (Germanischer Lloyd, 2010), 해외 선진국의 관련 연구 사례에서도 해당 파손기준식의 적용 예가 꾸준히 늘어나고 있는 추세이다
따라서 본 연구에서는 Hashin 및 Puck 파손이론에 기반한 적층 복합재료의 섬유 및 기지파손 평가를 위한 수치해석 방법을 제안한다. 즉, 섬유파손과 기지파손을 모두 평가할 수 있는 Hashin 및 Puck 파손이론을 바탕으로 범용유한요소해석 코드용 사용자정의 서브루틴을 구축하고 시리즈 파손해석을 수행한다. 나아가 문헌으로 심층 보고된 적층 복합재료의 섬유·기지파손의 결과와의 비교를 통해 제안 수치해석 방법의 적합성 및 유용성을 검증한다.
1980년 Hashin은 기존의 상호작용 기준식이 파손 판정 계산 시 복합재료를 단일적층구조로 가정하고 있는 이유로 섬유 파손이나 기지파손을 정밀하게 예측할 수 없음을 보인 바 있다 (Hashin, 1980).
Hashin 파손기준식을 나타내면 다음과 같다.
(1) Tensile fiber failure for
(2) Compressive fiber failure for
(3) Tensile matrix failure for
(4) Compressive matrix failure for
여기서
한편, Puck and Schürmann은 Hashin 파손기준식을 더욱 정밀하게 보완하였는데, 이를테면, 특정 조합하중에서 나타날 수 있는 기지파손을 세 가지로 분류하였으며, 기지파손으로 새롭게 생성되는 파단면에 대한 각도(angle of fracture plane) 결정기법을 제안하였다 (Puck & Schürmann, 1998; 2002).
Puck 파손기준식을 나타내면 다음과 같다.
(1) Fiber failure in tension
(2) Fiber failure in compression
(3) Matrix failure Mode A for transverse tension
(4) Matrix failure Mode B for moderate transverse compression
(5) Matrix failure Mode C for large transverse compression
여기서
Mode A, Mode B 및 Mode C의 기지파손모드를 섬유수직 방향응력-전단응력 평면 상에 나타내면 Fig. 1과 같다.
[Fig. 1] Graphical description for matrix failure Mode A, Mode B and Mode C (Puck & Schurmann, 1998)
한편, 임의의 조합하중에 놓인 적층 복합재료의 Puck 파손 기준식 적용을 위한 응력상태별 판단기준은 Table 1과 같다. 단, Table 1에서 는 수직-수직전단응력으로 인한 파단에 반하는 활성면의 파단 저항 파라미터,
[Table 1] Application condition of Puck failure criteria
Application condition of Puck failure criteria
한편, Puck 파손기준을 결정짓는 주요한 파라미터인 , , , 및
하지만, 식 (10), (11) 및 (13)의 주요 파라미터들은 미분항을 풀어야만 결정할 수 있다는 단점이 존재한다. 이에 따라 Puck 등은 실험적 및 현상학적 고찰을 통해 보다 간편한 재료 상수를 Table 2와 같이 제안하였다 (Puck, et al., 2002).
[Table 2] Recommendation of the inclination parameters
Recommendation of the inclination parameters
전술한 Hashin 및 Puck 파손기준식을 바탕으로 범용유한 요소해석코드 ABAQUS 전용 사용자정의 서브루틴을 구축한다. 본 연구에서는 UVARM 서브루틴을 통해 파손기준식을 계산하게 되는데, 여기서 UVARM이란 ABAQUS SOLVER를 통해 계산되어지는 특정 물리량을 사용자가 임의로 추출 및 가공함으로써 사용자가 원하는 새로운 변수를 도출해 낼 수 있도록 ABAQUS에서 제공하고 있는 서브루틴의 일종이다.
본 연구에서는 ABAQUS를 통해 각 증분마다 계산되어지는 적층 복합재료의 응력 및 변형률 정보를 바탕으로 식 (1)-(14)를 계산할 수 있으며 최종 계산 결과인 파손 분포를 ABAQUS 상에서 확인할 수 있다.
본 연구의 해석 대상은 탄소·유리섬유 및 에폭시수지로 구성된 적층 복합재료이다. 이미 여러 역학자들에 의해 임의의 조합하중에 놓인 적층 복합재료의 파손에 관한 연구가 실험적으로 수행된 바 있으며 본 연구에서는 Soden등에 의해 수행된 문헌을 참고하였다 (Soden, et al., 1998a; 1998b).
본 연구에서 사용된 재료는 E-glass 21xK43 유리섬유 및 LY556 에폭시(epoxy) 수지로 이루어진 적층 복합재료와 AS4 탄소섬유 및 3501-6 에폭시수지로 이루어진 적층 복합재료가 사용되었다. 대상 시험편의 플라이 적층배열, 재질 및 작용하중 등에 관한 상세 정보를 Table 3 (Soden, et al., 1998a; 1998b)에 나타낸다.
[Table 3] Details of laminates and loading conditions
Details of laminates and loading conditions
한편, 각 복합재료의 주요 재료물성치는 Table 4 및 Table 5 (Soden, et al., 1998a; 1998b)와 같다. 여기서 E-glass/LY556 적층 복합재료를 구성하는 각 플라이의 두께는, 90° 플라이의 경우, 0.172mm, +30° 및 −30° 플라이의 경우, 각각 0.414mm이며 따라서 플라이의 총 두께는 2mm이다. 또한 AS4/3501-6 적층 복합재료를 구성하는 각 플라이의 두께는 모든 섬유방향에 대해 0.1375mm이며 따라서 플라이의 총 두께는 1.1mm이다.
[Table 4] Material properties of E-glass/LY556 laminate
Material properties of E-glass/LY556 laminate
[Table 5] Material properties of AS4/3501-6 laminate
Material properties of AS4/3501-6 laminate
한편, 식 (5) 및 (6)에 사용된 변수
본 연구에서는 Table 3과 같은 다축 응력상태에 놓인 복합 재료 구조물의 일부를 추출하여 유한요소모델을 생성하였으며, 구조물의 일반적인 수준에서 고려할 수 있는 크기로 설정하였다. 즉, 유한요소모델의 길이 및 폭의 치수는 각각 300mm 및 300mm이다.
한편, 적층 복합재료에 작용되는 응력은 Soden 등에 의해 수행된 복합재료 시험조건과 동일한 크기·방향의 응력 (Table 3의 하중조건)을 부여하였다 (Soden, et al., 1998a; 1998b). Fig. 2는 Table 3의 Case B에 대한 하중 및 경계조건을 적용한 모습을 나타내고 있다.
Fig. 3은 E-glass/LY556 적층 복합재료의 ABAQUS 해석용 유한요소모델에 삽입된 적층배열과 각 플라이에서의 섬유 방향을 나타내고 있다.
본 절에서는 적층 복합재료의 초기 파손현상을 개발된 사용자 정의 서브루틴을 활용하여 전산해석을 수행함으로써 개발 코드의 적합성과 유용성을 검증한다.
Fig. 4에 보인 바와 같이, Soden, et al. (1998a; 1998b)에 의해 확보된 E-glass/LY556 및 AS4/3501-6 적층 복합재료의 초기 파손응력 포락선 (failure stress envelope)을 바탕으로 (실선으로 표기) 총 24가지의 전산해석 케이스를 임의로 선정하였다 (점으로 표기).
저자에 의해 임의로 선정된 24가지 해석 케이스의 하중 조건을 Table 6에 나타내었다. 이때, 적층 복합재료의 파손 및 비파손 상태의 정확한 평가가 이루어졌는지를 판단하기 위해 12가지의 파손 및 비파손 상태를 각각 선정하였다.
[Table 6] Failure analysis cases of composite laminates
Failure analysis cases of composite laminates
즉, Fig. 4와 Table 6에 보인 바와 같이, Case A-1, A-3, A-5, A-7, B-1, B-3, B-5, B-7, C-1, C-3, C-5 및 C-7은 다축응력 상태에서 초기파손응력 포락선의 내부에 위치하고 있으므로 파손이 발생하지 않는 경우이며, Case A-2, A-4, A-6, A-8, B-2, B-4, B-6, B-8, C-2, C-4, C-6 및 C-8은 외부에 위치하고 있으므로 초기 파손이 적층 복합재료의 임의의 부위에서 발생하는 경우이다.
전산해석 결과의 대표적 예로서 Fig. 5, Fig. 6 및 Fig. 7에 Case A-2, B-2 및 C-2의 조합하중에 놓인 대상 적층 복합재료의 Hashin 및 Puck 파손기준식 기반 전산해석 결과 (파손발생)를 보이고 있다. Fig. 5-Fig. 7에 보인 바와 같이 개발 ABAQUS 사용자정의 서브루틴을 통해 파손 분포 (UVARM1-UVARM4)를 손쉽게 파악할 수 있다. 그림의 적층 복합재료에서 파손이 일어난 부분, 즉, 파손지수가 1.0 이상인 영역은 회색으로 표시함으로써 복합재료 구조물의 파손여부를 쉽게 분석할 수 있다.
특히, 적층 복합재료에 섬유방향·섬유수직방향 이축하중이 동시에 작용하거나 섬유방향 압축하중과 섬유방향·섬유수직방향 전단하중이 동시에 작용하는 경우, 섬유파손은 발생하지 않는데 반해 기지파손이 발생하는 것을 관찰할 수 있다. 이는 기존 파손기준식으로는 판별할 수 없는 파손모드이며 따라서 보다 개선된 파손기준식의 필요성을 잘 보여주는 예라 할 수 있다.
Table 7은 제안 해석기법으로 계산된 Hashin·Puck 섬유 및 기지파손지수와 이의 문헌 결과 (Soden, et al., 1998a; 1998b)와의 비교를 나타내고 있다. Table 7에 보인 바와 같이, 본 연구를 통해 수행된 전산해석 결과가 문헌 결과와 79.2~91.7% 수준으로 일치하였으며, 이는 공학적으로 근접한 경향을 나타낸다고 볼 수 있다. 다만, 몇몇 Case에서 비파손을 파손으로 평가하는 경우가 있었으나 (Hashin의 경우, Case A-5, B-1, B-3, B-5, B-7, Puck의 경우, Case A-5, B-7) 보수적 설계 관점에서는 크게 문제되지 않는 결과라 판단된다.
[Table 7] Failure analysis results and their comparison to literature results
Failure analysis results and their comparison to literature results
본 연구에서는 Hashin 및 Puck 파손이론을 기반으로 적층 복합재료의 섬유·기지파손 평가용 수치해석 기법을 개발하였다. 나아가 파손 시나리오 전산해석 결과를 문헌으로 보고된 초기파손응력 포락선과 비교하였다.
본 연구의 주요한 결론을 기술하면 다음과 같다.
● 적층 복합재료의 정적파손기준으로 널리 활용되던 한계 기준식 및 상호작용 기준식을 보완할 수 있는 진보된 파손기준식을 활용하여 범용유한요소해석코드용 전산해석기법을 개발하였다
● 문헌으로 보고된 적층 복합재료의 초기파손응력 포락선과 의 비교를 통해 제안 전산해석기법을 검증하였으며 공학적으로 근접한 경향의 결과를 얻을 수 있었다.
● 본 연구를 통해 검증된 전산해석 기법 및 ABAQUS 사용자정의 서브루틴의 활용을 바탕으로 보다 복잡한 형상을 가지는 복합재료 구조물의 정적파손평가, 특히, 기지파손평가를 효과적으로 수행할 수 있을 것으로 예상된다.
● 풍력용 복합재료 블레이드, 상선용·특수선용 복합재료 프로펠러 등 산업계에서 그 개발이 활발히 이루어지고 있는 복합재료 구조물의 선급규정 기반 정적파손평가를 위한 유용한 도구로 활용될 수 있을 것으로 기대된다.
