응력 변환 (Stress Transformation)

변형률처럼 응력 텐서는 정의와 상관없이 좌표변환 또는 회전행렬의 전/후 곱을 통하여 변환(transform)할 수 있다.

2-D 응력 좌표 변환 (coordinate transforms)

좌표변환은 물체는 그대로 있고 좌표축이 회전한 경우이다. 주로 응력의 주방향(principal directions)을 찾는데 활용된다.

x, y 방향 힘의 평형은 다음과 같다.

σxxAcosθ+τxyAsinθ=σxxAcosθτxyAsinθτxyAcosθ+σyyAsinθ=σxxAsinθ+τxyAcosθ

면적 A를 소거하고 미지수 σxx,τxy 및 θ+90˚ 면수직응력 σyy를 구하면

σxx=σxxcos2θ+σyysin2θ+2τxysinθcosθσyy=σxxsin2θ+σyycos2θ2τxysinθcosθτxy=(σyyσxx)sinθcosθ+τxy(cos2θsin2θ)

행렬 표기법으로 쓰면

[σxxτxyτyxσyy]=[  cosθsinθsinθcosθ][σxxτxyτyxσyy][cosθsinθsinθ  cosθ]

앞 행렬을 Q로 쓰면 뒷 행렬은 그의 전치행렬이므로 다음과 같다.

σ=Q σ QT

[예제] 응력텐서 σ가 기준좌표계에서 다음과 같을 때 좌표계가 50˚ 회전했을 때 응력 텐서 σ'는

σ=[1223]=[  cos50sin50sin50cos50][1223][cos50sin50sin50  cos50]=[4.143   0.6380.6380.143]

좌표변환을 텐서표기법(tensor notation)으로 쓰면

σmn=QmiQnjσij

3-D 응력 좌표 변환 (coordinate transforms)

좌표축이 추가되었을뿐 원칙은 2-D와 동일하다.

σ=Q σ QTσmn=QmiQnjσijQij=cos(xi,xj)

행렬을 풀어 쓰면

[σxxτxyτxzτyxσyyτyzτzxτzyσzz]=[Q11Q12Q13Q21Q22Q23Q31Q32Q33][σxxτxyτxzτyzσyyτyzτzxτzyτzz][Q11Q21Q31Q12Q22Q32Q13Q23Q33]

회전 변환 (Rotations)

회전 변환은 좌표축을 유지하면서 물체를 회전한 경우이다. 회전 행렬 R극좌표 분해(polar decomposition)로부터 계산된다.

응력텐서의 회전은

σ=R σ RT

2-D에서 회전행렬은 변환행렬의 전치(R=QT)이므로 각도 θ 만큼 강체 회전의 경우

[σxxτxyτyxσyy]=[cosθsinθsinθ  cosθ][σxxτxyτyxσyy][  cosθsinθsinθcosθ]

변환 후 각 응력성분을 풀어 쓰면

σxx=σxxcos2θ+σyysin2θ2τxysinθcosθσyy=σxxsin2θ+σyycos2θ+2τxysinθcosθτxy=(σxxσyy)sinθcosθ+τxy(cos2θsin2θ)

3-D에서는

[σxxτxyτxzτyxσyyτyzτzxτzyσzz]=[R11R12R13R21R22R23R31R32R33][σxxτxyτxzτyzσyyτyzτzxτzyτzz][R11R21R31R12R22R32R13R23R33]

[예제] 위의 좌표계 변환 예제를 이번에는 각도 50˚의 강체회전에 적용한다.

σ=[cos50sin50sin50  cos50][1223][  cos50sin50sin50cos50]=[  0.2041.3321.332  3.796]

좌표 및 회전 역변환 (Transforming and rotating back)

좌표 변환식 양변의 앞에 QT, 뒤에 Q를 곱하면 QT=Q1 이므로

QT σ Q=QT Q σ QT Q=σ 

같은 방법으로 물체가 회전 후 원래 각도로 복귀하는 변환은

σ=RT σ R

출처 http://www.continuummechanics.org

댓글

이 블로그의 인기 게시물

전단응력 (Shear Stress)

절대압력과 계기압력

엑셀 상자그림(Box Plot) 그리기