[发明专利]基于GiD-Code_bright软件的批量修改单元材料方法

说明书

本发明公开了基于GiD‑Code_bright软件的批量修改单元材料方法，包括以下步骤：S1：使用GID‑Code_Bright建立几何模型，选定土体材料，设置边界条件和时间步，并完成计算；S2：提取软件数据；S3：读取软件数据；S4：根据提取结果计算变化值；S5：利用源文件生成不土体材料，替换不同单元格内土体材料；S6：重新计算，整理数据结果。本发明的有益效果是：通过GID‑Code_Bright建立几何模型，选定土体材料，通过软件读取数据相比于人工更便捷，减少了人工工作量，而且利用源文件生成不土体材料，替换不同单元格内土体材料，其参数精准度高，减少人为误差，计算结果更接近工程实际。

技术领域

本发明涉及岩土工程、矿山开采、近海工程，水利工程，防灾减灾工程技术领域，特别是基于GiD-Code_bright软件的批量修改单元材料方法。

背景技术

目前在有限元计算中，针对各种类型的岩土体本构，在应力应变响应的过程中土体性质和形状均会发生改变，并且这种改变的程度在各个单元格内并不相同，因此应当根据解决问题的实际需要，更改不同时间步中不同离散单元的本构参数。但是在有限元软件中附着材料过程中都是以面域进行附着，GID软件虽然提供了可以改变单个单元格材料的功能，但是计算改变参数时，提取原有单元格计算参数数据庞大，在对网格材料修改时人工修改工作量大，不贴合实际。经过发明人长期研究，发明了基于GiD-Code_bright软件的批量修改单元材料方法。

发明内容

本发明的目的在于克服现有技术的缺点，提供基于GiD-Code_bright软件的批量修改单元材料方法。

本发明的目的通过以下技术方案来实现：基于GiD-Code_bright软件的批量修改单元材料方法，包括以下步骤：

S1：使用GID-Code_Bright建立几何模型，选定土体材料，设置边界条件和时间步，并完成计算；

S2：提取软件数据；

S3：读取软件数据；

S4：根据提取结果计算变化值；

S5：利用源文件生成不土体材料，替换不同单元格内土体材料；

S6：重新计算，整理数据结果。

优选的，步骤S2中，利用“res_read.m”文件，读取结果文件中名为“工程名字.post.msh”和“工程名字.post.res”的文件，将所有网格、节点坐标以及各节点、高斯点物理量数据存储在“result.mat”子程序库中。

优选的，步骤S3中，利用“res_save_nodes.m、res_save_Gauss.m”提取节点、高斯点至matlab中指定矩阵中。

优选的，结果文件中，位移、水压力的物理信息存储在网格节点上，使用“res_save_nodes.m”程序读取，而其他物理信息储存在高斯点上，使用“res_save_Gauss.m”文件读取。

优选的，步骤S4中，不同离散单元内塑性应变量不同，因此本构数值变化不同，而数据提取后是存储在矩阵中，将矩阵数值带入经验公式中计算参数该变量，由此得到不同离散单元内不同本构数值。

优选的，步骤S5中，由于不同单元内本构参数不同，根据单元数量生成相应数量的材料，使用“replace.m”程序读取“工程文件名.mat”中的本构材料数据，确定复制材料的编号后，将该编号材料的数据按照单元数量，使用“replace.m”程序进行复制，再利用“replace.m”程序将所述步骤S3中的数据替换至新材料中，对新材料进行编号，生成新的“工程文件名.mat”文件，完成新材料的生成。

优选的，步骤S5中，使用“replace.m”程序读取“工程文件名.msh”中的网格信息，再利用其对材料编号进行替换，生成新的“工程文件名.msh”文件完成新材料的替换。