大跨度鋼桁架梁橋的地震響應(yīng)分析

1、大跨度橋梁結(jié)構(gòu)理論大跨度鋼桁架梁橋的地震響應(yīng)分析（南京長江大橋） 目錄1.緒論31.1課題研究背景31.2橋梁震害原因31.3橋梁抗震響應(yīng)主要影響因素42.大跨度橋梁抗震分析與地震波的輸入42.1動態(tài)時程分析法（直接動力分析理論）42.1.1時程分析法計(jì)算原理52.1.2時程分析法優(yōu)缺點(diǎn)52.2地震波的特性、輸入與選取52.2.1地震波的幅值62.2.2地震波的選取62.2.3地震波的輸入63.大跨度鋼桁架橋梁的數(shù)值分析63.1工程概況63.2該橋的基本設(shè)計(jì)參數(shù)73.3約束與連接的簡化處理73.4南京長江大橋南京聯(lián)主桁ANSYS建模73.4.1模型模態(tài)分析求解83.4.2模型模態(tài)計(jì)算結(jié)果103

2、.4.3模型模態(tài)結(jié)果分析104.大跨度鋼桁架橋梁地震響應(yīng)分析114.1地震波的選擇114.2時程分析134.3橫向地震波單向輸入下節(jié)點(diǎn)位移響應(yīng)分析141.緒論橋梁是鐵路工程，抗震防災(zāi)，危機(jī)管理系統(tǒng)的重要組成部分，在抗震救災(zāi)中有很重要的位置。為了減少地震災(zāi)害形成的破壞，不僅要對橋梁做好抗震加固措施，也要在橋梁設(shè)計(jì)上采取相應(yīng)措施以滿足抗震要求，提高橋梁的抗震能力，因此對響應(yīng)的橋梁地震響應(yīng)分析是必要的。 本文通過查閱相關(guān)文獻(xiàn)，對橋梁抗震的現(xiàn)狀，問題等作了分析總結(jié)，并采用大型通用軟件ANSYS對南京長江大橋中的一聯(lián)鋼桁架橋進(jìn)行了地震響應(yīng)分析，得出大跨度鋼桁架橋的動力特性和該橋體系的一些抗震性能，為以后

3、進(jìn)一步研究打下基礎(chǔ)。具體包括以下幾個方面：（1）通過查閱資料，闡述了對橋梁進(jìn)行抗震性能分析的現(xiàn)狀，存在的問題，以及對橋梁的減輕震害的一些設(shè)計(jì)思想，并總結(jié)了有關(guān)抗震性能分析的相關(guān)理論，為后期進(jìn)行橋梁結(jié)構(gòu)的抗震分析作基礎(chǔ)。（2）采用有限元軟件ANSYS對南京長江大橋其中一聯(lián)鋼桁架橋進(jìn)行模態(tài)分析，得到了橋梁結(jié)構(gòu)的自振頻率和振型圖，分析研究得到了其有關(guān)振動特性。（3）采用有限元軟件ANSYS對南京長江大橋其中一聯(lián)鋼桁架橋進(jìn)行地震波的瞬態(tài)分析，研究其分別在三種不同作用下該橋的抗震性能。通過分析研究，得到了橋梁跨中控制截面上下節(jié)點(diǎn)的位移響應(yīng)曲線，以及位移峰值，對其進(jìn)行對比分析，得出該橋體系的一些抗震性能為

4、后期抗震設(shè)計(jì)提供一些參考。1.1課題研究背景地震是一種不可控制的自然災(zāi)害，常常會造成嚴(yán)重的人員傷亡及建筑物破壞，給人類社會帶來嚴(yán)重災(zāi)難，造成巨大的經(jīng)濟(jì)損失。近年來已經(jīng)進(jìn)入多地震活躍期，地震活動非常頻繁。每年全世界有500百次地震的可能性，地震造成了嚴(yán)重的危害。為我們大家所熟悉的1976年，唐山大地震的震級為7.6級，轉(zhuǎn)眼瞬間，人數(shù)過百萬的重工業(yè)城市毀掉，給人們的安危以及國家?guī)碇卮髣?chuàng)傷，傷亡40余萬人，直接經(jīng)濟(jì)損失高達(dá)54億元。2008年，汶川大地震震級7.8級，對基礎(chǔ)設(shè)施造成嚴(yán)重?fù)p失，橋梁損傷達(dá)2900多座，人員傷亡和經(jīng)濟(jì)損失巨大。因此必須做好防震減災(zāi)工作。本文通過對一個大跨度鋼桁架梁橋的有

