*********************************************************************** * * RSTHF.dgibi * * Random Set Theory analysis with Finite Elements (RS-FEM) * *********************************************************************** 'OPTI' 'DIME' 2 'MODE' 'AXIS' 'ELEM' 'QUA4' ; 'OPTI' 'TRAC' 'PSC'; GRAPH = VRAI ; R = 5.00 ; RM = 25.00 ; HSA = 5.00 ; HSI = 10.00 ; HCL = 15.00 ; P1 ='POIN' 0.0 0.0 ; P2 ='POIN' 0.0 (-1.0 * HSA) ; P3 ='POIN' 0.0 (-1.0 * (HSA + HSI)) ; P4 ='POIN' 0.0 (-1.0 * (HSA + HSI + HCL)) ; P5 ='POIN' R 0.0 ; RH = 0.50 ; RV = 0.50 ; NVSA = ('ENTI' (HSA / RV) ) ; NVSI = ('ENTI' (HSI / RV) ) ; NVCL = ('ENTI' (HCL / RV) ) ; D1 = DROI NVSA P1 P2 ; D2 = DROI NVSI P2 P3 ; D3 = DROI NVCL P3 P4 ; NR = ( 'ENTI' (R / RH) ) ; NRM = ( 'ENTI' ((RM - R) / RH) ) ; NM = ( 'ENTI' (RM / RH) ) ; 'MESS' 'NR=' NR ; 'MESS' 'NRM=' NRM ; 'MESS' 'NM=' NM ; SRSA = D1 'TRAN' NR ('POIN' R 0.0) 'COUL' 'ROUG' ; D11 ='COTE' 3 SRSA ; SMSA = D11 'TRAN' NRM ('POIN' (RM - R) 0.0) 'COUL' 'ROUG' ; SRSI = D2 'TRAN' NR ('POIN' R 0.0) 'COUL' 'BLEU' ; D21 ='COTE' 3 SRSI ; SMSI = D21 'TRAN' NRM ('POIN' (RM - R) 0.0) 'COUL' 'BLEU' ; SRCL = D3 'TRAN' NR ('POIN' R 0.0) 'COUL' 'VERT' ; D31 ='COTE' 3 SRCL ; SMCL = D31 'TRAN' NRM ('POIN' (RM - R) 0.0) 'COUL' 'VERT' ; DP ='COTE' 4 SRSA 'COUL' 'NOIR' ; ST = (SRSA 'ET' SMSA 'ET' SRSI 'ET' SMSI 'ET' SRCL 'ET' SMCL) ; ELIM 0.1 ST ; CST = CONT ST COUL NOIR ; D4 ='COTE' 2 SMCL ; D5 ='COTE' 3 SMCL ; D6 ='COTE' 3 SMSI ; D7 ='COTE' 3 SMSA ; MOCO = 'MODE' DP MECANIQUE ELASTIQUE ISOTROPE COQ2 ; MOSA = 'MODE' (SRSA 'ET' SMSA) MECANIQUE ELASTIQUE ISOTROPE ; MOSI = 'MODE' (SRSI 'ET' SMSI) MECANIQUE ELASTIQUE ISOTROPE ; MOCL = 'MODE' (SRCL 'ET' SMCL) MECANIQUE ELASTIQUE ISOTROPE ; MOTOT= ( MOCO 'ET' MOSA 'ET' MOSI 'ET' MOCL ) ; BLR1 ='BLOQ' UR (D1 'ET' D2 'ET' D3) ; BLR2 ='BLOQ' UR (D5 'ET' D6 'ET' D7) ; BLZ1 ='BLOQ' UR UZ D4 ; BLT = BLR1 'ET' BLR2 'ET' BLZ1 ; CACO ='CARA' MOCO 'EPAI' 1.0 ; MACO ='MATE' MOCO 'YOUN' 3.0e9 'NU' 0.33 ; RS = OBJET @RSTH ; RST ='TABL' ; RST.'PRS' ='TABL' ; RST.'PRS'.'MIN' ='PROG' 1.0e5 ; RST.'PRS'.'MAX' ='PROG' 2.0e5 ; RST.'PRS'.'CPB' ='PROG' 1.0 ; RST.'PRS'.'MINS' = RS%'SCV' RST.'PRS'.'MIN' ; RST.'PRS'.'MAXS' = RS%'SCV' RST.'PRS'.'MAX' ; RST.'PRS'.'CPBS' = RS%'SCS' RST.'PRS'.'CPB' ; PRSN = ('EVOL' 'MANU' 'MIN' RST.'PRS'.'MINS' 'CPB' RST.'PRS'.'CPBS' ) ; PRSX = ('EVOL' 'MANU' 'MAX' RST.'PRS'.'MAXS' 'CPB' RST.'PRS'.'CPBS' ) ; RST.'ESA' ='TABL' ; RST.'ESA'.'MIN' ='PROG' 3.0e7 3.2e7; RST.'ESA'.'MAX' ='PROG' 3.6e7 3.8e7; RST.'ESA'.'CPB' ='PROG' 0.5 0.5 ; RST.'ESA'.'MINS' = RS%'SCV' RST.'ESA'.'MIN' ; RST.'ESA'.'MAXS' = RS%'SCV' RST.'ESA'.'MAX' ; RST.'ESA'.'CPBS' = RS%'SCS' RST.'ESA'.'CPB' ; ESAN = ('EVOL' 'MANU' 'MIN' RST.'ESA'.'MINS' 'CPB' RST.'ESA'.'CPBS' ) ; ESAX = ('EVOL' 'MANU' 'MAX' RST.'ESA'.'MAXS' 'CPB' RST.'ESA'.'CPBS' ) ; RST.'NSA' ='TABL' ; RST.'NSA'.'MIN' ='PROG' 0.30 ; RST.'NSA'.'MAX' ='PROG' 0.35 ; RST.'NSA'.'CPB' ='PROG' 1.00 ; RST.'NSA'.'MINS' = RS%'SCV' RST.'NSA'.'MIN' ; RST.'NSA'.'MAXS' = RS%'SCV' RST.'NSA'.'MAX' ; RST.'NSA'.'CPBS' = RS%'SCS' RST.'NSA'.'CPB' ; NSAN = ('EVOL' 'MANU' 'MIN' RST.'NSA'.'MINS' 'CPB' RST.'NSA'.'CPBS' ) ; NSAX = ('EVOL' 'MANU' 'MAX' RST.'NSA'.'MAXS' 'CPB' RST.'NSA'.'CPBS' ) ; RST.'ESI' ='TABL' ; RST.'ESI'.'MIN' ='PROG' 2.0e7 2.2e7 ; RST.'ESI'.'MAX' ='PROG' 2.6e7 2.8e7 ; RST.'ESI'.'CPB' ='PROG' 0.3 0.7 ; RST.'ESI'.'MINS' = RS%'SCV' RST.'ESI'.'MIN' ; RST.'ESI'.'MAXS' = RS%'SCV' RST.'ESI'.'MAX' ; RST.'ESI'.'CPBS' = RS%'SCS' RST.'ESI'.'CPB' ; ESIN = ('EVOL' 'MANU' 'MIN' RST.'ESI'.'MINS' 'CPB' RST.'ESI'.'CPBS' ) ; ESIX = ('EVOL' 'MANU' 'MAX' RST.'ESI'.'MAXS' 'CPB' RST.'ESI'.'CPBS' ) ; RST.'NSI' ='TABL' ; RST.'NSI'.'MIN' ='PROG' 0.35 ; RST.'NSI'.'MAX' ='PROG' 0.40 ; RST.'NSI'.'CPB' ='PROG' 1.00 ; RST.'NSI'.'MINS' = RS%'SCV' RST.'NSI'.'MIN' ; RST.'NSI'.'MAXS' = RS%'SCV' RST.'NSI'.'MAX' ; RST.'NSI'.'CPBS' = RS%'SCS' RST.'NSI'.'CPB' ; NSIN = ('EVOL' 'MANU' 'MIN' RST.'NSI'.'MINS' 'CPB' RST.'NSI'.'CPBS' ) ; NSIX = ('EVOL' 'MANU' 'MAX' RST.'NSI'.'MAXS' 'CPB' RST.'NSI'.'CPBS' ) ; RST.'ECL' ='TABL' ; RST.'ECL'.'MIN' ='PROG' 1.0e7 1.2e7 ; RST.'ECL'.'MAX' ='PROG' 1.6e7 1.8e7 ; RST.'ECL'.'CPB' ='PROG' 0.2 0.8 ; RST.'ECL'.'MINS' = RS%'SCV' RST.'ECL'.'MIN' ; RST.'ECL'.'MAXS' = RS%'SCV' RST.'ECL'.'MAX' ; RST.'ECL'.'CPBS' = RS%'SCS' RST.'ECL'.'CPB' ; ECLN = ('EVOL' 'MANU' 'MIN' RST.'ECL'.'MINS' 'CPB' RST.'ECL'.'CPBS' ) ; ECLX = ('EVOL' 'MANU' 'MAX' RST.'ECL'.'MAXS' 'CPB' RST.'ECL'.'CPBS' ) ; RST.'NCL' ='TABL' ; RST.'NCL'.'MIN' ='PROG' 0.40 ; RST.'NCL'.'MAX' ='PROG' 0.45 ; RST.'NCL'.'CPB' ='PROG' 1.00 ; RST.'NCL'.'MINS' = RS%'SCV' RST.'NCL'.'MIN' ; RST.'NCL'.'MAXS' = RS%'SCV' RST.'NCL'.'MAX' ; RST.'NCL'.'CPBS' = RS%'SCS' RST.'NCL'.'CPB' ; NCLN = ('EVOL' 'MANU' 'MIN' RST.'NCL'.'MINS' 'CPB' RST.'NCL'.'CPBS' ) ; NCLX = ('EVOL' 'MANU' 'MAX' RST.'NCL'.'MAXS' 'CPB' RST.'NCL'.'CPBS' ) ; RS%'RST' RST ; JX = RST.'PRS'.'CX' ; 'REPE' J JX ; 'SI' ('EGA' (@MOD &J 100) 0 ) ; 'MESS' 'J=' &J ; 'FINS' ; RS%'RSV' &J 0 ; CHP ='PRES' 'COQU' MOCO RST.'PRS'.'VN' 'NORM' ; MASA ='MATE' MOSA 'YOUN' RST.'ESA'.'VN' 'NU' RST.'NSA'.'VN' ; MASI ='MATE' MOSI 'YOUN' RST.'ESI'.'VN' 'NU' RST.'NSI'.'VN' ; MACL ='MATE' MOCL 'YOUN' RST.'ECL'.'VN' 'NU' RST.'NCL'.'VN' ; MAT = CACO 'ET' MACO 'ET' MASA 'ET' MASI 'ET' MACL ; RG1 ='RIGI' MOTOT MAT ; RGT = RG1 'ET' BLT ; DPT ='RESO' RGT CHP ; DZT ='EXCO' DPT 'UZ' 'UZ' ; DZDP ='REDU' DZT DP ; DZDPN ='MINI' DZDP ; RS%'RSR' DZDPN &J ; 'FIN' J ; VCHP ='VECT' CHP 3.0e-6 'FORC' 'COUL' 'NOIR' ; DFT ='DEFO' DZT ST VCHP ; 'SI' GRAPH ; 'TRAC' 'FACE' (ST 'ET' DP) 'TITR' 'Random-set Theory - Foundation settlement analysis - Mesh' ; 'TRAC' DZT DFT CST 'TITR' 'Vertical displacements and deformed state [m] - Example' ; 'DESS' ( PRSN 'ET' PRSX ) 'TITX' 'Foundation pressure load [N/m^{2}]' 'POSX' 'CENT' 'XBOR' 0.0e5 3.0e5 'XGRA' 2.0e4 'TITY' 'CPB [1]' 'POSY' 'CENT' 'YBOR' 0.0 1.0 'YGRA' 0.1 'TITR' 'Random Set Theory - Foundation settlement analysis' 'GRIL' 'POIN' 'GRIS' ; 'DESS' ( ESAN 'ET' ESAX ) 'TITX' 'Sand - E [N/m^{2}]' 'POSX' 'CENT' 'XBOR' 5.0e6 5.0e7 'XGRA' 5.0e6 'TITY' 'CPB [1]' 'POSY' 'CENT' 'YBOR' 0.0 1.0 'YGRA' 0.1 'TITR' 'Random Set Theory - Foundation settlement analysis' 'GRIL' 'POIN' 'GRIS' ; 'DESS' ( NSAN 'ET' NSAX ) 'TITX' 'Sand - Nu [1]' 'POSX' 'CENT' 'XGRA' 0.05 'XBOR' 0.0 0.5 'TITY' 'CPB [1]' 'POSY' 'CENT' 'YBOR' 0.0 1.0 'YGRA' 0.1 'TITR' 'Random Set Theory - Foundation settlement analysis' 'GRIL' 'POIN' 'GRIS' ; 'DESS' ( ESIN 'ET' ESIX ) 'TITX' 'Silt - E [N/m^{2}]' 'POSX' 'CENT' 'XBOR' 5.0e6 5.0e7 'XGRA' 5.0e6 'TITY' 'CPB [1]' 'POSY' 'CENT' 'YBOR' 0.0 1.0 'YGRA' 0.1 'TITR' 'Random Set Theory - Foundation settlement analysis' 'GRIL' 'POIN' 'GRIS' ; 'DESS' ( NSIN 'ET' NSIX ) 'TITX' 'Silt - Nu [1]' 'POSX' 'CENT' 'XBOR' 0.0 0.5 'XGRA' 0.05 'TITY' 'CPB [1]' 'POSY' 'CENT' 'YBOR' 0.0 1.0 'YGRA' 0.1 'TITR' 'Random Set Theory - Foundation settlement analysis' 'GRIL' 'POIN' 'GRIS' ; 'DESS' ( ECLN 'ET' ECLX ) 'TITX' 'Clay - E [N/m^{2}]' 'POSX' 'CENT' 'XBOR' 5.0e6 5.0e7 'XGRA' 5.0e6 'TITY' 'CPB [1]' 'POSY' 'CENT' 'YBOR' 0.0 1.0 'YGRA' 0.1 'TITR' 'Random Set Theory - Foundation settlement analysis' 'GRIL' 'POIN' 'GRIS' ; 'DESS' ( NCLN 'ET' NCLX ) 'TITX' 'Clay - Nu [1]' 'POSX' 'CENT' 'XBOR' 0.0 0.5 'XGRA' 0.05 'TITY' 'CPB [1]' 'POSY' 'CENT' 'YBOR' 0.0 1.0 'YGRA' 0.1 'TITR' 'Random Set Theory - Foundation settlement analysis' 'GRIL' 'POIN' 'GRIS' ; 'DESS' ( RS.'RT'.'EN' 'ET' RS.'RT'.'EX' ) 'TITX' 'u_{z} [m]' 'POSX' 'CENT' 'XBOR' -4.5e-2 0.0 'XGRA' 0.005 'TITY' 'CPB [1]' 'POSY' 'CENT' 'YBOR' 0.0 1.0 'YGRA' 0.1 'TITR' 'Random Set Theory - Foundation settlement analysis' 'GRIL' 'POIN' 'GRIS' ; 'FINS' ; 'FIN' ;