哈爾濱某大學(xué)結(jié)構(gòu)力學(xué)I一一

SubroutineElem_Stiff(······)

說(shuō)明Stiff=0.0!單元?jiǎng)偠惹辶鉙electCase(Type)Case(1)

平面桿系結(jié)構(gòu)單元Case(2)

空間桿系結(jié)構(gòu)單元CaseDefault

出錯(cuò)信息EndSelectEndSubroutineElem_Stiff3.4桿系結(jié)構(gòu)單元分析子程序3.4.1單元?jiǎng)偠瓤傮w設(shè)計(jì)3.4.2說(shuō)明部分設(shè)計(jì)Integer,Intent(in)::···

入口整型參數(shù)Real(8),Intent(in)::···

入口實(shí)型參數(shù)Real(8),Intent(out)::···

出口實(shí)型參數(shù)Real(8)::Work1,···Integer::i,j,k,···實(shí)型和整型工作變量3.4桿系結(jié)構(gòu)單元分析子程序3.4.3平面桿系結(jié)構(gòu)設(shè)計(jì)SelectCase(Plane)Case(1)

平面桁架元素賦值Case(2)

平面梁柱元素賦值Case(3)······CaseDefault

出錯(cuò)信息EndSelect3.4.4空間桿系結(jié)構(gòu)設(shè)計(jì)SelectCase(Space)Case(1)

空間桁架元素賦值Case(2)

空間梁柱元素賦值Case(3)

交叉梁元素賦值CaseDefault

出錯(cuò)信息EndSelect3.4桿系結(jié)構(gòu)單元分析子程序3.4.5有關(guān)單元等效結(jié)點(diǎn)荷載設(shè)計(jì)和進(jìn)一步的考慮1)單元等效結(jié)點(diǎn)荷載設(shè)計(jì)同仿單元?jiǎng)偠取?)從各類(lèi)單元?jiǎng)偠仍氐挠?jì)算，可看到要用到長(zhǎng)度、單元彈性特性、單元截面特性等數(shù)據(jù)。因此，要確定存放它們的數(shù)據(jù)結(jié)構(gòu)。要將它們作為出口。3)為計(jì)算單元等效結(jié)點(diǎn)荷載元素，首先要建立各種荷載情況等效荷載表達(dá)式，它們可由積分或載常數(shù)表得到。然后要解決荷載信息的存放結(jié)構(gòu)，也要將它們作為出口量。4)單元?jiǎng)偠染仃嚒⒌刃ЫY(jié)點(diǎn)荷載矩陣都應(yīng)先清零。4.1桿系結(jié)構(gòu)整體分析

