1、譯 文學 院： 船舶與建筑工程 專 業： 船舶與海洋工程 學 號： 0845511233 姓 名： 徐為兵 指導教師： 施興華 江 蘇 科 技 大 學2012年 04 月 02 日非線性框架的形狀靈敏度的可靠性分析 Terje Haukaas , Michael H. Scott摘要 提出了一個統一的,綜合的形狀靈敏度處理方法，其中包括構件節點坐標，橫截面的變化特性，以及球形無彈性框架結構形狀參數。一個新奇之處在于對無論是位移和力為基礎的非線性梁柱行為的有限元構想的幾何不確定性的考慮。形狀靈敏度方程使得不確定性幾何缺陷對結構可靠性評估的相對影響有全面調查。 為此，有限元的可靠性分析用于復雜結構

2、模型，在此情況下一些重要方法可以被應用。這里統一的方法是基于直接的鑒別方法，包括框架的有限元確定性的兼容性關系的變化，以及構件的橫截面幾何特性，以獲得完整的形狀靈敏度方程。形狀靈敏度分析方程，應用在OpenSees的軟件框架。大量的涉及鋼結構，鋼筋混凝土結構的例子表明幾何缺陷可能對結構可靠性評估產生重大影響。2006年愛思唯爾有限公司保留所有權利。關鍵詞：形狀靈敏度直接微分法;幾何缺陷;結構可靠性;梁列;非線性分析;OpenSees。1. 導言 結構工程中的許多應用需要計算關于輸入參數的結構響應數量的梯度,這被稱為響應靈敏度分析。響應靈敏度分析最常應用于問題的優化，如因受到限制最小化結構成本和

3、為系統識別的目的最小化實測和數值響應之間的差異。但在在結構可靠性分析中第一階和第二階可靠性方法提出另一種優化問題（一階可靠性方法和二階可靠性方法）。這些方法依賴于“最有可能的失效點”的測定，這是隨機變量約束優化問題的解決方案。作為一種副產品，一階可靠性方法分析為根據其對結構可靠性的相對影響來對不確定參數排名提供了重要方法。重要性在于此方法解決了由于不同的單位個體反應靈敏度無法直接比較的問題。它也強調，響應靈敏度作為一個獨立產品在結構設計上是有用的，因為它們表明結構的靈敏度響應了大量設計參數改變。響應靈敏度在基于梯度的優化算法的應用中有三個要求：效率，準確性，一致性。響應靈敏度的高效的計算需要基

4、于梯度算法可與自由梯度方法競爭，從計算成本方面來說包括響應面方法。在高維空間進行優化變量時這是特別重要的，在這種情況下，通過有限的差異反復運行以獲得梯度，是不可行的。在優化算法中為了避免收斂問題，精度是必要的，由于梯度上即使是很小的誤差可能執行不良。一致性要求源于利用近似的數值模型以獲得結構的響應。因為它在優化問題上采用近似響應，它是所要求的近似響應的梯度。因此，追求“精確”梯度理論邊界值問題是沒有意義的。事實上，這將導致函數值和梯度之間的不一致。 有兩種方法可用來獲取響應的靈敏度：有限差分法（FDMS）和直接分化法（DDM）。有限差分方法采用被擾亂的參數值進行結構分析的重新運行來估計響應靈敏

5、度。因此，它是一個計算效率低下的方法。此外，FDMS的準確性遭受關注。對于非線性問題選擇參數擾動值不是一個簡單的任務。如果擾動太小了，會引入四舍五入的誤差;而如果擾動過大，局部的非線性可能會導致對靈敏度的不準確估計。然而，FDMS滿足一致性的要求，因為它是在有限差分方程中采用的近似響應。DDM提供一個有吸引力的替代物來替代FDMS。在獲得和執行分析靈敏度方程的一次性成本內應獲得有限元響應算法，高效，準確，一致的反應靈敏度。沒有有限差分計算是發生在DDM內的，相反，響應方程與普通的響應計算一起經計算機分析鑒別和實施。一些研究人員對發展分析方程作出了貢獻，包括Choiand Santos2， Ts

