[物理学与PDEs]第5章第6节 弹性静力学方程组的定解问题
5. 6 弹性静力学方程组的定解问题
5. 6. 1 线性弹性静力学方程组
1. 线性弹性静力学方程组 $$\bee\label{5_6_1_le} -\sum_{j,k,l}a_{ijkl}\cfrac{\p ^2u_k}{\p x_j\p x_l}=\rho_0b_i,\quad i=1,2,3. \eee$$
2. (Korn 不等式) 设 $\Omega\subset{\bf R}^3$ 为有界区域, 则 $$\bex \exists\ C_0>0,\st \int_\Omega |{\bf E}|^2\rd x \geq C_0\sen{{\bf u}}_{H^1(\Omega)}^2,\quad \forall\ {\bf u}\in H^1_0(\Omega). \eex$$
3. 若 ${\bf A}=(a_{ijkl})$ 满足稳定性条件, 则 \eqref{5_6_1_le} 在有界区域上有弱解的存在性、唯一性及正则性.
5. 6. 2 非线性弹性静力学方程组
1. 非线性弹性静力学方程组 $$\bee\label{5_6_2_nle} -\sum_{j,k,l}a_{ijkl}(\n{\bf y})\cfrac{\p ^2y_k}{\p x_j\p x_l}=\rho_0b_i,\quad i=1,2,3. \eee$$
2. $I=\sev{\ba{ccc} n_1&n_2&n_3\\ \cfrac{\p y_2}{\p x_1}&\cfrac{\p y_2}{\p x_2}&\cfrac{\p y_2}{\p x_3}\\ \cfrac{\p y_3}{\p x_1}&\cfrac{\p y_3}{\p x_2}&\cfrac{\p y_3}{\p x_3} \ea}$ 仅与 $y_2,y_3$ 的切向偏导数有关. 事实上, $$\beex \bea I&=-\cfrac{\p y_2}{\p x_1}\sex{n_2\cfrac{\p y_3}{\p x_3}-n_3\cfrac{\p y_3}{\p x_2}}\\ &\quad+\cfrac{\p y_2}{\p x_2}\sex{n_1\cfrac{\p y_3}{\p x_3} -n_3\cfrac{\p y_3}{\p x_1}}\\ &\quad-\cfrac{\p y_2}{\p x_3}\sex{n_1\cfrac{\p y_3}{\p x_2}-n_2\cfrac{\p y_3}{\p x_1}}. \eea \eeex$$ 记 $$\beex \bea {\bf a}&=(0,-n_3,n_2)^T,\\ {\bf b}&=(-n_3,0,n_1)^T,\\ {\bf c}&=(-n_2,n_1,0)^T,\\ {\bf d}&=\sex{ -n_2\cfrac{\p y_3}{\p x_3}+n_3\cfrac{\p y_3}{\p x_2}, n_1\cfrac{\p y_3}{\p x_3} -n_3\cfrac{\p y_3}{\p x_1},-n_1\cfrac{\p y_3}{\p x_2}+n_2\cfrac{\p y_3}{\p x_1}}, \eea \eeex$$ 则直接计算有 $$\bex {\bf a}\cdot{\bf n}={\bf b}\cdot{\bf n}={\bf b}\cdot{\bf n}={\bf d}\cdot{\bf n}=0. \eex$$ 于是 $$\bex I=\n_{{\bf d}}y_2,\quad {\bf d}=(-\n_{{\bf a}}y_3,-\n_{{\bf b}}y_3,-\n_{{\bf c}}y_3). \eex$$
3. \eqref{5_6_2_nle} 可赋以的边界条件是 ${\bf y}={\bf y}_0$, ${\bf P}{\bf n}={\bf\sigma}({\bf y},\n_{{\bf \tau}}{\bf y})$, 或其组合.
4. \eqref{5_6_2_nle} 的解的唯一性一般不成立.
5. 问题: $$\bee\label{5_6_2_bvp} \bea -\sum_{j,k,l}a_{ijkl}\cfrac{\p ^2u_k}{\p x_j\p x_l}=\rho_0b_i,\quad i=1,2,3,&\quad\mbox{in }\Omega;\\ {\bf y}|_{\vGa_1}={\bf h},\ {\bf P}{\bf n}={\bf \sigma},&\quad\mbox{on }\p\Omega \eea \eee$$ 的解的存在性.
(1) 化 \eqref{5_6_2_bvp} 为泛函 (弹性体总势能) 的极小值问题.
(2) 通过下方有界、极小化序列、弱下半连续等方法求得极小值问题的解.
(3) 弱解的正则性 (问题).
6. 多凸的概念: $\hat W({\bf F})$ 为多凸的 (polyconvex), 如果 $\hat W({\bf F})$ 能表为 ${\bf F}$ 的子行列式的凸函数.
7. 函数 $g({\bf F})=\tr({\bf F}^T{\bf F})$, $\forall\ {\bf F}\in \bbR^{3\times 3}$ 是严格凸的.
[物理学与PDEs]第5章第6节 弹性静力学方程组的定解问题的更多相关文章
- [物理学与PDEs]第5章第5节 弹性动力学方程组及其数学结构
5.5.1 线性弹性动力学方程组 1. 线性弹性动力学方程组 $$\beex \bea 0&=\rho_0\cfrac{\p{\bf v}}{\p t}-\Div_x{\bf P}-\r ...
- [物理学与PDEs]第4章第2节 反应流体力学方程组 2.4 反应流体力学方程组的数学结构
1. 粘性热传导反应流体力学方程组是拟线性对称双曲 - 抛物耦合组. 2. 理想反应流体力学方程组是一阶拟线性对称双曲组 (取 ${\bf u},p,S,Z$ 为未知函数). 3. 右端项具有间 ...
- [物理学与PDEs]第4章第2节 反应流体力学方程组 2.3 混合气体状态方程
1. 记号与假设 (1) 已燃气体的化学能为 $0$. (2) 单位质量的未燃气体的化学能为 $g_0>0$. 2. 对多方气体 (理想气体当 $T$ 不高时可近似认为), $$\bex ...
- [物理学与PDEs]第4章第2节 反应流体力学方程组 2.2 反应流体力学方程组形式的化约
1. 粘性热传导反应流体力学方程组 $$\beex \bea \cfrac{\rd \rho}{\rd t}&+\rho \Div{\bf u}=0,\\ \cfrac{\rd Z}{\rd ...
- [物理学与PDEs]第4章第2节 反应流体力学方程组 2.1 粘性热传导反应流体力学方程组
1. 记号: $Z=Z(t,{\bf x})$ 表示未燃气体在微团中所占的百分比 ($Z=1$ 表示完全未燃烧; $Z=0$ 表示完全燃烧). 2. 物理化学 (1) 燃烧过程中, 通过化学反应 ...
- [物理学与PDEs]第3章第4节 磁流体力学方程组的数学结构
1. 在流体存在粘性.热传导及 $\sigma\neq \infty$ 时, 磁流体力学方程组是一个拟线性对称双曲 - 抛物耦合组. 2. 在流体存在粘性.热传导但 $\sigma=\infty$ ...
- [物理学与PDEs]第3章第2节 磁流体力学方程组 2.4 不可压情形的磁流体力学方程组
不可压情形的磁流体力学方程组 $$\beex \bea \cfrac{\rd {\bf H}}{\rd t}-({\bf H}\cdot\n){\bf u}&=\cfrac{1}{\sigma ...
- [物理学与PDEs]第3章第2节 磁流体力学方程组 2.3 磁流体力学方程组
1. 磁流体力学方程组 $$\beex \bea \cfrac{\p {\bf H}}{\p t} &-\rot({\bf u}\times{\bf H})=\cfrac{1}{\sigma ...
- [物理学与PDEs]第3章第2节 磁流体力学方程组 2.2 考虑到电磁场的存在对流体力学方程组的修正
1. 连续性方程 $$\bex \cfrac{\p \rho}{\p t}+\Div(\rho{\bf u})=0. \eex$$ 2. 动量守恒方程 $$\bex \cfrac{\p }{\p ...
随机推荐
- GIL:全局解释器锁 VS 用户程序锁
既然有了GIL锁,CPython还要多线程干什么? ''' GIL:全局解释器锁的来历 四核:同一时刻真正有四个任务在运行,多核的意义在于此 单核:看上去是并发的,因为进行了上下文切换,单核永远是串行 ...
- C#基础知识之泛型
泛型在c#中有很重要的位置,对于写出高可读性,高性能的代码有着关键的作用. 其实官方文档说明的很详细,我这边算是做个记录吧 一.什么是泛型? 泛型是 2.0 版 C# 语言和公共语言运行库 (CLR) ...
- 转://Oracle数据库升级后保障SQL性能退化浅谈
一.数据库升级后保障手段 为了保障从10.2.0.4版本升级到11.2.0.4版本更加平稳,我们事先采用了oracle性能分析器(SQL Performance Analyzer)来预测数据库的关键S ...
- php7 的yum源
yum源默认的版本太低了,手动安装有一些麻烦,想采用Yum更新安装的可以使用下面的方案: 1.检查当前安装的PHP包 yum list installed | grep php 如果有安装的PHP包, ...
- SegNet 理解与文章结构
SegNet: A Deep Convolutional Encoder-Decoder Architecture for Image Segmentation 发表于2016年,作者 Vijay B ...
- 使用jprofiler分析dump文件一个实例
3 .dump 线上文件栈 [root@yszyz10a153 ~]# jmap -dump:live,format=b,file=heap201712.hropf 72947 参考:https:/ ...
- 云端安装MQTT服务器
如果自己下载的3.1版本的MQTT, 安装步骤参考 https://developer.emqx.io/docs/emq/v3/cn/install.html 配置用户名和密码第一种是用http ht ...
- stm32之不定长接收
使用STM32CUBE_MAX配置工程,可以简化编程工作量,但是这样我们会遇到一些麻烦,比如利用串口接收不知道长度的数据的时候,我们可能会无从下手,前段时间看到他人程序中的串口不定长接收,此次特意总结 ...
- 横线和文字一排,文字居中显示vertical-align: middle;
<!DOCTYPE html><html> <head> <meta charset="UTF-8"> <title>& ...
- python+opencv读取视频,调用摄像头
引用 import cv2 import numpy 创建摄像头对象 cap = cv2.VideoCapture("videoTest/test1.mp4") #参数为视频文件目 ...