matlab的欧拉方法代码时间积分方案
这是一个有限元研究项目,其目的是针对梁单元在瞬态分析中研究时域中不同数值积分方案的阻尼效应。
特别是对于纵向自由度的冲击载荷。
研究的主要集成方案是Backward-Euler,Newmark方法,居中方案,一些Runge-kutta方案和广义alpha方案。
将所进行的仿真与商业代码Abaqus进行比较,并研究了不同的属性,例如准确性,稳定性,鲁棒性。这些方法和主要的有限元方法都是在Matlab中实现和编程的。
还对图形用户界面进行了编程,以允许用户轻松选择数值积分方法并查看位移结果。
GUI显示了在2m光束的末端进行3500
N冲击模拟的结果,持续时间为1e-4秒。
这些参数可以在
Newmark.m
和
dynamics.m
请享用
!
2022-05-11 16:24:51
88KB
系统开源
1