6、ay and Arora24，Liu and Der Kiureghian15，Zhang and Der Kiureghian25，Kleiber et al. 13Conte et al. 3，Roth and Grigoriu20， Scott et al. 22和Haukaas and Der Kiureghian9。 DDM的效率比FDMS要高，因為重復運行響應分析是沒有必要的。準確性可保證和響應響應在同一精度，因為相同的方程的解被用來獲得響應和響應靈敏度。一致性是他們已經在空間和時間上被有限元程序離散之后后， “通過區分響應方程實現的，因此，DDM是首選的計算響應靈敏度的方法。在本

7、文中，OpenSees的軟件框架16被擴大和應用于形狀靈敏度分析。OpenSees（地震工程模擬開放系統）是一個專門為地震工程分析而開發的開放源代碼的，面向對象的，通用的有限元代碼。OpenSees一開始是作為太平洋地震工程技術研究中心（PEER）試驗臺模擬的計算平臺，并一直被國家科學基金會贊助的 George E. Brown Jr. Network小網絡地震工程模擬（NEES）采用。由 Haukaas and Der Kiureghian8 對Open-Sees 的響應靈敏度和可靠性分析能力進行研發擴展這讓分析師可以表征輸入的參數作為隨機變量并計算結構響應事件的概率。這就是所謂的有限元可靠

8、性分析（FERA），與所謂的重點對響應的二次矩統計的隨機有限元方法不同。相比之下，可靠性的分析，特別是FERA適合于計算罕見的響應事件概率。這滿足了在性能為主的工程方面不斷增長的需求，以以一定概率的方式評估在激烈荷的罕見事件中結構行為。在OpenSees中響應靈敏度的實現是基于DDM的。實現分為了一個總體框架和具體對象的實現。后者反映了OpenSees被組織成元素，區域和實物的事實。靈敏度計算的框架，以及對選定的元素，區域和實物具體的實現已經到位。這包括與節點坐標相關的靈敏度，而以往的研究表明可能是結構可靠性不確定性的一個重要來源，尤其是考慮到非線性結構的行為9。在本文中，DDM形狀靈敏度方程

9、包括的相關的反應靈敏度有：（1）節點坐標，（2）全球性的結構或構件的形狀參數，（3）尺寸和纖維離散橫截面細節。特別重要的是在位移和力為基礎的構想中統一的梁柱元素的形狀靈敏度方程的發展。提出包含各級形狀靈敏度的梯度計算（結構，元素和截面）而且他們在OpenSees的應用允許在一個可靠性分析中包含廣泛的不確定的幾何缺陷。涉及鋼結構和鋼筋混凝土結構的兩個數值例子讓我們進一步了解了不確定幾何缺陷相對其他不確定結構特性的重要性。2.響應靈敏度在有限元可靠性分析中的應用本文對響應靈敏度的需要源于結構可靠性分析。為了獲得確切的可靠性評估，采用復雜的結構模型模擬結構性能。顯然，由于模型和輸入參數中的不確定性這

10、種預測只能在一定的概率意義上進行。這促使利用FERA獲得響應事件的概率預測。事實上，新興的基于性能的工程方法被設想應用在可靠性框架6,17。FERA的主要目標是獲得罕見的用戶定義的響應事件的性能函數的概率。分析的一個重要副產品是參數的重要性并根據它們的相對重要性設法排名變量。響應靈敏度就如下面的描述一樣代表分析中的一個重要組成部分。單一執行的情況下的可靠性問題函數是用多重積分表示的其中p是待求的概率，g是為了求解概率識別響應事件的執行函數，f（H）是隨機變量的聯合概率密度函數，這是收集在矢量H中的。在FERA中執行函數在響應數量U = U（H）方面被同有限有限元分析區分開來。隨機變量通常用邊緣

