基于遗传算法的某带缝重力坝性模量反分析02.docVIP

基于遗传算法的某带缝重力坝性模量反分析02.doc

  1. 1、本文档共6页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
基于遗传算法的某带缝重力坝性模量反分析02

基于遗传算法的某带缝重力坝弹性模量反分析 摘 要:根据某混凝土重力坝坝顶水平位移观测资料,采用逐步回归方法分离出坝顶水平位移水压分量,以大坝和基岩的弹性模量识别反问题作为优化问题,利用改进的遗传算法,并通过采用接触单元模拟坝体纵缝进行有限元反演分析。数值计算结果表明,本文计算方法能够较好的反映坝体结构的实际状况,参数反分析结果具有较高程度的合理性与可靠性。 关键词:重力坝;纵缝;接触单元;优化问题;遗传算法;反分析 1 引 言 对于水工结构实际工程,由于大坝及基岩工作条件复杂,在设计和复核计算中,对荷载、计算参数、边界条件等因素还难于进行精确模拟,这使得目前水工结构计算还难于做到完全与工程实际相吻合,有时甚至会有较大差别。因此,利用原位监测资料反分析大坝及基岩的实际运行状况,从很大程度上可以用来评价大坝和基岩的实际安全度;从而,一方面可对病险坝的处理提供科学依据,另一方面使得大坝在确保工程安全的前提下,充分发挥工程的效益[1]。 就混凝土重力坝而言,施工中一般沿坝轴线设置一条或多条纵向施工缝,以便于施工期混凝土浇注施工和满足温控设计要求,同时能降低运行期坝体所承受的地基不均匀沉陷及环境变化导致的坝体温度变形的影响。然而,由于坝体纵缝的存在,使得坝体变成一个不连续体系统,为解决这样一个复杂系统的接触问题,计算中应该采用合理有效的分析方法或计算理论,使得数值计算能够合理反映工程实际状况。 因此,本文针对坝体纵缝采用接触单元进行二维有限元数值模拟,提出利用坝顶水平位移的原型监测资料进行重力坝混凝土和岩基弹性模量的反演分析,其中将反分析作为优化问题,通过采用改进的遗传算法来进行求解。 2 接触单元及算法 目前,对于求解接触问题,一般采用基于连续介质力学有限元理论所建立的特殊界面单元模型如联结单元模型、节理单元模型和薄层单元模型来模拟,另外就是采用基于非连续介质理论或块体理论所建立的离散单元法和非连续变形分析方法[2]。 这两类方法属于惩罚弹簧法(罚函数法),由于弹簧刚度系数难于确定,因而不能得到结构比较真实的变形、应力及接触特性。而拉格朗日乘子方法和改进的拉格朗日方法可以隐式地得到结构的接触应力。本文在连续介质有限单元法中,采用接触对单元并结合改进的拉格朗日方法对大坝的纵缝进行较为合理的模拟。 2.1 接触模式 相互接触的两部分实体,在接触面的两侧一面定义为目标面单元,另一面定义为接触面单元,它们为无厚度的两层单元(亦称为接触对),覆盖在接缝两侧的实体单元上(如图1-a)。如果不考虑接触限制(即不施加接触协调条件),受压时将必然会导致互相侵入或穿透(如图1-b),这与实际情况不符。因此,必须控制接缝的开度。数学上,施加无穿透的约束是在计算过程中追踪物体的运动轨迹,一旦探测出发生接触(接触面和目标面发生接触),便将接触所需的运动约束(即法向无相对运动、切向可滑动)和节点力(即法向压力和切向摩擦力)作为边界条件直接施加在产生接触的节点上[3]。 但由于数值计算中,要精确描述节点恰好在一个接触断面上很困难。为此,可以采用两种方法来定义接触协调条件:一、罚函数法(如图1-c),用一个弹簧施加接触协调条件,弹簧刚度(或接触刚度)称为罚参数;二、拉格朗日乘子法(如图1-d),增加一个附加自由度(接触压力),以满足不侵入条件。 而将上述两种方法结合起来施加接触协调条件,并采用许可侵入量来进行描述,称为改进的拉格朗日方法(如图1-e)。 a)接触对表示 b)无接触协调条件 c)罚函数法接触协调条件 d)拉格朗日乘子法接触协调条件 e)改进的拉格朗日法接触协调条件 图1 接触模型状态示意图 在本文的有限元计算中采用改进的拉格朗日方法,接触追踪向量定义为 (1) 其中:为接触压力;τ为接触面剪切应力。 2.2 接触压力定义 接触压力的计算考虑如下情况: (2) 其中,为法向接触刚度;为第i次迭代的拉格朗日乘子,可由下式进行定义: (3) 其中,为容差,用来定义接触许可侵入量。 2.3 摩擦定义 图2 摩擦模型 本文计算采用简单库仑粘结-滑移摩擦模型对摩擦进行模拟,计算式表示为: (4) (5

文档评论(0)

qiwqpu54 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档