The nonlinear behavior of reinforced fly ash concrete (RFAC) beams until the ultimate failure is highly a multifaceted phenomenon due to the contention of heterogeneous material properties and the cracking behavior of concrete. This paper represents an exploration of nonlinear finite element analysis of reinforced concrete beams with the inclusion of fly ash under flexural loading. Finite element modelling of RFAC beams is carried out using a distinct reinforcement modelling method. The capability of the model to capture the critical crack regions, loads and deflections for various loadings in RFAC beams has been evaluated. Comparison is made between experimental results and finite element modelling with respect to crack formation and the ultimate capacity for flexural loading and mid-span deflection. The achieved results in the present study indicate acceptable approximation with those in the previous investigations.
Citation: Kong Fah Tee, Sayedali Mostofizadeh. Numerical and experimental investigation of concrete with various dosages of fly ash[J]. AIMS Materials Science, 2021, 8(4): 587-607. doi: 10.3934/matersci.2021036
Related Papers:
[1]
Jun Hu, Jie Wu, Mengzhe Wang .
Research on VIKOR group decision making using WOWA operator based on interval Pythagorean triangular fuzzy numbers. AIMS Mathematics, 2023, 8(11): 26237-26259.
doi: 10.3934/math.20231338
[2]
Manar A. Alqudah, Artion Kashuri, Pshtiwan Othman Mohammed, Muhammad Raees, Thabet Abdeljawad, Matloob Anwar, Y. S. Hamed .
On modified convex interval valued functions and related inclusions via the interval valued generalized fractional integrals in extended interval space. AIMS Mathematics, 2021, 6(5): 4638-4663.
doi: 10.3934/math.2021273
[3]
Yanhong Su, Zengtai Gong, Na Qin .
Complex interval-value intuitionistic fuzzy sets: Quaternion number representation, correlation coefficient and applications. AIMS Mathematics, 2024, 9(8): 19943-19966.
doi: 10.3934/math.2024973
[4]
Tatjana Grbić, Slavica Medić, Nataša Duraković, Sandra Buhmiler, Slaviša Dumnić, Janja Jerebic .
Liapounoff type inequality for pseudo-integral of interval-valued function. AIMS Mathematics, 2022, 7(4): 5444-5462.
doi: 10.3934/math.2022302
[5]
Mustafa Ekici .
On an axiomatization of the grey Banzhaf value. AIMS Mathematics, 2023, 8(12): 30405-30418.
doi: 10.3934/math.20231552
[6]
Muhammad Bilal Khan, Muhammad Aslam Noor, Thabet Abdeljawad, Bahaaeldin Abdalla, Ali Althobaiti .
Some fuzzy-interval integral inequalities for harmonically convex fuzzy-interval-valued functions. AIMS Mathematics, 2022, 7(1): 349-370.
doi: 10.3934/math.2022024
[7]
Waqar Afzal, Mujahid Abbas, Jongsuk Ro, Khalil Hadi Hakami, Hamad Zogan .
An analysis of fractional integral calculus and inequalities by means of coordinated center-radius order relations. AIMS Mathematics, 2024, 9(11): 31087-31118.
doi: 10.3934/math.20241499
[8]
Le Fu, Jingxuan Chen, Xuanchen Li, Chunfeng Suo .
Novel information measures considering the closest crisp set on fuzzy multi-attribute decision making. AIMS Mathematics, 2025, 10(2): 2974-2997.
doi: 10.3934/math.2025138
[9]
Li Li, Mengjing Hao .
Interval-valued Pythagorean fuzzy entropy and its application to multi-criterion group decision-making. AIMS Mathematics, 2024, 9(5): 12511-12528.
doi: 10.3934/math.2024612
[10]
Scala Riccardo, Schimperna Giulio .
On the viscous Cahn-Hilliard equation with singular potential and inertial term. AIMS Mathematics, 2016, 1(1): 64-76.
doi: 10.3934/Math.2016.1.64
Abstract
The nonlinear behavior of reinforced fly ash concrete (RFAC) beams until the ultimate failure is highly a multifaceted phenomenon due to the contention of heterogeneous material properties and the cracking behavior of concrete. This paper represents an exploration of nonlinear finite element analysis of reinforced concrete beams with the inclusion of fly ash under flexural loading. Finite element modelling of RFAC beams is carried out using a distinct reinforcement modelling method. The capability of the model to capture the critical crack regions, loads and deflections for various loadings in RFAC beams has been evaluated. Comparison is made between experimental results and finite element modelling with respect to crack formation and the ultimate capacity for flexural loading and mid-span deflection. The achieved results in the present study indicate acceptable approximation with those in the previous investigations.
1.
Introduction
Random numbers are very important in statistical, probability theory, and mathematical analysis in such complex cases, where the real numbers are difficult to record. The random numbers are generated from the uniform distribution when an interval is defined for their selection of random numbers. The random numbers are generated in sequence and depict the behavior of the real data. In addition, random data can be used for estimation and forecasting purposes. According to [1], "The method is based on running the model many times as in random sampling. For each sample, random variates are generated on each input variable; computations are run through the model yielding random outcomes on each output variable. Since each input is random, the outcomes are random. In the same way, they generated thousands of such samples and achieved thousands of outcomes for each output variable. In order to carry out this method, a large stream of random numbers was needed". To generate random numbers, a random generator is applied. The random numbers have no specific pattern and are generated from the chance process. Nowadays, the latest computer can be used to generate random numbers using a well-defined algorithm; see [1]. Bang et al. [2] investigated normality using random-number generating. Schulz et al. [3] presented a pattern-based approach. Tanyer [4] generated random numbers from uniform sampling. Kaya and Tuncer [5] proposed a method to generate biological random numbers. Tanackov et al. [6] presented a method to generate random numbers from the exponential distribution. Jacak et al. [7] presented the methods to generate pseudorandom numbers. More methods can be seen in [8,9,10].
The neutrosophic statistical distributions were found to be more efficient than the distributions under classical statistics. The neutrosophic distributions can be applied to analyze the data that is given in neutrosophic numbers. Sherwani et al. [11] proposed neutrosophic normal distribution. Duan et al. [12] worked on neutrosophic exponential distribution. Aliev et al. [13] generated Z-random numbers from linear programming. Gao and Ralescu [14] studied the convergence of random numbers generated under an uncertain environment. More information on random numbers generators can be seen in [15,16,17,18]. In recent works, Aslam [19] introduced a truncated variable algorithm for generating random variates from the neutrosophic DUS-Weibull distribution. Additionally, in another study [20], novel methods incorporating sine-cosine and convolution techniques were introduced to generate random numbers within the framework of neutrosophy. Albassam et al. [21] showcased probability/cumulative density function plots and elucidated the characteristics of the neutrosophic Weibull distribution as introduced by [22]. The estimation and application of the neutrosophic Weibull distribution was also presented by [21].
In [22], the Weibull distribution was introduced within the realm of neutrosophic statistics, offering a more inclusive perspective compared to its traditional counterpart in classical statistics. [21] further examined the properties of the neutrosophic Weibull distribution introduced by [22]. Despite an extensive review of existing literature, no prior research has been identified regarding the development of algorithms for generating random numbers using both the neutrosophic uniform and Weibull distributions. This paper aims to bridge this gap by presenting innovative random number generators tailored specifically for the neutrosophic uniform distribution and the neutrosophic Weibull distribution. The subsequent sections will provide detailed explanations of the algorithms devised to generate random numbers for these distributions. Additionally, the paper will feature multiple tables showcasing sets of random numbers across various degrees of indeterminacy. Upon thorough analysis, the results reveal a noticeable decline in random numbers as the degree of indeterminacy increases.
2.
Neutrosophic uniform distribution
Let xNU=xNL+xNUIxNU;IxNUϵ[IxLU,IxUU] be a neutrosophic random variable that follows the neutrosophic uniform distribution. Note that the first part xNL denotes the determinate part, xNUIxNU the indeterminate part, and IxNUϵ[IxLU,IxUU] the degree of indeterminacy. Suppose f(xNU)=f(xLU)+f(xUU)INU;INUϵ[ILU,IUU] presents the neutrosophic probability density function (npdf) of neutrosophic uniform distribution (NUD). Note that the npdf of NUD is based on two parts. The first part xNL, f(xLU) denotes the determinate part and presents the probability density function (pdf) of uniform distribution under classical statistics. The second part xNUIxNU, f(xUU)INU denotes the indeterminate part and IxNUϵ[IxLU,IxUU], INUϵ[ILU,IUU] are the measures of indeterminacy associated with neutrosophic random variable and the uniform distribution. The npdf of the uniform distribution by following [22] is given as
Note here that the npdf of uniform distribution is a generalization of pdf of the uniform distribution. The neutrosophic uniform distribution reduces to the classical uniform distribution when IxUU = 0. The neutrosophic cumulative distribution function (ncdf) of the neutrosophic uniform distribution is given by
Note that the first part presents the cumulative distribution function (cdf) of the uniform distribution under classical statistics, and the second part is the indeterminate part associated with ncdf. The ncdf reduces to cdf when IxUU = 0. The simplified form of ncdf of the Uniform distribution when L=U=S can be written as
where α and β are the scale and shape parameters of the Weibull distribution. The npdf of the Weibull distribution reduces to pdf of the Weibull distribution when INS=0. The ncdf of the Weibull distribution is expressed by
F(xNSW)=1−{e−(xNSWα)β(1+INW)}+INW;INWϵ[ILW,IUW].
(8)
The ncdf of the Weibull distribution reduces to cdf of the Weibull distribution under classical statistics when INW = 0. The neutrosophic mean of the Weibull distribution is given as [22]
μNW=αΓ(1+1/β)(1+INW);INWϵ[ILW,IUW].
(9)
The neutrosophic median of the Weibull distribution is given by
˜μNW=α(ln(2))1/β(1+INW);INWϵ[ILW,IUW].
(10)
4.
Simulation methodology
This section presents the methodology to generate random variates from the proposed neutrosophic uniform distribution and the neutrosophic Weibull distribution. Let uNϵ[uL,uU] be a neutrosophic random uniform from uN∼UN([0,0],[1,1]). The neutrosophic random numbers from NUD and NWD will be obtained as follows:
Step-4: From the routine, the first value of xNSW will be generated.
Step-5: Repeat the routine k times to generate k random numbers from NWD.
4.1. Examples
To illustrate the proposed simulation methods, two examples will be discussed in this section.
4.1.1. Example 1
Suppose that xNSU is a neutrosophic uniform random variable with parameters ([20,20],[30,30]) and a random variate xNSU under indeterminacy is needed. To generate a random number from NUD, the following steps have been carried out.
Step-1: Generate a uniform random number uN=0.05 from uN∼UN([0,0],[1,1]).
Step-2: Fix the values of INS=0.1.
Step-3: Generate values of xNSU using the expression xNSU=20+(0.05(1+0.1))(30−20)=20.5.
Step-4: From the routine, the first value of xNSU=20.5 will be generated.
Step-5: Repeat the routine k times to generate k random numbers from NUD.
4.1.2. Example 2
Step-1: Generate a uniform random number uN=0.30 from uN∼UN([0,0],[1,1]).
Step-2: Fix the values of INS=0.20, α=5, and β=0.5.
Step-3: Generate values of xNSW using the expression xNSW=5[−ln(1−(uN−INW)1+INW)]10.5=0.04.
Step-4: From the routine, the first value of xNSW=0.04 will be generated.
Step-5: Repeat the routine k times to generate k random numbers from NWD.
5.
Simulation study
In this section, random numbers are generated by simulation using the above-mentioned algorithms for NUD and NWD. To generate random numbers from NUD, several uniform numbers are generated from uN∼UN([0,0],[1,1]) and placed in Tables 1 and 2. In Tables 1 and 2, several values of INS are considered to generate random numbers from the NUD. Table 1 is depicted by assuming that NUD has the parameters aNS = 10 and bNS = 20 and Table 2 is shown by assuming that NUD has the parameters aNS = 20 and bNS = 30. From Tables 1 and 2, the following trends can be noted in random numbers generated from NUD.
Table 1.
Random numbers from NUD when aNS = 10 and bNS = 20.
1) For fixed INS, aNS = 10 and bNS = 20, as the values of u increase from 0.05 to 0.95, there is an increasing trend in random numbers.
2) For fixed u, aNS = 10 and bNS = 20, as the values of INS increase from 0 to 1.1, there is a decreasing trend in random numbers.
3) For fixed values of u and INS, as the values of aNS and bNS increases, there is an increasing trend in random numbers.
The random numbers for NWD are generated using the algorithm discussed in the last section. The random numbers for various values of u, INS, α, and β are considered. The random numbers when α=5 and β=0 are shown in Table 3. The random numbers when α=5 and β=1 are shown in Table 4. The random numbers when α=5 and β=2 are shown in Table 5.
Table 3.
Random numbers from NUD when α=5 and β=0.5.
In this section, the performance of simulations using classical simulation and neutrosophic simulation will be discussed using the random numbers from the NUD and the NWD distribution. As explained earlier, the proposed simulation method under neutrosophy will be reduced to the classical simulation method under classical statistics when no uncertainty is found in the data. To study the behavior of random numbers, random numbers from NUD when INS=1.1, aNS = 20, and bNS = 30 are considered and depicted in Figure 3. In Figure 3, it can be seen that the curve of random numbers from the classical simulation is higher than the curve of random numbers from the neutrosophic simulation. From Figure 3, it is clear that the proposed neutrosophic simulation method gives smaller values of random numbers than the random numbers generated by the neutrosophic simulation method. The random numbers from NWD when INS=0.9, α=5, and β=0.5 are considered and their curves are shown in Figure 4. From Figure 4, it can be seen that random numbers generated by neutrosophic simulation are smaller than the random numbers generated by the classical simulation method under classical statistics. The random numbers generated by the neutrosophic simulation are close to zero. The random numbers from NWD when INS=0.1, α=5, and β=1 (exponential distribution) are considered and their curves are shown in Figure 5. From Figure 5, it can be seen that the curve of random numbers generated by neutrosophic simulation is lower than the curve of random numbers generated by the classical simulation method under classical statistics. The random numbers from NWD when INS=0.1, α=5, and β=2 are considered and their curves are shown in Figure 6. From Figure 6, it can be seen that the curve of random numbers generated by neutrosophic simulation is lower than the curve of random numbers generated by the classical simulation method under classical statistics. From Figures 4–6, it can be concluded that the proposed simulation gives smaller values of random numbers as compared to the classical simulation method under classical statistics.
Figure 3.
Random numbers behavior from NUD when INS=1.1, aNS = 20, and bNS = 30.
The simulation method under neutrosophic statistics and classical methods was discussed in the last sections. From Tables 1 and 2, it can be seen that random numbers from the NUD can be generated when INS<1, INS=1 and INS>1. On the other hand, the random numbers from the NWD can be generated for INS<1, INS=1, and INS>1 when the shape parameter β<1. From Table 4 and 5, it can be noted that for several cases, the NWD generates negative results or random numbers do not exist. Based on the simulation studies, it can be concluded that the NWD generates random numbers INS<1, INS=1, and INS>1 only when β<1. To generate random numbers from NWD when β≥1, the following expression will be used
xNSW=α[−ln(1−uN+INW1+INW)]1β;1−uN+INW≥0.
8.
Concluding remarks
In this paper, we initially introduced the NUD and presented a novel method for generating random numbers from both NUD and the NWD. We also introduced algorithms for generating random numbers within the context of neutrosophy. These algorithms were applied to generate random numbers from both distributions using various parameters. We conducted an extensive discussion on the behavior of these random numbers, observing that random numbers generated under neutrosophy tend to be smaller than those generated under uncertain environments. It is worth noting that generating random numbers from computers is a common practice. Tables 1–5 within this paper offer valuable insights into how the degree of determinacy influences random number generation. Additionally, these tables can be utilized for simulation purposes in fields marked by uncertainty, such as reliability, environmental studies, and medical science. From our study, we conclude that the proposed method for generating random numbers from NUD and NWD can be effectively applied in complex scenarios. In future research, exploring the statistical properties of the proposed NUD would be advantageous. Additionally, investigating the proposed algorithm utilizing the accept-reject method could be pursued as a future research avenue. Moreover, there is potential to develop algorithms using other statistical distributions for further investigation.
Use of AI tools declaration
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
Acknowledgements
The authors are deeply thankful to the editor and reviewers for their valuable suggestions to improve the quality and presentation of the paper.
Conflict of interest
The authors declare no conflicts of interest.
References
[1]
Ismail AH, Kusbiantoro A, Chin SC, et al. (2020) Pozzolanic reactivity and strength activity index of mortar containing palm oil clinker pretreated with hydrochloric acid. J Clean Prod 242: 118565. doi: 10.1016/j.jclepro.2019.118565
[2]
Tee KF, Mostofizadeh S (2021) A mini review on properties of Portland cement concrete with geopolymer materials as partial or entire replacement. Infrastructures 6: 26. doi: 10.3390/infrastructures6020026
[3]
Mostofizadeh S, Tee KF (2019) Evaluation of impact of fly ash on the improvement on type II concrete strength. Proceedings of the 39th Cement and Concrete Science Conference, 190-193.
[4]
Gharehbaghi K, Tee KF, Gharehbaghi S (2021) Review of geopolymer concrete: A structural integrity evaluation. IJFE (in press).
[5]
Samad AAA, Mohamad N, Ali AZM, et al. (2017) Trends and development of green concrete made from agricultural and construction waste, Advanced Concrete Technology and Green Materials, Johor: Penerbit UTHM.
[6]
Tee KF, Agba IF, Samad AAA (2019) Experimental and numerical study of green concrete. Int J Forensic Eng 4: 238-254. doi: 10.1504/IJFE.2019.106677
[7]
Samad AAA, Hadipramana J, Mohamad N, et al. (2018) Development of green concrete from agricultural and construction waste, In: Sayigh A, Transition Towards 100% Renewable Energy, Springer, Cham, 399-410.
[8]
Golewski GL (2020) Changes in the fracture toughness under mode II loading of low calcium fly ash (LCFA) concrete depending on ages. Materials 13: 5241. doi: 10.3390/ma13225241
[9]
Golewski GL, Gil DM (2021) Studies of fracture toughness in concretes containing fly ash and silica fume in the first 28 Days of curing. Materials 14: 319. doi: 10.3390/ma14020319
[10]
Golewski GL (2021) The beneficial effect of the addition of fly ash on reduction of the size of microcracks in the ITZ of concrete composites under dynamic loading. Energies 14: 668. doi: 10.3390/en14030668
[11]
Lisantono A, Wigroho H, Purba R (2017) Shear behavior of high-volume fly ash concrete as replacement of Portland cement in RC beam. Procedia Eng 171: 80-87. doi: 10.1016/j.proeng.2017.01.312
[12]
Annapurna D, Kishore R, Sree MU (2016) Comparative study of experimental and analytical results of geopolymer concrete. IJCIET 7: 211-219.
[13]
Karnati V (2016) Flexural response of reinforced concrete beams using various cementitious materials [Master's Thesis]. University of Toledo.
[14]
Dahmani L, Khennane A, Kaci S (2010) Crack identification in reinforced concrete beams using ANSYS software. Strength Mater 42: 232-240. doi: 10.1007/s11223-010-9212-6
[15]
Vasudevan G, Kothandaraman S, Azhagarsamy S (2013) Study on non-linear flexural behavior of reinforced concrete beams using ANSYS by discrete reinforcement modeling. Strength Mater 45: 231-241. doi: 10.1007/s11223-013-9452-3
[16]
Ansys Fluent 14.0: User's guide, 2011. Available from: https: //www.scribd.com/doc/140163383/Ansys-Fluent-14-0-Users-Guide.
[17]
Nguyen K, Ahn N, Le T, et al. (2016) Theoretical and experimental study on mechanical properties and flexural strength of fly ash-geopolymer concrete. Constr Build Mater 106: 65-77. doi: 10.1016/j.conbuildmat.2015.12.033
[18]
Popovics S (1970) A review of stress-strain relationships for concrete. ACI J 67: 243-248.
[19]
Isojeh B, El-Zeghayar M, Vecchio F (2017) Simplified constitutive model for fatigue behavior of concrete in compression. J Mater Civil Eng 29: 04017028. doi: 10.1061/(ASCE)MT.1943-5533.0001863
[20]
Buckhouse ER (1997) External flexural reinforcement of existing reinforced concrete beams using bolted steel channels [Master's Thesis]. Milwaukee: Marquette University.
[21]
Halahla A (2019) Identification of crack in reinforced concrete beam subjected to static load using non-linear finite element analysis. CEJ 5: 1631-1646. doi: 10.28991/cej-2019-03091359
[22]
Kachlakev D, Miller TR, Yim S, et al. (2001) Finite element modeling of concrete structures strengthened with FRP laminates. FHWA-OR-RD-01-17, Oregon Department of Transportation Research Section
[23]
Desai YM, Mufti AA, Tadros G (2002) User manual for FEM PUNCH, Version 2.0. ISIS Canada.
[24]
Dawari V, Vesmawala G (2013) Modal curvature and modal flexibility methods for honeycomb damage identification in reinforced concrete beams. Procedia Eng 51: 119-124. doi: 10.1016/j.proeng.2013.01.018
[25]
Chahrour A, Soudki K (2005) Flexural response of reinforced concrete beams strengthened with end-anchored partially bonded carbon fiber-reinforced polymer strips. J Compos Constr 9: 170-177. doi: 10.1061/(ASCE)1090-0268(2005)9:2(170)
[26]
Ramazanianpour A, Malhotra V (1995) Effect of curing on the compressive strength, resistance to Chloride-ion penetration and porosity of concrete incorporating slag, fly ash or silica fume. Cement Concrete Comp 17: 125-133. doi: 10.1016/0958-9465(95)00005-W
Getahun MA, Shitote SA, Gariy ZCA (2018) Experimental investigation on engineering properties of concrete incorporating reclaimed asphalt pavement and rice husk ash. Buildings 8: 115. doi: 10.3390/buildings8090115
[29]
Xie J, Kayali O (2016) Effect of superplasticizer on workability enhancement of class F and class C fly ash-based geopolymers. Constr Build Mater 122: 36-42. doi: 10.1016/j.conbuildmat.2016.06.067
[30]
Nematollahi B, Sanjayan J (2014) Effect of different superplasticizers and activator combinations on workability and strength of fly ash based geopolymer. Mater Design 57: 667-672. doi: 10.1016/j.matdes.2014.01.064
[31]
Alrefaei Y, Wang Y, Dai J (2019) The effectiveness of different superplasticizers in ambient cured one-part alkali activated pastes. Cement Concrete Compos 97: 166-174. doi: 10.1016/j.cemconcomp.2018.12.027
Kong Fah Tee, Sayedali Mostofizadeh. Numerical and experimental investigation of concrete with various dosages of fly ash[J]. AIMS Materials Science, 2021, 8(4): 587-607. doi: 10.3934/matersci.2021036
Kong Fah Tee, Sayedali Mostofizadeh. Numerical and experimental investigation of concrete with various dosages of fly ash[J]. AIMS Materials Science, 2021, 8(4): 587-607. doi: 10.3934/matersci.2021036