基于有限条法的长跨斜拉桥地震分析.docx

基于集成有限条法的大跨斜拉桥的地震分析摘要：本文提出了一种基于集成有限条法和时程分析法的用于大跨斜拉桥进行地震分析的方法。该方法降低了结构集成时间、减少了结构的自由度从而提高了计算分析的速度。现有的有限条法仅对桥面进行建模，而将其他结构构件由边界条件代替。使用集成有限条法，可以对大跨斜拉桥的所有构件建模。因此，可以考虑包括桥面，桥墩和拉索之间的相互作用的真实的动态响应。为了验证本文提出方法的有效性，对一座大跨斜拉桥进行了在均匀和非均匀地震荷载下的地震反应分析。结果表明，集成有限条法可以成功应用于大跨度斜拉桥的地震分析，并且可以显著地降低分析时间。关键词：集成有限条法;大跨斜拉桥;地震分析; 时程分析1前言为了满足经济和社会发展对安全和高效的运输系统的需要，世界各国普遍大力发展了大跨度桥梁。如今，由于大跨度斜拉桥的经济性，与悬索桥比起来引起了更多的关注。由于刚度低，现代大跨度斜拉桥对动荷载非常敏感。因此，解决抗风和抗震设计是关键。大跨度桥梁在地震荷载和风荷载作用下的结构响应是非常复杂的，需要有与精确分析和相应的计算手段。大多数桥梁结构分析是通过有限元软件实现的。过去几十年，有限元法的准静态数值分析在研究和设计大跨桥梁中得到广泛应用（Cook等人，1989）。随着跨度的增加和轻质复合材料的应用，这些结构的动态特性变得更加复杂，使得常规有限元法不能进行正确分析。由Y. K. Cheung在20世纪70年代（Cheung和Cheung 1971）所提出的有限条法作为一种半分析数值分析，已经显示了分析桥梁结构的巨大潜力。然而，由于一些限制，这种方法不像有限元法那样流行。例如，有限条法只能处理棱柱形构件且为简支约束支撑的桥梁。此外，该方法不能处理具有剪力、内部支撑和集中荷载作用的情况。尽管有限条方法在结构分析中有很多实际应用，但都局限于桥梁的上部结构的地震分析中。桥梁的桥墩和拉索通常用边界条件代替。一种考虑桥墩和拉索影响的方法是将桥墩和拉索用其他单元（如边界单元）建模，通过迭代获得桥面和桥墩之间的相互作用。然而，该方法仅能进行静态和准静态分析。为了对复杂桥梁结构进行动态分析，Cheung等人2009; Shen等人2013; Naderian等人2015）提出了集成有限条法（IFSM）使得有限条法可用于大跨斜拉桥的整体动态分析。Shen等人比较了集成有限条法和有限元法分析所需的计算时间。在相似的建模和网格划分条件下，当通过集成有限条法对中等跨度桥梁进行动态分析时，用时274分钟，而有限元法则用时298分钟，集成有限条法大约提高10的计算速度。并且，在集成有限条法中每条的单元数可以进一步降低，从而提高计算效率。当对具有大量自由度的大跨斜拉桥进行动态分析时，与有限元法比所节约的时间将更加可观。集成有限条法可以对具有各种类型的约束支撑和内部支撑的大跨斜拉桥进行快速和准确的地震分析。在大跨桥梁的动态分析和设计方面，应特别考虑在整个桥梁寿命期内强震和强风的作用。近几十年来，在大跨结构工程中地震随机分析方面取得了重大进展。在有限条动力分析领域，Cheung和Cheung（1971）首先提出了用于具有不同边界条件的薄壁结构的自由振动分析的有限条法。Dawe（2002）基于薄板理论使用常规和样条有限条法对复合层压板，棱柱板和壳结构的进行了广泛的研究。该研究使用有限条法分析了单跨和多跨结构的自然频率，并评估了平板对动态荷载的瞬态响应。Wang和Zhang（2005）介绍了一种用于复合板自由振动分析的分层B样条有限条法。在该方法中，复合层叠板沿厚度方向分成多个层。在风致动力响应方面，Lau等人将有限条法与有限元法相结合应用到了大跨斜拉桥的颤振分析中。Li（1988）和Cheung等对任意形状的板和壳进行了自由振动分析。Cheung（2009）和Shen（2013）在频域法中引入了用于板梁和箱梁桥动态分析的集成有限条法。基于集成有限条法，本文对大跨斜拉桥进行了地震荷载作用下的时程分析。首先，简要介绍了大跨度斜拉桥的刚度，质量和阻尼矩阵的推导。使用从集成有限条法导出的刚度和质量矩阵，通过使用Newmark法进行了时程分析。通过比较有限元法获得的固有频率和从现场调查中报告的固有频率，对集成有限条法在大跨斜拉桥动态分析中解的精度进行了评估。此外，研究了均匀，不均匀和变化的地震激励对桥的结构特性的影响。数值计算结果表明，大跨度斜拉桥动力分析中集成有限条法的收敛性和精度非常高，证实了该数值技术的效率和能力。2位移函数通过在纵向上应用B3样条函数（Prenter 1975）获得在集成有限条法中。横向的条带的位移与有限元法相同，其中三次多项式用于表示条的横向形状函数。假设桥面是一个板壳结构，在分析中包括面内和面外自由度。平板样条有限条如图1所示，其中节点线的每个节点具有四个自由度，三个平移u，v和w，以及一个转动。平板条的总势能可以从平面内变形的和平面外变形的代数和得到。图1 平板样条有限条以为中心的平板条的位移参数向量由下式给出（1）集成有限条法是基于使用不等间隔的B3样条函数，通过该函数得到过渡段并将不同的结构构件连接在一起。此外，将支撑位置和集中载荷作用的位置与节点线上的节点重合，从而能得到更准确的结果。此外，引入不等间隔的内部结使得能够通过更紧密地间隔结点来描述在高应力梯度的区域中或在几何形状突变处的准确响应。以为中心的不等间隔的B3样条函数可以表示为(2)其中：(3)膜位移函数u和v以及弯曲位移函数w由横向多项式和纵向B3样条的乘积表示为如下所示：（4）（5）（6）其中r 为节点线上的纵向截面的总数。在上述等式中采用的横向形函数是用于垂直位移变化和用于平面内位移的线性内插的三次Hermite多项式函数，如下：（7）其中X=x/b. 到为一维行矩阵，并且每个矩阵具有（m + 3）个局部B3样条。、和分别与节线i和j的位移u和v相关，、和与位移w有关。桥墩和塔柱用与平板条类似的柱条来建模，并模拟为在一端固定，用于提供支撑边界条件，并且另一端为自由的垂直条，如图2所示。事实上，柱状条全局z方向与相应方向上由平面内刚度控制的平面壳体样条有限条中的局部y方向类似，而全局方向类似于壳体样条有限条的局部z方向。利用假定的位移函数和预设边界条件，可以通过使用传统的有限元概念来得到柱条的形函数。图2垂直条有时，一维柱条更适用于桥墩和塔柱的建模，因为它只有一条节点线并且类似于有限元法中使用的梁单元。属于节点线的每个节点具有三个平动自由度。因此，可以将位移函数定义为（8）（9）（10）对于位移-应变关系，因为在一维条中剪力和扭矩非常低而仅考虑在垂直和横向方向上的弯曲和轴向应力。为了在有限条中对拉索进行建模，开发了索条。通常，拉索只能承受轴向拉力，因此索条中只定义了轴向应力的应变-位移关系。在有限条法中所有条都必须预设边界条件。转换部分的建模利用有限元法中的单元和节点的概念，很容易建立其它构件的单元来将与桥面板连接在一起。然而，纵向方向上的概念单元与普通有限条法或样条有限条法中定义的不同。为了解决这个问题，在集成有限条法种开发了一种特殊的转换部分，用于连接复杂结构的不同构件，例如大跨斜拉桥。转换部分是通过使用不等间隔的B3样条函数来建立的。此外，支座可以用转换部分的特殊边界条件来建模。连接两个不同构件的典型转换部分如图3所示。不失一般性，假设正常和转换部分的宽度分别为H和h。例如，垂直线可以称为桥墩条上的节线，水平线可以称桥面板条上的节线。垂直线和水平线分别在桥面板条和桥墩条的结3和8处重叠。为了对允许转动但限制平动的固定支撑进行建模，结3和8应当具有相同的位移值以实现兼容性。为了在结3和8处具有相同的位移，h / H应当无限小。使用在样条有限条中建立的转换部分，满足结构的不同构件的位移协调性。图3转换部分（图是错的）刚度，质量和阻尼矩阵在集成有限条法中定义了大跨斜拉桥中所有结构构件的位移函数之后，可以通过使用标准有限元法计算不同构件的刚度k和质量m矩阵，如下所示：（11）（12）其中为密度，t 为条的厚度; D和B分别为弹性矩阵和应变矩阵; N为形函数矩阵。类似于有限元法，在建模过程中将集成有限条法条内部位移转换为节点位移，但是由于集成有限条法的半分析属性，所需节点的数量显着减少。因此，可以使用标准形成总刚的方法形成总体刚度矩阵K和总体质量矩阵M。有不同的方法来评估阻尼矩阵C。本文采用了经典的瑞利阻尼，并将桥梁的阻尼矩阵C表示为桥梁的刚度和质量矩阵的函数，如下所示：（13）其中和是瑞利阻尼因子，其可以通过具有与两个特定频率相关联的两个结构阻尼比j来研究。结构的第n阶模态阻尼比j由下式给出（14）假设对于两种不同模态具有相同的阻尼比j，可以计算和。数值研究（分析）一般情况本文的研究对象是香港的汲水门大桥，如图4所示。它是一个超大跨斜拉桥结构，总长750米，主跨度为430米。桥结构包括两个高度为145和133m;混凝土塔架，四个分别为24.8,32.9,16.9和26.8m的混凝土桥墩，176根拉索; 数个支座和一个桥面板。为了研究桥梁在地震荷载下的响应，建立了有限条模型。图4 汲水门桥桥梁的有限条法建模汲水门桥的主跨由复合钢混凝土结构组成，边跨由混凝土箱梁制成。桥梁模型所用的材料性质如表1所示，Zhang等人（2001），几何性质如图5所示。桥面板由四块壳样条条组成，所有条分为几个部分。这些部分的长度沿着桥面板变化，使得结的位置与拉索的的端部一致。表1桥的材料性能力学参数主跨边跨弹性模量 (kPa)2.00 x 1083.00 x 107密度 (kg/m3)3,8803,630泊松比0.30.2转动惯量（垂直） (m4)191363转动惯量（横向） (m4)2,5305,560图5桥面板的几何性质：（a）顶视图; （b）前视图塔柱的网格划分如图6所示。塔柱结构采用1D柱条单元。在每个塔柱中有三个支柱，其连接塔柱的上部，中间和下部。连杆梁是刚性的，并且所有六个自由度都被约束。然而，在1D柱条中，每个结仅具有三个平动自由度，而不考虑转动自由度。桥的四个墩也由1D柱条模型建模。塔柱和墩的边界条件都固定在与地面的连接处。176条拉索假定为线弹性单元，并且由有限条模型中的缆条建模。所有拉索通过将热荷载转换成在条的两端施加的等效拉力来施加预应力。补图图6集成方法的验证为了验证集成有限条法用于大跨斜拉桥的动态分析的性能，获得了桥的前10个固有频率（在表2中列出）。前两阶振动模态发生在塔柱，第三阶是甲板的垂直弯曲，接下来的三阶模态是桥面板的竖向弯曲模态、横向弯曲模态和扭转模态。将有限条结果与通过有限元法获得的结果和对桥面板振动的现场测量测试结果进行了比较。其中有限元分析是通过SAP2000软件进行的。在SAP模型中，桥面板由壳单元建模，墩和塔柱由梁单元建模，拉索由桁架单元建模。从表2中可以清楚地看出，在结果之间存在非常好的一致性，这在动态分析方面证实了集成有限条法的准确性。当使用集成有限条法时，发现固有频率的值略高，特别是桥面板的第一横向模态和第一扭转模态。当将从集成有限条法获得的固有频率与来自有限元分析的固有频率进行比较时的最大误差百分比为4.8，而当将相同结果与现场测量测试数据进行比较时，最大误差百分比小于10。一般来说，现场测量的自然频率的值略会低于集成有限条法和有限元法获得的值（或数值模拟的值）。大跨斜拉桥的自然频率对桥梁动态性能有重要影响。因此，固有频率的高一致性非常重要; 此外，结果之间良好的一致性证实了集成有限条法能够获得桥梁的动态特性。因此，可以非常有效对大跨斜拉桥进行地震分析。地震分析使用集成有限条法对汲水门大桥的从几个角度进行了地震分析研究。并将一些结果与有限元分析结果进行对比验证。在有限元SAP模型中，当进行动态分析时采用直接积分法。研究的关键点如图7所示。图7研究的关键点大跨斜拉桥的地震波效应当桥梁的跨度较小时，可以假定所有支撑在地震荷载下均匀地移动，或将受到的地震激励视为均匀激励。相比之下，各种空间效应，如行波效应，部分相干效应和局部场地效应，可能在大跨桥梁的地震反应中发挥重要作用。行波效应是一种动态现象，其中地震波在不同的时间到达不同的支撑。部分相干效应是指由于多重反射和折射波传播通过高度不均匀的土壤介质而导致的地震波相干性的损失。此外，由于不同的局部土壤条件，地震动的振幅和频率含量的变化被称为局部效应（Saxena et al。2000）。行波效应对大跨度桥梁的地震反应特别重要，而部分相干效应的重要性相对较小（De Silva 2005）。均匀激励下的地震反应作为最简单的情况，在相同方向上的所有自由度在相同的地面运动下振动。选择一条人工正弦波和一条天然波，并研究对集成有限条法的影响。在正弦波地面加速度的作用下，获得桥东塔的F点处的纵向加速度和位移输出结果（分别示于图8和图9），并将它们与有限元时程分析结果进行了对比。在均匀激励下的时程分析中，可以看出集成有限条法的结果与有限元的结果完全一致。从曲线可以看出，位移和加速度的幅度都在增加，前10s的增长速率高于随后的时间间隔。然而，从有限条时程分析得到的加速度和位移响应的幅度比有限元结果略高。图8 桥东塔的F点处的纵向加速度图9 桥东塔的F点处的纵向位移天然波选择集集地震波，并加给桥的横向的所有自由度。在集成有限条法和有限元中使用相同的分析方法，计算了桥面板中的点B处的横向加速度和位移，结果如图10和图11所示。与正弦波激励（图8和图9）相反，在集集地震下的地震反应并不均匀。如图11和12所示，在15s到20s之间的周期中遇到最大位移和最大加速度。值得注意的是，在这段时期之前的位移量非常低，这与同一时间间隔的集集地震的低加速度输入一致。同样，用有限元和集成有限条法的得到的时程分析的数值结果非常一致。图10B处的横向加速度图11 B处的横向位移非均匀激励下的地震反应当考虑大跨桥梁的不均匀激励时，应在不同支撑处输入不同的运动。使用集成有限条法和时域方法时，研究了行波效应和变化的地震激励对桥的影响。行波效应行波效应解释了地震动的时间延迟，因为地震波不同时到达桥基础，特别是对于大跨度桥梁这种大型结构，其中桥墩和塔柱之间的距离很长。在非垂直波在不同时间到达地面处的不同点并且在这些点处的运动之间产生时移的情况下，发生行波效应或异相效应。采用集集地震波作为输入的天然波;然而，对于参数研究，假设地震波的速度在纵向方向上介于500m / s到1000m / s之间。假设地震波从第一个西部桥墩开始传播，之后地面运动将逐一传递到后面的塔柱和桥墩。在有限条模型中，这是通过在各种支撑处输入具有不同时间滞后的地震波来实现的。还采用具有速度无限大的地震波，这意味着当在不同支撑之间传输地震波时没有时间延迟，这等于均匀激励。当施加非均匀激励时，在具有不同波速的情况下，位于西塔柱顶部的点E的加速度和位移响应为如图12至图17所示。由地震曲线可以看出，当桥处于均匀激励作用时，得到最大响应（图14和17）。此外，对于所有研究的速度，在15s到20s之间的周期中确定位移和加速度的最大幅度。可以得出结论，在均匀激励下，桥梁的位移和加速度高于非均匀激励的位移和加速度。随着地震波的速度增加，加速度的振幅增大。如图15和16所示，当速度从500增加到1000m / s时，在点E处记录的最大位移减小。图12 速度为500m / s时西塔柱点E处的纵向加速度图13 速度为1000m / s时西塔柱点E处的纵向加速度图14 速度为无穷大时西塔柱点E处的纵向加速度图15 速度为500m / s时西塔柱点E处的纵向位移图16速度为1000m / s时西塔柱点E处的纵向位移图17速度为无穷大时西塔柱点E处的纵向位移地震激励在时程分析中，数值结果高度依赖于所选择的地面加速度记录。为了研究不同地震激励下大跨斜拉桥支撑处的动力响应，使用台湾强运动阵列（SMART-1）的地震记录。SMART-1是由地球科学研究所（台湾）和加利福尼亚大学伯克利分校建造的强运动地震仪的密集数字阵列。SMART-1中的地震记录是在地震发生期间在不同地方观测到的原始数据。因此，记录已经包括了部分相干效应，衰减效应和场地效应的影响。使用在距离为1,000米的两个地方观察到的加速度记录，输入1和输入2。两个加速度记录的垂直和纵向分量如图18和图19所示。图18 地震波的垂直分量：（a）输入1; （b）输入2图19地震波的纵向分量：（a）输入1; （b）输入2进行三种情况分析。对于前两种情况，输入1和输入2分别施加在整个桥梁的支撑处。因此，这两种情况对应于均匀激励。在第三种情况下，输入1施加于两个西桥墩和西塔柱，而输入2应用于两个东桥墩和东塔柱。获得了在桥面板中部的点A处的垂直位移响应，以及在东塔柱顶部的点G处的纵向位移响应，并且示于图20-22和图23-25中。因此，可以比较不同案例研究中的地震反应，可以看到，在均匀激励（图20,21,23和24）下，点A和G处的位移数高于非均匀激励（图22和25）。在几乎所有情况下，在从开始到t = 5s的时间内目标点处的位移为零。此外，最大位移发生在10s至15s之间。看起来在所有三种情况下，位移曲线收敛到零。图20 A点的垂直位移响应：输入1和输入1图21A点的垂直位移响应：输入2和输入2图22A点的垂直位移响应：输入1和输入2图23 G点的垂直位移响应：输入1和输入1图24G点的垂直位移响应：输入2和输入2图25G点的垂直位移响应：输入1和输入2结论当使用传统的有限条法进行大跨斜拉桥动态分析时，仅对上部结构进行建模，而将支撑和拉索用边界条件代替，不能考虑桥面（上部结构）和其他构件间的相互作用。而在本文提出的集成有限条法中，可以对整个桥梁进行统一建模，可以考虑结构构件之间的结构相互作用。袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羁膂莈蚅袇膁蒀袀螃膀薂蚃肂腿节衿羈腿莄蚂袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羁膂莈蚅袇膁蒀袀螃膀薂蚃肂腿节衿羈腿莄蚂袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羁膂莈蚅袇膁蒀袀螃膀薂蚃肂腿节衿羈腿莄蚂袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羁膂莈蚅袇膁蒀袀螃膀薂蚃肂腿节衿羈腿莄蚂袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羁膂莈蚅袇膁蒀袀螃膀薂蚃肂腿节衿羈腿莄蚂袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