谱元方法课程.pptVIP

*Tele-mail：glqin@mail.xjtu.edu.cn*标准单元类包括了在各求解单元中的常量。类中的数据及方法主要有：单元内网格的划分(含两方向的节点数及总节点数)、标准插值节点、插值函数的导数矩阵及有关标准插值函数的积分。标准单元类定义，classCStdElement{protected://标准单元类的属性网格属性，标准插值节点，插值函数的导数，插值函数积分；public://标准单元类方法标准单元类的构造，Chebyshev多项式及其导数；插值函数的积分，属性设置及取值；}第62页，共94页，星期日，2025年，2月5日*Tele-mail：glqin@mail.xjtu.edu.cn*单元类主要数据及方法包含元素序号、单元内节点划分(和标准单元内划分相对应)、标准单元类、节点类和总体坐标与局部坐标的变换关系。单元类定义，classCElement{protected://单元类的属性单元及网格属性，标准单元类，坐标变换，单元矩阵；public://单元类方法单元构造，标准单元类的构造，节点类构造；元素矩阵的形成，属性设置及取值；}第63页，共94页，星期日，2025年，2月5日*Tele-mail：glqin@mail.xjtu.edu.cn*总体区域类主要数据及方法包含区域内单元数、总节点数、节点类、标准单元类、总体矩阵等。总体区域类定义，classCDomain{protected://整体区域类的属性区域的属性，单元及网格属性，单元类，标准单元类；总体矩阵，对角元位置，边界节点序号；public://整体区域类方法单元构造，标准单元类的构造，节点类构造；元素矩阵的形成，属性设置及取值；}第64页，共94页，星期日，2025年，2月5日*Tele-mail：glqin@mail.xjtu.edu.cn*6.2谱元方法程序设计中的几个问题插值函数积分的存储节点的编码及总体矩阵的形成对称带状矩阵的一维存储元素的编码影响元素和影响节点对角元的位置及一维数组的长度总刚度矩阵第一类边界条件的处理第65页，共94页，星期日，2025年，2月5日*Tele-mail：glqin@mail.xjtu.edu.cn*插值函数积分的存储在形成元素矩阵时，要计算如下积分，可以看出，在(i,j,k)变化时，不必占用(Nx+1)3个存储单元，这样既节约了存储空间又节省了计算时间。Bijk所占的空间可由下式计算，Bijk的检索按如下方式进行，先将i,j,k按由大到小排列，设ijk，则Bijk在数组中的位置为，第66页，共94页，星期日，2025年，2月5日*Tele-mail：glqin@mail.xjtu.edu.cn*3.7偏微分方程元素矩阵的形成具体考察Helmholtz方程，其中，作以下坐标变换在标准正方形单元内，使用伪谱Chebyshev谱逼近u(ξ,η)，则元素内线性方程组式变为，第30页，共94页，星期日，2025年，2月5日*Tele-mail：glqin@mail.xjtu.edu.cn*总体坐标系与局部坐标系第31页，共94页，星期日，2025年，2月5日*Tele-mail：glqin@mail.xjtu.edu.cn*元素矩阵第32页，共94页，星期日，2025年，2月5日*Tele-mail：glqin@mail.xjtu.edu.cn*3.8极坐标系下环形区域的谱元方法在极坐标系下，Helmholtz方程表达为，将极坐标系下环形段单元，变换为直角坐标系下的正方形单元，坐标变换关系如下，元素矩阵的形式为，第33页，共94页，星期日，2025年，2月5日*Tele-mail：glqin@mail.xjtu.edu.cn*极坐标坐标系与局部坐标系第34页，共94页，星期日，2025年，2月5日*Tele-mail：glqin@mail.xjtu.edu.cn*元素矩阵式中各变量的表达式如下，第35页，共9