#### Abstract

This paper comprehensively studies the effects of linear and nonlinear creep and shrinkage on shear strength and deformation of reinforced concrete (RC) deep beams through a time-dependent reliability analysis framework. The three-dimensional finite element-based program of ABAQUS was used to model RC deep beams. A viscoelastic approach was adopted to model linear creep and shrinkage by a user-defined material model developed and implemented in user subroutines UMAT and UEXPAN. The software was initially examined using two experimental results for short/long-term behavior of shallow and deep concrete beams. In the nonlinear range, creep was taken into account by the affinity hypothesis theorem to consider the effect of high-level sustained loading. Due to the complicated and time-consuming nature of the finite element method (FEM), adaptive neuro-fuzzy inference system (ANFIS) and Monte Carlo simulation (MCS) replaced these complex analyses. Finally, based on the numerical results obtained from the analyses of case study samples, it was concluded that, in a serviceability limit state, for RC deep beams, the probability of failure was reduced to less than one-fifth that for the shallow beams. On the contrary, in a strength limit state, a safety factor of about 1.7∼1.8 could be considered for the effect of sustained high-level loading on the shear strength of RC deep beams.

#### 1. Introduction

The role of RC deep beams in bearing huge gravitational loads (such as above the entrances to high-rise RC buildings with tubular structural systems) makes these members a vital structural element. Therefore, they have been the subject of extensive research studies to understand their performance within the linear and nonlinear ranges of behavior [1]. However, although a couple of studies have addressed the effects of linear/nonlinear creep and shrinkage on shallow beams and columns [2, 3], so far, these effects on the behavior of concrete deep beams have remained understudied [4]. Creep and shrinkage are affected by many factors, including concrete mix [5, 6], loading level [7, 8], concrete strength [9], loading time, geometric properties, construction method, and environmental conditions [10]. Because of the significant violating effect of the spreading strain on cross sections of shallow beams in plane fashion, prediction of the effects of creep and shrinkage on time-dependent reliability of deep concrete beam seems noteworthy to all professional designers [4]. In the present study, the scope of the reliability analysis is extended to linear and nonlinear creep and shrinkage in RC deep beams. Many researchers have conducted time-dependent reliability analysis of concrete structures under induced defects [11], corrosion [12], or their combination [13, 14]. Consideration of creep and shrinkage in this type of analysis is a more recent trend in the literature. Stewart [15] developed a probabilistic model to estimate reliability of shallow concrete beams in the serviceability limit state. Guo et al. [16] presented a time-dependent reliability assessment for PSC box-girder bridge including creep, shrinkage, and corrosion. Ma and Wang [17] implemented a time-dependent reliability for CFST arch bridge (Yajisha Bridge) considering creep, shrinkage, and corrosion.

Previous research has mainly focused on the study of linear/nonlinear creep and shrinkage in shallow RC beams, while their effects on deep beams are absent in the literature. Given this limitation, the present research focuses on three major developments as follows:(1)In case of RC deep beams, in addition to the serviceability limit state, the strength limit state is also considered(2)Effect of nonlinear creep is captured in both serviceability and strength limit states(3)Time-dependent reliability analysis is extended to RC deep beams for the probabilistic nature of deterioration over time