11、概率分布和相關系數指定。方程（1）的解析解是得不到的;然而，如FORM和SORM抽樣技術等方法提供了近似的解。特別需要關注的是在有限元可靠性分析中FORM，其次是高效的重要抽樣可用來糾正潛在的非線性。這是分析策略是有益的，因為它需要對執行函數相對較少的評估，實際上較少的執行有限元分析。此外，FORM分析產生參數重要性的方式，會被用在這個文章中。在FORM中， 方程 (1)中的積分邊界 g = 0在不相關的標準正態隨機變量被轉換的空間Y = Y(H)中被用超平面近似。非線性執行函數中，近似的理想點是界面上的g = 0的點即最接近Y空間的原點的點。這一點，被稱為最可能的失效點（MPP）和用Y*表示

12、 ，是約束優化問題的解解決這個優化問題最有效的算法利用執行函數的梯度，即。執行函數的微分的鏈式法則服從導數是容易得到的，因為g是響應量U的一個簡單的代數函數。意味著需要計算響應梯度，這是本文重點，矩陣是概率變換的雅可比行列式。這項工作要用到納塔夫變換14，因為所需的雅可比行列式是已經在OpenSees提供的。這個變換是Rosenblatt變換19,10一個有吸引力的替代品，因為它允許各種概率分布類型更廣泛的的相關值。當每次對執行函數進行評估時可靠性算法和有限元分析模塊之間的交流組成利用實現隨機變量H和返回值U更新有限元模型與 。已知經測定的Y*，概率p形成取決于其中U是標準正態分布的積分函數，

13、是在 FORM中定義的可靠性指數為。抽樣分布中心在Y*的重要性抽樣可能隨后被執行，因為與抽樣中心在隨機變量的平均中實現的蒙特卡羅抽樣法相比，它是一種高效的方案。本文的發展允許除隨機材料和負載變量之外在節點坐標和構件截面幾何特性方面的不完善的表征作為隨機變量。特別要關注的是對這些與其他不確定性來源相關的變量的重要性的調查。顯然，由于變量屬于不同維度向量的分量不能用于此目的。相反，FORM的重要性措施得到了運用，在此法中分量有統一的維度。這些方法的基礎是由Hohenbichler和Rackwitz11和Bjerager和Krenk1提出。FERA中重要性方法的應用是由Haukaas 和Der Ki

14、ureghian 9提出，下面的向量對隨機變量排名進行排名： 是由方程（3），是概率變換在MPP的雅可比行列式，意味著參數矩陣中的每個元素的平方根。有人指出，會減少到這是在可靠性分析中隨機變量之間沒有表現出相關性時眾所周知的“alpha-vector”的規模化版本。此外，-vector可以通常由測量得到，因此。中的元素被解釋為單個隨機變量對結構的可靠性的貢獻。而且，負（正）分量表示相應的隨機變量起到類似電阻（負載）變量的的作用。3.高次方的響應靈敏度方程 方程（3）中執行函數的梯度意味著需要計算的結構響應的靈敏度。為了通過DDM計算響應靈敏度的，控制結構響應方程被分化了。為了描述結構模型的材料

15、，幾何特性，負荷參數參數H的設置，全球靜態平衡方程有這樣的形式 其中Pr是結構內部抵抗力的矢量。內力可能明確地取決于H，也隱含通過位移響應矢量U。矢量Pf代表作用于結構的外部載荷。在方程（6）中忽略了慣性力和阻尼力，因為本文認為動態平衡的影響與元素和部分構想是獨立的。從（6）式到動態平衡和相應響應靈敏度的計算的推廣 是簡單的5。為了在結構水平上構造響應靈敏度方程，（6）式與從向量H選擇任何相關的參數有區別的： 微分的鏈式法則應考慮到Pr對顯性和隱性的依賴。矢量是反作用力對固定位移的條件導數。矢量外部負載的導數，只有當參數代表作用于結構的負載時，它才是非零的。重排式。方程（7）為節點響應靈敏度1

