首页/文章/ 详情

“力学大师”系列之伽辽金:力之大者,为国为民

28天前浏览342
关键词:力学,伽辽金法,有限元,加权残差法,苏俄

全文约7000字

“朝为田舍郎,暮登天子堂”,读书作为普通人逆天改命的重要机会,一直以来被赋予了很多功利性。然而除却这个目的,学习本身,何尝不是身心淬炼,升华蜕变的必经之路呢?作为理论与工程兼备的学科,力学的学习更是如此。学到深处,便胸怀大志,腹有良策,有包藏宇宙之机,吞吐天地之志。伽辽金靠着一股劲硬从最底层杀出来,饥贫、蹲监、战争都被他硬扛过去,将相本无种,男儿当自强!



只要是搞有限元算法,无论是CFD还是CAE,抑或是电磁,伽辽金都是绕不开的人物。伽辽金提出的伽辽金法与变分法、有限元法深度绑定,是将连续介质力学的偏微分方程转化为计算机可识别的代数方程组的核心桥梁。

不同于之前介绍的几位力学大师年纪轻轻便名满天下,伽辽金40岁前的生活极为坎坷。我现在刚刚30出头,在查阅伽辽金人生履历的时候非常感慨,他40岁前几个关键的节点如果换成我,大概率是扛不过去的,肯定泯然众人了。

他的一生,剧本密集到离谱:寒门抄写员、半工半读的穷学生、职业革命家、沙皇监狱里的政治犯、大学教授、苏联科学院院士、国家超级工程总顾问,最后还披上了一身苏维埃中将的军服。

用现在的话说,这是一个把“考研、进厂、上岸、踩缝纫机、评职称、拿帽子”全流程都刷了一遍的人。

铁汉子百炼成钢,哪怕到了晚年,力学人的钢筋铁骨仍然在伽辽金身上。开发了这么久的有限元软件,终于有机会好好学习一下这位大师的生平。

限于本人的水平,内容难免有错漏。

1. 寒门的出路

1871年3月4日(俄历2月20日),伽辽金出生在俄罗斯帝国维捷布斯克省的波洛茨克(今属白俄罗斯)。

在那个阶级固化、排犹情绪严重的沙皇俄国时代,贫困与边缘化的身份往往意味着一个天才尚未启蒙便要夭折。伽辽金的家境极其拮据,虽然在当地拥有一处简陋的旧房产,但微薄的收入根本无法维持一家人的生计。

父亲观念极端保守,反对儿子去上俄语学校。这一条几乎直接封死了正常升学通道——在那个年代的俄国,不接受俄语体系的正规中学教育,基本就告别了高等院校。

1883年,为了糊口和补贴家用,年仅12岁的伽辽金辍学打工,在波洛茨克当地的孤儿法庭当抄写员。

当时的法庭抄写是一份极度枯燥的工作:在昏暗的煤油灯下,日复一日地誊抄法律文书、判决书和账册,字迹不能涂改,不能潦草。这段童工经历塑造了他一辈子标志性的书写习惯:字写得小巧、紧凑、没有任何多余装饰,但异常清晰易读。这一经历也使他间接养成了严谨细致的作风。

今天讨论“寒门难出贵子”的时候,大家争论的是学区房、辅导班和信息差。伽辽金的起点是:没有学区,没有辅导,没有信息,只有一盏煤油灯和抄不完的案卷。他白天抄文书糊口,晚上刻苦自学数学、物理和外语。

很多时候即便是天才,也需要学习的动机(契机)。爱因斯坦从小就在父亲和叔叔的熏陶下,对数学非常感兴趣。牛顿则是在剑桥大学研究生舅舅的鼓励下学习。我没有查到资料能说明,像伽辽金这样出身最底层的人,究竟是什么让他对学习充满热情。或许,就是天生的热爱学习,又或许小小年纪的他当时就明白了:读书改变命运。

按照当时的学制,普通学校毕业生想报考高等工科院校,还得补考一个追加学年的课程。1893年,22岁的伽辽金只身跑到明斯克,以校外生的身份参加了公立中学的毕业考试——可以理解为沙俄版的社会考生自考。他一次通过全部科目,同年秋天考入圣彼得堡国立技术学院。

人生终于看到了一些希望。

2. 毕业生与地下党