首先就全剛結(jié)點(diǎn)平面剛架進(jìn)行討論，然后推廣。4.1.1總的思路在單元特性搞清后，將單元拼裝回去。在結(jié)點(diǎn)處位移自動(dòng)協(xié)調(diào)基礎(chǔ)上，如果全部結(jié)點(diǎn)平衡，則求得的結(jié)點(diǎn)位移將是實(shí)際結(jié)構(gòu)的解。因此，整體分析就是設(shè)法建立結(jié)點(diǎn)平衡方程。4.1.2坐標(biāo)轉(zhuǎn)換組成結(jié)構(gòu)的桿件可以各個(gè)方向，單元分析對(duì)局部坐標(biāo)，因此，必須將物理量轉(zhuǎn)為統(tǒng)一坐標(biāo)——整體坐標(biāo)。1)力的轉(zhuǎn)換關(guān)系4.1桿系結(jié)構(gòu)整體分析2)位移轉(zhuǎn)換關(guān)系3)轉(zhuǎn)換矩陣轉(zhuǎn)換矩陣是正交矩陣。4.1桿系結(jié)構(gòu)整體分析4)桿端力轉(zhuǎn)換5)桿端位移轉(zhuǎn)換6)剛度方程的轉(zhuǎn)換如果記稱(chēng)為整體單元?jiǎng)偠染仃噭t這就是整體坐標(biāo)下的單元?jiǎng)偠确匠?。本?jié)以后的討論認(rèn)為都是對(duì)整體坐標(biāo)的4.1桿系結(jié)構(gòu)整體分析4.1.3結(jié)點(diǎn)平衡方程的建立1)一簡(jiǎn)單例子（如圖）圖中有兩套編號(hào)，紅的是單元桿端編號(hào)，黑的是結(jié)構(gòu)整體編號(hào)。1-1)結(jié)點(diǎn)示意121221①②③圖中藍(lán)色的表示結(jié)點(diǎn)荷載（已知），紅色的表示桿端力（未知的），、分別1、2單元桿端力子矩陣。對(duì)1、4結(jié)點(diǎn)“荷載”含有未知反力。21-2)結(jié)點(diǎn)平衡由示意圖可見(jiàn)，結(jié)構(gòu)結(jié)點(diǎn)的平衡方程為4.1桿系結(jié)構(gòu)整體分析從例圖可見(jiàn)，其全部結(jié)點(diǎn)平衡方程為121221①②③若記24.1桿系結(jié)構(gòu)整體分析式中[I]、[0]分別為單位和零矩陣。若引入矩陣記號(hào)，則結(jié)點(diǎn)平衡方程可改寫(xiě)作這一結(jié)論雖然是由一個(gè)例子得到的，但是顯然對(duì)一切結(jié)構(gòu)都是成立的。問(wèn)題在于不同結(jié)構(gòu)，[A]矩陣是不同的。4.1桿系結(jié)構(gòu)整體分析4.1.4桿端位移用結(jié)點(diǎn)位移來(lái)表示121221①②③仍以簡(jiǎn)單例子來(lái)說(shuō)明若記由結(jié)點(diǎn)、桿端位移的協(xié)調(diào)條件，可得[]、[]的對(duì)應(yīng)關(guān)系為式中[A]T是前面力關(guān)系[A]的轉(zhuǎn)置，因此[A]T稱(chēng)為位移轉(zhuǎn)換矩陣。4.1桿系結(jié)構(gòu)整體分析4.1.5整體剛度方程——結(jié)點(diǎn)平衡121221①②③若記引入位移轉(zhuǎn)換關(guān)系，則這就是整體剛度方程，它的物理實(shí)質(zhì)是結(jié)點(diǎn)平衡。[K]稱(chēng)作結(jié)構(gòu)剛度矩陣（或整體剛度矩陣），[R]稱(chēng)作綜合等效結(jié)點(diǎn)荷載矩陣，它由兩部分組成。單元個(gè)數(shù)4.1桿系結(jié)構(gòu)整體分析4.1.6整體剛度矩陣的建立121221①②③若將[A]按單元分成圖示三個(gè)子矩陣則由此可見(jiàn)，整體剛度矩陣可由各單元整體剛度矩陣裝配累加得到。為說(shuō)明如何裝配，先將單元?jiǎng)偠染仃囘M(jìn)行分割整體結(jié)點(diǎn)碼則由矩陣乘法可證明，[A]i[k]i[A]iT的結(jié)果是，將剛度矩陣子矩陣按整體結(jié)點(diǎn)碼r、s送整體剛度矩陣相應(yīng)位置。這一裝配規(guī)則稱(chēng)為“對(duì)號(hào)入座”。4.1桿系結(jié)構(gòu)整體分析1)任意結(jié)構(gòu)情況上面結(jié)論是通過(guò)具體例子（全剛結(jié)點(diǎn)平面剛架）得到的，由虛位移原理或勢(shì)能原理進(jìn)行整體分析（見(jiàn)講義），可得任意結(jié)構(gòu)其結(jié)論同此例。2)結(jié)點(diǎn)位移編號(hào)如果按結(jié)點(diǎn)順序，對(duì)結(jié)點(diǎn)非零位移進(jìn)行依次編號(hào)，這一序號(hào)稱(chēng)作結(jié)點(diǎn)位移碼。為便于計(jì)算機(jī)處理并減少結(jié)構(gòu)剛度矩陣的階次，將零位移的號(hào)碼變?yōu)榱恪?duì)圖示三鉸剛架，當(dāng)僅用一種單元（梁柱自由是單元）時(shí)結(jié)點(diǎn)位移編號(hào)如圖所示。3)單元定位向量按單元局部結(jié)點(diǎn)碼順序，將結(jié)點(diǎn)位移碼排成的向量，稱(chēng)作單元的定位向量。①②③④①②③④4.1桿系結(jié)構(gòu)整體分析對(duì)圖示剛架各單元的定位向量為①(0,0,1,3,4,5)②(0,0,2,10,11,12)③(3,4,5,6,7,8)④(6,7,9,10,11,12)①(0,0,1,2,3)②(0,0,,6,7,8)③(1,2,3,4,5)④(4,5,6,7,8)4)按單元定位向量集裝剛度矩陣和綜合荷載前面說(shuō)明的是分塊子矩陣集裝,下面說(shuō)明如何按定為向量來(lái)集裝.如果如圖所是采用各種不同的單元（一端有鉸），則定位向量為①②③④①②③④如何獲得帶鉸的單元?jiǎng)偠染仃嚭偷刃Ш奢d矩陣定位向量①②③④①②③④4.1桿系結(jié)構(gòu)整體分析4-1)剛度集裝以3單元為例來(lái)說(shuō)明定位向量單元局部位移碼根據(jù)單元局部位移碼和定位向量的對(duì)應(yīng)關(guān)系用定位向量位移碼送元素。根據(jù)單元局部位移碼和定位向量的對(duì)應(yīng)關(guān)系用定位向量位移碼送元素，定位向量元素為零時(shí)不送。①②③④①②③④4.1桿系結(jié)構(gòu)整體分析4-2)荷載集裝以4單元為例來(lái)說(shuō)明定位向量局部位移碼此結(jié)論同樣適用于剛度集裝4.1桿系結(jié)構(gòu)整體分析4.1.7整體分析總結(jié)1)對(duì)局部坐標(biāo)和整體坐標(biāo)不一致的單元，要對(duì)剛度、荷載進(jìn)行坐標(biāo)轉(zhuǎn)換。2)需對(duì)“結(jié)構(gòu)”進(jìn)行結(jié)點(diǎn)、位移的局部和整體編號(hào)。4)集裝所得整體剛度矩陣是對(duì)稱(chēng)、帶狀稀疏矩陣，當(dāng)支撐條件能限制剛體位移時(shí)，矩陣非奇異。3)根據(jù)單元局部位移碼和定位向量的對(duì)應(yīng)關(guān)系用定位向量位移碼送元素，定位向量元素為零時(shí)不送。據(jù)此可集裝、累加得到整體剛度矩陣。5)綜合荷載由兩部分組成，因此首先要將直接作用結(jié)點(diǎn)的荷載按結(jié)點(diǎn)位移碼送入，如果還有單元等效荷載，再按定位向量集裝、累加。4.1桿系結(jié)構(gòu)整體分析8)如果有某位移碼方向彈性支撐，需進(jìn)行將彈簧剛度送入位移碼對(duì)應(yīng)的對(duì)角線元素位置累加。9)如果有某位移碼方向已知支撐位移，需進(jìn)行將“邊界條件處理”。具體做法以后介紹。7)整體剛度方程實(shí)質(zhì)是全部結(jié)點(diǎn)的平衡條件。6)剛度矩陣帶狀稀疏，其帶寬取決于結(jié)點(diǎn)、位移編碼。最大半帶寬=定位向量中最大元素差+1。4.1桿系結(jié)構(gòu)整體分析2.5.8邊界條件的處理10)當(dāng)用虛位移或勢(shì)能原理作整體分析時(shí)（勢(shì)能為例），應(yīng)變能為單元應(yīng)變能之和，外力勢(shì)能為單元外力勢(shì)能之和+結(jié)點(diǎn)外力勢(shì)能。全部桿端力的外力勢(shì)能彼此抵消。1)乘大數(shù)法設(shè)第i個(gè)位移為已知值a，N=108或更大的數(shù)。乘大數(shù)法是將剛度矩陣Kii改為NKii，將Ri改為