16、3給出了一個方程的線性系統： 這里是切線剛度矩陣。對于向量H中的每個參數，方程（8）中等式右側和方程組的解的集合給出了相應的節點響應靈敏度矢量。方程（8）的線性形式和切線剛度矩陣的再用有利于提高DDM計算上的效率。 方程（8）式中的靈敏度方程。需要力向量p對于每個元素的導數的集合，對于固定節點位移：這里是集合程序， u是元素位移矢量。有人指出該集合應所有包含非彈性材料響應的元素，不管是否對應個別元素的參數。存在于向量U和u之間的一到一的映射是獨立于的，因此不需要微分。利用DDM計算響應靈敏度的關鍵是從結構模型的每個元素中獲得。 的計算取決于非彈性材料的響應的元素方程式，正如接下來的部分的描述：

17、4.元素的平衡和運動方程的概述 平衡和協調方程的概述需要合成元素貢獻的抵御力矢量如圖1所示，展示的是Filippou和Fenves4的結論。圖1左側的方程代表不同層次之間的內力平衡而右側的代表每個級別變形之間兼容性的關系。圖1的中間列表示在每個級別鏈接力和變形的要素關系。在全球體系中，一個元素抵抗力和節點位移都分別包含在向量p和u中。在一個基本的系統中構想梁柱元素，與剛體位移模式無關，其中向量v收集元素變形和向量q收集相應的力。對于小位移，節點位移和元素變形之間兼容性的關系是線性的，正如矩陣向量產品描述的那樣。矩陣A描述力和位移在球形坐標系統和基本坐標系統之間的轉變的定義如球形坐標系中元素的定

18、位規定的那樣。 相對梯度關系給出基本系統和球形系統之間中的元素力的等式關系。 圖1所示的和b矩陣在元素的基本系統內描述平衡和兼容性關系，這在本文后面討論。部分兼容性矩陣，將部分變形和在橫截面區的任何點的材料應變聯系起來。 圖1 梁柱元素的主方程式5. 靈敏度推導的統一方法 如上式（9），通過DDM進行響應靈敏度分析要求對結構模型中每個元素的p對于固定的u的條件導數進行計算。當參數代表了材料的一種性能，而且采用標準的有限元方程時，在給定位移和應變13的條件下為了確定，直接求平衡方程和材料本構關系是足夠的。此外，Haukaas和Der kiureghian7承認，為獲得正確的靈敏度結果，代表幾何參

19、數時，給定位移的條件下，還必須求出運動學關系。 雖然上述方法提供了基于位移的有限元素的正確的結果，但是它不適用于基于力的方程，其中未指定位移場。在此采用的獲得形狀靈敏度方程的方法是合并元素力向量的全導數與連接圖1中每個級別的平衡和兼容性關系的導數。在代表的材料參數的情況下，最初是由Scott等 22發展的，本文中導數被擴展到代表幾何參數的情形。這方法基本上由圖1中每個級別相結合的四個方程組成。總的來說，這些方程是：（一）平衡方程的導數。（二）運動兼容性方程的導數。（三）力矢量包括本構關系的在目前的水平的全導數。（四）力矢量在下一水平連接較低水平的全導數。 這種方法被用來證明球形元素的力的條件導

20、數。球形坐標系統和基本坐標系統的平衡和兼容性關系的導數式是：同時，球形和基本系統下力的全導數式是：首先，式(10a)通過代入球形和基本力被按照方程式定義展開。式(10c) 和 (10d)分別對應于：其中和。然后將此式與(10b)合并然后，元素剛度矩陣的定義，通過消去和相關的項來簡化表達式，元素力的條件導數的最終表達式是式(13)右邊的第一和第三項表示在元素力的作用下元素長度變化的靈敏度，正如矩陣所表述的那樣。這個矩陣直接計算元素的方向余弦并且對于那些不表示元素的節點坐標的任何參數來說它是等于零的。第二項表示與下一水平的連接，基本坐標系中計算是必需的。條件導數依賴于非線性材料行為的元素構想，正如