The objective of this study is developing a comprehensive three-dimensional finite element model of RC deep beams in order to carry out time-dependent reliability analysis considering both linear and nonlinear creep and shrinkage. To build the prediction model, artificial neural networks (ANNs) have been developed to solve complex and more elaborated problems [18]. Unfortunately, the data required to construct ANN cannot be gathered directly from experimental results. Indeed, there are no data available in the literature obtained from the results of experimental evidence pertaining to the behavior of concrete deep beams subjected to linear/nonlinear creep and shrinkage effects. Thus, train, test, and checking data should be obtained indirectly using a verified finite element method (FEM). Because of modeling and computational costs, multiple applications of FEM are not practical and an artificial intelligent-based modeling technique such as neuro-fuzzy inference system (ANFIS) may be used instead of these time-consuming and complex analyses. So far, researchers have successfully used ANN and ANFIS in civil engineering problems [19, 20]. A few research studies can be found in the literature which have focused on applying ANN and ANFIS to the prediction of strength and displacement of concrete structures including creep and shrinkage. Bal and Buyle-Bodin [21] applied a ANN to estimate the dimensional changes due to creep in 186 tests extracted from the RILEM (International Union of Laboratories and Experts in Construction Materials). Bilir et al. [22] used ANFIS for predicting the cracks widths of a specific mortar considering shrinkage as an important parameter in concrete behavior. Pendharkar et al. [23] presented a closed-form equation for estimating long-term deflection of high-rise concrete composite frames using the neural network. Al-Zewainy et al. [24] provided a statistical approach for predicting the long-term deflection of RC beams using ANN. Although both tools have widely been utilized with a satisfactory performance, researches have shown that, for many applications, ANFIS exhibits a greater performance than ANN and regression-based models do [25, 26]. Thus, in this work, the ANFIS technique is used to predict the maximum deflection and shear strength of RC deep beams in relation to a twofold performance level of serviceability and ultimate shear strength limit states. The Monte Carlo simulation method is utilized to attain the probability of failure for each limit state condition.

#### 2. Analytical Modeling of Creep and Shrinkage of Concrete

##### 2.1. Linear and Nonlinear Concrete Creep and Shrinkage in Time-Dependent Deflection

Creep and shrinkage are expressed in terms of strain. At time , the total strain of an uncracked concrete specimen under uniaxial loading has several components including elastic strain , creep strain , and shrinkage strain . It is usually assumed (although not completely true) that the strains are independent of each other, so they can be calculated separately and then summed up together. Therefore, the concrete strain at a given point and time can be written as follows [2]:where is the instantaneous strain, represents the creep strain, and shows the drying shrinkage strain. Creep strain is classified into two categories: linear and nonlinear ranges. At service level stresses in structures, which are less than about , and in the absence of cracking, the creep strain depends on stress linearly; however, when the stress exceeds this level, nonlinear creep happens.

Previous studies have claimed that linear creep does not considerably affect shear and flexural strengths of shallow concrete beams [27]. As stated above, shallow and deep beams are different regarding strain distribution in the depths of their cross sections. Furthermore, the effects of strain distribution should be examined on the long-term behavior of RC deep beams. In the present study, the long-term behavior of deep beams in the serviceability limit state is first addressed using the accurate FEM approach through user subroutines UMAT and UEXPAN added to ABAQUS finite element software. Subroutine UMAT is associated with simulating the behavior of concrete under time-dependent factors (creep). Subroutine UEXPAN is a behavioral model of material under volume expansion or contraction used for simulation of shrinkage in concrete.

In linear creep theory, creep strain can be defined as follows [2]:where is the creep coefficient, which is the ratio of the creep strain to the elastic strain at the start of loading.

When a concrete specimen is under a permanent stress exceeding , nonlinear creep occurs [2]. Previous studies have indicated that the viscoelastic approach cannot be used for predicting the concrete behavior at high levels of stress, and it is necessary to make some corrections. Affinity hypothesis can be applied for generalizing linear creep to nonlinear creep from the serviceability point of view; it is used within the stress range of to [28]. In this hypothesis, the linear creep factor can be multiplied by a coefficient to consider the effects of nonlinear creep on displacements of beams. This hypothesis assumes that the linear and nonlinear creep strains are related through the actual stress ratio of within the stress range of to . The proportional coefficient can be estimated by the following equation [28]:

Previous research has suggested that the coefficient could be properly utilized for nonlinear creep problems [28]. The proportional coefficient proposed by Bažant and Kim [29] can be applied in the stress range exceeding , as defined in the following equation:where and . Parameter models the failure effects of concrete at high-stress levels. In this study, considering the effects of nonlinear creep, displacements of deep beams can be obtained using the proportional coefficient and outcome of subroutine UMAT for linear creep.

##### 2.2. Rate-of-Creep Method and Volume Expansion to Consider Linear/Nonlinear Creep and Shrinkage