5、限元模擬分析，運(yùn)用ANSYS分析軟件對其合理建模，并對其模態(tài)進(jìn)行分析，得到自振頻率和振型，從而對其進(jìn)行三個不同激勵方向的地震波的瞬態(tài)分析，通過對比分析，求的該橋在不同地震響應(yīng)下，橋梁結(jié)構(gòu)桿件的位移變形，得到其抗震性能。1.2橋梁震害原因地震是不可控制和準(zhǔn)確預(yù)測的，所以從橋梁的設(shè)計(jì)開始時就設(shè)置減輕震害的一些措施很重要。通過對橋梁震害的研究和分析，得知其發(fā)生原因，由此建立并選擇合理的抗震設(shè)計(jì)方法，進(jìn)而采用相應(yīng)的加固方式。通過查閱相關(guān)地震災(zāi)害記錄，我們看到的橋梁損傷的原因主要有以下四類：（1）在對橋梁選址時，所選的場地對抗震不利，地震時候引起了地基失效或者地基變形。（2）橋梁抗震設(shè)防低于地震的強(qiáng)度。

6、（3）存在著一些人為因素，對工程進(jìn)行設(shè)計(jì)或者施工時候出現(xiàn)問題，應(yīng)該加大力度以避免此類事件的發(fā)生。（4）橋梁結(jié)構(gòu)本身抗震能力不足。通過對橋梁地震災(zāi)害的分析，我們發(fā)現(xiàn)，科學(xué)合理的抗震結(jié)構(gòu)體系，支撐結(jié)構(gòu)的合理優(yōu)化以及減隔震的一些方法能降低地震對工程造成的危害。因此我們要擇優(yōu)選取結(jié)構(gòu)體系，對關(guān)鍵部位進(jìn)行優(yōu)化處理，加強(qiáng)減隔震的加固方法，這對橋梁抗震有很大的意義。1.3橋梁抗震響應(yīng)主要影響因素對橋梁結(jié)構(gòu)進(jìn)行地震反應(yīng)分析，屬抗震動力學(xué)問題范疇。一般動力學(xué)問題都需要通過三個要素來說明，即輸入激勵，系統(tǒng)以及輸出反應(yīng)。兩個已知的元素是和未知的第三元素，是求解該類問題的過程。因此，對橋梁結(jié)構(gòu)的地震反應(yīng)分析，是把橋梁

7、結(jié)構(gòu)作為一個系統(tǒng)的。一般會采用有限元方法，將結(jié)構(gòu)分為許多離散的單元問題，在各個節(jié)點(diǎn)處連接起來成為一個集合整體；而把地震地面運(yùn)動看成是對系統(tǒng)的輸入；地震反應(yīng)作為一個系統(tǒng)的輸出。這樣，地震反應(yīng)分析就是已知地震輸入和結(jié)構(gòu)系統(tǒng)求地震反應(yīng)的問題。在進(jìn)行地震反應(yīng)分析時，準(zhǔn)確無誤地預(yù)測出地震對橋梁結(jié)構(gòu)的影響是進(jìn)行抗震設(shè)計(jì)的基礎(chǔ)。對大跨度橋梁的地震響應(yīng)能產(chǎn)生影響的因素有很多，主要包括多點(diǎn)激勵，行波效應(yīng)，多種非線性因素以及阻尼的問題和樁基-土-結(jié)構(gòu)的相互作用等等。2.大跨度橋梁抗震分析與地震波的輸入2.1動態(tài)時程分析法（直接動力分析理論）20世紀(jì)60年代抗震分析方法又出來一種新的方法，動態(tài)時程分析法。它被用于對

8、工程以及超高層建筑的抗震性能研究。發(fā)展至20世紀(jì)80年代時，國家已經(jīng)開始普遍采用此法，并將它作為抗震設(shè)計(jì)規(guī)范的分析法之一。時程分析方法又被叫做直接動力分析方法，采用逐步積分的方法對工程結(jié)構(gòu)進(jìn)行求解的運(yùn)動微分方程。采用該方法，工程結(jié)構(gòu)桿件的內(nèi)力，位移的時程關(guān)系可以通過各個粒子隨著時間變化的速度和加速度及位移的動態(tài)響應(yīng)。而且通過時程分析方法，結(jié)構(gòu)地震響應(yīng)能夠被準(zhǔn)確真實(shí)地表達(dá)，補(bǔ)充了反應(yīng)譜法的缺陷，從而可以對反應(yīng)譜法進(jìn)行彌補(bǔ)和對比。時程分析方法在理論上與實(shí)際應(yīng)用當(dāng)中分別被稱作步步積分法與動態(tài)設(shè)計(jì)。在工程結(jié)構(gòu)抗震分析，可以用來解決運(yùn)動輸入地面加速度記錄結(jié)構(gòu)方程的積分表達(dá)式。地震作用輸進(jìn)去時程分析方法以