进了大学不等于脱贫。没有家庭供养,伽辽金的大学生活是一场体力与意志的拉力赛:白天上课,晚上兼 职,当过绘图员,也当过私人家庭教师,教数学。放到现在,这就是一个同时打着三份工还要保证学业的大学生。

人不会总是走背字,在学业上,他遇到了对他一生学术道路产生深远影响的导师——俄罗斯工程力学奠基人之一、技术学院院长、圣彼得堡力学学派的掌门、维克多·勒沃维奇·基尔皮乔夫(Viktor Lvovich Kirpichov)。基尔皮乔夫的主张非常明确:工程结构分析必须和数学物理理论深度结合,反对工程师只会套经验公式。

在基尔皮乔夫门下,伽辽金既练了锻造、铸造、机械加工这些实打实的工程手艺,又在数学物理方程和弹性力学上打下了硬底子。

当时民主主义思想广泛传播,列宁已经将马克思主义带入了沙俄。圣彼得堡并不仅仅是学术的圣地,它同样是革命风暴的策源地。当时的圣彼得堡学生群体思想极其活跃,面对沙皇政权的腐 败与工人的悲惨境遇,许多青年学生纷纷投身政治。伽辽金作为一名来自社会底层的苦学生,对工人群体的苦难有着天然的同理心。

1899年,临近毕业的伽辽金因为参加学生罢 课被遣送回原籍,随身只发了一张“通行证”(到地方后须向警察报到)。好在当局很快允许被遣返的学生回来,这年年底,他以校外生身份完成学业,拿到工程师资格。

同样是1899年,他加入了俄国社会民主工党(RSDLP)——后来苏联共产党的前身。别人毕业拿双证,他毕业拿双身份:工程师,外加一个地下党员。这一选择也改变了他此后的命运轨迹。

3. 工程师、革命家、监狱与学术

毕业后刚参加工作的前五年非常重要,这个时间基本决定了人未来的事业走向。有的人开始柴米油盐,不再学习新技术钻研新方法,踏实过日子。有些人开始理论联系工程实际,对力学钻得更深。当然这五年也是大家结婚生子、买车买房、升职加薪的关键期。

不同的是,革命年代,更富传奇。

毕业后的头三年,伽辽金在哈尔科夫的俄国蒸汽机车与机械制造联合公司工厂当工程师,同时给工人夜校讲课。

1903年底,他远赴远东,参加中东铁路(东清铁路,横贯中国东北)的建设。只干了半年,就因能力出色被调回圣彼得堡,出任北方机械与锅炉厂的技术部主任,负责大型锅炉和承重金属结构的计算。

如果剧本照此发展,他大概率会成为帝俄工业界一位优秀、务实、逐渐发福的总工程师。但1905年来了。

这一年,俄国爆发第一次革命。列宁在伦敦主持召开了俄国社会民主工党第三次代表大会,制定了党在民主革命中的策略路线。11月,列宁由国外回到圣彼得堡,直接领导布尔什维克的工作。

伽辽金参与组织圣彼得堡的工程师联盟,因组织工程师罢 工被捕,关了一个多月。1906年8月5日(俄历7月23日),警察包围了乌杰利纳亚车站附近阿列克谢耶夫街13号的一栋房子,工党圣彼得堡委员会的成员几乎被一锅端。被捕现场有人向警察开了枪——按常理这是掉脑袋的重罪,但1907年3月26日圣彼得堡法 院的判决却轻得出奇:19名被告中,1人判两年,伽辽金等8人各判一年半,其余释放。

伽辽金被关进了圣彼得堡最臭名昭著的“十字架”监狱,人生上升的势头被时代打断。

时也命也,谁都没想到,这堵高墙成了他学术生涯真正的孵化器。其实不少大佬都有类似的经历,文王拘而演周易,牛顿在疫情爆发时窝在家乡搞出万有引力定律,霍金坐轮椅后开始研究大爆炸理论,胡子哥在监狱里写《我的奋斗》...

在狱中简陋的监房里,伽辽金申请获得了纸笔与力学书籍。远离了工厂的噪音与革命的喧嚣,伽辽金将全部精力倾注于弹性力学中最深奥的问题之一——受压杆件的纵向弯曲与稳定性。

欧拉早在18世纪就给出了理想压杆临界载荷的公式,但在实际工程中,结构往往由多跨、变截面、带有复杂约束的纵横交错杆系组成,极难获得解析解。

伽辽金在狱中克服了没有计算机、没有计算器的困难,完全依靠极其精密的代数手算,开创性地将欧拉纵向弯曲理论推广至多跨复杂组合柱系统。

