【连续介质力学】连续体运动学

最新推荐文章于 2024-02-20 16:07:52 发布

LAWKAWAI

最新推荐文章于 2024-02-20 16:07:52 发布

阅读量685

点赞数

分类专栏：连续介质力学文章标签：笔记

本文链接：https://blog.csdn.net/weixin_38362786/article/details/131086254

版权

连续介质力学专栏收录该内容

16 篇文章 8 订阅

订阅专栏

简介

一个质点从t=0出发，随着时间有不同的构形，运动的描述是运动学。

所以需要建立运动方程来表征连续体是如何演化及其性质（例如位移、速度、加速度、质量密度、温度等）如何随时间变化

初始构形或者参考构形： $\mathcal{B}_0$

当前构形： $\mathcal{B}_t$

在这里插入图片描述
从描述单个质点的运动开始

研究在变形当中质点之间的相对距离

定义变形和应变张量，在这之前先定义连续性质

连续介质

任何介质都有质量

如果连续体在任何地方性质一致，被称为齐次的

考虑某个点P为中心的球作为初始构形：
在这里插入图片描述

体积： $\Delta V_0$
质量： $\Delta m$
那么，质量密度：
在这里插入图片描述
质点：一个微小体积元，所具有的确定的性质，例如质量密度、速度、温度等

点：空间中的一个点

质点迹线：单个质点运动过程中的轨迹
在这里插入图片描述

运动的类型

刚体运动：保持原来的构形，表征为质点之间的距离保持一致；可以分类为：位移和旋转

带变形的运动：表征为质点之间的距离会发生变化

刚体运动

建立运动方程，直角坐标系 $OX_1X_2X_3$ 在刚体上，所以任意质点的位置向量在运动过程中是不会发生改变的

另一个直角坐标系 $ox_1x_2x_3$ 表示正交基 $(\hat e_1, \hat e_2, \hat e_3)$

在这里插入图片描述
质点Ｐ在坐标系 $\hat e_i$ 和 $\hat I_i$ 的位置向量为 $\vec X$ 和 $\vec x$
那么：
$\vec x = \vec c + \vec X$
其中 $\vec c(t)$ 与时间相关，描述的是坐标系 $\hat I_i$ 的位移运动
下标表示：

在 $ox_1x_2x_3$ 的分量可以点乘 $\hat e_i$ 得到：
在这里插入图片描述
其中 $\hat I_j \cdot \hat e_i = a_{ji}$ 表示从坐标系 $\hat I_i$ 到坐标系 $\hat e_i$ 的变换矩阵

并且 $a_{ik} a_{kj} = \delta_{ij}$ ， $a_{ji}$ 是一个正交矩阵

考虑 $Q_{ij} = a_{ji}$ ，那么：
$\boxed{\vec x = \vec c + Q \cdot \vec X} 刚体运动方程$

NOTE：分量的变换定律和正交变换密切相关，但它们有着完全不同的含义

问题2.1 一个连续体是以b为边的正方形，受到一个刚体运动，这个刚体运动是以逆时针旋转30°的旋转运动，求出运动方程，并求出质点D的新的位置

在这里插入图片描述

$\vec c = \vec 0$ ，所以空间坐标系和质点坐标系重叠了
在这里插入图片描述

构形的类型

参考构形或者初始构形： 在时间 $\equiv t_0$ ，即未变形之前的构形，质点P的位置向量为 $\vec X^P$

当前构形或者变形构形： 在时间 $t = t$ 的构形

运动可以表征为一个双射函数 $\varphi$ ，这可以保证存在逆函数 $\varphi ^{-1}$

在这里插入图片描述

质量密度

当前构形的质量密度：
在这里插入图片描述
质量密度是一个标量场，是一个关于位置和时间的函数： $\rho=\rho(\vec x, t)$

运动的描述

质点和空间坐标系

质量体 $\mathcal B_0$ 在初始构形

在任意时间，质量体会占据一个新的空间 $\mathcal B_t$ ，考察质点P

在这里插入图片描述
质点P在参考构形，时间t=0的位置向量为：

定义质点坐标系：
在这里插入图片描述
质点P在当前构形位于位置 P’, 位置向量为：

其中，定义了空间坐标系：
在这里插入图片描述

位移向量

位移向量定义为当前构形的位置向量 $\vec x$ 和参考构形的位置向量 $\vec X$ 的差：
在这里插入图片描述

速度向量

定义为位置向量的变化率：
在这里插入图片描述

加速度向量

