有限元分析報告

1、精選優(yōu)質(zhì)文檔-傾情為你奉上西安市新城區(qū)某公司科研辦公樓結構設計有限元分析報告撰 寫 人： 王 平 班 級： 工程力學1203 學 號： 指導教師： 張衛(wèi)喜 2016年6月15日專心-專注-專業(yè)目錄1 工程概況工程名稱：西安市新城區(qū)某公司科研辦公樓；建筑所在地：西安市；建設規(guī)模：總建筑面積約4700m2，主體結構6層，無地下室。結構總高度22.5m，底層結構高度4.5m，其余層結構高度為3.6m，幾何模型圖如圖1所示；抗震設防烈度：抗震設防烈度為8度，設計基本地震加速度值0.2g,第一組。場地類別為類，特征周期為0.35s。周期折減系數(shù)為0.75。建筑設計使用年限：50年。結構重要性等級：二級。

2、圖1 框架幾何模型圖2 分析依據(jù)框架結構是由梁、板、柱以剛接相連接而成，構成承重體系的結構，即由梁、板、柱組成框架共同抵抗使用過程中出現(xiàn)的水平荷載和豎直荷載。本設計報告采用ANSYS有限元軟件分析。根據(jù)框架結構體系特點，本結構分析主要依據(jù)以下國家規(guī)范：1國家標準：建筑結構荷載規(guī)范（GB50009-2012).北京：中國建筑工業(yè)出版社.2012；2國家標準：建筑抗震設計規(guī)范（GB50011-2010).北京：中國建筑工業(yè)出版社.2010；3國家標準：混凝土結構設計規(guī)范（GB50010-2010).北京：中國建筑工業(yè)出版社.2010；4建筑、勘察等技術文件。3 荷載與計算工況3.1荷載簡化及荷載組

3、合1）恒載：含柱、梁、板自重；2）梁間活載以換算線荷載施加（含墻、門、窗等），如圖2所示；3）板面活荷載按建筑功能分類取值。圖2 框架結構梁間荷載示意圖3.2 邊界條件1）柱底固結；2）板四周與梁固結，如圖3所示。圖3 框架結構板邊界示意圖3.3 工況工況一：模態(tài)分析；工況一：分析結構可變荷載組合下的剛度與強度；（注:由于永久荷載組合與可變荷載組合相比較小，故在此不進行校核）工況二：對結構在水平地震力和恒載組合作用下進行分析。3.3.1模態(tài)分析給模型加約束之后，進行前6階模態(tài)分析。3.3.2 靜力分析在可變荷載效應控制的荷載作用下，對框架進行剛度、強度校核。4 有限元模型4.1 基本假定1）假

4、設混凝土為線彈性材料；2）鋼筋混凝土取整體分析，計鋼筋貢獻；3）節(jié)點假定為剛接，底層柱與基礎剛接，樓板與梁四邊剛接。4.2 力學模型1）主要受力構件為柱、梁、板；2）抗側(cè)力構件為柱；3）水平聯(lián)系構件為框架梁（考慮板的貢獻）；4）柱為壓彎構件，框架梁含彎曲、軸向變形、扭轉(zhuǎn)變形。4.3 主要物理參數(shù)取值表1 主要物理參數(shù)樓層混凝土等級彈性模量（N/m2）泊松比板厚（mm）柱（b×h）（mm）橫梁(b×h)縱梁(b×h)（mm）次梁（b×h）（mm）邊跨梁（b×h）（mm）中跨梁（b×h）（mm）1C303×10100.21006

5、00×600300×600300×400300×600250×50025500×5006120查閱相關資料得：鋼筋混凝土有限元模型根據(jù)鋼筋的處理方式主要分為三種，即分離式、組合式、整體式。分離式模型把鋼筋和混凝土作為不同的單元處理，即混凝土和鋼筋各自被劃分為足夠小的單元。考慮到鋼筋是一細長材料，通?？梢院雎云淇辜簦话沅摻羁勺鳛榫€單元來處理?；炷羷t作為實體單元來處理。分離式模型中在鋼筋和混凝土之間可以插入連接單元來模擬鋼筋與混凝土的粘結和滑移。當鋼筋和混凝土粘結較好，可以認為兩者無滑移時，采用組合式模型。組合式模型中最常用兩種方式第

