【土木建筑】3桿系靜力分析4

1、subroutine elem_stiff() 說明說明 stiff=0.0 ! 單元?jiǎng)偠惹辶銌卧獎(jiǎng)偠惹辶?select case (type) case (1) 平面桿系結(jié)構(gòu)單元平面桿系結(jié)構(gòu)單元 case (2) 空間桿系結(jié)構(gòu)單元空間桿系結(jié)構(gòu)單元 case default 出錯(cuò)信息出錯(cuò)信息 end selectend subroutine elem_stiff integer,intent (in) : 入口整型參數(shù)入口整型參數(shù) real(8),intent(in) : 入口實(shí)型參數(shù)入口實(shí)型參數(shù) real(8),intent(out) : 出口實(shí)型參數(shù)出口實(shí)型參數(shù) real(8) : wor

2、k1, integer : i,j,k, 實(shí)型和整型工作變量實(shí)型和整型工作變量 select case (plane) case (1) 平面桁架元素賦值平面桁架元素賦值 case (2) 平面梁柱元素賦值平面梁柱元素賦值 case (3) case default 出錯(cuò)信息出錯(cuò)信息 end select select case (space) case (1) 空間桁架元素賦值空間桁架元素賦值 case (2) 空間梁柱元素賦值空間梁柱元素賦值 case (3) 交叉梁元素賦值交叉梁元素賦值 case default 出錯(cuò)信息出錯(cuò)信息 end select xyyxxfyfxfyfydxdy

3、dxd yxyxcossinsincosffff yxyxcossin-sincosdddd 10000cssc 00sincostsc xyyxxfyfxfyfydxdydxd t1tt )(teeeepftpf eedtd xyyxxfyfxfyfydxdydxd eeeeeeeedtktdktpftpft )( eeeedkpf tktkeet eeeedkpf 1234 2f 1f2 2dp 12f 21f21orj 各各桿桿交交ijifp1234 314d32223d21122d111dfpffpffpfp tt4dt3dt2dt1ddppppp tt32t31t22t21t12t1

4、1fffffff 2 2dp 12f 21f fiiiiiip 000000000000000000d fap d123 000000000000000000iiiiiia1234 tt4t1 tt32t31t22t21t12t11dddddd ta 1234 31kkdiagk eqdpkf;fap tt3t2t1eqpppp )(eqdpkap eqeqpap 記記 teqdakapp kaka;rpp teqd記記 kr 則則 31ttiiiiakaaka3211234 000000000000000000iiiiiia iikkkkk 22211211rssr(0,0,1)(0,0,2

5、)(3,4,5)(6,7,8)(6,7,9)(10,11,12)1234 56(0,0,1)(0,0,2)(3,4,5)(6,7,8)(6,7,9)(10,11,12)1234 56(0,0,1)(0,0,2)(3,4,5)(6,7,8)(6,7,9)(10,11,12)1234 56(0,0)(0,0)(1,2,3)(4,5)(6,7,8)12345(0,0,1)(0,0,2)(3,4,5)(6,7,8) (6,7,9)(10,11,12)1234 56(0,0)(0,0)(1,2,3)(4,5)(6,7,8)12345 55545352514544434241353433323125242

6、3222115141312113k87654543215432187654(0,0,1)(0,0,2)(3,4,5)(6,7,8) (6,7,9)(10,11,12)1234 56(0,0)(0,0)(1,2,3)(4,5)(6,7,8)1234587600 4543214 pppppp54321 nininnnjninninijiiiinjinjirrrrkkkkkkkkkkkkkkkkkkkk2121212122222211111211 akraakrakrkkkkkkkkkkkkniniininnnjnnnjnj2211212122222111121100010000 xy nnnini

7、niiinikkkkkkkkk111111136101453606672 1) 有關(guān)有關(guān) f90 語(yǔ)法語(yǔ)法導(dǎo)出類型 新特性新特性 有關(guān)有關(guān) f90 語(yǔ)法語(yǔ)法 eedtd 有關(guān)有關(guān) 有關(guān)有關(guān) 有關(guān)有關(guān) 054301054301054301 332313322212312111a. 32223222.有關(guān)有關(guān) 有關(guān)有關(guān) indx (2)(1)(3)24135i = 6i = 4i = 310 kn4 kn/m4m4mea= 109 n!=!=!=!= 621,.,j , i ij)ie(ijkk 0)(0)()ie(j)ie(iji 4) nn a tldla =單位上三角陣對(duì)角陣 bylxldl

