专利名称：一种混凝土裂缝扩展仿真方法
技术领域：
本发明涉及一种混凝土裂缝扩展仿真方法，属于计算机仿真技术领域。
背景技术：
计算机仿真技术是利用计算机对系统的结构、功能和行为等进行动态性逼真模仿 的技术，目前已广泛应用于各个工程技术领域。混凝土裂缝扩展仿真作为计算机仿真技术 在混凝土断裂性能研究中的具体应用，对于建筑、桥梁、水利等的工程设计和施工具有重要
眉、ο现有混凝土裂缝扩展仿真方法大多采用非线性迭代的方法，其一般流程为计算 等效荷载列阵一组装等效刚度矩阵一求解控制方程一计算各积分点应变一计算各积分点 名义应力一修正各积分点当前弹性模量、本构矩阵一计算各积分点当前实际应力一判定各 积分点应力状态一判定开裂点一对所有开裂点进行状态标记、应力释放、弹模修正等一计 算不平衡力一重复以上步骤直至满足收敛标准。该迭代流程在处理混凝土开裂问题时，通 常设定超出抗拉强度的积分点/单元应力不分前后、均标记开裂、进行应力释放的计算策 略，该计算策略与混凝土实际开裂情况不吻合，导致计算结果与实际情况误差较大，更无法 精确估算裂缝出现的位置、裂缝间距以及裂缝宽度等，因此不适用于对裂缝分布规律以及 开裂全过程的精细化数值模拟。此外，目前常用的混凝土开裂模型有分离裂缝模型或弥散 裂缝模型，分离裂缝模型将裂缝设置在单元边界上，一旦出现裂缝就调整模型单元网格，该 模型概念清晰，可实现混凝土结构裂缝分布位置、裂缝宽度以及裂缝开展全过程仿真模拟。 但在实际使用时，裂缝模拟过程复杂繁琐、需要不断调整网格，对计算机技术的要求较高， 因而受到了诸多限制。弥散裂缝模型设定当单元某一特征代表点的应力超过了混凝土材料 的抗拉强度，就认为单元开裂，裂缝分布于整个单元内部，裂缝既连续、微小且是彼此平行 的。由于出现裂缝后只需对材料本构矩阵加以调整，不需要增加节点和重新划分单元，易于 有限元程序实现，但该模型将裂缝弥散至整个单元，同样无法精确模拟混凝土裂缝出现、扩 展全过程，获取结构裂缝宽度。针对目前混凝土裂缝扩展仿真方法中存在的不足，需采用一种改进、优化的混凝 土裂缝扩展模拟方法，以期更精确、合理地进行混凝土开裂精细化数值仿真。
发明内容
本发明所要解决的技术问题在于克服现有技术对所有可能开裂积分点同时进行 应力释放，从而导致的仿真结果与混凝土实际开裂情况误差较大的缺陷，提供一种混凝土 裂缝扩展仿真方法，该方法对弥散裂缝模型进行了修正、完善，通过最可能开裂点的判定， 进行开裂小迭代循环，确定混凝土裂缝出现、发展过程。该方法可精确反映混凝土裂缝出 现、扩展全过程以及裂缝出现的精确范围。本发明采用如下技术方案解决上述技术问题
一种混凝土裂缝扩展仿真方法，包括根据各积分点应力状态判定可能开裂的积分点，并对可能开裂的积分点进行应力释放、刚度修正的过程；该过程按照以下步骤
首先根据各积分点应力状态判定可能开裂的积分点；然后根据客户的关注要点从可能 开裂的积分点中选择最可能开裂点；在保持其余积分点状态不变的同时，对最可能开裂点 进行应力释放、刚度修正；进行位移、应力重分布，并再次判定最可能开裂点，直至没有新开 裂点出现。所述根据各积分点应力状态判定可能开裂的积分点是以各积分点的主拉应力、断 裂能、损伤因子中的一种作为判定指标。本发明优选主拉应力作为判定指标，具体按照以下 方法
判断各积分点的主拉应力是否大于混凝土抗拉强度，如是，则判定该积分点为可能开 裂的积分点。所述选择最可能开裂的积分点可分别或综合参考主拉应力、断裂能、损伤变量以 及材料强度等参数进行选择，本发明具体按照以下方法
按照下式计算各可能开裂的积分点的开裂风险因子，并选择开裂风险因子的值最大的 积分点为最可能开裂点，
权利要求
1.一种混凝土裂缝扩展仿真方法，包括根据各积分点应力状态判定可能开裂的积分 点，并对可能开裂的积分点进行应力释放、刚度修正的过程；其特征在于，该过程按照以下 步骤首先根据各积分点应力状态判定可能开裂的积分点；然后根据客户的关注要点从可能 开裂的积分点中选择最可能开裂点；在保持其余积分点状态不变的同时，对最可能开裂点 进行应力释放、刚度修正；进行位移、应力重分布，并再次判定最可能开裂点，直至没有新开 裂点出现。
2.如权利要求1所述混凝土裂缝扩展仿真方法，其特征在于，所述根据各积分点应力 状态判定可能开裂的积分点是以各积分点的主拉应力、断裂能、损伤因子中的一种作为判 定指标。
3.如权利要求2所述混凝土裂缝扩展仿真方法，其特征在于，所述根据各积分点应力 状态判定可能开裂的积分点具体按照以下方法判断各积分点的主拉应力是否大于混凝土抗拉强度，如是，则判定该积分点为可 能开裂的积分点。
4.如权利要求1所述混凝土裂缝扩展仿真方法，其特征在于，所述选择最可能开裂的 积分点具体按照以下方法按照下式计算各可能开裂的积分点的开裂风险因子，并选择开裂风险因子的值最大的 积分点为最可能开裂点，式中，y为积分点开裂风险因子V力积分点主拉应力，单位为MPa ； J；为积分点抗拉 强度，单位为Mpa。
5.如权利要求1所述混凝土裂缝扩展仿真方法，其特征在于，所述在保持其余积分点 状态不变的同时，对最可能开裂点进行应力释放、刚度修正，具体按照以下步骤 步骤1、查找与最可能开裂点相关的单元，记为Ei (i=l，2……)； 步骤2、读取上一迭代步中E i的弹性模量、泊松比，计算E i单元上迭代步单元刚度矩阵；步骤3、根据E i单元当前的弹性模量、泊松比，计算E i单元本迭代步单元刚度矩阵； 步骤4、得到本迭代步E i单元刚度矩阵的改变量；步骤5、由总体刚度指示矩阵，将步骤中得到的单元刚度矩阵的改变量迭代至总体刚度 矩阵中；步骤6、计算因积分点开裂所产生的应力释放值。
全文摘要
本发明公开了一种混凝土裂缝扩展仿真方法，属于计算机仿真技术领域。本发明包括根据各积分点应力状态判定可能开裂的积分点，并对可能开裂的积分点进行应力释放、刚度修正的过程；该过程按照以下步骤首先根据各积分点应力状态判定可能开裂的积分点；然后根据客户的关注要点从可能开裂的积分点中选择最可能开裂点；在保持其余积分点状态不变的同时，对最可能开裂点进行应力释放、刚度修正；进行位移、应力重分布，并再次判定最可能开裂点，直至没有新开裂点出现。本发明可精确反映混凝土裂缝出现、扩展全过程以及裂缝出现的精确范围，既保留了弥散裂缝模型不需要重新划分网格、计算简便的优点，又克服了弥散裂缝模型与实际开裂情况不符的缺点。
文档编号G06F17/50GK102129491SQ201110047579
公开日2011年7月20日 申请日期2011年2月28日 优先权日2011年2月28日
发明者刘兰珣, 刘婷, 吴胜兴, 王潘绣, 赵海涛, 黄冬辉 申请人:河海大学