6、一種是分層式，即在橫截面上分成許多混凝土層和若干鋼筋層，并對截面的應變做出某些假定，主要適用于桿件系統(tǒng)的分層組合式；另一種組合方式是采用等參數(shù)單元，使用帶鋼筋膜的單元來進行模擬，主要適用于二維或三維結構的帶膜組合式或帶鋼筋的組合式。整體式模型也是在考慮鋼筋與混凝土之間無滑移的情況下使用的。假定鋼筋彌散于整個單元中，并把單元視為連續(xù)均勻材料，由此求出的是綜合了混凝土與鋼筋兩種材料的單元剛度矩陣。在ANSYS中可以考慮的鋼筋混凝土模型一般是兩種，即分離式模型和整體式模型。分離式模型雖然在理論上可行，但在進行復雜結構的分析時，對計算機的性能要求較高，計算效率較低。因此本文空間結構計算模型采用的是整體

7、式有限單元模型來處理鋼筋混凝土結構。4.4單元選取4.4.1 梁、柱單元的選取梁、柱采用空間Timoshenko單元，擬合其彎曲、拉壓及扭轉(zhuǎn)受力狀態(tài)。Beam 188是一個二節(jié)點的三維線性梁, 適用于分析細長的梁, 其元素是基于Timoshenko 梁理論的，具有扭切變形效果。Beam 188 在每個節(jié)點上有6或7個自由度，（自由度）數(shù)目的變化是由KEYOPT(1)來控制的。當 KEYOPT(1) = 0時 (默認), 每節(jié)點有6個自由度。 分別是沿x,y,z的位移及繞其的轉(zhuǎn)動。 當 KEYOPT(1) = 1時,會添加第七個自由度 (翹曲量) ，Beam 188單元示意圖如圖4所示。圖4 T

8、imoshenko梁單元示意圖4.4.2 板單元的選取板采用彈性殼單元。Shell 181適用于薄到中等厚度的殼結構。該單元有四個節(jié)點,單元每個節(jié)點有六個自由度,分別為沿節(jié)點X,Y,Z方向的平動及繞節(jié)點X,Y,Z軸的轉(zhuǎn)動。退化的三角形選項用于網(wǎng)格生成的過渡單元。Shell 181單元具有應力剛化及大變形功能。該單元有強大的非線性功能,并有截面數(shù)據(jù)定義,分析,可視化等功能,還能定義復合材料多層殼；Shell 181單元示意圖如圖5所示。圖5 彈性殼單元示意圖4.5分網(wǎng)與有限元模型4.5.1 建模思路經(jīng)查閱ANSYS相關書籍得知Beam188 單元適合于分析從細長到中等粗短的梁結構，該單元基于鐵木

9、新科梁結構理論，并考慮了剪切變形的影響。有大量研究也表明該單元可以很好的模擬型鋼梁與柱，計算結果與實際比較接近，那是因為鋼材的材質(zhì)比較均勻。但框架混凝土梁與柱材質(zhì)并不均勻，混凝土內(nèi)部含各種鋼筋，故采用整體式模型，并把單元視為連續(xù)均勻材料，那么Beam188便可以很好的模擬框架混凝土梁與柱了。Shell 181適用于薄到中等厚度的殼結構.該單元有四個節(jié)點,單元每個節(jié)點有六個自由度,分別為沿節(jié)點X,Y,Z方向的平動及繞節(jié)點X,Y,Z軸的轉(zhuǎn)動。在此處剛好可以模擬板結構。綜上所述，得出以下建模思路：用Beam188、Shell181單元建立空間框架結構模型，與PKPM計算結果進行比較；4.5.2 分網(wǎng)

