점탄성의 기억과 에너지 소산

같은 변형에서도 응력이 달라지는 이유

고무줄을 빠르게 늘린 다음 길이를 유지해보자. 길이는 거의 같아도 필요한 힘이 시간에 따라 달라질 수 있다. 천천히 늘렸을 때의 힘 기록도 다를 수 있다. 현재 변형률 하나만으로 응력을 정하는 관계는 이런 관측을 모두 설명하지 못한다. 재료가 과거 변형의 어떤 정보를 기억하는지 구성방정식에 포함해야 한다.

여기서는 스프링과 점성 요소인 대시폿을 두 방식으로 연결한다. 이들은 작동 원리를 이해하는 기본 모델이며 모든 고분자나 조직의 보편적 설명은 아니다. 스프링은 회수 가능한 기계 에너지를 저장한다. 대시폿은 기계 일을 모델에 명시하지 않은 열적 저장소로 넘긴다. 따라서 응력의 변화가 곧 에너지 소멸이나 능동 재료의 증거는 아니다.

가정과 단위, 지정한 실험

공간적으로 균일한 작은 일축 공학 변형률을 고려한다. 시편의 운동을 푸는 대신 변형 이력을 외부에서 지정한다. 탄성계수 E와 점도 eta는 양의 상수다. 온도 의존성, 소성, 관성, 공간 기울기, 큰 변형은 제외한다. 전체 변형률과 내부 탄성 변형률은 무차원이며 작은 범위로 제한한다. 응력 단위는 Pa, 점도는 Pa s, 에너지 밀도는 J/m³다.

$$\varepsilon(t)=\varepsilon_0\sin(\omega t),\qquad \dot\varepsilon=\varepsilon_0\omega\cos(\omega t),\qquad [E]=\mathrm{Pa},\quad[\eta]=\mathrm{Pa\,s}.$$

진폭의 부호를 바꾸면 변형 방향의 약속이 뒤집힌다. 각주파수는 음수가 아니다. 이 실험은 전체 변형률0에서 시작하되 초기 변형률 속도를 지정한 매끄러운 이력이다. 순간적으로 변형률을 뛰게 하는 계단 입력과는 다르다. 입력의 모양과 초기 내부 상태를 먼저 정해야 응력 변화의 원인을 올바르게 비교할 수 있다.

병렬 연결에서 Kelvin–Voigt를 유도한다

병렬인 두 요소에는 같은 변형률이 걸리고 응력을 더한다. 스프링 응력은 변형률에, 대시폿 응력은 변형률 속도에 비례한다. 둘을 합하면 현재의 두 양으로 결정되는 순간 관계를 얻는다.

$$\sigma_{KV}=E\varepsilon+\eta\dot\varepsilon.$$

사인 실험의 시작점에서 변형률은0이지만 속도는0일 필요가 없다. 따라서 Kelvin–Voigt 응력은 처음부터 양수나 음수일 수 있다. 숨은 기억이 갑자기 나타난 것이 아니라 대시폿이 현재 속도에 반응하는 것이다. 이상적인 불연속 변형률 계단은 대시폿에 충격 응력을 요구하므로 같은 식에서 유한한 보통 완화 곡선으로 그리면 안 된다. 매끄러운 시작 실험과 불연속 하중의 수학적 이상화를 구분하자.

직렬 연결에서 Maxwell을 유도한다

직렬인 두 요소는 같은 응력을 받고 변형률을 더한다. 스프링 변형률은 응력을 E로 나눈 값이고, 대시폿 변형률 속도는 응력을 eta로 나눈 값이다. 전체 변형률을 미분하면 응력의 일차 진화방정식을 얻는다.

$$\dot\varepsilon=\frac{\dot\sigma}{E}+\frac\sigma\eta, \qquad \dot\sigma+\frac\sigma\tau=E\dot\varepsilon, \qquad \tau=\frac\eta E.$$

완화 시간 tau는 점성과 탄성의 상대적인 시간척도다. Kelvin–Voigt와 달리 현재 전체 변형률과 속도만으로 Maxwell 응력을 정할 수 없다. 초기 스프링 변형률 같은 내부 상태도 필요하다. 처음 전체 변형률이0인데 스프링이 미리 변형되어 있다면 대시폿에는 반대 부호의 초기 변형률이 있다. 전체 길이가 같다는 사실만으로 내부 에너지까지 같다고 단정할 수 없는 이유다.

