結(jié)構(gòu)力學(xué)基礎(chǔ)概念：力法：力法與位移法對比分析

結(jié)構(gòu)力學(xué)基礎(chǔ)概念：力法：力法與位移法對比分析1結(jié)構(gòu)力學(xué)基礎(chǔ)概念：力法與位移法對比分析1.1緒論1.1.1結(jié)構(gòu)力學(xué)的基本概念結(jié)構(gòu)力學(xué)是研究結(jié)構(gòu)在各種外力作用下變形、應(yīng)力和穩(wěn)定性的一門學(xué)科。它主要關(guān)注結(jié)構(gòu)的強度、剛度和穩(wěn)定性，是土木工程、機械工程、航空航天工程等領(lǐng)域的基礎(chǔ)。結(jié)構(gòu)力學(xué)中的結(jié)構(gòu)可以是橋梁、建筑物、飛機機翼等，這些結(jié)構(gòu)在設(shè)計時需要考慮其在不同載荷下的響應(yīng)，以確保安全性和功能性。1.1.2力法與位移法的簡介在結(jié)構(gòu)分析中，力法和位移法是兩種基本的分析方法，它們分別從力和位移的角度來求解結(jié)構(gòu)的內(nèi)力和變形。1.1.2.1力法力法，也稱為力矩分配法，是一種基于結(jié)構(gòu)的平衡條件來求解未知力的方法。它適用于超靜定結(jié)構(gòu)，即結(jié)構(gòu)的未知力數(shù)目多于平衡方程數(shù)目的結(jié)構(gòu)。力法的基本思想是，通過引入多余未知力（稱為力法未知量），將超靜定問題轉(zhuǎn)化為靜定問題，然后通過滿足變形協(xié)調(diào)條件來求解這些未知力。1.1.2.2位移法位移法，也稱為位移分析法，是一種基于結(jié)構(gòu)的變形條件來求解未知位移的方法。它同樣適用于超靜定結(jié)構(gòu)，但與力法不同，位移法直接求解結(jié)構(gòu)的位移，然后通過位移與內(nèi)力的關(guān)系來計算內(nèi)力。位移法的基本思想是，通過引入未知位移（稱為位移法未知量），將結(jié)構(gòu)的變形問題轉(zhuǎn)化為一組線性方程，然后求解這些未知位移。1.2力法與位移法的對比分析1.2.1力法的原理與應(yīng)用1.2.1.1原理力法的原理基于結(jié)構(gòu)的平衡條件和變形協(xié)調(diào)條件。在超靜定結(jié)構(gòu)中，除了滿足靜力平衡條件外，還必須滿足變形協(xié)調(diào)條件，即結(jié)構(gòu)在未知力作用下的變形必須與實際的變形相協(xié)調(diào)。力法通過引入多余未知力，將超靜定結(jié)構(gòu)轉(zhuǎn)化為一系列靜定結(jié)構(gòu)，然后通過求解變形協(xié)調(diào)方程來確定這些未知力。1.2.1.2應(yīng)用示例假設(shè)有一個連續(xù)梁，兩端固定，中間有一個支座，梁上作用有均布載荷。這個結(jié)構(gòu)是超靜定的，因為它有三個支座反力，但只有兩個平衡方程（水平和垂直方向的力平衡）。為了使用力法求解，我們可以引入一個多余未知力，假設(shè)為中間支座的豎向反力，然后將結(jié)構(gòu)分解為兩個靜定梁，分別計算它們在均布載荷和假設(shè)反力作用下的內(nèi)力和變形。最后，通過變形協(xié)調(diào)條件，即中間支座處的撓度必須相等，來求解這個多余未知力。1.2.2位移法的原理與應(yīng)用1.2.2.1原理位移法的原理基于結(jié)構(gòu)的變形條件和位移與內(nèi)力的關(guān)系。在超靜定結(jié)構(gòu)中，位移法通過引入未知位移，將結(jié)構(gòu)的變形問題轉(zhuǎn)化為一組線性方程。這些方程描述了結(jié)構(gòu)在各種外力作用下的位移響應(yīng)。一旦未知位移被求解，就可以通過位移與內(nèi)力的關(guān)系來計算結(jié)構(gòu)的內(nèi)力。1.2.2.2應(yīng)用示例考慮一個框架結(jié)構(gòu)，由多個梁和柱組成，受到外部載荷的作用。使用位移法，我們首先確定結(jié)構(gòu)的關(guān)鍵位移，如節(jié)點的水平和豎向位移，然后建立這些位移與結(jié)構(gòu)內(nèi)力之間的關(guān)系。通過求解位移方程，我們可以得到節(jié)點的位移，進而計算出梁和柱的內(nèi)力。在實際計算中，位移法通常與矩陣方法結(jié)合使用，形成有限元分析的基礎(chǔ)。1.2.3力法與位移法的對比1.2.3.1解題思路力法：從力的角度出發(fā)，通過滿足變形協(xié)調(diào)條件來求解未知力。位移法：從位移的角度出發(fā)，通過滿足位移與內(nèi)力的關(guān)系來求解未知位移。1.2.3.2計算復(fù)雜性力法：通常需要計算結(jié)構(gòu)在各種載荷下的變形，計算量可能較大。位移法：直接求解未知位移，然后通過位移計算內(nèi)力，計算過程更為直接，但需要建立復(fù)雜的位移與內(nèi)力關(guān)系。1.2.3.3適用性力法：適用于結(jié)構(gòu)簡單、未知力數(shù)目較少的情況。位移法：適用于結(jié)構(gòu)復(fù)雜、未知位移數(shù)目較多的情況，尤其是現(xiàn)代結(jié)構(gòu)分析軟件中廣泛采用的有限元方法。1.2.3.4精度與可靠性力法：在處理線性問題時精度較高，但對于非線性問題的處理能力有限。位移法：能夠更好地處理非線性問題，如大變形、材料非線性等，因此在現(xiàn)代工程分析中更為常用。1.2.4結(jié)論力法和位移法各有優(yōu)勢，選擇哪種方法取決于結(jié)構(gòu)的復(fù)雜性、未知量的類型以及分析的目的。在實際工程中，位移法由于其計算的直接性和對非線性問題的處理能力，被更廣泛地采用。然而，對于一些特定的結(jié)構(gòu)和問題，力法可能提供更簡單、更直觀的解決方案。1.3代碼示例以下是一個使用Python進行位移法分析的簡單示例，計算一個簡支梁在集中載荷作用下的位移和內(nèi)力。這個例子使用了numpy庫來進行矩陣運算。importnumpyasnp