9、及結(jié)構(gòu)所在場地的地震波，自始至終，步步積分，直至地震作用終止。2.1.1時程分析法計(jì)算原理當(dāng)橋梁結(jié)構(gòu)的各個支撐處受到地面震動作用時，運(yùn)動方程為：MsMsbMsbTMbXsXb+CsCsbCsbTCbXsXb+KsKsbKsbTKbXsXb=0Pb其中Xb各支座的地面強(qiáng)迫位移 Xs結(jié)構(gòu)中所有非支座節(jié)點(diǎn)的位移 Pb地面作用于各支座的力，即支座的支反力M、C、K節(jié)點(diǎn)的質(zhì)量矩陣，阻尼矩陣和剛度矩陣當(dāng)采用集中質(zhì)量的離散化模型時，Msb=0。 絕對位移 Xx可分解為擬靜位移 Ys和動態(tài)相對位移 Yr ，忽略 Xb所產(chǎn)生的阻尼力，則運(yùn)動方程可簡化為：MsYr+CsYr+KsYr=MsKs-1KsbXb節(jié)點(diǎn)的

10、動力位移和擬精力位移的和為節(jié)點(diǎn)的總位移。每一步單元內(nèi)力運(yùn)用上式運(yùn)動方程，采用逐步積分法，節(jié)點(diǎn)的動力位移與加速度增量即可被計(jì)算出。2.1.2時程分析法優(yōu)缺點(diǎn)進(jìn)行時程分析可以得到數(shù)值上較為精確的分析結(jié)果，但是存在著一些參數(shù)難以確定的問題，因而本質(zhì)仍然比較模糊。其他問題如：輸入地震動；簡化結(jié)構(gòu)分析模型是否與實(shí)際相符；結(jié)構(gòu)-基礎(chǔ)-土相互作用問題；結(jié)構(gòu)構(gòu)件的非線性動力特性和屈服后的行為；數(shù)值積分的精度及穩(wěn)定性等都有待于解決。時程分析方法的缺點(diǎn)是很顯然的，計(jì)算量過于龐大，而且由于是用確定性的時間歷程來模擬尚未發(fā)生的地震，理論上應(yīng)取許多條地面運(yùn)動加速度曲線作為樣本分別進(jìn)行計(jì)算后進(jìn)行統(tǒng)計(jì)分析比較合理。但受到計(jì)

11、算效率低的制約，目前在工程上只能用很少幾條樣本曲線進(jìn)行結(jié)構(gòu)分析統(tǒng)計(jì)。往往難于保證得到可靠的統(tǒng)計(jì)量。2.2地震波的特性、輸入與選取為了得到橋梁在地震中每時每刻的地震反應(yīng)，需要對橋梁進(jìn)行動力時程分析，而動力時程分析是以輸入合理的地震波為基礎(chǔ)的，因此，選擇合理的地震波是對橋梁進(jìn)行動力時程分析和橋梁抗震設(shè)計(jì)的基礎(chǔ)。但地震波的選取受多種因素的影響，包括該橋的動力特性，地震帶和場址類型等。而且不同的橋在同一地點(diǎn)或在不同地點(diǎn)的同一座橋，地震波的選擇也是不一樣的。2.2.1地震波的幅值地震動峰值加速度（PGA）的大小可以反映地震波幅值。地震動幅值有水平方向和垂直方向的，地震動 水平方向的幅值比較大，因此水平地

