我需要知道abaqus子程序中“积分点”的定义。我是abaqus软件的新手,我在等你的帮助
OP现在问了这个问题已经过去了2.5年,所以我的回答可能更适合那些在这里关注过链接的人,希望能提供一些见解。因为FEM编程是特殊的,0我将尝试回答这个问题,而不是将其标记为偏离主题。无论如何,我的一些答案一般适用于FEM,有些是针对Abaqus的。
快速检查:如果您只是要求用于积分点的(通常或标准)位置的特定数值,则答案取决于它。幸运的是,标准值可广泛用于各种元素(参见下面的资源)。
但是,我假设您要求编写用户元素(UEL)子例程,但还不熟悉元素的表达方式或集成点是什么。
答案:在基于标准位移的FEM中,单个有限元的本构响应通常通过在元素上或元素内的一个或多个点处的数值积分(也称为正交)来获得。这些点的位置和位置取决于元素类型,某些性能权衡等,以及所使用的特定集成技术。我见过的用于连续(固体)有限元的集成技术包括:
根据我的经验,迄今为止的标准实践是使用高斯积分或简化积分方法(通常是高斯积分的变化)。在高斯求积法中,积分点的位置取自被称为高斯点的元素内的特殊(“最佳”)点,已被证明可为给定的计算费用水平提供高水平的可靠精确解 - 至少对于用于许多等参数有限元的典型多项式函数。在某些情况下,已发现其他集成技术具有竞争力1但高斯求积法肯定是黄金标准。还有其他我不熟悉的技术。
实用建议:假设采用等参数公式,在UEL中使用“元素形状函数”和由节点自由度定义的主要场变量(具有固体力学焦点,这些通常是位移)来计算元素应变,应力等,在每个集成点。如果这对您没有意义,请参阅下面的资源。
请注意,如果需要节点(或任何其他点)处的应力,则必须从积分点推断它们,再次使用形状函数,或直接在节点处计算/积分。
建议的资源:请:如果您正在编写用户子程序,您应该已经知道集成点是什么。对不起,但就是这样。在尝试编写UEL之前,您必须至少了解基础知识。
也就是说,我觉得你对FEA / FEM的编程很感兴趣。如果你有动力,但没有大学,你可以参加一两个FEM课程,那么有大量的资源,从大规模开放在线课程(MOOCs)到大量的教科书 - 我通常推荐anything written by Zienkiewicz。对于一个可读但又“可靠”的介绍,强调固体力学,我喜欢库克斯等人的Concepts and Applications of Finite Element Analysis, 4th Edition(又名“库克书”)。祝好运!
0在提出正确的问题之前,您通常需要大量的背景知识。
1 Trefethen,2008,“高斯正交优于Clenshaw-Curtis?”,DOI 10.1137 / 060659831
你的问题不是很清楚。
你的意思是在python环境中?你有shell元素的截面点,这是通过你的shell部分设置的槽厚度。积分点的数量取决于您的元素类型。
您可以在Abaqus脚本手册中找到大量信息。例如
http://www.tu-chemnitz.de/projekt/abq_hilfe/docs/v6.12/books/cmd/default.htm
FEM中的积分点,主要变量得到解决。要时刻铭记在心。在Abaqus的用户子程序中,计算在每个集成点进行。记住这一点并继续前进。如果您不满意,请查看任何FEM书籍,以了解集成点的定义/说明。它不依赖于子程序。
积分点是元素中的节点值之一。例如,八节点C3D8R连续体砖元素具有八个积分点,在砖的每个角落处有一个积分点。
同样在子程序中的其他变量,如状态变量,SVARS,或存储在积分点,因此如果你的元素说4 SVARS你需要跟踪,那么整个8节点元素中将有8 * 4 = 32个SVARS。
我希望这回答了你的问题。