本文選題：土 + 結(jié)構(gòu)相互作用��； 參考：《天津大學(xué)》2015年博士論文

【摘要】：近來年,對土—結(jié)構(gòu)相互作用問題的研究廣度和深度逐漸擴(kuò)大,各類復(fù)雜模型越來越龐大且精細(xì)化,導(dǎo)致在有限元分析中單元節(jié)點(diǎn)的自由度數(shù)可謂以十萬百萬計(jì),需高配置的計(jì)算設(shè)備和長時(shí)間的計(jì)算時(shí)間才可完成分析,當(dāng)需多個(gè)不同工況的比較時(shí)則更成倍地增加工作量,亟待一種可行的既滿足計(jì)算精度又節(jié)省計(jì)算成本的分析計(jì)算方法。動力子結(jié)構(gòu)法是將土—結(jié)構(gòu)相互作用體系分成若干子系統(tǒng),先分別求出各子系統(tǒng)的單體反應(yīng),然后再聯(lián)系各單體的反應(yīng),使其滿足相互作用的條件,從而得到整個(gè)體系的反應(yīng),大大減少計(jì)算量,是一種理想的計(jì)算方法,但該方法僅能解決線性問題。本文根據(jù)課題組新近提出的能解決局部非線性問題的線性—非線性混合約束模態(tài)綜合法的思想,在該方法的具體實(shí)施方面進(jìn)行了研究,充分利用現(xiàn)有商業(yè)軟件的優(yōu)勢,將該方法的實(shí)施與ANSYS軟件相結(jié)合,對實(shí)施過程中諸如結(jié)構(gòu)受力全過程中剛度矩陣的提取方法;彈塑性模型中等效彈性模量的計(jì)算;一致粘彈性邊界問題;以及非線性的接觸問題等進(jìn)行了深入的分析與研究,使混合約束模態(tài)綜合法能解決復(fù)雜狀況下的土—結(jié)構(gòu)相互作用體系,開辟了動力子結(jié)構(gòu)法的研究領(lǐng)域。本文主要的創(chuàng)新工作與成果有:1.混合約束模態(tài)綜合法在ANSYS的實(shí)施,首先要解決剛度矩陣提取中的對號入座問題,為此在深入分析ANSYS二次開發(fā)四種方法的基礎(chǔ)上,決定采取FORTRAN語言的二次開發(fā)工具UPFs進(jìn)行矩陣提取的開發(fā)工作,自編程序?qū)崿F(xiàn)了按照節(jié)點(diǎn)編號由小到大的順序進(jìn)行矩陣的提取。2.非線性子結(jié)構(gòu)在塑性階段中的剛度矩陣的提取工作是本文的重點(diǎn)難點(diǎn)之一。為了在非線性時(shí)程分析中生成某一時(shí)刻的彈塑性剛度矩陣,提出了將非線性塑性區(qū)域“等效”為線彈性區(qū)域,即在塑性階段上的小段Δt時(shí)間內(nèi)利用分段等效線性化手段,對結(jié)構(gòu)各個(gè)進(jìn)入塑性階段的單元進(jìn)行等效處理,并推導(dǎo)了等效彈性模量的公式,完成了對塑性單元的彈性模量和剪切模量等效工作;利用了ANSYS多點(diǎn)重啟動的方法,恢復(fù)結(jié)構(gòu)上一時(shí)刻的計(jì)算結(jié)果后能繼續(xù)進(jìn)行下一時(shí)刻分析;提出了ANSYS與MATLAB程序間的相互調(diào)用和控制的操作方法,從而解決了非線性子結(jié)構(gòu)在塑性階段中的剛度矩陣的有序進(jìn)行。通過算例模型進(jìn)行一般性分析,包括Pushover靜力非線性分析和簡諧波動力時(shí)程非線性分析求得響應(yīng)時(shí)程曲線,與ANSYS直接計(jì)算法所得到的響應(yīng)時(shí)程曲線比對吻合良好,驗(yàn)證了在彈塑性階段二次開發(fā)的程序文件所提取的矩陣的正確可行。3.建立了考慮局部非線性的地基土—剪力墻結(jié)構(gòu)相互作用的模型,然后根據(jù)結(jié)構(gòu)存在局部塑性區(qū)域的特點(diǎn),將整體體系中未進(jìn)入非線性階段的區(qū)域劃分為若干個(gè)線性子結(jié)構(gòu),而將進(jìn)入塑性階段的局部區(qū)域獨(dú)立劃分為若干個(gè)非線性子結(jié)構(gòu)。提出了混合約束模態(tài)綜合法的詳細(xì)分析流程并編寫了開發(fā)程序,著重說明了在ANSYS和MATLAB二者程序之間的實(shí)現(xiàn)過程。4.在土—結(jié)構(gòu)相互作用分析模型中,在地基土邊界引入一致粘彈性邊界條件來模擬半無限域的土體,推導(dǎo)了帶有一致粘彈性邊界的子結(jié)構(gòu)方程,論證了勢能判據(jù)截?cái)鄿?zhǔn)則同樣適用于帶有一致粘彈性邊界的線性子結(jié)構(gòu),證明了采用混合約束模態(tài)綜合法同樣能處理復(fù)雜的邊界問題。5.通過選擇合理的目標(biāo)面單元和接觸面單元、迭代計(jì)算接觸剛度和設(shè)置摩擦系數(shù),將高度非線性的接觸問題引入混合約束模態(tài)綜合法中。提出了箱型基礎(chǔ)與周圍土體接觸關(guān)系的處理方法,將接觸單元一分為二,目標(biāo)面單元劃入剛度相對較高的箱型基礎(chǔ)的線性子結(jié)構(gòu)內(nèi),接觸面單元劃入剛度相對較低的周圍土體的非線性子結(jié)構(gòu)內(nèi),推導(dǎo)出了箱型基礎(chǔ)與周圍土體具有接觸關(guān)系的線性與非線性模態(tài)綜合方程,拓寬了混合約束模態(tài)綜合法的應(yīng)用領(lǐng)域。6.通過算例計(jì)算對線性和非線性土體區(qū)域給出了推薦范圍,當(dāng)設(shè)置一致粘彈性邊界時(shí),非線性子結(jié)構(gòu)的區(qū)域?yàn)?倍的基礎(chǔ)寬度,兩側(cè)線性子結(jié)構(gòu)區(qū)域各為1.5倍的基礎(chǔ)寬度,使得模型總體自由度數(shù)比使用一般自由邊界子結(jié)構(gòu)土體單元自由度數(shù)減少很多,從而提高了計(jì)算效率。7.進(jìn)行了混合約束模態(tài)綜合法擴(kuò)展到三維土—結(jié)構(gòu)相互作用分析,重點(diǎn)研究了與二維問題的不同之處,重新推導(dǎo)了適合于三維問題的一致粘彈性邊界條件的子結(jié)構(gòu)方程;重新推導(dǎo)了箱型基礎(chǔ)與周圍土體具有接觸關(guān)系的線性與非線性模態(tài)綜合方程,以及三維模型等效彈性模量的公式,結(jié)果表明混合約束模態(tài)綜合法能解決三維土—結(jié)構(gòu)相互作用復(fù)雜問題。8.在混合約束模態(tài)綜合法的實(shí)施中,通過逐步提取的結(jié)構(gòu)體系全程階段等效彈塑性剛度矩陣,進(jìn)行特征值方程計(jì)算得到了結(jié)構(gòu)前n階主要頻率和振型在全時(shí)段變化情況,從而了解結(jié)構(gòu)剛度在彈塑性動靜力分析中的變化規(guī)律,可為利用結(jié)構(gòu)動力特性對結(jié)構(gòu)體系進(jìn)行深入分析研究提供支撐。
[Abstract]:In recent years, the breadth and depth of the soil structure interaction have been expanded and the complex models are becoming more and more large and fine. The free degree of unit nodes in the finite element analysis can be described as one hundred thousand million. It needs high configuration of computing equipment and long time calculation time to complete the analysis, when many different workers are needed. In comparison with the situation, the workload is more multiplied, and a feasible analytical method is needed to satisfy both the calculation precision and the cost of calculation. The dynamic substructure method divides the soil structure interaction system into several subsystems, first to find the single reaction of each subsystem, and then to contact the reaction of the various monomers, so that it satisfies the phase. The condition of interaction, thus obtaining the reaction of the whole system, greatly reducing the amount of calculation, is an ideal calculation method, but this method can only solve the linear problem. In this paper, the idea of the linear nonlinear mixed constraint mode synthesis method which can solve local nonlinear problems recently proposed by the project group is the concrete implementation of this method. In order to make full use of the advantages of the existing commercial software, the implementation of the method is combined with the ANSYS software, and the method of extracting the stiffness matrix in the whole process of structural stress, the calculation of the equivalent elastic modulus in the elastoplastic model, the uniform viscoelastic boundary boundary problem and the nonlinear contact problem are carried out in the implementation process. In the analysis and research, the mixed constrained modal synthesis method can solve the soil structure interaction system under complex conditions and open up the research field of the dynamic substructure method. The main innovation work and achievements of this paper are as follows: 1. the implementation of the hybrid constrained modal synthesis method in the ANSYS, first of all, we should solve the problem of the number of seats in the stiffness matrix extraction. On the basis of the in-depth analysis of the four methods of the two development of ANSYS, the two development tool UPFs of the FORTRAN language is adopted to develop the matrix extraction. The self compiled program realizes the extraction of the stiffness matrix of the.2. nonlinear substructure in the plastic phase according to the node number from small to large. In order to generate the elastoplastic stiffness matrix of a certain moment in the nonlinear time history analysis, it is proposed that the nonlinear plastic region is "equivalent" to the linear elastic region, that is, the piecewise equivalent linearization method is used in the small segment delta T time of the plastic stage, and the units that enter the plastic stage are carried out. The formula of the equivalent elastic modulus is derived, and the equivalent work of the elastic modulus and the shear modulus of the plastic unit is completed. Using the method of ANSYS multi point restart, the next time analysis can be carried out after the recovery of the calculation results at the time of the structure, and the operation of mutual call and control between the ANSYS and the MATLAB program is put forward. This method solves the order of the stiffness matrix of the nonlinear substructure in the plastic stage. Through the example model, the general analysis is carried out, including the Pushover static nonlinear analysis and the harmonic dynamic time history nonlinear analysis to obtain the response time history curve, which is in good agreement with the response time curve comparison obtained by the ANSYS direct calculation method. To verify the correctness and feasibility of the matrix extracted from the program documents developed at the two time of the elastoplastic stage, the model of the soil shear wall structure interaction with local nonlinearity is established by.3.. Then, according to the characteristics of the local plastic region in the structure, the region which has not entered the nonlinear stage in the whole system is divided into several lines. The partial region that enters the plastic stage is divided into several nonlinear substructures independently. The detailed analysis process of the mixed constrained modal synthesis method is proposed and the development program is written. The implementation process between the ANSYS and the MATLAB two programs is emphasized, and the.4. is in the soil structure interaction analysis model, and the foundation is in the foundation. The soil boundary is introduced into the uniform viscoelastic boundary condition to simulate the soil in the semi infinite domain. The substructure equation with the uniform viscoelastic boundary is derived. It is proved that the potential energy criterion truncation criterion is also applicable to the linear substructure with the uniform viscoelastic boundary. It is proved that the mixed constrained modal synthesis method can also deal with complex boundary questions. By selecting the reasonable target surface element and the contact surface element,.5., the contact stiffness and the friction coefficient are iteratively calculated, and the highly nonlinear contact problem is introduced into the mixed constrained modal synthesis method. The contact relationship between the box type foundation and the surrounding soil is proposed. The contact element is divided into two parts and the target surface unit is divided into the stiffness phase. In the linear substructure of the higher box foundation, the contact surface element is inserted into the nonlinear substructure of the surrounding soil with relatively low stiffness. The linear and nonlinear modal synthesis equations of the contact relationship between the box type foundation and the surrounding soil are derived, and the application field of the mixed constrained modal synthesis method,.6., is widened by the calculation example. The range of linear and nonlinear soil is recommended. When the uniform viscoelastic boundary is set, the region of the nonlinear substructure is 2 times the base width and the basic width of the linear substructure of the two sides is 1.5 times, making the total degree of freedom of the model less than the free degree of the soil unit using the general free boundary substructure. The hybrid constrained modal synthesis method is extended to the three-dimensional soil structure interaction analysis, and the difference between the.7. and the two dimensional problem is studied. The substructure equation of the uniform viscoelastic boundary condition suitable for three dimensional problems is rederived, and the contact relation between the box foundation and the surrounding soil is re derived. The linear and nonlinear modal synthesis equations of the system and the formula of the equivalent elastic modulus of the three-dimensional model show that the mixed constrained modal synthesis method can solve the complex problem of the three-dimensional soil structure interaction..8., in the implementation of the mixed constrained modal synthesis method, can be obtained by gradually extracting the equivalent elastoplastic stiffness matrix in the whole process stage of the structure. The eigenvalue equation is used to calculate the changes of the main frequencies and modes of the n order before the structure, so as to understand the change law of the structural stiffness in the elastoplastic dynamic static analysis, which can provide support for the in-depth analysis of the structural system by using the dynamic characteristics of the structure.
【學(xué)位授予單位】：天津大學(xué)
【學(xué)位級別】：博士
【學(xué)位授予年份】：2015
【分類號】：TU470