在这本长达130页、完全在狱中手写完成的巨幅手稿中,伽辽金不仅发展了纵向弯曲理论,更展现出了他未来将数学工具与工程实际完美接轨的独特天赋。

1908年刑满出狱。1909年,这篇诞生于铁窗之下的论文,在《圣彼得堡实用技术学院学报》上正式发表。这篇论文逻辑严密、计算精确,能够直接应用于桥梁骨架、大型厂房铁架以及船舶结构的稳定性设计,立刻在俄罗斯力学界引发了轰动。

同年,受圣彼得堡国立技术学院邀请,38岁的伽辽金正式进入该校担任结构力学讲师,开启了他长达三十余年的大学教授生涯。

伽辽金的人生路,终于开始好走起来。

4. 伽辽金法:把论文写在祖国的大地上

4.1 工程中思考的近似思想

伽辽金深刻认识到,力学绝不能闭门造车,必须深刻理解当时世界上最先进的材料、工艺与结构形式。

从1909年到1914年第一次世界大战爆发前夕,伽辽金利用假期与学术休假,频繁自费前往欧洲工业强国进行学术考察与技术调研。

他先后访问了德国、奥地利、瑞士、比利时、瑞典等国的顶尖重工业企业、造船厂、大型桥梁工地以及水力发电站。他深入研究了欧洲当时刚刚兴起的无梁楼盖结构、薄板结构以及大型电站钢结构框架。

不同于普通的走马观花,伽辽金每到一个城市,最喜欢去的地方不是博物馆,而是正在热火朝天建设中的建筑工地。

他详细记录欧洲工程师在处理复杂受力构件时的实际做法。这段游历让他敏锐地察觉到了横亘在“理论”与“工程实践”之间的巨大鸿沟:

“经典弹性力学理论已经发展到了极高的水平,建立了一系列精美的偏微分方程;然而,这些方程只有在极少数边界条件极其简单的理想情况下才能求出精确的解析解。而在实际工程中,工程师们面对的往往是奇形怪状的非对称板、带有无数加强筋的桥梁骨架。面对这些偏微分方程,欧洲的工程师们往往只能依靠经验公式或过大的安全系数。”

“必须找到一种简便、通用、精度可控的偏微分方程近似积分方法!” 这一想法在伽辽金的脑海中深深扎根。

1913年,圣彼得堡电厂计划建造一座具有前沿水平的大跨度重型金属框架结构,要求该建筑必须能够承受巨大的机械动力载荷和极高的温度变化。

伽辽金作为首席咨询顾问,亲自设计并计算了该电厂的金属骨架。他将自己多年来研究的金属框架屈曲理论、多层构件内力分配方法悉数应用其中。该结构成为俄罗斯首批承受超重载荷的特种金属框架之一,被欧洲工程界誉为结构工程的经典范例。

4.2 伽辽金方法

变分法中有一个著名的里茨法。

1908年,德国数学家沃尔特·里茨(Walther Ritz)提出了一种近似求解变分问题的方法:把未知解表示为一组基函数的线性组合,然后通过让系统总能量最小来定出组合系数。里茨法漂亮、强大,但有一个前提——你的问题必须能写出能量泛函。一旦碰上写不出能量表达式的物理问题(非保守系统、耗散问题、各种复杂偏微分方程),里茨法直接熄火。

1913年,俄国造船学家布勃诺夫在评论克雷洛夫的著作时,提出过一个类似的做法:直接把级数展开式代入微分方程,乘以坐标函数再积分。但他本人只把这个技巧当成里茨法的一个变体,没有多想。

1915年,伽辽金在《工程师通报》上发表了一篇论文,题目朴实无华:《杆件与平板:关于杆件与平板弹性平衡某些问题的级数解法》。

这篇论文,就是后来“伽辽金方法”的出生证。

它的逻辑干净到近乎暴力,拢共三步。

第一步,猜。假设微分方程的解,可以用一组事先选好的、满足边界条件的基函数的线性组合来逼近——系数待定。

第二步,算误差。把这个近似解直接代回原微分方程。因为是猜的,方程两边当然不相等,多出来的那一截,叫“残差”(Residual)。

第三步,消灭误差。伽辽金要求:残差和每一个基函数都正交——残差乘以每个基函数、在整个求解区域上积分,结果必须等于零。有几个待定系数,就有几个正交条件,正好凑成一个代数方程组。解出系数,完事。