Under constant environmental factors, stress is expressed by a compliance function defined as the strain at time caused by a unit uniaxial stress applied at time . The relationship between stress and strain in concrete can be expressed using the Kelvin chain model. The model includes a number of Kelvin units called . Each Kelvin unit includes a spring with a hardness of parallel to a damper with the viscosity of . Note that parameter is known as retardation time [30]. In the step by step analysis, the values of and can be assumed constant in each time interval whose proportional magnitude would be adopted in the next step. The outline of the computer program consisting of implementation of rate-of-creep (UMAT) and shrinkage method (UEXPAN) in ABAQUS is provided in Algorithm 1 [30].

| ||||||||||||

#### 3. Uncertainty Modeling of Shrinkage and Creep

Different material models can be considered in analysis for estimating creep and shrinkage effects on concrete structures including ACI209-08 [31], CEB-FIP90 [32], the Japan Society of Civil Engineers (JSCE) [33], the Gardner and Lockman (GL2000) [34], and B3 [35]. The most frequently used model in codes is the one suggested by ACI Committee 209. In this code, the shrinkage strain may be calculated using the following equation [31]:where is the shrinkage strain at time (age of concrete); denotes the ultimate shrinkage strain (at the time of infinity) which depends on experimental results, while coefficient depends on the curing condition; is the age of concrete at start of drying; shows the uncertainty factor of the shrinkage model. Indeed, creep coefficient may be written as follows [31]:where is the creep coefficient at time (age of concrete), represents the ultimate creep coefficient (at the time of infinity) which depends on experimental results, is the observation (current) time (days), and is the uncertainty factor of the creep model. For standard conditions and at ambient relative humidity of 40%, the average value suggested by the ACI-209 standard for the ultimate shrinkage strain is . In the absence of specific data, the average value recommended for is calculated using the following equation [31]:where is an empirical factor.

In this study, coefficients of variations of shrinkage and creep properties are set to 54.2% and 51.7%, respectively, based on the previous studies [35, 36].

#### 4. Deep Beams’ Description and FE Modeling

ABAQUS suggests several models to simulate the nonlinear behavior of concrete. The concrete damaged plasticity (CDP) model is different from others in that it is able to separately simulate compressive and tensile behaviors of concrete and introduce reduction of elastic modulus. In addition to static loads, the model is also capable of modeling cyclic loads. Note that the behavioral model uses Drucker–Prager failure criteria with Lubliner et al.’s yield function [37]. For the serviceability limit state, subroutines UMAT and UEXPAN are used, while for the strength limit state, the CDP model along with a proposed method, here called the “fixed volume contraction scheme,” type of simulating shrinkage, is employed to evaluate the probability of failure in each target age.

Due to lack of experimental data related to creep and shrinkage of deep concrete beams, for gathering train and test data, required for time-dependent reliability analysis, the finite element simulation method was used. In the beginning, the software was initially examined using the experimental results for long-term deflection of shallow concrete beams by Gilbert [38]. Then, the results of the software were compared with the final strength of deep concrete beams determined in the laboratory work by Salamy et al. [39]. According to the two-step verifications, it was concluded that the software could manage to properly model deep concrete beams and capture the effects of creep and shrinkage. In the Gilbert’s test [38], 6 conventional concrete beams with different geometries and loading conditions underwent long-term loading. All specimens were simply supported beams with a span of 3.5 m under a constant load for 400 days. The applied force varied from 11.8 kN to 27 kN. The analytical results obtained by implementing subroutines UMAT and UEXPAN were compared with the test results. The geometrical characteristics of the tested beams are presented in Table 1 and Figure 1.

As observed in Figure 2, the proposed finite element analytical model shows good accuracy concerning changes in reinforcement and loading on beams. The model does not indicate a notable deviation in the changes of parameters; this properly matches experimental results. As the creep model of ACI does not consider drying creep, analytical results show 1% to 6% deviation from laboratory results after 400 days of loading.

**(a)**

**(b)**

**(c)**

In the second stage of verification, the experimental work by Salamy et al. [39] is addressed. Six simply supported RC deep beams are selected and modeled with differences in loading, physical, and geometrical properties. The geometrical properties of these 6 beams are provided in Figure 3 and Table 2.

In the following, a bullet-point summary for RC deep beams modeling is shown:(i)Geometry modeling using ABAQUS(ii)Material modeling using the concrete damaged plasticity (CDP) model or UMAT and UEXPAN subroutines(iii)Analysis of RC deep beams as well as the postprocessing procedure(iv)Verification of results using experimental data