21、接下來的基于位移和基于力的梁柱構想部分所展示的那樣。6.基于位移的元素的梯度計算 基于位移的梁柱元素的構想遵循標準有限元分析程序，其中元素位移場是從節點位移12,26篡改過來的。兼容的部分變形是從元素變形改過來的，在截面和基本力之間有弱平衡，如圖1表明的那樣。線性軸向位移和立體的 在基于位移的元素的實現中，在標準化的定義域內= 1, 1。用數字的集合來對平衡關系進行評估（尤其是兩點高斯集合）。x和的定義域的轉換是按照式進行的，其中L是元素的長度。因此，轉換的雅克比式是而且平衡關系的數字集合是 其中是集合點的個數，是第i個點的下標，是相關聯的集合重量。對于高斯集合來說，點和重量都是確定的，因此它

22、們的導數都將是等于零。 通過將x和轉換式替換到是可以對式(15)進一步簡化的，此時平衡關系變為：其中標準的插入矩陣為式(17)中的形式與元素的長度是無關的，因此對于所有參數來說都將等于零。為了根據前文所建立的程序得到條件導數，關于的元素平衡和兼容性關系被分化了：基本的和部分的力矢量的全導數是 在式（18a）中是由和對于的獨立性組成的。類似于球形系統中的元素的力條件導數的推導，式 （18c）和（18d）的基本和部分力的導數分別被插入式（18a）。然后，將式（18b）與結果表達式和基本系統中單元剛度矩陣的定義合并，可將涉及的項消除。這個過程在位移為基礎的構想中產生了相對于條件導數的下列方程：其中是

23、部分剛度矩陣。向量，是根據截面的基本響應的梯度計算的。這在本文后面將進行討論。導數，是式（14）中插補矩陣的一個簡單的縮放，其中只有依賴于的項是1/L的共同的因素：元素長度的導數，是通過對描述元素方向的方向余弦進行微分獲得的。當不對應于其中一個元素節點的一個坐標時此導數等于零。通過將式（20）代入式（19）可進一步簡化基本力的條件導數： 方程（20）和（21）的形式是特定于線性軸和三次埃爾米特橫向位移領域的假設，因為它可能將插值矩陣標準化為式（17）中的。這個標準化使得產生一個有效的數值執行結果，因為在基本系統中它只需要式（21）中的一項來說明形狀靈敏度。涉及的導數的項將出現在條件導數中，當假

24、設位移場沒有被元素長度標準化時。 7.基于力要素的梯度計算 在力的方程式中23，積分形式所表現的是協調關系而非平衡關系。平衡和協調方程是s = bq和v=，分別表示在圖1中。根據基本系統中的靜力平衡，矩陣b從元素結束力內插部分力： 由于通過式（22）中的元素長度對x坐標標準化，力插值矩陣不不依賴于任何參數，因此，其導數等于零。協調關系是靠數值積分來進行評估的Neuenhofer和Filippou18為基于力的元素開發了一種狀態確定的程序，它繞過必須滿足式（23）中協調關系所需要的內部迭代，而沿元素在每節執行平衡。GaussLobatto正交基于力的元素的實現的標準，因為它將積分點置于元素的端點

25、，在沒有單元載荷時此處的彎矩是最大的。在基于力的方程式中為了確立條件導數，將平衡和協調關系對進行求導：基本力矢量和截面力矢量的全導數是式（24a）和（24b）利用了和b相對于的獨立性。合并方程得到的一個表達式的過程概念上類似基于位移的列方程式，但稍微有點復雜。首先，式（24c）和（24d）的基本和截面力的導數被插入式（24a）中。重新整理最終表達式，給出了的方程： 其中，是截面撓度矩陣。將這個表達式與式（24b）合并，然后從元素撓度矩陣的定義出發，和恒等式，消去與相關的項，條件導數的最終表達式為 從式（23）的元素相容性關系出發，可進一步簡化條件導數：重要的是要注意式（21）和（27）對于基于

26、位移和力方程組功能的等價，分別的，其中只有一項是用來說明元素基本系統的形狀靈敏度。當剪切力的插值出現時，另外的涉及的項會在式（27）中出現；然而，通常情況下，忽略剪切效應，方程（27）的形式造成導致一個有效的數值實施，因為是零。切面力的條件導數，取決于本構模型在沿元素的各積分點，如下面的部分所討論的。8.在截面和物質水平上的梯度計算 沿元素的的每一個截面響應在截面變形方面被定義，e和相應的截面力，或應力合力，s，如圖1所示。不管元素方程組，截面力的條件導數，需要確定元素對球形抵抗載荷力向量梯度的貢獻，如式（21）和（27）。這可以通過以下方式獲得導數一個封閉形式的壓力，由此產生的可塑性直接微分

27、關系或材料的數值積分超過截面的應力。在前者情況下，問題簡化為求一個截面本構律的解析梯度方程，而在后一種情況下，截面的兼容和平衡方程，和，分別地，必須進行求解。截面的積分是通過一個用戶定義的纖維數量用數值積分評估，：其中假想水平截面的截面兼容性矩陣仍保持水平，包含纖維的位置，。在第i個光纖位置的材料應力是，是相應的纖維面積。和前面的章節的程序相同，條件截面力的條件導數是通過對截面的平衡和協調方程對微分獲得的： 截面力矢量和材料應力的全導數是： 其中，是第i個纖維位置的應變。截面力和材料應力的替代物從式（29c）和（29d）求導而得，分別代入（29a）隨后結合式（29a）和（29b），其次通過截面

28、剛度矩陣消除關于的項，給出截面力的條件導數其中是材料剛度。導數和分別對應位在第i個截面纖維置和大小的變化。對于構件橫截面的尺寸和細節來說，這些項在計算結構響應敏感度是很重要的。關于光纖離散寬的法蘭部分和的計算在附錄I中有說明。類似的方法進行鋼筋混凝土部分中顯得尤為重要以確定鋼筋的位置和數額的響應靈敏度，如在在本文末尾的數值例子中說明。 在計算結構響應靈敏度的剩余任務是獲得各處纖維位置的物質響應。許多參考文獻提供特定本構關系的響應靈敏度的詳細推導，包括Zhang and Der Kiureghian25，Kleiber等13，Roth和Grigoriu20，和Haukaas and Der Ki

29、ureghian8。OpenSees的實現包括的彈性模型和為無彈性行為的鋼和混凝土許多數字單軸材料模型。9.關于球形參數的反應梯度 在上一節的推導包括元素節點坐標的反應靈敏度。然而，它經常獲得敏感全球的形狀參數，其中包括離散到幾個有限的元素框架構件的終點坐標，以及描述結構全球幾何缺陷的參數。因此，一個全球性的形狀參數的攝動擾動幾個節點的坐標。 為了獲得與一個全球形狀參數相關的響應靈敏度，需要獨立的節點坐標之間明確的關系，以及形狀參數。對于一個形狀參數的情況下，記為，這可能與構件節點坐標重合，這種關系表示為其中是取決于球形參數的節點坐標的集合。 比如， 作為一個例子，在2維空間中假想有一個有n個

30、節點的多元素直線構件。該構件的坐標兩端均表示為和。對于所求的反應靈敏度這四個參數，任何一個都可以是參數的。為了建立式（31）明確的形式，其他構件的節點坐標用端點坐標表示： 其中n1是被離散的構件的標號。 為了獲得球形參數的響應靈敏度，對節點響應靈敏度應用微分的鏈式法則： 矩陣在前面部分是從導數得到的，是對式(31)求導得到的。以式（32）為例，當響應量u的靈敏度方面的要求構件，節點1的x坐標式（33）變為當表示一個球形結構的形狀缺陷時，矢量 是通過計算每個節點處由于單位球形缺陷而產生的缺陷建立的。因此，被解釋為一個影響系數的矢量。球形結構的形狀缺陷，以及多元構件的端點坐標是不確定量。因此，在結

31、構可靠性分析中它們應被視作隨機變量，就像接下來的鋼筋和混凝土框架的數值例子說明的那樣。10.數值算例 這篇文章推導出的響應靈敏度方程已經被應用在OpenSees中而且已經被有線差異計算所證實。為了研究與其他不確定的結構參數相關的不確定幾何缺陷的重要性，靜力推覆可靠性分析被應用在兩種結構中：鋼架結構和鋼筋混凝土框架結構。靜力推覆分析是地震工程中能力評估流行的分析方法，因為它的評估結構以達到目標位移的需求的能力。因此，本文的分析是為了著重介紹結構能力的不確定性評估。 10.1.鋼結構的有限元可靠性分析 圖2中第一個結構認為是三開間，三層鋼結構框架。每個構件都離散成四個基于位移的元素來代表沿構件長度