【相似文獻(xiàn)】

相關(guān)期刊論文 前10條

1 劉開祥;;我所模態(tài)綜合分析技術(shù)應(yīng)用研究取得重大成果[J];交通部上海船舶運(yùn)輸科學(xué)研究所學(xué)報(bào);1982年01期

2 殷學(xué)綱,董明,黃河;混合界面模態(tài)綜合法的模態(tài)修正方法[J];計(jì)算結(jié)構(gòu)力學(xué)及其應(yīng)用;1987年04期

3 邱向軍;剩余柔度模態(tài)與附著模態(tài)在模態(tài)綜合法中的等價(jià)性[J];航空學(xué)報(bào);1987年03期

4 關(guān)銀發(fā),毛清華;工程中動態(tài)特性的模態(tài)綜合解法和試驗(yàn)研究[J];機(jī)械設(shè)計(jì)與研究;1989年02期

5 伍旭強(qiáng),朱繼梅;實(shí)驗(yàn)?zāi)Ｐ偷哪B(tài)綜合技術(shù)[J];上海力學(xué);1992年03期

6 張德文;部件試驗(yàn)?zāi)B(tài)綜合的動柔度法[J];強(qiáng)度與環(huán)境;2004年03期

7 董興建;孟光;;基于虛擬單元的模態(tài)綜合法[J];振動工程學(xué)報(bào);2009年06期