#定義結(jié)構(gòu)參數(shù)

L=4.0#梁的長度

E=200e9#材料的彈性模量

I=0.1#梁的截面慣性矩

P=10000.0#集中載荷

#定義剛度矩陣

k=(E*I)/(L**3)*np.array([[12,6*L,-12,6*L],

[6*L,4*L**2,-6*L,2*L**2],

[-12,-6*L,12,-6*L],

[6*L,2*L**2,-6*L,4*L**2]])

#定義載荷向量

f=np.array([0,P,0,0])

#定義邊界條件

bc=np.array([1,0,1,0])#1表示固定，0表示自由

#應(yīng)用邊界條件

k_mod=k[np.ix_(bc==0,bc==0)]

f_mod=f[bc==0]

#求解位移

u=np.linalg.solve(k_mod,f_mod)

#計算內(nèi)力

M=np.dot(k,u)

#輸出結(jié)果

print("位移:",u)

print("內(nèi)力:",M)1.3.1代碼解釋定義結(jié)構(gòu)參數(shù)：包括梁的長度L、材料的彈性模量E、梁的截面慣性矩I和作用在梁上的集中載荷P。定義剛度矩陣：k矩陣描述了梁在各種外力作用下的變形特性。它是根據(jù)梁的長度、材料屬性和截面形狀計算出來的。定義載荷向量：f向量表示作用在梁上的外力，這里假設(shè)梁的一端受到集中載荷P的作用。定義邊界條件：bc向量表示梁的支承情況，1表示固定支座，0表示自由端。應(yīng)用邊界條件：通過np.ix_函數(shù)和bc向量，從k矩陣和f向量中去除固定支座的影響，得到修改后的剛度矩陣k_mod和載荷向量f_mod。求解位移：使用np.linalg.solve函數(shù)求解修改后的剛度矩陣k_mod和載荷向量f_mod，得到位移向量u。計算內(nèi)力：通過np.dot函數(shù)計算剛度矩陣k和位移向量u的點積，得到內(nèi)力向量M。輸出結(jié)果：打印出位移向量u和內(nèi)力向量M。這個例子展示了位移法的基本流程，即通過建立剛度矩陣和載荷向量，應(yīng)用邊界條件，求解位移，然后計算內(nèi)力。在實際工程分析中，剛度矩陣和載荷向量的構(gòu)建會更加復(fù)雜，可能需要考慮多個自由度和多種載荷類型。2力法原理2.1力法的基本方程力法，作為結(jié)構(gòu)力學(xué)中解決超靜定結(jié)構(gòu)問題的一種方法，其核心在于將結(jié)構(gòu)的超靜定問題轉(zhuǎn)化為靜定問題，通過求解未知的多余約束力來確定結(jié)構(gòu)的內(nèi)力和變形。力法的基本方程是基于最小勢能原理或最小余能原理建立的，它描述了結(jié)構(gòu)在未知多余約束力作用下的變形協(xié)調(diào)條件。2.1.1最小勢能原理對于彈性結(jié)構(gòu)，當(dāng)結(jié)構(gòu)處于平衡狀態(tài)且滿足變形協(xié)調(diào)條件時，其總勢能為最小?？倓菽苡蓛?nèi)部勢能和外部勢能組成。內(nèi)部勢能是由于結(jié)構(gòu)內(nèi)部變形產(chǎn)生的能量，而外部勢能則是由于外部荷載作用于結(jié)構(gòu)上產(chǎn)生的能量。力法的基本方程就是通過最小化總勢能來求解未知的多余約束力。2.1.2最小余能原理最小余能原理是力法的另一種表述，它基于能量守恒原理。余能是結(jié)構(gòu)在未知多余約束力作用下，內(nèi)部勢能與外部勢能之差。當(dāng)結(jié)構(gòu)達到平衡狀態(tài)且滿足變形協(xié)調(diào)條件時，余能為最小。力法的基本方程通過最小化余能來確定未知的多余約束力。2.1.3方程的建立力法的基本方程通常表示為：δ其中，δij是結(jié)構(gòu)在第i個多余約束力作用下，第j個多余約束力的位移影響系數(shù)；Fi是第i個多余約束力；R2.2力法的解題步驟力法解題步驟主要包括以下幾點：確定超靜定次數(shù)：首先，需要識別結(jié)構(gòu)的超靜定次數(shù)，即結(jié)構(gòu)中多余約束的數(shù)量。選取基本結(jié)構(gòu)：在結(jié)構(gòu)中去除多余約束，形成一個靜定的基本結(jié)構(gòu)。建立基本方程：根據(jù)力法的基本方程，建立未知多余約束力的方程組。求解未知力：通過求解方程組，得到未知的多余約束力。計算內(nèi)力和變形：利用得到的多余約束力，計算結(jié)構(gòu)的內(nèi)力和變形。2.2.1示例：求解一個超靜定梁假設(shè)我們有一個兩端固定的超靜定梁，受到中間點的集中荷載作用。該梁的超靜定次數(shù)為1，即存在一個未知的多余約束力。2.2.1.1步驟1：確定超靜定次數(shù)超靜定次數(shù)為1，因為兩端固定提供了4個約束（2個垂直力，2個彎矩），而靜定梁只需要2個約束（一個垂直力，一個彎矩）。2.2.1.2步驟2：選取基本結(jié)構(gòu)將一端的固定約束改為鉸支，形成一個靜定的基本結(jié)構(gòu)。2.2.1.3步驟3：建立基本方程假設(shè)未知的多余約束力為F，作用于梁的一端。根據(jù)最小勢能原理，建立基本方程：δ其中，δ11是梁在F作用下，梁另一端的位移影響系數(shù)；R2.2.1.4步驟4：求解未知力通過計算δ11和R1，可以求解出2.2.1.5步驟5：計算內(nèi)力和變形利用得到的F，可以計算出梁的內(nèi)力分布和變形情況。通過以上步驟，我們可以系統(tǒng)地應(yīng)用力法來解決超靜定結(jié)構(gòu)問題，確定結(jié)構(gòu)的內(nèi)力和變形。力法在解決復(fù)雜結(jié)構(gòu)問題時，尤其在計算機輔助設(shè)計中，提供了強大的工具和方法。3結(jié)構(gòu)力學(xué)基礎(chǔ)概念：位移法原理3.1位移法的基本方程位移法是結(jié)構(gòu)力學(xué)中一種基于位移的分析方法，它以結(jié)構(gòu)的位移作為基本未知量，通過建立位移與內(nèi)力之間的關(guān)系，進而求解結(jié)構(gòu)的內(nèi)力和位移。位移法的基本方程通常由平衡方程、變形協(xié)調(diào)方程和物理方程組成。3.1.1平衡方程平衡方程描述了結(jié)構(gòu)在任意截面處的內(nèi)力平衡條件。對于一個平面框架結(jié)構(gòu)，平衡方程可以表示為：在節(jié)點處，水平方向和垂直方向的力平衡；在節(jié)點處，力矩平衡。例如，對于一個簡單的平面框架節(jié)點，其平衡方程可以表示為：∑3.1.2變形協(xié)調(diào)方程變形協(xié)調(diào)方程確保了結(jié)構(gòu)各部分之間的位移連續(xù)性。在位移法中，這些方程將結(jié)構(gòu)的位移聯(lián)系起來，確保了結(jié)構(gòu)的整體性和連續(xù)性。3.1.3物理方程物理方程描述了內(nèi)力與位移之間的關(guān)系，通?；诓牧系膽?yīng)力-應(yīng)變關(guān)系。對于線彈性材料，物理方程可以簡化為胡克定律，即內(nèi)力與位移成正比。3.2位移法的解題步驟位移法的解題步驟主要包括以下幾點：確定基本未知量：選擇結(jié)構(gòu)的獨立位移作為基本未知量，通常包括節(jié)點的水平位移、垂直位移和轉(zhuǎn)角位移。建立剛度矩陣：根據(jù)結(jié)構(gòu)的幾何和材料特性，建立結(jié)構(gòu)的剛度矩陣。剛度矩陣反映了結(jié)構(gòu)對位移的抵抗能力。應(yīng)用邊界條件：將結(jié)構(gòu)的邊界條件（如固定支座、鉸支座等）應(yīng)用到剛度矩陣中，形成修改后的剛度矩陣。求解位移：利用修改后的剛度矩陣和荷載向量，通過求解線性方程組得到結(jié)構(gòu)的基本未知位移。計算內(nèi)力：根據(jù)求得的位移，利用物理方程計算結(jié)構(gòu)各部分的內(nèi)力。3.2.1示例：平面框架的位移法分析假設(shè)我們有一個簡單的平面框架，由兩根梁和一個節(jié)點組成，如圖所示：A

|

|

|

|

B

|

|

|

|

C節(jié)點B受到水平荷載P的作用，節(jié)點A和C為固定支座。我們使用位移法來分析這個結(jié)構(gòu)。3.2.1.1步驟1：確定基本未知量在這個例子中，我們選擇節(jié)點B的水平位移作為基本未知量。3.2.1.2步驟2：建立剛度矩陣假設(shè)每根梁的剛度為k，剛度矩陣可以表示為：k但是，由于我們只關(guān)心節(jié)點B的水平位移，因此可以簡化為：k3.2.1.3步驟3：應(yīng)用邊界條件節(jié)點A和C為固定支座，因此它們的位移為零。在剛度矩陣中，這相當(dāng)于將與固定支座相關(guān)的行和列刪除。3.2.1.4步驟4：求解位移利用簡化后的剛度矩陣和荷載向量，我們可以通過以下方程求解節(jié)點B的水平位移：k其中，Δ是節(jié)點B的水平位移，P是作用在節(jié)點B的水平荷載。3.2.1.5步驟5：計算內(nèi)力一旦我們求得了節(jié)點B的水平位移，就可以利用物理方程計算梁AB和BC的內(nèi)力。例如，梁AB的內(nèi)力可以表示為：F3.2.2Python代碼示例下面是一個使用Python和NumPy庫來求解上述平面框架位移的簡單示例：importnumpyasnp

#定義剛度和荷載

k=1000#假設(shè)剛度為1000N/m

P=500#假設(shè)水平荷載為500N

#求解位移

Delta=P/k

#計算內(nèi)力

F_AB=k*Delta

F_BC=-k*Delta

print(f"節(jié)點B的水平位移為:{Delta}m")

print(f"梁AB的內(nèi)力為:{F_AB}N")

print(f"梁BC的內(nèi)力為:{F_BC}N")在這個示例中，我們首先定義了剛度k和荷載P的值。然后，我們使用荷載P除以剛度k來求解節(jié)點B的水平位移Δ。最后，我們利用剛度k和位移Δ來計算梁AB和BC的內(nèi)力。通過位移法，我們可以有效地分析結(jié)構(gòu)在荷載作用下的響應(yīng)，為結(jié)構(gòu)設(shè)計和優(yōu)化提供重要的理論基礎(chǔ)。4力法與位移法對比分析4.1方法的適用范圍在結(jié)構(gòu)力學(xué)中，力法和位移法是兩種基本的分析方法，它們各自適用于不同類型的結(jié)構(gòu)分析。4.1.1力法力法，也稱為間接法，主要適用于超靜定結(jié)構(gòu)的分析。這種方法通過設(shè)定結(jié)構(gòu)的多余未知力作為基本未知量，然后根據(jù)變形協(xié)調(diào)條件建立方程，求解這些未知力。力法的優(yōu)點在于它能夠直接處理結(jié)構(gòu)的超靜定問題，而不需要預(yù)先確定位移。然而，它的缺點是對于復(fù)雜結(jié)構(gòu)，建立和求解方程可能較為繁瑣。4.1.1.1適用范圍示例連續(xù)梁：連續(xù)梁由于存在多個支座，通常具有多個超靜定度，力法可以有效地解決這類問題。框架結(jié)構(gòu)：對于框架結(jié)構(gòu)，尤其是多層多跨框架，力法能夠處理結(jié)構(gòu)的復(fù)雜超靜定情況。4.1.2位移法位移法，也稱為直接法，主要適用于靜定和超靜定結(jié)構(gòu)的分析。這種方法通過設(shè)定結(jié)構(gòu)的位移作為基本未知量，然后根據(jù)平衡條件和變形條件建立方程，求解這些未知位移。位移法的優(yōu)點在于它能夠直接處理結(jié)構(gòu)的變形和位移，對于現(xiàn)代計算機輔助設(shè)計軟件來說，位移法更為適用。然而，它的缺點是對于某些特定的超靜定結(jié)構(gòu)，可能需要較多的位移未知量，從而增加計算的復(fù)雜度。4.1.2.1適用范圍示例桁架結(jié)構(gòu)：桁架結(jié)構(gòu)的分析通常采用位移法，因為桁架的節(jié)點位移是結(jié)構(gòu)分析的關(guān)鍵。有限元分析：在有限元分析中，位移法是主要的分析方法，它能夠處理各種復(fù)雜的結(jié)構(gòu)和材料問題。4.2計算效率與精度分析4.2.1力法力法的計算效率和精度主要取決于結(jié)構(gòu)的超靜定度和未知力的數(shù)量。對于超靜定度較低的結(jié)構(gòu)，力法的計算效率較高，因為未知力的數(shù)量較少。然而，對于超靜定度較高的結(jié)構(gòu)，力法的計算效率會降低，因為需要求解的方程數(shù)量增加，這可能導(dǎo)致計算時間的顯著增加。在精度方面，力法能夠提供較高的精度，因為它直接處理結(jié)構(gòu)的力平衡和變形協(xié)調(diào)條件。4.2.2位移法位移法的計算效率和精度則更多地依賴于位移未知量的數(shù)量和所采用的單元類型。對于復(fù)雜結(jié)構(gòu)，位移法可能需要大量的位移未知量，這會增加計算的復(fù)雜度和時間。然而，現(xiàn)代計算機技術(shù)的發(fā)展使得位移法在計算效率上有了顯著的提升，尤其是在有限元分析中。在精度方面，位移法同樣能夠提供較高的精度，特別是當(dāng)采用高階單元時，它能夠更準(zhǔn)確地模擬結(jié)構(gòu)的變形和應(yīng)力分布。4.2.3對比分析計算效率：對于超靜定度較低的結(jié)構(gòu)，力法可能更為高效；而對于復(fù)雜結(jié)構(gòu)，尤其是需要考慮非線性問題的結(jié)構(gòu)，位移法（尤其是有限元分析）通常更為高效。精度：力法和位移法在精度上各有優(yōu)勢，但總體來說，位移法由于能夠更直接地處理結(jié)構(gòu)的變形，因此在大多數(shù)情況下能夠提供更高的精度。適用性：位移法的適用范圍更廣，它不僅適用于超靜定結(jié)構(gòu)，也適用于靜定結(jié)構(gòu)，以及各種復(fù)雜的結(jié)構(gòu)和材料問題。4.2.4結(jié)論選擇力法還是位移法進行結(jié)構(gòu)分析，應(yīng)根據(jù)結(jié)構(gòu)的類型、超靜定度、所需精度以及計算資源來決定。在實際工程應(yīng)用中，位移法由于其廣泛的適用性和較高的計算效率，通常被優(yōu)先考慮。然而，對于某些特定的超靜定結(jié)構(gòu)，力法可能仍然是一個更優(yōu)的選擇。5實例分析5.1力法解決連續(xù)梁問題5.1.1原理力法，也稱為力矩分配法，是一種解決結(jié)構(gòu)超靜定問題的方法。在連續(xù)梁的情況下，超靜定意味著梁的支座提供的約束力超過了靜力平衡所需的最小約束力。力法通過設(shè)定未知的多余約束力（如支座反力或力矩）為零，然后逐步修正這些力，直到滿足變形協(xié)調(diào)條件，即結(jié)構(gòu)在多余約束力作用下的變形與實際變形相匹配。5.1.2內(nèi)容考慮一個簡單的三跨連續(xù)梁，兩端固定，中間支座為鉸接。假設(shè)梁的剛度為EI，長度為L，中間支座處的荷載為P。此梁有三個多余約束力：兩端的支座反力和中間支座的力矩。5.1.2.1步驟設(shè)定未知的多余約束力：假設(shè)兩端的支座反力為零，中間支座的力矩為零。計算結(jié)構(gòu)在假設(shè)力下的變形：使用梁的彎曲方程和邊界條件，計算梁在假設(shè)力下的撓度和轉(zhuǎn)角。修正多余約束力：根據(jù)計算出的變形與實際變形的差值，修正多余約束力，直到滿足變形協(xié)調(diào)條件。重新計算結(jié)構(gòu)響應(yīng)：使用修正后的多余約束力，重新計算梁的內(nèi)力和變形。5.1.2.2示例假設(shè)我們有以下連續(xù)梁的參數(shù)：梁的剛度EI=1000kN·m2梁的長度L=5m中間支座處的荷載P=100kN我們可以使用Python和NumPy庫來計算修正后的多余約束力：importnumpyasnp

