Nonlinear interface effects in multilayered structures: vibro-acoustic modeling and experimental analysis

1. 연구 배경 및 문제 정의 (Problem)

• 배경: 운송, 토목, 전자, 생체 의료 등 다양한 분야에서 다층 구조물 (적층판, 샌드위치 패널 등) 이 널리 사용되고 있습니다. 이러한 구조물의 동적 거동은 층 간의 인터페이스 (접합부) 상태에 크게 의존합니다.
• 문제점: 실제 제조 공정에서는 접착제 결함, 부분적 박리, 얇은 유연한 중간층 등으로 인해 불완전한 인터페이스 (imperfect interfaces) 가 발생합니다.
• 기존 연구의 한계: 기존 모델들은 주로 선형적인 스프링 - 슬라이딩 인터페이스를 가정하거나, 저주파수 영역에서의 거동만 설명하는 경우가 많습니다. 그러나 진동 레벨 (excitation level) 에 따라 인터페이스의 기계적 특성이 변하는 비선형 거동을 정량적으로 분석하고 모델링하는 방법은 아직 미흡합니다.
• 목표: 다양한 진동 진폭 하에서 다층 구조물의 비선형 인터페이스 거동을 실험적으로 규명하고, 이를 설명할 수 있는 진폭 의존성 (amplitude-dependent) 을 가진 등가 인터페이스 파라미터를 식별하는 방법론을 개발하는 것.

2. 방법론 (Methodology)

가. 이론적 모델링 (Theoretical Modeling)

• Zig-Zag 형식화 (Zig-Zag Formulation): Guyader 와 Lesueur 가 제안한 모델을 기반으로 하여, 각 층을 독립적인 운동학 필드로 기술하면서도 인터페이스 조건을 통해 층 간 결합을 수행합니다. 이는 자유도를 줄이면서도 각 층의 운동을 독립적으로 묘사할 수 있게 합니다.
• 불완전 인터페이스 모델: 인터페이스를 선형 스프링 층으로 모델링하여, 인접 층 간의 전단 응력 연속성과 변위 불연속성 (jump) 을 연결합니다.
  • 비선형 인터페이스 거동을 처리하기 위해 조각별 선형 근사 (piecewise linear approximation) 를 사용했습니다.
  • 인터페이스의 비선형성을 나타내는 등가 인터페이스 파라미터 $B$ (인터페이스 순응도, compliance) 를 도입했습니다. $B=0$은 완전 결합, $B \to \infty$는 완전 박리를 의미합니다.
• 등가 Kirchhoff-Love 판 모델: 다층 구조를 주파수에 의존하는 등가 굽힘 강성 ($D$) 을 가진 단일 등가 판으로 변환하여 전체 구조의 동적 응답을 유도했습니다.
• 시그모이드 기반 모델 (Sigmoid-based Model): 굽힘 강성 $D(f)$ 의 주파수 의존성을 명확한 식으로 표현하기 위해 Arasan 등이 제안한 시그모이드 함수를 적용했습니다.
  • 저주파수 한계 ($D_{low}$) 와 고주파수 한계 ($D_{high}$) 사이의 전이 영역을 정의하며, 이 전이 주파수 ($f_T$) 와 기울기 ($R$) 를 통해 인터페이스 파라미터 $B$의 영향을 정량화합니다.

나. 실험적 분석 (Experimental Analysis)

• 시편: 유리/에폭시/유리 (Glass/Epoxy/Glass) 3 층 빔 구조물 (두께: 유리 3mm, 에폭시 0.35mm).
• 실험 장치:
  • 진동 여기: 양단에 부착된 압전 버저 (Piezoelectric buzzers) 를 사용하여 1 kHz ~ 50 kHz 대역의 선형 주파수 칩 (chirp) 신호를 인가.
  • 측정: 스캐닝 레이저 진동계 (Scanning Laser Vibrometer) 를 사용하여 빔의 횡방향 변위장 (transverse displacement field) 을 측정.
  • 조건: 7 가지 다른 진폭 레벨 ($A_0$: 0.1V ~ 3.0V) 에서 데이터를 수집.