8 李興泉;鄧兆祥;龐劍;李英強(qiáng);;模態(tài)綜合的界面自由度里茲變換方法[J];華南理工大學(xué)學(xué)報(bào)(自然科學(xué)版);2012年12期

9 朱德懋;結(jié)構(gòu)動力分析的模態(tài)綜合法[J];南京航空航天大學(xué)學(xué)報(bào);1978年04期

10 王文亮;杜作潤;陳康元;;模態(tài)綜合技術(shù)短評和一種新的改進(jìn)[J];航空學(xué)報(bào);1979年03期

相關(guān)會議論文 前10條

1 董興建;;基于虛擬單元的模態(tài)綜合法[A];第四屆全國動力學(xué)與控制青年學(xué)者研討會論文摘要集[C];2010年

2 方明霞;;包含軟子結(jié)構(gòu)的混合界面模態(tài)綜合法建立轎車非線性動力學(xué)模型[A];中國力學(xué)學(xué)會學(xué)術(shù)大會'2005論文摘要集（下）[C];2005年

3 鄭鋼鐵;劉明輝;;子結(jié)構(gòu)模態(tài)綜合的研究進(jìn)展與應(yīng)用[A];第十屆全國振動理論及應(yīng)用學(xué)術(shù)會議論文集（2011）下冊[C];2011年

