边界元法与无网格法无网格法概论优秀PPT.ppt

《边界元法与无网格法无网格法概论优秀PPT.ppt》由会员分享，可在线阅读，更多相关《边界元法与无网格法无网格法概论优秀PPT.ppt（46页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、边界元法与无网格边界元法与无网格法无网格法概论法无网格法概论第一页，本课件共有46页参考文献参考文献|张雄，刘岩.无网格法，清华大学出版社，2004|刘更，刘天祥，谢琴.无网格法及其应用.西北工业大学出版社，2005.|G.R.Liu,Y.T.Gu(王建明、周学军).无网格法理论及程序设计,山东大学出版社，2007.|S.N.Atluri,S.P.Shen.The Meshless Local Petrov-Galerkin Method,Tech Science Press,2002.第二页，本课件共有46页有限元法存在的某些困难有限元法存在的某些困难|冲压成型：网格畸变|裂纹动态扩展：网格

2、重分|高速碰撞：网格畸变|奇异性问题：解析函数|自适应问题：网格重分(h)、近似函数(p)|应变局部化：网格重分|薄壳问题：近似函数高阶连续性问题|复杂三维结构有限元网格的生成第三页，本课件共有46页无网格法无网格法|直接利用分布在求解域中的离散点离散点离散点离散点来构造近似函数的一种求解偏微分方程的数值方法。不需要借助于网格。网格法（有限元法、边界元）无网格法对某些特殊问题，无网格法很有效。第四页，本课件共有46页无网格法概论无网格法概论|无网格法的研究历史无网格法的研究历史|全域插值函数全域插值函数|典型无网格法典型无网格法第五页，本课件共有46页无网格法的研究历史无网格法的研究历史|七十

3、年代：非规则网格有限差分法|1977年：Smoothed particle hydrodynamics SPH归一化光滑函数算法 分片试验不稳定的起因及稳定化方案 克服零能模态的具体方案MLSPH水下爆炸仿真模拟、高速碰撞等第六页，本课件共有46页无网格法的研究历史无网格法的研究历史|1992年：Diffuse element(Nayroles等)|1994年：Element Free Galerkin(Belytschko)动态裂纹扩展数值模拟三维撞击、流体晃动分析 板壳分析EFG和有限元、边界元法耦合相变问题；扩散问题|1995年：Reproducing Kernel Particle M

4、ethod（W K Liu）多尺度分析、自适应分析结构动力学、流体动力学动态断裂和局部化金属加工成形中厚梁板、微电子机械系统纳米管起皱|节理岩体 2000|边界条件 2001|质点积分 2006实质上与EFG等价！第七页，本课件共有46页无网格法的研究历史无网格法的研究历史|1996年：Finite Point Method（Onate等）流体动力学 弹塑性分析弹塑性分析|1996年：Hp Clouds(Oden等)铁摩辛柯梁问题厚板的弯曲问题基于云团法的新型hp有限元Hp无网格云团法|1996年：PUFEM和GFEM（Babuska等）动态裂纹扩展问题|1998年：Local boundar

5、y integral equation method(LBIE)和Meshless local Petrov-Galerkin法(MLPG)(Atluri,Sladek)第八页，本课件共有46页无网格法的研究历史无网格法的研究历史|将无网格法的思想引入有限元法中PUFEM Babuska，1996GFEM DuarteXFEM Belytschko流形元法（石根华）动态裂纹扩展动态裂纹扩展节理岩体节理岩体应变局部化应变局部化网格连续网格连续近似函数不连续近似函数不连续第九页，本课件共有46页无网格法类型无网格法类型|2000年：紧支径向基函数配点法紧支径向基函数配点法|2001年：最小二乘配点

6、无网格法最小二乘配点无网格法|2001年：加权最小二乘无网格法加权最小二乘无网格法|2003年：伽辽金最小二乘无网格法伽辽金最小二乘无网格法|2003年：伽辽金配点无网格法伽辽金配点无网格法|2004年：边界弱形式配点法边界弱形式配点法|2005年：物质点有限元法物质点有限元法|2006年：质点积分无网格伽辽金法质点积分无网格伽辽金法|2009年：冲击爆炸三维物质点法数值仿真软件冲击爆炸三维物质点法数值仿真软件MPM3D第十页，本课件共有46页无网格法活动无网格法活动|1996年 Computer Methods in Engineering Mechanics and Engineering

