第八章 复合材料结构分析有限元法.ppt

符号规定 四节点 ( hg=1) 八节点( hg=2) 参 考 资 料 高阶理论有限元 付晓华，陈浩然. 复合材料多层厚板精化高阶理论及其有限元法，复合材料学报，1992，vol. No. 2； 陈浩然，温玄玲，含分层损伤轴对称层合壳的高阶理论有限元分析，大连理工大学学报，vol. 39，No. 5. P601. 张依芬，温玄玲，陈浩然，含分层损伤复合材料层合板的精化有限元法，大连理工大学学报，vol. 37，NO. 4 P398，1997. Global-Local有限元法 Shi yibing, Chen Haoran. Mixed finite element for thicker laminate shell analogies of composite laminates. Composite structures. Vol.20. No.3.1992 Hamilton理论为基础有限元法 Chen Haoran, Yang Zhenglin etc. process-induced stress analogical composite laminates using semi-analytical Hamilton method. Composite Structures, Vol.41, No. 9. P49-55. (998) 拟协调元层合板有限元法 陈浩然，衣翃， 复合材料层合板拟协调单元，复合材料学报，Vol. 4，No.1 1987 陈浩然，衣翃，复合材料层合板壳通用单元和等网格加筋板壳的局部稳定性，航空学报，vol. 9, NO.3 Pa177 §7.4 考虑剪切效应的分项插值单元 （1）内力势能→应变能→刚度阵 （43） 若板是薄板或忽略剪切效应，则 阵必须是奇异阵，即剪切能量为零。而面内，弯矩耦合刚度阵为非奇异阵，则总单元刚度阵也必须是非奇异阵。 设 为广义位移, 在单元内位移采用分项插值，且等参插值。 （44） 其中 为 i 节点位移向量，即 设应变位移阵为 （45） 其中 将（45）代入（43），并对节点位移向量 变分得 (46) 其中 （2）刚度阵奇异性的判别 a、秩的判别规则 在刚度阵求解时必须采用数值积分，一般采用高斯积分方法。然而采用不同的高斯积分方式所得到的刚度阵，是具有不同性质的特征。 考虑剪切效应的板单元，要求剪切刚度阵为奇异，而其它刚度阵必须非奇异，总单元刚度阵必须非奇异。 判别刚度阵奇异与非奇异，首先可判断它是缺秩或满秩。若缺秩为奇异。令K为n阶方阵，其满秩，则其秩＞n。 ① 矩阵相乘的秩 若 则秩 ② 矩阵相加的秩 若 则秩 b、单刚的秩 其中 为权系数, D为弹性阵，其d×d，秩为d。 2D d = 3 3D d = 6 轴对称 d = 5 B 应变—节点位移向量，其尺寸 , 为节点未知数（自由度） 故 秩为d 为高斯积分阶数。 故 , 秩为 若结构由M个单元组成则总刚度阵的秩为: （47） 为总刚度非奇异性判别准则。（M为单元总数） （3）数例 以四与八节点平面等参元为例 d=3 每节点两个自由度 (黑字为非奇异, 红字为奇异) 192 非奇异 2 16 [(5×9)+(4×5)] ×2-34=96 48 非奇异 1 16 25×2-18=32 c 24 非奇异 2 2 13×2-3=23 6 奇异 1 2 6×2-3=9 b 12 奇异 2 1 8×2-3=13 3 奇异 1 1 4×2-3=5 a M n M n 二次单元（八节点） 线性单元（四节点） (a) (b) (c) 奇异 奇异 奇异 非奇异 非奇异 非奇异 （4）简易刚度阵秩的判别准则 构造考虑剪切效应板单元，若仅判别单刚奇异性，可采用如下简易判别条件，即： 奇异 （48） J 为增加一个单元形成的自由度的增加 d 为弹性阵的秩 为高斯积分阶数 其假定原结构单元总刚是满秩，在此基础上再增加一单元，去判别此单元增加后引起秩的变化。 该方法计算