公开(公告)号：CN111475914A

公开(公告)日：2020-07-31

申请号：CN202010146716.8

申请日：2020-03-05

Applicant: 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院)

Inventor： 雷卫东 ,  秦晓飞 ,  周海 ,  陈锐

IPC: G06F30/20 ,  G06F119/14

Abstract: 本发明涉及动力学分析领域，特别涉及一种基于时域边界元法分析半无限域双孔洞动力问题的数值方法。包括以下步骤：S1.根据半无限域双孔洞问题建立位移边界积分方程、应力边界积分方程；S2.将所解问题定义域的边界划分成若干个边界元，并引入无穷边界单元，求解出无穷边界单元的插值函数；S3.根据步骤S2划分的边界元及其单元类型，对边界积分方程进行离散；S4.将源点P遍历所有边界节点后，基于算法程序求解出边界积分方程的位移和动力。采用边界元法分析，可以降低求解问题的维度，只需在边界上划分网格，自动满足无穷远处的边界条件，无需施加人工边界。与其他数值方法相比，本发明提供的方法，显著地减少了计算工作量，提高了计算精度和效率。

公开(公告)号：CN111475973A

公开(公告)日：2020-07-31

申请号：CN202010146889.X

申请日：2020-03-05

Applicant: 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院)

Inventor： 雷卫东 ,  秦晓飞 ,  刘志坚 ,  陈锐

IPC: G06F30/23 ,  G06F119/14

Abstract: 本发明涉及一种时域边界元与Newmark精细直接积分法的分域耦合算法。包括以下步骤：S1.求出并已知第n时间步的所有结点位移和结点力；S2.假定通过有限元控制方程求得 再根据Newmark精细直接积分控制方程，由 求得 S3.利用位移协调条件，通过求得 S4.将 代入时域边界元控制方程得到 和S5.通过空间上力转换矩阵，由 求得 S6.根据收敛条件 判断 是否收敛，ε为指定的收敛系数；若满足收敛条件，以 作为公共界面上真实的有限元子域的等效结点力，然后进行下一时间步的计算；若不满足，则执行步骤S7；S7.取重复步骤S2～步骤S6，在公共界面上对未知量进行不断的迭代更新，将收敛后的位移和力作为公共界面上结点的真实位移和力，然后进行下一时间步的计算。

公开(公告)号：CN111159951A

公开(公告)日：2020-05-15

申请号：CN201911408608.7

申请日：2019-12-31

Applicant: 哈尔滨工业大学(深圳)

Inventor： 雷卫东 ,  秦晓飞 ,  段志遥 ,  陈锐

IPC: G06F30/23 ,  G06F30/13

Abstract: 本发明提供了一种基于ABAQUS有限元与边界元的耦合方法,包括：S1、ABAQUS根据输入文件调用ABAQUS/Explicit分析计算生成结果文件Job-1.odb；S2、运行odbHistoryOutput.py把位移和坐标数据写入到TXT文本；S3、读取TXT文本的数据生成边界坐标矩阵和节点位移向量；S4、调用SolveStrategye进行时域边界元分析计算；S5、根据时域边界元分析计算将边界的面力转换成节点力；S6、利用平衡原理f=-Mp生成节点力矩阵并写入TXT文本；S7、运行Inpmodify.py生成下一个时间步耦合的输入文件。数据读写处理速度快，方便快捷。

公开(公告)号：CN107480338A

公开(公告)日：2017-12-15

申请号：CN201710586507.3

申请日：2017-07-18

Applicant: 哈尔滨工业大学深圳研究生院

Inventor： 雷卫东 ,  刘志坚 ,  李宏军 ,  秦晓飞

IPC: G06F17/50

CPC classification number: G06F17/5018 ,  G06F2217/78

Abstract: 本发明提供了一种基于时空积分域处理时域边界元奇异积分的方法，包括以下步骤：首先，通过奇异分离法处理奇异子矩阵元素，采用线性单元进行离散时，将所有变量转变为对的函数，将形函数 、，及线元 用 来表示。本发明的有益效果是：对于奇异积分的奇异部分积分的处理不再采用非奇异积分处理时采用的先积时间后空间的积分顺序，变换积分顺序后重新进行积分，计算量大大减少，积分结果也变得简便很多，可提高时域边界元处理弹性动力学问题的计算精度和效率。

公开(公告)号：CN111159951B

公开(公告)日：2023-08-18

申请号：CN201911408608.7

申请日：2019-12-31

Applicant: 哈尔滨工业大学(深圳)

Inventor： 雷卫东 ,  吴秉臻 ,  秦晓飞 ,  段志遥 ,  陈锐

IPC: G06F30/23 ,  G06F30/13

Abstract: 本发明提供了一种基于ABAQUS有限元与边界元的耦合方法,包括：S1、ABAQUS根据输入文件调用ABAQUS/Explicit分析计算生成结果文件Job‑1.odb；S2、运行odbHistoryOutput.py把位移和坐标数据写入到TXT文本；S3、读取TXT文本的数据生成边界坐标矩阵和节点位移向量；S4、调用SolveStrategye进行时域边界元分析计算；S5、根据时域边界元分析计算将边界的面力转换成节点力；S6、利用平衡原理f=‑Mp生成节点力矩阵并写入TXT文本；S7、运行Inpmodify.py生成下一个时间步耦合的输入文件。数据读写处理速度快，方便快捷。