首页/文章/ 详情

数值积分:高斯积分点和权重推导

9月前浏览779

一、概述

  1. 适用场景:数值积分计算、有限元分析、工程力学等需高精度积分的场景

  2. 核心范围:以区间 [-1,1](标准区间)的 1-4 阶高斯积分(高斯 - 勒让德积分)为例讲解推导过程,更高阶的积分参数以此类推。

  3. 内容结构:基础理论→分阶推导→结果汇总→区间转换→应用注意事项

二、高斯积分基础理论

1. 基本公式

标准区间 [-1,1] 的高斯积分公式为:

   

  •        :积分阶数(积分点个数)

  •        :积分点([-1,1] 内的坐标,为        次勒让德多项式的根)

  •        :积分权重(非负系数,通过代数精度条件求解)

  • 核心优势:        阶积分可精确计算 2n-1 次多项式(代数精度为 2n-1)

2. 关键支撑:勒让德多项式

(1)定义与表达式

勒让德多项式    是区间 [-1,1] 上的正交多项式,前 5 次表达式为:

  •        (0 次)

  •        (1 次)

  •        (2 次)

  •        (3 次)

  •        (4 次)

(2)核心性质(积分点推导依据)

  1. 根的分布:        的        个根均在开区间 (-1,1) 内,且关于原点对称(        

  2. 正交性:对任意次数        的多项式        ,满足        

三、1-4 阶高斯积分分阶推导

1. 1 阶高斯积分(n=1)

(1)目标

求 1 个积分点    和 1 个权重    ,满足代数精度 1(精确积分 1 次多项式)

(2)步骤 1:求积分点(1 次勒让德多项式的根)

1 次勒让德多项式:    

令    ,解得:    

(3)步骤 2:求权重(代数精度条件)

取    (0 次多项式,必被精确积分),代入积分公式:

   

解得:    

(4)结果

  • 积分点:        

  • 权重:        

  • 代数精度:1

2. 2 阶高斯积分(n=2)

(1)目标

求 2 个积分点    和 2 个权重    ,满足代数精度 3(精确积分 3 次多项式)

(2)步骤 1:求积分点(2 次勒让德多项式的根)

2 次勒让德多项式:    

令    ,即    ,解得:

   

(对称分布,符合勒让德多项式根的性质)

(3)步骤 2:求权重(列方程组求解)

取 2 个低次多项式    、    (次数 < 2,必被精确积分),列方程:

  1.        

  2.        

代入    、    ,化简方程 2:

   

代入方程 1:    

(4)结果

  • 积分点:        ,        

  • 权重:        

  • 代数精度:3

3. 3 阶高斯积分(n=3)

(1)目标

求 3 个积分点    和 3 个权重    ,满足代数精度 5(精确积分 5 次多项式)

(2)步骤 1:求积分点(3 次勒让德多项式的根)

3 次勒让德多项式:    

令    ,因式分解得:    ,解得:

   

(对称分布,含原点)

(3)步骤 2:求权重(列方程组求解)

取 3 个低次多项式    、    、    (次数 < 3),列方程:

  1.        

  2.        

  3.        

利用对称性(    )得    ,且    ,简化方程:

  • 方程 2:        (自动满足)

  • 方程 3:        ,代入        

           

代入方程 1:    

(4)结果

  • 积分点:        ,        ,        

  • 权重:        ,        

  • 代数精度:5

4. 4 阶高斯积分(n=4)

(1)目标

求 4 个积分点    和 4 个权重    ,满足代数精度 7(精确积分 7 次多项式)

(2)步骤 1:求积分点(4 次勒让德多项式的根)

4 次勒让德多项式:    

令    ,设    ,转化为二次方程:    

用求根公式    (    ),解得:

   

因此积分点为:

   

   

(3)步骤 2:求权重(列方程组求解)

取 4 个低次多项式    、    、    、    (次数 < 4),利用对称性(    )得    ,简化方程:

  1.        

  2.        

代入    、    ,解方程:

由方程 1 得    ,代入方程 2:

   

化简计算得:

   

(4)结果

  • 积分点:        ,        ,        ,        

  • 权重:        ,        

  • 代数精度:7

四、1-4 阶高斯积分结果汇总表

阶数 n      
积分点          ([-1,1] 区间,精确值 / 近似值)      
积分权重          (精确值 / 近似值)      
代数精度      
1      
                     
 /                 
1      
2      
             
 /                 
             
 /                 
3      
3      
             
;           /                 
             
 /           ;           /                 
5      
4      
             
 /           、                
             
 /           、                
7      

五、区间转换:从 [-1,1] 到任意区间 [a,b]

1. 转换公式

设原积分区间为 [a,b],通过变量替换    (    ),则    ,积分公式变为:

   

  •        :[-1,1] 区间的高斯积分点(查表获取)

  •        :[-1,1] 区间的高斯权重(查表获取)

2. 示例(2 阶积分,a=0,b=2)

  • 转换关系:        ,        

  • 积分公式:        

  • 代入 2 阶点和权重:        

六、应用注意事项

  1. 精度选择:若被积函数为低次多项式(如≤3 次),2 阶积分即可精确计算;若为高次或复杂函数,需提高阶数(如 4 阶及以上)

  2. 权重性质:所有权重均为正数,避免数值计算中的正负抵消误差

  3. 计算精度:实际应用中积分点和权重建议保留 5-6 位小数,平衡精度与计算量

来源:有限元先生
理论
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-11-13
最近编辑:9月前
外太空土豆儿
博士 我们穷极一生,究竟在追寻什么?
获赞 46粉丝 48文章 117课程 0
点赞
收藏
作者推荐

ABAQUS UEL二次开发帖子大放送

概述  为方便阅读,总结了 公 众号发过的有关于UEL二次开发的帖子,如果觉得有用,欢迎转发、推荐。帖子总结基于多面体比例边界有限元法的ABAQUS用户自定义单元 摘要:采用ABAQUS平台提供的UEL二次开发接口,将比例边界有限元法(The Scaled Boundary Finite Element Method,SBFEM)嵌入到ABAQUS中,开发一种任意多面体自定义单元,该自定义单元没有面和节点数量的限制,相比较经典的结构化有限元网格,该种单元非常适用于复杂几何模型的离散,尤其在涉及到多种尺度的几何模型中,任意多面体单元能够快速的完成尺寸由小到大的过度,降低了有限元模型的自由度。蒙皮单元在UEL二次开发中的作用 摘要:为了用户能够在abaqus中实现自己的有限元理论,abaqus提供了二次开发接口uel,uel子程序设计用来完成一个单元的核心计算,包括单元对运动方程左端贡献:刚度矩阵、质量矩阵和阻尼矩阵,以及对方程右端不平衡力的贡献,这里要注意区分“方程右端量”和“右端不平衡力”的区别。  uel是一款功能强大的子程序,但同时也是一把双刃剑。用户自定义单元无法直接与abaqus自带的功能完美结合,比如给用户自定义单元施加复杂荷载,面力、体力等,也无法给用户自定义单元定义surface,也就无法与abaqus自带的接触非线性模型相结合。  帖子分享了蒙皮单元在uel二次开发中几个常见的作用,包括可视化、施加复杂荷载和定义接触属性,重在讲解思路。基于ABAQUS UEL二次开发的比例边界有限法 摘要:分享一篇关于ABAQUS二次开发的文章,这是比例边界有限法与商业软件结合的第一篇文章,来自于CMAME期刊,很有意义,文章配套的有开源程序,方便学习。CPE8/CPS8用户自定义单元二次开发(1)CPE8/CPS8用户自定义单元二次开发(2) 摘要:有限元法编写二维四边形二阶用户自定义单元,采用abaqus提供的uel二次开发平台,包括单元的刚度矩阵K、质量矩阵M和阻尼矩阵C,uel需要向abaqus主程序输出AMARTIX、RHS等矩阵。    编写的uel子程序适用于静动力计算、频域计算等等。开源线弹性静力小变形UEL子程序分享  摘要:分享一个适用于线弹性静力小变形的UEL子程序,程序来源于github,也可私信我获取全部程序。该UEL可计算的单元类型有:CPE4 C3D8 C3D20 CPE6 CPE8。该程序采用UVARM子程序给用户自定义单元进行可视化。在计算之前首先在自定义单元的初始位置嵌套一层abaqus自带的单元,并且将材料属性设置为无限小,因此嵌套的单元对模型的刚度无贡献,在UEL中计算单元应变和应力,并传递参数到UVARM子程序中。C3D8 BBAR UEL用户自定义单元二次开发(1)C3D8 BBAR UEL用户自定义单元二次开发(2) 摘要:开发了适用于静力通用、频率分析和动力隐式(固定增量步长和自适应增量步长均可)的三维八节点线性UEL,即ABAQUS自带的C3D8单元,该UEL考虑了B-BAR修正,避免剪切锁死。  采用编写的UEL,分别设置了静力通用分析步、频率分析和动力隐式分析步,将计算结果与ABAQUS对比,位移、速度和加速度与ABAQUS均保持一致,说明该UEL复现了一小部分C3D8单元的计算功能。C3D20 UEL用户自定义单元开发(1)C3D20 UEL用户自定义单元开发(2) 摘要:开发了适用于静力通用计算的三维二十节点(C3D20)的用户自定义单元,在挖孔悬臂梁受剪切荷载算例中,位移计算结果与ABAQUS自带单元保持一致。对比刚度矩阵,与abaqus保持一致。C3D8/C3D20 UEL在KOYNA混凝土坝抗震计算中的应用 摘要:以Koyna混凝土坝为对象进行地震响应计算。将自编的八节点UEL和二十节点UEL应用到计算中。分别进行了混凝土坝模态计算和地震时程计算。其中,在模态计算中共设置四种计算工况,分别为:ABAQUS-C3D8、UEL-C3D8、ABAQUS-C3D20、UEL-C3D20。在地震时程计算中设置两种计算工况,分别为:ABAQUS-C3D8、UEL-C3D8。计算结果表明,自编UEL与ABAQUS自带单元结果一致。一种新思路用于实现ABAQUS用户自定义单元 摘要:一种新思路用于实现ABAQUS自定义单元,不需要在INP文件中采用“USER ELEMENT”等关键字定义单元,更重要的是,也不需要在UEL接口中进行繁琐的FORTRAN编程,只需要在INP文件中直接导入相应的矩阵即可。目前这种自定义单元适用的计算有:static, frequency extraction, modal dynamic, mode-based steady-state dynamics, complex eigenvalue extraction, and subspace-based steady-state dynamics。更多的功能还在探索之中。  这种自定义单元方式为多种CAE软件协同二次开发提供了一种可能,众所周知,ABAQUS的非线性方程组求解能力是行业翘楚,这时候如果能得到描述方程组的关键矩阵,可直接导入。 如固体力学的运动方程中的KK、MM和CC,直接导入ABAQUS便可以求解,相较于UEL二次开发,该方法节省了向ABAQUS主程序输出关键矩阵的编程工作,至于 KK、MM和CC,可以自己采用高级语言编程生成,也可以从其他软件中导出。再比如,COMSOL以其多场耦合计算功能著称,采用该方法就可以将COMSOl的多场耦合功能与ABAQUS结合, 本质是数学中的方程组求解。  这种技术可用于解决用户自定义单元在ABAUS中无法可视化的问题 。用户自定义单元的可视化,要么是将数据导出到第三方软件进行处理,要么是采用UMAT套一层单元进行可视化,这两种方法都需要大量的编程工作,采用该技术可避免上述工作。UEL如何与接触非线性模型结合 摘要:分别采用经典有限元法和比例边界有限法在abaqus中编制了C3D8单元的UEL程序,计算结果均与abaqus保持同样精度。  在此基础上将用户自定义单元与abaqus自带的接触非线性模型相结合,并完成了单元的位移、速度和加速度的可视化。来源:有限元先生

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