三角形单元有限元法程序设计.pptxVIP

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;4、引入约束条件;5、线性方程组求解;6、单元应力;例2：简支梁，梁高3m,跨度18m，厚度1m，承受均布荷载10N/m2。已知 按平面应力问题进行计算。;网格划分;考察点y(m);6.2 提高计算精度的方法 ;平均法整理单元应力;（3）网格合理布局 根据应力梯度使网格的布局合理化。即在梯度大的区域网格密些，梯度小的区域应稀些。密、稀网格之间应逐步过渡。 （4）改用高阶单元 受集中力的悬臂梁，采用128个三结点三角形常应变单元，以及3个八结点四边形高阶单元结果。由图可见，采用高阶元的计算精度比常应变元高得多。;带圆孔方板的网格划分;6.3 通用有限分析软件;习 题