一维波动方程推导过程

应力波反射法检测基桩原理

1.1 基桩动测技术的发展及国内外研究现状

一百年以前，动力打桩公式 1865年B.de Saint Venant提出一维波动方程 50年代后期A.Smith提出了波动方程在桩基中应用的差分数值 解法，它把锤一桩一土系统简化为质量块、弹簧和阻尼器模型 从而使波动方程打桩分析进入实用阶段。

1967年美国G.G.Goble等人发表了“关于桩承载力的动测研究”一文， 1975年发表了“根据动测确定桩的承载力”研究报告 1970年以后，美国己把动力试桩技术用于实际工程 1977年PDI公司开始生产以PDA(Pile Driving Analyzer)打桩分析仪 采用波动方程程序(Case Pile Wave-equation Analysis program/contimuous,简CAPWAPC程序)对桩的侧阻分布、端阻和桩身缺陷

进行实测波形的拟合法分析。

方便、快捷、一定的准确度被各国接受 要求较高的人员素质、专业理论知识、 丰富的工程经验 缺乏与静荷载试验在桩周分层摩阻力和端阻力方面对比。

1.2.1 一维杆的纵向波动方程

一根材质均匀的等截面弹性杆，长度为L，截

面积为A，弹性模量为E，体密度为ρ 。若杆变

形时符合平截面假定，在

化工热力学

3.4 酶促反应动力学 酶促反应动力学酶促反应动力学(kinetics of enzymecatalyzed reactions)是研究酶促反应速度 及其影响因素的科学。酶促反应的影响因素 主要包括酶的浓度、底物的浓度、pH、温度、 抑制剂和激活剂等。

化工热力学

一. 酶浓度的影响在一定温度和pH下，酶 促反应在底物浓度大于 100 Km时，速度与酶的浓 度呈正比。 酶浓度对速度的影响机 理：酶浓度增加，[ES]也 增加，而V=k3[ES],故反 应速度增加。

化工热力学

二. 温度对酶促反应速度的影响 酶促反应与其它化学反应一样，随温度的增加，反应 速度加快。化学反应中温度每增加10℃反应速度增加的 倍数称为温度系数Q10。一般的化学反应的Q10为2~3,而 酶促反应的Q10为1~2。 在一定范围内，反应速度达到最大时对应的温度称为 该酶促反应的最适温度(optimum temperature Tm ).一 般动物组织中的酶其最适温度为35~40℃,植物与微生物 中的酶其最适温度为30~60℃,少数酶可达60℃以上，如 细菌淀粉水解酶的最适温度90℃以上。

化工热力学

温度对酶促反应速度的影响机理：

1. 温度影响反应体系中的活化分子

第七章 一维波动方程的傅里叶解 小结及习题答案

第二篇 数学物理方程

——物理问题中的二阶线性偏微分方程及其解法

Abstracts:1、根据物理问题导出数理方程—偏微分方程；

2、给定数理方程的附加条件：初始条件、边界条件、物理条件 (自然条件，连接条件)，从而与数理方程一起构成定解问题； 3、方程齐次化；

４、数理方程的线性导致解的叠加。

一、数理方程的来源和分类（状态描述、变化规律）

1、来源

I．质点力学：牛顿第二定律F?mr

??弦?2u(r,t)???a2?2u(r,t)?0(波动方程);?杆 振动：2?弹性体力学(弹性定律)?t?膜连续体力学? ??????流体力学：质量守恒律：???(?v)?0;?t??v?1?热力学物态方程：?(v??)v?p?f?0(Euler eq.).??t????II.麦克斯韦方程

?D?d???d????D??;????E?dl???B?ds???E?B;???????B?d??0???B?0;?H?dl???(j?D)?ds???H?j?D. ??E???u,B???A,u,A满足波动方程。???Lorenz力公式?力学方程；Maxwell eqs.+电导定律?电报方程。III. 热力学统计

第二篇 数学物理方程

——物理问题中的二阶线性偏微分方程及其解法

Abstracts:1、根据物理问题导出数理方程—偏微分方程；