4 路觀平;;連拱壩交叉模態(tài)綜合法的完備性[A];第八屆全國振動理論及應(yīng)用學(xué)術(shù)會議論文集摘要[C];2003年

5 王慶洋;曹登慶;鄒元杰;;基于子結(jié)構(gòu)模態(tài)綜合法的航天器復(fù)雜結(jié)構(gòu)固有特性研究[A];中國力學(xué)大會——2013論文摘要集[C];2013年

6 董興建;孟光;;基于剩余動柔度的自由界面模態(tài)綜合法[A];第八屆全國動力學(xué)與控制學(xué)術(shù)會議論文集[C];2008年

7 張惠僑;朱昌明;;用部件模態(tài)綜合法分析浮吊臂架的動態(tài)特性[A];中國機(jī)械工程學(xué)會物料搬運(yùn)學(xué)會第二屆年會論文集（三）--金屬結(jié)構(gòu)[C];1984年

8 陳國平;韋勇;;非比例阻尼結(jié)構(gòu)線性振動系統(tǒng)的子結(jié)構(gòu)模態(tài)綜合[A];第八屆全國振動理論及應(yīng)用學(xué)術(shù)會議論文集摘要[C];2003年

9 朱昊;郭軍;;基于動態(tài)子結(jié)構(gòu)的柔性附件伸展研究[A];北京力學(xué)會第18屆學(xué)術(shù)年會論文集[C];2012年

