程序设计 有限元课件技术方案.pptVIP

第2章 有限元程序设计方法;2、单元分析 ;总刚存贮;等带宽形式;方阵存贮和半带宽存贮地址关系;对 称;; 变带宽存贮：按列存贮方式。从左到右，逐列存放；对每一列，先存主对角线元素，然后由下而上顺序存放，直到顶线下第一个元素为止。为避免混淆，我们把存贮[K]的一维数组称为[A]。 实现变带宽存贮的关键问题是：总刚中元素Kij在一维数组A中的地址是什么？为此，需要知道主元Kii在A中的位置和相应列高hi。 主元位置：采用一个一维数组MAXA存主元在A中位置。 MAXA =[1,2,4,6,10,12,16,18，22]。; 从第j列的主对角线元素起到该列上方第一个非零元素为止，所含元素的个数称为第j列的列高，记为hj ；如果把第j列上方第1个非零元素的行号记为mj，则第j列的列高为 hj = j - mj + 1 其实，hj就是第j行的左带宽，因而必有 UBW= max(hj) j=1,2, …,N;1; 连接数组： 1号单元： LM=[0, 0, 1, 0, 0, 2] 2号单元： LM=[0, 0, 2, 3, 4, 5] 3号单元： LM=[3, 4, 5, 6, 7, 8] 4单元： LM=[0, 0, 1, 6, 7, 8]; a) 如果 ID（i, j ) = 0 则表明j号节点第i个自由度受有约束。;主元在一维数组[A]中的地址; MAXA(j) = h1 +…+ hj-2 - hj-1 + 1 =（h1 +…+ hj-2 + 1）+ hj-1 = MAXA（j - 1）+ hj-1 因为永远有 MAXA(1)= 1， MAXA(2)= 2 故计算主元地址的公式可写为 MAXA(j+1)= MAXA(j)+ hj （2-1） 式中， j = 2，3，…，N； hj——刚度矩阵[K]第j列的列高。;Ki,j在一维数组[A]中的地址 ;4、引入约束条件;5、线性方程组求解;6、单元应力;例2：简支梁，梁高3m,跨度18m，厚度1m，承受均布荷载10N/m2。已知 按平面应力问题进行计算。;网格划分;考察点y(m);2.2 提高计算精度的方法 ;平均法整理单元应力;（3）网格合理布局 根据应力梯度使网格的布局合理化。即在梯度大的区域网格密些，梯度小的区域应稀些。密、稀网格之间应逐步过渡。 （4）改用高阶单元 受集中力的悬臂梁，采用128个三结点三角形常应变单元，以及3个八结点四边形高阶单元结果。由图可见，采用高阶元的计算精度比常应变元高得多。;带圆孔方板的网格划分;;2.3 通用有限分析软件;习 题