For concrete, the C3D8R element, which is a linear brick element for general purposes, is used, and for rebar, a 3D beam element, namely, B31, is used [40]. As can be seen in Table 3, there is an excellent consistency between the analytical and experimental results.

As a conclusion, the two-stage verifications can demonstrate capability of the analytical modelling to simulate time-effect deformation and predict strength of the RC deep beams.

#### 5. Definition of Input Variables

The properties of the input variables assumed in this work are summarized as follows. Shrinkage and creep model uncertainty factors (, ) suggested by Bazant and Baweja [35], compressive strength of concrete () according to Yang [36], compression reinforcement ratio (), shear span-to-effective depth ratio (), age of concrete in the first loading (), and relative humidity () are assumed to be random variables following a normal distribution [36]. It is also assumed that these input variables are independent of each other. Note that sensitivity analysis is conducted based on the two parameters of shear span-to-effective depth ratio () and sustained load . Creep and shrinkage parameters such as , , , , , and are randomly altered for each analysis based on Gaussian distribution fuzzy membership functions. The depth of beam changes incrementally from 0.25 m to 2.0 m with the same and sustained load . The change of is applied with a rise in with a constant value of *a*. For each time of dependent reliability analysis, three load levels of 40%, 70%, and 90% (for the strength limit state only) are chosen for the shear/flexural strength. Table 4 summarizes the boundary range of input random variables for ANFIS modeling.

#### 6. Definition of the Limit State

A structure must be designed to simultaneously satisfy both serviceability and strength requirements. Two limit states for serviceability and shear strength are considered in this study. The serviceability limit state is usually dealing with excessive static deflections that can cause occurrence of visible cracks in concrete beams. In the serviceability limit state, maximum deflections of deep beams are compared to allowable beam deflection limits. Accordingly, the serviceability limit state function can be written as follows:where is the allowable deflection limit and represents the maximum midspan deflection of deep concrete beams.

The deflection limit in equation (8), which takes into account creep and shrinkage effects, is usually unclear. Unlike for shallow beams, in design codes, there is no explicitly prescribed deflection limit for the effect of long-term loading on RC deep beams. Thus, midspan deflection limits for shallow beams such as and can be considered along with two stricter limitations of and for sensitive structures, where is the beam span [17, 41].

In the ultimate strength reliability assessment, the probability of shear failure of deep concrete beams (in the case of shallow beams, flexural failure can be conceived) in their lifetime is addressed. In order to assess failure probability, it is necessary to know the margin between two random variables of resistances and applied load . Thus, the ultimate limit state function can be expressed as follows:where in general, and are considered as time-independent output variables.

The sign function describes the state of structure. For positive value of , the state of structure is safe; defines the domain of failure events; and conventionally defines the limit state surface, i.e., the boundary between the two mentioned domains.

#### 7. Architecture of ANFIS Modeling

ANFIS (adaptive neuro-fuzzy inference system) is a well-known artificial neural network based on a fuzzy concept for modeling complex systems. ANFIS represents Sugeno–Tsukamoto fuzzy models in which a hybrid learning algorithm is used. ANFIS was first proposed by Jang [42]. ANFIS is a type of artificial neural network that combines the advantages of both fuzzy inference and ANN (artificial neural network) methods. Advantages of ANN can be categorized as follows [43]:(i)Using a self-learning process for finding a reliable solution(ii)A simplified algorithm without any complex mathematical formulation(iii)Ability to solve moderate complex nonlinear problems(iv)Robustness of algorithm in processing structural analysis

Despite its merit, some disadvantages of this method can be highlighted as follows:(i)Training process significantly impacts the computational cost(ii)Effect of the ambiguity of information is not considered in this method

Some advantages of ANFIS compared to ANN can be introduced as follows [44]:(i)Combination of powerful mathematical concept of the artificial neural network and human knowledge-based fuzzy logic system. This provides a robust framework for solving complicated problems.(ii)The ability to simulate fully nonlinear complex systems.(iii)Rapid learning and reducing computational costs.

Therefore, according to the mentioned points, ANFIS is used for predicting long-term deflection of RC deep beams.