里茨法是从“能量最低”出发的,只有守规矩的保守系统才有能量可谈;伽辽金方法直接从微分方程本身出发,根本不问你有没有能量泛函。

这一处理,把方法的适用范围从弹性静力学踹向了几乎整个数学物理:流体力学、传热、电磁场、声学,统统通用。普通的学者解决一个问题,大师一出手就是一类问题。

他还用虚位移(或然位移)原理给方法配上了力学解释,并在论文里一口气算了一大批杆件与平板的弯曲、稳定算例。

4.3 有限元

伽辽金法简直就是给有限元法量身定制的。在我们建立了控制方程以后,对于每一个离散的单元都可以按如下流程处理:

(1) 根据单元类型,确定插值函数。此时速度、压力等变量,都可以用权函数表达。

(2) 采用伽辽金方法,权函数=插值函数,控制方程与权函数相乘,积分取0。

(3) 在每个单元域内,方程转换为权函数的积分形式,最终形成单元矩阵。

我们在开发的流动求解器、传热求解器中,都使用了这套技术。

静界有限元与您共研工业软件技术

4.4 大国工程师

1917年,十月革命一声炮响,改变了整个国家的命运,也彻底改变了伽辽金的命运。

由于早年的革命经历和顶尖的学术水平,已是中年的伽辽金迅速获得重用:

1920年:伽辽金被任命为圣彼得堡理工学院结构力学系主任。

1923年—1929年:出任圣彼得堡理工学院土木工程学院院长。

1924年—1929年:兼任列宁格勒国立大学结构力学教授、列宁格勒铁道运输工程学院弹性力学教授。

1928年:当选为苏联科学院通讯院士。

1935年:当选为苏联科学院院士。

这个阶段他的学术产出同样是巅峰。1930年,他在三维弹性力学上取得突破,证明各向同性弹性体的位移可以用一个向量函数来表示——这就是后来教科书上标准的“伽辽金向量”,成为求解三维弹性力学问题的通用工具之一。1933年,他出版了371页的专著《弹性薄板》,汇总多年研究薄板弯曲的解析与数值成果。在没有计算机的年代,这本书就是结构工程师案头的“计算手册”。

在工程界,他是苏联工业化背后的“定海神针”:沃尔霍夫水电站、康多波加纸浆厂、国家钢铁设计院、有色金属设计院的技术委员会里都有他;他还编制了梯形截面水坝与挡土墙的应力计算表。第一个五年计划的明珠——第聂伯河水电站,当时欧洲最大的水电站,他参与了设计咨询;1932年电站竣工,他出任国家验收委员会委员。

1936年4月,苏联政府一纸命令,任命他为苏维埃宫钢结构骨架与墙体楼板初步设计鉴定委员会主席。

苏维埃宫规划高度415米,顶部还要立一尊巨型列宁雕像,含尖塔接近500米,建成即是当时全球最高、也最浮夸的建筑。

这个审查任务政治意味和技术难度双高,伽辽金带着团队对钢框架承载方案做了极为严苛的论证。

后来地基都打好了,二战来了,项目无限期搁置。

5. 力学将军

就在红色帝国全力发展的时候,希特勒闪击波兰,开始了征服欧洲的进程。为了避免两线作战,希特勒与斯大林签署《苏德互不侵犯条约》以稳住苏联。

可是体量与资源庞大的苏联是希特勒眼中的一块大肥肉。1940年,德军绕开固若金汤的马奇诺防线,一举拿下法国,“西线无战事”,终于可以回头攻苏。于是,1940年12月18日,希特勒正式签署入侵苏联的首个作战命令——第21号命令,代号“巴巴罗萨”闪击列宁格勒(圣彼得堡)。

不久,德军兵临列宁格勒城下,开始了惨绝人寰的“列宁格勒大封锁”。

此时的伽辽金院士已是古稀之年的老人,且兼任列宁格勒军事工程技术大学结构力学系主任。政府原本安排他第一批乘坐飞机撤离至后方安全区域,但伽辽金斩钉截铁地拒绝了撤离安排,坚持留在战火纷飞的列宁格勒。

1941年夏天,列宁格勒市成立了列宁格勒防御工事建设指导委员会。委员会中虽有数位院士,但唯有伽辽金拥有极其丰富的土木、建筑与结构工程实战经验。伽辽金事实上成为了整个列宁格勒防御工事建设的总工程师。

