【正文】 面高度；——截面有效高度，；——混凝土抗壓強(qiáng)度設(shè)計(jì)值， 、——鋼筋強(qiáng)度， ——受拉、受壓區(qū)鋼筋面積，由規(guī)范鋼筋表查得； 、——自鋼筋的重心分別至截面最近邊緣的距離；構(gòu)件大偏壓破壞的概率為： 構(gòu)件小偏壓破壞的概率為： clear。在Elem table item at node I欄后面的下拉菜單中選取FI，在Elem table item at node J欄后面的下拉菜單中選取FJ，在Optional scale factor后面欄中輸入1，在Items to be plotted on欄后面選擇Deformed shape單選按鈕，最后單擊OK按鈕，得到隧道襯砌支護(hù)結(jié)構(gòu)的軸力圖，如圖330所示。單擊Add按鈕，打開(kāi)一個(gè)Define Additional Element Table Items對(duì)話框，如圖326所示。選擇菜單Main Menu：Preprocessor→Modeling→Delete→Nodes命令，打開(kāi)一個(gè)刪除節(jié)點(diǎn)選取對(duì)話框，在輸入欄中輸入11111111111111112122，單擊OK按鈕。隧道本身是受豎向均布?jí)毫Γ鐖D319，在此將豎向均布荷載q置換為每個(gè)節(jié)點(diǎn)的節(jié)點(diǎn)荷載，節(jié)點(diǎn)荷載等于均布荷載的集度q乘上節(jié)點(diǎn)相鄰梁?jiǎn)卧酵队暗?/2，再乘上結(jié)構(gòu)的計(jì)算寬度b，這里，b=1。在DOFs tu be constrained欄后面選取UX，UY，接著在Apply as欄后面的下拉菜單中選取Constant value選項(xiàng)，并在Displacement value欄后面輸入0值，然后單擊OK按鈕，完成對(duì)彈簧節(jié)點(diǎn)位移的約束。 element divisions欄后面輸入1，再單擊OK按鈕。創(chuàng)建隧道襯砌支護(hù)線模型圖38隧道襯砌支護(hù)線模型選擇菜單Main Menu：Preprocessor→Modeling→Create→Lines→Lines→straight line命令，然后依次選中1，2兩點(diǎn)，再依次選中2，3兩點(diǎn)，直到82，1兩點(diǎn)，單擊OK按鈕。這是因?yàn)樗淼酪r砌支護(hù)仰拱和腰部以及頂部的厚度不同，所以要建立多個(gè)BEAM3實(shí)常數(shù)。⑵ 選擇菜單Utility Menu：File→Change Title命令，出現(xiàn)Change Title對(duì)話框，在輸入欄中輸入Tunnel Support Structural Analysis，單擊OK按鈕關(guān)閉該對(duì)話框。圖129結(jié)構(gòu)軸力圖（單位：N）⑶ 列出主要數(shù)據(jù)1）列表顯示單元的彎矩，剪力和軸力 圖130 List Element Table Data對(duì)話框選擇菜單Main Menu：General Postproc→List Results→Element Table Data命令，打開(kāi)一個(gè)List Element Table Data對(duì)話框，如圖130所示。在User label for item欄后面輸入MI，在Item，Comp Results data item欄后面的左邊下拉菜單中選取By sequence num，并在右欄輸入6，然后單擊Apply按鈕；再次在User label for item欄后面輸入MJ，在Item，Comp Results data item欄后面的左邊下拉菜單中選取By sequence num，并在右欄輸入12，然后圖125定義單元數(shù)據(jù)表對(duì)話框單擊Apply按鈕；使用同樣的方法依次輸入FI，1；FJ，7。3）第二次求解選擇菜單Main Menu：Solution→Solve→Current LS命令，打開(kāi)一個(gè)求解選項(xiàng)信息和當(dāng)前求解載荷步對(duì)話框，接受默認(rèn)設(shè)置，單擊OK按鈕，開(kāi)始求解計(jì)算，直到出現(xiàn)一個(gè)Solution is done提示欄，單擊Close按鈕關(guān)閉提示欄。第八步：后處理（對(duì)計(jì)算結(jié)果進(jìn)行分析）⑴ 計(jì)算分析修改模型1）查看隧道襯砌支護(hù)結(jié)構(gòu)變形圖選擇菜單Main Menu：General Postproc→Plot Results→Deformed Shape命令，彈出Plot Deformed Shape對(duì)話框，如圖120所示，選中Def+underformed，單擊OK按鈕，將出現(xiàn)隧道襯砌支護(hù)結(jié)構(gòu)變形圖，如圖121所示。選擇菜單Main Menu：Solution→Define Loads→Apply→Structural→Force/Moment→On Nodes命令，打開(kāi)節(jié)點(diǎn)位置施加載荷對(duì)話框，以41點(diǎn)為例，選擇41點(diǎn)，打開(kāi)打開(kāi)Apply F/M on Nodes對(duì)話框。這里左側(cè)加至39點(diǎn)，右側(cè)加至45點(diǎn)，再根據(jù)判斷去掉多余彈簧，由于此處數(shù)據(jù)量較大，建議用命令流添加。 element divisions欄后面輸入1，再單擊OK按鈕。得出二襯的中心線如圖16所示。按相同步驟添加其他幾個(gè)。在【/FILNAM】Enter new Jobname輸入框中輸入工作名稱(chēng)Support，單擊OK按鈕關(guān)閉該對(duì)話框。⑵ 選擇菜單Utility Menu：File→Change Title命令，出現(xiàn)Change Title對(duì)話框，在輸入欄中輸入Tunnel Support Structural Analysis，單擊OK按鈕關(guān)閉該對(duì)話框。第四步：定義材料屬性⑴ 選擇菜單Main Menu：Preprocessor→Material Props→Material Models命令，出現(xiàn)Define Material Behavior對(duì)話框。(3)保存幾何模型文件選擇菜單Utility Menu：File→Save as命令，打開(kāi)Save Database對(duì)話框，在Save Database ，單擊OK按鈕。⑹ 劃分網(wǎng)格圖111二襯中心線節(jié)點(diǎn)圖在網(wǎng)格劃分工具欄中單擊Mesh按鈕，打開(kāi)一個(gè)對(duì)話框，單擊Pick All按鈕，生成82個(gè)梁?jiǎn)卧?，如圖111所示。圖113彈性系數(shù)原理添加彈簧后的圖像如圖114。在Direction of force/mom欄后面的下拉菜單中選取FY，在Force/moment 。圖121初次分析計(jì)算隧道襯砌支護(hù)結(jié)構(gòu)變形圖圖120查看變形圖對(duì)話框 由上圖可知有些彈簧為受拉彈簧，這意味著該處土體受拉，這顯然不符合常理，需要去除彈簧并重新計(jì)算，直到結(jié)構(gòu)變形圖中沒(méi)有受拉彈簧為止。4）查看第二次分析計(jì)算結(jié)構(gòu)變形圖選擇菜單Main Menu：General Postproc→Plot Results→Deformed Shape命令，彈出Plot Deformed Shape對(duì)話框，選中Def+underformed，單擊OK按鈕，將出現(xiàn)第二次分析計(jì)算的隧道襯砌支護(hù)結(jié)構(gòu)變形圖。最后得到定義好后的單元數(shù)據(jù)表對(duì)話框，如圖126所示。在Items to be listed欄后面的下拉菜單中選擇MI、MJ、FI、FJ選項(xiàng)，然后單擊OK按鈕，打開(kāi)單元數(shù)據(jù)表文件，如圖131所示。(3)首先選擇主菜單parametersangular units,在AFUN一欄選擇Degrees DEG，使ansys可以辨認(rèn)角度。圖35定義隧道腰部和頂部BEAM3實(shí)常數(shù)2對(duì)話框按相同步驟添加其他幾個(gè)。得出如圖38所示2 保存幾何模型文件選擇菜單Utility Menu：File→Save as命令，打開(kāi)Save Database對(duì)話框，在Save Database ，單擊OK按鈕。3 劃分網(wǎng)格圖313隧道支護(hù)單元圖在網(wǎng)格劃分工具欄中單擊Mesh按鈕，打開(kāi)一個(gè)對(duì)話框，單擊Pick All按鈕，生成82個(gè)梁?jiǎn)卧鐖D313所示。圖317橫向豎向約束另外給最低點(diǎn)1點(diǎn)施加水平方向約束。例如，第1點(diǎn)的節(jié)點(diǎn)荷載p為，第2節(jié)點(diǎn)的節(jié)點(diǎn)荷載p2為：圖219施加荷載圖解同時(shí)點(diǎn)2所受的力為1桿的j端里和2桿的i端力，兩者相疊加。3）第二次求解選擇菜單Main Menu：Solution→Solve→Current LS命令，打開(kāi)一個(gè)求解選項(xiàng)信息和當(dāng)前求解載荷步對(duì)話框，接受默認(rèn)設(shè)置，單擊OK按鈕，開(kāi)始求解計(jì)算，直到出現(xiàn)一個(gè)Solution is done提示欄，單擊Close按鈕關(guān)閉提示欄。在User label for item欄后面輸入MI，在Item，Comp Results data item欄后面的左邊下拉菜單中選取By sequence num，并在右欄輸入6，然后單擊Apply按鈕；再次在User label for item欄后面輸入MJ，在Item，Comp Results data item欄后面的左邊下拉菜單中選取By sequence num，并在右欄輸入12，然后單擊Apply按鈕；使用同樣的方法依次輸入FI，1；FJ，7。⑶ 列出主要數(shù)據(jù)1）列表顯示單元的彎矩，剪力和軸力 選擇菜單Main Menu：General Postproc→List Results→Element Table Data命令，打開(kāi)一個(gè)List Element Table Data對(duì)話框，如圖331所示。close all。ES=。h0=h01(i)。syms b f h0 N fy1 As1 ES fy As as as1 M h MUX sigMAXb=1。7 個(gè)變量分別為，1軸力；2fy；3fy1；4厚度h0；5混凝土抗壓強(qiáng)度f(wàn)；6彎矩；7，h。normX=eps。1/2*(x(1,1)x(2,1)*As1+x(3,1)*As)^2/b/x(5,1)^2。close all。MUX=[476463 479323 484990 493049 503047 514320 526069 537564 547851 556438 563134 567558 569375 568332 565896 565777 564040 559800 552502 541680 529327 542166 553763 546051 526246 531923 537976 542756 545846 547415 547476+08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08+08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +08 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 +07 60874 60984 26098 101203 105690 ]。normX=eps。...... x(4,1)*(x(1,1)x(2,1)*As1+(33/10000*ES+x(3,1))*As)/(x(5,1)+5/4*(33/10000*ES+x(3,1))*As/x(4,1))x(4,1)*x(5,1)*(x(1,1)x(2,1)*As1+(33/10000*ES+x(3,1))*As)/(x(5,1)+5/4*(33/10000*ES+x(3,1))*As/x(4,1))^21/2*(x(1,1)x(2,1)*As1+(33/10000*ES+x(3,1))*As)^2/(x(5,1)+5/4*(33/10000*ES+x(3,1))*As/x(4,1))^2+x(5,1)*(x(1,1)x(2,1)*As1+(33/10000*ES+x(3,1))*As)^2/(x(5,1)+5/4*(33/10000*ES+x(3,1))*As/x(4,1))^3。命令流如下%%%%%%%%%%%%%%%%%%%%%%%%%%% clear。fy1=。h=h01(i)。 %%%%%%%%%%%%存放bbeta數(shù)據(jù)的矩陣%%%%%end%%%%%%%%%將該命令流文件從matlab中打開(kāi)，修改方法同鋼筋混凝土，此處函數(shù)為，只有3個(gè)變量。M1=[465 413 235 169 953 2276 4228 6693 9084 10101 7430 4044 608 2938 6217 8527 8976 6594 2519 967 633 552 193 736 2516 5206 8125 9357 5447 8245 16475 20104 19978 16992 11983 5793 861 7282 12900 17263 20014 20956]。clc。1/2*x(1,1)]。g=x(5,1)*x(4,1)*(x(1,1)x(2,1)*As1+(*ES+x(3,1))*As)/(x(5,1)+(*ES+x(3))*As/(*x(4,1)))*x(5,1)*(x(1,1)x(2,1)*As1+(*ES+x(3,1))*As)^2....../(x(5,1)+(*ES+x(3))*As/(*x(4,1)))^2+x(2,1)*As1*(x(4,1)as1)x(1,1)*(+x(6,1)/x(1,1)