位移与应变

贡献者： ACertainUser

本文存在未完成的内容。

　　¹ 在应力中我们讨论了如何表示材料中各处的受力，现在我们讨论如何表示材料中各处的变形。此处我们只探讨小变形的情况，即材料变形的程度非常轻微。

1. 位移

图 1：材料变形后，材料中的每一点的位置发生改变

图 2：上图的叠加

　　材料变形后，材料中每一点都移动了。例如，原先位于 $x_{0}$ 的点在变形后运动到了 $a_{0}$ 位置，原先位于 $x_{1}$ 的点运动到了 $a_{1}$ 位置...如果我们在弹性力学中暂且假定 “变形前后，原先的一个点不会分裂为两个，也不会有两个点合并为一个”，那么我们就可以在点变形前的位置 $x$ 与变形后的位置 $a$ 之间建立一种映射 $x \to a,$ 亦即定义一个函数。 $a = a (x) .$

位移函数

　　为了更好地表示 $x$ 处点位置的变化，我们定义位移函数

\begin{matrix} (1) & u (x) = a (x) - x, \end{matrix}

由此，点

x

的位置变化也可以表述为

\begin{matrix} (2) & x + u (x) = a (x) . \end{matrix}

　　 $u$ 是一个矢量函数，也被称为位移场。它的 3 个分量都是关于点坐标的函数，分别表示点在 $x_{1}, x_{2}, x_{3}$ ²各方向上的位移。 $u (x) = (\begin{matrix} u_{1} (x) \\ u_{2} (x) \\ u_{3} (x) \end{matrix}) = (\begin{matrix} u_{1} (x_{1}, x_{2}, x_{3}) \\ u_{2} (x_{1}, x_{2}, x_{3}) \\ u_{3} (x_{1}, x_{2}, x_{3}) \end{matrix}) .$ 由于变形前后点不合并、分裂的性质， $u (x)$ 是一个 “好函数”，单值、连续且多阶可导。

2. 应变

　　有时我们使用应变 $ε$ 来描述材料的变形。类似于应力，应变也定义在材料中的每一个微元处，可记为一个 $3 \times 3$ 的矩阵（二阶张量）.应变矩阵也是对称矩阵 $ε = ε^{T}$ ，有 $6$ 个独立变量。

\begin{matrix} (3) & ε = (\begin{matrix} ε_{11} & ε_{12} & ε_{13} \\ ε_{21} & ε_{22} & ε_{23} \\ ε_{31} & ε_{32} & ε_{33} \end{matrix}) . \end{matrix}

　　应变角标的含义与应力的类似，第一个角标表示面的法方向，第二个角标表示变形的方向。我们可将应变分为两类，正应变（两角标相同）与切应变（两角标不同）。

正应变

图 3：正应变示意图。改绘自冯西桥的《弹性力学》

　　正应变与微元的长度变化有关。如图 3 所示， $x_{1}$ 方向上，微元的原本长度为 $l_{0} = d x_{1}$ ，变形后长度为 $l = d x_{1} + \frac{\partial u_{1}}{\partial x_{1}} d x_{1}$ ，长度变化 $Δ l = l - l_{0} = \frac{\partial u_{1}}{\partial x_{1}} d x_{1}$ ，那么定义正应变

\begin{matrix} (4) & ε_{11} = \frac{Δ l}{l_{0}} = \frac{\partial u_{1}}{\partial x_{1}} . \end{matrix}

切应变

图 4：切应变示意图。改绘自冯西桥的《弹性力学》

　　而切应变与微元的角度变化有关。如图 4 所示，定义切应变

\begin{matrix} (5) & ε_{12} = \frac{1}{2} (90^{\circ} - α) = \frac{1}{2} (\frac{\partial u_{1}}{\partial x_{2}} + \frac{\partial u_{2}}{\partial x_{1}}), \end{matrix}

有时也使用工程切应变

γ

。别紧张，他们只是相差一个系数

\begin{matrix} (6) & γ_{12} = 2 ε_{12} = 90^{\circ} - α . \end{matrix}

3. 应变几何方程

　　推广上一节的结论，我们共可以得到 $6$ 个联系应变与位移的独立方程。这 $6$ 个方程被称为应变几何方程。

\begin{matrix} (7) & ε_{i j} = \frac{1}{2} (\frac{\partial u_{i}}{\partial x_{j}} + \frac{\partial u_{j}}{\partial x_{i}}) (i, j = 1, 2, 3) . \end{matrix}

4. 应变协调方程

5. 画廊：经典的应变类型

　　图中黑色为原微元体，红色为变形后的微元体。

图 5：

u_{1} = 1.5

。各应变均为

0

，微元形状不变、只发生平移。

图 6：

u_{1} = 0.2 x_{1}

。

ε_{11} = 0.2

（其余项为

0

），微元发生拉伸变形

图 7：

u_{1} = 0.2 x_{2}

ε_{12} = 0.1

,微元发生简单的剪切变形

图 8：

u_{1} = 0.2 x_{2}, u_{2} = 0.2 x_{1}

ε_{12} = 0.2

，微元发生剪切变形

图 9：同上，但是是俯视图

6. 附录：应变模拟器

　　用以绘制上述应变示意图的 Matlab 代码

%定义位移函数，你可以设置自己的线性位移函数
u1 = @(x,y,z) -0.2*x + 0.2*y;
u2 = @(x,y,z) 0;
u3 = @(x,y,z) 0;

