【正文】 話框，如圖222所示，選中Def+underformed，單擊OK按鈕，將出現隧道襯砌支護結構變形圖，如圖323所示。3）第二次求解選擇菜單Main Menu：Solution→Solve→Current LS命令，打開一個求解選項信息和當前求解載荷步對話框，接受默認設置，單擊OK按鈕，開始求解計算，直到出現一個Solution is done提示欄，單擊Close按鈕關閉提示欄。圖324最終隧道結構變形圖5）保存計算結果選擇菜單Utility Menu：File→Save as命令，打開Save Database對話框，在Save Database ，單擊OK按鈕。在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。在Elem table item at node I欄后面的下拉菜單中選取MI，在Elem table item at node J欄后面的下拉菜單中選取MJ，在Optional scale factor后面欄中輸入1，在Items to be plotted on欄后面選擇Deformed shape單選按鈕，最后單擊OK按鈕，得到隧道襯砌支護結構的彎矩圖，如圖329所示。⑶ 列出主要數據1）列表顯示單元的彎矩，剪力和軸力 選擇菜單Main Menu：General Postproc→List Results→Element Table Data命令，打開一個List Element Table Data對話框，如圖331所示。（2） 生成可靠行文件選擇Utility MenuFileWrite DB Log File進入指定文件選項，在write database log to選項填入名稱，暫可取erchen82，如圖336圖336 保存文件（3） 進入概率設計階段圖337 assign 選擇Main MenuProb DesignAnalysis FileAssign彈出Assign Deterministic Model File選項，在analysis ，如圖337， 選擇Main MenuProb DesignProb DefinitnsRandom Input進入變異性分析，如圖338所示，單擊add，添加變異系數，以E為例，在define a random variable選項卡中選中定義好的變量E，下方選擇gauss高斯分布，單擊OK，彈出Quantify Gaussian distribution如圖339，在Mean ，standard *=，單擊ok退出，同樣方法依次添加ceyalixishu，k，houdu1—houdu7的均值和標準差。close all。as=。ES=。由于該隧道為變截面隧道，所以其厚度有變異性，將其列出，142節(jié)點對應的厚度。h0=h01(i)。圖42 matlab命令流選擇在左側選中“判斷大小偏心”，如圖42打開該命令后，選中所有命令，按F9運行，得到結果如圖43g0即為大偏心破壞截面，其余為小偏心，將小于0的截面按照大偏心破壞計算可靠度指標β。syms b f h0 N fy1 As1 ES fy As as as1 M h MUX sigMAXb=1。ES=。7 個變量分別為，1軸力；2fy；3fy1；4厚度h0；5混凝土抗壓強度f；6彎矩；7，h。FinData=zeros(1,42)。normX=eps。gX=[x(4,1)(x(1,1)x(2,1)*As1+x(3,1)*As)/b/x(5,1)1/501/2*x(7,1)+as。1/2*(x(1,1)x(2,1)*As1+x(3,1)*As)^2/b/x(5,1)^2。alphaX=gs/norm(gs)。close all。ES=。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 ]。As=As12(i)。normX=eps。gX=[x(4,1)*x(5,1)/(x(5,1)+5/4*(33/10000*ES+x(3,1))*As/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(7)+as。...... 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。alphaX=gs/norm(gs)。命令流如下%%%%%%%%%%%%%%%%%%%%%%%%%%% clear。f=20e6。fy1=。h01=[ ]。h=h01(i)。 抗壓破壞由公式（45）為功能函數。 %%%%%%%%%%%%存放bbeta數據的矩陣%%%%%end%%%%%%%%%將該命令流文件從matlab中打開，修改方法同鋼筋混凝土，此處函數為，只有3個變量。for i=1:42N=N1(i)。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]。as1=0。clc。*(muXx))/norm(gs)FinData(1,i)=bbeta。1/2*x(1,1)]。...... x(4,1)*x(5,1)*As/(x(5,1)+5/4*(33/10000*ES+x(3,1))*As/x(4,1))9/4*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))^2*As+5/4*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*As/x(4,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)+x(7,1)/2as)。sigmaX=sigMAX(:,i)。FinData=zeros(1,31)。As11=[ ]。syms b f h0 N fy1 As1 ES fy As as as1 M h MUX sigMAXb=1。*(muXx))/norm(gs)FinData(1,i)=bbeta。1/2*x(1,1)]。As*x(4,1)(x(1,1)x(2,1)*As1+x(3,1)*As)/b/x(5,1)*As。g=(x(1,1)x(2,1)*As1+x(3,1)*As)*x(4,1)(x(1,1)x(2,1)*As1+x(3,1)*As)^2/2/b/x(5,1)+x(2,1)*As1*(x(4,1)as1)x(1,1)*(x(6,1)/x(1,1)+x(7,1)/2as)。sigmaX=sigMAX(:,i)。 SigMAX為標準差，為均值乘以變異系數所得，順序同上。As1=。ES=。close all。FinData(1,i)=g。42個單元的厚度FinData=zeros(1,42)。As1=。fy=3e8。清空對話框syms b f h0 N fy1 As1 ES fy As as as1 M h 變量符號，此處為不變常量b=1。如圖342，單擊ok開始循環(huán)。圖331 List Element Table Data對話框圖332 單元表數據第九步 變異性分析圖333 get函數（1） 提取軸力彎矩結果選擇菜單parametersGet Scalar Data,彈出Get Scalar Data對話框如圖333 ，選擇Results data里的element results選項，單擊ok進入get Element Results Data選項卡，在Name of parameter to be difined填入n1，Element number N 填1，Results data to be retrieved選擇By sequence num，后面填SMISC，1提取單元1的軸力，如圖334所示：圖334 get Element Results Data由于本模型有82個單元，所以以相同方式提取82次。6）畫結構軸力圖圖330結構軸力圖（單位：N）選擇菜單Main Menu：General Postproc→Plot Results→Contour Plot→Line Element Results命令，打開一個Plot LineElement Results對話框。圖327定義好后的單元數據表對話框3）設置彎矩分布標題選擇菜單Utility Menu：File→Change Title命令，出現Change Title對話框，在輸入欄中輸入BENDING MOMENT distribution，單擊OK按鈕關閉該對話框。2）將節(jié)點彎矩、剪力、軸力制表圖326定義單元數據表對話框選擇菜單Main Menu：General Postproc→Element Table→Define Table命令，打開一個Element Table Data對話框，如圖325所示。經過2次計算，得到沒有受拉彈簧時的隧道結構模型。2）刪除受拉彈簧單元選擇菜單Main Menu：Preprocessor→Modeling→Delete→Elements命令，打開一個刪除單元選取對話框，在輸入欄中輸入11111111111111112122