基于Fourier谱元法拱坝结构静动力分析探究.docVIP

基于Fourier谱元法拱坝结构静动力分析探究 　　简介：将谱方法和有限元理论结合形成的谱元法，同时具有谱方法任意阶收敛和有限元可操作性的优点。将Fourier积分函数运用于谱元法对拱坝结构计算中，具有快速、高效、精确的求解优势，对提高数值计算精度和谱元法在工程中的实用性，具有较高的参考价值。 　　关键词：拱坝结构静动力分析、谱元法、Fourier积分函数 　　 　　Research on static and dynamic analysis of dam structure 　　based on Fourier spectral element method 　　Tang Bo 　　School of Civil Engineering, Tong ji University, Shang Hai, China, 200092 　　Abstract：The spectral element method originated from the finite element method (FEM) and the spectral method has the advantages of excellent operability and convergence of arbitrary order. The dam structural calculation using spectral element method with Fourier integral function can solve equations fast, efficiently and accurately. Therefore, it has a large value in improving the numerical accuracy and the practicality of spectral method in practical engineering. 　　Keywards：Static and dynamic analysis of dam structure, spectral element method, Fourier integral function 　　 　　1、引言 　　 结构静、动力分析的高速、有效和准确性，一直以来都是研究人员关注的问题，在结构计算理论也不断发展的基础上，新的计算方法也不断被提出，例如数值差分法、有限元方法等。本文在前人的基础上，利用谱方法的优点，对结构各种荷载激励下的静、动力计算分析方法做出了初步探究。 　　2、基本理论 　　 我们知道，一般静、动力问题的基本理论方程为： 　　 几何方程： (1) 　　为应变向量，为微分算子，为形函数，为几何矩阵，为节点位移向量，为拱坝结构位移矩阵。 　　 物理方程：(2) 　　 为结构应力向量，为刚度矩阵。 　　 力平衡方程：(3) 　　 是微分算子，是该点在单位体积的体积力向量，和分别为质量密度和阻尼系数。 　　 力的边界条件：（4-1） 　　为边界上已知弹性体单位面积上作用的面积力量；为边界外法线方向余弦形成的列向量； 　　 几何边界条件：（4-2） 　　 为边界上已知的位移向量。 　　3、基本动力方程 　　 由有限元法可知，由达朗贝尔原??可知结构的基本动力方程为： 　　（5） 　　单元刚度矩阵：（5-1） 　　单元质量矩阵： （5-2） 　　 单元阻尼矩阵：（5-3） 　　 矩阵形式表达为： 　　 　　 上式中，为单元荷载向量。 　　 则对整个结构可知：（6） 　　 分别为结构刚度矩阵，阻尼矩阵，质量矩阵，荷载矩阵。 　　4、利用Fourier谱元法形成计算形函数矩阵 　　 对于任意函数，采用fourier逼近则有： 　　 由最佳平方逼近可知得： 　　 　　 　　 对于在离散点 给定值的函数，可以类似的求出在离散点的离散Fourier系数和，即得到如下表达式： 　　 （）（7） 　　 利用三角变换可得：（8） 　　 其中： 　　 在实际运用中，采用如下点集： ，则代入上式可得： 　　 （9） 　　 其中： 　　 则可以得到： 　　 对于函数导数的计算可以通过对的导数计算，配置点处的导数为： 　　 （10） 　　 对三维函数，在标准的正方体单元内，分别定义三个方向上的节点系，其中： 　　 　　 则有： 　　 则的傅里叶谱展开式为： 　　 （11） 　　 其中：，， 　　 显然： 　　 若共i个点，则（13）式可表示为： 　　 　　 得到位移形函数再代入下面的计算中。 　　5、刚度矩阵、质量（阻尼）矩阵、荷载向量推导