A(1,:)=[0,0,0];
A(2,:)=[1,0,0];
A(3,:)=[1,1,0];
A(4,:)=[0,1,0];
A(5,:)=[0,1,1];
A(6,:)=[0,0,1];
A(7,:)=[1,0,1];
A(8,:)=[1,1,1];

hold on
axis equal

%绘制顶点
%for i = 1:8
%    scatter3(A(i,1),A(i,2),A(i,3),'k');
%endfor

view(-30,60)
xlabel('Axis 1','fontsize',15)
ylabel('Axis 2','fontsize',15)
zlabel('Axis 3','fontsize',15)

line([A(1,1), A(2,1)],[A(1,2), A(2,2)],[A(1,3), A(2,3)],'color','k');
line([A(1,1), A(4,1)],[A(1,2), A(4,2)],[A(1,3), A(4,3)],'color','k');
line([A(1,1), A(6,1)],[A(1,2), A(6,2)],[A(1,3), A(6,3)],'color','k');

line([A(2,1), A(7,1)],[A(2,2), A(7,2)],[A(2,3), A(7,3)],'color','k');
line([A(2,1), A(3,1)],[A(2,2), A(3,2)],[A(2,3), A(3,3)],'color','k');

line([A(3,1), A(4,1)],[A(3,2), A(4,2)],[A(3,3), A(4,3)],'color','k');
line([A(3,1), A(8,1)],[A(3,2), A(8,2)],[A(3,3), A(8,3)],'color','k');

line([A(4,1), A(5,1)],[A(4,2), A(5,2)],[A(4,3), A(5,3)],'color','k');

line([A(6,1), A(5,1)],[A(6,2), A(5,2)],[A(6,3), A(5,3)],'color','k');
line([A(6,1), A(7,1)],[A(6,2), A(7,2)],[A(6,3), A(7,3)],'color','k');

line([A(8,1), A(5,1)],[A(8,2), A(5,2)],[A(8,3), A(5,3)],'color','k');
line([A(8,1), A(7,1)],[A(8,2), A(7,2)],[A(8,3), A(7,3)],'color','k');

for i = 1:8
    B(i,1) = u1(A(i,1),A(i,2),A(i,3))+A(i,1);
    B(i,2) = u2(A(i,1),A(i,2),A(i,3))+A(i,2);
    B(i,3) = u3(A(i,1),A(i,2),A(i,3))+A(i,3);
endfor

%for i = 1:8
%    scatter3(B(i,1),B(i,2),B(i,3),'r');
%endfor

line([B(1,1), B(2,1)],[B(1,2), B(2,2)],[B(1,3), B(2,3)],'color','r');
line([B(1,1), B(4,1)],[B(1,2), B(4,2)],[B(1,3), B(4,3)],'color','r');
line([B(1,1), B(6,1)],[B(1,2), B(6,2)],[B(1,3), B(6,3)],'color','r');

line([B(2,1), B(7,1)],[B(2,2), B(7,2)],[B(2,3), B(7,3)],'color','r');
line([B(2,1), B(3,1)],[B(2,2), B(3,2)],[B(2,3), B(3,3)],'color','r');

line([B(3,1), B(4,1)],[B(3,2), B(4,2)],[B(3,3), B(4,3)],'color','r');
line([B(3,1), B(8,1)],[B(3,2), B(8,2)],[B(3,3), B(8,3)],'color','r');

line([B(4,1), B(5,1)],[B(4,2), B(5,2)],[B(4,3), B(5,3)],'color','r');

line([B(6,1), B(5,1)],[B(6,2), B(5,2)],[B(6,3), B(5,3)],'color','r');
line([B(6,1), B(7,1)],[B(6,2), B(7,2)],[B(6,3), B(7,3)],'color','r');

line([B(8,1), B(5,1)],[B(8,2), B(5,2)],[B(8,3), B(5,3)],'color','r');
line([B(8,1), B(7,1)],[B(8,2), B(7,2)],[B(8,3), B(7,3)],'color','r');

图 10：简单的剪切变形

　　另一个以点阵代替正方体的实现的 Matlab 代码：

clc
clear

function s = SHPERE(x,y,z,r)
  [X,Y,Z] = sphere(50);
  X = r*X+x;
  Y = r*Y+y;
  Z = r*Z+z;
  s = surf(X,Y,Z);
end

u1 = @(x,y,z) 0.2*y;
u2 = @(x,y,z) -0.05*y;
u3 = @(x,y,z) 0;

hold on
axis equal

for x = 0:3
  for y = 0:3
    for z = 0:3
      a = u1(x,y,z)+x;
      b = u2(x,y,z)+y;
      c = u3(x,y,z)+z;
      SPHERE(a,b,c,0.5);
    end
  end
end

1. ^ 本文参考了冯西桥的《弹性力学》课程与陆明万的《弹性理论基础》
2. ^ 为了方便表示，以 $x_{1}, x_{2}, x_{3}$ 轴代指 $x, y, z$ 轴

致读者：小时百科一直以来坚持所有内容免费无广告，这导致我们处于严重的亏损状态。长此以往很可能会最终导致我们不得不选择大量广告以及内容付费等。因此，我们请求广大读者热心打赏 ，使网站得以健康发展。如果看到这条信息的每位读者能慷慨打赏 20 元，我们一周就能脱离亏损，并在接下来的一年里向所有读者继续免费提供优质内容。但遗憾的是只有不到 1% 的读者愿意捐款，他们的付出帮助了 99% 的读者免费获取知识，我们在此表示感谢。