申请/专利权人：北京大学;北京大学重庆大数据研究院

申请日：2024-02-06

公开（公告）日：2024-03-19

公开（公告）号：CN117725803A

主分类号：G06F30/23

分类号：G06F30/23;G06F30/17;G06T17/20;G06F111/04;G06F113/24;G06F119/14

优先权：

专利状态码：在审-实质审查的生效

法律状态：2024.04.05#实质审查的生效;2024.03.19#公开

摘要：本申请公开了一种基于混合有限元空间求解板弯曲问题的方法。包括：按照预设的处理方式对不同类型的边界条件进行处理，得到载荷向量；将用于求解二维线弹性问题的内蕴混合有限元空间基函数在积分点处的值按照拉格朗日插值点的类型进行分类并组合，得到用于求解板弯曲问题的内蕴混合有限元空间基函数在积分点处的值；基于用于求解板弯曲问题的内蕴混合有限元空间基函数在各个积分点处的值构造系数矩阵；基于系数矩阵和载荷向量构造线性代数方程组并求解，直接得到弯矩与位移结果。本申请解决了相关技术在分析薄板弯曲问题时，采用传统有限元法求解对应的偏微分方程无法在网格数量较少的情况下提升计算精度的技术问题。

主权项：1.一种基于混合有限元空间求解板弯曲问题的方法，其特征在于，包括如下求解步骤：步骤S1，获取单元数据、材料参数、边界条件，其中，单元数据包括：单元形状、节点坐标、单元数量、单元维度，边界条件包括：应力约束、弯矩约束、位移约束、转角约束，材料参数包括：板的厚度、杨氏模量、泊松比、密度；步骤S2，确定边界条件的类型，并按照预设的处理方式对不同类型的边界条件进行处理，得到载荷向量b；步骤S3，获取当前单元的积分点、积分权重、单元面积，依据积分点计算k阶拉格朗日基函数在积分点处的值和k-2阶拉格朗日基函数在积分点处的值，其中，k为大于等于3的整数；步骤S4，获取当前单元的二维对称矩阵空间的张量基，并结合步骤S3得到的k阶拉格朗日基函数在积分点处的值，得到用于求解二维线弹性问题的内蕴混合有限元空间基函数在积分点处的值；步骤S5，将步骤S4得到的用于求解二维线弹性问题的内蕴混合有限元空间基函数在积分点处的值按照拉格朗日插值点的类型进行分类，将分类后的用于求解二维线弹性问题的内蕴混合有限元空间基函数在积分点处的值按照一定条件进行组合，得到用于求解板弯曲问题的内蕴混合有限元空间基函数在积分点处的值；步骤S6，确定当前单元的用于求解板弯曲问题的内蕴混合有限元空间基函数的整体自由度和k-2阶拉格朗日基函数的整体自由度；步骤S7，利用步骤S5所得的用于求解板弯曲问题的内蕴混合有限元空间基函数在各个积分点处的值计算单元能量矩阵，并结合步骤S6所得的用于求解板弯曲问题的内蕴混合有限元空间基函数的整体自由度，构造系数矩阵A；步骤S8，计算当前单元的k阶拉格朗日基函数在积分点处的二阶导数，并与步骤S4所得的当前单元的二维对称矩阵空间的张量基相乘，得到用于求解二维线弹性问题的内蕴混合有限元空间基函数在积分点处的divDiv值；步骤S9，将步骤S8得到的用于求解二维线弹性问题的内蕴混合有限元空间基函数在积分点处的divDiv值按照一定条件进行组合，得到用于求解板弯曲问题的内蕴混合有限元空间基函数在积分点处的divDiv值；步骤S10，利用步骤S3得到的k-2阶拉格朗日基函数在积分点处的值与步骤S9得到的用于求解板弯曲问题的内蕴混合有限元空间基函数在积分点处的divDiv值计算单元divDiv算子矩阵，并将其追加至系数矩阵A中；步骤S11，对循环步骤S3至步骤S10所得的系数矩阵A和步骤S2所得的载荷向量b进行处理，得到最终的线性代数方程组Ax=b；步骤S12，求解线性代数方程组Ax=b，直接得到弯矩与位移结果。

全文数据：

权利要求：

百度查询： 北京大学;北京大学重庆大数据研究院 基于混合有限元空间求解板弯曲问题的方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD