记得很早之前,在做某单位的活儿的时候,来了一个博士,在闲聊时告诉我,“瞬态动力学我们都用欧拉法就够了,什么龙格库塔什么Newmark方法都是浮云!”,令人印象深刻。时间一长,竟也被潜移默化,没有深究个中深意,总觉得欧拉法精度差,那就把时间步取密一点不就好了? 最近需要自己处理一些瞬态动力学过程,遇到了这样那样的问题。不得不老老实实翻书,清华的《计算动力学》一开始就从中心差分法的稳定性开始分析,而且有言:“欧拉法一般只能用于起步或与其他方法配合使用”。这些不得不使我好奇被博士尊为"就够了"的欧拉法未什么这么不受待见,自己用书上的方法推导了一下欧拉法的稳定性,与其他方法相对比,才知道欧拉法精度之差,之难以稳定…… 总结了一下,写了个小程序做了三个算例,算是在这个问题上彻底给自己一个交代,下次遇到学术问题,不要全信他人,哪怕他是博士……呵呵 [attach]49510[/attach] |
GMT+8, 2025-2-18 02:28 , Processed in 0.050553 second(s), 15 queries , Gzip On.
Powered by Discuz! X3.4
Copyright © 2001-2021, Tencent Cloud.