32、曲率的非線性分布。一纖維離散代表框架構件的寬法蘭鋼截面的響應，如圖2所示。每個法蘭有兩個纖維和網格有十個纖維。每個應力應變行為都用圖5a中所示的單軸材料模型表示，它有個3個材料參數：(1)彈性模數E；（2）屈服強度；(3)二階剛度比。 結構模型的所有材料和幾何參數都被認為是不確定的。各個構件的截面尺寸d，和被模型化平均值分別為250，2，250和20毫米，和2%的變異系數(cov)。每個的彈性模數E，是一個對數隨機變量，其平均值為200,000MPa，5% cov，而且和其他構件的彈性模數的相關系數為0.6。每個構件鋼材的屈服強度，是一個均值為300MP構件a的對數正態分布隨機變量，10cov

33、，與其他構件的相關系數為0.6。每個構件的剛度比是均值為0.02的對數正態分布隨機變量，10cov，與其他構件的相關系數為0.6。總共，21個構件×7個參數=147隨機變量表示材料和截面的幾何參數。此外，16個連接節點每個節點的兩個坐標被認為是不相關的正態隨機變量，給出總計179個隨機變量。縱坐標的標準偏差是10毫米，而水平坐標的標準偏差在10毫米到25毫米之間。該高度的標準偏差的變化是由于幾何缺陷的存在而造成建筑潛在的一個全球性的搖擺。應當注意到，假定每個構件都保持筆直，內部構件節點的位置由式（32）描述。為了調查幾何缺陷的重要性相對于其他結構參數確定載荷應用于結構。外部連接的重力

34、荷載為50kN，內部的為100kN。橫向的載荷隨高度變化，最大值為400kN。圖2 鋼框架結構。顯示節點號和元素號（括號內） 為了評估結構需求的側向位移，為了獲得總漂移最大值超過3的概率需要進行有限元可靠性分析。為了這次響應事件功能函數是表1 鋼結構框架中最重要的25個參數的排名 其中節點13的水平位移。當獲得它的功能函數及其梯度之后，可得到FORM中的MPP。即，隨機變量的不同的實現五個循環的有限元分析是必需的。由此產生的可靠性指標，是2.01，這意味著有0.022的概率超過3的目標漂移。圖3顯示的位移響應在節點13的響應與平均實現和隨機的變量的MPP實現。正如所料，平均觀察為溫和的非線性，

35、同時顯著屈服發生在失效的位移。本文特別關注的是根據方程（5）中的的重要方法排名隨機變量。179個隨機變量中最重要的 25個變量如表1所示。列構件的屈服強度，構件12除外，位居第34位，排名在最重要的參數內。這強調了一個事實，在這個例子中，屈服發生在柱中。網格深度，d，和法蘭寬度，表1中從12到25的構件，表明缺陷在幾何尺寸中的重要性，盡管cov只有2%。右列中的節點4,8,和12的水平坐標排在50個最重要的隨機變量中，因此表明缺陷在結構的球形幾何中相對較高的重要性。 10.2鋼筋混凝土結構的有限元可靠性分析調查對鋼架的觀察對鋼筋混凝土是否也是有效的。為此，圖4中對兩開間，兩層的結構進行了可靠性

36、分析。構件以一個單一的基于力的元素代表捕獲的曲率變化，由于軸向和力矩力量的相互作用。 如圖4中橫截面的每個纖維是仿照單軸材料模型。混凝土材料纖維的核心和蓋是用單軸模型與改良的KentPark骨干曲線描述的21并用零張力線性卸載/重新加載，分別如圖5b和c所示，。在前面的例子中使用的雙線性模型代表對加固鋼的應力應變響應。圖3 平均載荷位移響應和鋼結構隨機變量MPP的實現 圖4 鋼筋混凝土框架結構。顯示節點編號和元素編號（括號內）。認為所有的材料和幾何參數是不確定的。的橫截面尺寸的b和h是不相關的正常隨機變量且圖4中給出了平均值和5的cov。每個鋼筋的面積，是一種平均值如圖4的正態隨機變量。 和2

