【连续介质力学】涉及积分的定理

2024-03-14 18:10

本文主要是介绍【连续介质力学】涉及积分的定理,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

涉及积分的定理

分部积分

分部积分:
∫ a b u ( x ) v ′ ( x ) d x = u ( x ) v ( x ) ∣ a b − ∫ a b v ( x ) u ′ ( x ) d x \int_a^bu(x)v'(x)dx=u(x)v(x)|_a^b-\int_a^bv(x)u'(x)dx abu(x)v(x)dx=u(x)v(x)ababv(x)u(x)dx

其中, v ′ ( x ) = d v d x v'(x) = \frac{dv}{dx} v(x)=dxdv

散度定理

给定一个体积为V的域B,边界为S,那么应用在向量场的散度定理(也叫Green’s Theorem, 格林定理)为:
∫ V ∇ x ⃗ ⋅ v ⃗ d V = ∫ S v ⃗ ⋅ n ^ d S = ∫ S v ⃗ ⋅ d S ⃗ ∫ V v i , i d V = ∫ S v i n ^ i d S = ∫ S v i d S i \int_V \nabla_{\vec x}\cdot \vec vd V = \int_S \vec v\cdot \hat ndS = \int _S \vec v\cdot d\vec S \\ \int_V v_{i,i}dV=\int_Sv_i\hat n_idS=\int_Sv_idS_i Vx v dV=Sv n^dS=Sv dS Vvi,idV=Svin^idS=SvidSi
其中, n ^ \hat n n^ 是向外垂直于表面S的
在这里插入图片描述
令T是在域B的二阶张量场,那么应用散度定理在这个场:
∫ V ∇ x ⃗ ⋅ T d V = ∫ S T ⋅ n ^ d S = ∫ S T ⋅ d S ⃗ ∫ V T i j , j d V = ∫ S T i j n ^ j d S = ∫ S T i j d S j \int_V \nabla_{\vec x}\cdot Td V = \int_S T\cdot \hat ndS = \int _S T\cdot d\vec S \\ \int_V T_{ij,j}dV=\int_ST_{ij}\hat n_jdS=\int_ST_{ij}dS_j Vx TdV=STn^dS=STdS VTij,jdV=STijn^jdS=STijdSj

通过利用散度定理,也可以证明:
在这里插入图片描述
其中,假设 δ i k , j = 0 i k j \delta_{ik,j}=0_{ikj} δik,j=0ikj,另外,因为 x k , j = δ k j x_{k,j}=\delta_{kj} xk,j=δkj,得到:
在这里插入图片描述
给定一个定义在域B的二阶张量,以下成立:
在这里插入图片描述
因此,可以证明:
在这里插入图片描述

问题1.47 令 Ω \Omega Ω 是一个边界为 Γ \Gamma Γ的域,m是一个二阶张量场, ω \omega ω是一个标量场,证明:

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

路径的独立性

连接两个点A和B的曲线称为从A到B的路径
在这里插入图片描述
建立线积分与路径无关的条件:
b ⃗ ( x ⃗ ) \vec b(\vec x) b (x )是一个连续的向量场,那么积分 ∫ C b ⃗ ⋅ d r ⃗ \int_{C}\vec b\cdot d\vec r Cb dr 与路径无关,当且仅当 b ⃗ \vec b b 是保守场,这意味着存在一个标量场 ϕ \phi ϕ使得 b ⃗ = ∇ x ⃗ ϕ \vec b = \nabla_{\vec x} \phi b =x ϕ,也就是向量场是某个标量场的梯度:
∫ A B b ⃗ ⋅ d r ⃗ = ∫ A B ∇ x ⃗ ϕ ⋅ d r ⃗ ∫ A B ( b 1 e ^ 1 + b 2 e ^ 2 + b 3 e ^ 3 ) ⋅ d r ⃗ = ∫ A B ( ∂ ϕ ∂ x 1 e ^ 1 + ∂ ϕ ∂ x 2 e ^ 2 + ∂ ϕ ∂ x 3 e ^ 3 ) ⋅ d r ⃗ \int_A^B\vec b\cdot d\vec r=\int_A^B\nabla_{\vec x}\phi \cdot d\vec r \\ \int_A^B(b_1\hat e_1+b_2\hat e_2+b_3\hat e_3)\cdot d\vec r=\int_A^B(\frac{\partial \phi}{\partial x_1}\hat e_1+\frac{\partial \phi}{\partial x_2}\hat e_2+\frac{\partial \phi}{\partial x_3}\hat e_3)\cdot d\vec r ABb dr =ABx ϕdr AB(b1e^1+b2e^2+b3e^3)dr =AB(x1ϕe^1+x2ϕe^2+x3ϕe^3)dr

