带孔钢板的单元矩阵法求解
将四个钢板的拉伸转化为四个弹簧之间的组合与作用。已知每个钢板的相对刚度表达式为 K_i=\dfrac{A_iE}{l_i} ,根据之前弹簧的受力分析可以得到弹簧单元的刚度矩阵为:
$$
K^{(i)} = \begin{bmatrix} k^{(i)}_{eq}&-k^{(i)}_{eq}\\ -k^{(i)}_{eq}&k^{(i)}_{eq} \end{bmatrix} \tag{1}
$$
分别写出单元1~4的在总刚度矩阵的位置:
$$
K^{(1G)} = \begin{bmatrix} k_{1}&-k_{1}&0&0\\ -k_{1}&k_{1}&0&0\\ 0&0&0&0\\ 0&0&0&0\\ \end{bmatrix} \begin{matrix} u_1\\ u_2\\ u_3\\ u_4\\ \end{matrix} \tag{2}
$$
$$
K^{(2G)} = \begin{bmatrix} 0&0&0&0\\ 0&k_2&-k_2&0\\ 0&-k_2&k_2&0\\ 0&0&0&0\\ \end{bmatrix} \begin{matrix} u_1\\ u_2\\ u_3\\ u_4\\ \end{matrix} \tag{3}
$$
$$
K^{(3G)} = \begin{bmatrix} 0&0&0&0\\ 0&k_3&-k_3&0\\ 0&-k_3&k_3&0\\ 0&0&0&0\\ \end{bmatrix} \begin{matrix} u_1\\ u_2\\ u_3\\ u_4\\ \end{matrix} \tag{4}
$$
$$
K^{(4G)} = \begin{bmatrix} 0&0&0&0\\ 0&0&0&0\\ 0&0&k_4&-k_4\\ 0&0&-k_4&k_4\\ \end{bmatrix} \begin{matrix} u_1\\ u_2\\ u_3\\ u_4\\ \end{matrix} \tag{5}
$$
将上述的每个单元的刚度组合起来即相加得到总刚度矩阵:
$$
K^{(G)} =K^{(1G)}+K^{(2G)}+K^{(3G)}+K^{(4G)}\\ =\begin{bmatrix} k_{1}&-k_{1}&0&0\\ -k_{1}&k_1+k_2+k_3&-k_2-k_3&0\\ 0&-k_2-k_3&k_2+k_3+k_4&-k_4\\ 0&0&-k_4&k_4\\ \end{bmatrix}
$$