[发明专利]基于脉冲响应函数的弹性连接子结构综合方法

权利要求书

1.一种基于脉冲响应函数的弹性连接子结构综合方法，其特征在于，包括以下步骤：

步骤1：通过数值积分方法Newmark法或者试验方法获得每个子结构的脉冲响应函数矩阵H^(s)(t)，其中s表示子结构的序号，H^(s)(t)矩阵的元素表示系统第s个子结构的第j个自由度在脉冲激励下第i个自由度的位移响应；

步骤2：弹性连接件的动力学描述，包括以下步骤：

步骤2.1：系统位移向量的分块：

系统的位移向量用全体子结构的位移向量表示为：

u=u(1)Tu(2)T...u(Ns)TT;---(1)]]>

式中，N_s为系统中子结构的个数，T表示矩阵转置。根据子结构间的连接关系，将子结构的位移向量u^(s)分块为界面自由度和内部自由度

u(s)=uc(s)Tui(s)TT;---(2)]]>

并令：

uc=uc(1)Tuc(2)T...uc(Ns)TT;---(3)]]>

步骤2.2：确定界面位移相容条件：

子结构间所有弹性连接件的自由度u_e全部分布于各个连接件的边界上，位移相容条件要求u_e＝u_c；通过下式求得布尔矩阵B：

Bu=Σs=1NsB(s)u(s)=uc;---(4)]]>

则界面位移相容条件最终表示为：

Σs=1NsB(s)u(s)=ue;---(5)]]>

步骤2.3：确定界面力相容条件：

连接件自由度u_e对全体子结构边界的作用力向量用λ表示，则系统第s个子结构受到的来自连接件的作用力g^(s)为：

g(s)=B(s)Tλ;---(6)]]>

界面力相容条件(即牛顿第三定律)要求子结构边界对连接件自由度u_e的作用力向量为λ_e＝-λ；

步骤2.4：建立连接件的运动方程：

令所有连接件的总体质量、阻尼和刚度矩阵分别为Μ_e、C_e和Κ_e，则所有连接件的运动方程为：

Meu··e+Ceu·e+Keue=λe;---(7)]]>

式中：和u_e分别是连接件自由度的加速度、速度和位移向量；

步骤3：建立子结构的运动方程，包括以下步骤：

步骤3.1：建立时间连续形式的子结构运动方程：

由Duhamel积分可知，系统第s个子结构的运动方程为：

u(s)(t)=∫0tH(s)(t-τ)[f(s)(τ)+g(s)(τ)]dτ;---(8)]]>

式中，f^(s)为子结构所受的外载荷；

步骤3.2：将步骤3.1中的子结构运动方程进行时间离散，得：

un(s)=Σi=0n-1Hn-i(s)(fi(s)+fi+1(s)+B(s)Tλi+B(s)Tλi+1)dt2;---(9)]]>

式中，dt为积分步长，角标代表时刻(如u_n＝u(ndt))；

步骤3.3：利用Newmark方法表示子结构的速度和加速度

u··n(s)=1βdt2(un(s)-un-1(s))-1βdtu·n-1(s)-(12β-1)u··n-1(s);---(10)]]>

u·n(s)=γβdt(un(s)-un-1(s))+(1-γβ)u·n-1(s)+(1-γ2β)u··n-1(s)dt;---(11)]]>

式中，γ和β是Newmark法的无量纲参数；

步骤4：综合求解式(5)、(7)、(9)、(10)和(11)得各个子结构位移响应和子结构间界面力的时间递推公式为：

un(s)=u~n-1(s)+H1(s)B(s)Tλndt2λn=G-1(p~n-1-Σs=1NsK‾eB(s)u~n-1(s));---(12)]]>

式中：

u~n-1(s)=Σi=0n-2Hn-i(s)[fi(s)+fi+1(s)+B(s)T(λi+λi+1)]dt2+H1(s)(fn-1(s)+fn(s)+B(s)Tλn-1)dt2;---(13)]]>

G=dt2K‾eΣs=1NsB(s)H1(s)B(s)T+I;---(14)]]>

K‾e=1βdt2Me+γβdtCe+Ke;---(15)]]>

p~n-1=M~eΣs=1NsB(s)u··n-1(s)+C~eΣs=1NsB(s)u·n-1(s)+K~eΣs=1NsB(s)un-1(s);---(16)]]>

M~e=(12β-1)Me+(γ2β-1)Cedt;---(17)]]>

C~e=1βdtMe+(γβ-1)Ce;---(18)]]>

K~e=1βdt2Me+γβdtCe;---(19)]]>

其中I是单位矩阵；

求出各子结构的位移后，各子结构的加速度及速度响应可以根据式(10)和式(11)分别求出。

2.一种基于脉冲响应函数的弹性连接子结构算法，其特征在于，包括以下步骤：

步骤1：定义基本未知量，分别是系统的第s个子结构的位移函数列向量u^(s)(t)、速度函数列向量加速度函数列向量以及连接件对系统的作用力函数列向量λ(t)；

步骤2：通过Newmark法或者试验方法获得每个子结构的脉冲响应函数矩阵H^(s)(t)，其中s表示子结构的序号；H^(s)(t)矩阵的元素表示系统第s个子结构的第j个自由度在脉冲激励下第i个自由度的位移响应；

步骤3：根据子结构之间的连接关系，求得布尔矩阵B^(s)，B^(s)的作用是从第s个子结构的自由度中筛选出第s个子结构的界面自由度，该界面自由度指该子结构与弹性连接件相连接的自由度；

步骤4：根据系统结构及参数计算得出系统弹性连接件的质量矩阵Μ_e、阻尼矩阵C_e和刚度矩阵Κ_e；

步骤5：初始化基本未知量，即t＝0的初始时刻系统不受外力作用，步骤1中的基本未知量全部为零向量；

步骤6：求解系统在外力f(t)作用下的响应，步骤1中四组基本未知量的时间递推公式如下：

un(s)=u~n-1(s)+H1(s)B(s)Tλndt2u·n(s)=γβdt(un(s)-un-1(s))+(1-γβ)u·n-1(s)+(1-γ2β)u··n-1(s)dtu··n(s)=1βdt2(un(s)-un-1(s))-1βdtu·n-1(s)-(12β-1)u··n-1(s)λn=G-1(p~n-1-Σs=1NsK‾eB(s)u~n-1(s));]]>

其中，dt为确定的积分步长，下角标代表时刻，如β和γ是Newmark法的无量纲参数，G、以及分别由下式确定；

u~n-1(s)=Σi=0n-2Hn-i(s)[fi(s)+fi+1(s)+B(s)T(λi+λi+1)]dt2+H1(s)(fn-1(s)+fn(s)+B(s)Tλn-1)dt2;]]>

G=dt2K‾eΣs=1NsB(s)H1(s)B(s)T+I;]]>

K‾e=1βdt2Me+γβdtCe+Ke;]]>

p~n-1=M~eΣs=1NsB(s)u··n-1(s)+C~eΣs=1NsB(s)u·n-1(s)+K~eΣs=1NsB(s)un-1(s);]]>

其中，

M~e=(12β-1)Me+(γ2β-1)Cedt;]]>

C~e=1βdtMe+(γβ-1)Ce;]]>

K~e=1βdt2Me+γβdtCe.]]>