응력 특이점이란 — 메시를 줄여도 응력이 수렴하지 않는 이유
유한요소해석에서 메시를 세밀하게 나눌수록 특정 지점의 응력이 계속 증가하는 경우가 있습니다. 이는 해석 오류가 아니라 응력 특이점(Stress Singularity)이라는 수학적 현상입니다. 이를 구분하지 못하면 안전한 구조를 불합격 처리하게 됩니다.
현상
일반적인 FEA에서는 메시를 세밀화할수록 결과가 참값에 수렴합니다. 요소 크기를 절반으로 줄였을 때 최대응력 변화가 몇 퍼센트 이내라면 수렴했다고 판단합니다.
그런데 어떤 지점에서는 메시를 아무리 줄여도 응력이 계속 올라갑니다. 요소 크기를 1/2로 줄이면 응력이 1.3배, 다시 1/2로 줄이면 또 1.3배. 수렴하지 않고 발산합니다.
원인
이는 해석 모델이 물리적으로 불가능한 형상이나 조건을 포함하고 있기 때문입니다. 탄성 이론상 해당 지점의 응력이 무한대가 되도록 정의되어 있어, 요소를 아무리 작게 만들어도 유한한 값으로 수렴할 수 없습니다.
주요 발생 원인
| 원인 | 설명 |
|---|---|
| 날카로운 내부 모서리 | 필렛 반경이 0인 재진입 코너. 실제 부품에는 반드시 가공 반경이 존재하지만 CAD 모델에서 생략된 경우. |
| 점 하중 | 면적이 0인 점에 유한한 힘을 가하면 응력 = 힘 / 0 이 되어 무한대가 됩니다. |
| 점 구속 | 단일 절점을 고정하는 경계조건. 반력이 한 점에 집중되어 동일한 문제가 발생합니다. |
| 재료 물성 급변 경계 | 강성이 크게 다른 두 재료가 만나는 모서리. |
진단 방법
최대응력 지점이 특이점인지 아닌지를 판별하는 절차는 다음과 같습니다.
- 해당 부위의 메시를 단계적으로 세밀화합니다 (예: 요소 크기 1.0 → 0.5 → 0.25 → 0.125)
- 각 단계의 최대응력을 기록합니다
- 응력 변화 추이를 확인합니다
수렴하면 실제 응력 집중이며, 그 값을 판정에 사용할 수 있습니다.
계속 증가하면 특이점이며, 해당 값은 물리적 의미가 없습니다.
특이점의 응력값은 "틀린 값"이 아니라 "의미 없는 값"입니다
수치 자체는 그 메시 조건에서 정확히 계산된 값입니다. 다만 메시를 바꾸면 값이 달라지므로, 구조의 안전성을 판정하는 기준으로 쓸 수 없다는 뜻입니다.
대응 방법
1. 형상 수정
가장 근본적인 방법입니다. 실제 부품에 존재하는 필렛 반경을 모델에 반영하면 특이점이 사라지고 응력이 수렴합니다. 실제 가공 반경을 알 수 있다면 우선 적용합니다.
2. 하중·구속 조건 개선
점 하중을 면 하중으로, 점 구속을 면 구속 또는 Remote Displacement로 변경합니다. 실제 하중 전달 방식에 가깝게 모델링하는 것이 원칙입니다.
3. 응력 선형화
형상 수정이 어려운 경우 적용합니다. 단면을 가로지르는 경로를 설정하고 그 경로를 따라 응력을 막응력(membrane)과 굽힘응력(bending) 성분으로 분해합니다. 특이점의 국부 피크는 이 분해 과정에서 제거되며, 남은 성분을 해당 기준의 허용값과 비교합니다.
이 방법은 압력용기 및 원자력 기기 평가에서 표준적으로 사용되며, 관련 코드에 절차가 규정되어 있습니다.
4. 평가 위치 조정
특이점에서 일정 거리 떨어진 지점의 응력을 평가값으로 사용하는 방법입니다. 다만 거리 선정의 근거를 명확히 제시해야 하며, 임의로 정하면 타당성을 인정받기 어렵습니다.
실무 권고
- 최대응력이 허용값을 초과하면, 먼저 그 지점이 특이점인지 확인합니다
- 메시 수렴성 검토 결과를 보고서에 포함합니다
- 응력 선형화를 적용한 경우 경로 설정 근거를 명시합니다
- 특이점을 근거로 설계를 변경하기 전에 진단을 먼저 수행합니다
불필요한 설계 보강은 비용과 중량 증가로 이어집니다. 특이점 여부를 정확히 판별하는 것만으로도 상당한 설계 여유를 확보할 수 있습니다.