2、给定数理方程的附加条件：初始条件、边界条件、物理条件

(自然条件，连接条件)，从而与数理方程一起构成定解问题；

3、方程齐次化；

４、数理方程的线性导致解的叠加。

一、数理方程的来源和分类（状态描述、变化规律）

1、来源

I ．质点力学：牛顿第二定律F mr = 连续体力学2222()(,)(,)0(()0;v 1()0(Euler eq.).u r t a u r t t v t v v p f t ρρρ?????-?=?????????+??=????-?+??=+=?????

弹性定律弦弹性体力学

杆 振动：波动方程);膜流体力学：质量守恒律：热力学物态方程： II.麦克斯韦方程

;;00;().,,,D D E l B s E B B B H l j D s H j D E u B A u A σρτρσ??=???=?=????=????=???=?=+????=+??=-?=?????????

???????????d d d d d d d 满足波动方程。Lorenz 力公式力学方程；Maxwell eqs.+电导定律电

本文详细推导了高等传热学的基本方程：连续性方程、动量方程和能量方程

方程推导

1.导热微分方程

x方向导入微元体的热流量为 x

x+dx方向导出微元体的热流量为： Tdydz x

x dx x x Tdx x ( dydz)dx x x x

同理可得y、z方向的导入、导出热流量。

根据能量守恒：导入微元体的总热流量+微元体内的生

成热=导出微元体的总热流量+微元体内能的增加 微元体内能的增加：dU c

Tdxdydz 微元体内的生成热：qdxdydz T T T T

经整理有： c z z q x x y y

该式可在（1）导热系数为常数；（2）导热系数为常数，无内热源（3）导热系数为常数、稳态（4）导热系数为常数、无内热源、稳态等情况下简化 T 1 T 1 T T

圆柱坐标系： c r 2 q r r r r z z

T 1 2 T 1 T 1 T

球坐标系： c r sin q r2 r r r2sin r2si

完全弹性 碰撞的速度公式推导过程

完全弹性碰撞的速度公式推导过程完全弹性碰撞的速度公式是怎么推导的无从得知，书上没讲，很多资料也没有讲，我想多半是为了不要影响思维的连贯性，所以将之省略了。我开始以为不复杂，就是上标下标看着烦人，所以就打算试着推导一下。谁知这个推导并没有想象中那么简单。第一次因为上下标搞混了，推导了半天没结果就放一边了。第二次仔细地推导，花了更多的时间，结果还是一塌糊涂。我终于明白书上为什么没有把这个推导过程放在书里了，的确是太复杂，学习的时候多半会干扰对碰撞本身的关注。但是这么放弃也有点不甘心，就又花了些时间，第三次准备将其推导出来。闲人可以看看，我也是放假闲着没事推导的，实在是很复杂很恐怖的推导。我自己都不想再看，因为象那样用常规的方式根本就推导不出来! 动量守恒定律: MpVp'+MqVq'=MpVp+MqVq(1-1) 动能守恒: (1/2)MpVp'2+(1/2)MqVq'2=(1/2)MpVp2+(1/2)MqVq2(1-2) 前两次推导吃了亏，所以第三次推导前仔细看了看书上结果公式的特点。有这样几个地方需要注意: 1、撞击后有两个速度，我们需要求的结果分别是这两个速度; 2、任一撞后的速度公式中，不能有另一个待求

从正态总体N（μ1，σ）和N（μ2，σ）中分别抽取含量为n1和n2的样本，两样本均数差值X1 －X

2 服从正态分布

N（μ1－μ2，?），其中

X1－X2?＝X1－X2?2(＋1n11) ① n2其中①式中σX1 －

X2 为两样本均数差值的标准误，其估计值为

n?n11SX－X＝SC2(＋)＝SC2(12) ② 12n1n2n1?n2其中②式中SC2为两样本合并的方差，其计算公式为：

?XSc?21?(X1)2/n1??X22?(?X2)2/n2n1?n2?22 ，则可用公式

③

如已计算出S1 和 S12③ 计算出

SX－X＝S2x?S2x2＝S21/n1?S22/n2④