초기 기억을 명시한다

초기 탄성 변형률을 e-zero라 하자. 동차항은 초기 기억을 전달하고 적분항은 이후의 하중 반응을 전달한다. 지수 커널은 오래된 변형률 속도보다 최근의 이력에 더 큰 가중치를 준다.

$$\sigma_0=Ee_0,\qquad \sigma(t)=\sigma_0e^{-t/\tau} +E\int_0^t e^{-(t-s)/\tau}\dot\varepsilon(s)\,ds.$$

E와 eta, 현재 지정한 변형이 같아도 초기 내부 상태가 다른 시편의 응력은 다를 수 있다. 두 응력의 차이는 지수적으로 감소한다. 초기 기억항을 제거하면 실제 장 자체가 바뀐다. 반면 실제 장을 유지하고 방정식에서 완화항을 누락하는 대조는 다른 질문이다. 첫 비교는 초기조건의 역할을, 둘째 비교는 방정식 항의 필요성을 검사한다.

시작 과도와 정상 주기 반응을 분리한다

사인 이력에 대해 무차원 주파수 x를 정의하자. 직접 적분하면 두 주기 계수와 감쇠 과도를 얻는다. a와 b에는 변형률 진폭이 들어 있으므로 탄성률이 아니라 응력 단위의 계수다.

$$x=\omega\tau,\qquad a=\frac{E\varepsilon_0x^2}{1+x^2},\qquad b=\frac{E\varepsilon_0x}{1+x^2},\qquad \sigma=a\sin(\omega t)+b\cos(\omega t)+(\sigma_0-b)e^{-t/\tau}.$$

순수한 주기 상태의 초기 응력은 b다. 초기 응력을0으로 두면, 선언된 초기 응력에서 오는 동차 기여가0이어도 일반적으로 시작 과도는 남는다. 따라서 초기 기억0과 이미 정상 주기 상태라는 말은 다르다. 실제 응력, 주기 응력, 과도를 따로 그려본 뒤 응력–변형률 경로를 해석하자. 처음부터 정상 주기 공식을 사용하면 시작할 때의 에너지 저장을 놓칠 수 있다.

저장 에너지와 소산, 기계 전력

Kelvin–Voigt의 스프링은 전체 변형률에 따라 에너지를 저장한다. Maxwell에서는 전체 변형률 중 스프링 몫만 회수 가능한 에너지를 저장한다. 대시폿의 소산률은 두 경우 모두 음수가 아니다.

$$\psi_{KV}=\frac12E\varepsilon^2,\quad \mathcal D_{KV}=\eta\dot\varepsilon^2,\qquad \psi_M=\frac{\sigma^2}{2E},\quad \mathcal D_M=\frac{\sigma^2}{\eta}.$$

실제 에너지를 미분하고 구성식을 사용하면 국소 전력 수지를 얻는다. 스프링이 저장 에너지를 되돌려줄 때 입력 전력은 음수일 수 있다. 수동성의 비음 조건은 모든 순간의 입력 전력이 아니라 소산에 적용된다.

$$\sigma\dot\varepsilon=\dot\psi+\mathcal D, \qquad [\sigma\dot\varepsilon]=[\mathcal D]=\mathrm{W/m^3}.$$

소산 에너지는 기계적 부분계에서 빠져나간다. 여기서는 시편의 온도, 열용량, 열 경계 교환을 계산하지 않는다. 기계 소산을 열 생성으로 해석하려면 추가 물리 설정이 필요하며, 이 구성식만으로 완전한 열방정식을 얻은 것은 아니다. 에너지 밀도와 전력을 같은 숫자처럼 비교하지 말고 J/m³와 W/m³를 구분하자.

실제 관측 시간 구간을 적분한다

유한한 실험에서는 지정한 시간 구간의 일을 기록한다. 실제 전력과 실제 소산을 따로 적분하고 저장 변화와 비교한다. 물리적 기준 체적을 곱할 때도 모든 밀도 항에 같은 체적을 사용하고 그 변환의 가정을 확인해야 한다.

$$W_{[t_0,t_1]}=\int_{t_0}^{t_1}\sigma\dot\varepsilon\,dt =\psi(t_1)-\psi(t_0)+\int_{t_0}^{t_1}\mathcal D\,dt.$$