In the following, a bullet-form ANFIS architecture designed for 7 input variables (depicted in Table 4) using the Gaussian membership function is shown:(i)Selecting the proper input data(ii)Fuzzification step: assigning the membership function(iii)Implication of rules for the establishing the fuzzy network(iv)Normalization and weighted outputs(v)Defuzzification as well as aggregation(vi)Output layer including midspan deflection or shear strength of RC deep beams

In ANFIS, a hybrid method (a hybrid algorithm adjusts the resultant parameters in a forward pass and the premise parameters in a backward pass) and gradient descent scheme along with the least-squares method aim to minimize the data error.

#### 8. Verification of ANFIS

The ANFIS model developed in this research is first trained by 100 input-output datasets and further tested by 25 datasets. Due to the large number of input parameters and despite the satisfactory results in the first stage, for final training of the network, 280 input-output datasets are used. The 100 trained data for estimating the deflection of RC deep beams are reported in Table 5.

To enhance the accuracy of the ANFIS model, a wide range of input parameters are selected and applied to the model. These datasets are created through FEM analyses, which have been verified by experimental results. After testing different learning algorithms with 480 different epochs, the model reaches the minimum error size (Figure 4). On the contrary, the maximum convergence in the model is reached after 1800 epochs. To evaluate the adequacy of the model developed in ANFIS, three norms are used. These norms are root-mean-square error (RMSE), mean squared error (MSE), and mean absolute deviation (MAD) [45, 46].

Figure 4 depicts the training process of the developed ANFIS model to establish the fuzzy relationship between input and output variables, including the maximum deflection of RC deep beams. ANFIS model minimizes the error by increasing epoch’s number to stabilize the network.

In addition, Table 6 lists the results of sensitivity analysis of ANFIS at the serviceability level for the training step based on root-mean-square error (RMSE), mean squared error (MSE), and mean absolute deviation (MAD).

Figure 5 demonstrates a good consistency between target and output data, confirming the efficiency of the ANFIS model.

**(a)**

**(b)**

In the case of the strength limit state, the performance of the model is evaluated using two different error functions. The root-mean-square error (RMSE) and mean absolute percentage error (MAPE) [46], for the strength limit state of RC deep beams, are represented in Table 7.

As presented in Figure 5 and Tables 6 and 7, it can be concluded that the ANFIS model is capable enough to predict deflection and shear strength of RC deep beams. Comparison of figures and errors suggests that subclustering the ANFIS model along with the hybrid optimization approach has been close to FEM and a few experimental results. Figure 6 shows the sample results of training data in time-dependent ANFIS models for minimizing the error size by a proper number of epochs to stabilize sufficient ANFIS models.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

#### 9. Monte Carlo Simulation

When failure surface cannot be expressed analytically, the Monte Carlo Simulation (MCS) method is useful in reliability analysis. Problems subjected to complex geometry or complex phenomena such as creep and shrinkage with a large number of basic variables need an efficient approach including MCS. Unlike analytical methods such as first-order reliability method (FORM) and second-order reliability method (SORM), which suffer correlations of variables and give approximate solutions to most practical problems with high numbers of input variables, MSC is able to yield more reliable solutions with adequate realizations [47].

Note that the ANFIS model is developed and MCS is adopted for each target age. Therefore, to estimate the probability of failure of concrete beams for each point, a massive calculation is required. In the following, a bullet-point summary for time-dependent reliability analysis of RC beams is shown.(i)Finite element modeling (FEM)(ii)Train and test data gathering using FEM results(iii)ANFIS network setup and gathering of all datasets(iv)Monte Carlo Simulation (MCS) and probability of failure.

##### 9.1. Serviceability Limit States

By performing structural reliability analysis at a certain age with four different allowable values of the midspan deflection of a beam with *L* = 3.0 m, the probability of RC deep beam failure at the serviceability limit state is estimated. In analyses, linear creep, nonlinear creep (which is dependent on the level of loading), and shrinkage effects in the serviceability limit state are addressed. For deflection limits of , , , and , the probability of failure in different ages of construction is assessed [17, 41]. As noted earlier, the probability of failure is first calculated for a constant length by changing from 0.5 to 4. The corresponding curves are presented in Figure 7.