Na。請(qǐng)考慮為什麼這樣做能使邊界條件得到滿足？2)置換法（劃零置1）設(shè)第i個(gè)位移為已知值a。4.1桿系結(jié)構(gòu)整體分析上述置換工作量大一些，顯然可看出邊界條件得到精確滿足。4.1桿系結(jié)社構(gòu)整體宅分析★3)關(guān)于斜叢邊界的短處理如圖賢示意都的斜扛支座耗情況主，有卻多種御處理孝方案亡。3-1撐)通過(guò)單具元的坐材標(biāo)轉(zhuǎn)換睛來(lái)處理xy圖示賠有斜償支座懂單元劣，r結(jié)點(diǎn)窩處以?xún)A孔角-來(lái)進(jìn)跌行坐頃標(biāo)轉(zhuǎn)元換，紐奉也即在r結(jié)點(diǎn)處皆整體坐制標(biāo)為圖斑示xy。r3-遺2)通過(guò)增雞加一個(gè)洽單元來(lái)霧處理圖示有槽斜支座伴單元，r結(jié)點(diǎn)處補(bǔ)沿y方向把增加強(qiáng)一個(gè)剛結(jié)的枕單元，此單元有“無(wú)炎窮大絕”的江抗拉紅剛度撒、但沒(méi)有抗類(lèi)彎剛度適。單元戰(zhàn)長(zhǎng)度可艱任意。3-鴉3)對(duì)整體浮剛度矩丑陣進(jìn)行薄處理（認(rèn)參見(jiàn)匡采文起教薪材）最大半帶寬2.5桿系結(jié)軌構(gòu)整體嗎分析2.抽5.鞏10總剛度聽(tīng)矩陣的悅存儲(chǔ)與蹤蝶對(duì)應(yīng)解回法因?yàn)猷徔倓偘染叵ш噷?duì)拍稱(chēng)、噸帶狀疾稀疏梳，利附用這歲一特點(diǎn)可減今少存儲(chǔ)苗、提高孤解算效隊(duì)率。零元素零元脫素非零元素最大半帶寬主對(duì)角線元素等半帶存儲(chǔ)零元素非零元素變帶寬一維存儲(chǔ)帶寬是變的到P牌.5祥52.5桿系支結(jié)構(gòu)讀整體業(yè)分析目前一腰般都用單變帶寬砍存儲(chǔ)，蠟下面結(jié)添合程序籌說(shuō)明存儲(chǔ)引和解償法。咐首先醬介紹呼一些F9干0的語(yǔ)法稻。定義進(jìn)導(dǎo)出型類(lèi)型導(dǎo)出類(lèi)辰型——結(jié)點(diǎn)typ光e:緞:typ洪_(kāi)Jo輝intrea旁l色::口x,y！坐標(biāo)in槽te壩ge鉛r亦::尸G服DO蕩F(死3)！整體混位移忌碼end驕ty惱pety傳p_洋Jo殊in衫t1)有關(guān)F90語(yǔ)法導(dǎo)出類(lèi)攝型新特憤性用結(jié)禁點(diǎn)導(dǎo)哭出類(lèi)肥型作州為成奏員導(dǎo)勒出單近元類(lèi)爪型:typ厘e(cuò):忙:typ社_El表eme算ntint嬸ege觀r陶::輸Joi膝ntN副o(2夠)！結(jié)點(diǎn)編姥號(hào)ty頸pe郊(t眨yp冰_J搜oi極nt)牽::右No孟de(案2)！結(jié)點(diǎn)int蝕ege逝r范::滴Glb午DOF黎(6)認(rèn)!定位盟向量re污al弄::秩E賄A,絞EIend怎ty仙pety始p_丙El達(dá)em畝en妙ttyp掏e(ty晨p_鋪El助em趨en漢t)過(guò)::瀉Ele亡m(5及)！定義5個(gè)單斥元類(lèi)休型…!對(duì)單良元i的端孤點(diǎn)j的x,y凳,GD區(qū)OF(譯1:3攔)的賦值Ele份m(i艦)%N庭ode潛(j)=Joi以nt(腔Ele蹤蝶m(i縫)%J膝oin廢tNo碰(j))…由導(dǎo)出沉類(lèi)型定百義新類(lèi)葬型由導(dǎo)出白類(lèi)型定肅義變量re肺al沖:蠟:互A(笨5)詠,B福(5口,1辟0)水,C股(5從)B=0窩.0!對(duì)B清零A=蝦1.汽0!對(duì)A賦1:勾A出(i琴)=般1.料0,伙i絞=1引,5C=圈A+招2!數(shù)組貌與標(biāo)慣量運(yùn)誤算:A堤(1:執(zhí)5)+忠(/搶2,2合,2,企2,2柴/)A=綱C+貪A!數(shù)組昆與數(shù)亞組運(yùn)等算（鹿同形貼）C=sq證rt喝(A)!數(shù)組徑的函捉數(shù)運(yùn)殃算:C寇(i)贈(zèng)=sq醒rt頌(A稠(i),頌i=1沖,5數(shù)組內(nèi)塵部函數(shù)：do矛t_誼pr摸od程uc床t(存ve謠ct全or圾_a僅,v萄ec膽to夫r_賣(mài)b)!點(diǎn)積如：dot召_pr鄉(xiāng)豐odu航ct(鳳(/1襯,2,穗3/)沾,(/焦2,3甚,4/妹))的值為20（待價(jià)續(xù)）有關(guān)F9水0語(yǔ)法數(shù)組運(yùn)真算與賦業(yè)值：ma塞tm恩ul糾(m紅at壞ri弓x_姑a,吹ma碰tr劉ix睜_b)矩陣次相乘如：lo集cE盟Di濱sp=mat僑mul僚(T,度glb粱EDi矩sp)tra兔nsp落ose船(ma告tri揀x)矩陣轉(zhuǎn)泉置如：glb犁EDi淺sp=ma食tm儲(chǔ)ul虎(t擇ra冠ns占po鼓se浙(T亮),晚lo好cE萌Di極sp)siz礦e(a死rra個(gè)y,d洪im)求數(shù)組橋第di俊m維的漢長(zhǎng)度dim為可昌選變框元:si內(nèi)ze捷(a住,d糟im兵=2調(diào))若arr繪ay為一維栽時(shí)，可小不用di姻m。sum糠(ar叛ray襪,di初m,m席ask嘴)數(shù)組元揚(yáng)素求和di蘇m,匠ma岔sk為可鹽選變侵元；mas歐k=條件表犬達(dá)式sum軟(a(鑼1:1爽0))對(duì)a的1到10元素求容和su軍m(午a(突1:淡10悼),階ma伶sk追=a折>0摸)對(duì)a(1劑:10臺(tái))中大于0的元咬素求包和（續(xù)）有關(guān)F9見(jiàn)0語(yǔ)法whe位re結(jié)構(gòu)新特鐵性例whe幅re乓(C>科0)C=燦0A=柴B*賊Dend目wh臭erewhe摸re轟(C>焰0)A恢=壇Ben檢d跡wh辣er搞e定義wh滲er禾e(數(shù)組關(guān)六系表達(dá)防式)數(shù)組齒賦值姨語(yǔ)句…els逆ewh鎮(zhèn)ere數(shù)組泛賦值社語(yǔ)句...end蒼wh駕ere規(guī)則:1)同形數(shù)犁組;2擺)不許環(huán)嵌套;3材)最多二洗分叉有關(guān)F9嫂0語(yǔ)法cy肯cl矛e和exi汗t語(yǔ)句新特或性用在do循環(huán)穗中cy喪cl肢e——作下一規(guī)個(gè)循環(huán)收步exi廟t——跳出澇循環(huán)墳，執(zhí)有行en鑼d接do后一哈條語(yǔ)眨句等效例do...if倒(.洲no停t.爛co艷nd寄1)區(qū)th覆en...if(跨con本d2)got宏o5..止.end舟ifen遞d售do5柜..鞏.用法do..遙.if骨(咸co仙nd孕1)票c頭yc透le...if馳(co克nd2彎)e甩xit...end探do...有關(guān)F9大0語(yǔ)法數(shù)組捏構(gòu)造漢函數(shù)spr您ead語(yǔ)法spr杰ead歪(數(shù)組名,di巡壽m,吼nc萄op叨ie該s)將數(shù)組刪沿dim維方丘向復(fù)慚制nc吉op咽ie異s形成新張數(shù)組di律m,粗nc艇op死ie摔s—整型、暑位置變倚元、關(guān)掀鍵字變斗元（若無(wú)按位衫置引津用，汽可略奸關(guān)鍵秩字）例:（僅輝限一其維數(shù)盒組）1）spr輸ead圍(on焰e,d君im=泄1,n窩cop艦ies鐵=3)spr銳ead葬(on草e,1粒,3)spr救ead浴(on顆e,n羊cop犧ies=3,趙dim嚷=1)[1,1,1]或[1,1,1]T2）ELo討cVe盤(pán)c(1跨:6)舟=(/漿1,腦0,3絞,4,規(guī)5,0裂/)spr暮ead懂(EL左ocV低ec,逐dim=1,串nco察pie皇s=3舉)3）spr惜ead眾(A(吊2,2拳:),取dim具=1,垂nco厘pie梨s=2叔)如果di燥m=景2呢?有關(guān)F9丸0語(yǔ)法指針po慨in糖te岡rpo趁in叉te搭r是變量搬的屬性眼，可以旬指向相社同類(lèi)型謎的變量;被指向敘的目標(biāo)塊必須具寄有ta袍rg明et屬性或poi乎nte課r屬性可以將話指針變注量理解半為別名賠、稱(chēng)號(hào)re州al閑,t鋸ar經(jīng)ge考t棒:李:肚a,并b,殃ED憲is焰p(鞏6)!可被略指針蜘所指re竿al視,p廁oi怪nt陶er誘:禿:養(yǎng)p1厭,p侄2!稱(chēng)號(hào)補(bǔ)：班彼長(zhǎng)、芝課代聯(lián)表！p1丘,p物2是指針嶼，可以采指向?qū)嵳n型數(shù)據(jù)re永al嚇,p含oi造nt夾er堪:樸:他G(晉:)!先進(jìn)集承體！G是指君針，少可以舍指向落一維開(kāi)實(shí)型論數(shù)組指針融是一釀種“稱(chēng)號(hào)”，上丈述聲抗明語(yǔ)珠句建制立了“稱(chēng)號(hào)驕”，但并培未“授意予”某個(gè)嗚變量續(xù)這個(gè)性稱(chēng)號(hào)貌，因但此是典指向“空”，并未增占用內(nèi)僅存a=性3.賴(lài)0p1牙=>哭a!p身1指向a!稱(chēng)號(hào)p1授予a，a的數(shù)據(jù)盾有兩個(gè)功名：固示定名a和流動(dòng)緣瑞名p1逃;既可用p1也可盤(pán)用a（p1丙—班長(zhǎng)，a立—張三）a=飾4.榨0!者a的值變館為4.咬0,p1也變米為4.0p1編=職>朱b!班長(zhǎng)換再人了G丑=>ED扭is鵲p!先進(jìn)集漸體有了鹽得主！ED芝is齒p(:)的長(zhǎng)度停就是G(:飯)的長(zhǎng)淋度，賺用G和用EDi睛sp同樣效匹果又如：rea晃l,t若arg容et伐::慰a,brea菌l,p糕oin毀ter名::鉗p,qa=指3.衛(wèi)14b=驕2.愚0p究=>去a!丹p哨=沿a頓=討3.工14q=瞎>b!沿q岔=秧b膛=鑒2.啊0q壩=橫p!晝(q指向頂?shù)臄?shù)奔據(jù)b)=崇(p指向的每數(shù)據(jù)a梁)！即：嶺所有=3籮.14指針集可以晃指向奪有名映的數(shù)態(tài)據(jù)區(qū)厭，也捆可以拼指向原無(wú)名堡的數(shù)戚據(jù)區(qū)rea舅l,p護(hù)oin肺ter筒::洗b(:迷,:)in杠te帖ge計(jì)r罩::新nre騙ad鼠(*粒,*拌)牲nal拜lo都ca畫(huà)te法(津b(鈔n,演n)啟)!指針蔽指向閑了一逼個(gè)剛角開(kāi)辟乳的數(shù)敘組!以下可瀉以當(dāng)作陽(yáng)常規(guī)數(shù)福組用b(奸1,喪1)織=情1腿.1b(1機(jī),2)泉=索1.2...dea廉llo脾cat拿e(b嗚)有關(guān)F90語(yǔ)法用指封針建堅(jiān)立動(dòng)普態(tài)數(shù)倡組指針三與al摘lo像ca堅(jiān)ta納bl白e數(shù)組的描區(qū)別具備al管lo雄ca禍ta停bl丟e數(shù)組的咳所有功驕能還可以洪用在導(dǎo)到出類(lèi)型株中，例匯如整體坦剛度矩彎陣的變挖帶寬存境儲(chǔ)：ty茫pe消:親:ty綢p_潑Kc衛(wèi)ol!整體剛角度矩陣K的列rea材l(8踏),p尿oin旦ter予::筍ro妙w(:燙)!該列俘的行判元素en夾d睬ty積pe..擴(kuò).typ蠟e(ty播p_若Kc剛ol間),惑al睬lo跪ca堵ta尼bl謠e::Kc臉ol(:)…all抓oca諷te紛(Kc擾ol真(N疼Gl巴bD刮OF))!分配了NG就lb暴DO活F列...all倍oca階te灘(Kc浴ol(厭5)%萄row搏(3:彩5))!第5列只用3至5行(1奧)NE晝le吵m,裙NJ疾oi害nt,NG松lb計(jì)DO面F,NJL患o(jì)ad,NE永Lo浪ad單元數(shù)畫(huà)結(jié)點(diǎn)逼數(shù)總雨自由度貴數(shù)結(jié)凝點(diǎn)荷載肚數(shù)單箭元荷載雷數(shù)[Jo什int然-結(jié)點(diǎn)]…NJ雨oi誤nt行(2)綱Jo尚int壤%X,摘Joi謎nt%那Y,這GDO介F(1盆:3)X坐標(biāo)Y坐標(biāo)睜結(jié)嶺點(diǎn)位移旺碼[El烘em蠟-單元]…NE肝le果m行(3)歌Jo醋int加No1設(shè),Jo體int峰No2卷,E晨A,E償I起點(diǎn)沃號(hào)外終點(diǎn)品號(hào)次剛度[JL俱oa鵝d-結(jié)點(diǎn)荷紙載]…NJ惕Lo植ad行(4)Jo交in易tN還o,Lod右DOF,Lod的Val作用點(diǎn)援號(hào)局懷部位移陰碼荷翠載值[EL本oa胞d-單元荷哥載]…NE悟Lo競(jìng)ad行(5)El游em慢No,Ind攔x,印Po爸s,Lod憶Val單元號(hào)睡類(lèi)型躁號(hào)位僅置荷勇載值Indx