37、的cov。每個構件涉及構件到具體的厚度是一個平均厚度為75毫米和10cov的正態隨機變量。 圖5 梁和柱蓋區域（a）鋼鐵（b）在自由混凝土，以及柱核心區域的（c）有約束混凝土材料模型。每個的加強筋的彈性模量E，是一個對數正態分布隨機變量的均值為20MPa，5cov，與其他構件E的相關系數0.6。加強筋屈服強度，是對數正態分布隨機變量的每個構件構件平均420MPa，10的cov，與其他構件的相關系數0.6。剛度比，每個構件的加強筋是一個對數正態分布的隨機平均0.05，10的cov，和其他構件的的相關系數0.6。圖5b和c中的所有混凝土材料參數是對數正態分布隨機變量，平均值如圖所示和10的cov。

38、混凝土的強度參數，和，相關性為0.6，以及相應的應變，和。每個節點的兩個坐標被視為不相關的正態隨機變量。垂直坐標被分配了一個10毫米的標準偏差，而橫向坐標分配標準從10毫米到至多20毫米的偏差。這個總也有142個隨機變量可靠性分析。重力荷載850kN=1700千牛，而橫向載荷和。所有應用負載是確定的。這個例子中的執行函數被定義為確定最大位移將超過2偏離的概率：其中是節點7的側向位移。一個FORM有限元的可靠性分析在九個迭代收斂可靠性指數2.76與相應的概率0.0029。隨機變量的最可能失效點即發生重大非線性處，如圖6所示。表2中列出了25個最重要的隨機變量。值得注意的是，它被發現的構件的橫截面

39、的高度h毗鄰的內在聯系是最重要的參數，其次是相同構件的加強筋的屈服強度。這些結果表明，幾何缺陷在鋼筋混凝土結構的可靠性的重大影響。人們還注意到，蓋混凝土的深度也有相當的重要性。這一發現可進一步證明在鋼筋混凝土結構蓋的大量分散的調查。正如前面的例子，用水平節點坐標表示的全球結構形狀缺陷在25個最重要的隨機變量中也淫節點2，3，5也有相當的重要性。表2 在鋼筋混凝土例子中25個最重要的參數排名11. 總結 一個全面和統一的通過直接微分的方法來處理響應靈敏度方程被發現了。幾何，材料和負荷參數都包括在雙方的基于力和基于位移的非彈性梁柱響應方程的建立。全球形狀參數的靈敏度方程解釋了被離散成多個有限元素的

40、結構構件的幾何缺陷。分析方程已經被應用，并被OpenSees軟件驗證。三段式三層鋼框架的有限元可靠性分析展示橫斷面構件尺寸，特別是截面深度和法蘭寬度，有很高的重要性。對于兩段，兩層鋼筋混凝土結構，截面深度，以及混凝土蓋的厚度，都是與材料參數相關的重要參數。在每個在每個例子中，相對于其他結構參數的幾何缺陷的重要性排名表示著不確定的幾何參數對可靠性評估的重要影響，即使分散在概率分布中很小。致謝這個工作和在OpenSees中相應的軟件應用在國家科學基金(NSF)的EEC-9701568號文件和加拿大自然科學和工程研究協會(NSERC)的發現授權的規定下已由太平洋地震工程研究中心提供支持。他們的支持是非常值得感謝的。 附錄 1 寬法蘭截面形狀靈敏度方程 正如方程（30）表明的那樣，為了解釋定義一個離散纖維橫截面的幾何參數中的變化有必要計算各個纖維位置的導數和。第i個纖維位置處的截面兼容矩陣的導數是。如圖7所示，截面是由以下幾個集合參數定義的，總截面高d，腹板厚度，面板寬度，和面板厚度。從腹板的高度