Python3之彈性力學 應力張量1

2022-03-02 12:33:39 字數 2612 閱讀 3290

已知某點的應力張量為:

\[\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...