10、對梁和柱均設置為按份數(shù)進行劃分，份數(shù)取為20；對板進行劃分時，考慮到板的邊界與梁重合，故分網(wǎng)大小與梁相同，有限元模型圖如圖4所示。圖6 有限元模型圖5 靜力分析 5.1模態(tài)結果對框架結構進行模態(tài)分析，提取前六階振型圖，并統(tǒng)計前六階的頻率，周期和振型特點。表2 ANSYS分析結構動力參數(shù)階數(shù)周期（s)頻率（Hz）振幅特點10.8121.2309平動，橫向20.7561.3219平動，縱向30.7121.4046扭轉(zhuǎn)40.2603.8421平動，橫向50.2444.0930平動，縱向60.2294.3756扭轉(zhuǎn) 表3 PKPM計算結構動力參數(shù)階數(shù)周期（s)頻率（Hz）振幅特點10.9051.104

11、5平動，橫向20.8511.1758扭轉(zhuǎn)30.8421.1884平動，縱向40.2953.3944平動，橫向50.2773.6088平動，扭轉(zhuǎn)60.2763.6271平動，扭轉(zhuǎn)對表2及表3的數(shù)據(jù)進行分析可知，對于PKPM計算結果，結構在二階模態(tài)時發(fā)生扭轉(zhuǎn)變形；對于ANSYS計算結果，結構在三階模態(tài)時出現(xiàn)扭轉(zhuǎn)變形。ANSYS計算各階模態(tài)位移圖如圖813所示。MNMXXYZ 0.131E-03.261E-03.392E-03.522E-03.653E-03.783E-03.914E-03.NODAL SOLUTIONSTEP=1SUB =1FREQ=1.2309USUM (AVG)RSYS=0DM

12、X =.SMX =.圖7 第一階振形MNMXXYZ 0.928E-04.186E-03.278E-03.371E-03.464E-03.557E-03.649E-03.742E-03.835E-03NODAL SOLUTIONSTEP=1SUB =2FREQ=1.32186USUM (AVG)RSYS=0DMX =.835E-03SMX =.835E-03圖8 第二階振形MNMXXYZ 0.162E-03.324E-03.486E-03.648E-03.811E-03.973E-03.NODAL SOLUTIONSTEP=1SUB =3FREQ=1.40461USUM (AVG)RSYS=0D

13、MX =.SMX =.圖9 第三階振形MNMXXYZ 0.124E-03.248E-03.372E-03.496E-03.620E-03.744E-03.868E-03.992E-03.NODAL SOLUTIONSTEP=1SUB =4FREQ=3.84214USUM (AVG)RSYS=0DMX =.SMX =.圖10 第四階振形MNMXXYZ 0.910E-04.182E-03.273E-03.364E-03.455E-03.546E-03.637E-03.728E-03.819E-03NODAL SOLUTIONSTEP=1SUB =5FREQ=4.09304USUM (AVG)RSY

14、S=0DMX =.819E-03SMX =.819E-03圖11 第五階振形MNMXXYZ 0.156E-03.312E-03.467E-03.623E-03.779E-03.935E-03.00109.NODAL SOLUTIONSTEP=1SUB =6FREQ=4.37558USUM (AVG)RSYS=0DMX =.SMX =.圖12 第六階振形5.2靜力分析結果5.2.1 工況二MNMXXYZ 0.NODAL SOLUTIONSTEP=1SUB =1TIME=1USUM (AVG)RSYS=0DMX =.SMX =.圖13 可變荷載組合下框架結構變形圖由圖13可得，在可變荷載組合作用下

15、框架結構最大變形發(fā)生在屋面樓板，最大撓度為11.5mm。5.2.2 工況三MNMXXYZ 0.NODAL SOLUTIONSTEP=2SUB =1TIME=2UZ (AVG)RSYS=0DMX =.SMX =.圖14 整體框架位移圖MNMXXYZ 0.NODAL SOLUTIONSTEP=2SUB =1TIME=2UZ (AVG)RSYS=0DMX =.SMX =.圖15 單榀框架位移圖由上圖14和圖15可知，在只有地震力作用下，框架水平最大位移為36.634mm，PKPM電算結果為27.80mm。5.2.3層間位移角由ANSYS計算的各層位移及層間位移角與手算、PKPM電算的各層位移及層間位