定义为速度向量的变化率：
在这里插入图片描述

拉格朗日描述和欧拉描述

连续的性质：质量密度、温度、速度、加速度是质点固有的，这些性质可以随着时间变化，连续运动可以表征为一个双射的函数 $\varphi$ 和逆函数 $\varphi^{-1}$

这保证了我们可以联系当前构形和参考构形之间的连续性质

运动的拉格朗日描述

用质点坐标系 $\vec X$ 表示的质点运动：
在这里插入图片描述

这个运动方程被称为拉格朗日运动描述

这个方程给出了在当前时刻 t 的位置 $\vec x$ ，其中质点在时刻 $t = t_0$ 占据着位置 $\vec X$

上面的方程如果应用于某个质点P, 则该方程描述的是这个质点的迹线

运动的欧拉描述

用空间坐标系 $\vec x$ 表示的质点运动：
在这里插入图片描述
这个方程给出了在初始时刻 $t_0$ 的初始位置 $\vec X$ ，该质点在当前时刻 $t$ 具有坐标 $x_1, x_2, x_3)$

在这里插入图片描述
在 $t = 0$ 时刻有： $\vec x^P(\vec X, t=0)=\vec X^P$

质点P的迹线：
$\vec x(\vec X^P,t_0)=\vec X^P \rightarrow \vec x(\vec X^P,t_1)=\vec X^{P'} \rightarrow \vec x(\vec X^P,t_2)=\vec X^{P''}$

不同时刻在点P的质点：
$\vec X(\vec x^P,t_0)=\vec X^P \rightarrow \vec X(\vec x^{S'},t_1)=\vec X^{S} \rightarrow \vec X(\vec x^{Q''},t_2)=\vec X^Q$

其中 $\vec X(\vec x^{S'},t_1)=\vec X^{S}$ 表示在时刻 $t_1$ 位于点 $P = S ‘$ 的质点为S，其初始位置为 $\vec X^S$

存在逆的充分必要条件是：雅可比行列式不为0
在这里插入图片描述
**不可穿透性公理：**两个粒子不能同时占据同一位置。如后所述，当雅可比行列式为正时，这个条件是确定的

拉格朗日变量和欧拉变量

在连续体 $\mathcal B$ 上的物理量 $\mathcal Z$ 可以表达成：

拉格朗日形式： $(Z(\vec X, t))$

欧拉形式： $z(\vec x, t)$

在这里插入图片描述

问题2.2 考虑以下拉格朗日描述的运动方程：

在这里插入图片描述

质点时间导数

随时间变化的导数： $\frac{D}{Dt}$

跟着质点P运动并且记录随着时间变化的温度

在这里插入图片描述
如果属性是由拉格朗日描述：
$\theta = \theta(X_1, X_2, X_3, t)$

这种情况下，质点的时间导数表示为：

$\boxed{\dot \theta (\vec X, t)\equiv\frac{D\theta(\vec X, t)}{Dt}=\frac{d\theta (\vec X, t)}{dt}}$

这个属性是由质点坐标系描述的，意味着该属性与同一个质点运动过程相连接

如果属性是由欧拉描述：
$\theta = \theta(x_1, x_2, x_3, t)$

观察不再跟踪质点P，而是固定在某个点 $x_1, x_2, x_3)$ ，观察经过的质点
在这里插入图片描述
在时刻 $t_1$ : 观测到质点Q的属性
在时刻 $t_2$ : 观测到质点R的属性
在时刻 $t_3$ : 观测到质点P的属性

需要强调的是质点的时间导数与质点的固有属性关于时间的导数有关，例如，这是关于同一个质点的

然而，一个固定在某个点的观察者，只能得到当地的变化率的信息
为了得到完整的信息，我们需要知道质点沿着迹线其属性是怎么变化的，而这个额外的信息被称为对流的变化率，这与质量输运有关

所以，为了得到质点的时间导数，必须考虑以下两项：

当地变化率
对流变化率
所以：

质点速度： $\vec v(\vec X, t)=\dot {\vec x}(\vec X, t)$

欧拉描述： $\vec v(\vec X(\vec x, t), t)=\vec v(\vec x, t)$

可以定义欧拉描述下的时间导数算子： $*(\vec x, t)$
$\boxed{\frac{D * (\vec x, t)}{Dt}=\frac{\partial *(\vec x, t)}{\partial t}+\nabla_{\vec x}*(\vec x, t)\cdot \vec v(\vec x, t)}欧拉描述的质点时间导数$