- 博客(3)
- 收藏
- 关注
原创 「Python与地震工程」单自由度体系求解之通用微分方程数值解法(基于odeint)
借助Python求解此方程最快捷(编程者的工作量小,并非最快速)的方法就是将方程降阶为状态空间形式调用scipy.integrate.odeint()函数。需要注意的是,对于地震激励,需要把离散的地震波“包装”成一个连续函数才能被odeint识别。单自由度体系分析是结构动力学和地震工程学习的启蒙知识,也是复杂结构动力分析的基石,具有重要的作用和研究价值。本文将探讨如何用Python的科学计算模块来实现单自由度体系的动力响应分析。经典的弹性有阻尼单自由度体系在任意激励P(t)作用下的运动方程可写为。
2023-03-26 11:22:14
203
原创 「Python与地震工程」单自由度体系求解之频域法(傅里叶变换法)
图中, p(t)↔P(ω) 为激励的Fourier变换对;u(t)↔U(ω) 为响应的Fourier变换对;H(iω) 为频域响应函数,对于地震作用下的单自由度体系。「Python与地震工程」单自由度体系求解之频域法(傅里叶变换法)频域法(傅里叶变换法)求解体系响应的基本流程。转载本文请注明出处。
2023-03-26 10:47:19
223
原创 「Python与地震工程」单自由度体系求解之Newmark-β法
「Python与地震工程」单自由度体系求解之Newmark-β法原理Newmark-β法是地震工程领域最经典的逐步积分算法。Nathan Mortimore Newmark推导过程请查阅结构动力学或地震工程学教材,此处仅简单列出逐步递推公式。已知第 步响应,则第 步位移响应可按下式计算:其中则第 i+1 步速度、加速度响应为对于地震激励程序代码import numpy as npimpor
2023-03-07 14:14:49
1075
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人