시작 실험의 첫 주기는 보통 처음과 끝의 응력이 다르다. 응력–변형률 평면의 경로가 닫히지 않으며 저장 에너지 변화도0일 필요가 없다. 공유한 해석식에서 계산한 작은 수지 잔차는 진단이다. 구현을 검증하려면 별도의 미분, convolution, 시간 적분이 필요하다. 특히 모든 항을 잔차0이 되도록 정의하는 방식은 잘못된 방정식을 검출하기 어렵다.

저장·손실 탄성률은 주기 상태를 기술한다

사인 변형률을 위상 기준으로 삼자. 같은 위상의 계수는 저장 탄성률, cosine 계수는 손실 탄성률을 정의한다. 복소 지수의 부호 약속만 바꾸고 위상 해석을 그대로 두면 겉보기 손실 부호가 뒤집힐 수 있다.

$$\sigma_{\rm steady}=\varepsilon_0[E'\sin(\omega t)+E''\cos(\omega t)],\quad E'_{KV}=E,\ E''_{KV}=\eta\omega,\quad E'_M=\frac{Ex^2}{1+x^2},\ E''_M=\frac{Ex}{1+x^2}.$$

주파수가 양수인 정상 주기에서는 저장 에너지가 원래 값으로 돌아오므로 한 주기의 일이 소산과 같다. 이는 체적당 에너지이며 전력이나 무차원 도형 면적이 아니다.

$$W_{\rm steady\ cycle}=\oint\sigma\,d\varepsilon=\pi\varepsilon_0^2E''\ge0.$$

주파수0에서 지금의 사인 이력은 항등적으로0이므로 유한한 주기가 없다. 미리 변형된 Maxwell 스프링은 외부 전력0에서도 완화할 수 있다. 아주 작은 양의 주파수는 주기를 교육용 관측 시간이나 수치 표현 범위 밖으로 보낼 수 있다. 둘을 주파수0의 유한한 주기로 바꾸어 표시하면 안 된다.

손으로 계산하는 기준 예제

E=2Pa, eta=1Pa s, 진폭0.002, 각주파수1rad/s, 초기 탄성 변형률0을 택하자. 계산하기 쉬운 교육용 수치이며 특정 재료에 맞춘 계수는 아니다.

$$\tau=0.5\ \mathrm{s},\quad x=0.5,\quad a=0.0008\ \mathrm{Pa},\quad b=0.0016\ \mathrm{Pa},\quad E'_M=0.4\ \mathrm{Pa},\quad E''_M=0.8\ \mathrm{Pa}.$$

Maxwell 응력은 Pa 단위로 0.0008 sin(t)+0.0016 cos(t)−0.0016 exp(−2t)다. 정상 주기 손실은 약1.0053×10의−5승 J/m³이다. 실제 첫 주기 일은 약8.7731×10의−6승 J/m³이며 그 일부가 끝점에서 저장된 채 남는다. 같은 실험의 Kelvin–Voigt 저장·손실 탄성률은 각각2Pa와1Pa다. 같은 계수를 사용해도 요소의 연결과 내부 상태에 따라 전력 분배가 달라진다.

작은 변화와 상쇄, 3차 시작 변화

거의 같은 현재 응력에서 초기 응력을 직접 빼면 실제 변화를 잃을 수 있다. 초기 응력 변화율과 구동 cosine이1에서 벗어나는 몫을 분리하면 더 안정적으로 계산할 수 있다. 초기 응력을 eta 곱하기 초기 변형률 속도로 맞추면 1차 미분이 사라지고 사인 원점에서는 2차 미분도 사라진다.

$$\sigma_0=\eta\varepsilon_0\omega \quad\Longrightarrow\quad \sigma(t)-\sigma_0=-\frac{E\varepsilon_0\omega^3}{6}t^3+O(t^4).$$

표시 응력이 그대로여도 이 3차 변화는 비영이다. expm1 형태, 부호 있는 로그, 정규화한 적분 좌표를 사용하면 반올림된 끝점 뺄셈에서 사라지는 정보를 보존할 수 있다. 선형 표현 범위 아래의 비영 관측은 로그 크기와 함께 제공 불가로 표시하고 정확한0으로 강제하지 않는다. 수치 주기 적분에서도 큰 탄성 기여가 상쇄되며 아주 작은 손실을 가릴 수 있으므로 해석적인 정상 주기 식과 비교해야 한다.

3차라는 결론은 화면의 소수뿐 아니라 실제 입력에서 초기 등식이 성립해야 한다. 예를 들어 진폭0.001과 각주파수0.7은 이진 부동소수점으로 저장되며, 두 값의 정확한 곱은 따로 입력한 초기 탄성 변형률0.0007과 같지 않을 수 있다. 탄성계수와 점도가 각각1이면 초기 응력 변화율은 약−2.27×10의−20승 Pa/s다. 구동 항 둘을 각각 먼저 반올림하면 이 작은 차이를 지우고3차 변화만 남기는 오류가 생긴다. 구현에서는 시간0도 포함해 입력의 차와 비를 유지한 뒤 초기 변화율을 관측한다. 정확히 조율한3차 실험과 작은1차 기여가 남은 인접 실험을 구분하자. 표시 자릿수나 같은 식을 공유한 작은 잔차만으로 초기 등식이 성립한다고 판단하지 않는다.

실험 순서와 풀이

먼저 같은 하중에서 시작 상태와 정상 주기 초기 상태를 비교하자. 다음에는 주파수를 바꾸며 응력 위상, 저장·손실 탄성률을 관찰하자. 마지막으로 실제 시간 구간의 일·저장 변화·소산을 함께 읽자. 학습된 구성식 연구에서는 내부 상태 누락, 잘못된 위상 약속, 낮은 적합 오차 뒤에 숨은 에너지 위반을 이런 비교로 찾을 수 있다.

한 주파수에서 위상과 진폭을 모두 측정해도 모델 구조를 구분하기에는 부족할 수 있다. 적절한 양의 계수를 택하면 두 기본 모델 모두 그 주파수의 양의 저장·손실 쌍을 재현할 수 있기 때문이다. 여러 주파수와 시작 과도·완화 이력은 추가적인 판별 정보를 준다. 한 곡선을 맞춘 다음에는 초기 상태, 단위, 응력 정의를 명시한 채 다른 하중 이력에서 예측을 검증해야 한다. 이는 계수 추정과 모델 구조의 식별을 구분하는 예다.

연습 1 — 병렬과 직렬

각 연결에서 두 요소가 공유하는 양은 무엇인가?

해설: 병렬은 변형률을 공유하고 응력을 더한다. 직렬은 응력을 공유하고 변형률을 더한다. E와 eta가 같아도 이 기하·힘 조건이 달라 서로 다른 구성방정식을 얻는다.

연습 2 — 초기 응력

처음 전체 변형률이0인데 Maxwell 응력이 비영일 수 있는가?

해설: 가능하다. 스프링 변형률이 e-zero, 대시폿 변형률이 그 음수이면 합은0이다. 스프링 응력은 E e-zero이므로 내부 초기조건을 명시해야 한다.

연습 3 — 음의 전력

기계 전력이 음수이면 능동 모델이라는 증거인가?

해설: 아니다. 스프링이 저장 에너지를 되돌려줄 수 있다. 소산의 비음성과 실제 구간에서 입력 일=저장 변화+소산 적분이라는 수지를 함께 확인한다.

연습 4 — 첫 주기의 일

초기 응력0에서 시작한 첫 주기 일에 정상 주기 식을 바로 써도 되는가?

해설: 일반적으로 안 된다. 과도가 남아 끝점 응력이 초기값과 다르고 저장 에너지가 변한다. 실제 응력장으로 일을 적분해야 한다.

연습 5 — 주파수0

미리 변형된 Maxwell 모델에 이 실험의 주파수0을 주면 어떻게 되는가?

해설: 지정한 변형률과 속도는0이다. 초기 응력은 지수 완화하고 저장 에너지는 줄며 소산은 양수다. 주기 손실을 붙일 유한한 주기는 없다.

연습 6 — 숨은 3차 변화

처음 응력 변화율이0인데 작은 양의 시간에 응력이 바뀌는 이유는?

해설: 1차 미분0이 함수 전체의 상수를 뜻하지 않는다. 조율한 초기조건에서는 선도 변화가3차다. 작은 변화를 반올림 응력의 차이와 별도로 보존해야 한다.

자료와 한계

A. F. Bower의 Applied Mechanics of Solids §3.6은 작은 변형의 점탄성, 하중 이력과 기본 기계 모델을 설명한다. 여기서는 명시한 일축 약속으로 매끄러운 이력과 에너지 관계를 직접 유도했다. 두 요소 예제는 모든 재료의 보정된 모델이나 일반적인 큰 변형 점탄성 법칙이 아니다.