Dynamic.py 文件源码

python
阅读 23 收藏 0 点赞 0 评论 0

项目:StructEngPy 作者: zhuoju36 项目源码 文件源码
def spectrum_analysis(model:Model.fem_model,n,spec):
    """
    sepctrum analysis\n
    n: number of modes to use\n
    spec: a list of tuples (period,acceleration response)
    """
    freq,mode=eigen_mode(model,n)
    M_=np.dot(mode.T,model.M)
    M_=np.dot(M_,mode)
    K_=np.dot(mode.T,model.K)
    K_=np.dot(K_,mode)
    C_=np.dot(mode.T,model.C)
    C_=np.dot(C_,mode)
    d_=[]
    for (m_,k_,c_) in zip(M_.diag(),K_.diag(),C_.diag()):
        sdof=SDOFSystem(m_,k_)
        T=sdof.omega_d()
        d_.append(np.interp(T,spec[0],spec[1]*m_))
    d=np.dot(d_,mode)
    #CQC
    return d
评论列表
文章目录


问题


面经


文章

微信
公众号

扫码关注公众号