8、 t tldla byl yxld t ulldla t nnnnijiinjnjaaaaaaaaaaaaaa2121222221111211 iidjiijll t 1111213231210nnijllllll nnjiiinndldldlddldldldd02223222221113111211111jiiiijldu ijaijl nknkjkkkikkjikijldlula11 iklikiklikik10j , ik )j , i(kikkkjkikijdllamin11)(1111da 1111ad 211112lda 121121adl 222211122dlda 221112

9、222ldad 12a22aiijiikkkjkikijdldlla 11iiikkkjkikijjiddllal 11121 j,ijjjkkkjkjjddla 112n.,j32 112jkkkjkjjjjdladiiiikkkjkikijjiddllkl 11121 j,ii11 ik 121jikkkjkjjjjdlkdii1ji1)max(111jii ,ii n.,j32 776657565546544436534333275242322241311100000kkkkdkldklldkllldlld nnnnijbbbyyyllllll2121213231211111022121

10、byyl iikikikbyyl 11kikikiiylby 11n,i21 12122ylby 11by 2bkiikikiiylby 111i iyiyibib nnnnnndydydyxxxlllll222111212321312111110nnnndyx 1nx ydxl1t nndddd11122111iiinikkkiidyxlx 1 nikkkiiiiixldyx1nni)n(i)n(ixlyy 1iiinidyy )1(nx)1( nnnyx)1(11 nnnyxnx1 nx 21 1 )(111) 1()(,n , njyxj , iixlyyjjjjjijiji pdkf

11、pdtk dtk dkmnq1f2f3f肚松衯宸&愮鐝d)? $?d悡!餯怉 扈鋹a 嘬貑 d?噡1/2001騫寸15鏈?-crm鍦氱敤.files/imgr_logo.gif 冣杁9/erp鏂噡1/2001騫寸15鏈?-crm鍦氱敤.files/logo.gif 冟?a/erp鏂噡1/2001騫寸15鏈?-crm鍦敤.files/logo_compute.gif 冡?疘1/erp鏂噡1/2001騫寸15鏈?-瀵瑰啿鍔涢噺.files/ 9/erp鏂噡1/2001騫寸15鏈?-瀵瑰啿鍔涢噺.files/0830.gif冧塖阇8/erp鏂噡1/2001騫寸15鏈?-瀵瑰啿鍔涢噺.file

12、s/4-2.gif 冨篰飁/erp鏂噡1/2001騫寸15鏈?-瀵瑰啿鍔涢噺.files/imgr_logo.gif 冨杁9/erp鏂噡1/2001騫寸15鏈?-瀵瑰啿鍔涢噺.files/logo.gif 冦旿?a/erp鏂噡1/2001騫寸15鏈?-瀵瑰啿鍔涢噺.files/logo_compute.gif 冧?疘4/erp鏂噡1/2001騫寸15鏈?-灝忕櫧榧犲拰erp.files/ /erp鏂噡1/2001騫寸15鏈?-灝忕櫧榧犲拰erp.files/0830.gif 冪阇?/erp鏂噡1/2001騫寸15鏈?-灝忕櫧榧犲拰erp.files/biaoshi.gif 冭賏?;/erp鏂噡

13、1/2001騫寸15鏈?-灝忕櫧榧犲拰erp.files/cio.gif 冩?=/erp鏂噡1/2001騫寸15鏈?-灝忕櫧榧犲拰erp.files/email.gif冩?a/erp鏂噡1/2001騫寸15鏈?-灝忕櫧榧犲拰erp.files/imgr_logo.gif 冭聖杁/erp鏂噡1/2001騫寸15鏈?-灝忕櫧榧犲拰erp.files/logo.gif 冩碝?d/erp鏂噡1/2001騫寸15鏈?-灝忕櫧榧犲拰erp.files/logo_compute.gif 冭?疘:/erp鏂噡1/2001騫寸17鏈?-鏂版湇鍔粡嫻庢誕鐜?files/ f/erp鏂噡1/2001騫寸17鏈?-鏂

14、版湇鍔粡嫻庢誕鐜?files/astd_oct.gif 冮烜?e/erp鏂噡1/2001騫寸17鏈?-鏂版湇鍔粡嫻庢誕鐜?files/biaoshi.gif 凅?c/erp鏂噡1/2001騫寸17鏈?-鏂版湇鍔粡嫻庢誕鐜?files/email.gif 冮?g/erp鏂噡1/2001騫寸17鏈?- 潕肚松衯?雎?lián)]=?牓2d?h糖 =m 痞 ?(瘣滸?( h?棧哐鋃$*?z 3蓚w 傯= %u旄 冏?塧 ?腦?u zyy h?鵓?鎂+x脲燴哠k?q 梽桃 羀徬 ? w? 錭x?c?e ?gei鎃+, nit噈鲃 c ?u锳 9 硄 輕儍l8踇縎; w岐?u rh ?wx鵛 8 吷y 兞 戀w 六

15、 堲?-猴w兂.+嬼? 墎 娮 ?嗠鸊gij倭r j ?侜! ?u?唻蕛? 凓 | ? 縖鴰嬔r 責(zé) | a伾 鼐*?% 磣嬓冴 苻 高 hv霼?m | 兪 t?4撠 ?聝 ?飌z0?t蕞筍 訌?3髀続c?p? t 壉鄪mr +?p ; 8齚 p 蕞驅(qū)_嬝3? ?%?佝較#?u嬻芺婩 *?莀x _ ? 萄阯k詴qm? x拱曃 ?鳡?+腰+?艃? ?u? 葭?橗?岶 %媗祆 ?呿祴?峙 % ms?叫 篂?x ?$ ?f笂斚汶弸邏s?苾铏p蔭淥? 呉坡y詢生蓯籀 塡塠= ?鬋 昆?sz ?|?匡?(別 ? 鯇z0 l 慶琶 5艭危?內(nèi)? g?a? j? ?$孅?0lh 効譱蘑 ? r?;n?囏a杸) g? ?t$撈-0? c+僿?v胐顑o 2 ? ?刻 n稍嶀?寃?泲? ?a礙鵠yzf猏虓贋用 拷t,v3 b勵(lì) =桽= ?st傄u烶輤 腙鯧0_儖?絅= 計(jì)裌 巹4飪秞 ? 塀墕銒t隿+廢 榰n 82 )撅 ?ug?儅導(dǎo)鴣 鮁 e 5?o軸c冾:_)5隠瑾.觠0峗e ?王?謷呻 u?蓌 x0?艣專 ?u?v=zt望 ?v ;3蒛h8#?d?z痣w兠 兯凔?黃t ? 彄貖 ?ty?嗥