多级弹簧-质量系统瞬态分析（基于Newmark）

最新推荐文章于 2023-10-21 18:45:31 发布

battlestar

最新推荐文章于 2023-10-21 18:45:31 发布

阅读量1.2k

点赞数 4

分类专栏： Python学习 CAE学习文章标签： Python 弹簧阻尼隐式 Newmark

本文链接：https://blog.csdn.net/battlestar/article/details/87926139

版权

该程序主要用于实现多自由度下动态系统的隐式分析，计算模型参见多级弹簧-质量系统的瞬态计算（基于中心差分法）。下面是具体的Python程序。

#coding:utf-8
"""
Created on Tue Nov 20 12:41:40 2018

@author: yujin.wang
"""
import time
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import sys
sys.setrecursionlimit(10000) # 10000 is an example, try with different values
################Force###############
def CenterDiff(t,endtime,m,k,fext,u,delta,du0):
	fint = k.dot(u)
	ddu = m.I.dot((fext - fint))
	du = du0 + ddu * delta
	u = u + du*delta
	t += delta
	t_list.append(t);u1_list.append(u[0,0]);u2_list.append(u[1,0])
	if t > endtime:
		return t_list,u1_list,u2_list
	else:
		return CenterDiff(t,endtime,m,k,fext,u,delta,du)
	
def NewmarkSolver(t,endT,deltaT,TOL,u,du,ddu,k,c,m,iternum):
	t += deltaT
	a1 = m/(beta*deltaT*deltaT) + gamma *c/(beta*deltaT)
	a2 =