[CIVIL NX] 수화열 해석에서 시간이산계수는 어떻게 설정하나요? – MIDAS Support

질문

수화열 해석에서 시간이산계수(Temporal discretization Factor, Φ)는 어떻게 설정하나요?

답변

수화열 해석은 시간에 따라 변화하는 콘크리트의 온도를 계산하기 위해 해석 시간을 여러 개의 시간 단계(Time Step)로 나누어 순차적으로 해석합니다.
이때 시간 이산 계수(Φ)는 현재 시간 단계와 다음 시간 단계의 상태를 해석에 어느 정도 반영할 것인지를 결정하는 계수입니다. 예를 들어 현재 온도가 20℃, 다음 시간 단계의 온도가 30℃라고 가정하면 다음과 같이 이해할 수 있습니다.

Φ = 0.0 → 현재 단계의 온도인 20℃를 기준으로 계산
[Forward difference, 1 × 20℃ + 0 × 30℃ = 20℃ ]

Φ = 0.5 → 현재와 다음 단계 온도의 평균인 25℃를 기준으로 계산
[Crank-Nicolson method, 0.5 × 20℃ + 0.5 × 30℃ = 25℃ ]

Φ = 1.0 → 다음 단계의 온도인 30℃를 기준으로 계산
[Backward difference, 0 × 20℃ + 1 × 30℃ = 30℃ ]

위의 20℃, 25℃, 30℃는 시간 이산 계수의 개념을 쉽게 이해하기 위해 온도만을 이용하여 단순화한 예시입니다.

실제 수화열 해석에서는 온도 자체를 단순히 가중 평균하는 것이 아니라, 현재 시간 단계와 다음 시간 단계의 열전달 방정식에 각각 가중치를 적용하여 계산합니다.

여기서 Crank–Nicolson method(크랭크–니콜슨 방법)는 시간에 따른 온도 변화를 계산하기 위해 사용하는 대표적인 시간 적분법(Time Integration Method)입니다. 현재 시간 단계와 다음 시간 단계의 열전달 방정식을 동일한 비율인 50:50으로 고려하여 다음 시간 단계의 해를 계산하는 방법입니다.

Φ 값이 증가할수록 다음 시간 단계의 상태를 더 많이 반영하는 방식으로 해석이 수행됩니다.

일반적으로 Φ가 0에 가까울수록 현재 시간 단계의 상태를 중심으로 계산하고, Φ가 1에 가까울수록 다음 시간 단계의 상태를 중심으로 계산합니다. 특히 Φ = 1.0은 수치적으로 안정적인 특성이 있어 시간에 따라 온도 변화가 지속적으로 발생하는 수화열 해석에서 일반적으로 적용됩니다.

CIVIL NX HYPER-S에서는 수치적으로 안정적인 특성이 있는 Φ = 1.0를 고정값으로 설정하고 있습니다.