Data Loading...
점탄성 및 점소성(크리프)에 대한 이해와 활용 Flipbook PDF
점탄성 및 점소성(크리프)에 대한 이해와 활용
396 Views
164 Downloads
FLIP PDF 1.49MB
점탄성, 점소성(크리프)의 재료 정의 및 활용 방법 황정필 | 태성에스엔이
Contents
• 점탄성(Viscoplasticity) 정의 • 점소성(Viscoplasticity) 정의 • 크리프(Creep) 정의
• 크리프(Creep) 이론 • 크리프(Creep) 이론 모델 • Ansys WB에서 크리프 재료 활용 • APDL에서의 Curve Fitting
점탄성(Viscoelasticity) 정의 점탄성(Viscoelasticity) 외부의 하중에 의해 고체 성질의 탄성과 액체 성질의 점성이 동시에 나타나는 현상 –
탄성(Elasticity) 외부 하중을 제거하였을 때 원래의 형태로 완전히 복원되는 성질
𝜎 = 𝐸𝜀 𝐹 = −𝑘𝑥 –
𝜎 : 수직 응력 𝐸 : 탄성계수 𝜀 : 수직 변형률
훅(Hook)의 법칙
점성 (Viscosity) 외부 하중이 형상 변형에 모두 사용되어 복원되지 않는 성질
𝜏 = 𝜂𝛾‘ 𝑑𝛾 𝜏 = 𝜂 𝑑𝑡
𝜏 : 전단 응력 𝜂 : 점성계수 𝛾 : 전단 변형률
뉴턴(Newton)의 법칙
고체 이상 고체 (완전 탄성)
→ 변형이 시간 지연 없이 바로 반응 → 변형된 비율만이 중요함(𝜀) → 변형이 에너지로 저장됨
액체 실제 물체(점탄성)
→ 짧은 시간 : 탄성 성질 중요 → 긴 시간 : 점성 성질 중요
이상 액체 (완전 점성)
→ 변형의 시간 변화가 중요 시간 지연이 있고 장시간의 변형에서 두드러짐(𝛾‘) → 변형이 에너지로 저장 안됨
점탄성(Viscoelasticity) 정의 재료의 거동 비교 하중 조건
𝜎 𝜎0
Loading 𝝈
탄성 (Elasticity) 𝜀
= 𝝈𝟎
𝜀0 Unloading
𝝈 = 𝑬𝜺
𝝈=𝟎 t1
Maxwell model Generalize Maxwell model Kelvin model 표준 선형 model (SLM)
t
t1
t
Spring과 Dashpot의 연결 상태에 따라 점탄성에서 사용하는 Equation이 달라짐 • • • •
※ 동적 특성 고려 안함
점성 (Viscosity) 𝜀 𝜀0 𝜏 = 𝜂
𝑑𝛾 𝑑𝑡 t
t1
점탄성 (Viscoelasticity) 𝜀 𝜀0
𝝈 = 𝒇(𝜺, 𝒕) t1
t
점탄성(Viscoelasticity) 정의 점탄소성의 구분 항복 변형률을 기점으로 선형 점탄성, 비선형 점탄성 영역의 재료로 구분 –
선형 점탄성 영역(Linear Viscoelastic Region; LVR) 항복 이하의 영역에서 점성과 탄성의 양면성을 지닌 재료
–
비선형 점탄성 영역(Non Linear Viscoelastic Region; Non-LVR) 점성과 탄성의 양면성을 지닌 재료로 항복을 초과하는 영역
Modulus & Stress
LVR
Non-LVR
Modulus Stress
점탄성 (Viscoelasticity)
점소성 (Viscoplasticity)
Yield Strain
Strain
보통 재료의 LVR 범위는 변형량이 0.3% 안이며, 고무와 같은 재료의 경우 3% 이상도 가능
점소성(Viscoplasticity) 정의 재료의 특성 비교 –
탄성 (Elasticity) 외부 하중이 제거되었을 때 완전하게 원래의 모양으로 되돌아오는 성질 (탄성을 Spring으로 표현, 대부분 변형량이 적음)
𝜎 = 𝐸𝜀 (Hook’s law) –
탄소성 (Elasto-Plasticity)
외부 하중의 증가에 따라 변형이 점차 탄성영역에서 소성영역으로 이행하는 성질 (탄성을 Spring, 소성을 Slider로 표현, 영구변형 존재)
–
점탄성 (Viscoelasticity) 외부 하중에 의해 고체의 탄성과 유체의 점성 특징이 동시에 나타나는 성질 (점성의 Dashpot와 탄성의 Spring을 함께 표현, 항복점 이하의 영역)
점소성(Viscoplasticity) 정의 재료의 특성 비교 –
점소성(Viscoplasticity) 외부 하중에 의해 고체의 탄소성과 유체의 점성 특징이 동시에 나타나는 성질 (점성의 Dashpot와 탄소성의 Spring과 Slider를 함께 표현, 항복점 이상의 영역)
𝒇 = 𝒎𝒙” + 𝒄𝒙’ + 𝒌𝒙 𝒇 = 𝒌(𝒕)𝒙 [𝑲] = ( 𝒆𝒍𝒂𝒔𝒕𝒊𝒄 + 𝒗𝒊𝒔𝒄𝒐 + [𝒑𝒍𝒂𝒔𝒕𝒊𝒄]) 점성에 관한 함수 적용
• 유변학에서의 점소성(Viscoplasticity)
점성과 탄성의 양면성을 지닌 재료로 항복을 초과하는 영역의 재료(비선형 점탄성 영역(Non Linear Viscoelastic Region; Non-LVR)), 재료의 탄성률이 변형률에 따라 값이 달라짐 sx
ex
점소성(Viscoplasticity) 정의 변형률 속도에 따른 거동 차이 (Fate Dependent) 소성변형에서 응력과 변형률의 관계가 변형률 속도(strain rate)에 따라 달라질 때
sx
e 1'
e
' 2
sx
e 3'
ex
e 1'
e
' 2
e e 1'
T1
' 3
e
' 2
e 3'
T2
e 1' e 2' e 3'
빨리 당긴 경우 중간 느리게 당긴 경우
T2 > T1 ex
• 점소성(Viscoplasticity)의 종류 → 크리프(CREEP) → RATE → ANAND ※ CREEP과 점소성은 Engineering의 관점에서는 차이가 있지만 재료 관점에서 보면 동일
크리프(CREEP) 정의 크리프(CREEP)와 응력완화(Stress Relaxation) 𝜀
e
하중 제거
𝜀0
𝝈 = 𝒇(𝜺, 𝒕) t1
– –
creep ①
②
recover ③ relaxation
t
t
크리프(CREEP) Creep은 그래프와 같이 일정 응력 아래서 strain이 증가함(선형이며 exponential creep을 보여줌) 응력완화(Stress Relaxation) Stress Relaxation은 그래프와 같이 일정 변형률 아래서 응력이 점진적으로 감소함
크리프(CREEP)
응력완화(Stress Relaxation)
크리프(CREEP) 정의 크리프(CREEP) 외력이 일정하게 유지되어 있을 때, 시간의 흐름에 따라 재료의 변형이 증대하는 현상으로 항복점 이하의 낮은 응력 상태에서도 파손이 발생할 수 있음 (일반적으로 고온 상태의 재료에서 많이 발생하며, 일부 재료에서 상온상태에서도 이러한 현상이 발생함) ①②③
e
1차
2차
1차 크리프
3차
기울기 ① > ② > ③ 짧은 시간 발생, 변형률 속도 감소 ③ ② ①
크리프(CREEP) 거동 함수
2차 크리프 기울기 ① < ② < ③ 파단까지 변형률 속도가 급격히 증가 Necking, Damage와 같은 파손과 연관
※ 크리프는 1차, 2차, 3차로 나눌 수 있으며 각 영역에서의 특징을 가지고 있음
크리프(CREEP) 이론 ANSYS에서 사용되는 크리프(CREEP) 이론 크리프의 거동을 정의하기 위해 변형률 속도(Strain Rate) 함수를 계산하여 재료의 거동 특성을 표현함 점점 감소
e
일정
점점 증가
e
de e dt y ax b
Time
Time
2
y ' 2ax
ecr f1 s f 2 e f 3 t f 4 T
ecr As B T C t D e E 상수 D의 값에 따라 Primary, Secondary를 표현
※ 변형률 속도(strain rate)는 재료의 종류에 따라 응력, 온도, 시간 및 변형률에 따른 함수로 결정
크리프(CREEP) 이론 ANSYS에서 사용되는 크리프(CREEP) 이론 ANSYS Mechanical에서 Elastic, creep, plastic strain은 모두 현재 응력상태에서 평가되지만, 각각 독립적으로 계산되며 서로 관계를 가지고 있지 않음
[𝑲] = ( 𝒆𝒍𝒂𝒔𝒕𝒊𝒄 + 𝒑𝒍𝒂𝒔𝒕𝒊𝒄 + 𝒗𝒊𝒔𝒄𝒐 )
e eel e pl ecr
ecr As B T C t D e E
※ 변형률 속도(strain rate)는 재료의 종류에 따라 응력, 온도, 시간 및 변형률에 따른 함수로 결정 –
Stress dependency
ecr s n
ecr e
–
Temperature-dependency
–
Time-hardening
ecr t m
–
Strain-hardening
ecr e n
Q RT
크리프(CREEP) 이론 ANSYS에서 사용되는 크리프(CREEP) 이론 –
Time-hardening term 크리프 변형률 속도는 크리프 시작 순간 부터 오직 시간에 의해서만 결정된다고 가정 (곡선이 위/아래로 이동) 응력이 σ1 에서 σ2로 변하면, creep 변형률은 점A에서 점B로 계산
ecr t –
ec
s1 s2
A
B
m
Strain-hardening term 크리프 변형률 속도가 재료의 현재 변형률 에 대해서만 좌우된다고 가정 (곡선이 좌/우로 이동) 응력이 σ1 에서 σ2로 변하면, CREEP 변형률은 점A에서 점B로 계산
ecr e n
t ec
s1 A
s2 B
t
2차 CREEP(Secondary)은 time- 이나 strain-hardening을 동반하지 않으며, 2차 CREEP에서 CREEP strain rate는 보통 일정한 값을 이용
크리프(CREEP) 이론 ANSYS에서 사용되는 크리프(CREEP) 이론 –
Temperature-dependency term CREEP 효과가 열적으로 활성화(thermally activated)되며, 온도 의존도가 보통 Arrhenius law를 통해 표현
ecr e –
Q RT
Q : 활성화에너지 (activation energy) R : 일반기체 상수 T : 절대온도
Stress dependency term CREEP strain은 보통 응력에 의존되는데, 특히 전이 CREEP과 연관, 정상상태의 CREEP 거동(2차 CREEP)은 다양한 방법으로 표현 가능
ecr s n
Norton’s law
ecr e Cs
exponential law
ecr sinh As
hyperbolic sine law
크리프(CREEP) 이론 Implicit와 Explicit Solution
Impact Velocity (m/s)
Strain Rate (/s)
Effect
12000
> 108
Vaporization of colliding solids
Implicit
Explicit
Implicit 크리프는 아래의 소성 모델(rate-independent plasticity)과 조합할 수 있음 • BISO, MISO, and NLISO with CREEP (creep with isotropic hardening) • BKIN with CREEP (creep with kinematic hardening) • HILL and CREEP (anisotropic creep) • BISO, MISO, and NLISO with HILL and CREEP (anisotropic creep with isotropic hardening) • BKIN with HILL and CREEP (anisotropic creep with kinematic hardening)
크리프(CREEP) 이론 모델 ANSYS WB에서 사용되는 크리프 이론 모델(Implicit) Creep Equation Description
Type
Strain Hardening
Primary
Time Hardening
Primary
Generalized Exponential
Primary
Generalized Graham
Primary
Generalized Blackburn
Primary
Modified Time Hardening
Primary
Modified Strain Hardening
Primary
Generalized Garofalo (Hyperbolic sine)
Secondary
Exponential Form
Secondary
Norton
Secondary
Time Hardening
Both
Rational Polynomial
Both
Generalized Time Hardening
Both
Primary
ANSYS에서 제공하는 크리프 이론 모델은 총 13개가 있으며, 각각의 이론 모델에서 사용되는 함수의 특성에 따라 크리프의 거동이 달라짐
크리프(CREEP) 이론 모델 ANSYS WB에서 사용되는 크리프 이론 모델(Implicit) –
ecr As B T C t D e E
Strain Hardening
εcr C1σ C2 ε C3 e C4 /T
–
Modified Strain Hardening
C εcr C1σ C2 C3 1ε 3
• • • •
1 C3 1 C4 /T
e
Strain-hardening term과 Norton 법칙 포함 상수 C3 는 보통 음의 값으로 e이 증가하면서 크리프 변형률 속도가 감소하게 되는 1차 크리프 표현 가능 2차 creep효과 표현가능(e이 증가함에 따라, 크리프 변형률 속도가 거의 일정할 수 있기 때문) Arrhenius equation(반응속도상수와 온도의 관계를 나타내는 식)을 포함
크리프(CREEP) 이론 모델 ANSYS WB에서 사용되는 크리프 이론 모델(Implicit) –
Time Hardening
ecr As B T C t D e E
εcr C1σ C2 t C3 e C4 /T –
Generalized Time Hardening r
εcr ft e
C6
T
f C1s C2s 2 C3s 3 r C4 C5s –
Modified Time Hardening
C1σ C2 t C3 1e C4 /T εcr C3 1 • • •
Time-hardening term과 Arrhenius equation, Norton’s law을 포함 시간에 대한 exponential 항은 1차 creep에 대해 감소하는 creep strain rate을 표현 (보통 -0.5와 -1.0 사이의 값 사용) creep strain rate이 일정한 2차 creep에서 주요 부분 모사 가능
크리프(CREEP) 이론 모델 ANSYS WB에서 사용되는 크리프 이론 모델(Implicit) –
Generalized Blackburn
ecr e Cs
exponential law
εcr C1eC2σ 1 e rt C6 eC7 σ t r C3 s C4 –
C5
Generalized Graham
εcr C1σ C2 t C3 C4t C5 C6t C7 e C8 /T –
Generalized Exponential
εcr C1σ C2 re rt r C5σ C3 e C4 /T • •
Time-hardening creep 법칙에서 파생된 함수들 포함 Generalized Blackburn은 Norton’s law 대신 exponential 법칙을 사용하여 Generalized Exponential처럼 time-hardening 항에 대해 exponential 형태를 포함
크리프(CREEP) 이론 모델 ANSYS WB에서 사용되는 크리프 이론 모델(Implicit) –
Combined Time Hardening
C1σ C2 t C3 1e C4 /T εcr C5σ C6 teC7 /T C3 1 –
Rational Polynomial
εc t cpt εc εmt 1 pt
εcr C1
c C7 εmC8 σ C9 • •
•
εm C210C3σ σ C4 p C10εmC11 σ C12
1차 및 2차 creep 효과를 직접적으로 모사 가능 만약 Time Hardening (TBOPT=11)에서 시간미분을 하면, time-hardening (TBOPT=2)과Norton’s law (TBOPT=10)를 모두 포함 가능 Rational Polynomial 형식은 보통 원자력 분야의 강철에 대해 사용
크리프(CREEP) 이론 모델 ANSYS WB에서 사용되는 크리프 이론 모델(Implicit) –
Generalized Garofalo
εcr C1 sinh C2 σ C3 e C4 /T –
Exponential Form
εcr C1e σ/C 2 e C3 /T –
Norton
εcr C1σ C2 e C3 /T •
creep strain rate의 시간이나 변형률 의존성을 포함하지 않기 때문에, 2차 creep에 알맞음(즉, creep strain rate가 일정)
크리프(CREEP) Curve Fitting 크리프(CREEP) Curve Fitting 해석 MAPDL에서는 실제 재료의 시험데이터를 기반으로 13개의 크리프 이론 모델에 사용되는 상수를 계산함 (구해진 상수 값을 이용하여 해석 수행 가능)
계산 데이터 및 시험데이터 비교
CREEP 모델 상수 확인
※ 사용할 크리프 이론 모델에 따라서 필요한 시험데이터의 종류가 결정됨 ※ 반대로, 시험데이터와 크리프 거동에 맞추어 크리프 이론 모델을 선택 해야함 (1차 크리프인지 2차 크리프인지 선택)
APDL에서 Curve Fitting하기 크리프(CREEP) Curve Fitting MAPDL에서 크리프 Curve Fitting 해석을 위해 Text file로 작성된 시험데이터를 준비해야 함 – – – –
파일은 space나 tab을 통해 열을 구분해야 함 헤더(파일의 처음 몇 줄)는 사용되는 변수에 대한 정보를 포함해야함 사용 가능한 각각의 데이터 열에서, 해더에 /#이나 var를 입력하여 어떠한 변수들이 열과 연관되어 있는지를 지정함 만약, 변수가 상수라고 가정하면 헤더에 /var, value 를 사용함 Abbreviation used in text data file
/seqv,200 /temp,100
Time
time
Equivalent Creep Strain
creq
Equivalent Creep Strain Rate
dcreq
Equivalent Stress
seqv
0.00406109 0.000181314
Temperature
temp
0.00664691 0.000165303
/1,creq /2,dcreq 0.00215869 0.000203055
0.0102068 0.000152217
※ 예를 들어, creep 테스트에서 등가응력이 200이고 온도가 100으로 고정된 상태에서의 결과가 있다면 오 른쪽 테이블과 같이 입력됨
0.0151416 0.000140946 0.0220102 0.000130945
ANSYS WB에서 크리프 재료 활용하기 해석 목적 간단한 예를 통한 크리프 재료 정의 방법 확인 및 거동 특성 확인
• 해석 모델 정의 및 경계조건 설정 –
해석 모델 길이 1m의 정육면체 모델 사용
–
재료 물성 크리프 및 소성 모델 정의
–
경계 조건 3면의 Normal 방향 구속 단일 방향으로 일정한 하중 적용
–
해석 결과 시간에 따른 재료의 거동 확인 각각의 탄성, 소성, 크리프 결과 확인
ANSYS WB에서 크리프 재료 활용하기 크리프 및 소성 재료 설정 Bilinear Isotropic Hardening 재료와 Strain Hardening의 크리프 재료를 함께 정의 BISO 소성 재료 정의
크리프 재료 정의
–
정의한 상수의 단위에 맞춰서 WB에서도 해석을 수행해야 함 (단위 주의!)
ANSYS WB에서 크리프 재료 활용하기 해석 물성 적용 및 해석 설정 –
–
물성 정의
–
해석 설정
–
경계 조건 설정
해석 설정에서 Auto Time Stepping 설정을 활성화하고 Substep 수를 충분히 설정하는 것이 좋음
ANSYS WB에서 크리프 재료 활용하기 해석 결과 간단한 예를 통한 크리프 재료 정의 방법 확인 및 거동 특성 확인
0.43645 m/m 0.00150 m/m 0.13465 m/m 0.30030 m/m Strain (m/m)
Time (sec)
[𝑲] = ( 𝒆𝒍𝒂𝒔𝒕𝒊𝒄 + 𝒑𝒍𝒂𝒔𝒕𝒊𝒄 + 𝒗𝒊𝒔𝒄𝒐 ) Total Strain = Elastic + Plastic + Creep Strain
요약 • 점탄성 점탄성은 고체의 탄성, 유체의 점성 특성을 함께 가진 재료로 항복 변형률 이하에서 재료 정의 - Spring과 Dashpot로 나타낼 수 있음 - 유변학에서 선형 점탄성 영역(LVR)에 포함
• 점소성 점소성은 고체의 탄소성, 유체의 점성 특성을 함께 가진 재료로 항복 변형률 이상에서 재료 정의 - Spring과 Dashpot, Slider로 나타낼 수 있음 - 유변학에서 비선형 점탄성 영역(Non-LVR)에 포함 - 변형률 속도(strain rate) 함수로 시간, 변형률, 응력, 온도에 대한 변수로 함수 정의
• 크리프 -
크리프 이론 모델의 정의 점소성(크리프) 및 소성 모델을 조합하거나 점탄성 모델과 초탄성 재료 조합가능 Curve Fitting Creep 예제 모델
※ 각 재료에 대한 정의 및 특성에 대해 확인하였으며, ANSYS 정의 방법에 대해 확인 하였음
감사합니다. [email protected]
※ 본 자료의 모든 콘텐츠의 저작권은 소프트웨어 개발사와 ㈜태성에스엔이에 있으므로 무단 전재 및 변형, 배포할 수 없습니다.
29