• 데이터 처리:
  • CFAT (Corrected Force Analysis Technique): 측정된 변위장을 기반으로 등가 굽힘 강성 ($D$) 을 역산출하는 방법.
  • 2D FFT 및 피팅: 측정된 변위장에서 종방향 굽힘 모드 (bending mode) 의 파수 ($k_b$) 를 분리하고, 이를 시그모이드 모델에 피팅하여 실험적 전이 주파수 ($f_T^{exp}$) 를 도출.
  • 파라미터 식별: 도출된 $f_T^{exp}$와 이론적 모델 ($f_T(B)$) 을 비교하여 각 진폭 레벨에 해당하는 인터페이스 파라미터 $B$를 수치적으로 역해결 (root-finding) 함.

3. 주요 기여 (Key Contributions)

1. 비선형 인터페이스 파라미터의 실험적 식별: 진동 진폭의 변화에 따라 인터페이스 순응도 파라미터 $B$가 어떻게 변하는지를 최초로 정량적으로 규명했습니다.
2. 시그모이드 기반 모델의 적용: 복잡한 다층 구조의 주파수 의존적 굽힘 강성 거동을 간결한 시그모이드 함수로 모델링하여, 실험 데이터로부터 인터페이스 특성을 효율적으로 추출하는 방법론을 제시했습니다.
3. CFAT 와 이론 모델의 결합: 실험적으로 측정된 동적 특성을 이론적 모델과 직접 비교하여, 인터페이스의 비선형성이 전체 구조의 동적 특성 (특히 전이 주파수) 에 미치는 영향을 명확히 입증했습니다.

4. 주요 결과 (Results)

• 진폭 의존성 확인: 진동 진폭이 증가함에 따라 식별된 인터페이스 파라미터 $B$가 약 2 배까지 증가하는 것을 확인했습니다. 이는 진동 레벨이 커질수록 인터페이스의 결합이 약화되거나 (슬립 증가), 비선형적 거동이 강화됨을 의미합니다.
• 시스템 연화 (Softening) 현상: $B$의 증가는 등가 굽힘 강성 ($D$) 및 동적 영률 ($E$) 의 감소로 이어집니다. 특히 전이 주파수 영역 ($f_T$) 에서 주파수 하향 이동 (약 3 kHz 감소) 이 관측되었으며, 이는 시스템이 진동 레벨이 높아질수록 '연화 (softening)'되는 비선형 거동을 보임을 의미합니다.
• 모델 정확도: 제안된 시그모이드 기반 모델은 Guyader 의 참조 모델과 매우 높은 정확도 (상대 오차 최대 1.2%) 로 일치하며, 실험적으로 측정된 CFAT 결과와도 7 kHz ~ 30 kHz 대역에서 잘 부합했습니다.
• 주파수 영역별 영향: 저주파수 영역에서는 인터페이스 품질의 영향이 미미하나, 고주파수 영역 및 전이 영역에서는 인터페이스 파라미터 변화가 구조 동적 특성에 큰 영향을 미쳤습니다.

5. 의의 및 결론 (Significance & Conclusion)

• 학술적 의의: 제조 결함으로 인한 불완전 인터페이스가 다층 구조물의 등가 물성 (굽힘 강성, 영률) 에 진폭 의존적인 비선형 변화를 일으킨다는 사실을 실험적으로 증명했습니다. 이는 기존 선형 모델로는 설명할 수 없었던 현상을 규명한 것입니다.
• 실용적 의의:
  • 다층 구조물의 동적 설계 시, 작동 진동 레벨에 따른 인터페이스 특성의 변화를 고려해야 함을 시사합니다.
  • 제안된 방법론은 접착 품질 평가, 비파괴 검사, 그리고 3D 프린팅 등 첨가제 제조 기술을 통한 정밀한 인터페이스 제어 구조물의 설계에 활용될 수 있습니다.
• 향후 과제: 다양한 접착 품질을 가진 시편을 추가 테스트하여 인터페이스 파라미터와 접착 품질 간의 상관관계를 규명하고, 능동 제어 가능한 인터페이스를 가진 다층 구조물 설계로 연구 범위를 확장할 계획입니다.

핵심 요약: 이 연구는 유리/에폭시/유리 3 층 빔을 대상으로, 진동 진폭이 커질수록 인터페이스의 비선형성 (슬립 증가) 이 강화되어 구조의 등가 강성이 감소하고 전이 주파수가 낮아지는 현상을 실험 및 이론적 모델 (시그모이드 함수 기반) 을 통해 정량적으로 규명했습니다. 이는 다층 구조물의 동적 거동 예측 및 설계에 있어 인터페이스 비선형성 고려의 중요성을 강조합니다.