1在H0：μ1＝μ2＝0的条件下，t的计算公式为：

t?|X1?X2|SX1?X2，ν＝n1?n2?2⑤

例3-3 测得14名慢性支气管炎病人与11名健康人的尿中17酮类固醇（u mol/24h）排出量如下，试比较两组人的尿中17酮类固醇的排出量有无不同。

病人X1：10.05 18.75 18.99 15.94 13.96 17.22 14.69 15.10 9.42

波动方程数值模拟技术及其应用

作者姓名: 陈睿 专业班级: 2008050603指导教师: 熊晓军

摘 要

波动方程数值模拟技术在地震勘探中的应用非常广泛，特别是对于碳酸盐岩这一类重要的油气储集层。

本文主要介绍了声学波动方程的基本理论，相位移波动方程数值模拟方法，相位移加插值波动方程数值模拟方法的原理，并且采用相位移加插值的方法进行实际碳酸盐岩模型的数值模拟，根据实际区域的地质剖面猜测初始的地震模型，通过波动方程对该猜测的初始模型进行正演与偏移，再把通过偏移的地震剖面与实际的地震记录剖面对比，反复调整其中的相关参数，更新地质剖面，从而获得更加正确的地质解释模型。对比地质模型与原始的地震资料，从而确定了猜测的正确性，为该地区以后的储层预测、地震资料解释提供了一定的参考价值。

综上的论述，本次研究为相同地震、地质条件下礁滩储层的波场特征认识积累了一些经验，为准确地进行礁滩储层预测奠定了一定的基础。

关键词：相位移 波动方程 数值模拟 偏移

I

Numerical Simulation Technology Of Wave

Equation And Its Application

Abstract:The numer

第五次作业（前三题写在作业纸上）

一、用有限差分方法求解一维非定常热传导方程，初始条件和边界条件见说明.pdf文件，热扩散系数α=const，

?T?2T??2 ?t?x1. 用Tylaor展开法推导出FTCS格式的差分方程

2. 讨论该方程的相容性和稳定性，并说明稳定性要求对求解差分方程的影响。 3. 说明该方程的类型和定解条件，如何在程序中实现这些定解条件。

4. 编写M文件求解上述方程，并用适当的文字对程序做出说明。（部分由网络搜索得

到，添加，修改后得到。） function rechuandaopde

%以下所用数据，除了t的范围我根据题目要求取到了20000，其余均从pdf中得来 a=0.00001;%a的取值 xspan=[0 1];%x的取值范围 tspan=[0 20000];%t的取值范围

ngrid=[100 10];%分割的份数，前面的是t轴的，后面的是x轴的 f=@(x)0;%初值

g1=@(t)100;%边界条件一 g2=@(t)100;%边界条件二

[T,x,t]=pdesolution(a,f,g1,g2,xspan,tspan,ngrid);%计算所调用的函数 [x,t]=meshgrid(x,t);

mesh

第五次作业（前三题写在作业纸上）

一、用有限差分方法求解一维非定常热传导方程，初始条件和边界条件见说明.pdf文件，热扩散系数α=const，

?T?2T??2 ?t?x1. 用Tylaor展开法推导出FTCS格式的差分方程

2. 讨论该方程的相容性和稳定性，并说明稳定性要求对求解差分方程的影响。 3. 说明该方程的类型和定解条件，如何在程序中实现这些定解条件。

4. 编写M文件求解上述方程，并用适当的文字对程序做出说明。（部分由网络搜索得

到，添加，修改后得到。） function rechuandaopde

%以下所用数据，除了t的范围我根据题目要求取到了20000，其余均从pdf中得来 a=0.00001;%a的取值 xspan=[0 1];%x的取值范围 tspan=[0 20000];%t的取值范围

ngrid=[100 10];%分割的份数，前面的是t轴的，后面的是x轴的 f=@(x)0;%初值

g1=@(t)100;%边界条件一 g2=@(t)100;%边界条件二

[T,x,t]=pdesolution(a,f,g1,g2,xspan,tspan,ngrid);%计算所调用的函数 [x,t]=meshgrid(x,t);

mesh