類(lèi)型

pos1均布長(zhǎng)度2集中位置3...2)某程序蓮輸入數(shù)態(tài)據(jù)說(shuō)明3,5爹,7,游1,10,0鞋,0即,0,董00,載4,算1深,2頃,34,4康,4踐,5,的64,4半,4線,5,撞74,0挪,0夫,0,希01,2遷,1懇.0E遲9,季162,3萌,1曲.0E擱9,混245,攔4,友1哀.0蜂E9剛,防122,傻1,向10匯.0槐E3!結(jié)點(diǎn)2，自歐由度1，值為10E跌32,或1,扎4,哥-4徹.0察E3！單元2，均匯布，吊長(zhǎng)4米，值大為-4餡E32-1尖)數(shù)據(jù)親文件捐例子宰：(2)(1)(3)24135i=6i=4i=310kN4kN/m4m4mEA=恢109N(1編)(2)坐標(biāo)位移宣碼(3贈(zèng))(4)(5)結(jié)點(diǎn)煙號(hào)EA,曾EIrea薄d(5釀,*)NEl罰em,延NJo捐int尖,NG宴lbD墻OF,斧NJL飽oad磚,NE繁Loa腫dall及oca濤te建(Jo耗in疊t(徑NJ回oi換nt))all庭oca朽te排(Ele觸m(N鞏Ele糖m))all數(shù)oca林te拳(JLo泰ad(域NJL詢(xún)oad))all配oca貞te芹(EL腎oa帥d(經(jīng)NE盆Lo萄a(bǔ)d))..錄.rea惑d(5扯,*)蹄(J字oin勻t(i揉),i狠=1,舒NJo摸int歉)re剝ad脖(5繞,*畜)驗(yàn)(Ele欺m(i關(guān)e)%襪Joi針ntN撞o,E裝lem厲(ie求)%E左A,&ELe壟m(i嗽e)%籮EI,出ie=1,晨NEl遞em)if桐(NJ后Lo咳ad>0)羽re汗ad(狹5,*京)(JL傷oa性d(活i)棚,i=1宏,N矛JL堆oa農(nóng)d)if巡壽(NE勁Lo捕ad>0)冒re岔ad(若5,*果)(ELo雁ad(為i),俘i=1,供NEL炒oad恒)2-酸2)程序讀鐮入計(jì)算病所需數(shù)跌據(jù)：2-3樓)求始行床碼和分扯配帶寬邪子程序!=春==紫==就==改==猶==坊==血==皮==隨==慕==杰==援==故==壯==返==拾==賤=sub塑rou仁tin老eSet宮Mat津Ban寇d(Kc幸ol,E掙lem左)!接口里簡(jiǎn)單!=劃==幼==鴿==鴨==歲==坐==進(jìn)==前==燒==斑==姥==民==遵==斷==傷==循==扛=typ致e(t牧yp_雕Kco麥l),延int哥ent曾(inou序t)缸:室:Kco判l(wèi)(:)！總剛列ty湖pe擴(kuò)(t劫yp款_E疏le躬me皆nt予),粱in持te裕nt掠(i晚n))湖:圓:濫El敏em叮(:霉)！單元int重ege沾r(i燈kin圍d)夕::扒m檢in羊DO樣F,盞EL樓oc舉Ve掠c(開(kāi)6)in該te賴(lài)ge噸r(位ik于in盟d)詢(xún)::諸R似ow耍1(菜si截ze哪(K輸co縮慧l,繼di遠(yuǎn)m=諸1)扣)！Row懼1為自練動(dòng)數(shù)畜組，揉子程嶺序結(jié)碗束后摘自動(dòng)乒釋放債。！這脈樣做壤可使敬接口候簡(jiǎn)單快，減灣少了溜數(shù)組賣(mài)的控慈制變籠量。int助ege肯r(i悶kin妻d)溜::ie每,j飾,N階Gl浩bD售OF做,N殃El哈emNG宿lb己DO抬F=siz濫e(K居col區(qū),di初m=1)！使接口王簡(jiǎn)單NE扔le申m=siz日e(E序lem齒,di住m=1岔)Ro熱w1踢=NG吩lb營(yíng)DO植F!先設(shè)所科有始行滴碼同終悲行碼!確定各寫(xiě)列始行僑碼向量doie=1列,N漿el盾em!對(duì)單元穗循環(huán)!確定摘定位陡向量EL陸oc缺Ve柔c(:龍)很=El此em守(i摩e)捷%G錫lb冰DO習(xí)F(:)!尋找蝦定位網(wǎng)向量左中大品于零罷的最桂小值min嬸DOF=min宿val手(EL索ocV釘ec,皮mas城k=ELo省cVe衰c>0膨)!屏蔽徒定位用向量烏中小庸于等薯于零泥的wh姻er被e隱(ELo點(diǎn)cVe父c>腔0)!尋找Ro番w1匆(E餓Lo嗎cV拔ec多)和mi逆nD感OF中的最紅小值Row盾1(E咬Loc等Vec粱)=蒼mi宣n(R南ow1慕(EL誰(shuí)ocV偏ec)蘇,&mi揉nD鉗OF)end槐wh澆ereen默d梁do!為各列創(chuàng)的帶寬猶分配空武間do脊j看=1嚇,N洽Gl槽bD比OF!對(duì)總自氣由度數(shù)呀循環(huán)all味oca鍋te醬(Kc淹ol(卵j)%憂row河(Ro薪w1(讀j):腥j)!給Kc擔(dān)ol傭(j事)%奶ro螞w從Ro驕w1難(j績(jī))到j(luò)個(gè)空陜間Kco刪l(j涉)%r巷ow=z市ero!總剛元田素清零end蘿doret灶urnend勤su抵bro急uti炎neSet披Mat棒Ban娛d3)整體剛游度矩陣邊的集成doie=1,不NEl椅em…!計(jì)算單阻剛EK蒼(6勁,6種),沫EL鈔oc蜘Ve劉c(產(chǎn)6)do反j=1看,6!對(duì)單元截逐列集巖成JGD耍OF涼=EL裕oc測(cè)Ve艘c(談j)!取出位監(jiān)移碼if戒(JG言DOF神=丹0)飛cyc游le!作下一碼循環(huán)步whe快re鄙(ELo厲cVe醋c>0臥.娃And翠.ELo遠(yuǎn)cVe虧c<=禿J悼GD增OF斬)!位移煮碼非僑零同旁時(shí)小荒于第j個(gè)局孟部碼欠對(duì)應(yīng)爪的位頂移碼Kco尚l(wèi)(J雨GDO鬧F)%殃row閘(EL迎ocV申ec)葬=構(gòu)&Kc師o(wú)l墨(J剖GD軋OF隸)%灶ro驢w(蔬EL衡oc趨Ve衛(wèi)c)答+你EK胖(:把,j幕)!集成en排d暫wh能er犧eend灰doen挖d隆do局部碼位移碼4)變帶寬猴矩陣的張分解求鋤解4-勵(lì)1)LD田L(fēng)T分解法藝求解[A]{x}=若對(duì)稱(chēng)正定，則可分解為=單位歌上三姨角陣對(duì)角傍陣原方慣程變肌為求解步政驟：1.分解：2.解y：3.解x：LD臉LT分解法Gau烤ss消去前消去處理縮慧右端昏項(xiàng)回代(不同的b只做弱一次螞分解)存放粗：主對(duì)角—上三壁角—時(shí)不必后求和(上三怎角:i<j)不動(dòng)存在處存在處(第j列系數(shù))上三堂角:i<j4-育2)分解膠一般慣情況歷：逐列忠分解對(duì)角線:i=j新>1(第j列系鼻數(shù))第i列中第1個(gè)非岸零元壩素的嚼行碼禍為：4-乖3)變列寬夜存貯的降分解修腎正：第j列中魯?shù)?個(gè)非零末元素的謙行碼為封：則分解調(diào)順序4-藥4)仍F擱90實(shí)現(xiàn)!三角躲分解dia絲式g(1霞:nc甘ol)計(jì)=好(/鈔(Kc摸ol繡(j涌)%雷ro建w(鍋j)季,j=1,灣nco駛l)芬/)do屑j=懶2,幼nco正lrow韻1=菊l(shuí)b閉oun膏d(K饑col鵲(j)宏%ro刮w,1欠)!i航_1較jdo雨i伶=銜r(shí)更ow偽1,倆j-府1ro胡w_斧1州=夏ma截x(消ro晝w1驢,l姑bo多un彈d(殼Kc古ol臂(i境)%配ro煙w,折1)述)!攝i_阿1s歸=奮su歐m(罰di竭ag妨(r性ow豪_1譽(yù):i槐-1癢)戰(zhàn)*K摧co振l(園i)似%r序ow互(r濾ow束_1們:i擇-1雅)房誠(chéng)*袖&Kc甲o(hù)l紐奉(j堆)%尤ro貢w(脆ro火w_饞1:防i-遙1)演)!求和詞部分Kco便l(j汁)%r池ow(抬i)=笑(Kco賢l(j忠)%r抄ow(沉i)-燙s)/緊dia閱g(i)end君do！第j列系數(shù)殃分解完低畢s=廈su辭m(d柱iag壩(ro才w1:臭j-1掘)*K歪col決(j)普%ro肉w(r屈ow1部:j-廉1)*血*2)dia樓g(j)=dia蜻g(j)-輪s!第j列的主騰對(duì)角元趙素end顯do4-5最)一般揮情況