12、震的剪切運(yùn)動對工程結(jié)構(gòu)的危害比較大?？偧铀俣确逯蛋l(fā)生在高頻率，高頻率的振動，并不損壞橋，可以對橋梁造成危害的主要是從基礎(chǔ)的運(yùn)動速度，尤其是基礎(chǔ)位移?；诖耍岢隽擞行Х逯档母拍睢Ｓ纱伺缮拿~有有效峰值加速度、有效峰值速度、持續(xù)加速度、持續(xù)速度、均方根加速度、有效反應(yīng)譜加速度等。2.2.2地震波的選取地震波幅值的選取原則是選擇的典型強(qiáng)震地震波記錄與擬建場地的抗震設(shè)防烈度，對應(yīng)的加速度峰值大小相同，如果選擇的強(qiáng)震地震波的加速度峰值與抗震設(shè)防對應(yīng)的加速度峰值不同，則調(diào)整強(qiáng)震記錄到加速度峰值相當(dāng)。在選擇強(qiáng)震記錄時，最大峰值加速度應(yīng)符合橋梁所在地區(qū)的烈度要求外，場地條件也應(yīng)盡量接近，也就是地震加速度時

13、程的主要周期應(yīng)盡量接近橋址區(qū)場地的卓越周期。2.2.3地震波的輸入地震輸入可采用一致激勵和多點(diǎn)激勵模式，使用于各類橋梁線形和非線性的地震反應(yīng)分析。當(dāng)大跨度橋梁結(jié)構(gòu)的跨度與地震波的波長處于相同數(shù)量級時，應(yīng)考慮地震輸入時間和空間的變異性，為了能夠比較合理地預(yù)測未來地震作用下橋梁結(jié)構(gòu)的地震反應(yīng)。目前通常采用兩種方法去選擇輸入地震波，第一種方法稱為類比地震波方法，第二種方法稱為人工地震波方法。如果橋梁結(jié)構(gòu)是在指定橋址的情況下進(jìn)行地震分析的，則對地震波進(jìn)行輸入應(yīng)選擇第一種方法，即類比地震波方法，但必須先對橋梁地震站點(diǎn)討論和基本的安全評價工作。此時便可以選擇相似的地震記錄作為地震波進(jìn)行輸入，這就要求所選地

14、震波與所建橋梁的地址環(huán)境相類似，并且具有類似的地震級別的地震動記載。一般情況下，對橋梁進(jìn)行地震分析時需采用若干條相似地震波，依據(jù)分別進(jìn)行的時程分析的計(jì)算結(jié)果，從而更科學(xué)合理對該橋結(jié)構(gòu)進(jìn)行綜合抗震性能評價 。3.大跨度鋼桁架橋梁的數(shù)值分析3.1工程概況南京長江大橋是我國自行設(shè)計(jì)建造的首座特大鐵路，公路兩用橋梁，由正橋和引橋兩部分組成。其中正橋分上下兩層，上層為公路，四車道，公路行車道寬15m，兩側(cè)各設(shè)凈寬2.5m的人行道，總寬為19.5m。下層為雙線鐵路，兩線路中心距為4m.正橋上部結(jié)構(gòu)由10孔鋼梁組成，浦口岸第一孔為跨徑128m的簡支鋼絎梁，其余9孔為三聯(lián)三等跨160m的連續(xù)鋼桁梁。主桁為平弦

15、菱形桁，節(jié)間距8m，桁高16m,下加勁桁弦高14m；兩主桁中心距為14m.3.2該橋的基本設(shè)計(jì)參數(shù)該橋主桁桿件有31種不同類型的截面，眾多類型的縱橫向連接系以及鐵路，公路兩種不同橋面系等，結(jié)構(gòu)體系十分復(fù)雜。為了能夠較準(zhǔn)確地分析，根據(jù)大橋的結(jié)構(gòu)特點(diǎn)，選取了南京聯(lián)主桁（160m三跨連續(xù)鋼桁梁）為研究對象進(jìn)行有限元模型建模，同時根據(jù)相關(guān)技術(shù)資料確定了大橋各主要構(gòu)件的截面性質(zhì)和材料性質(zhì)，如下表所示：南京長江大橋鋼桁梁桿件截面特性截面性質(zhì)截面積（m2）豎彎慣性矩（m4）橫彎慣性矩（m4）扭轉(zhuǎn)慣性矩（m4）數(shù)值0.343310-20.199840.128610-50.125870.128610-50.01

16、78390.266610-30.13685南京長江大橋材料特性構(gòu)件名稱所有材料彈性模量（pa）泊松比密度（kgm3）連續(xù)鋼桁梁I級鋼2.110110.287.85103公路橋面輕質(zhì)混凝土-2.01033.3約束與連接的簡化處理在建模時，三跨連續(xù)鋼桁梁橋的約束條件通過限制模型中的相應(yīng)節(jié)點(diǎn)自由度來實(shí)現(xiàn)。具體做法是，按照設(shè)計(jì)圖紙，對四個支撐處都施加豎向，繞豎向和橫向的轉(zhuǎn)動約束，對中跨的兩個支撐施加縱向約束，四個支撐處在同一側(cè)分別選取一個支撐點(diǎn)進(jìn)行橫向約束。3.4南京長江大橋南京聯(lián)主桁ANSYS建模建模假設(shè)基于以下假設(shè)的基礎(chǔ)上形成了有限元模型：（1）假設(shè)橋面板和縱梁的形心在同一水平面上；（2）假設(shè)各節(jié)