10 李洪亮;陳予恕;;模態(tài)綜合法在雙轉(zhuǎn)子系統(tǒng)中的應(yīng)用[A];第十三屆全國非線性振動暨第十屆全國非線性動力學(xué)和運(yùn)動穩(wěn)定性學(xué)術(shù)會議摘要集[C];2011年

相關(guān)重要報(bào)紙文章 前1條

1 ;航天CAE技術(shù)（九）[N];中國航天報(bào);2006年

相關(guān)博士學(xué)位論文 前1條

1 張海順;局部非線性土與結(jié)構(gòu)相互作用的混合約束模態(tài)綜合法研究[D];天津大學(xué);2015年

相關(guān)碩士學(xué)位論文 前9條

1 顏誠;基于模態(tài)綜合法2MW風(fēng)電增速箱柔體動力學(xué)分析及優(yōu)化設(shè)計(jì)[D];湘潭大學(xué);2015年

2 李英強(qiáng);基于模態(tài)綜合的汽車結(jié)構(gòu)優(yōu)化研究[D];重慶大學(xué);2012年

3 李光明;基于DMAP的模態(tài)綜合研究及實(shí)現(xiàn)[D];華中科技大學(xué);2012年

4 溫爭鳴;自由界面模態(tài)綜合法的研究與系統(tǒng)實(shí)現(xiàn)[D];華中科技大學(xué);2014年

5 安方;固定界面模態(tài)綜合法的理論分析[D];南京航空航天大學(xué);2008年

6 孫運(yùn)見;基于SiPESC平臺的模態(tài)綜合法研究及實(shí)現(xiàn)[D];大連理工大學(xué);2012年

7 袁昭旭;復(fù)雜航天器結(jié)構(gòu)動力學(xué)兩類精確建模方法研究及軟件開發(fā)[D];哈爾濱工業(yè)大學(xué);2014年

8 鄧俊;試驗(yàn)?zāi)B(tài)和解析模態(tài)相結(jié)合的復(fù)雜結(jié)構(gòu)模態(tài)綜合技術(shù)研究[D];西北工業(yè)大學(xué);2002年

9 王慶洋;有阻尼組合結(jié)構(gòu)時(shí)頻域參數(shù)辨識方法研究[D];哈爾濱工業(yè)大學(xué);2013年

，

本文編號：2101904