《相似模拟试验和数值模拟在岩土工程中的应用及实际案例》由会员分享，可在线阅读，更多相关《相似模拟试验和数值模拟在岩土工程中的应用及实际案例（26页珍藏版）》请在金锄头文库上搜索。
1、相似模拟与模型试验在岩土工程中的应用相似模拟与其它一样是社会生产发展的必然产物。由于社会生产的不断发展，岩土工程 所提出的问题日益复杂和繁琐。用数学方法很难得到精确的解析解，只能作一些假设与简化 再求解，因而带来一些误差。于是人们不得不通过实验的方法来探求那些靠数学方法无法研 究的复杂现象的规律性。但是直接的实验的方法有很大的局限性，其实验的结果只能推广到 与实验条件完全相同的实际问题中去，这种实验方法常常只能得出个别量的表面规律性关 系，难以抓住现象的内在本质。相似模拟正是为解决这些问题而产生的，它不直接的研 究自然现象或过程的本身，而是研究与这些自然现象或过程相似的模型， 它是理论与实际密
2、切相结合的科学研究方法，是解决一些比较复杂的生产工程问题的一种有 效方法。一、相似模拟与模型试验的方要研究内容它是研究自然界相似现象的一门科学。它提供了相似判断的方法。并用于指导模型试验， 整理试验结果，并把试验结果用于原型的理论基础。二、相似常数设c表示相似常数，x表示原型中的物理量，X表示模型中的物理量，贝I：xc =斗ix i其中c 表示第i个物理量所对应的相似常数。 i物理量包含于现象之中。而表示现象的物理量，一般都不是孤立的，互不关联的，而是 处在自然规律所决定的一定关系中，所以说各种相似常数之间也是相互关联的。在许多的情 况下这种关联表现为数学方程的形式。下面举例说明：设两个物体受
3、力与运动相似贝它们的质点的运动方程和力学方程均可用同一方程描述，即：原型的运动方程与物理方程dsdvv =f = mdtdt模型的运动方程与物理方程dsdvv =f = m dtdt因为两个物体的现象相似，其对应物理量互成比例，即 st 十c厂cts st t ，联合得到ccv t = c = 1cccf t 二 c 二 1 ccmv由，可以说明，各相似常数不是任意选择的，它们之间是相互关联的。三、相似三定理1. 相似第一定理相似第一定理是指出两个相似物体之间物理量的关系，具体可以归纳为二点。一、 相似现象可以用完全相同的方程组来表示。二、用来表征这些现象的一切物理量在空间相对 应的各点在时间
4、上相对应的各瞬间各自互成一定比例关系。2. 相似第二定理相似第二定理描述了物理体系中各个物理量之间的关系，相似准则之间的函数关系。兀关系式（准则方程）f （兀,兀,兀）=01 1 2 n兀关系式的性质 对于彼此相似的现象，兀关系式相同。 兀关系式中的兀项在模型试验中有自变项与应变项之分。自变项是由单值条件的物理量所组成的定性准则，应变项是包含非单值条件的物理量的非定性准则。 若能做到原型与模型中的自变兀项相等，由应变兀项与自变兀项之间的关系式可 以得到应变兀项，然后推广到原型中去，作为工程设计的各种参数。3. 相似第三定理相似第三定理是解决两个同类物理现象满足什么样的条件才能相似的问题。第一条
5、件：由于相似现象服从同一的自然规律，因此，可被完全相同的方程能所描述 第二条件：具有相同的文字方程式，其单值条件相似，并且从单值条件导出的相似准则 的数值相等。所谓的单值条件是指从一群现象中，根据某一个现象的特性，把这个具体的现象从一群 现象中区分出来的那些条件，单值条件中的物理量又称为单值量。单值条件包括几何条件、 物理条件、边界条件和初始条件。4. 相似三定理之间的关系相似第一和第二定理是从现象已经相似这一基础上出发来考虑问题，第一定理说明了相 似现象各物理量之间的关系，并以相似准则的形式表示出来。第二定理指出了各相似准则之 间的关系，便于将一现象的实验结果推广到其它现象。相似第三定理直接
6、同代表具体现象的 单值条件相联系，并且强调了单值量相似，所以显于出了科学上的严密性，是构成现象相似 的充要条件。是一切模型试验应遵守的理论指导原则。但是在一些复杂的现象中，很难确定现象的单值条件，仅能借经验判断何为系统最主要 的参量，或者虽然知道单值量，但是很难做到模型和原型由单值量组成的某些相似准则在数 值上的一致，这使得相似第三定理真正的实行，并因而使模型试验结果带来近似的性质。 一、 同类相似与异类相似同类相似是指相似的物体是同类物质，模型与原型的全部物理量相等，物理本质一致， 区别在于各物理量的大小比例不同。异类相似是指相似的物体不同类。仅因为对应量都遵循 相同的方程式，具有数学上的相
7、似性。五、相似准则的导出方法相似准则的导出方法有三种：定律分析法，方程分析法和因次分析法。从理论上说，三 种方法可以得到同样的结果，只是用不同的方法对物理现象作数学上的描述。但是作为三种 不同的方法，又有各自的适用条件。1. 三种方法的介绍定律分析：这种方法是建立在全部现象的物理定律已知的基础上的，通过剔除次要因素， 从而推算出数量足够的，反映现象实质的兀项。这种方法的缺点上：1) 流于就事论事，看不出现象的变化过程和内在联系，故作为一种方法，缺乏典型意 义2) 由于必须找出所有的物理定理，所以对于未能掌握其全部机理的，较为复杂的物理 现象，运用这种方法是不可能的，甚至无法找到近似解3) 常常
8、有一些物理定理，对于所讨论的问题表面上看去关系不密切，但又不宜于妄加 剔除，而必须通过实验找出各个定律间的制约关系，决定其重要因素，这实际问题 的解决带来不便。优点：对于模型制作有指导性意义。方程分析法：根据已知现象的微分或积分方程推出兀项。此方法的的优点：1) 结构严密，能反映出现象的本质，故可望得到问题的可靠性结论2) 分析程序明确步骤易于检查3) 各种成份的地位一览无遗，有利于推断，比较和校验缺点：对现象的机理不清楚，没有建立方程的问题，无法解决因次分析法：是根据正确选定参量，通过因次分析法考察各参量的因次，求出和兀定理 一致的函数关系式，并据此进行相似现象的推广。因次分析法的优点，对于
9、一切机理尚未彻 底弄清，规律也未充分掌握的现象来说，尤其明显。它能帮助人们快速地通过相似性实验核 定所选参量的正确性，并在此基础上不断加深人们对现象机理和规律性的认识。以上各种方法，日前应用最广泛的是因次分析法，但是也不排除将各种方法结合使用的 可能性。六、相似准则导出方法的解题步骤1. 三种方法的解题步骤1) 定律分析法的步骤 分析现象，抓住主要矛盾，排除次要因素 写出主要矛盾的物理表达式 作等效变化，转化为具有相同因次的物理量 两两作比值，求出相似准则兀2) 方程分析法通常的方程分析法有：相似转换法和积分类比法相似转换法的步骤 写出现象的基本微分方程 写出全部的单值条件，并令其二现象相似
10、将微分方程按不同现象写出 进行相似转换 求出相似准则兀积分类比法的步骤 写出现象的基本微分方程和全部的单值条件 用方程的任一项，除其它各项 进行积分类比转换，求出相应的准则3) 因次分析法因次分析法一般分为两种：指数分析法和矩阵分析法。这两种方法的基本原理一样，运算步骤稍有不同。指数分析法主要用于现象的物理量较少的情况，而矩阵分析法主要用于现 象物理量较多的情况。指数分析法 列出相似准则的表达式 根据方程两边因次相等列出物理量参数的方程K个 设物理量有M个，任选其中的MK个物理量为已知量 将这MK个物理量，依次用MK个单位向量代入方程，得到MK组解 把这MK组解代入相似准则的表达式中，可以得出
11、MK个独立的相似准则 矩阵分析法矩阵分析法与指数分析法的基本原理一样，矩阵分析法把线性方程组的求解用矩阵的求 解来代替。其运算步骤不再此重复。2. 证明指数分析法解出的独立兀项的广泛代表意义例设某现象由5个物理量Al, A2, A3, A4, A5组成，这5个基本物理的独立因次为L，M， N物理量的表达式Ai 二 M 卩7y;i 二 123,4,5相似准则的表达式兀=A1x A2 y A3 zA4uA5v因为兀项为零，故有对于 L a X +a Y + a Z + a U +a V 二 012345对于 M P X + p Y +p Z +p U + p V 二 012345对于 T y X
12、+y y + y z +y U +y V 二 012345固定U，V这两个参数，设U=0, V=1则可以得出一组解，设为X = X1, Y=Y1, Z = Z1, 但若设U=0, V=N则方程得出另一组解，设为X = X2, Y=Y2, Z = Z2这两组解之间存在着 如下关系，即：111X1 二 X 2 Y1 二 Y 2 Z1 =Z 2NNN由上式可知，这个相似准则和前一个相似准则只差方次关系，又因为相似准则可以通过加、减、乘、除、幂运算等进行相互变换，故这两个相似准则实为同一个无因次量群。设U=1, V = 0则可以得出一组解，设为X=X3, Y=Y3, Z=Z3，但若设U=1, V=1则
13、 方程得出另一组解，设为X = X4, Y=Y4, Z = Z4这两组解之间存在着如下关系,即：X4=Y4 =Y1 + Y32Z4=Z1 + Z 3 2故U=1, V=1的相似准则可以用U=0, V=1和U=1, v = o的相似准则表示，所以说 U=0, V=1和U=1, V = 0的相似准则可以表示U, V为任何实数的相似准则。3. 三种方法解题1）定律分析法已知一个简支梁受有大小为4KN/M均布荷载，简支梁的跨度为4M,截面的高为0.5M,宽为0.4M，跨中截面的最大正应力为480 KN /M 2,求当梁的跨度为2M,截面尺寸相同受均布荷载为2KN/M时的跨中截面的最大正应力。跨中弯矩的
14、公式M=ql 28最大正应力公式6M5 =bh2bh25 解：由最大正应力公式可以推出M =-6ql2又因为M 8所以兀=ql 2bh25ql2 _ q l 2= m m -bh25 b h 25m m m又因为截面的尺寸相同所以可以简化为ql2ql2m m5m所以 5 = $m m 5 =60 KN / M 几何方程mql 22）方程分析法以弹性力学中的极坐标的平面应力问题为例说明1. 写出现象的基本微分方程1）静力学平衡方程d51 比 5 -5d +吩 +4 + f = 0OpP 如 p P1 05氏2t4 + P4+ P4+ f = 0p 040pp 4c -iscc&c -isccc -1 plplUc -iduc -i6pcc-i8c -i E=c -1(8 / +Tpe-c-1pc -18810二兀E8_ E8二兀p81 dup p dedudupep dedp3）物理方程=-(8p8（52(1+ p)pepe4）边界条件（2 个）p = 8 l +t mp pe另外一个类似2. 写出全部的单值条件，并令其二现象相似 1）几何单值条件相似P=c= cP l 式中： c c c 分别表示长度相似常数，应变相似常数和位移相似常数 lU