因此:
b 1 = ∂ ϕ ∂ x 1 ; b 2 = ∂ ϕ ∂ x 2 ; b 3 = ∂ ϕ ∂ x 3 ; b_1 = \frac{\partial \phi}{\partial x_1}; \quad b_2 = \frac{\partial \phi}{\partial x_2}; \quad b_3 = \frac{\partial \phi}{\partial x_3}; b1=x1ϕ;b2=x2ϕ;b3=x3ϕ;

由于场是保守的,所以 b ⃗ \vec b b 的旋度为0:
在这里插入图片描述
因此,得出结论:
在这里插入图片描述

因此,如果以上条件不满足,则场不是保守的

Kelvin-Stokes定理

令S是一个曲面, F ⃗ ( x ⃗ , t ) \vec F(\vec x, t) F (x ,t)是一个向量场,根据Kelvin-Stokes定理,有:
∮ Γ F ⃗ ⋅ d Γ ⃗ = ∫ Ω ( ∇ ⃗ x ⃗ ∧ F ⃗ ) ⋅ d S ⃗ = ∫ Ω ( ∇ ⃗ x ⃗ ∧ F ⃗ ) ⋅ n ^ d S \boxed{\oint_{\Gamma}\vec F \cdot d\vec \Gamma=\int_{\Omega}(\vec \nabla_{\vec x}\wedge \vec F)\cdot d\vec S=\int_{\Omega}(\vec \nabla_{\vec x}\wedge \vec F)\cdot \hat n dS} ΓF dΓ =Ω( x F )dS =Ω( x F )n^dS

在这里插入图片描述

如果 p ^ \hat p p^表示切向于边界 Γ \Gamma Γ的单位向量,那么Stoke’s定理为:
∮ Γ F ⃗ ⋅ p ^ d Γ = ∫ Ω ( ∇ ⃗ x ⃗ ∧ F ⃗ ) ⋅ d S ⃗ = ∫ Ω ( ∇ ⃗ x ⃗ ∧ F ⃗ ) ⋅ n ^ d S \oint_{\Gamma}\vec F\cdot \hat pd \Gamma=\int_{\Omega}(\vec \nabla_{\vec x}\wedge \vec F)\cdot d\vec S=\int_{\Omega}(\vec \nabla_{\vec x}\wedge \vec F)\cdot \hat ndS ΓF p^dΓ=Ω( x F )dS =Ω( x F )n^dS

参考笛卡尔坐标系下的表示:
F ⃗ = F 1 e ^ 1 + F 2 e ^ 2 + F 3 e ^ 3 \vec F= F_1 \hat e_1+ F_2 \hat e_2+ F_3 \hat e_3 F =F1e^1+F2e^2+F3e^3
d S ⃗ = d S 1 e ^ 1 + d S 2 e ^ 2 + d S 3 e ^ 3 d\vec S = dS_1\hat e_1+dS_2\hat e_2+dS_3\hat e_3 dS =dS1e^1+dS2e^2+dS3e^3
d Γ ⃗ = d x 1 e ^ 1 + d x 2 e ^ 2 + d x 3 e ^ 3 d\vec \Gamma = dx_1\hat e_1+dx_2\hat e_2+dx_3\hat e_3 dΓ =dx1e^1+dx2e^2+dx3e^3

