已知某點的應力張量為:
\[\left[
\begin
\sigma_ &\tau_ &\tau_\\
\tau_ &\sigma_ &\tau_\\
\tau_ &\tau_ &\sigma_
\end
\right]
=\left[
\begin
0 &1 &2\\
1 & \sigma_ & 1\\
2 &1 &0
\end
\right]
\]並已知經過該點的某一平面上的應力向量為零向量,求 \(\sigma_y\) 和主應力?
由題意,存在某個微分面(單位法向量為 \(\boldsymbol\)),其上的應力向量 \(\boldsymbol=\boldsymbol\),即
\[\boldsymbol=\boldsymbol\cdot\boldsymbol=
\left[\begin
0 & 1 & 2\\
1 & \sigma_ & 1\\
2 & 1 & 0
\end
\right]
\left[
\begin
n_1\\
n_2\\
n_3\end
\right]
=\left[
\begin
0\\0\\
0\end
\right]
\]線性方程組存在非零解,必然行列式為零,即
\[\left|\begin
0 & 1 & 2\\
1 & \sigma_ & 1\\
2 & 1 & 0
\end
\right|
= 0 + 2 + 2 -4\sigma_y - 0 - 0
= 0\]
求得 \(\sigma_y = 1\)。
於是,應力張量為
\[\left[
\begin
\sigma_ & \tau_ & \tau_\\
\tau_ & \sigma_ & \tau_\\
\tau_ & \tau_ & \sigma_
\end
\right]
=\left[
\begin
0 & 1 & 2\\
1 & 1 & 1\\
2 & 1 & 0
\end
\right]
\]求主應力,即為求應力張量的特徵值。
\[\left|\,\boldsymbol-\sigma\boldsymbol \,\right| = 0\]或
\[\left|
\begin
-\sigma & 1 & 2\\
1 & 1-\sigma & 1\\
2 & 1 & -\sigma
\end
\right|
=(1-\sigma)\sigma^2 + 2 + 2 - 4(1-\sigma) + \sigma + \sigma = 0
\]整理得
\[-\sigma^3 + \sigma^2 + 6\sigma = -\sigma(\sigma-3)(\sigma+2) = 0
\]主應力
得到三個主應力分別為
\[\left\
\sigma_1 & = & 3\\
\sigma_2 & = & 0\\
\sigma_3 & = & -2
\end
\right.
\]
from sympy import init_printing, matrix
init_printing(use_unicode=true)
# 生成矩陣物件
sigma = matrix([[0, 1, 2], [1, 1, 1], [2, 1, 0]])
sigma
\[\left[\begin0 & 1 & 2\\1 & 1 & 1\\2 & 1 & 0\end\right]
\]\[\left\
\sigma_1 & = & 3\\
\sigma_2 & = & 0\\
\sigma_3 & = & -2
\end
\right.
\]
sigma.eigenvals() # 求特徵值
\[\left \
\]
sigma.eigenvects()
\[\left [ \left ( -2, \quad 1, \quad \left [ \left[\begin-1\\0\\1\end\right]\right ]\right ), \quad \left ( 0, \quad 1, \quad \left [ \left[\begin1\\-2\\1\end\right]\right ]\right ), \quad \left ( 3, \quad 1, \quad \left [ \left[\begin1\\1\\1\end\right]\right ]\right )\right ]
\]
Python3之彈性力學 應力張量2
已知某應力張量的分量為 sigma 3,quad sigma sigma 1,quad sigma sigma 0,quad sigma 2 求 1 全部主應力 2 最大主應力對應的主方向 3 求方向向量為 boldsymbol left 0,dfrac dfrac right 的斜面上的正應力 s...
python3語言 Python3學習之語言基礎3
三 判斷 迴圈語句,函式,命名空間,作用域 1 python3 條件控制 python中if語句的一般形式如下所示 python 中用 elif 代替了 else if,所以if語句的關鍵字為 if elif else。注意 1 每個條件後面要使用冒號 表示接下來是滿足條件後要執行的語句塊。2 使用...
python3語言 Python3學習之語言基礎1
一 python3入門,資料型別,字串 python 中的變數不需要宣告。每個變數在使用前都必須賦值,變數賦值以後該變數才會被建立。usr bin python3 counter 10 整型變數 miles 3.14 浮點型變數 name bob 字串 print counter print mil...