**(a)**

**(b)**

**(c)**

**(d)**

As observed in Figure 7, creep and shrinkage have less significant effects on deflection of RC deep beams. The shear behavior seems to play a dominant role in time-dependent deflection of RC deep beams. For beams with and plus exceedance probability for limitation at load level of , the beam strength is equal to 0.0332 and 0.183, respectively. A fivefold increase can be observed for the shallow beam. Note that failure probability of 0.183 has occurred at limitation, which is a rigorous limitation for beams. A sudden drop of failure probability in case of and is also observed. This decline can be interpreted as a margin between the two types of RC beams (i.e., deep and shallow). Shear-dominated behavior cannot occur for values greater than 2, while flexural-dominated behavior can be seen for the ratios of and . For the load level of of the beam strength ratio, it is observed that nonlinear creep significantly affects the final deflection of RC beams. These results confirm that, in the nonlinear creep range, the affinity hypothesis theorem must be used for correct estimation of the long-term deflection of RC beams. As demonstrated in Figure 7(b), after 3 years, about 85% of the final failure probability is developed. After several years, under sustained load, the rate of change in is very small. Design engineers and insurance companies can take this conclusion into account in their design and insurance premium calculation.

##### 9.2. Strength Limit States

At a strength limit state level, the concrete damaged plasticity model is used to model RC deep beams where a constant volume reduction scheme, a kind of simulating shrinkage, is employed to model shrinkage and to adapt to the concrete model. Verification of this method was carried out by Hyodo et al. [48]. Three RC beams with effective depths of 250, 500, and 1000 mm have been selected and exposed to drying shrinkage. The geometry and properties of the beams are outlined in Tables 8 and 9. *b* is the beam width; *h* is the beam height; *L* is the beam length; *c* represents the distance between two loading points; *a* shows the shear span; *d* denotes the effective depth; *f*_{c} refers to the compressive strength; *f*_{t} reflects the splitting tensile strength; *E*_{c} is Young’s modulus.

The time-dependent shear strengths of the beams calculated based on the experimental test and the proposed model in this paper are shown in Table 10. For all beams, the proposed model predicts the shear strengths of RC beams accurately. At this step, the verification of software is completed.

In the first stage of the analyses, linear creep and shrinkage are applied to the FEM model to obtain test and train data for developing ANFIS. As mentioned earlier, the collapse probability of beams is first calculated with a constant length and altering from 0.5 to 4. Figure 8 reveals the relationship between collapse probability (shear or flexural failure) and type of the RC beam. It is clear that shrinkage and linear creep within the low-stress level range do not have any major effect on the collapse probability of RC beams.

As shown in Figure 8, over time, the collapse probability will grow mildly. Previous research investigations have also demonstrated that linear creep does not considerably affect the shear and flexural strengths of shallow beams, confirming the results shown in Figure 8 [27]. For both types of beams, low values of collapse probability can be achieved, as observed in this figure. For instance, with linear creep and shrinkage assumption, for , a maximum collapse probability of 0.0114 is obtained, which is very small and can be neglected by design engineers and insurance companies. Note that no tangible difference can be observed between the rates of collapse over time for RC deep and shallow beams.

In the second stage, the collapse probability for RC beams at the load level of 70% of the beam strength is investigated in the presence of creep and shrinkage effects. Figure 9 indicates the results for beams with different shear span-to-depth ratios subjected to the load level of 70% of the beam strength. It is noteworthy that, within this range of stress, the linear creep should be substituted by a nonlinear creep. Furthermore, a damage index considering stiffness degradation along with the creep damage model is employed here.

To this end, the following nonlinear damage index (equation (10)) recommended by Mazzotti and Savoia [49] is used (Figure 10):