力之大者,为国为民。在敌机的狂轰滥炸与炮火下,伽辽金亲自带队奔赴现场,指导防空掩体、地下指挥所、反坦克壕沟与重型混凝土碉堡的修建。他利用自己的壳体与板结构理论,计算并优化了加固混凝土掩体抵御德军重型炮弹轰炸的最小壁厚与配筋方案,挽救了无数守军与平民的生命。

为了表彰伽辽金在战争期间对国防工程的卓越贡献,苏联最高苏维埃授予伽辽金工程兵中将军衔。

在列宁格勒被德军完全合围前夕,这位身体已经严重透支的国家级院士,才在苏联最高统帅部的严令下,随科学院力学研究所紧急撤离,疏散到了后方的莫斯科。

来到莫斯科后,伽辽金甚至顾不上调养自己在围城期间因饥饿和严寒受损的身体,立刻无缝接入了苏联科学院军事工程委员会的工作。

长期的战争高压、列宁格勒封锁期间的饥寒与繁重的国防科研工作,严重摧残了老人的健康。

1945年5月9日,苏联赢得了卫国战争的胜利。然而仅仅两个月后,1945年7月12日,伽辽金因病在莫斯科逝世,享年74岁。

6. 后记

伽辽金离开这个世界时,真正意义上的电子计算机尚未问世。他一生中的所有复杂数值计算,都是依靠他那精细的手笔在纸上一步步完成的。然而,他留下的数学思想,却在计算机时代绽放出极其璀璨的光芒。

20世纪50年代,随着喷气式飞机、核反应堆与航天器的兴起,工程界对复杂三维结构求解的需求呈爆炸式增长。

以雷·克拉夫、奥利弗·齐恩基维奇为代表的有限元先驱们发现,伽辽金方法提供了构建有限元方程最优雅、最通用的数学框架。

变分法的局限在于只能用于有能量泛函的弹性力学问题。

而伽辽金法可以无缝推广到CFD、传热学、电磁场以及声学等几乎所有连续介质物理领域。

今日所有主流商业有限元与多物理场仿真软件,其底层核心求解器无一例外都包含着伽辽金法的代码,并演化成一个庞大的算法家族:里茨-伽辽金法、彼得罗夫-伽辽金法、间断伽辽金法...


来源:静界有限元
静力学复合材料通用航天船舶建筑电子声学理论爆炸材料控制工厂
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-07-28
最近编辑:28天前
静界有限元
博士 签名征集中
获赞 27粉丝 4文章 58课程 0
点赞
收藏
作者推荐

从水滴在流场中的运动求解到气象武器