17、點(diǎn)滿足理想連接條件；（3）忽略橋面板的垂度變化；（4）建模坐標(biāo)定義為，順橋向?yàn)?X方向，橫橋向?yàn)閅方向，豎向?yàn)閆向。采用ANSYS軟件建立該連續(xù)鋼桁梁橋有限元模型時，單元選取如下：主桁桿件，鐵路縱梁，橫梁，上下平橫向連接系等采用三維梁單元BEAM4模擬，公路縱梁，公路橋面混凝土以及鐵路混凝土軌枕，無縫鋼軌等對結(jié)構(gòu)剛度貢獻(xiàn)較小的部分簡化為質(zhì)量單元MASS21，質(zhì)量單元加在上，下平聯(lián)的相關(guān)節(jié)點(diǎn)上。所建立的有限元模型共有7960個單元，5227個節(jié)點(diǎn)，如圖所示。3.4.1模型模態(tài)分析求解該大跨度鋼桁梁橋的前十階振型圖如下： 一階振型圖 二階振型圖 三階振型圖 四階振型圖五階振型圖 六階振型圖七階振型

18、圖 八階振型圖九階振型圖 十階振型圖3.4.2模型模態(tài)計(jì)算結(jié)果通過ANSYS軟件對該大跨度鋼桁梁橋進(jìn)行模態(tài)分析，由于低階振型對結(jié)構(gòu)起控制作用，所以只提取前十階振型進(jìn)行分析。其自振頻率和周期如下表所示：階次自振頻率f/Hz自振周期1 0.15289 6.54042 0.17727 5.64113 0.22345 4.47534 0.36890 2.71085 0.56478 1.77066 0.62091 1.61057 0.69155 1.44608 1.2485 0.80099 1.3697 0.730110 1.4913 0.6706 模態(tài)分析是動力學(xué)分析中很重要的基礎(chǔ)分析，是其他動力學(xué)分

19、析的前提，通常用于研究結(jié)構(gòu)的振動性質(zhì)，從而計(jì)算結(jié)構(gòu)的固有頻率及主振型。這些參數(shù)是重要的結(jié)構(gòu)參數(shù)，并成為其他類型動力學(xué)分析的基礎(chǔ)。3.4.3模型模態(tài)結(jié)果分析 通過以上分析發(fā)現(xiàn)該大跨度鋼桁梁橋以橫橋向振動為主，故對橫橋向的地震響應(yīng)最為敏感。4.大跨度鋼桁架橋梁地震響應(yīng)分析4.1地震波的選擇橋址區(qū)大地構(gòu)造單元地處揚(yáng)子準(zhǔn)地臺下寧蕪斷裂帶中的寧蕪斷陷西北緣的江浦坳陷。江浦坳陷是寧蕪斷陷盆地內(nèi)的次級構(gòu)造，形成晚，經(jīng)歷的構(gòu)造運(yùn)動少，地質(zhì)構(gòu)造簡單，斷裂不甚發(fā)育。根據(jù)本橋工程場區(qū)地震危險(xiǎn)性評價報(bào)告50年超越概率10%的地震基本烈度為度。場址地震危險(xiǎn)性分析結(jié)果50年超越概率10%3%2%地震烈度（度）6.57.0

20、7.1基巖水平峰值加速度0.0890.1270.140我國規(guī)范對公路橋梁的抗震有關(guān)規(guī)定（1）公路橋梁抗震設(shè)計(jì)規(guī)范規(guī)定：“按規(guī)范進(jìn)行抗震設(shè)計(jì)的橋梁，其設(shè)防目標(biāo)是，當(dāng)遭受橋梁設(shè)計(jì)基準(zhǔn)期內(nèi)發(fā)生概率較高的多遇地震影響時，一般不受損壞或不需修理可繼續(xù)使用，當(dāng)遭受橋梁設(shè)計(jì)基準(zhǔn)期內(nèi)發(fā)生概率較低的罕遇地震影響時，應(yīng)保證不致倒塌或產(chǎn)生嚴(yán)重結(jié)構(gòu)損傷，經(jīng)加固修復(fù)后仍可繼續(xù)使用?！保?）抗震設(shè)防烈度為6度及以上地區(qū)的公路橋梁，必須進(jìn)行抗震設(shè)計(jì)。（3）公路橋梁抗震設(shè)防烈度和設(shè)計(jì)基本地震加速度取值的對應(yīng)關(guān)系，應(yīng)符合規(guī)定。設(shè)計(jì)基本地震加速度值為0.15g和0.30g地區(qū)內(nèi)的公路橋梁，除規(guī)范另有規(guī)定外，應(yīng)分別按抗震設(shè)防烈度7