F ⃗ \vec F F 的旋度的分量为:
在这里插入图片描述
那么,Stoke‘s定理表示成以上分量形式:
∮ Γ F 1 d x 1 + F 2 d x 2 + F 3 d x 3 = ∫ Ω ( ∂ F 3 ∂ x 2 − ∂ F 2 ∂ x 3 ) d S 1 + ( ∂ F 1 ∂ x 3 − ∂ F 3 ∂ x 1 ) d S 2 + ( ∂ F 2 ∂ x 1 − ∂ F 1 ∂ x 2 ) d S 3 \oint_{\Gamma}F_1dx_1+F_2dx_2+F_3dx_3 \\ =\int_{\Omega}(\frac{\partial F_3}{\partial x_2}-\frac{\partial F_2}{\partial x_3})dS_1+(\frac{\partial F_1}{\partial x_3}-\frac{\partial F_3}{\partial x_1})dS_2+(\frac{\partial F_2}{\partial x_1}-\frac{\partial F_1}{\partial x_2})dS_3 ΓF1dx1+F2dx2+F3dx3=Ω(x2F3x3F2)dS1+(x3F1x1F3)dS2+(x1F2x2F1)dS3

特殊情况:曲面S是就是平面 Ω \Omega Ω,上式仍然成立
在这里插入图片描述
如果 Ω \Omega Ω x 1 − x 2 x_1-x_2 x1x2平面,那么上式变为:
∮ Γ F ⃗ ⋅ d Γ = ∫ Ω ( ∇ ⃗ x ⃗ ∧ F ⃗ ) ⋅ e ^ 3 d S \oint_{\Gamma}\vec F \cdot d \Gamma=\int_{\Omega}(\vec \nabla_{\vec x}\wedge \vec F)\cdot \hat e_3 dS ΓF dΓ=Ω( x F )e^3dS

这就是大家都知道的平面上的Stoke’s定理,即格林公式
张量分量:
∮ Γ F 1 d x 1 + F 2 d x 2 = ∫ Ω ( ∂ F 2 ∂ x 1 − ∂ F 1 ∂ x 2 ) d S 3 \oint_{\Gamma}F_1dx_1+F_2dx_2=\int_{\Omega}(\frac{\partial F_2}{\partial x_1}-\frac{\partial F_1}{\partial x_2}) dS_3 ΓF1dx1+F2dx2=Ω(x1F2x2F1)dS3

在这里插入图片描述

格林公式

F ⃗ \vec F F 是一个向量场,应用散度定理,有:
∫ V ∇ x ⃗ ⋅ F ⃗ = ∫ S F ⃗ ⋅ n ^ d S \int_V \nabla_{\vec x}\cdot \vec F=\int_S \vec F \cdot \hat n dS Vx F =SF n^dS

根据:
在这里插入图片描述
F ⃗ = ϕ ∇ x ⃗ ψ \vec F = \phi \nabla_{\vec x}\psi F =ϕx ψ,代入上面两个等式:
在这里插入图片描述
这就是格林第一公式

且有:
在这里插入图片描述
这就是格林第二公式

问题1.48 令 b ⃗ \vec b b 是一个向量场,定义为 b ⃗ = ∇ ⃗ x ⃗ ∧ v ⃗ \vec b =\vec \nabla_{\vec x}\wedge \vec v b = x v ,证明:

在这里插入图片描述
参考教材:
Eduardo W.V. Chaves, Notes On Continuum Mechanics

这篇关于【连续介质力学】涉及积分的定理的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/809269

相关文章

uva 1342 欧拉定理(计算几何模板)

题意: 给几个点,把这几个点用直线连起来,求这些直线把平面分成了几个。 解析: 欧拉定理: 顶点数 + 面数 - 边数= 2。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#inc

STL经典案例(四)——实验室预约综合管理系统(项目涉及知识点很全面,内容有点多,耐心看完会有收获的!)

项目干货满满,内容有点过多,看起来可能会有点卡。系统提示读完超过俩小时,建议分多篇发布,我觉得分篇就不完整了,失去了这个项目的灵魂 一、需求分析 高校实验室预约管理系统包括三种不同身份:管理员、实验室教师、学生 管理员:给学生和实验室教师创建账号并分发 实验室教师:审核学生的预约申请 学生:申请使用实验室 高校实验室包括:超景深实验室(可容纳10人)、大数据实验室(可容纳20人)、物联网实验

