各向异性材料中的条单元法（幻灯片）精选.pptVIP

用条元法解决各向异性固体材料中的压力波问题 1.1 简介 1.2 系统方程 1.3 静应力问题的条元法 1.4 动态问题的条元法 1.5 备注 1.1 简介 我们要解决的问题是有压力波在各向异性的薄板中的扩散问题，对于这类问题来说，要想得到一个解析解是很困难的，而诸如有限元法，边界元法等数值解法却可以用来解决此类问题。 目前来说，有限元法已经又了很大发展，但是，由于对单元划分的细致性，它的存储和计算都需要大量的空间。边界元法需要的存储空间相对较少，但在各向异性和不均一的固体材料中使用边界元法却是很困难的，因为对这些材料来说，我们很难找到简单的格林函数。于是就有了条元法，从存储需求上，它具有明显的优越性。 1.2 系统方程 我们只考虑平面上的问题，先定义一个区域Df： x=±a, -∞≤y≤+∞ 它包括带有边界s1,s2,s3,s4的区域Dp。 同时，我们假定材料的各向异性被定义在x-y平面内，于是相应的位移为u,v。 1.2.1 应变-位移关系 其中 是应变矢量 U是位移矢量 特别符号L为 1.2.2 应力-应变关系 应力和应变的关系可作如下表示: 其中 c是各向异性材料的非主轴刚度系数矩阵。 1.2.3 变换方程式 对于不受外力作用的各向异性材料来说，其系统方程如下： 将前述的应力与位移的关系代入，可得 其中 于是作用在平面y方向的压力可表达为： 作用在平面x方向的压力可表达为： 1.2.4 条元法方程 将Df区域划分为N个无限条带单元，对任一单元可作如下假定： 其中N(x)是单元的形态方程矩阵，V(y)是单元的表面振幅的位移矢量。 同时我们假定表面谐波的势能为 是振幅矢量，有 将Df区域的所有单元累加，得到Df区域的近似方程 假定： 代入得， 令Pt=0，可以得到关于k的特征值方程 解特征方程，我们可以得到2M个特征值kj(j=1,2, ‥,2M)，将特征向量 叠加，转移方程可写为 通过在边界S3，S4上的节点，常量C可以用位移矢量表示 通过前述的应力矢量方程，我们可以得到单元内任一点的应力： 于是 将所有条元矩阵的R1，R2叠加，就可以得到所有节点线的应力矢量 通过上式及前述的位移矢量方程，最终可得 其中K是刚度矩阵。上式给出了边界S3,S4上的拉力和位移的关系，但对于 给定的Rbt，通过上式我们就能解得Vbt。 1.3 条元法解决静应力问题 在静应力问题中，时间和频率为0。在半无限区域中，应力问题从起点处开始。当x≥a时，由于对称性，我们只考虑1/4区域，其边界条件为 x=0 y=0 于是线性载荷就可以用三角分布来近似，忽略掉时间部分，其静应力为 1.4 动态问题的条元法 对于无界部分的动态问题，我们需要用无反射的边界来减少区域的大小。对于弹性波，有三个无反射的边界，其条件的一般形式为 第一种情况 第二种情况 第三种情况 在条元法中，上述三种无反射边界条件都能生效。但是，研究发现，上述三种无反射边界条件在半空间区域问题中的运行结果都不是十分令人满意。于是，一种新的粘弹性无反射边界被构思出来。方法是： 在x=a处附加少量假想的粘弹性单元，逐渐增加单元在x方向上的弹性，这些粘弹性单元能有效地吸收来自任一角度的外来波。对对k个粘弹性单元，其材料常量矩阵为： 其中C是弹性常量矩阵。在下面的计算中，α0，β分别选作0.01和1.624。在条元法模型中，在0≤x≤15λ0中共有120个单元，其中100个是弹性单元，20 个是附加在右表面的粘弹性单元。 1.4.1 2D空间的谐波问题 应用条元法解决无边界2D表面的谐波线性载荷问题，仅考虑其1/4区域， 因此，线性载荷可用三角应力分布来近似，而时间部分可以忽略。 1.4.2 半区域的谐波问题 对于在各向异性的半空间材料的普通时域线性载荷来说，没有完全形式的精确弹性力学解法。通过傅里叶变换可以得到积分形式的解： 其中 kt,kl是横向和纵向波的波数，在靠近两极时ζ=±k