直接施加本质边界条件的FE/EFG耦合算法及其数值实现  

A FE/EFG coupled algorithm directly specifying essential boundary condition

作　　者：柳军[1] 严波[1,2] 卢岳川 孙英学 

机构地区：[1]重庆大学工程力学系,重庆400044 [2]核反应堆系统设计技术重点实验室,成都610041

出　　处：《计算力学学报》2011年第5期766-772,共7页Chinese Journal of Computational Mechanics

基　　金：核反应堆系统设计技术国家级重点实验室基金(ZDS-A-0908)资助项目

摘　　要：提出一种可以直接施加本质边界条件的有限元与无网格Galerkin(FE/EFG)耦合算法。将问题域分成FE和EFG两种类型的子域,采用转换矩阵耦合两子域的交界面;通过另一转换矩阵将无网格区域本质边界上的名义位移转换成真实位移,从而可在其上直接施加本质边界条件;采用二次转换实现两种转换矩阵之间的协调。提出全域统一采用单元思想和分块矩阵法计算刚度矩阵,提高了程序的通用性和刚度矩阵的计算速度。编制了相应计算程序,通过数值算例验证了方法和数值实施过程的正确性和有效性。A finite element and meshfree Galerkin （FE/EFG） coupled algorithm based on transformation matrices is presented. By means of two transformation matrices, the FE sub-domain and EFG sub-domain can be coupled effectively and the essential boundary condition can be imposed directly on the EFG subdomain boundary. Two transformations are needed to deal with the nodes on the essential boundary and in the influence domain of the nodes on the interface between the two sub-domains. Furthermore,a partitioned-matrix method is proposed to calculate the transferred stiffness matrix which is stored in one-dimension variable bandwidth. The system equation is solved by the LU decomposition method. The cornputer program of the FE/EFG coupling method for three dimensional problems is developed. The efficiency of the algorithm presented in this papre is demonstrated by two numerical examples,a plate with a hole and a plate with a centre-crack under uniform tension.

关 键 词：有限元法 EFG法 耦合算法 本质边界条件 数值实现 

分 类 号：O346.1[理学—固体力学]