混凝土结构非线性分析课程报告

混凝土结构非线性分析课程报告目 录1、结构非线性分析简介111 结构线性分析与非线性分析的区别112 非线性行为的原因113 非线性结构有限元分析中应注意的问题31. 4 钢筋混凝土结构非线性分析的意义42、混凝土结构非线性相关研究52. 1基于ABAQUS 纤维梁单元的钢筋混凝土柱受力破坏全过程数值模拟53、预应力混凝土结构非线性相关研究831 预应力混凝土结构非线性有限元分析83. 2预应力混凝土结构组合式非线性分析模型114、桥梁结构非线性相关研究1341 钢筋混凝土纤维梁柱单元实用模拟平台1342 预应力混凝土薄壁高墩刚构桥梁极限承载力分析165、总结与展望18参考文献201、结构非线性分析简介11 结构线性分析与非线性分析的区别线性分析在结构方面就是指应力应变曲线刚开始的弹性部分，也就是没有达到应力屈服点的结构分析 非线性分析包括状态非线性，几何非线性，以及材料非线性结构线性分析与非线性分析的区别类型线性分析非线性分析结构属性结构属性（刚度、阻尼等）在分析中是恒定的结构属性随时间、变形和荷载而变化。并与用户定义的属性、荷载及指定的分析参数有关初始状态分析从零应力状态开始，即使是用到了先前的非线性分析的刚度分析可以从一个先前的非线性分析继续，初始状态为先前分析的所有荷载、变形、应力等结构响应所有的结构响应，如位移、内力、应力等直接与荷载成正比。不同线性分析结果可以叠加因结构属性可能发生变化，而且可能有初始非零应力状态，所以响应与荷载可能不成正比。不同的非线性分析结果一般不能叠加两者之间的区别其实很多，不过两个关键，一个是材料定义的时候不同，（材料属性根据需要设置，静力学分析一般只要弹性模量和泊松比，如果考虑体载荷或动力学分析还需要定义密度）。另一个就是在求解设置选项的时候不同，因为非线性一般存在收敛困难的问题。12 非线性行为的原因引起结构非线性的原因很多，主要可分为以下3种类型。（1）状态变化（包括接触）许多普通结构表现出一种与状态相关的非线性行为。例如，一根只能拉伸的电缆可能是松弛的，也可能是绷紧的；轴承套可能是接触的，也可能是不接触的；冻土可能是冻结的，也可能是融化的。这些系统的刚度由于系统状态的改变而突然变化。状态改变或许和载荷直接有关（如在电缆情况中），也可能是由某种外部原因引起的（如在冻土中的紊乱热力学条件）。接触是一种很普遍的非线性行为，接触是状态变化非线性类型中一个特殊而重要的子集。（2）几何非线性结构如果经受大变形，其变化的几何形状可能会引起结构的非线性响应。如图5.2所示的钓鱼杆，在轻微的载荷作用下，会产生很大的变形。随着垂向载荷的增加，杆不断弯曲导致动力臂明显减少，致使杆在较高载荷下刚度不断增加。 （3）材料非线性结构非线性中，最典型的分析是材料非线性（由于材料本身非线性的应力-应变关系导致的结构响应非线性叫材料非线性。除了材料本身固有的应力-应变关系外，加载过程的不同，机构所处的环境的变化均可导致材料的应力应变的非线性）材料非线性包括弹塑性分析，蠕变分析，超弹性分析，弹塑性分析。13 非线性结构有限元分析中应注意的问题1. 因为非线性，查看某个单独荷载的作用结果无意义。应将各种荷载放在同一种工况下进行分析。2. 在做几何非线性分析时，尽量不要使用释放梁端约束功能。建议在需要释放梁端约束的节点位置建立两个节点，节点间用弹性连接连接，在需要释放约束的方向不输入刚度值即可。几何非线性分析中不推荐使用释放梁端约束的原因如下： a. 几何非线性分析中的单元几何刚度是使用节点的坐标计算的(不断修正)。 b. 释放梁端约束后，因为梁单元的端部节点和另一个单元的节点共享一个节点，单元的位移和节点的位移会有不一致的问题，从而造成几何非线性分析不容易收敛。3. 几何非线性分析和P-delta分析不能同时进行，几何非线性分析属于大位移分析，P-delta分析属于小位移分析，做几何非线性分析时不必再做P-delta分析。如果用户这两个分析都要做，建议另存模型后分别分析。4. 几何非线性分析中采用的方法有TL法、UL法、CR法等。5.几何非线性分析中适用的单元有桁架(包括索)、梁单元、平面应力单元、板单元，如果与其他单元(如实体单元)混合使用时，只能考虑其他单元的刚度效应，不能考虑其他单元的几何非线性效应。6.施工阶段分析可同时考虑非线性累加模型和收缩和徐变分析，即非线性分析的累加模型可以考虑收缩和徐变。1. 4 钢筋混凝土结构非线性分析的意义 (1)由于钢筋和混凝土的抗拉强度相差很大，钢筋混凝土结构在正常使用状态下，大部分受弯构件部已经开裂而进入非线性状态，但钢筋并未屈服仍在弹性状态下工作，因此，作为一个结构或构件来说，必然是在非线性状态下工作，这时用弹性分析方法求得纳结构内力和变形就不能反映结构的实际工作状态。(2)混凝土和钢筋在一个结构中共同工作的条件是两者之间的变形协调，没有相对滑移，但实际上，这种条件并不能完全满足，特别是在反复荷载下，光圆钢筋与混凝土之间的粘结往往会破坏，某些情况下，会导致变形过大，而传统的线弹性结构分桥不能反映这些现象。(3)在钢筋混凝土结构的设计中，在内力分析时，往往按弹性计算，而在构件截面设计时，却按极限状态进行计算，其结果是内力分析和截面设计的结果都不能反映结构的实际受力状态，造成了钢筋混凝土结构内力分析和截面设计的严重脱节。(4)与其他任何形式的结构一样，节点和连接是保证钢筋混凝土结构能作为一个复杂体系承受外力的基本条件，而传统的弹性结构分析时将节点理想化为刚接(例如框架)或者铵接(例如衍架)均不能反映节点的复杂受力状态和变形情况，从而难以为设计提供正确的信息。(5)在长期荷载作用下，混凝土会产生一定的徐变变形，这时，结构的内力和变形就发生了变化，按弹性分析求得的内力和变形就不能反映实际情况了。2、混凝土结构非线性相关研究2. 1基于ABAQUS 纤维梁单元的钢筋混凝土柱受力破坏全过程数值模拟研究机构：沈阳建筑大学 华夏幸福基业股份有限公司项目简介1、关键词: 钢筋混凝土柱; 纤维梁单元; 单轴本构模型; 破坏准则; 单元消除; 显式动力2、研究方法：基于有限元软件ABAQUS 的显式求解模块Explicit，研究了钢筋混凝土梁柱构件纤维梁单元的构型及其单元消除技术，构建了包含材料破坏准则的钢筋、混凝土单轴本构关系，并利用用户材料接口VUMAT 编制了相应的计算子程序。为验证上述方法的合理性和正确性，本文对轴压与侧向往复加载下钢筋混凝土柱构件的受力破坏全过程进行了数值模拟。3、实例分析本文首先模拟了钢筋混凝土柱试件TP74 在两种不同工况下的受力破坏全过程。有限元建模时将构件等分为5 个单元，并将保护层混凝土、钢筋的箱型截面各划分为16 个纤维，核心混凝土的矩形截面划分为25 个纤维。图10 为在试验加载幅值60mm 范围内，试件计算所得与实测滞回曲线的对比。可以看出，计算曲线与实测曲线基本吻合，能够较好地模拟钢筋混凝土柱的承载力、刚度变化以及滞回捏拢行为，但在加载位移幅值较大时计算曲线较实测曲线略显丰满，且对强度退化现象反映不足。这可能是由于下述原因造成: 试验中每级荷载均加载三次，计算时只加载一次，因此对混凝土的累积破坏效应反映不足; 计算中未考虑钢筋的黏结滑移与屈曲; 计算中钢筋本构关系与实际的钢筋应力-应变关系有一定的差别。总体而言，基于B31 梁单元并采用本文建立的混凝土、钢筋单轴本构关系能够较为准确地模拟构件在压-弯荷载共同作用下的非线性性能发展过程。4、结论( 1) 采用基于梁单元建立的混凝土、钢筋单轴本构关系，可以较为准确地模拟钢筋混凝土柱构件在压-弯荷载共同作用下的非线性性能发展过程。( 2) 单元消除( 杀死) 技术能够简单、合理地描述构件由连续体变为非连续体的断裂破坏过程，是采用有限单元法模拟结构倒塌过程时的一种有效手段。( 3) 本文构建的材料单轴本构关系和破坏准则能够细致地描述梁柱构件受力过程中截面的内力变化情况，准确地模拟单个纤维的破坏与构件截面的渐进破坏过程。( 4) 对于非线性性态发展较为强烈的结构构件的模拟分析，考虑材料失效( 破坏准则) 可以更为准确、细致地描述因局部材料破坏造成的构件承载能力下降、耗能能力衰减等现象。( 5) 本文方法对于钢筋混凝土柱构件的模拟结果与其在真实受力状态下的破坏形态相吻合，较好地描述了钢筋混凝土柱构件的断裂破坏过程以及残余构件的运动状态。本文方法可用于钢筋混凝土框架结构在强震下倒塌全过程的模拟分析。3、预应力混凝土结构非线性相关研究31 预应力混凝土结构非线性有限元分析研究机构：同济大学结构工程与防灾研究所1、关键词:预应力混凝土;有限元法;非线性分析2、研究方法：提出了能反映预应力钢筋与混凝土间相对滑移的新型单元模型, 该模型可将无粘结预应力筋对结构的作用直接反映在单元模型内, 从而免去迭代试算过程。将该模型与混凝土、钢筋模型相结合, 所建立的材料和几何非线性分析程序能模拟预应力钢筋张拉, 进行预应力混凝土结构从开裂至破坏的全过程分析。文章最后给出了预应力混凝土连续梁的有限元分析实例, 对所提模型的正确性和可应用性进行了验证。3、实例分析计算过程首先模拟预应力筋张拉和锚固滑移, 在预应力筋左端加T =70kN 张力, 然后在反方向加T 的6 %, 模拟锚固滑移损失。张拉完成后的混凝土梁变形及预应力筋应力分别如图3a 和图5所示, 预应力筋应力值在1033MPa 至1087MPa 之间, 试验值为1041MPa 。张拉完成后, 对预应力筋进行灌浆粘合, 预应力筋单元换成有粘结单元。张拉完成至试验加载大致为28 天,根据计算钢筋松弛、混凝土收缩和徐变所引起的预应力筋应力值损失在10 %左右。混凝土梁的最终变形和破坏形态分别如图3b 及图3c 所示, 图4 为荷载与跨中挠度关系曲线。图1.试验模型图2.计算模型图3. 分析模型受力变形和破坏形态 图4.荷载-跨中挠度 图5.有粘结预应力筋应力分布图6.无粘结预应力筋应力分布4、结论本文提出了预应力钢筋与混凝土间相对滑移和相互作用的预应力筋单元模型, 可用来模拟预应力钢筋的施工张拉和无粘结预应力混凝土的计算分析, 免去了以往模型中预应力传递计算的迭代试算过程。结合该模型及混凝土、钢筋单元模型建立的预应力混凝土结构非线性分析有限元模型能反映结构的材料和几何非线性, 以及预应力混凝土结构所特有的预应力筋松弛、混凝土收缩和徐变等因素;可对由张拉时的摩擦、随时间变化的预应力筋松弛和混凝土徐变所造成的预应力损失进行较为细致的分析计算;可模拟有粘结或无粘结预应力混凝土结构从张拉施工、加载后开裂至破坏的全过程;可得到预应力筋的应力变化和应力分布的较为详细的分析结果。本文的分析实例所得到计算结果与试验吻合良好, 验证了所提模型的正确性和可应用性。3. 2预应力混凝土结构组合式非线性分析模型研究机构：东南大学混凝土&预应力混凝土教育部重点实验室， 华东建筑设计研究院有限公司1.关键词：预应力混凝土；非线性分析；组合式模型；纤维截面；粘结滑移；低周疲劳；抗震性能2、研究路径：建立基于纤维截面的预应力混凝土结构组合式分析模型。开发了一种粘结单元模拟预应力筋与混凝土间的粘结滑移，采用基于纤维截面的粘结单元模拟普通纵筋与混凝土的粘结滑移，模型还考虑了材料低周疲劳特性对结构损伤的贡献。利用3 组预应力混凝土框架试验对模型进行验证。3.实例分析：利用本文的组合式模型对试验进行分析，得到的各层跨中竖向位移与竖向力关系如图7 所示。通过图7 可以看到：理论计算得到的承载力略大于试验结果，但两者总体上吻合较好，组合模型能够准确地预测预应力混凝土框架在竖向荷载下的受力性能。4结论文章建立了基于纤维截面的预应力混凝土结构组合式分析模型，综合考虑了钢筋粘结滑移，材料低周疲劳等影响结构非线性性能的主要因素，并通过3 个有代表性的试验进行了验证，得到以下结论：1、模型不仅能够准确地模拟预应力混凝土结构的竖向变形能力和滞回特性，也能够准确地预测预应力混凝土结构在地震作用下的动力性能。2、模型的建立是基于杆系纤维截面单元，因而具有计算效率高，收敛性好的优点，能够作为预应力混凝土结构非线性分析的有效手段。3、模型没有考虑节点区域的剪切变形，无法反映节点破坏模式。鉴于节点损伤对结构抗震性能具有重要的影响，模型仍需加以完善。4、桥梁结构非线性相关研究41 钢筋混凝土纤维梁柱单元实用模拟平台研究机构：天津大学建筑工程学院，滨海土木工程结构与安全教育部重点实验室1、关键词：钢筋混凝土；梁柱单元；纤维模型；非线性行为；滞回性能；约束效应2、研究路径：将基于ABAQUS 建立钢筋混凝土精细化纤维梁柱单元模拟平台FENAP，编制相应的材料库，开发多种钢和混凝土材料本构模型。同时利用FENAP 平台分别模拟梁、柱构件在反复加载条件下的滞回性能，并将计算结果与OpenSees 所得结果进行对比分析，以期得到一些有益的结论。3、算例分析：模拟过程中，采用位移控制、反复循环加载。FENAP 与OpenSees 计算得到的墩顶力-位移、墩底截面弯矩-曲率以及纤维的应力-应变关系对比曲线如图8(a)图8(c)所示。4、结论本文建立了一种基于ABAQUS 的钢筋混凝土精细化纤维梁柱单元的模拟平台FENAP，通过对不同梁柱构件滞回性能的数值模拟，得到以下结论：(1) FENAP 平台界面友好，方便实用，能有效模拟构件的刚度退化、强度退化等损伤效应以及轴力和弯矩的多维耦合效应等复杂非线性动力行为，且可考虑箍筋对混凝土的约束作用。(2) 在FENAP 平台中建立了独立的单轴材料库，添加了多种材料本构模型，能模拟钢和混凝土材料的多种非线性行为，且设定了可扩展的材料库接口，方便后期添加更多的用户材料本构模型。(3) FENAP 平台可有效模拟桥梁构件的多种复杂非线性行为，且计算效率和求解精度高，为长大桥梁地震灾变过程模拟提供了一种实用分析手段。42 预应力混凝土薄壁高墩刚构桥梁极限承载力分析研究机构：内蒙古大学交通学院.交通运输部科学研究院.中交通力建设股份有限公司，1、关键词：桥梁工程；高墩刚构桥梁；理论分析；极限承载力；虚拟层合单元；非线性2、研究方法采用实体退化的虚拟层合板单元，综合考虑其结构几何非线性和钢筋、混凝土材料非线性效应，同时计入混凝土开裂后的刚度折减，在虚拟层合单元有限元程序的基础上，引入破坏准则（混凝土材料）和屈服模型（钢筋）及正交分布裂缝模型，构建预应力混凝土结构空间非线性有限元极限承载力分析程序。为验证该程序对钢筋混凝土结构的有效性，对一预应力空心板梁的加载破坏试验过程进行全程模拟。利用本文方法对预应力混凝土高墩刚构桥梁的极限承载力进行了分析，并与通用有限元程序的计算结果进行比较。3、算例分析 4、结论（）采用实体退化的虚拟层合板单元，引入破坏准则（混凝土材料）和屈服模型（钢筋）及正交分布裂缝模型，构建了预应力混凝土结构空间非线性有限元极限承载力分析程序，模拟了一简支空心板梁破坏试验过程，在预测其极限承载力和结构的整体性能方面具有较好精度。（）本文方法可有效克服由于混凝土开裂后刚度为的单元数目较多时造成的整体刚度矩阵“假性奇异”引起的计算不收敛现象。同时，由于模型单元数目少，可避免由于单元划分过细引起的裂缝应力局部失真，较通用有限元程序可更合理地模拟较复杂混凝土高墩刚构桥的极限承载力及破坏过程，计算结果符合对钢筋混凝土结构破坏特征的常规认知，具有合理性。（）本文方法相对常规通用有限元程序的主要优势是对结构开裂后的刚度矩阵构成进行了更有效地处理，使之可以较准确地预测结构出现裂缝后的受力行为。但本文未对单元内部块体划分对计算精度的影响进行深入分析，需要进一步深入研究。5、总结与展望混凝土结构是土木工程中应用最广泛的一种结构。这种结构在建造和使用过程中的安全性、经济性和技术合理性有着重要的社会和经济意义。然而，由于钢筋混凝土是由两种力学性能差别较大的材料组成，受力极为复杂，因而目前除少数结构(如钢筋混凝土连续梁、框架等)外，人们仍用线弹性理论来分析钢筋混凝土结构的内力，而以极限状态设计法确定构件的承载能力，使结构内力分析和构件截面承载力计算结果均不能反映结构的实际受力状态。近年来随着土木工程学科的成熟和有限元法的不断发展，国内外许多学者采用数值模拟技术进行结构破坏模式与破坏行为预测等混凝土结构非线性行为的研究，学者们在研究中积累了复杂和多样的混凝土本构模型，与此同时也研究出了多种能够反映混凝土构件特点的单元模式，并且发展了较先进的迭代控制方法，我们越来越强调结构整体非线性分析。此外，随着设计理论和施工技术的不断发展，土木结构不断的朝着更大高度和更大跨度的方向发展，许多超级工程的模拟需要将非线性分析与有限元法结合起来，我们除了在混凝土的本构关系的表述和试验研究方法而继续进行更深入的研究之外，钢筋混凝土结构非线性有限元分析也应进一步向实用方向发展，努力把现有的分析方法与工程设计结合起来。为造福广大人民群众和促进学科领域前沿的不断发展而努力奋斗。参考文献1、王强,宋雪迪,郝中华,王占飞.基于ABAQUS纤维梁单元的钢筋混凝土柱受力破坏全过程数值模拟J. 土木工程学报. 2014(12)2、吴晓涵.预应力混凝土结构非线性有限元分析J. 建筑结构学报. 2007(04)3、于琦,孟少平,吴京,郑开启.预应力混凝土结构组合式非线性分析模型J. 工程力学. 2011(11)4、禚一,李忠献.钢筋混凝土纤维梁柱单元实用模拟平台J.工程力学.2011(04)5、李海青,杨万里,高璇.预应力混凝土薄壁高墩刚构桥梁极限承载力分析J. 中国公路学报. 2013(06)袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羁膂莈蚅袇膁蒀袀螃膀薂蚃肂腿节衿羈腿莄蚂袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羁膂莈蚅袇膁蒀袀螃膀薂蚃肂腿节衿羈腿莄蚂袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羁膂莈蚅袇膁蒀袀螃膀薂蚃肂腿节衿羈腿莄蚂袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羁膂莈蚅袇膁蒀袀螃膀薂蚃肂腿节衿羈腿莄蚂袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羁膂莈蚅袇膁蒀袀螃膀薂蚃肂腿节衿羈腿莄蚂袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