1. 引论

有限的单元，无限的能力

基本就是说说有什么用，不过需要吐槽一下，那些视频看上去很炫酷，但是在没学过FEM的人看来就平平无奇，也不知道干什么用的。

引论

1.1 力学的分类：质点、刚体、变形体的力学

以卫星为例子，描述轨道只要简化为质点就可以，使用三个坐标就可以标示，如果描述姿态还需要再加3个坐标，如果考虑电池板等柔性元件就更加复杂。

想起一个大佬说的，做有限元分析，不是考虑加什么东西，而是考虑在不影响结果的情况下，可以从现实世界的模型里减去什么。

下面更加具体的分析：

质点力学：初中学过的
刚体力学（理论力学）：物体形状不变，但是质量分布不能简化为一个点
变形体的力学
- 简单形状
  - 简单形状，简单组合：材料力学，小变形，形状简单
  - 简单形状，复杂组合：结构力学，各种复杂桁架的受理分析（土木应该会经常用吧）
  - 复杂形状：弹塑性力学，需要使用微元法分析了，预感这是这门课解决的问题，未必能用解析法解析了。

1.2 变形体力学的要点

变形体力学的要点，曾老师设计了三个同种材料单拉实验，三个实验的受力-位移曲线都不一样，引导出描述材料中的不变量：弹性模量E的定义：

$E = \frac{\sigma}{\epsilon}$

其中：

$\sigma = F/A$
$\epsilon=u/L$
F 为试样受拉力
A 为试样横截面积
u 为受拉力F以后，试样的形变量
L 为试样长度

1.3 微分方程求解的方法

这是我踩的第一个坑，我想了10分钟才明白为什么是u的二阶导数。这里并不是位移-速度-加速度。首先视频里说“长度上受到均布载荷P”，我里所应当的认为P是压强，应变也就是F/A，其实不是。

我们如果在自由端施加F的拉力，的确小杆子内部的应变处处相等，但是这样的简单载荷是不需要微分方程求解的，所以给定了一个复杂的情况，这里的P的单位其实是N/m，杆子内部的受力不均匀，是一个从左到右逐渐加码的载荷。

所以，应变的大小其实是 $\epsilon(x) = \frac{du}{dx}$ ，那么应力就是 $\sigma(x) = E \epsilon(x)$ ，截面受力 $F(x) = EA \epsilon(x)$ 考虑到P是均匀加负荷，所以F的导数等于P，也就是力从0开始逐渐变大到LP（从右到左看的话）。