基于小型振动台试验探究地震损伤对砖塔动力特性的影响

基于小型振动台试验探究地震损伤对砖塔动力特性的影响一、引言1.1研究背景与意义砖塔作为一种具有深厚历史文化底蕴的古建筑，承载着丰富的历史信息与文化价值，是人类文明的瑰宝。它们不仅是建筑艺术的杰出代表，体现了不同历史时期的建筑风格和工艺水平，更是历史变迁、文化传承的见证者，反映了当时的社会、经济、宗教等多方面的情况，具有不可替代的文物价值。然而，砖塔多为高耸结构，高宽比较大，材料多为砖石，其抗拉、抗剪强度相对较低。在地震等自然灾害作用下，砖塔容易遭受严重破坏，甚至倒塌。如1976年唐山大地震，大量砖石结构建筑包括砖塔遭受重创，许多塔体出现裂缝、倾斜，部分甚至完全坍塌，造成了不可挽回的文化损失。2008年汶川地震同样对灾区及周边地区的砖塔带来了巨大冲击，众多珍贵的砖塔建筑面临严重的损毁威胁，这些地震灾害的实例表明，砖塔在地震作用下的安全问题亟待关注。为了更好地保护这些珍贵的文化遗产，深入了解地震损伤对砖塔动力特性的影响至关重要。动力特性作为结构的固有属性，反映了结构在动力荷载作用下的振动特征，包括自振频率、振型和阻尼比等参数。地震损伤会使砖塔的结构材料性能劣化、构件连接松动、结构整体性降低，从而导致其动力特性发生改变。通过研究地震损伤对砖塔动力特性的影响，可以揭示砖塔在地震作用下的损伤演化规律和破坏机制，为砖塔的抗震性能评估、保护加固以及抗震设计提供科学依据。小型振动台试验是研究结构在地震作用下动力响应和破坏机制的重要手段之一。它能够在实验室条件下模拟不同地震波的作用，对砖塔模型进行加载试验，实时测量和记录模型在地震作用下的加速度、位移、应变等响应数据。与现场原位测试相比，小型振动台试验具有可重复性强、试验条件可控、成本相对较低等优点，可以系统地研究不同地震参数、结构形式和材料特性等因素对砖塔动力特性的影响。通过对试验数据的分析，可以深入了解砖塔在地震过程中的力学行为，验证和完善理论分析模型，为实际砖塔的抗震保护提供可靠的试验支持。因此，开展地震损伤对砖塔动力特性影响的小型振动台试验研究具有重要的科学意义和工程应用价值。1.2研究现状1.2.1砖石古塔抗震性能研究砖石古塔作为珍贵的历史文化遗产，其抗震性能一直是研究的重点。国内外学者通过震害调查、理论分析、数值模拟和试验研究等多种方法，对砖石古塔的抗震性能进行了广泛而深入的探讨。在震害调查方面，许多学者对不同地区、不同类型的砖石古塔在地震中的破坏情况进行了详细记录和分析。如潘毅等人对汶川地震中砖石结构古塔的震害进行了调查，发现地震对砖石古塔的破坏形式主要包括塔身开裂、倾斜、倒塌等，且破坏程度与地震烈度、结构形式、材料性能等因素密切相关。袁建力研究了砖石古塔震害程度与地震烈度的对应关系，为评估砖石古塔在不同地震烈度下的安全性提供了重要参考。理论分析上，学者们运用结构力学、材料力学等知识，建立了各种理论模型来分析砖石古塔的抗震性能。张永亮等探讨了应用振型分解反应谱法分析不同高度砖石古塔地震性能的适用性，为砖石古塔的抗震设计和分析提供了理论依据。李胜才等应用显式积分法分析了古塔结构破坏的大转动与非线性大位移以及不连续位移场，反演了地震作用下古塔的损伤机制。数值模拟技术的发展为砖石古塔抗震性能研究提供了有力工具。陈平等利用有限元分析软件对崇文塔进行了模拟分析，探讨了造成崇文塔损害的原因。芦苇等应用时程分析法分析和评估了不同基础刚度条件下万寿寺塔的抗震性能。通过数值模拟，可以对砖石古塔在地震作用下的应力、应变分布以及动力响应进行详细分析，深入了解其抗震性能和破坏机理。1.2.2振动台试验研究振动台试验能够在实验室条件下模拟地震作用，直观地研究结构的动力响应和破坏过程，对于砖石古塔抗震性能研究具有重要意义。谢启芳等对西安钟楼设计制作了1∶6的振动台试验模型，研究了钟楼模型在不同地震波作用下的抗震性能。Kim等对一个石结构古塔制作了足尺的振动台试验，研究石结构古塔的地震响应规律。赵祥等设计并制作1∶10的古塔模型进行振动台试验，研究了模型地震激励前后的动力反应。宋晓滨等设计制作了1∶5的木塔模型振动台试验，研究了不同等级地震作用下结构的动力响应规律。卢俊龙等设计制作了1∶8的玄奘塔模型，并进行了振动台试验，得到内填黏土的实心古塔在输入单向、双向以及三向地震波时，结构的地震响应规律和破坏机制。这些振动台试验研究，为深入了解砖石古塔在地震作用下的动力特性、破坏模式和抗震性能提供了大量的试验数据和宝贵的经验，也为数值模拟和理论分析的验证提供了依据。1.2.3砖砌体结构数值模拟研究砖砌体结构是砖石古塔的主要结构形式，对砖砌体结构的数值模拟研究有助于更好地理解砖石古塔的力学行为。在数值模拟方法上，常用的有离散元法（DEM）、有限元法（FEM）、格子Boltzmann方法（LBM）以及混合方法等。离散元法通过计算颗粒间的相互作用力来模拟砌体单元间的相互作用，能够处理不规则形状的颗粒，适用于较大的形变变形分析场景。有限元法是一种较为成熟的数值模拟方法，常用于分析结构中的静力和动力响应，在砌体匀质化研究中主要用于研究砌体的永久变形和塑性变形，虽然其需要进行网格划分，计算系统较为复杂，且网格精度影响计算精度，但模拟结果较为精确，能够更好地模拟颗粒的变形过程。格子Boltzmann方法是一种较新的数值计算方法，由于其使用的离散方程，能够模拟流体的宏观和微观行为，在研究砌体匀质化中的流体力学过程时表现出了很高的精度和可靠性。一些研究人员尝试将不同方法进行融合，如DEM-FEM方法，通过将离散元方法和有限元方法融合，能够更好地模拟砖墙的弹性和塑性变形过程，提高了研究的真实性和准确性。许多学者运用这些数值模拟方法对砖砌体结构在不同荷载作用下的力学性能进行了研究。通过建立合理的数值模型，可以模拟砖砌体结构的开裂、破坏等现象，分析其受力性能和变形特性，为砖石古塔的抗震性能研究提供理论支持。1.2.4研究不足尽管目前在砖石古塔抗震性能、振动台试验以及砖砌体结构数值模拟等方面取得了一定的研究成果，但仍存在一些不足之处。在砖石古塔抗震性能研究中，虽然对震害现象和破坏机理有了一定的认识，但对于一些复杂因素的影响，如材料的劣化、结构的损伤累积以及不同结构形式之间的相互作用等，还缺乏深入系统的研究。现有理论模型和分析方法在准确性和适用性方面仍有待提高，难以全面准确地描述砖石古塔在地震作用下的复杂力学行为。振动台试验研究方面，由于试验条件和模型制作的限制，一些试验结果可能存在一定的局限性。例如，试验模型与实际砖石古塔在材料性能、结构尺寸和连接方式等方面可能存在差异，导致试验结果不能完全反映实际情况。此外，目前振动台试验多集中在少数典型古塔，对于不同地区、不同类型砖石古塔的研究还不够全面。在砖砌体结构数值模拟研究中，数值模型的建立和参数选取对模拟结果的准确性影响较大，但目前还缺乏统一的标准和方法。不同数值模拟方法都有其各自的优缺点和适用范围，如何根据具体问题选择合适的方法以及如何将不同方法进行有效结合，还需要进一步探索和研究。综上所述，针对地震损伤对砖塔动力特性影响的研究，仍有必要通过更深入的小型振动台试验以及更完善的数值模拟分析，全面系统地揭示其内在规律，为砖塔的抗震保护提供更加科学可靠的依据。1.3研究内容与方法1.3.1研究内容砖塔模型设计与制作：根据实际砖塔的结构特点和尺寸，按照相似理论设计并制作1:X比例的小型砖塔模型。在模型制作过程中，选用与实际砖塔材料性能相似的材料，如采用特制的小型砖块和模拟的砂浆，确保模型能够真实反映实际砖塔的力学性能。同时，对模型的尺寸、材料性能等参数进行精确测量和记录，为后续试验提供准确的数据基础。振动台试验方案设计与实施：制定详细的振动台试验方案，选择合适的地震波，如El-Centro波、Taft波等，并根据不同的地震烈度对地震波进行调幅，设置多个加载工况。在试验过程中，利用传感器实时测量砖塔模型在地震作用下的加速度、位移、应变等响应数据，记录模型的裂缝开展、破坏形态等现象。同时，在每次加载前后，采用环境激励法测试模型的动力特性，包括自振频率、振型和阻尼比等参数，以研究地震损伤对砖塔动力特性的影响规律。试验数据处理与分析：对振动台试验获取的数据进行整理和分析，绘制加速度时程曲线、位移时程曲线、应变时程曲线等，分析砖塔模型在不同地震作用下的动力响应特征。通过对比不同加载工况下模型的动力特性参数，研究地震损伤对砖塔自振频率、振型和阻尼比的影响规律。采用损伤识别方法，如基于动力特性变化的损伤识别方法、基于应变模态的损伤识别方法等，对砖塔模型的损伤程度和损伤位置进行识别和评估。数值模拟分析：利用有限元分析软件，建立砖塔模型的数值模型，考虑材料的非线性、接触非线性等因素，对砖塔在地震作用下的动力响应和损伤演化过程进行数值模拟。将数值模拟结果与振动台试验结果进行对比分析，验证数值模型的准确性和可靠性。通过数值模拟，进一步研究不同因素，如结构形式、材料性能、损伤程度等，对砖塔动力特性和抗震性能的影响，为砖塔的抗震设计和保护加固提供理论依据。地震损伤对砖塔动力特性影响规律总结：综合振动台试验和数值模拟分析结果，总结地震损伤对砖塔动力特性的影响规律，包括自振频率随损伤程度的变化规律、振型的改变特征、阻尼比与损伤的关系等。建立地震损伤与砖塔动力特性参数之间的定量关系模型，为实际砖塔的抗震性能评估提供科学的方法和依据。根据研究结果，提出针对砖塔的抗震保护建议和措施，如合理的加固方法、抗震设计改进建议等，以提高砖塔在地震作用下的安全性和稳定性。1.3.2研究方法试验研究方法：通过小型振动台试验，直接获取砖塔模型在地震作用下的动力响应数据和破坏现象，为研究地震损伤对砖塔动力特性的影响提供第一手资料。试验过程中，严格控制试验条件，确保试验数据的准确性和可靠性。同时，采用多种测量手段和仪器，如加速度传感器、位移传感器、应变片等，对模型的响应进行全面测量，以便深入分析砖塔在地震作用下的力学行为。数值模拟方法：运用有限元分析软件进行数值模拟，能够对砖塔在复杂地震作用下的力学性能进行详细分析，弥补试验研究的局限性。通过建立合理的数值模型，可以模拟不同工况下砖塔的动力响应和损伤演化过程，分析各种因素对砖塔动力特性的影响。数值模拟还可以对试验结果进行验证和补充，为研究提供更全面的信息。对比分析方法：将振动台试验结果与数值模拟结果进行对比分析，相互验证和补充，以提高研究结果的准确性和可靠性。对比分析不同加载工况下砖塔模型的动力响应和动力特性参数的变化，以及不同损伤程度下砖塔的力学性能差异，从而深入揭示地震损伤对砖塔动力特性的影响规律。同时，将研究结果与已有文献资料进行对比，分析异同点，进一步完善研究成果。1.4技术路线本研究的技术路线如图1-1所示，首先广泛收集国内外砖石古塔抗震性能相关研究资料，深入分析现有研究成果与不足，明确研究方向。依据相似理论设计并制作砖塔模型，对模型材料性能和尺寸进行精确测定。制定全面的振动台试验方案，选用多种地震波并进行不同工况加载。在试验过程中，运用传感器实时采集砖塔模型在地震作用下的加速度、位移、应变等响应数据，详细记录裂缝开展和破坏形态。利用环境激励法测试模型加载前后的动力特性参数。对试验数据进行整理和分析，绘制时程曲线，研究动力响应特征和损伤对动力特性的影响规律，采用损伤识别方法评估损伤程度和位置。利用有限元软件建立砖塔数值模型，模拟其在地震作用下的动力响应和损伤演化过程，将模拟结果与试验结果对比验证，分析不同因素对砖塔动力特性和抗震性能的影响。综合试验和模拟结果，总结地震损伤对砖塔动力特性的影响规律，建立定量关系模型，提出抗震保护建议和措施。图1-1技术路线图二、振动台试验模型研制2.1砖石古塔结构形态与地震损伤特征砖石古塔作为我国古代建筑的瑰宝，具有丰富多样的结构形态，主要包括楼阁式、密檐式、亭阁式、覆钵式和金刚宝座式等类型，其中楼阁式和密檐式砖塔较为常见。楼阁式砖塔的建筑形式源自中国传统建筑中的楼阁，其特点是层数较多，一般为奇数层，常见的有五层、七层、九层等。各层之间有明显的塔身和塔檐，塔檐较为宽大，多采用斗拱结构支撑。塔身内部通常设有楼梯，可供人登临远眺，如西安大雁塔，通高64.517米，底层边长25.5米，共七层，采用四方楼阁式结构，体现了唐代建筑的雄伟与精致。楼阁式砖塔的结构较为复杂，其内部空间较大，各层之间的连接方式对结构的整体性和稳定性影响较大。在地震作用下，由于其较高的高度和较大的体量，楼阁式砖塔更容易受到地震力的作用，结构的动力响应较为复杂。密檐式砖塔从楼阁式塔发展而来，与楼阁式塔相比，其塔檐更加密集，层数较多，但内部空间相对较小。密檐式砖塔的第一层塔身特别高，设有门窗，多雕刻佛像或佛经故事，而上部各层塔檐紧密相连，几乎没有实际的楼层空间。许多密檐式砖塔为实心塔，也有的中空但内部没有阶梯，一般不做登临之用，如西安小雁塔，现存13层，高43.30米，南北各层辟有券门，单壁较厚形成筒体结构，整体轮廓呈现为秀丽的卷刹。密檐式砖塔的结构整体性较强，但其塔檐在地震作用下容易受到破坏，导致塔身的稳定性下降。由于其内部空间较小，结构的受力分布相对较为集中，在地震作用下更容易出现应力集中的现象。在地震作用下，砖石古塔容易遭受各种损伤，其损伤形态主要包括裂缝、倾斜、倒塌等。裂缝是砖石古塔最常见的损伤形式之一，可分为水平裂缝、垂直裂缝和斜裂缝等。水平裂缝通常出现在塔身的薄弱部位，如各层塔檐与塔身的连接处、门窗洞口周围等，这是由于地震作用下塔身的水平剪力过大，导致砌体结构的水平抗剪能力不足而产生的。垂直裂缝一般是由于塔身的竖向压力分布不均匀，或者在地震作用下塔身发生扭转，导致砌体结构的竖向承载能力受到破坏而形成的。斜裂缝则多发生在塔身的墙角处或应力集中区域，是由于地震作用下的主拉应力超过了砌体结构的抗拉强度所致。例如，在汶川地震中，许多砖石古塔出现了不同程度的裂缝，部分古塔的裂缝宽度甚至达到了数厘米，严重影响了结构的安全性。倾斜也是砖石古塔在地震作用下常见的损伤现象，主要是由于地震作用下塔身的重心偏移，导致结构的稳定性下降。倾斜可能是由于地基不均匀沉降、塔身结构的不对称性以及地震力的偏心作用等原因引起的。当古塔发生倾斜后，其结构的受力状态将发生改变，进一步加剧了结构的损伤程度。如果倾斜角度过大，古塔可能会发生倒塌，造成严重的破坏。如著名的虎丘塔，由于地基不均匀沉降等原因，塔身逐渐倾斜，成为了“中国的比萨斜塔”。在强烈地震作用下，砖石古塔可能会发生倒塌，这是最严重的损伤形态。倒塌可能是由于塔身结构的整体性被严重破坏，无法承受地震力的作用，或者是由于地基失效，导致塔身失去支撑而倒塌。倒塌不仅会对古塔本身造成毁灭性的破坏，还可能对周围的人员和建筑物造成威胁。在历史上的多次地震中，都有砖石古塔倒塌的案例，如1976年唐山大地震中，许多砖石古塔在地震中倒塌，造成了不可挽回的文化损失。了解砖石古塔的结构形态和地震损伤特征，对于开展振动台试验研究具有重要的指导意义。在试验模型设计和制作过程中，需要充分考虑这些因素，以确保试验模型能够真实反映实际砖石古塔的结构性能和地震响应特征。通过对地震损伤特征的分析，还可以为砖石古塔的抗震性能评估和保护加固提供科学依据。2.2试验模型原型塔的基本特征本试验选取某典型砖塔作为研究对象，该砖塔为楼阁式砖塔，具有较高的历史文化价值和建筑艺术价值，其结构特征在楼阁式砖塔中具有一定的代表性。该砖塔通高X米，共Y层，平面呈正方形。底层边长为Z米，自下而上各层边长以一定比例逐渐递减，形成了稳定而优美的轮廓。塔身整体为砖石结构，采用青砖和石灰砂浆砌筑而成。青砖的尺寸为长a毫米、宽b毫米、厚c毫米，其抗压强度经检测达到XMPa，具有较好的耐久性和力学性能。石灰砂浆由石灰、砂和水按一定比例配制而成，其粘结强度和抗压强度分别为YMPa和ZMPa，保证了砖塔结构的整体性和稳定性。在结构构造方面，砖塔每层均设有塔檐，塔檐采用叠涩出挑的方式，由多层青砖向外悬挑而成，最大悬挑长度为d毫米。塔檐的坡度为i，有利于排水和保护塔身。每层塔檐下设有斗拱，斗拱为砖制，其形式为x，斗拱的设置不仅增加了塔檐的支撑力，还具有装饰性，体现了古代建筑的艺术风格。砖塔内部设有楼梯，楼梯采用砖石结构，宽度为e毫米，坡度为j。楼梯沿着塔身内壁盘旋而上，连接各层，可供人登临。每层内部空间为正方形，边长为f毫米，高度为g毫米，内部空间较为宽敞，可用于供奉佛像或存放物品。该砖塔的基础采用砖石基础，基础埋深为h米，基础底面尺寸为长m米、宽n米。基础采用分层砌筑的方式，每层厚度为k毫米，基础的抗压强度和稳定性能够满足砖塔的承载要求。该砖塔的基本特征为试验模型的设计提供了重要依据，在模型设计过程中，将严格按照相似理论，对砖塔的尺寸、材料性能、结构构造等进行相似模拟，确保试验模型能够真实反映原型砖塔的动力特性和抗震性能。2.3振动台模型试验相似关系相似理论是指导模型试验的重要理论基础，它为确定模型与原型之间的相似程度、等级等提供了依据。在振动台模型试验中，相似理论确保模型能够准确反映原型在地震作用下的动力特性和响应。相似理论主要包含三个重要定理：相似第一定理：两个相似的系统，单值条件相同，其相似判据的数值也相同。这意味着如果模型和原型相似，那么它们在几何形状、材料性质、边界条件和初始条件等单值条件上应保持一致，并且由这些条件导出的相似判据，如无量纲量群，在模型和原型中的数值也相等。相似第一定理规定了现象相似的必要条件，即相似现象能为文字上完全相同的方程组所描述，用来表征这些现象的一切物理量在空间相对应的各点和在时间上相对应的瞬间各自互成一定的比例关系，且各相似常数值要服从于某种自然规律的约束。相似第二定理：当一现象由n个物理量的函数关系来表示，且这些物理量中含有m种基本量纲时，则能得到(n-m)个相似判据。该定理也被称为π定理，它确定了描述相似现象的方程组的解的一般数学结构，即任何定解问题的积分结果都可以表示成准数方程的形式。通过π定理，可以将描述物理现象的众多物理量转换为少数几个独立的相似准数，这些相似准数之间的函数关系反映了物理现象的本质规律。相似第三定理：凡具有同一特性的现象，当单值条件（系统的几何性质、介质的物理性质、起始条件和边界条件等）彼此相似，且由单值条件的物理量所组成的相似判据在数值上相等时，则这些现象必定相似。相似第三定理明确了模型满足什么条件时现象才能相似，它是模型试验所必须遵循的法则，规定了现象相似的充分条件。推导相似比例关系通常采用量纲分析法，以确定模型与原型之间各物理量的相似比。在砖石古塔振动台模型试验中，主要考虑的基本物理量有几何尺寸（L）、弹性模量（E）、密度（ρ）、时间（t）、力（F）等。根据相似理论，首先确定几何相似比C_L。结合试验条件和实际需求，选取几何相似比为1:X。即模型的几何尺寸是原型的1/X，如原型塔的高度为H，模型塔的高度则为H/X。基于几何相似比，可进一步推导其他物理量的相似比：弹性模量相似比：根据相似理论，弹性模量相似比C_E与几何相似比C_L和应力相似比C_σ相关。在振动台试验中，通常假设模型与原型材料的应力应变关系相似，即应力相似比等于弹性模量相似比，所以C_E=C_σ。密度相似比：考虑模型与原型在重力和惯性力作用下的相似性，密度相似比C_ρ与几何相似比C_L和加速度相似比C_a有关。在地震作用下，加速度相似比通常取为1，即模型与原型在地震作用下的加速度响应相同，所以C_ρ=C_σ/(C_LC_a)。时间相似比：时间相似比C_t与几何相似比C_L和速度相似比C_v相关。速度相似比又与加速度相似比和时间相似比有关，通过推导可得C_t=\sqrt{C_L/C_a}。力相似比：力相似比C_F与质量相似比C_m和加速度相似比C_a相关，而质量相似比又与密度相似比和几何相似比有关，经推导C_F=C_ρC_L^3C_a。通过上述相似比例关系的推导，建立了模型与原型之间物理量的对应关系，为振动台模型试验的设计和实施提供了重要依据。在模型制作和试验过程中，严格按照这些相似比来确定模型的材料、尺寸以及加载制度等参数，能够保证模型试验结果能够准确反映原型砖塔在地震作用下的动力特性和响应。2.4模型设计在进行砖塔模型设计时，为简化制作过程，在保持主要结构特征和力学性能相似的前提下，对原型塔进行了一定程度的简化。去除了一些对整体动力特性影响较小的装饰性构件，如塔身上的小型雕刻、局部彩绘等，以减少模型制作的复杂性，同时确保这些简化不会对模型的动力特性和地震响应产生显著影响。在模型材料选择方面，遵循相似理论和试验要求，力求选择与原型砖塔材料性能相似且适合模型制作的材料。对于砖块材料，考虑到原型砖塔采用青砖，经过调研和试验，选用了一种特制的小型青砖，其尺寸按照几何相似比进行缩放，抗压强度和弹性模量等力学性能指标与原型青砖相似。通过对多批次砖块的抽样检测，得到该小型青砖的平均抗压强度为XMPa，弹性模量为YMPa。砌筑模型所用的砂浆采用了由石灰、砂和少量添加剂按特定比例配制而成的材料。其中，石灰与砂的比例经过多次试验优化，以确保砂浆的粘结强度和抗压强度与原型塔的砂浆性能相似。添加剂的加入旨在改善砂浆的工作性能和耐久性。经测试，该砂浆的粘结强度达到ZMPa，抗压强度为WMPa。为了准确模拟砖塔的结构，对砖块和砌筑模具进行了精心设计。根据相似比确定砖块的尺寸为长a毫米、宽b毫米、厚c毫米，确保砖块之间的连接和受力方式与原型塔相似。砌筑模具采用定制的钢模具，具有高精度的尺寸和良好的脱模性能，能够保证砖块的成型质量和砌筑精度。在砌筑过程中，严格控制砂浆的厚度和饱满度，使其均匀分布在砖块之间，确保模型的整体性和力学性能。在模型设计过程中，充分考虑了各种因素对模型动力特性的影响，通过合理选择材料和设计模具，确保了模型能够真实反映原型砖塔的力学性能和地震响应特征。2.5试验模型制备在进行模型材料制备时，砖块原料的混合与处理是关键步骤。将选取的优质黏土、适量的页岩以及特定添加剂按比例混合，以确保砖块具有合适的物理性能。在混合过程中，使用专业的搅拌设备，充分搅拌均匀，使各种原料充分融合。随后，对混合原料进行陈化处理，将其放置在特定环境中一段时间，让水分均匀分布，提高原料的可塑性和稳定性。采用定制的钢模具进行砖块成型，将处理好的原料放入模具中，通过压力机施加一定压力，使砖块成型。成型后的砖块进行自然养护或蒸汽养护，使其强度逐渐增长。在养护过程中，严格控制养护温度和湿度，确保砖块的质量。砂浆制备方面，按照设计比例将石灰、砂和添加剂混合，加入适量的水进行搅拌。搅拌过程中，控制搅拌速度和时间，保证砂浆的均匀性和和易性。搅拌完成后，对砂浆的性能进行检测，如粘结强度、抗压强度等，确保其符合试验要求。模型制作流程如下：首先搭建模型制作平台，确保平台水平、稳定。在平台上铺设防潮层，防止水分对模型造成影响。然后，根据设计尺寸在平台上放线，确定模型的位置和尺寸。从模型底部开始，按照设计要求进行底层砖块的砌筑。在砌筑过程中，使用水平仪和铅垂线确保砖块的水平度和垂直度。在砖块之间均匀涂抹砂浆，厚度控制在合理范围内，确保砂浆饱满，使砖块之间的连接牢固。每砌筑一层，都要对模型的尺寸、垂直度进行检查和调整，确保模型的准确性。在模型的关键部位，如墙角、门窗洞口等，采取加强措施，如增加砖块数量、使用特殊的连接方式等，以提高模型的整体稳定性。随着模型的逐层砌筑，在适当位置安装预先准备好的传感器预埋件，确保预埋件位置准确，固定牢固。继续向上砌筑，直至完成整个模型的主体结构。模型主体完成后，进行塔檐和斗拱等附属结构的安装。塔檐采用预制构件，按照设计位置和角度进行安装，使用砂浆或其他连接材料确保连接牢固。斗拱结构较为复杂，先在地面进行组装调试，然后整体安装到模型上，保证斗拱的位置和角度准确，与主体结构连接紧密。完成模型的主体和附属结构施工后，对模型进行整体检查和修整。检查模型的外观是否符合设计要求，对表面不平整、砂浆不饱满等问题进行处理。对模型的内部结构进行检查，确保内部连接牢固，无松动现象。最后，对模型进行编号和标识，记录模型的相关信息，便于后续试验和数据分析。模型制备流程如图2-1所示。图2-1模型制备流程图三、砖塔模型的振动台试验3.1试验设备本次试验选用了型号为XX的电动式振动台，该振动台具有高精度、高可靠性和宽频带的特点，能够满足模拟多种地震波的试验需求。其主要参数如下：最大承载重量为X吨，足以支撑砖塔模型及相关附属设备的重量；频率范围为0.1-100Hz，可覆盖常见地震波的频率成分；最大位移为±Ymm，最大加速度为±Zg，能够模拟不同强度的地震激励。电动式振动台的工作原理基于电磁感应定律。其主要由振动台体、功率放大器、信号发生器和控制系统等部分组成。信号发生器产生模拟地震波的电信号，该信号经过功率放大器放大后，输入到振动台体的励磁线圈中。励磁线圈在磁场中受到电磁力的作用，产生与输入电信号相对应的机械振动。振动台体通过台面将振动传递给放置在其上的砖塔模型，从而实现对地震作用的模拟。在振动过程中，控制系统实时监测和调整振动台的输出参数，确保试验的准确性和稳定性。为了全面测量砖塔模型在地震作用下的动力响应，试验采用了多个加速度传感器。这些传感器型号为XX，具有高精度、高灵敏度和宽频响的性能特点。其测量范围为±Xg，能够满足不同地震强度下的加速度测量需求；灵敏度为YmV/g，分辨率高，能够精确测量微小的加速度变化；频率响应范围为0.5-1000Hz，可准确捕捉模型在地震作用下的高频振动响应。加速度传感器基于压电效应工作，当传感器受到加速度作用时，内部的压电材料会产生电荷，电荷的大小与加速度成正比。通过测量电荷的大小，即可得到模型的加速度响应。高速摄像机在试验中用于记录砖塔模型的变形和破坏过程。选用的高速摄像机型号为XX，帧率可达Xfps，能够清晰捕捉模型在地震作用下的瞬间变形和裂缝开展情况；分辨率为Y×Z像素，可提供高清晰度的图像，便于后续对模型的破坏形态进行分析。在试验过程中，将高速摄像机固定在合适的位置，使其能够完整拍摄到砖塔模型的各个部位。通过高速摄像机拍摄的视频，可直观地观察模型在不同地震作用下的响应特征，为研究模型的破坏机制提供重要依据。试验还配备了专业的数据采集与分析软件，如XX软件。该软件具备强大的数据采集功能，能够同时采集多个传感器的数据，并以高速、高精度的方式进行存储。在数据分析方面，软件提供了丰富的工具和算法，可对采集到的数据进行滤波、积分、频谱分析等处理。通过该软件，能够快速准确地得到砖塔模型的加速度时程曲线、速度时程曲线、位移时程曲线以及频率响应等参数，为深入分析模型的动力特性和地震响应提供了便利。3.2模型的定位和组装在将砖塔模型放置于振动台台面之前，需要对振动台台面进行清洁和检查，确保台面上无杂物、无损伤，避免影响模型的放置和试验结果。使用水平仪对台面进行水平度检测，若台面不水平，需通过调整振动台的支撑脚来使其达到水平状态，保证模型在振动过程中受力均匀。依据预先设计的试验方案，在振动台台面上精确确定砖塔模型的放置位置。采用全站仪等测量仪器，按照坐标定位的方法，在台面上标记出模型的中心位置以及各边缘的位置。定位过程中，需考虑振动台的振动方向和模型的受力特点，确保模型的主要受力方向与振动台的振动方向一致，以准确模拟地震作用下的受力情况。例如，如果试验主要研究砖塔在水平方向的地震响应，应将模型的长轴方向与振动台的水平振动方向平行放置。为了确保砖塔模型在振动过程中稳固，不发生位移或晃动，采用了特制的固定装置。在模型底部与振动台台面之间设置高强度的螺栓连接，通过在模型底部预留的螺栓孔和台面上对应的安装孔，使用螺栓将模型与台面紧密固定。同时，在模型的周围增设支撑框架，支撑框架采用型钢制作，具有足够的强度和刚度。支撑框架与模型之间通过橡胶垫等柔性材料连接，既能提供稳定的支撑，又能避免对模型造成损伤。在固定过程中，使用扭矩扳手严格控制螺栓的拧紧力矩，确保每个螺栓的紧固程度一致，使模型在振动过程中能够与振动台台面协同运动。模型的组装从基础部分开始，按照由下而上的顺序进行。在组装过程中，严格遵循设计图纸和施工规范，确保每个构件的位置准确无误。对于砖块之间的连接，使用设计好的砂浆进行砌筑，保证砂浆的饱满度和粘结强度。在砌筑过程中，使用水平尺和铅垂线随时检查砖块的水平度和垂直度，偏差控制在允许范围内。对于塔檐、斗拱等附属结构，先在地面进行预组装，检查各构件之间的配合精度，确保无误后再进行整体安装。安装时，使用专用的连接件将附属结构与主体结构牢固连接，如采用榫卯连接或螺栓连接等方式。在组装过程中，避免对模型造成碰撞和损伤，防止影响模型的结构性能。在模型组装完成后，对整个模型进行全面检查，包括结构的完整性、构件的连接牢固性、各部分尺寸的准确性等。对模型的外观进行检查，查看是否有砖块松动、砂浆开裂等现象，如有问题及时进行修复和处理。通过全面细致的检查和调试，确保模型的组装质量满足试验要求，为后续的振动台试验提供可靠的保障。3.3测点布置为全面、准确地获取砖塔模型在地震作用下的动力响应信息，依据研究目的和砖塔的结构特点，在关键部位合理布置测点。测点布置主要集中在塔基、塔身和塔顶等部位，旨在监测模型在地震过程中的加速度、位移和应变等参数，为后续的数据分析和损伤评估提供丰富的数据支持。在塔基部位，布置加速度传感器和应变片。加速度传感器用于测量塔基在地震作用下的加速度响应，了解地震波输入的特性和能量传递情况。应变片则用于监测塔基在地震作用下的应力应变状态，判断塔基是否出现开裂、破损等损伤情况。塔基作为砖塔与基础的连接部位，承受着塔身传来的全部荷载，在地震作用下，塔基的受力状态直接影响着砖塔的整体稳定性。通过在塔基布置测点，可以及时捕捉到塔基的动力响应和损伤特征，为分析砖塔的地震破坏机制提供重要依据。塔身是砖塔的主要承重结构，也是地震作用下最容易出现损伤的部位。在塔身的不同高度和不同方位布置多个加速度传感器和位移计。加速度传感器可以测量塔身不同部位在地震作用下的加速度响应，分析加速度沿塔身高度的分布规律，了解地震作用下塔身的振动特性。位移计则用于测量塔身的水平位移和竖向位移，研究塔身的变形情况。在塔身的薄弱部位，如门窗洞口周围、各层塔檐与塔身的连接处等，增设应变片，监测这些部位的应力应变变化，判断是否出现裂缝等损伤。由于塔身较高，在地震作用下，塔身不同部位的动力响应和受力状态存在差异，通过在塔身合理布置测点，可以全面掌握塔身的动力特性和损伤演化过程。塔顶作为砖塔的最高部位，在地震作用下的位移和加速度响应较大，对砖塔的整体稳定性影响也较大。在塔顶布置加速度传感器和位移计，测量塔顶在地震作用下的加速度和位移响应。塔顶的动力响应能够反映砖塔整体的振动特性和地震作用的强烈程度。通过监测塔顶的加速度和位移，可以评估砖塔在地震作用下的稳定性，为砖塔的抗震性能评估提供重要数据。测点布置的依据是基于对砖塔结构力学性能和地震响应特征的分析。通过对砖塔在地震作用下的受力分析和振动理论研究，确定了塔基、塔身和塔顶等关键部位在地震响应中的重要性。这些部位的动力响应能够直接反映砖塔的整体性能和损伤情况。同时，考虑到传感器的测量精度和可靠性，以及试验操作的便利性，合理选择了传感器的类型和布置位置。例如，加速度传感器选择灵敏度高、频率响应范围宽的型号，以准确测量地震作用下的加速度变化；位移计选择精度高、量程合适的型号，以满足塔身和塔顶位移测量的需求。应变片则根据不同部位的受力特点，选择合适的规格和粘贴位置，确保能够准确测量应力应变的变化。测点布置的作用在于能够实时、准确地获取砖塔模型在地震作用下的动力响应数据。这些数据是研究地震损伤对砖塔动力特性影响的基础，通过对数据的分析，可以深入了解砖塔在地震过程中的力学行为，揭示地震损伤对砖塔自振频率、振型和阻尼比等动力特性参数的影响规律。同时，根据测点数据可以判断砖塔模型的损伤程度和损伤位置，为砖塔的抗震性能评估和保护加固提供科学依据。3.4地震波选择与处理在振动台试验中，地震波的选择对试验结果有着至关重要的影响。合适的地震波应能够真实地模拟实际地震的特性，从而准确地反映砖塔在地震作用下的动力响应。本次试验选用了El-Centro波作为输入地震波。El-Centro波是1940年美国加利福尼亚州ImperialValley地震时在El-Centro台站记录到的地震波，该地震的震级为7.1级，属于中强地震。El-Centro波具有丰富的频率成分，涵盖了从低频到高频的多个频段，能够较好地模拟不同类型地震的频谱特性。其加速度时程曲线具有明显的峰值和复杂的波动特征，包含了多个强烈的震动周期，能够对砖塔模型产生较为全面的地震作用。许多学者在相关结构的振动台试验研究中，也常常选用El-Centro波作为输入地震波，其试验结果具有较好的代表性和参考价值。因此，选择El-Centro波能够为本次砖塔模型的振动台试验提供可靠的地震激励。在选用El-Centro波后，需要对其进行必要的处理，以满足试验要求。首先是滤波处理，由于实际记录的地震波中可能包含噪声和其他干扰信号，这些信号会影响试验结果的准确性。采用带通滤波器对El-Centro波进行滤波处理，去除噪声和低频干扰信号，保留有效频率范围内的地震波信号。通过滤波处理，可以提高地震波信号的质量，使试验结果更加准确可靠。其次是调幅处理，根据试验方案中设定的不同地震烈度，对El-Centro波的峰值加速度进行调整。通过调幅处理，可以模拟不同强度的地震作用，研究砖塔模型在不同地震烈度下的动力响应和损伤演化规律。在调幅过程中，依据相似理论，确保模型所受的地震作用与原型在相应地震烈度下的受力情况相似。例如，将El-Centro波的峰值加速度分别调整为0.1g、0.2g、0.4g等，以模拟不同地震烈度下的地震作用。通过对地震波的调幅处理，可以系统地研究砖塔模型在不同地震强度下的抗震性能，为砖塔的抗震保护提供更全面的依据。3.5电磁式振动台非线性迭代补偿反馈在砖塔模型的振动台试验中，电磁式振动台的非线性问题对试验结果的准确性有着显著影响。振动台系统中的非线性因素较为复杂，主要包括台面运动过程中的摩擦非线性、功率放大器的饱和非线性以及传感器的非线性特性等。这些非线性因素会导致振动台输出的实际振动信号与输入的控制信号之间存在偏差，使得试验结果无法真实反映砖塔在理想地震作用下的动力响应。为解决振动台的非线性问题，采用非线性迭代补偿反馈方法。该方法的基本原理是基于振动台系统的逆动力学模型，通过对系统非线性特性的分析和建模，实时计算出系统的非线性误差，并将该误差作为补偿信号反馈到控制系统中，对输入信号进行修正，从而实现对振动台输出的精确控制。具体而言，在每次振动试验过程中，控制系统首先根据设定的地震波信号向振动台发送初始控制信号。传感器实时采集振动台台面的实际振动响应数据，包括加速度、位移等信息，并将这些数据反馈给控制系统。控制系统利用预先建立的非线性模型，对传感器反馈的数据进行分析和处理，计算出由于系统非线性导致的实际振动响应与期望振动响应之间的误差。根据计算得到的误差，控制系统生成相应的补偿信号，该补偿信号包含了对非线性因素的修正信息。将补偿信号叠加到下一次的输入控制信号中，对输入信号进行调整和补偿。通过这样不断地迭代计算和补偿，逐渐减小振动台输出的误差，使振动台的实际振动响应尽可能接近期望的地震波信号。以某一次试验为例，在试验初期，由于振动台系统的非线性，实际输出的加速度信号与期望的El-Centro波加速度信号存在较大偏差，加速度峰值误差达到了X%。经过第一次迭代补偿后，误差减小到了Y%。随着迭代次数的增加，误差逐渐收敛，经过N次迭代后，加速度峰值误差降低到了Z%，有效地提高了振动控制精度，使试验结果更加准确可靠。通过非线性迭代补偿反馈方法，能够显著提高振动台的振动控制精度，减小系统非线性对试验结果的影响。这为砖塔模型在振动台试验中获取准确的动力响应数据提供了有力保障，有助于更精确地研究地震损伤对砖塔动力特性的影响。3.6输入期望信号的频谱分析对经过滤波和调幅处理后的El-Centro波进行频谱分析，有助于深入了解其频率成分和幅值分布情况，为砖塔模型振动台试验提供更为精准的参考依据。运用快速傅里叶变换（FFT）算法对处理后的地震波进行频谱分析。该算法能够高效地将时域信号转换为频域信号，从而清晰地呈现出地震波在不同频率下的幅值分布。在MATLAB软件平台上，通过调用相关函数，对处理后的El-Centro波加速度时程数据进行快速傅里叶变换，得到其频谱图，如图3-1所示。图3-1El-Centro波频谱图从频谱图中可以看出，El-Centro波的频率成分丰富，涵盖了从低频到高频的多个频段。在低频段，0-5Hz范围内，幅值相对较大，表明该频段的能量较为集中。其中，在1Hz左右出现了一个明显的幅值峰值，这意味着在这个频率附近，地震波的能量较为突出。低频成分主要与地震波的长周期振动特性相关，它对砖塔结构的整体振动响应和大变形起着重要作用。在地震作用下，砖塔的整体晃动和较大幅度的位移变形往往与低频成分的激励密切相关。在中高频段，5-20Hz范围内，虽然幅值相对低频段有所减小，但频率成分更为复杂。该频段包含了多个幅值较小的峰值，表明地震波在这个频率区间内具有多种不同频率的振动成分。中高频成分主要影响砖塔结构的局部振动响应，如塔檐、斗拱等附属结构的振动以及塔身局部的应力集中区域的响应。这些部位在中高频地震波的作用下，更容易出现局部的开裂、破损等损伤。通过对El-Centro波频谱分析结果的研究，可以为砖塔模型振动台试验提供以下重要参考：在试验过程中，根据地震波的主要频率成分，合理设置传感器的频率响应范围，确保能够准确捕捉到砖塔模型在不同频率下的振动响应。例如，对于主要关注砖塔整体振动响应的测点，选择频率响应范围能够覆盖低频段的传感器；而对于监测塔檐、斗拱等局部结构响应的测点，则需要选择频率响应范围更宽，能够涵盖中高频段的传感器。同时，根据频谱分析结果，结合砖塔的结构特点和固有频率，分析地震波与砖塔结构之间的共振可能性。如果地震波的某些频率成分与砖塔的固有频率相近，可能会引发共振现象，导致砖塔结构的振动响应显著增大，从而加剧结构的损伤。通过这种分析，可以提前采取措施，如调整地震波的输入方式或对砖塔结构进行适当的加固，以避免共振的发生。频谱分析结果还可以为后续的数值模拟分析提供参考，在建立砖塔数值模型时，根据地震波的频率特性和幅值分布，合理设置模型的参数和边界条件，使数值模拟结果更加准确地反映实际情况。3.7采样频率与同步系数的对应关系在砖塔模型振动台试验的数据采集过程中，采样频率的选择至关重要，它直接影响到所采集数据的准确性和完整性，进而影响对砖塔动力特性的分析结果。采样频率是指单位时间内采集数据的次数，单位为赫兹（Hz）。根据奈奎斯特采样定理，为了能够准确地还原原始信号，采样频率应至少为信号最高频率的两倍。在砖塔模型振动台试验中，砖塔在地震作用下的振动响应包含了丰富的频率成分。通过对地震波的频谱分析以及砖塔结构动力学理论分析可知，砖塔振动响应的频率范围较宽，其中高频成分主要与砖塔的局部振动相关，如塔檐、斗拱等部位的振动；低频成分则主要与砖塔的整体振动相关。为了全面准确地捕捉砖塔在地震作用下的振动响应，需要合理选择采样频率。在本试验中，经过前期的理论分析和预试验，确定了采样频率为1000Hz。选择这一采样频率主要基于以下考虑：一方面，1000Hz的采样频率能够满足奈奎斯特采样定理对于砖塔振动响应信号的要求，确保能够准确捕捉到信号中的高频成分，从而精确还原砖塔在地震作用下的真实振动情况。另一方面，考虑到试验数据采集系统的性能和存储能力，1000Hz的采样频率既能够保证采集到足够详细的数据，又不会给数据存储和后续处理带来过大的负担。同步系数是指在多通道数据采集过程中，各通道数据之间的时间同步精度。在砖塔模型振动台试验中，由于需要同时采集多个测点的加速度、位移、应变等数据，各通道数据的同步性对于准确分析砖塔的动力特性至关重要。如果各通道数据不同步，将会导致数据分析结果出现偏差，无法准确反映砖塔在地震作用下的真实受力状态和变形情况。为了确保各通道数据的同步性，在试验中采取了一系列措施。选用了具有高精度同步功能的数据采集系统，该系统能够保证各通道数据在采集过程中的时间同步误差控制在极小的范围内。在试验前，对数据采集系统进行了严格的校准和调试，确保各通道的采样起始时间一致。同时，在试验过程中，通过同步触发信号来启动各通道的数据采集，进一步提高了数据的同步精度。通过实际试验和数据分析，研究了采样频率与同步系数之间的对应关系。在不同的采样频率下，对同步系数进行了测试和分析。结果表明，随着采样频率的提高，同步系数的精度也相应提高。当采样频率较低时，由于数据采集的时间间隔较大，同步系数的误差相对较大；而当采样频率提高到一定程度后，同步系数的误差明显减小，各通道数据的同步性得到了显著改善。在本试验选择的1000Hz采样频率下，同步系数达到了0.99以上，满足了试验对于数据同步性的要求。这意味着在该采样频率下，各通道采集的数据能够在时间上保持高度的一致性，为准确分析砖塔的动力特性提供了可靠的数据基础。3.8试验测试方法与加载方案在动力特性测试方法的选择上，本试验采用自由振动法。自由振动法是一种较为常用且有效的动力特性测试方法，其原理基于结构动力学理论。当砖塔模型在外界激励下产生自由振动时，其振动过程遵循一定的动力学方程。根据牛顿第二定律，结构的振动方程可以表示为M\ddot{u}+C\dot{u}+Ku=0，其中M为质量矩阵，C为阻尼矩阵，K为刚度矩阵，\ddot{u}、\dot{u}、u分别为加速度、速度和位移响应。在自由振动状态下，模型仅在自身惯性力和阻尼力的作用下振动，通过测量模型在自由振动过程中的响应，如加速度、位移等，就可以分析得到其动力特性参数。在试验中，首先对砖塔模型施加一个初始激励，使其产生自由振动。可以采用多种方式施加初始激励，如使用小型锤击装置对模型进行瞬间敲击，或者利用电磁激振器对模型施加一个短暂的冲击力。在模型产生自由振动后，利用布置在关键部位的加速度传感器和位移传感器，实时采集模型的振动响应数据。通过对采集到的数据进行处理和分析，如采用快速傅里叶变换（FFT）等算法，可以得到模型的自振频率、振型等动力特性参数。自由振动法的优点在于操作相对简单，能够较为直观地反映结构的固有动力特性，而且不需要复杂的加载设备和试验条件。同时，该方法对模型的损伤较小，不会对后续的地震加载试验产生较大影响。加载制度是振动台试验的关键环节，它直接影响到试验结果的准确性和可靠性。本次试验采用多工况加载制度，旨在全面研究砖塔模型在不同地震强度作用下的动力响应和损伤演化规律。根据《建筑抗震设计规范》（GB50011-2010）以及相关研究成果，确定了不同工况下的加载参数。试验共设置了X个工况，每个工况对应不同的地震烈度和峰值加速度。具体加载工况如表3-1所示。工况编号地震波峰值加速度（g）对应地震烈度1El-Centro波0.1Ⅶ度2El-Centro波0.2Ⅷ度3El-Centro波0.4Ⅸ度4El-Centro波0.6Ⅹ度5El-Centro波0.8Ⅺ度6El-Centro波1.0Ⅻ度表3-1加载工况表在每个工况加载前，先对砖塔模型进行环境激励下的动力特性测试，作为初始状态数据。然后按照设定的峰值加速度对El-Centro波进行调幅，将调幅后的地震波输入振动台，对砖塔模型进行加载。在加载过程中，利用传感器实时采集模型的加速度、位移、应变等响应数据，并通过高速摄像机记录模型的裂缝开展、破坏形态等现象。每个工况加载完成后，再次对模型进行环境激励下的动力特性测试，对比加载前后的动力特性参数，分析地震损伤对砖塔动力特性的影响。相邻工况之间设置一定的时间间隔，以便对模型进行检查和数据整理。同时，在整个试验过程中，密切关注模型的状态，确保试验的安全进行。通过这种多工况加载制度，可以系统地研究砖塔模型在不同地震强度下的动力响应和损伤演化规律，为深入了解地震损伤对砖塔动力特性的影响提供丰富的数据支持。3.9试验过程与试验现象在完成砖塔模型的定位、组装以及测点布置等准备工作后，正式启动振动台试验。首先开启振动台控制系统，对系统进行初始化设置，确保振动台处于正常工作状态。按照预先设定的加载方案，依次输入不同工况的地震波信号。在输入地震波前，再次检查各传感器与数据采集系统的连接是否正常，确保能够准确采集试验数据。在加载过程中，密切关注振动台的运行状态和砖塔模型的响应。随着地震波的输入，振动台台面开始按照设定的加速度和频率进行振动，砖塔模型在振动台的激励下产生动力响应。传感器实时采集砖塔模型在地震作用下的加速度、位移和应变等数据，并将数据传输至数据采集系统进行存储和初步处理。同时，高速摄像机全程记录砖塔模型的变形和破坏过程，以便后续对试验现象进行详细分析。在工况1，输入峰值加速度为0.1g的El-Centro波时，砖塔模型在地震作用下整体响应较小。通过传感器数据监测和高速摄像机观察，发现模型无明显裂缝产生，仅在塔顶部位出现了微小的位移，位移量约为Xmm。模型的加速度响应也相对较小，塔顶的最大加速度为Yg。此时，砖塔模型基本处于弹性阶段，结构性能未受到明显损伤。当进入工况2，峰值加速度提高到0.2g时，砖塔模型的动力响应有所增强。在模型的底部和塔身中部，开始出现一些细微的裂缝。这些裂缝主要为水平裂缝，宽度较窄，约为0.1-0.2mm。塔顶的位移明显增大，达到了Zmm，加速度也相应增加，塔顶最大加速度达到了Wg。模型的自振频率和阻尼比等动力特性参数开始发生变化，自振频率略有降低，阻尼比稍有增加。随着工况3，峰值加速度提升至0.4g，砖塔模型的损伤进一步发展。塔身的裂缝数量增多，宽度增大，部分裂缝开始贯通。在塔檐与塔身的连接处，出现了较为明显的裂缝，导致塔檐的局部稳定性受到影响。塔顶的位移继续增大，达到了Amm，加速度峰值也提高到了Bg。此时，模型的自振频率显著降低，阻尼比明显增大，表明模型的结构损伤程度在加剧。在工况4，峰值加速度为0.6g时，砖塔模型的破坏现象更加严重。塔身出现了多条贯通裂缝，部分砖块开始松动，塔体的整体性受到较大影响。塔顶的位移急剧增大，达到了Cmm，模型出现了明显的倾斜，倾斜角度约为D°。加速度响应也呈现出较大的波动，表明模型的结构刚度在迅速下降。当工况5，峰值加速度达到0.8g时，砖塔模型的损伤已非常严重。塔身的裂缝进一步扩展，部分区域的砖块脱落，塔体出现了较大的变形。塔顶的位移持续增大，达到了Emm，模型的倾斜角度进一步加大，约为F°。此时，模型的自振频率大幅降低，阻尼比急剧增大，结构接近破坏状态。在最后一个工况6，峰值加速度为1.0g时，砖塔模型最终发生倒塌。塔身大部分砖块散落，塔体完全失去承载能力。通过对试验过程和现象的详细记录与分析，可以清晰地了解砖塔模型在不同地震强度作用下的损伤演化过程和破坏机制，为后续的试验数据分析和理论研究提供了重要依据。3.10模型试验分析对砖塔模型在振动台试验中采集到的数据进行深入分析，能够揭示地震损伤对砖塔动力特性的影响规律，为评估砖塔的抗震性能提供重要依据。通过对试验数据的处理，计算得到砖塔模型在不同工况下的自振频率和阻尼比。自振频率是结构的固有属性，它反映了结构在自由振动状态下的振动快慢。阻尼比则表征了结构在振动过程中能量耗散的程度。在试验初期，工况1加载前，砖塔模型处于初始完好状态，其自振频率经计算为f_1Hz，阻尼比为\xi_1。随着地震作用的增强，在工况3加载后，模型出现了一定程度的损伤，此时自振频率降低至f_3Hz，阻尼比增大至\xi_3。当工况6加载完成，模型接近倒塌状态时，自振频率进一步下降到f_6Hz，阻尼比大幅增加到\xi_6。通过对不同工况下自振频率和阻尼比的对比分析，可以看出随着地震损伤的加剧，砖塔模型的自振频率逐渐降低，阻尼比逐渐增大。这是因为地震损伤导致砖塔结构的刚度下降，质量分布发生改变，从而使自振频率降低；同时，损伤的出现增加了结构内部的能量耗散机制，使得阻尼比增大。自振频率和阻尼比的变化能够直观地反映出砖塔结构的损伤程度，可作为评估砖塔抗震性能的重要指标。对砖塔模型在不同工况下的峰值加速度、速度和位移响应进行分析，有助于了解模型在地震作用下的动力响应特征。在不同工况下，砖塔模型不同部位的峰值加速度响应呈现出明显的变化规律。随着地震烈度的增加，即从工况1到工况6，各部位的峰值加速度逐渐增大。以塔顶部位为例，工况1时塔顶峰值加速度为a_{1t}g，到工况6时增大至a_{6t}g。同时，峰值加速度沿塔身高度的分布也发生了变化。在低烈度地震作用下，峰值加速度在塔身上部和下部的差异相对较小；而在高烈度地震作用下，塔顶部的峰值加速度明显大于底部，这表明在强烈地震作用下，塔顶部的动力响应更为剧烈，更容易受到损伤。速度响应方面，随着地震作用的增强，砖塔模型各部位的速度响应也逐渐增大。工况2时，塔身中部的峰值速度为v_{2m}m/s，到工况5时增大至v_{5m}m/s。速度响应的变化反映了模型在地震过程中的运动状态变化，较大的速度响应意味着模型在地震作用下的动能增加，结构受到的冲击力也相应增大。位移响应是衡量砖塔结构变形的重要指标。随着地震烈度的提高，砖塔模型的位移响应显著增大。工况3时，塔顶的水平位移为d_{3t}mm，而在工况6时，塔顶水平位移急剧增大到d_{6t}mm，且模型出现了明显的倾斜。通过对位移响应的分析，可以直观地了解砖塔模型在地震作用下的变形情况，位移的增大表明结构的损伤在不断积累，结构的稳定性逐渐降低。动力放大系数是结构在动力荷载作用下的响应与静态荷载作用下响应的比值，它反映了结构在地震作用下的动力放大效应。计算砖塔模型在不同工况下各测点的动力放大系数，发现随着地震烈度的增加，动力放大系数呈现出先增大后减小的趋势。在低烈度地震作用下，结构基本处于弹性阶段，动力放大系数相对较大。随着地震烈度的提高，结构逐渐进入非线性阶段，损伤不断发展，结构的刚度降低，动力放大系数逐渐减小。例如，在工况1时，某测点的动力放大系数为β_1，到工况4时增大至β_4，而在工况6时减小到β_6。动力放大系数的变化规律表明，在地震作用过程中，结构的动力响应并非简单地随地震力的增大而线性增大，而是受到结构自身损伤和非线性特性的影响。层间位移角是衡量结构层间变形的重要参数，它反映了结构在水平地震作用下各楼层之间的相对变形程度。根据试验数据计算砖塔模型各层的层间位移角，随着地震烈度的增加，各层的层间位移角逐渐增大。在工况2时，某层的层间位移角为\theta_2，到工况5时增大至\theta_5。当层间位移角超过一定限值时，结构可能会发生严重破坏甚至倒塌。通过对层间位移角的分析，可以评估砖塔结构在不同地震烈度下的安全性能，判断结构是否满足抗震设计要求。在本试验中，当工况6加载时，部分楼层的层间位移角超过了相关规范规定的限值，表明此时砖塔模型的结构已处于危险状态。通过对砖塔模型在振动台试验中的自振频率、阻尼比、峰值加速度、速度、位移响应、动力放大系数和层间位移角等参数的分析，可以全面了解地震损伤对砖塔动力特性的影响规律，为评估砖塔的抗震性能提供了有力的数据支持和理论依据。四、基于ABAQUS的数值模拟分析4.1ABAQUS有限单元法选取与简介ABAQUS是一款功能强大的工程模拟有限元软件，在结构力学、热传导、流体力学等众多领域都有广泛应用。它具备丰富的单元库，包含多种类型的单元，如实体单元、壳单元、梁单元等，能够满足不同结构形式的建模需求。在材料模型方面，ABAQUS提供了广泛的选择，涵盖线性弹性、非线性弹性、塑性、粘弹性、粘塑性等多种材料行为的模拟。这使得用户可以根据实际材料特性，准确地定义模型的材料属性，从而更真实地模拟结构在各种工况下的力学响应。ABAQUS采用的有限单元法，是一种将连续体离散化为有限个单元进行数值求解的方法。其基本原理是将复杂的连续求解区域划分成有限个单元，这些单元通过节点相互连接。在每个单元上，假设一个简单的函数（插值函数）来近似表示其位移分布规律。以二维问题为例，对于一个平面应力问题，将求解区域离散为三角形或四边形单元，在三角形单元中，通常假设位移函数为线性函数，通过单元节点的位移来确定整个单元的位移分布。利用弹塑性理论中的变分原理或其他方法，建立单元节点的力和位移之间的力学特性关系，得到一组以节点位移为未知量的代数方程组。通过求解这些方程组，可以得到节点的位移分量，进而利用插值函数确定单元集合体上的场函数，如应力、应变等。在求解过程中，ABAQUS会自动考虑单元之间的连接和相互作用，确保整个模型的力学行为得到准确模拟。选择ABAQUS软件进行砖塔在地震作用下的动力响应数值模拟，具有多方面的优势。ABAQUS强大的非线性分析能力能够准确模拟砖塔在地震作用下的复杂力学行为。地震作用下，砖塔结构会进入非线性阶段，出现材料非线性（如砖和砂浆的塑性变形）和几何非线性（如结构的大变形）。ABAQUS能够精确处理这些非线性问题，通过合理设置材料本构模型和非线性求解算法，准确模拟砖塔在地震过程中的损伤演化和破坏过程。ABAQUS拥有丰富的接触算法和界面模型。砖塔由砖块和砂浆组成，砖块之间通过砂浆粘结，在地震作用下，砖块与砂浆之间可能会出现脱粘、滑移等接触非线性行为。ABAQUS可以通过定义合适的接触对和接触属性，准确模拟这种接触非线性，考虑砖块与砂浆之间的相互作用，从而更真实地反映砖塔结构的力学性能。ABAQUS具备高效的计算性能和良好的并行计算能力。对于大型复杂的砖塔模型，计算量通常较大，ABAQUS能够充分利用计算机资源，通过并行计算技术加速计算过程，缩短计算时间。这使得在进行多工况、多参数的数值模拟分析时，能够快速得到计算结果，提高研究效率。ABAQUS还拥有友好的用户界面和强大的后处理功能。用户可以方便地进行模型的创建、参数设置和计算任务的提交。在计算完成后，后处理模块能够以直观的图形、图表等方式展示计算结果，如应力云图、位移云图、时程曲线等，便于用户对模拟结果进行分析和评估。4.2有限元模型的建立在建立砖塔有限元模型时，遵循简化原则对模型进行适当简化，以提高计算效率并突出主要结构特征。去除了一些对整体动力性能影响较小的细节部分，如塔身上微小的装饰雕刻等，但保留了对结构力学性能起关键作用的部分，如塔身、塔檐、斗拱和基础等。通过这种简化方式，既能保证模型能够准确反映砖塔的主要力学行为，又能避免因模型过于复杂而导致计算量过大，影响计算效率和结果的准确性。在模型构建过程中，采用实体单元对砖塔结构进行模拟。实体单元能够较好地模拟砖塔的三维结构特征和复杂的受力状态，精确地反映结构在地震作用下的应力应变分布情况。对于塔身、塔檐等主要结构部件，选用八节点六面体线性减缩积分单元（C3D8R）。该单元在保证计算精度的同时，能够有效减少计算量，提高计算效率。在划分网格时，根据砖塔结构的特点和分析需求，对不同部位采用不同的网格密度。对于结构的关键部位，如塔基与塔身的连接处、塔檐与塔身的连接处以及门窗洞口周围等容易出现应力集中的区域，采用较小的单元尺寸，加密网格，以提高计算精度，准确捕捉这些部位的应力应变变化。而对于结构相对简单、受力较为均匀的部位，则适当增大单元尺寸，降低网格密度，在不影响计算精度的前提下，减少计算量。通过合理的网格划分，既能保证模型的计算精度，又能提高计算效率，使模拟结果更加准确可靠。材料本构关系模型的选择对数值模拟结果的准确性至关重要。考虑到砖塔在地震作用下材料的非线性行为，选用混凝土塑性损伤（CDP）模型来描述砖和砂浆的力学性能。该模型能够较好地模拟材料在受拉和受压状态下的非线性力学行为，包括材料的开裂、损伤演化和塑性变形等现象。在CDP模型中，通过定义混凝土的单轴受压应力-应变关系、受拉应力-应变关系以及损伤参数等，来准确描述砖和砂浆的材料特性。对于砖的材料参数，根据试验测定的砖的抗压强度、抗拉强度、弹性模量等数据，确定其在CDP模型中的相应参数。例如，已知砖的抗压强度为f_{cu}，则在模型中设置砖的单轴受压峰值应力为f_{cu}，并根据相关规范和试验数据确定受压应力-应变曲线的形状参数。对于砂浆的材料参数，同样依据试验结果进行设置。同时，考虑到砖与砂浆之间的粘结作用，通过设置合适的界面参数来模拟两者之间的相互作用。例如，定义界面的粘结强度、摩擦系数等参数，以准确反映砖与砂浆在地震作用下的粘结和滑移行为。通过合理选择和设置材料本构关系模型及相关参数，能够更真实地模拟砖塔在地震作用下的力学响应和损伤演化过程。4.3模型模态分析在进行有限元模型的模态分析时，首先需确定合理的边界条件。根据砖塔实际的工作状态，将砖塔底部与基础的连接视为固定约束，即限制砖塔底部所有节点在X、Y、Z三个方向的平动位移和转动位移。这样的边界条件设置能够模拟砖塔在实际中基础对塔身的约束作用，确保模型在分析过程中的力学行为与实际情况相符。通过在ABAQUS软件中执行模态分析命令，设置合适的分析参数，如提取的模态阶数等，得到砖塔模型的自振频率和振型。经计算，得到了模型前X阶自振频率，结果如表4-1所示。模态阶数自振频率（Hz）1f12f23f34f45f5表4-1砖塔模型前X阶自振频率从表中数据可以看出，随着模态阶数的增加，自振频率逐渐增大。第一阶自振频率反映了砖塔整体的基本振动特性，其频率相对较低，表明砖塔在一阶振动时的振动周期较长，振动较为缓慢，主要表现为整体的侧移振动。高阶模态的自振频率较高，对应着砖塔结构的局部振动或复杂的振动形态。通过ABAQUS软件的后处理功能，还得到了砖塔模型各阶自振频率对应的振型。一阶振型呈现出整体弯曲的形态，砖塔的顶部位移最大，随着高度的降低，位移逐渐减小，这与砖塔在水平地震作用下的整体侧移变形特征相符。二阶振型则表现为在一阶振型的基础上，出现了一个反弯点，结构的振动形态更为复杂。高阶振型中，还出现了局部构件的振动，如塔檐的局部振动等。这些振型特征反映了砖塔结构在不同振动模式下的变形特点，为深入理解砖塔的动力特性提供了直观的依据。将有限元模型计算得到的自振频率和振型与试验结果进行对比，以验证有限元模型的准确性。在试验中，通过环境激励法或自由振动法等测试手段，得到了砖塔模型在初始状态下的自振频率和振型。对比结果表明，有限元模型计算得到的自振频率与试验值在一定误差范围内基本吻合。对于一阶自振频率，有限元计算值为f_{1FEM}Hz，试验值为f_{1Exp}Hz，相对误差为\Deltaf_1=\frac{|f_{1FEM}-f_{1Exp}|}{f_{1Exp}}\times100\%，误差在合理范围内。振型方面，有限元模型得到的各阶振型与试验观察到的振动形态基本一致，能够准确反映砖塔模型在不同振动模式下的变形特征。通过对比分析，验证了所建立的有限元模型能够较好地模拟砖塔的动力特性，为后续的地震响应分析提供了可靠的基础。4.4ABAQUS动力时程分析方法动力学基本方程是描述结构在动力荷载作用下运动规律的基础，对于砖塔这类结构，其动力学基本方程可表示为：M\ddot{u}+C\dot{u}+Ku=F(t)其中，M为结构的质量矩阵，它反映了结构各部分的质量分布情况，对于砖塔模型，质量矩阵的计算需要考虑砖塔的几何形状、材料密度等因素；C为阻尼矩阵，阻尼在结构振动中起到消耗能量的作用，常见的阻尼模型有瑞利阻尼等，在砖塔模型中，阻尼的取值会影响结构在地震作用下的响应幅值和衰减特性；K为刚度矩阵，它体现了结构抵抗变形的能力，砖塔结构的刚度与砖块和砂浆的材料性能、结构的几何形状以及各部分之间的连接方式密切相关；\ddot{u}、\dot{u}、u分别为加速度、速度和位移响应向量，它们描述了结构在动力荷载作用下的运动状态；F(t)为动力荷载向量，在砖塔的地震响应分析中，F(t)通常是根据输入的地震波来确定的，它反映了地震作用对砖塔结构的激励。在ABAQUS中，求解动力学基本方程的时间积分方法主要有隐式算法和显式算法，二者各有特点。隐式算法在求解过程中，考虑了当前时刻以及前一时刻的状态变量来计算下一时刻的响应。以Newmark法为例，它是一种常用的隐式时间积分方法。在Newmark法中，假设加速度在时间步长内按线性变化，通过对动力学基本方程进行离散化处理，得到关于位移、速度和加速度的线性方程组。在每个时间步内，需要求解一个大型的线性方程组，以得到结构的响应。隐式算法的优点是计算稳定性好，对于大多数结构动力学问题都能得到较为准确的结果。它可以处理较大的时间步长，对于一些变化较为缓慢的动力荷载，能够有效地减少计算量。然而，隐式算法在求解线性方程组时，需要较大的计算资源和内存空间，尤其是对于大型复杂结构，计算成本较高。在砖塔的地震响应分析中，如果砖塔结构较为复杂，节点和单元数量较多，使用隐式算法可能会导致计算时间过长，对计算机硬件要求较高。显式算法则是基于当前时刻的状态变量来直接计算下一时刻的响应。以中心差分法为例，它是一种显式时间积分方法。在中心差分法中，加速度通过位移和速度的差分来近似计算，不需要求解联立方程组。显式算法的时间步长通常由结构的最小特征尺寸和材料的波速决定，为了保证计算的稳定性，时间步长必须足够小。显式算法的优点是计算效率高，每个时间步的计算成本较低，适用于求解高速动态问题和涉及大量接触非线性的问题。在砖塔的地震响应分析中，如果关注砖塔在地震作用下的瞬间响应和破坏过程，显式算法能够快速得到结果。然而，显式算法由于时间步长小，需要进行大量的时间步计算，对于长时间的动力分析，计算量会显著增加。而且，显式算法在处理一些复杂的接触和摩擦问题时，可能会出现数值振荡等不稳定现象。在ABAQUS中进行砖塔的动力时程分析时，需要根据具体问题的特点和要求，合理选择时间积分方法。如果砖塔结构较为简单，对计算精度要求较高，且计算资源充足，隐式算法可能是更好的选择。而如果关注砖塔在地震作用下的快速响应和破坏过程，或者砖塔结构存在大量的接触非线性，显式算法则更具优势。在实际分析中，也可以通过对比两种算法的计算结果，来验证分析的准确性和可靠性。4.5数值模拟结果在ABAQUS中完成动力时程分析计算后，对模拟结果进行深入分析，以全面了解砖塔在地震作用下的力学行为和损伤演化过程。首先，计算模型的阻尼系数。阻尼系数的准确计算对于理解结构在地震作用下的能量耗散机制至关重要。采用瑞利阻尼模型来计算阻尼系数，该模型假设阻尼矩阵是质量矩阵和刚度矩阵的线性组合，即C=αM+βK，其中α和β为瑞利阻尼系数。通过对砖塔模型的前两阶自振频率ω_1和ω_2以及对应的阻尼比\xi_1和\xi_2进行计算，利用公式\xi_i=\frac{\alpha}{2\omeg