## 3 Seismic Parameters Model Development  
Different criteria such as stiffness, energy absorption, ductility, and seismic factors are discussed herein [51; 52]. All parameters were computed numerically using the load-displacement curve of each group, and the results are compared. The simplified procedures of converting a full graph into a bilinear graph are shown in Figure 4. This figure shows the procedures that calculate the parameters needed for the finite element analysis [53]. As can be seen in this figure, the first step is to add the load-displacement diagram. This diagram is added based on finite element results. The second step is to find out the first plastic hinge. The plastic hinge usually occurs after the crack tip and at the time of crack propagation. It is worth mentioning that the plastic hinge is automatically calculated with the developed script. In step three, the developed script is defined to calculate the yield and failure load and corresponding displacement. In the next step, the time and \"k\" factor are defined. In step five, the data of step three are employed to calculate the ductility factor. Step six is scripted to choose a different method of computing factor of \"R\". In the final steps (seven and eight), the accuracy of all previous steps is determined, based on two main conditions (see step eight). In fact, these eight steps are subjected to recursion to achieve the required accuracy. In Figure 5, the approach of the parametrical study is shown. In general, the concept of finding these parameters is carried out by converting the load-displacement diagram into a bilinear diagram. In fact, two linear parts of the bilinear diagram contributed to the calculations of the parameters based on a boundary between the linear and nonlinear parts. In Figure 6, the parametrical approach of the bilinear diagram of group FS (group F with steel bars) is presented. In order to create a two-linear graph, the first plastic hinge and its corresponding displacement were considered. It should be mentioned that the failure criteria are the maximum load and displacement. After finding the yield force and its corresponding displacement by employing linear analysis of Equation (2), the boundary of linear and nonlinear parameters is generated. Further, finding a linear analysis parameter requires the lateral force which can be computed by \(C\) and \(W\) in which \(C\) = ABI/R and \(W\) is the weight of the RC element [51; 52] (A: Design base acceleration (in relation to gravity acceleration g), B: Response coefficient of the building obtained from the design response spectrum, I: Importance factor and R is the Behavior coefficient of the building). \[V=C\cdot W \tag{2}\] \begin{table}
\begin{tabular}{c c c c c c} \hline \hline
**Group of** & **Steel (S)** & **CFRP Bars** & **GFRB Bars** & **Number of** & **Optimized** \\
**Specimens** & **(C)** & **(G)** & **(G)** & **Dangitudinal** & **CFRP with** \\  & & & & **Bars \#25** & **\#12 Bars** \\ \hline A & AS & AC & AG & 6 & AO \\ \hline B & BS & BC & BG & 6 & BO \\ \hline C & CS & CC & CG & 8 & CO \\ \hline D & DS & DC & DG & 12 & DO \\ \hline E & ES & EC & EG & 4 & EO \\ \hline F & FS & FC & FG & 4 & FO \\ \hline G & GS & GC & GG & 26 & GO \\ \hline H & HS & HC & HG & 24 & HO \\ \hline \hline \end{tabular}
\end{table} 
Table 4: Model’s Name. Figure 4: Flowchart for Seismic Calculation. Figure 5: The general response of a structure. On the other hand, there are some conditions that are needed to confirm and control the obtained values. One of the most important conditions is controlled by the area under the graph. It indicates that linear and nonlinear area of the bilinear graph that has been computed. It should be noted that the area under the curve of linear and nonlinear parts must be almost the same. A small difference in the computed area may still be considered as high accuracy of the computation.