#定義參數(shù)

EI=1000#梁的剛度，單位：kN·m2

L=5#梁的長度，單位：m

P=100#中間支座處的荷載，單位：kN

#計算中間支座的力矩修正值

#假設(shè)兩端支座反力為零，中間支座力矩為零

#初始計算中間支座的撓度

defdeflection_at_midspan(EI,L,P):

"""計算中間支座處的撓度"""

return(P*L**3)/(48*EI)

#計算修正值

defmoment_correction(deflection,L):

"""計算力矩修正值"""

return(deflection*4*EI)/L**2

#初始撓度

initial_deflection=deflection_at_midspan(EI,L,P)

#力矩修正值

moment_correction_value=moment_correction(initial_deflection,L)

print(f"修正后的中間支座力矩為：{moment_correction_value}kN·m")5.1.3解釋在上述代碼中，我們首先定義了梁的剛度、長度和荷載。然后，我們計算了在假設(shè)力（即兩端支座反力和中間支座力矩為零）下的中間支座撓度。接著，我們使用力法的原理，計算了修正中間支座力矩所需的力矩值。這個修正值將用于調(diào)整結(jié)構(gòu)，使其滿足變形協(xié)調(diào)條件。5.2位移法解決框架結(jié)構(gòu)問題5.2.1原理位移法是一種基于結(jié)構(gòu)位移的分析方法，適用于框架結(jié)構(gòu)。它通過設(shè)定結(jié)構(gòu)的關(guān)鍵位移（如節(jié)點的轉(zhuǎn)角和線位移）為未知數(shù)，然后建立這些位移與結(jié)構(gòu)內(nèi)力之間的關(guān)系，最終求解這些未知位移。位移法的核心是利用結(jié)構(gòu)的剛度矩陣，該矩陣描述了結(jié)構(gòu)在單位位移下的內(nèi)力響應(yīng)。5.2.2內(nèi)容考慮一個簡單的兩層框架結(jié)構(gòu)，每層有兩個柱子和一個橫梁。假設(shè)柱子和橫梁的剛度分別為EI和EI/12，層高為H，柱子的寬度為B。此框架有四個關(guān)鍵位移：兩個節(jié)點的轉(zhuǎn)角和兩個節(jié)點的線位移。5.2.2.1步驟建立剛度矩陣：根據(jù)柱子和橫梁的剛度，建立框架的剛度矩陣。應(yīng)用邊界條件：將已知的位移（如固定支座的位移為零）應(yīng)用于剛度矩陣。求解未知位移：使用線性代數(shù)方法求解關(guān)鍵位移。計算內(nèi)力：使用求得的位移，計算框架的內(nèi)力。5.2.2.2示例假設(shè)我們有以下框架結(jié)構(gòu)的參數(shù)：柱子和橫梁的剛度EI=1000kN·m2層高H=3m柱子的寬度B=1m我們可以使用Python和SciPy庫來求解關(guān)鍵位移：fromscipy.linalgimportsolve