The computed collapse probability for both deep and shallow beams is revealed in Figure 10. As can be seen, the effect of shrinkage and nonlinear creep on the shear strength of RC deep beams is considerable. In the case of shallow beams, the shrinkage and nonlinear creep do not have prominent effect on flexural strength of shallow beams. If a volume reduction scheme without the nonlinear damage index is applied, 75% of the total collapse probability is achieved. It can be deduced that, at moderate sustained load levels, the greatest share of collapse probability belongs to the shrinkage. For instance, in the case of deep beams is up to 5 times larger than that in shallow beams after 1000 days. As a conclusion, within the load level range of 70% of the beam strength, a safety factor of 1.25 is applicable to the shear design of RC deep beams considering shrinkage and creep effects. In other words, a coefficient of 0.75 can be considered in calculating the shear strength of RC deep beams in order to capture the effects of shrinkage and nonlinear creep in the design. For shallow beams, is equivalent to 0.05, confirming the ineffectiveness of shrinkage and creep on their flexural strength.

In the ultimate load step, a constant load of 90% of the beam strength is applied to test the effect of shrinkage and creep at a high-stress level. The results are shown in Figure 11.

These results confirm that shrinkage and nonlinear creep play an important role in shear strength of RC deep beams. Shrinkage along with creep accelerates the failure of deep beams. According to Figure 11, the collapse probability increases when *a*/*d* decreases. It seems that, in high-level sustained loading, the role of nonlinear creep can hardly be distinguished from that of shrinkage, as both of them cause collapse in a very high-level loading. In shallow beams, shrinkage has no effect on flexural strength, and all beams are damaged by nonlinear creep. Based on the results shown in Figure 11, a safety factor between 1.7 and 1.8 can be considered for preventing shear failure in RC deep beams subjected to the high loading level. On the contrary, a strength reduction coefficient of 0.2 can be used for calculating the shear strength of RC deep beams. Sudden drops because of transition from shallow to deep beams are again notable. It is because of the difference between the two types of beams. In beams with *a*/*d* = 0.5–2, the principal failure is shear-dominant, while in beams with *a/d* values greater than 2, flexural failure may occur. In the shear-dominated case, shrinkage and nonlinear creep play an important role.

#### 10. Conclusion

This paper presented a time-dependent reliability analysis of deep concrete beams with a special attention to linear/nonlinear creep and shrinkage. Because of the lack of experimental results for RC deep beams considering creep and shrinkage, FEM analysis was carried out by ABAQUS software with ANFIS employed to obtain the data required for time-dependent reliability analysis. Monte Carlo simulation was used to determine the probabilities of failure. The results indicated that ANFIS was a powerful tool in the prediction of deflection as well as shear and flexural strengths of RC beams. In addition, the following conclusions could be drawn:(1)It was found that shear-dominated behavior of the deep beam played an important role in the failure probability of RC beams in the serviceability limit state. Reduction in *a*/*d* resulted in a higher safety margin in the final long-term deflection. Also, more than fivefold increase was observed in the failure probability for shallow beams. Thus, RC beams with a higher *a*/*d* ratio in the presence of creep and shrinkage exhibited far lower reliability levels than deep beams in the serviceability limit state. Generally, within the linear creep range, exceedance probability of the beam deflection limit was as low as about 0.18. The results of the analyses indicated that, by considering the shrinkage effect, a safety factor within 1.2∼1.3 could be applied. For a practical calculation, a reduction coefficient between 0.8 and 0.7 could be considered for the shear strength of RC deep beams.(2)In the strength limit state, shrinkage and nonlinear creep effects led to greater collapse probability, especially in RC deep beams. In comparison with shallow beams, the effect of shrinkage was more pronounced (shrinkage increased the probability of failure by more than 4 times). One reason could be the difference in failure patterns of shallow and deep beams. When the concrete beam was affected by shrinkage, development of cracks in full depth could cause reduction in the shear strength. Thus, the collapse probability rose in the presence of shrinkage. The analysis results indicated that a safety factor of 1.7∼1.8 was required to prevent the risk of shear failure in a high-level sustained loading of RC deep beams. For medium and low load levels, it was not necessary to consider an additional safety factor to prevent flexural or shear failure.

Based on the current study, the following prospective for future researches can be suggested:(1)The time-dependent reliability analysis has been verified using two-stage verifications, but the specific experimental evidence on long-term behavior of RC deep beams can be carried out to enhance this investigation(2)More investigations shall be carried out to study the effect of other parameters such as temperature on time-dependent reliability of concrete structures(3)Proposing a new formulation for long-term deflection of RC deep beams may be a potential alternative for future exploration

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.