平面剛架單元相關問題,急求?

1.請用勢能駐值原理推導剛架單元剛度矩陣

2.推導平面剛架單元的形函數矩陣?它的力學意義是什麼?在計算跨間荷載引起的等效節點荷載時,如何應用?


剛架單元如下圖所示,一個單元有兩個節點,每個節點有兩個自由度(撓度w和轉角phi),一共四個自由度,設為a_1, a_2, a_3, a_4,如下圖所示。

1.總勢能Pi =int_0^{l_e}frac{1}{2} EI(w(x))^2dx-int_0^{l_e}qwdx-Q_ia_i

用Rayleigh Ritz, 設approximation function為:

	ilde{w}(x)=N_ia_i=N_1a_1+N_2a_2+N_3a_3+N_4a_4

可推得:

	ilde{w}

其中N_i為形函數。

將近似函數帶入原泛函數Pi。得到:

Pi =int_0^{l_e}frac{1}{2} EI(N_i

這樣Pi就只是a_1, a_2, a_3, a_4的函數了。用駐值原理得到:

frac{partial Pi}{partial a_j}= int_0^{l_e}EIN_i

所以剛度矩陣的元素K_{ij}=N_i

2.推到形函數需要從近似函數需滿足邊界條件推導。邊界條件如下:

	ilde{w}(0)=a_1,	ilde{w}(l_e)=a_3;  	ilde{w}	ilde{w}

設形函數N_i為polynomial, 即N_i(x)=c_0+c_1 x+c_2x^2+c_3x^3,帶入到邊界條件中可求得:

N_1=1-frac{3}{l_e^2}x^2+frac{2}{l_e^3}x^3

N_2=x-frac{2}{l_e}x^2+frac{1}{l_e^2}x^3

N_3=frac{3}{l_e^2}x^2-frac{2}{l_e^3}x^3

N_4=-frac{1}{l_e}x^2+frac{1}{l_e^2}x^3

形函數可應用於將跨間載荷轉化為等效節點載荷,轉化公式如下:

Q_i=int N_i q dx, Q_i為等效節點載荷,q為跨間載荷,N_i為形函數。原理是兩者虛功相等。

舉個例子:

對於這個例子,三個形函數分別為N_1=(x-1)x/2N_2=1-x^3N_3=(1+x)x/2,

帶入公式可以得到Q_1=Q_3=q/3, Q_2=4q/3.


推薦閱讀:

有限元方法的核心思想是什麼?
零基礎如何學習abaqus?
abaqus有限元分析中,運動耦合的控制點到底怎麼理解?
有限元中,分析複雜裝配體,考慮計算機計算能力的基礎之上,不論是不是薄板件,都用實體單元類型,可以嗎?
說一說radioss、optistruct、nastran、abaqus的優缺點。各用於什麼場合?

TAG:有限元分析FEA |