中图分类号:TV222 文献标志码:A
黄河水资源的合理调配对流域内工农业生产和居民生活用水至关重要。黄河下游地区通过引黄涵闸对黄河水进行有效引用,引黄涵闸的合理设计直接影响水资源利用效率与工程安全。传统的引黄涵闸闸室设计方法耗时费力,当面对复杂的设计参数与工况时,易出现计算误差。随着计算机技术飞速发展,利用编程语言辅助工程设计成为提升设计效率与精度的重要途径。AutoLISP语言是一种嵌入AutoCAD软件中的编程语言,具有强大的图形处理与数据计算能力[1],将其应用于黄河下游引黄涵闸闸室设计,有望显著优化设计流程,提高设计质量,为水利工程建设注入新的活力。
1项目背景
濮阳段黄河下游引黄涵闸设计项目面临复杂的地理环境与多样的水利需求。该区域地势平坦,黄河水位受季节和上游来水影响波动较大,同时需要兼顾周边农业灌溉、工业用水和生态补水等任务。工程区出露地层为第四系全新统人工填土( Q4f )和第四系全新统河流冲洪积层(
),岩性主要为堤身土、粉细砂、砂壤土、粉质壤土和粉质黏土。闸室和涵洞段地基基础位于第2层的壤土层上,该层存在渗透变形、地震液化、抗滑稳定、沉降变形、抗冲刷掏刷、基坑开挖边坡稳定性以及基坑排降水等工程地质问题。基础持力层为第2层壤土层,第4层粉砂层可作为较好的桩端持力层。传统设计手段难以快速、准确地处理大量设计数据与复杂工况。AutoLISP语言开发用于该项目具有显著优势,其与AutoCAD软件紧密集成,能够充分利用AutoCAD强大的绘图功能,便捷、高效地绘制引黄涵闸相关图纸。通过编写AutoLISP程序,可以实现设计参数自动输入、图形快速生成以及复杂计算自动化,极大缩短设计周期[2。例如,当绘制三维模型时,利用AutoLISP语言能够迅速根据输入参数构建精确的模型,而无须手动逐一绘制各个部件。濮阳段黄河下游引黄涵闸设计参数详情见表1。
2黄河下游引黄涵闸闸室设计
2.1设计流程
黄河下游引黄涵闸闸室设计流程如下:首先,收集并整理详细的设计参数,包括地质、水文等基础数据。其次,根据这些参数,在AutoCAD环境下编写AutoLISP程序,绘制精确的三维模型。再次,程序自动识别不同的计算工况及其对应的荷载组合。最后,分别对涵闸渗流稳定与闸室稳定进行计算、分析,根据计算结果,对设计进行优化调整(黄河下游引黄涵闸闸室设计流程如图1所示)。
2.2输入设计参数
在数据录入环节,AutoLISP程序运用严谨的数据校验算法。以水位数据为例,程序预先设定合理水位区间,该区间根据黄河下游长期水文监测数据的统计学分析得出,包括不同季节、年份的水位波动范围。当设计人员输入水位数据时,程序会实时比较输入值与预设区间,一旦发现输入值超出范围,立即触发基于异常处理机制的提示弹窗,保证数据在初始阶段的准确性。对于地质数据,例如地基承载力特征值与土壤内摩擦角等,AutoLISP程序构建了专门的数据结构。采用哈希表结合链表的数据存储方式,将地质参数的类别作为哈希表的键值,链表用于存储同类别下的详细参数信息[3]。当后续构建计算模型调用数据时,这种数据结构设计能够大幅提升检索效率,降低时间复杂度。
在计算闸室基础稳定性过程中,程序能够快速定位并获取所需的地基承载力特征值,为复杂的力学计算提供精准数据支持。此外,程序还对输入参数进行数据类型匹配校验。地震加速度应为数值型数据,如果输入字符型数据,那么程序会根据数据类型识别算法及时提示设计人员修正,避免数据类型错误导致后续计算错误。这些输入设计参数不仅决定了闸室主体结构的选型与尺寸设计,而且在涵闸渗流稳定与闸室稳定计算中,其会作为核心变量参与复杂的数学模型运算,任何微小偏差都可能在历经多轮计算迭代后导致设计结果显著偏离,因此输入参数的精确性与可靠性至关重要。
表1濮阳段黄河下游引黄涵闸设计参数

图1黄河下游引黄涵闸闸室设计流程

2.3绘制三维模型
获取精准无误的设计参数后,AutoLISP语言基于强大、卓越的图形绘制能力,在AutoCAD平台上迅速、精准地构建引黄涵闸闸室的三维模型[4。程序根据输入的闸底板参数,精准确定底板的长、宽、高尺寸以及空间位置。在此基础上,按照既定顺序依次绘制闸墩、闸门等核心部件。为了增强模型的直观可视化效果,程序会赋予不同部件高度仿真的材质与鲜明颜色,将闸墩模拟为混凝土材质外观,呈现出粗糙、坚固的质感;将闸门设置为金属材质外观,展现光泽与冷硬感。同时,借助AutoLISP丰富的三维坐标变换函数,设计人员能够便捷地对模型进行全方位旋转、灵活缩放等操作,可从任意角度细致观察模型,全面检查设计的合理性(引黄涵闸闸室三维模型如图2所示)。
2.4识别计算工况及荷载
在实际运行历程中,黄河下游引黄涵闸需要应对正常运行、检修维护和洪水侵袭等多种复杂工况。AutoLISP程序能够根据输入的特征水位和其他相关关键参数,运用智能算法自动、精准识别不同计算工况。在每种工况下,程序进一步基于预先构建的荷载计算模型,综合确定作用在闸室上的荷载组合。这些荷载包括水压力、结构自重、土压力和地震力等。在正常运行工况下,水压力与闸室自重占据主导地位;在洪水工况下,水压力会急剧攀升,同时需要充分考量地震力可能带来的叠加影响。程序利用严谨的数学运算,准确计算各荷载在不同工况下的数值和作用方向,并将其无缝对接为后续渗流稳定与闸室稳定计算的核心输人条件。
2.5计算涵闸渗流稳定
引黄涵闸的渗流稳定是保障工程长期安全运行的核心要素。利用AutoLISP编写的渗流计算程序,以达西定律为理论根基,结合先进的有限元分析方法,对涵闸地基及其周边土体的渗流状况进行精细化模拟计算。程序将闸室及其周边土体精细划分为众多计算单元,进行反复迭代运算,求解各单元的水头分布情况。并在计算进程中充分考虑不同土层渗透系数的显著差异,运用专业算法对计算结果进行科学修正[5。该段引黄涵闸渗流稳定性结果见表2。根据计算输出的各分段阻力系数、修正前/后的水头损失以及各部位水头数据,如果某部位水头异常偏高或渗流坡降逾越允许阈值,那么表明存在渗流不稳定风险,此时需要及时调整设计方案,例如增设防渗墙、优化地基加固处理方式等。
图2引黄涵闸闸室三维模型

表2该段引黄涵闸渗流稳定性结果(单位:m)

2.6计算闸室稳定
闸室稳定计算对保障黄河下游引黄涵闸安全运行具有至关重要的作用,主要包括抗滑稳定与抗倾稳定2个核心。在抗滑稳定计算中,关键计算过程如公式(1)所示。

式中: Ks 为抗滑稳定安全系数,是衡量闸室抗滑动能力的关键指标; f 为基底摩擦系数,该系数与地基土的性质、闸底板和地基的基础情况紧密相关,一般需要根据地址勘察数据与工程经验确定; ΣG 为作用在闸室上的竖向荷载总和,包括闸室自身结构的重力、上部附属设施重力和可能存在的其他竖向荷载; ΣH 为作用在闸室的水平荷载总和。
抗倾稳定计算过程如公式(2)所示。

式中: Kt 为抗倾稳定安全系数,反映了闸室抵抗倾覆的性能; ΣM#L 为抗倾力矩总和,其由竖向荷载对倾覆轴产生的力矩和部分水平荷载形成的抗倾力矩构成;
为倾翻力矩总和,源于水平荷载对倾覆轴的力矩作用。
利用AutoLISP程序,根据识别出的计算工况和荷载,并紧密结合上述公式进行闸室稳定计算。该段闸室稳定计算结果见表3。
表3闸室稳定计算结果

表4不同工况下AutoLISP程序与Excel的计算结果比较

3结果与分析
针对濮阳段引黄涵闸,运用AutoLISP编写程序,对6种不同工况条件下的闸室稳定情况展开计算分析。这6种工况包括正常运行、检修、小洪水、大洪水、地震单独作用以及洪水与地震组合作用工况,全面模拟引黄涵闸在实际运行中可能面临的复杂场景。
在正常运行工况下,AutoLISP程序根据输入的正常水位、结构自重等参数,代入抗滑稳定公式(1)与抗倾稳定公式(2)进行运算。程序快速检索与调用数据结构中存储的参数,在极短时间内完成计算,输出抗滑稳定安全系数 Ks1 与抗倾稳定安全系数 Kt1 。在检修工况下,考虑部分设备拆除、无水压力等情况,程序同样能够迅速调整计算参数,精准得出相应工况下的稳定系数 Ks2 与 Kt2? 在洪水工况下,无论是小洪水还是大洪水,程序均能根据水位骤升带来的水压力变化,自动更新荷载计算模型。在计算水压力过程中,程序运用流体力学原理并结合实际水位高度与闸室尺寸,准确计算出不同位置水压力的大小与方向,进而精确计算抗滑与抗倾稳定系数 Ks3 / Kt3 与 Ks4 Kt4 。在地震单独作用工况下,程序根据输入的地震参数,例如地震加速度等,运用动力学原理计算地震力对闸室的作用,得出稳定系数 Kss 与 Kt5 。在洪水与地震组合作用工况下,程序综合考虑水压力与地震力的叠加影响,进行更复杂的耦合计算,输出稳定系数 Ks6 与 Kt6 。不同工况下AutoLISP程序与Excel的计算结果比较见表4。
将上述AutoLISP程序计算结果整理成详细的程序输出报告,报告中不仅记录了每种工况下的最终稳定计算结果,而且包括各工况下的荷载取值、参数设定等中间计算过程。同时,使用传统的Excel建表方式对相同数据进行处理和计算。在Excel中,设计人员需要手动输入大量公式来构建计算模型,在复杂工况下,公式嵌套层数多、容易出现逻辑错误,例如在洪水与地震组合工况下,Excel需要多次手动调整公式,以考虑不同荷载的组合。而AutoLISP程序会根据预设的计算逻辑,自动完成参数更新与计算。在计算效率方面,AutoLISP程序具有高效的算法与数据处理能力,计算6种工况的闸室稳定情况仅需数秒,而Excel完成同样计算可能需要数分钟甚至更长时间,尤其在数据量增大、工况复杂程度较高工况下,差距更显著。在计算精度上,AutoLISP程序采用专业的工程计算模型,严格遵循力学原理与规范要求,其计算结果、理论分析与实际工程经验高度契合(见表4)。在洪水与地震组合作用的复杂工况下,Excel计算的抗滑稳定安全系数为1.22,而AutoLISP程序计算结果为1.25,与专业力学分析比较后可知,AutoLISP程序结果更接近真实值,充分表明AutoLISP程序在引黄涵闸闸室设计中计算闸室稳定的优势与可行性,能够为工程设计提供高效、精准的技术支撑,大幅提升设计质量与效率。
4结语
本文将AutoLISP语言应用于黄河下游引黄涵闸闸室设计,尤其是濮阳段项目,从项目背景分析到完整的设计流程实现,再到结果比较和分析,AutoLISP不仅能够快速处理复杂的设计参数与多样的计算工况,而且能够精确计算涵闸渗流稳定与闸室稳定,结果表明,AutoLISP语言功能强大,在提升设计效率与精度方面成效显著。
参考文献
[1]代敬龙.基于AutoLISP的选煤厂非标参数化设计研究煤矿机械,2024,45(11):214-216.
[2]张碧,于文涛,赵玉波.Autolisp在总图设计中的探索与应用].工程建设与设计,2024(9):177-180.
[3]陶文东.基于AutoLISP的AutoCAD机械典型零件二次开发[].现代制造技术与装备,2023,59(7):42-44.
[4]郑伟安,孟大鹏,王秀元,等.一种基于AutoLISP的地形图斜坡类要素处理程序的设计与实现.山东国土资源,2023,39(4):58-63.
[5]许冬亮.AutoCAD二次开发在大比例尺地形图更新中的应用研究.测绘与空间地理信息,2023,46(3):162-164.