7、 出版了无网格法专辑（139卷）|2000年 Computational Mechanics 出版了无网格法专辑（25卷，2-3期）|近年来许多著名数值方法国际会议都设置了无网格法的主题会。|许多著名有限元专家，如Zienkiewicz、Belytscho、Atluri、WK Liu、KJ Bathe等都对无网格法进行了深入研究。第十一页，本课件共有46页无网格法概论无网格法概论|无网格法的研究历史无网格法的研究历史|全域插值函数全域插值函数|典型无网格法典型无网格法第十二页，本课件共有46页点插值法点插值法MLS拟合（n m）：I=1,2,n线性函数：二次函数：配点：函数逼近：第十三页，本课

8、件共有46页紧支试函数紧支试函数|移动最小二乘法(Moving Least Square)|核函数近似(Kernel approximation)|重构核近似(Reproducing Kernel approx.)|单位分解法(Partition of Unity)|径向基函数（Radial basis functions）|点插值法（Point interpolation method）|自然邻接点插值|Kriging插值|非均匀有理B样条（NURBS）函数u(x)可以近似为大多数无网格法形函数不满足插值特性，即第十四页，本课件共有46页径向基函数径向基函数|一类以点x到节点xI的距离dI为

9、自变量的函数，也称为距离函数 中心位于节点xI的距离基函数MQ：RMQ：TPS：Gaussians：4th order spline radial basis function：紧支性：紧支性：当当第十五页，本课件共有46页径向基函数径向基函数第十六页，本课件共有46页径向基函数径向基函数|增加多项式基后，MQ插值正定|如果p中包含常数基和线性基，则插值具有一阶一性；|Wang等采用局部形式 径向基点插值法|Hermite型径向基函数插值第十七页，本课件共有46页移动最小二乘近似移动最小二乘近似|近似函数 基函数（多项式或其它已知函数）待定系数线性基：二次基：第十八页，本课件共有46页移动最小

10、二乘近似移动最小二乘近似 待定系数的确定待定系数的确定|有限元法有限元法 令uh(x)在单元节点i处等于函数u(x)在该节点处的函数值ui待定系数的个数必须等于单元自由度uh(xi)=u(xi)具有插值特性依赖于网格|MLS 使uh(x)在节点处的误差在加权最小二乘意义下取极小值精度高，并且可具有高阶连续性能够精确重构基中的任何函数计算量大uh(xi)u(xi)不具有插值特性(拟合)第十九页，本课件共有46页移动最小二乘近似移动最小二乘近似第二十页，本课件共有46页移动最小二乘近似移动最小二乘近似第二十一页，本课件共有46页移动最小二乘近似移动最小二乘近似|当基函数中最高阶完备多项式的阶数k=

11、0时，MLS形函数退化为为Shepard函数|如果在MLS近似中将权函数在域内取为1，在域外取为0，则MLS近似退化为标准的最小二乘近似|MLS近似可以精确重构包含在基底中的任何函数pi(x)，即对于裂纹扩展问题，基函数可以取为对于声场，基函数可以取为第二十二页，本课件共有46页广义移动最小二乘近似广义移动最小二乘近似|MLS要求近似函数在各节点处误差的平方和为最小，但对近似函数导数的误差没有任何约束|Atluri等人在分析欧拉梁时，要求近似函数及其导数在各点处误差的平方和最小可只将函数在面力边界各节点处的导数作为自变量，要求近似函数在所有节点处误差的平方和与近似函数导数在面力边界各点处误差的

12、平方和之和最小第二十三页，本课件共有46页无网格法概论无网格法概论|无网格法的研究历史无网格法的研究历史|全域插值函数全域插值函数|典型无网格法典型无网格法第二十四页，本课件共有46页典型无网格法典型无网格法第二十五页，本课件共有46页典型无网格法典型无网格法|伽辽金型无网格法伽辽金型无网格法|配点型无网格法配点型无网格法配点型无网格法配点型无网格法|基于局部弱形式和边界积分方程的无网格法基于局部弱形式和边界积分方程的无网格法基于局部弱形式和边界积分方程的无网格法基于局部弱形式和边界积分方程的无网格法|最小二乘无网格法最小二乘无网格法|物质点法物质点法|加权余量法加权余量法第二十六页，本课件共

13、有46页伽辽金型无网格法伽辽金型无网格法|EFGMLSGalerkin|RKPMRKGalerkin|PIMPIGalerkin第二十七页，本课件共有46页伽辽金型无网格法伽辽金型无网格法等效积分弱形式（虚位移原理）MLS近似：|计算量大|精度高，稳定性好|需要背景网格进行积分|系数矩阵对称|不易施加本质边界条件处理第二十八页，本课件共有46页伽辽金型无网格法伽辽金型无网格法 数值积分|背景网格积分第二十九页，本课件共有46页伽辽金型无网格法伽辽金型无网格法 数值积分第三十页，本课件共有46页伽辽金型无网格法伽辽金型无网格法 数值积分|节点积分稳定化方案零能模态第三十一页，本课件共有46页伽辽