16、移角如表4、表5和表6所示。表4 ANSYS計算各層位移、層間位移角層數(shù)絕對位移相對位移層間位移角636.6344.9891/721531.6455.2691/683426.3765.0591/711321.3175.7791/622215.5386.3791/56419.1596.8691/655表5 手算各層位移、層間位移角層數(shù)絕對位移相對位移層間位移角627.4242.4141/1491525.0103.6621/983421.3484.6721/771316.6765.4451/661211.2315.9811/60215.2505.2501/857表6 PKPM計算各層位移、層間位移

17、角層數(shù)絕對位移相對位移層間位移角627.802.501/1798525.303.301/1011422.004.651/744317.355.751/617211.606.451/55615.155.151/874由表4表6可知，ANSYS計算、PKPM電算以及手算結果的層間位移角均符合規(guī)范要求且相差不大。5.3 強度校核5.3.1 工況二可變荷載組合下框架結構的彎矩圖如圖17所示。XYZ -83010.6-41098.5813.62142725.784637.8LINE STRESSSTEP=1SUB =1TIME=1SMIS2 SMIS2MIN =-ELEM=840MAX =ELEM=40

18、418圖16 可變荷載組合下彎矩圖由圖16可知，在可變荷載組合下最大彎矩發(fā)生在梁跨中，為168.462kN·m。5.3.2 工況三水平地震力下框架結構的彎矩圖如圖18所示。XYZ -88636.4-1098066676.3LINE STRESSSTEP=1SUB =1TIME=1SMIS2 SMIS2MIN =-ELEM=MAX =ELEM=2841圖17 水平地震力下彎矩圖由圖可知，在水平地震荷載作用下，最大彎矩發(fā)生在一層柱底，為377.302kN·m。6 基于ANSYS、PKPM、手算的誤差分析6.1計算原理的不同6.1.1手算框架結構是多次超靜定結構，縱橫向的構件相互

19、制約、相互影響。如果要精確計算結構內(nèi)力，只能采用空間結構計算模型，而不能采用平面結構計算模型，所列的力法、位移法方程是很繁瑣的。為了簡化工作量，手算只計算一榀框架結構模型，而忽略其他方向，這樣計算是建立在結構規(guī)則的前提下的。建筑抗震設計規(guī)范里提出，建筑形體規(guī)則的參考指所有平面規(guī)則和豎向規(guī)則。平面規(guī)則有扭轉(zhuǎn)規(guī)則、凹凸規(guī)則；豎向規(guī)則則有側(cè)向剛度規(guī)則和樓層承載力突變不能太大。只有滿足這些條件，才可以簡化成平面結構。因此，可以在計算框架結構承受水平荷載作用下的內(nèi)力、位移時采用底部剪力分配法；在計算框架結構承受豎向荷載作用下的內(nèi)力、位移時采用彎矩分配法。6.1.2 ANSYS與PKPM兩者都是有限元軟件

20、，在有限元基本理論上都是一致的，表現(xiàn)如下：第一步：問題及求解域定義。根據(jù)實際問題近似確定求解域的物理性質(zhì)與幾何區(qū)域；第二步：求解域離散化。將求解域近似為具有不同有限大小和形狀且彼此相連的有限個單元組成的離散域，習慣性稱為有限元網(wǎng)格劃分；第三步：確定狀態(tài)變量及控制方法。一個具體的物理問題通?？梢杂靡唤M包含問題狀態(tài)變量邊界條件的微分方程表示，為適合有限元求解，通常將微分方程化為等價的泛函形式；第四步：單元推導。對單元構造一個適合的近似解，即推導有限單元的形式，其中包括選擇合理的單元坐標系，建立單元試函數(shù)，以某種方法給出單元各狀態(tài)變量的離散關系，從而形成單元矩陣；第五步：總裝求解。將單元總裝形成離散域的總矩陣方程（聯(lián)合方程組），反映對近似求解域的要求，即單元函數(shù)的連續(xù)性要滿足一定的連續(xù)條件；第六步：聯(lián)立方程組求解和結果解釋。有限元法最終導致聯(lián)立方程組。聯(lián)立方程組的求解可以用直接法、迭代法和隨機法。求解結果是單元節(jié)點處狀態(tài)變量的近似值。不同的是，兩者在單元庫、建模、后處理等子塊程序有較大差異。ANSYS擁有