渗流模型结果导入结构模型进行水土耦合计算.pdfVIP

渗流模型结果导入结构模型进行水土耦合计算 第一步：在ADINA-Thermal 模块进行渗流计算 01：读入渗流模型命令流文件 生成模型如下： 02：查看渗流结果并导出孔压结果文本文件 查看总水头云图如下： 定义孔压水头参数： 查看孔压水头云图如下： 查看孔压云图如下： 点击如下查看FLUX 参数云图如下： 则浸润面如下图所示： 把孔压大于零的结果提取为文本格式文件： 03：通过Excel 调整孔压结果文本文件 打开fe-pressure.txt 并删除一些不需要的内容： 删除一些不必要的内容删除一些不必要的内容 通过Excel 导入文本数据并进行分列： 把不需要的Node 列删除掉： ADINA-Structure 结构模型所需要的孔压应是负值，故在此需把孔压列乘以-1。 转化前 转化后 在两列中间施加“1”列并存为文本fe-pressure.txt。 第二步：在ADINA-Structures 模块进行结构计算 把渗流模型的网格复制到结构模块，设置材料、边界、荷载等进行结构计算。 初始阶段竖直Z 方向位移云图： 最终阶段竖直Z 方向位移云图： 初始阶段最大主压应力云图： 最终阶段最大主压应力云图： 此时最大主压应力在节点NODE 505 ，最大主压应力值为4020152Pa ，此时应为总 应力。 第三步：导入渗流的孔压进行水土耦合计算 01：在上步模型中读入fe-pressure.txt 作为节点位移荷载 打开单元组的Porous 的选项，以便于随后进行porous 的水土耦合计算。 在节点荷载处读入fe-pressure.txt 文本文件。 此时，模型及其荷载如下显示： 02：修改求解的dat 文件 通过UltraEdit 打开03_水土耦合计算.dat 求解文件。 点击 在 UltraEdit 里面打开列快模式，查找到 C*** PRESCRIBED DISPLACEMENTS ，此时需要把图示中的2 改为8，此操作的目的是把节点位移荷 载改为孔压荷载。 此时可按照上图的格式在Excel 内建立如下的数据，然后把其复制粘贴到UltraEdit 的数据中，这样就可以把其中的“2”换成“8”。 （注意：在选择列的时候不要多选了，不要把空格也选中进去） 读入POR 结果文件如下，可见荷载已经转变为了节点孔压荷载： 竖直Z 方向位移云图如下： 最大主压应力云图如下： 此时最大主压应力在节点NODE 505 ，最大主压应力值为3864664Pa ，此时应为有 效应力。 第四步：验证计算结果 查看渗流计算时节点505 的孔压结果。 可见对于节点NODE 505 有： 4020152Pa-3864664Pa=155488Pa ≈160819Pa，符合总应力-有效应力=孔压的要求。