应变局部化分析的嵌入强间断多尺度有限元法  

EMBEDDED STRONG DISCONTINUITY MODEL BASED MULTISCALE FINITE ELEMENT METHOD FOR STRAIN LOCALIZATION ANALYSIS

作　　者：卢梦凯 张洪武[1] 郑勇刚[1] Lu Mengkai;Zhang Hongwu;Zheng Yonggang(International Research Center for Computational Mechanics,State Key Laboratory of Structural Analysis for Industrial Equipment,Department of Engineering Mechanics,Dalian University of Technology,Dalian 116024,China)

机构地区：[1]大连理工大学,国际计算力学研究中心,工业装备结构分析国家重点实验室,运载工程与力学学部工程力学系,大连116024

出　　处：《力学学报》2017年第3期649-658,共10页Chinese Journal of Theoretical and Applied Mechanics

基　　金：国家自然科学基金(11232003,11672062);中央高校基本科研业务费(DUT14YQ217)资助项目

摘　　要：固体材料的应变局部化行为是导致结构破坏失效的重要因素之一,开展相关数值模拟分析对于结构安全性评估具有重要意义.然而由于材料的非均质和多尺度特性,采用传统数值方法进行求解时通常需要从最小特征尺度离散求解的结构,这将大幅度增加计算规模和成本.针对这一问题,本文提出了一种基于嵌入强间断模型的多尺度有限元方法.该方法从粗细两个尺度离散求解模型,首先在细尺度单元上引入嵌入强间断模型来描述单元间断特性,所附加的跳跃位移自由度则通过凝聚技术进行消除,从而保持细尺度单元刚度阵维度不变.其次,提出了一种增强多节点粗单元技术,其可根据局部化带与粗单元边界相交情况自适应动态地增加粗节点,新构造的增强数值基函数可以捕捉细尺度间断特性,完成物理信息从细单元到粗单元的准确传递以及宏观响应的快速分析;再次,在细尺度解的计算中,将细尺度解分解为降尺度解与单胞局部摄动解,从而消除弹塑性分析时单胞内部的不平衡力.最后,通过两个典型算例分析,并与完全采用细单元的嵌入有限元结果进行对比,验证了所提出算法的正确性与有效性.Strain localization is a common factor that may lead to the failure of solid structure and its numerical analysis becomes an important aspect for the structural safety evaluation.Due to the heterogeneity and multiscale nature,however,traditional numerical methods need to resolve the structure at the fine scale to obtain reasonable results,which increases drastically the computational scale and cost.To solve this problem,an embedded strong discontinuity model based multiscale finite element method is proposed here.In this method,both coarse and fine scale elements are used to represent the structure.The embedded strong discontinuity model is first introduced into the fine element to describe the discontinuity and the corresponding additional displacement jump degree of freedom on the elemental level can be eliminated with the condensation technique,which keeps the dimensions of the sti_ness matrix unchanged.Then,an enhanced multi-node coarse element technique is proposed,which can adaptively insert coarse nodes according to the intersection between the discontinuity line and coarse element boundary and thus guarantees the proper transformation of information between the fine and coarse elements.The problem can then be e_ectively solved on the coarse scale level.Moreover,a solution decomposition technique,in which the fine scale solution is decomposed into the downscaling and local perturbation solutions,is adopted to eliminate the unbalance forces within the unit cell in the elasto-plastic analysis.Finally,two representative examples are presented to demonstrate the accuracy and e_ectiveness of the proposed method through the comparisons with the results of the embedded finite element method.

关 键 词：多尺度有限元法 嵌入强间断模型 应变局部化 增强粗单元 弹塑性分析 

分 类 号：O242.21[理学—计算数学] O344.3[理学—数学]