不知火舞的被虐|伊人天伊人天天综合网|博洛尼亚天气|任你懆这里只有精品4|久久美日韩精品久久|掌中之物漫画免费阅读观看|0丨d老妇

<p>高斯勒讓德積分是有限元中最常見的數(shù)值積分方法之一，在有限元中有著廣泛的應(yīng)用。實際上，關(guān)于該積分方法的書籍和公眾號文章也已經(jīng)較多，本文主要是基于現(xiàn)有的教程，基本上把該方法的具體理論和使用重復(fù)了一遍，另外基于常用的數(shù)值計算庫PETSc描述下在PETSc中如何使用高斯勒讓德積分（本文后續(xù)都將高斯勒讓德積分簡稱為高斯積分）。

ETABS中有兩種彈塑性時程分析方法，分別是非線性模態(tài)分析法（FNA法）和非線性直接積分法。其中FNA適合于帶有少量非線性連接單元的結(jié)構(gòu)，計算速度快是其主要優(yōu)點，在減隔震分析中多被采用；而非線性直接積分法適用范圍更廣，適用于除時間相關(guān)效應(yīng)外的所有非線性行為，適用性強是其主要優(yōu)點，在大震彈塑性時程分析當(dāng)中多被采用。本文主要介紹非線性直接積分法的相關(guān)設(shè)置與應(yīng)用。

5.77350269E-01 -5.77350269E-01 1.00000000E-02 npt = 3 coords = -5.77350269E-01 5.77350269E-01 1.00000000E-02 npt = 4 coords = 5.77350269E-01 5.77350269E-01 1.00000000E-02 因此Abaqus中平面應(yīng)力單元高斯積分點的順序為

假設(shè)是母單元的自然坐標(biāo)系，是由高斯積分點控制的坐標(biāo)系（術(shù)語可能不專業(yè)），假設(shè)高斯積分方案為。

06-通過多重積分法求解線接觸彈性變形Fortran和MATLAB程序，程序請見下文附件及百度網(wǎng)盤鏈接

07-通過多重積分法求解點接觸彈性變形Fortran和MATLAB程序，程序請見下文附件及百度網(wǎng)盤鏈接

，由于被積函數(shù)較為復(fù)雜，如果代數(shù)積分進(jìn)行計算要消耗大量的計算資源，因此有限元理論中引入數(shù)值積分，即我們熟悉的高斯積分和hammer積分，對于直角坐標(biāo)系我們采用高斯積分，對于面積或者體積坐標(biāo)系我們采用hammer數(shù)值積分，具體這兩類積分的原理大家可以參考相關(guān)數(shù)值積分教材即可，這里我們只利用hammer數(shù)值積分的結(jié)論。

(21)4、高斯積分公式（20）中的單元剛度矩陣通過數(shù)值積分求得，本案例中的四節(jié)點和八節(jié)點四邊形等參單元均采用2*2個積分點的高斯積分即可求得精確結(jié)果。高斯積分點的坐標(biāo)具體如圖所示。

而等式左端積分就表示在整個流場中單位時間內(nèi)噴泉噴出或吸收了多少水量。所以說，散度的意義是場的有源性。</strong></p><p>&nbsp;</p><p><strong style="background-color: rgb(255, 218, 81);">3.高斯公式的理解</strong></p><p>我們已經(jīng)知道了高斯公式等式左端就是體積流量的變化，等式右端是什么呢？

使用隱式梯形積分法，構(gòu)建變電站電磁瞬態(tài)仿真模型。使用順序高斯消去法，獲取變電站工況特征并對其進(jìn)行模擬，將此模擬結(jié)果作為電磁暫態(tài)仿真模型的計算環(huán)境。構(gòu)建仿真實驗環(huán)節(jié)，使用歷史數(shù)據(jù)對實驗對象不同工況下的過電壓走向進(jìn)行分析，并將其作為對照數(shù)據(jù)。經(jīng)實驗結(jié)果證實，仿真結(jié)果與對照數(shù)據(jù)走向一致，說明此方法可有效避免仿真結(jié)果失真問題。

因此積分點上的應(yīng)變和應(yīng)力相對于節(jié)點上是更加準(zhǔn)確的。 高斯積分點探測法 在面對面的接觸計算時，相比節(jié)點探測法來說，使用高斯積分點探測方法可以得到更加準(zhǔn)確的計算結(jié)果。在面對面的接觸中，如使用罰函數(shù)接觸算法或增廣拉格朗日算法時，Workbench默認(rèn)的探測方法是高斯積分點法。

如果只有一個分布，它們就會用最大似然法來估計。但是由于有 K 個這樣的集群，并且概率密度被定義為所有這些 K 分布的密度的線性函數(shù)，即其中 πk是 k的混合系數(shù)th分配。為了通過最大對數(shù)似然法估計參數(shù)，請計算 p(X∣μ,Σ,π)。現(xiàn)在定義一個隨機(jī)變量 γk(X)，使得γk(X)=ρ(k∣X)。

高斯積分的原理簡單來說是：對于f(x)在a-b區(qū)間的積分，我們可以取a到b之間的n個點，總的積分值等于這些點處的函數(shù)值的加權(quán)和。圖中的積分點坐標(biāo)是等參坐標(biāo)。等參變換可以理解為一種歸一化，等參變換不影響本文的討論，為理解上的直觀起見，本文仍采取自然坐標(biāo)xy描述。以一個平面應(yīng)力問題的四節(jié)點矩形單元為例。單元的坐標(biāo)系建立在中心。

——勒讓德積分從上圖我們可以看出單元積分點的分布特別符合高斯——勒讓德積分域（-1，1），于是采用高斯——勒讓德積分，相應(yīng)的積分點，權(quán)系數(shù)數(shù)值分析教材中肯定會涉及，針對兩點高斯積分，積分點為 ，權(quán)系數(shù)為1。

</strong></p><h2>2.5 高斯積分</h2><p>上述積分的計算需要使用數(shù)值積分，一維高斯積分點的坐標(biāo)以及權(quán)系數(shù)如下：</p><p><img src="https://img.jishulink.com/msimage/202403/e2b09e5be1a98ebfad3d72b9372ef67f.png"></p><p>采用2x2x2的高斯積分點。

對于這樣一種線性單元，在構(gòu)造剛度矩陣的時候，需要進(jìn)行下式所示的積分。（四節(jié)點矩形單元應(yīng)該是8×8)其中B矩陣是單元形函數(shù)對空間坐標(biāo)的相關(guān)偏導(dǎo)，D矩陣是本構(gòu)矩陣。該積分中的被積矩陣（8×8）的每一個元素都是一個三元函數(shù)，其針對單元域的積分值成為一個剛度系數(shù)。如上單元在高斯積分方案下的減縮積分就是取被積函數(shù)在積分域中心點的函數(shù)值乘以2（曾攀04P178），實際上就是梯形積分公式。