21、度和8度的要求進(jìn)行抗震設(shè)計(jì)?？拐鹪O(shè)防烈度和設(shè)計(jì)基本地震加速度a的對應(yīng)關(guān)系抗震設(shè)防烈度6789設(shè)計(jì)基本地震加速度值0.05g0.10(0.15)g0.20(0.3)g0.40g選擇了1940，EI Centro site，270Deg地震波其peak=0.3569g,Duration=53.72sec.時程曲線如下所示：地震加速度數(shù)據(jù)時間步/s 加速度/(m/s*2) 時間步/s 加速度/(m/s*2)2.00E-020.05066901 2.80E-01 0.13991464.00E-020.041274758 3.00E-01 0.108932886.00E-020.031780518 3.

22、20E-01 0.0769529328.00E-020.022286376 3.40E-01 0.0454721961.00E-010.016489872 3.60E-01 0.013991461.20E-010.04377317 3.80E-01-0.0175892361.40E-010.071056468 4.00E-01-0.0241852241.60E-010.069957202 4.20E-01 0.043873131.80E-010.045672018 4.40E-01 0.111931682.00E-010.021386932 4.60E-01 0.169895742.20E-0

23、1 -0.002938197 4.80E-01 0.0419742822.40E-010.047371044 5.00E-01-0.0859473722.60E-010.097740202 5.20E-01-0.21386932地震加速度數(shù)據(jù)時間步/s 加速度/(m/s*2) 時間步/s 加速度/(m/s*2)5.40E-01 -0.10193764 1.36E+00 -0.53367375.60E-01 0.082149676 1.38E+00-0.73854765.80E-01 0.0982401 1.40E+00-0.99938446.00E-01 0.11892692 1.42E+00

24、 -1.26922746.20E-01 0.18188898 1.44E+00 -1.52906466.40E-01 0.24485006 1.46E+00 -1.6190096.60E-01 0.28682444 1.48E+00 -1.30920166.80E-01 0.1409142 1.50E+00 -1.00938047.00E-01 0.10293724 1.52E+00 -0.90144817.20E-01 0.0995386 1.54E+00 -0.968407587.40E-01 0.096440918 1.56E+00 -0.939425067.60E-01 0.03377

25、9326 1.58E+00 -0.910442547.80E-01 -0.077652456 1.60E+00 -1.03936848.00E-01 -0.18888422 1.62E+00 -1.18926928.20E-01 -0.31580696 1.64E+00 -0.788516828.40E-01 -0.33179664 1.66E+00 0.920436588.60E-01 -0.37776844 1.68E+00 2.34856028.80E-01 -0.49569674 1.70E+00 2.15867549.00E-01 -0.6026314 1.72E+00 1.9787

26、8669.20E-01 -0.51768304 1.74E+00 1.79889789.40E-01 -0.43373428 1.76E+00 1.6190099.60E-01 -0.34978552 1.78E+00 1.65898329.80E-01 -0.26483814 1.80E+00 1.69895741.00E+00 -0.36077916 1.82E+00 1.72894541.02E+00 -0.46671422 1.84E+00 1.64898721.04E+00 -0.57364888 1.86E+00 1.35917181.06E+00 -0.68957798 1.88

27、E+00 1.05935061.08E+00 -0.80650668 1.90E+00 0.76353271.10E+00 -0.75953528 1.92E+00 0.429736861.12E+00 -0.68258176 1.94E+00 0.0259841121.14E+00 -0.5086886 1.96E+00 -0.377768441.16E+00 -0.41974282 1.98E+00 -0.780521981.18E+00 -0.24784788 2.00E+00 -1.17927321.20E+00 0.30581292 2.02E+00 -1.58903081.22E+00 0.86047332 2.04E+00 -1.98878261.24E+00 0.67358732 2.06E+00 - 2.