ABAQUS应用于土石坝数值计算的技巧总结.doc

ABAQUS应用于土石坝数值计算的技巧总结 作者：dailinghui ABAQUS软件是国际著名的有限元通用软件，功能非常强大，但是在进行土石坝应力应变计算这一专业而又复杂的问题上，还存在一些局限性。因此，必须对软件和计算问题有深入的了解，才能找到解决办法，突破局限性。为了便于沟通和交流，本文对ABAQUS软件进行土石坝计算时的一些经验技巧加以总结和探讨。共分为初始建模、网格划分、地应力平衡、无限元节点调整、增量步设置、计算过程控制、结果后处理几个部分进行总结。 一、初始建模 1、ABAQUS CAE建模：ABAQUS软件的前处理功能比较强大，可以通过人性化的界面来建立几何模型。但是广大工程人士最为了解和熟悉的画图软件当属AutoCAD。ABAQUS软件有专门的接口，可以将格式为dxf的CAD文件导入进来。在导入的过程中，可以采用以下步骤来进行：（1）首先将原始图形（坝体剖面图）导到Skectch模块；（2）在Part模块中建立Part（坝体）的时候，由于ABAQUS只允许外边框存在，而不能有内部线，因此，可以先把Skectch模块中的坝体剖面图全部添加进来，只保留坝体的外边框，将其余的线条都删去；（3）对Part进行Partition分割操作：把坝体剖面图全部添加进来，只保留分割坝体所需要用到的线条，把其他多余线条全部删去，包括坝体的外边框。建第二个、第三个Part的时候，重复以上操作即可。采用上述方法建模的好处是，整个建模过程用到的所有线条，都是基于原始图形（坝体剖面图）的，因此可以避免线条有交叉、结点不重合等情况的发生。而且在Assembly模块里将各个Part进行组装的时候，不但不需要再调整各个Part位置，还能保证各个Part之间接触面完全重合，为建立正确的接触关系打下基础。 2、其他前处理软件建模：ABAQUS软件的前处理功能是基于“自上而下”的思路设计的，即先建出模型的整体，然后根据不同部位的要求切出各个小块。而其他有限元软件中还有另外一种建模思路——“自下而上”的思路，即首先按不同部位的要求建出各个小块，然后再组合成整体。专业的前处理软件有很多，比较知名的有Patran软件、Hypermesh软件等。Patran软件是基于“自下而上”的建模思路设计的，它可以为ABAQUS做专业的前处理，最终生成ABAQUS的输入文件，即inp文件。用Patran为ABAQUS作前处理时，建议在进行Equivalence操作（即将各个小块组合成整体）之后，对节点和单元进行一次Renumber操作，即对单元和节点重新编号。如果不进行这项操作，也不会对计算和最终结果产生影响，但是在对结果进行后处理时，将会遇到麻烦。例如用Tecplot软件（专业的后处理软件）进行后处理时，由于节点号和单元号不连续，将会导致无法单独挑出关心部位，如心墙。而进行了Renumber操作之后，就能避免这种情况的发生。 二、网格划分 在对沥青混凝土心墙坝划分网格时，由于心墙与整个坝体相比，厚度方向的尺寸显得非常小，但不能据此在心墙厚度方向上只划分一层网格。因为一层网格会导致整个心墙的应力应变全都受到接触面的影响，从而产生较大误差。如果采用其他后处理软件（如Tecplot）对ABAQUS的结果进行后处理时，两种材料交界面上的节点应力值会被强行平均，这种情况下，如果心墙厚度方向上只有一层网格，内部进行插值后，最终得到的整个心墙的应力值就是完全错误的。所以心墙厚度方向上至少应划分2层～4层网格。 三、地应力平衡 有深厚覆盖层存在时，坝基的地应力平衡就是必不可少的。进行地应力平衡时，首先只对坝基施加重力荷载，然后从结果文件中读出坝基各个单元的应力分量，将其保存到一个文本文档里。然后借助Excel软件读取应力文档进行格式调整和处理，将其转换为ABAQUS可以识别的格式。在这个过程需要注意一个问题，即Excel软件读取应力文档之后，默认的科学记数法的小数位数是2位，而ABAQUS的数据是采用6位有效数字的，即用科学记数法表达时，有5位小数。因此，在Excel里需要将小数调成5位，这样在进行地应力平衡时可以减小误差，得到更为精确的结果。 一般情况下，如果坝基的材料采用弹性模型或者摩尔库仑模型时，进行一次地应力平衡的结果可以使位移降至10-17m的数量级，效果非常好。但是坝基材料采用邓肯张模型时，一次地应力平衡的结果一般只能使位移降至10-4m的数量级，而且位移的最大值可能达到10-3m的数量级。这种情况下，如果对结果要求严格，可重复进行一次或两次地应力平衡，一般能使位移再降低1～2个数量级，而再进行更多次的地应力平衡，就不会有明显的效果了。 四、无限元节点调整 无限元在ABAQUS软件中的实现方法比较简单，虽然无法在ABAQUS CAE里直接划分，但可以通过在inp文件里