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