14、金型无网格法伽辽金型无网格法 位移边界条件的处理位移边界条件的处理|拉格朗日乘子法第三十二页，本课件共有46页伽辽金型无网格法伽辽金型无网格法 位移边界条件的处理位移边界条件的处理|修正变分原理第三十三页，本课件共有46页伽辽金型无网格法伽辽金型无网格法 位移边界条件的处理位移边界条件的处理|罚函数法第三十四页，本课件共有46页加权余量法加权余量法平衡方程应变-位移关系应力-应变关系 边界条件线弹性力学的控制方程线弹性力学的控制方程第三十五页，本课件共有46页加权余量法加权余量法对于复杂问题，只能求近似解加权余量法不要求余量在各点均为零，而要求余量的加权积分为零 平均意义上满足方程近似解：加权

15、余量法的物理意义：选取合适的待定参数强迫余量在某种平均意义下为零等效积分形式检验函数Test functions 试探函数(trial function)待定系数第三十六页，本课件共有46页加权余量法加权余量法如何选取权函数？取：选择不同的权函数，得到不同的加权余量法 Test function第三十七页，本课件共有46页加权余量法加权余量法1.配点法强迫余量在域内N个离散点上为零!2.子域法强迫余量在N个子域VI上的积分为零!第三十八页，本课件共有46页加权余量法加权余量法3.最小二乘法强迫余量的均方和为最小！4.伽辽金法在有限元法中主要采用伽辽金法第三十九页，本课件共有46页加权余量法加权

16、余量法5.局部彼得洛夫-伽辽金法余量在各子域s上消除，且检验函数和试探函数取自不同的函数空间第四十页，本课件共有46页加权余量法加权余量法|子域s以各节点为中心：球形、椭球、立方体在子域中积分，易于实现|检验函数MLS权函数Dirac 函数误差函数微分方程的基本解：等价于LBIE单位阶跃函数：无域内积分，效率高近似函数：与Galerkin法等价|最终得到的求解方程的系数矩阵一般是不对称第四十一页，本课件共有46页无网格法的应用无网格法的应用|天体物理|水下爆炸|高速碰撞|动态裂纹扩展|板壳分析|流体力学|相变问题|大变形|自适应|金属加工成型|塑性变形局部化|纳米管起皱第四十二页，本课件共有4

17、6页岩土数值模拟岩土数值模拟|岩土特性：强度和压力相关性，剪胀性，非关联塑性，以岩土特性：强度和压力相关性，剪胀性，非关联塑性，以及应变软化和等压屈服特性。及应变软化和等压屈服特性。|Drucker-Prager弹塑性本构模型弹塑性本构模型 屈服面与压力相关，剪胀性和非关联塑性。屈服面与压力相关，剪胀性和非关联塑性。屈服面：屈服面：f s和和f t平面上平面上 第四十三页，本课件共有46页岩土数值模拟岩土数值模拟|Drucker-Prager模型模型剪切失效剪切失效-屈服函数屈服函数 拉伸失效拉伸失效-屈服函数屈服函数 屈服面 剪切失效剪切失效-塑性势函数塑性势函数 非关联流动非关联流动拉伸失

18、效拉伸失效-塑性势函数塑性势函数 关联流动关联流动摩擦角：粘聚力：c剪胀角：第四十四页，本课件共有46页岩土数值模拟岩土数值模拟 边坡失效边坡失效|粘土边坡粘土边坡E=70MPa,=0.3,=2.1 g/cm3 =20;=00;c=10.0 kPaMPM 边坡滑移时间边坡滑移时间:5.8s 滑移距离滑移距离:7.5m CPU时间时间:6分钟分钟 LS-DYNA(FEM)滑移时间滑移时间:6.2s 滑移距离滑移距离:6.5m CPU时间时间:262分钟分钟 第四十五页，本课件共有46页岩土数值模拟岩土数值模拟 边坡失效边坡失效|沙土边坡沙土边坡=20;c=0.1 kPa;=00MPM 休休 止止 角角:15.5 滑移时间：滑移时间：8s 计算时间：计算时间：6分钟分钟LS-DYNA(FEM)休休 止止 角角:17.5 滑移时间滑移时间:8.5s 计算时间计算时间:19.8小时小时严重网严重网格畸变格畸变第四十六页，本课件共有46页