importnumpyasnp

#定義參數(shù)

EI=1000#柱子和橫梁的剛度，單位：kN·m2

H=3#層高，單位：m

B=1#柱子的寬度，單位：m

#建立剛度矩陣

#假設(shè)框架為兩層，每層有兩個柱子和一個橫梁

#剛度矩陣為8x8，其中前4行和列對應(yīng)第一層的位移，后4行和列對應(yīng)第二層的位移

#柱子的剛度為EI/H，橫梁的剛度為EI/12/B

defstiffness_matrix(EI,H,B):

"""建立框架的剛度矩陣"""

k_column=EI/H

k_beam=EI/(12*B)

K=np.zeros((8,8))

#柱子的剛度

K[0,0]=K[2,2]=K[4,4]=K[6,6]=2*k_column

K[1,1]=K[3,3]=K[5,5]=K[7,7]=k_column

K[0,2]=K[2,0]=K[4,6]=K[6,4]=-k_column

#橫梁的剛度

K[1,3]=K[3,1]=K[5,7]=K[7,5]=-k_beam

K[1,1]+=k_beam

K[3,3]+=k_beam

K[5,5]+=k_beam

K[7,7]+=k_beam

returnK

#應(yīng)用邊界條件

#假設(shè)框架底部固定，即前兩個位移為零

