结构力学数值方法：谱方法与随机结构分析教程.pdf

结构力学数值方法：谱方法与随机结构分析教程

1绪论

1.1结构力学数值方法概述

结构力学数值方法是解决复杂结构力学问题的有效工具，它通过将连续的

物理问题离散化，转化为一系列的代数方程组，从而可以利用计算机进行求解。

这些方法包括有限元法(FEM)、边界元法(BEM)、离散元法(DEM)、谱方法

(SpectralMethod)等。每种方法都有其适用范围和特点，其中谱方法以其高精度

和快速收敛性在解决某些特定问题时表现出色。

1.2谱方法的基本概念

谱方法是一种基于函数展开的数值方法，它将问题的解表示为一组正交基

函数的线性组合。与有限元方法不同，谱方法使用的是全局基函数，这意味着

基函数在整个求解域内都有定义，而不是仅在局部单元内。这种全局性使得谱

方法在处理光滑解时具有极高的精度，收敛速度通常比有限元方法快得多。

1.2.1示例：使用谱方法求解一维波动方程

假设我们有一维波动方程：

22

∂∂

2

=

∂∂

,

其中，是位移，是波速。我们使用谱方法求解此方程。

步骤：选择基函数

1.2.1.11

我们选择Chebyshev多项式作为基函数，因为它们在[-1,1]区间内具有良好

的正交性和光滑性。

1.2.1.2步骤2：离散化

将和离散化，使用Chebyshev-Gauss-Lobatto点进行空间离散化，使用时

间离散化方法（如Runge-Kutta方法）进行时间离散化。

1.2.1.3步骤3：求解

将离散化后的方程组求解，得到位移的数值解。

1

1.2.1.4代码示例

importnumpyasnp

importmatplotlib.pyplotasplt

fromscipy.specialimporteval_chebyt

#参数设置

N=32#Chebyshev多项式的阶数

c=1#波速

L=1#区间长度

T=1#时间长度

dt=0.01#时间步长

#Chebyshev-Gauss-Lobatto点

x=np.cos(np.pi*np.arange(N+1)/N)

#初始条件

u0=np.sin(np.pi*x)

u1=np.sin(2*np.pi*x)

#空间导数矩阵

D=np.zeros((N+1,N+1))

foriinrange(N+1):

forjinrange(N+1):

ifi==j:

D[i,j]=0

elifi==0ori==N:

D[i,j]=-1/(L*(x[i]-x[j])*np.sqrt(1-x[j]**2))

else:

D[i,j]=-1/(L*(x[i]-x[j])*(1-x[j]**2))

#时间离散化

t=np.arange(0,T+dt,dt)

u=np.zeros((len(t),N+1))

u[0,:]=u0

u[1,:]=u1

forninrange(1,len(t)-1):

u[n+1,:]=2*u[n,:]-u[n-1,:]+dt**2*c**2*D@D@u[n,:]

#绘制结果

plt.figure()

plt.plot(x,u[-1,:],o,label=SpectralMethod)

plt.plot(x,np.sin(np.pi*x)*np.cos(c*np.pi*t[-1]),label=AnalyticalSolution)