微积分-积分应用5.4(功)

术语“功”在日常语言中用来表示完成一项任务所需的总努力量。在物理学中,它有一个依赖于“力”概念的技术含义。直观上,你可以将力理解为对物体的推或拉——例如,一个书本在桌面上的水平推动,或者地球对球的向下拉力。一般来说,如果一个物体沿着一条直线运动,位置函数为 s ( t ) s(t) s(t),那么物体上的力 F F F(与运动方向相同)由牛顿第二运动定律给出,等于物体的质量 m m m 与其

Java验证辛钦大数定理

本实验通过程序模拟采集大量的样本数据来验证辛钦大数定理。   实验环境: 本实验采用Java语言编程,开发环境为Eclipse,图像生成使用JFreeChart类。   一,验证辛钦大数定理 由辛钦大数定理描述为: 辛钦大数定理(弱大数定理)  设随机变量序列 X1, X2, … 相互独立,服从同一分布,具有数学期望E(Xi) = μ, i = 1, 2, …, 则对于任意正数ε ,

CPC23三 K.(Lucas定理)

K.喵喵的神·数 Time Limit: 1 Sec Memory Limit: 128 MB Description 喵喵对组合数比较感兴趣,并且对计算组合数非常在行。同时为了追求有后宫的素质的生活,喵喵每天都要研究质数。 我们先来复习一下什么叫做组合数。对于正整数P、T 然后我们再来复习一下什么叫质数。质数就是素数,如果说正整数N的约数只有1和它本身,N

变速积分PID控制算法

变速积分PID控制算法 变速积分PID控制算法:变速积分PID的基本思想:变速积分的PID积分项表达式: 注:本文内容摘自《先进PID控制MATLAB仿真(第4版)》刘金琨 编著,研读此书受益匪浅,感谢作者! 变速积分PID控制算法: 在普通的PID控制算法中,由于积分系数 k i k_i ki​是常数,所以在整个控制过程中,积分增量不变。而系统对积分项的要求是,系统偏差大

梯形积分PID控制算法

梯形积分PID控制算法 梯形积分PID控制算法: 注:本文内容摘自《先进PID控制MATLAB仿真(第4版)》刘金琨 编著,研读此书受益匪浅,感谢作者! 梯形积分PID控制算法: 在PID控制律中积分项的作用是消除余差,为了减小余差,应提高积分项的运算精度,为此,可将矩形积分改为梯形积分。梯形积分的计算公式: ∫ 0 t e ( t ) d t = ∑ i = 0 k e

抗积分饱和PID控制算法

抗积分饱和PID控制算法 抗积分饱和PID控制算法:1.积分饱和现象:2.抗积分饱和算法: 注:本文内容摘自《先进PID控制MATLAB仿真(第4版)》刘金琨 编著,研读此书受益匪浅,感谢作者! 抗积分饱和PID控制算法: 1.积分饱和现象: 所谓积分饱和现象是指若系统存在一个方向偏差,PID控制器的输出由于积分作用的不断累加而加大,从而导致执行机构到达极限位置 X m

积分分离PID控制算法

积分分离PID控制算法 积分分离PID控制:积分分离控制基本思路:积分分离控制算法表示:积分分离式PID控制算法程序流程图: 注:本文内容摘自《先进PID控制MATLAB仿真(第4版)》刘金琨 编著,研读此书受益匪浅,感谢作者! 积分分离PID控制: 在普通的PID控制中引入积分环节的目的,主要为了消除静差,提高控制精度。但在过程启动、结束或大幅度增减设定时,短时间内系统输出

量化交易面试:什么是中心极限定理?

中心极限定理(Central Limit Theorem, CLT)是概率论和统计学中的一个重要定理,它描述了在一定条件下,独立随机变量的和的分布趋向于正态分布的性质。这个定理在量化交易和金融分析中具有重要的应用价值。以下是对中心极限定理的详细解释: 基本概念: 中心极限定理指出,当我们从一个具有任意分布的总体中抽取足够大的样本时,样本均值的分布将近似于正态分布,无论原始总体的分布是什么样的。