defapply_boundary_conditions(K):

"""應(yīng)用邊界條件"""

K[0,:]=K[:,0]=0

K[1,:]=K[:,1]=0

K[0,0]=K[1,1]=1

returnK

#求解未知位移

#假設(shè)框架頂部受到水平荷載F=100kN

defsolve_displacements(K,F):

"""求解未知位移"""

#F為荷載向量，前兩個位移已知為零，因此只考慮后六個位移

F=np.array([0,0,0,0,100,0,0,0])

#只保留未知位移對應(yīng)的行和列

K_reduced=K[2:,2:]

F_reduced=F[2:]

#求解未知位移

delta=solve(K_reduced,F_reduced)

#將已知位移和求得的未知位移合并

displacements=np.array([0,0,delta[0],delta[1],delta[2],delta[3],delta[4],delta[5]])

returndisplacements

#計算內(nèi)力

#使用求得的位移，計算框架的內(nèi)力

defcalculate_internal_forces(K,displacements):

"""計算內(nèi)力"""

#內(nèi)力向量為K乘以位移向量

forces=np.dot(K,displacements)

returnforces

#建立剛度矩陣

K=stiffness_matrix(EI,H,B)

#應(yīng)用邊界條件

K=apply_boundary_conditions(K)

#求解未知位移

displacements=solve_displacements(K,100)

#計算內(nèi)力

forces=calculate_internal_forces(K,displacements)

print(f"關(guān)鍵位移為：{displacements}")

print(f"內(nèi)力為：{forces}")5.2.3解釋在上述代碼中，我們首先定義了框架結(jié)構(gòu)的參數(shù)，包括柱子和橫梁的剛度、層高和柱子的寬度。然后，我們建立了框架的剛度矩陣，并應(yīng)用了底部