关键词:飞机结冰,水滴轨迹,结冰模型,朗缪尔,龙格库塔法,人工降雨,流体力学科学与技术的发展,是无数人智慧接力的结果。所谓愚公移山,子又生孙,孙又生子;子又有子,子又有孙,一代代前仆后继,才有今日之昌明科技。本篇尝试从人物切入,介绍一下飞机结冰计算中经典的水滴轨迹求解问题。这个问题目前在国内外已不是什么难题,相关专业的研究生毕业论文大都有相关介绍,站在前人肩膀上,只要通晓原理,很快就可以编写出程序求解。但是就是这样一个问题,涉及了数学方法、力学模型、试验技术、计算机技术,多位天才接力才得以解决。我在多年前开发结冰软件的时候,曾仔细做过相关研究,因此本打算写一篇算法介绍的文章。但在准备阶段,查阅资料的时候竟发现很多有趣的东西。为此才准备转变写法,以期分享这些发现。1. 龙格自牛顿建立经典力学与微积分手段以后,大量工程问题有了统一的分析框架,即受力分析、构建平衡方程、求解。然而很多问题构建方程容易,求解可就难了。有些有解析解的,数学家们发挥聪明才智给出了答案。但是没有解析解的,该如何处理呢?如果这些问题实际生活中碰不到也就罢了,可偏偏有的问题又很重要。尤其工业革命以后,大家需要研究天体运行以解释宇宙,研究大炮的弹道以提升军事。于是离散化的数值求解方法应运而生,最初欧拉提出了欧拉法。以弹道问题为例,通过牛顿第二定律,建立考虑重力、阻力、升力、推力、科里奥利力等等的控制方程后,可以知晓弹体的速度控制方程。这样弹道轨迹就变成了一个已知初值(初始位置)的预测未来的问题:如何计算出它下一秒在哪个位置,下下秒又在哪里。欧拉给出来的思路通俗易懂:把整个时间历程分成多步求解,假设这一秒(步)内的速度是恒定不变的,那么根据速度矢量,乘以单步的时间,再和初始位置矢量相加,就能预测下一秒的位置。这种思想(将时间分割为小步长,用当前点的切线斜率近似下一时刻的值)简单是简单了,可是精度很差,如盲人开车,容易越开越偏。自欧拉后,又过了将近150年,德国的数学家、物理学家卡尔·龙格决定解决这个问题。当时的德国数学界由克莱因等人推动的“纯数学至上”风气主导,对应用性的数学研究不够重视。龙格当时已经从纯数学研究转向光谱分析,欧拉法的一阶精度导致误差累积过快,无法满足光谱数据对精度的苛刻要求。卡尔·龙格 Carl Runge,1856 – 1927龙格的思路是:既然欧拉法是基于“当前”的斜率往前推一步,越走越偏,属于一锤子买卖。那么如果我假装往前迈个半步,然后用那个半步的斜率再回来修我当前斜率,如此来来回 回,四渡赤水呢?结果还真行,龙格的方法综合多个局部信息加权平均,得到比单点外推更精确的结果。龙格这种多阶段试探+综合修正的思想竟然暗合兵法,神似四渡赤水。不直接突围,而是将战略转移拆成多次渡河机动,每次渡河都是一次试探和一次修正,获取敌军部署、地形和时机信息,综合四次机动后的全局态势,确定战略方向,最终选择最优路线,避免了硬冲导致的被围歼风险。2. 库塔龙格的方法发表不久后,引起了马丁·威廉·库塔的注意。库塔当时正在攻读博士论文,其方向就是微分方程的数值解法。在龙格的基础上,库塔将这种“多点试探”的方法进行了高度的数学泛化。他不仅将方法扩展到了方程组,更通过严密的数学证明,推导出了今天工程界奉为圭臬的经典四阶龙格-库塔法(RK4)的系数矩阵。龙格库塔法在工程中应用极广泛,作为求解常微分方程初值问题的黄金标准,其应用场景几乎覆盖了所有需要随时间演化模拟的科学与工程领域。从弹道轨迹到病毒传播,从电路瞬态到游戏物理,它无处不在。值得一提的是,库塔不仅仅是数学家,也是空气动力学名宿。航空史上有著名的“库塔-儒可夫斯基定理”(解释了机翼是如何产生升力的),库塔就是其中的关键贡献者。库塔所处的时代正是飞机刚刚诞生(1903)的时候,数学、力学与航空工程互相融合促进。库塔一定想不到,他完成的龙格库塔法 会在半个多世纪以后用在飞机结冰上。马丁·威廉·库塔 Martin Wilhelm Kutta,1867 – 19443. 朗缪尔3.1 粒子调控就在库塔完善了龙格库塔法的同时,远在美国的20岁年轻人欧文·朗缪尔正在潜心电灯泡的研究。欧文·朗缪尔 Irving Langmuir 1881 – 1957做科研的人大都有自己的三板斧,不管问题是什么,都有相对固定的几个套路去处理。朗缪尔是研究物理化学的,两者的交叉天然的就会出现在微观粒子上。朗缪尔研究灯泡也是从这个角度切入,他发现,通过将灯泡填充以惰性气体(如氩气)可以显著延长钨丝的使用寿命。随后,这种在一个系统里面引入额外“粒子”观察反应,成了朗缪尔研究的主要手段。自此,他开始涉足表面化学领域,并发现当分子氢被引入到钨丝灯泡中时,会分解为原子氢并在灯泡表面形成一层厚度仅为一个原子的物质,这相当于早期的“原子制造”。单纯引入额外粒子已经不能满足他,还必须要让它们“听话”。继续沿着这个思路,朗缪尔又研究了等离子体,并成为等离子体的命名人。终于在1932年,由于表面化学的工作,朗缪尔被授予诺贝尔奖。有趣的是,朗缪尔荣获诺贝尔奖的工作是关于油膜的,通过控制界面分子排布来调控宏观界面性质,把粒子层面的应用算是玩明白了。3.2 数学模型与试验二次世界大战促使航空飞行器技术快速发展,轰炸机、运输机尺寸越来越大。在著名的驼峰航线,美军动用大量C-46/C-47运输机飞跃喜马拉雅山向中国重庆抗战后方运输物资。C-46运输机该航线气象条件恶劣,结冰事故频发。根据库塔-儒可夫斯基定理,飞机升力和阻力大小受到外形影响极大,机翼形状一定设计好,最忌“整形”。然而飞机机翼结冰后,外形发生改变,升力下降,阻力增加,造成飞机失速迎角提前,再加复杂地形和气流影响,因此事故不断。现在我们知道,为了应对机翼结冰,需要空气动力学设计和除防冰系统协同。当时美军首先想到的是如何让机翼表面不结冰,即防冰。这一下撞到了朗缪尔的研究领域,“通过控制界面分子排布来调控宏观界面性质”,如果在飞机表面镀膜呢?我想最初军方和朗缪尔的思路可能是这样。就这样朗缪尔开始了飞机结冰的研究。但是对于当时的朗缪尔来说,还有个更关键的问题要搞明白。高空之中,尤其巡航高度,温度通常在零下数十度。在地面这个温度环境中,水都结成了冰,冰撞在飞机上要么弹开要么破碎,怎会积聚呢?除非,它们在零下温度情况下,仍未结冰。虽然已年逾六十,朗缪尔研究手段和物理直觉仍旧狠辣,经过对高空云层的研究,他很快确定,在没有外来粒子(灰尘等)作为凝结核的情况下,云层中漂浮的水滴果然以液态形式存在,这种低于0℃的水滴被称作过冷水滴。飞机机翼表面刚好也是提供了凝结核,才产生了结冰。随后朗缪尔分两路推进研究,一路是理论。根据前面提到的经典力学的分析框架:受力分析、构建平衡方程、求解。朗缪尔把水滴等效成球体,同时考虑流场作用,由于机翼前缘等局部位置可以看成是圆柱的一部分,为此经过简化和处理,建立了圆柱体流场中水滴运动的控制方程。随后又通过一系列推导,由于圆柱体流场求解相对简易,在不借助龙格库塔法的基础上,终于实现了水滴轨迹的艰难求解。另一路则是试验,水滴轨迹和其尺寸关系极大。可是云层中的水滴尺寸究竟是多少?它们又是否一样大,如果不一样大,又有何分布规律?科研就是如此,一环扣一环。朗缪尔根据他推导的圆柱体绕流流场水滴轨迹规律发现,不同直径圆柱体获取的水滴收集量不同,如果设计一个含有多个直径的圆柱体进行测试,就可以根据结冰量的差异反推出云层里的水滴数据。这就是旋转多圆柱法。朗缪尔通过在华盛顿山结冰观测站的大量实验,归纳出了著名的 Langmuir A 到 G 等一系列水滴尺寸分布谱(其中 Langmuir D 谱最为经典)。直到今天,民航局的适航认证标准(如 FAA 的 Appendix C)以及各大结冰仿真软件(LEWICE, FENSAP-ICE)在设置云层过冷水滴的水动力学直径分布(MVD)时,依然在使用 Langmuir 分布 作为标准输入之一。旋转多圆柱法3.3 “科学狂人”经此一役,朗缪尔对大气科学产生了极大兴趣。他开始思考,高空中的过冷水滴没有依附的微粒,就无法打破表面张力结成冰晶,也就无法变重掉落。如果人为向云层提供凝结核,或者强行让它降温产生冰晶,就能触发连锁反应,逼迫整片云下雨。这就是人工降雨的思维起点了。为了验证这个设想,他带领团队开启了试验。引入外来粒子,这是朗缪尔从博士就用惯的套路。朗缪尔助手文森特·谢弗在实验室里用一个改装的家用冷柜来模拟高空低温云层。他在冷柜里吹气,让呼出的水分在低温下形成一片微型的“过冷云雾”。谢弗尝试往里面扔各种粉末(粉尘、滑石粉、土壤)充当凝结核,但过冷云雾毫无反应。直到 1946 年 7 月的一天,纽约天气酷热,实验室里的冷柜温度不够低。为了让冷柜快速降温,谢弗顺手把一大块干冰(固态二氧化碳,温度为 -78.5°C)扔进了冷柜。神迹发生了。就在干冰进入冷柜的瞬间,原本雾蒙蒙的过冷水滴瞬间消失,取而代之的是无数的微小冰晶在光线下闪烁。干冰强大的局部降温能力,强行把局部空气降到了 -40°C 以下。在这种极寒下,过冷水滴不需要任何凝结核,自身就会自发结晶(自发成核)。这时候团队意识到,飞机结冰的原因可能不止凝结核那么简单。过冷水滴本就不稳定,水滴与机翼的碰撞会改变张力状态,同样可以造成结冰的发生。后来朗缪尔的另一位同事伯纳德·冯内古特从化学文献中找到了更优方案。他发现碘化银(AgI)的晶体结构与冰相似,可以作为"冰核"让过冷水滴在其表面冻结,且不需要干冰那样的极低温。不管怎么说,这个冰箱里的试验验证了原理可行性。随后朗缪尔与美军开始策划真正人工降雨试验。1946年11月13日,谢弗乘坐一架小型飞机,从纽约州斯克内克塔迪起飞,在附近寻找到了一片过冷的层云。他从飞机上向约3英里长的云带播撒了约3磅(约1.4公斤)碎干冰,干冰撒入后不到5分钟,云滴就转化为冰晶,几分钟后,雨雪从云层底部倾泻而下,试验成功。现代无人机增雨1947年,在美军的支持下,朗缪尔主持了著名的“卷云计划”。他甚至曾疯狂地派遣 B-17 轰炸机飞入一个正在肆虐的大西洋飓风中,倾倒了 180 磅干冰,试图用人工降雨人为改变飓风的能量分布。虽然那次试验导致飓风诡异地大拐弯并席卷了佐治亚州,尽管有人认为是巧合,但是也引发了巨大的法律和政 治风波,同时也证明了朗缪尔这套理论对大气的恐怖影响力。正所谓“天有不测风云”,自古以来天气从来不以人的意志为转移。但是作为战争中非常关键的要素,任何作战计划都必须考虑天气。近年一直传闻美军已具备相当水平的气象武器,从朗缪尔的气象试验中可以看出,这些传闻绝非空穴来风。真是蝴蝶效应,本身只是一项飞机结冰研究,最后竟然演变成了气象武器的雏形。4. NACA与NASA与此同时,NACA(NASA的前身)在俄亥俄州组建了Lewis结冰研究团队,Rinaldo J. Brun、William Lewis等多位科学家,将朗缪尔圆柱流场结果扩展到翼型流场上。他们发表了多篇具有里程碑意义的 NACA 技术报告(如著名的 NACA Report 1215)。他们利用当时极其原始的模拟计算机(如机械微分分析仪)或者手算,画出了水滴在 NACA 0012 等经典翼型周围的运动轨迹。又过去近30年,此时CFD、数值模拟、计算机技术蓬勃发展,已经具备了大规模求解水滴轨迹的基本能力。NASA组建了 LEWICE 软件团队,试图真正把拉格朗日轨迹法、面元法与现代网格技术完美结合并写成工程软件。所谓的拉格朗日轨迹法就是前面的朗缪尔的粒子追踪技术。根据相对原理,假定机翼不动,远方的水滴粒子撞击过来。根据粒子受力平衡,创建微分方程组。再通过龙格库塔法时间推进求解。水滴在流场中,会受到流场速度的扰动。为简化,采用势流作为控制方程,面元法作为求解技术。面元法的基本思路是在翼面每个单元中心(控制点)布置源强,流场中的每个点的速度由这些单元和对应的源强叠加确定。面元法面元法有个很大好处,就是一旦求解出来翼型表面的源分布,空间中任意一个点的流速,都可以通过一个矩阵乘法完成计算。如果采用更复杂的CFD技术,只能获取流场网格点的速度,这样一来,需要时刻判断水滴所处单元位置,再通过单元点的速度进行水滴位置插值获取速度,如此计算量大增。LEWICE软件于80年代问世,90年代逐步公开,该软件通过一系列精妙设计,可以在数秒内完成冰形计算。LEWICE软件5. 后续我在《我是如何开发飞机结冰软件的》中详述了自己研发飞机结冰软件的历程。人做一个事情前,结果是未知的。越是难、越是耗时的事情,越需要发大愿行之。尽管前人已经做了诸多工作,很多资料也都详实完备,轮到自己而为,仍然走了很多弯路。在多年以后重新整理这个研究的历史,故事之曲折精彩令我瞠目结舌。随着国内民机的发展,属于我们的“结冰”传奇,或许才刚刚开始。来源:静界有限元

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