(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202210897651.X (22)申请日 2022.07.28 (71)申请人 天津城建大 学 地址 300000 天津市西青区津静 路26号 (72)发明人 李杰　苏纪浩　张学杰　穆全平　 李淼　包第啸　 (74)专利代理 机构 天津市君砚知识产权代理有 限公司 12 239 专利代理师 刘雅爽 (51)Int.Cl. G06F 30/23(2020.01) G06F 119/02(2020.01) (54)发明名称 基于LS-DYNA的非均质柱子模型及建立方法 (57)摘要 本发明公开了一种基于LS ‑DYNA的非均质柱 子模型及建立方法， 包括： 获取用户输入的修改 指令并解析； 获得指令类型和对应参数， 利用预 设的随机函数产生随机数， 根据随机数产生的新 数据对原数据进行替换， 形成新的K文件， 将新的 K文件存放在预设的指定位置； 获取LS ‑DYNA求解 器的位置， 计算CPU核数， 以及存放新的K文件的 位置， 生成DOS命令， 进行批量计算； 逐一提取计 算结果中的位移信息和合力信息， 生成爆炸冲击 荷载后的柱状图以及统计数据。 采用随机数对修 改的参数生成新的K文件的形式， 实现对参数的 批量修改， 采用DOS命令形式， 实现批量计算， 极 大的提高了 计算效率， 减轻了 计算负担 。 权利要求书2页 说明书6页 附图5页 CN 115186560 A 2022.10.14 CN 115186560 A 1.一种基于LS ‑DYNA的非均质柱子模型建立方法， 其特 征在于， 包括以下步骤： S1、 获取用户输入的修改指令， 并对修改指令进行解析； S2、 根据解析所述修改指令获得指令类型和对应参数， 利用预设的随机函数产生随机 数， 根据随机数产生的新数据对原数据进行替换， 并将替换后的新数据形成新的K文件， 将 新的K文件 存放在预设的指定位置； S3、 获取LS ‑DYNA求解器 的位置， 根据求解器位置计算CPU核数， 以及存放新的K文件的 位置， 生成DOS命令， 进行批量计算； S4、 逐一提取计算结果中的位移信息和合力信息， 根据用户输入的D值柱状图、 残余承 载力的水平坐标和垂直坐标， 以及坐标分度值， 生成爆炸冲击荷载后的柱状图以及统计数 据， 完成非均质柱子模型的建立。 2.根据权利要求1所述的基于LS ‑DYNA的非均质柱子模型建立方法， 其特征在于， 在S1 中， 所述用户输入的修改指令包括数据类型以及数据符合的分布类型， 所述数据类型包括 混凝土抗压强度、 钢筋形态、 屈服 强度、 直径、 划分块数的长度大小以及柱子模型的长度； 所 述钢筋形态包括纵筋或箍筋。 3.根据权利要求2所述的基于LS ‑DYNA的非均质柱子模型建立方法， 其特征在于， 在S2 中， 所述根据解析所述修改指令获得指令类型和对应参数， 利用预设的随机函数产生 随机 数， 还包括将输入的划分块数 的长度大小以及柱子模型 的长度， 带入空间混凝土抗压强度 分布公式中， 根据混凝土材料型号的不同， 计算得出不同的数据， 并把得出的不同数据带入 到抗压强度随机数产生 函数中。 4.根据权利要求2所述的基于LS ‑DYNA的非均质柱子模型建立方法， 其特征在于， 在S2 中， 所述根据随机数产生的新数据对原数据进行替换， 还包括根据K文件中， 命令行与修改 位置的对应关系， 相邻两参数预留的固定的字符位数， 通过计算要修改的参数所在的字符 位数， 直接对K文件的命令行中， 要修改参数的字符位数进行替换。 5.根据权利要求1所述的基于LS ‑DYNA的非均质柱子模型建立方法， 其特征在于， 在S2 中， 还包括所述将替换后的新数据形成新的K文件时， 将原K文件中未修改数据直接进行复 制， 将替换后的新数据与复制过来的未修改数据进行拼 合， 形成新的K文件。 6.根据权利要求5所述的基于LS ‑DYNA的非均质柱子模型建立方法， 其特征在于， 还包 括对新的K文件按照原文件名_n的命名规则进行命名， 其中， n 为自然数。 7.根据权利要求1所述的基于LS ‑DYNA的非均质柱子模型建立方法， 其特征在于， 在S4 中， 还包括所述逐一提取计算结果中的位移信息和合力信息， 生成位移时程曲线， 只提取时 程曲线中所需要的一段数据， 保存为xls文件， 采用xlrd、 xlwt、 xlutils专门用来处理xls文 件的库， 求 解出数据的均值， 变异系数， 标准差数据。 8.一种基于LS ‑DYNA的非均质柱子模型， 其特征在于， 包括第一数据获取模块、 第一执 行模块、 第二数据获取模块和第二执 行模块； 所述第一数据获取模块， 用于获取用户输入的修改指令， 并对修改指令进行解析； 所述第一执行模块， 用于根据解析所述修改指令获得指令类型和对应参数， 利用预设 的随机函数产生随机数， 根据随机数产生的新数据对原数据进行替换， 并将替换后的新数 据形成新的K文件， 将新的K文件 存放在预设的指定位置； 所述第二数据获取模块， 用于获取LS ‑DYNA求解器的位置， 根据求解器位置计算CPU核权　利　要　求　书 1/2 页 2 CN 115186560 A 2数， 以及存放 新的K文件的位置， 所述第二执行模块； 用于根据CPU核数以及存放新的K文件的位置， 生成DOS命令， 进行 批量计算； 逐一提取计算结果中的位移信息和合力信息， 根据用户输入的D值柱状图、 残余 承载力的水平坐标和垂直坐标， 以及坐标分度值， 生成爆炸冲击荷载后的柱状图以及统计 数据。 9.根据权利要求8所述的基于LS ‑DYNA的非均质柱子模型， 其特征在于， 所述第一数据 获取模块中， 所述用户输入的修改指令包括数据类型以及数据符合的分布类型， 所述数据 类型包括混凝土抗压强度、 钢筋形态、 屈服强度、 直径、 划分块数 的长度大小以及柱子模型 的长度； 所述钢筋形态包括纵筋或箍筋。权　利　要　求　书 2/2 页 3 CN 115186560 A 3