下载文档原格式
里兹法（Ritz method）是一种求解微分方程的数值方法。
以下是使用里兹法求解微分方程的一般步骤：
确定定解条件：里兹法需要确定微分方程的定解条件，包括边界条件和初始条件。
这些条件将用于构建能量积分方程。
选择基函数：在里兹法中，选择一组适当的基函数是非常关键的。
基函数应该满足边界条件，并且能够较好地近似微分方程的解。
常用的基函数包括多项式、三角函数等。
构建能量积分方程：使用选择的基函数，将微分方程转化为能量积分方程。
这个过程涉及到对微分方程进行乘积积分，并应用边界条件。
求解线性方程组：通过构建能量积分方程，得到一个线性方程组。
该方程组的解将给出微分方程的近似解。
使用数值方法（如高斯消元法、LU分解等）求解该线性方程组。
评估解的精度：得到近似解后，需要评估解的精度。
可以通过计算残差、误差估计等方法来评估解的精度，并根据需要改进基函数或增加基函数的数量来提高精度。
需要注意的是，里兹法是一种近似方法，得到的解是微分方程的近似解。
因此，在实际应用中，需要根据问题的要求和解的精度来选择合适的基函数和求解方法。
同时，里兹法的收敛性和误差分析也是需要考虑的重要方面。
里兹法求微分方程1. 引言微分方程是数学中的重要概念，广泛应用于物理、工程、经济等领域。
里兹法（Ritz method）是一种常用的数值解微分方程的方法，通过将待求函数表示为一组已知函数的线性组合，将微分方程转化为一个待求系数的优化问题。
本文将介绍里兹法求解微分方程的基本原理和具体步骤，并通过一个实例来说明其应用。
2. 里兹法的基本原理里兹法的基本思想是将待求函数表示为一组已知函数的线性组合，通过确定系数的值来近似满足微分方程。
这里的已知函数称为试探函数（trial functions），可以根据具体问题的特点选择不同的试探函数。
假设我们要解的微分方程为：L[u(x)]=f(x)其中，u(x)是待求函数，L是微分算子，f(x)是已知函数。
我们将待求函数表示为试探函数的线性组合：u(x)=∑c ini=1ϕi(x)其中，c i是待定系数，ϕi(x)是已知函数。
将这个表示代入微分方程中，得到：L[∑c ini=1ϕi(x)]=f(x)我们可以将方程两边乘以试探函数的权重函数w(x)，然后对方程两边进行积分，得到：∫L b a [∑c ini=1ϕi(x)]w(x)dx=∫fba(x)w(x)dx这里的a和b是微分方程的区间。
将积分运算符和微分算子交换顺序，得到：∑c i ni=1∫Lba[ϕi(x)]w(x)dx=∫fba(x)w(x)dx我们可以定义一个误差函数R(c1,c2,...,c n)，表示上式左边减去右边的结果。
我们的目标是找到一组系数c i，使得误差函数最小化。
这样，原先的微分方程就被转化为了一个优化问题。
3. 里兹法的求解步骤里兹法的求解步骤如下：步骤1：选择试探函数根据待求函数的特点，选择一组试探函数。
常见的试探函数包括多项式、三角函数、指数函数等。
步骤2：计算误差函数将试探函数的表示代入微分方程，并乘以权重函数w(x)，然后对方程两边进行积分，得到误差函数R(c1,c2,...,c n)。
里兹法里兹法是近似计算的经典方法。
它是势能驻值原理具体应用的典型范例。
里兹法的基本思路是：将无限自由度体系近似地用有限自由度体系来代替，应用势能原理求得代用体系的精确解，从而求得原体系的近似解。
里兹法的具体做法是：选择一组满足求解域位移边界条件的试函数作为实际问题的近似解。
显然，近似解的精度与试函数的选择有关，如果精确解包含在试函数族中，由里兹法将得到精确解。
例1试用里兹法求图1.4所示悬臂梁的挠度方程。
设梁线弹性，EI为常数。
解：1.选取可能位移状态此悬臂梁应满足的位移边界条共有两个，即在固定端A处的挠dv应为零。
符合这两个条件的可能位移状态有无限多个，度v和转角dx这是一个无限自由度体系。
在近似分析中，我们假设挠曲线为一个多项式v = a 1x 2+ a 2x 3 + … + a n x 1+n (1)由于在x=0处要满足v=dxdv=0的条件，故在式（1）中没有包含常数项和一次项。
式（1）中共有n 个任意参数a 1，a 2，…，a n ，只要这n 个参数定义了，梁的挠度方程也就确定了，所以梁的变形决定于这n 个任意参数。
这里决定变形状态的参数个数就是体系的自由度。
采用式（1）所表示的多项式，就相当于把原来的无限自由度体系近似地作为n 个自由度体系来看待。
2.按单自由度体系计算在式（1）中只取第一项，即挠度方程为v = a 1x 2 (2)这时把梁按单自由度体系计算，体系势能为∏ = U － P v B其中222110()(2)222llEI EIU v dx a dx EIla ''===⎰⎰21B v a l =因此 2112EIla Pl a ∏=- 由势能驻值原理有21140d EIla Pl da ∏=-= 求得14Pla EI=代入式（2），得24Plv x EI=B 点的挠度为34BPl v EI= 与精确解33BPl v EI=相比，误差较大。
3.按两个自由度体系计算为了提高精度，在式（2）中保留前两项v = a 1x 2 + a 2x 3 (3)即把梁当做两个自由度体系看待，这时22221211220(26)2(33)2l EI U a a x dx EIl a la a a l =+=++⎰ 222231122122(33)()EIl a la a l a P a l a l ∏=++-+势能驻值条件为21212312202(23)002(36)0EIl a la Pl a EIl a la Pl a ∂∏⎧=+-=⎪∂⎪⎨∂∏⎪=+-=⎪∂⎩由此可得12Pl a EI =,26Pa EI=-代入式(3),挠度方程为23(3)6Pv lx x EI=- （4）式（4）实际上就是挠度的精确解。
在线性多自由度系统振动中，振动问题归结为刚度矩阵和质量矩阵的广义特征值问题，缺点：当系统自由度较大时，求解计算工作量非常大。
本章介绍邓克利法，瑞利法，里茨法，传递矩阵法等计算方法，可作为实用的工程计算方法对系统的振动特性作近似计算。
1、邓克利法由邓克利（Dunkerley ）在实验确定多圆盘的横向振动固有频率时提出的，便于作为系统基频的计算公式 。
自由振动作用力方程：0KX XM =+ n R ∈X 左乘柔度矩阵F = K -1，位移方程：0X X FM =+ 定义D=FM 为系统的动力矩阵：0X XD =+ 作用力方程的特征值问题：φφM K 2ω= 位移方程的特征值问题：φφλ=D 特征值：22221n ωωω<<< ，n λλλ>>> 21 关系：2/1i i ωλ=位移方程的最大特征根：211/1ωλ=，对应着系统的第一阶固有频率。
位移方程的特征方程：0=-I D λ展开：0)()1(1111=++++---n n n nna a a λλλD tr d d d a nn -=+++-=)(22111例：022211211=--λλd d d d0)]()([)1(21122211221122=-++--d d d d d d λλ当 M 为对角阵时：)(FM D tr tr =∑==ni iii m f 1特征方程又可写为：0)())((21=---n λλλλλλ有：∑=-=ni i a 11λtrD -=∑=-=ni i ii m f 1∑∑===ni iii ni im f 11λ∑∑===ni i ii ni im f 1121ω如果只保留第 i 个质量，所得的单自由度系统的固有频率为：iii ii i m f m k 12==ω例：两自由度系统柔度矩阵：⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+=2111111111k k k k kF （1）只保留 m 1 时1111k f =，1121m k =ω（2）只保留 m 2 时122122111k k k f =+=，21222m k =ω将2i ω代入：22221121111nni iωωωω+++=∑=对于梁结构系统，第二阶及第二阶以上的固有频率通常远大于基频，因此左端可只保留基频项，有：22221211111nωωωω+++≈例：三自由度系统⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡0002223101220010*********x x x k x x x m 采用常规方法，固有频率：m k /3730.01=ω，m k /3213.12=ω，m k /0286.23=ω邓克利法当 m1 单独存在时：m k /21=ω 当 m2 单独存在时：k k k k k k 21212112=+=，m k /1222=ω当 m3 单独存在时：kk k k k 251111321123=++=，52123k k =，mk 523=ω代入邓克利法公式：22221211111nωωωω+++≈，mk /3535.01=ω2、瑞利法瑞利法是基于能量原理的一种近似方法，可用于计算系统的基频，算出的近似值为实际基频的上限，配合邓克利法算出的基频下限，可以估计实际基频的大致范围。
Ritz法1.Raleigh法的原理Ritz里茨法是由Raleigh瑞利法发展而来的，学习里茨法需要先理解瑞利法。
1.瑞利法的计算步骤？（1）选取满足位移边界条件容许函数（容许函数可理解为假设的振型函数）；（2）利用容许函数计算最大动能和最大应变能；（3）利用最大动能和最大应变能计算最低阶自由振动频率。
2.为什么计算最大动能和最大应变能？根据能量守恒，对于变形体的自由振动，如果不考虑能量转化为温度等其他能量，那么自由振动的最大动能应该与最大应变能相等。
如最简单的质量块——弹簧系统，给定初始位移后释放，质量块会在一定范围内来回振荡。
有这样一个物理事实，在位移最大处速度为零，位移最大（即变形最大），此时动能为零，变形能最大；与之相反，位移为零（即没有变形），变形能为零，此时速度最大即动能最大。
因为能量守恒，只考虑动能和变形能时，最大动能等于最大变形能。
3. 计算最大动能和最大应变能后，如何计算频率？（1）如果使用瑞利法计算质量块——弹簧系统的自由振动频率，有该系统的应变势能可表示为212U kx ，动能可表示为212T mx。
（2）式中x 为质量块的位移，对于该线性系统，存在解析解形式如 sin x A t ，其中A 为振幅，ω为频率，φ为初相位。
（3）由解析解的形式可知最大位移x A ，求导一次得最大速度为xA 。
（4）将最大位移与最大速度，代入能量中得最大应变能与最大动能 2222max max 111,222U kA T m A mA ，观察频率存在的形式。
（5）提取出频率，满足以下关系2221212kA k m mA （6）瑞利法只能求一阶自由振动频率。
4. 计算频率后如何计算振型？（1）如果是瑞利法，假设的容许函数中一般不含有待定系数，则容许函数就是振型；（2）如是里茨法，容许函数中含有待定系数，求出频率后，得到的特征向量即为待定系数，特征向量代回容许函数就是振型。
2. Ritz 法因为瑞利法只能计算最低阶的自由振动频率，并且由于容许函数的选择一般不是精确的振型，不能完全满足所有边界条件，所及计算出的一阶自由振动频率会偏大。
里茨－加廖金方法里茨－加廖金方法Ritz-Galerkin method求解数学物理方程的近似方法，主要用于椭圆型边值问题,W.里茨于1908年对此做了开创性的工作这类方法从变分原理出发，选定有限个试探函数,,…,，用它们的线性组合[438-04]构造近似解,[kg1]从而把问题归结为确定组合中的系数。
极小值原理表达物理基本定律的一种形式，其表达可概括如下：给出一个依赖于物理状态（为一函数）的变量()（数学上称为泛函），同时给出()的容许函数集，即一切容许取的物理状态,则真实的物理状态就是中使()达到极小值的函数。
例如弹性力学中著名的极小能量原理的表述是：弹性体在外力作用下的平衡位移使总势能达到极小。
这里，总势能是位移函数的泛函。
现讨论小变形的均匀弹性膜，膜在区域上的垂直位移用函数(,)表示,假定在边界上膜固定在坐标平面上，即[438-05]， (1)在外力(,)作用下弹性膜的总势能为[438-06](2)它的容许函数集是满足边界条件(1)并使积分(2)存在的全体函数。
由极小能量原理得出：弹性膜的平衡位移是中使总势能(2)达到极小的函数，即[438-07]。
(3)虚功原理又称虚位移原理，在力学中与极小能量原理同属变分原理。
它的一般表述是：平衡系的力对虚位移所作的虚功为零。
在弹性膜的例子中,如果仍以表示平衡位移，则虚功原理可表达为：对所有[kg1][kg1]，使[438-08] (4)成立。
变分问题(3)或(4)确定的平衡位移,也是泊松方程第一边值问题的解，即满足[438-09]式中为拉普拉斯算子。
里茨法从与微分方程问题等价的极小值原理出发，选择有限个试探函数[kg1],,…,，在它们的线性组合[439-01]中去找近似解。
一般而言，试探函数须属于容许函数集把[439-02]代入(2)，得到[439-03](5)式中[439-04]，[439-05]。
由多元微分学可推出极小解[439-06]的系数应满足[439-07](6)故问题归结为求解线代数方程组(6)。
微分方程的Ritz近似方程1. 引言微分方程是数学中的重要工具，用于描述自然界中许多现象和过程。
然而，解析求解微分方程并不总是可行的，特别是对于复杂的非线性微分方程或者边界条件不明确的问题。
为了克服这些困难，人们发展了一种近似求解微分方程的方法，其中Ritz近似方程就是其中之一。
Ritz近似法基于变分原理，通过将原问题转化为一个极值问题来求解微分方程。
在这个极值问题中，我们引入一个试验函数，并通过调整试验函数中的参数来最小化误差函数。
Ritz近似法在实际应用中取得了广泛的成功，并被广泛应用于结构力学、流体力学、电磁场等领域。
2. Ritz方法的基本原理Ritz方法是一种基于变分原理的求解偏微分方程近似解的方法。
其基本思想是，在给定边界条件下，在某个函数空间中寻找一个函数使得误差最小化。
具体而言，Ritz方法通过以下步骤进行：1.假设待求解的函数形式，并引入待定参数；2.将待求解的函数代入原微分方程，得到一个含有待定参数的积分表达式；3.引入一个试验函数，并确定试验函数中的参数；4.将试验函数代入积分表达式中，得到一个关于参数的函数；5.调整参数使得该函数取极值，即最小化误差；6.求解参数，得到近似解。
Ritz方法的关键在于选择合适的试验函数。
一般来说，试验函数需要满足边界条件，并且足够接近真实解。
常用的试验函数包括多项式、三角函数、指数函数等。
3. Ritz方法在常微分方程中的应用Ritz方法可以用于求解各种类型的常微分方程，例如一阶和二阶常微分方程。
下面以一阶常微分方程为例进行说明。
考虑一阶常微分方程：y′(x)=f(x,y(x))其中f(x,y)是已知的函数。
我们希望求解上述方程在某个区间[a,b]上的近似解。
为了使用Ritz方法，我们假设待求解的近似解为：y N(x)=∑c iNi=1ϕi(x)其中ϕi(x)是试验函数，c i是待定系数。
将上述形式的近似解代入微分方程中，得到：∑c i Ni=1ϕi′(x)=f(x,∑c iNi=1ϕi(x))对上述等式两边同时乘以一个试验函数ϕj(x)并在区间[a,b]上积分，可以得到：∫ϕj ba (x)∑c iNi=1ϕi′(x)dx=∫ϕjba(x)f(x,∑c iNi=1ϕi(x))dx利用分部积分法和试验函数的正交性质，可以将上述等式转化为一个线性方程组：∑A ijNi=1c i=b j其中A ij是由试验函数的导数和非线性项的乘积计算得到的系数矩阵，b j是由边界条件和非线性项计算得到的向量。
结构动力学--瑞利法和里兹法研究报告周靖丰 学号：15213736一、 瑞利法和里兹法简介弹性系统振动问题归结为对系统的刚度矩阵和质量矩阵的研究。
然而，从微分方程出发，研究弹性体的振动，除了一些简单的情况以外，要精确求解往往是不可能的，而工程中遇到的实际结构总是比较复杂，因此近似解法占有很重要的地位。
通过对截止模态的研究发现对低频率固有频率的研究具有重要的意义，这对工程实践具有重要意义。
瑞利法和里兹法基于能量守恒原理，可以用来计算振动系统固有频率的的近似值。
二、 推导过程和算例多自由度系统 ➢ 瑞利法讨论自由度为n 的保守系统。
设系统做某阶主振动，则系统的动能和势能为：1T 2=T ••X M X 1V 2=T X KX设解的形式为：sin()X A t ωθ=+，带入到上式中，求出动能和势能的最大值为：2max 12T ω=T A MA ， max 12V =T A KAA 是振幅矩阵，根据保守系统机械能守恒原理，最大动能与最大势能应该相等，则有以下固有频率计算公式：2()R ω==T T A KAA A MAR （A ）称为瑞利商。
根据模态的定义有：2()i R ω==(i)T (i)(i)(i)T (i)ΦK ΦΦΦM Φ任选一个列阵ψ作为假设模态，这不是实际模态，但是总能表示为简正模态的线性组合：jN =a nj j a =Φ∑N ΨΦ其中，a 为系数a j （j=1，2，…,n ）组成的列阵。
取A=ψ，则有：22121()nj jj njj a R aω====∑∑TT TN N TT T NN a ΦK Φa a ΛaΨ=a ΦM Φa a Ea这样得到的瑞利商不是系统的任一阶固有频率的平方，但必介于系统的最高和最低的固有频率的平方21ω和2n ω之间：221n R ωω≤≤Ψ（），若恰当选择系数使假设模态接近k)Φ（，其中除k a 以外的其他系数j a （j ≠k ）均为小量，令j j k a =a j=ε（1，...,k-1,k+1,...,n ）带入到瑞利商中可以得到：n2222kj k j j=1R =+-ωωωε∑Ψ（）（）因此若假设模态Ψ与第K 阶模态ΚΨ（）的差别为一阶小量，则瑞利商与第K 阶固有频率平方2K ω的差别为二阶小量。
结构动力学--瑞利法和里兹法研究报告周靖丰学号：15213736—、瑞利法和里兹法简介弹性系统振动问题归结为对系统的刚度矩阵和质量矩阵的研究。
然而，从微分方程出发，研究弹性体的振动，除了一些简单的情况以外，要精确求解往往是不可能的，而工程中遇到的实际结构总是比较复杂，因此近似解法占有很重要的地位。
通过对截止模态的研究发现对低频率固有频率的研究具有重要的意义，这对工程实践具有重要意义。
瑞利法和里兹法基于能量守恒原理，可以用来计算振动系统固有频率的的近似值。
二、推导过程和算例多自由度系统?瑞利法讨论自由度为n的保守系统。
设系统做某阶主振动，则系统的动能和势能为:1 ? T? 1T -X M X V -X T KX2 2设解的形式为：X As in( t ），带入到上式中，求出动能和势能的最大值为T max A T MA ，V max 1A T KA2 2A是振幅矩阵，根据保守系统机械能守恒原理，最大动能与最大势能应该相等，则有以下固有频率计算公式：R（A）称为瑞利商。
根据模态的定义有:R（①①）①①T K①① 2 ①①T M①①i任选一个列阵作为假设模态，这不是实际模态，但是总能表示为简正模态的线性组合:na j N =①N a j其中，a为系数a j （j=1 , 2, -；n）组成的列阵。
取A=，则有:R(A) A T KA A T MAna T①N K①N a R（旳a T①N M①』 a T A a2 2a.j jj 1a T Ea n2a j率的平方12和2之间：12R （屮）2，若恰当选择系数使假设模态接近①代，其中除a k以外的其他系数a j （j k）均为小量，令a j = j a（j=1，…,k-1,k+1,..., n ）带入到瑞利商中可以得到：nR（旳=2+ （2- ；）2j=i因此若假设模态屮与第K阶模态M K）的差别为一阶小量，则瑞利商与第K阶固有频率平方K的差别为二阶小量。
