## 5 Results and Discussions
### Tensile Tests Prediction  
The finite element model for the tensile specimen geometry was created in ABAQUS software. Three models were later defined by assigning local material references with respect to the printing trajectories. Three different simulations have been carried out to simulate the behavior of each of the specimens according to the suggested model. The part of the specimen pinched between the fixed jaws is encastered and a displacement is applied on a reference point coupled with the part pinched between the mobile jaws.  
Figure 7 compares the stress-strain curves obtained experimentally with those obtained by finite element simulations. The results show that for all three cases, the elastic behavior is well reproduced with a slight difference because of the simplifications described before. Thus, the model predicts well the stress evolution until the yielding point. By comparing the model response for the three different orientations, it is clear that the anisotropic yielding criterion taken in the model captures the yield point in the experimental curves. The transition between the elastic and the plastic domain for the \(0^{\circ}\) specimen occurs around  
\begin{table}
\begin{tabular}{c c c c} \hline \hline
**Ratio** & \(R_{L}\) & \(R_{T}\) & \(R_{LT}\)=\(R_{LT}\) \\ \hline Expression & \(\frac{\sigma_{L}^{0}}{e^{0}}\) & \(\frac{\sigma_{L}^{0}}{e^{0}}\) & \(\frac{\sigma_{T}^{0}}{e^{0}}=\frac{\sigma_{T}^{0}}{e^{0}}\) \\ Value & 1 & 0.565 & 0.857 \\ \hline \hline \end{tabular}
\end{table}
Table 4: Implementation of [PERSON]’s potential function in ABAQUS using the yield stress ratios for the transverse isotropic material.  
\begin{table}
\begin{tabular}{c c c c} \hline \hline \(\sigma_{limit}\)(MPa) & & \(m\) & \\ \hline Mean value & Std & Mean value & Std \\
6.58 & 3.41 & 221.6 & 96.19 \\ \hline \hline \end{tabular}
\end{table}
Table 5: Exponential law fitting results.  
the value of \(\sigma_{L}^{\emptyset}\). In addition, for this orientation, the saturation of the stress at the average tensile strength is noticed (see Figure 7a). The onset of plastic behavior observed for the 45\({}^{\circ}\) tensile sample is well simulated (see Figure 7b) by coupling [PERSON]'s criterion with the 0\({}^{\circ}\) flow curve. Therefore, the numerical approach is able to highlight the anisotropic strengths according to the filament orientation. For the sample with 90\({}^{\circ}\) (Figure 7c), orientation exhibits a brittle fracture. However, a plastic behavior appears in the numerical response. This behavior was intentionally discarded since it does not appear in the experimental tests. The elongation at fracture was managed by applying the appropriate corresponding displacement value without implementing a fracture initiation criterion. The average values reported in Table 2 of the fracture elongation were respected.  
The model allows good prediction of the elastic domain without any estimation of the air gap. The plastic component was defined in order to demonstrate the ability of [PERSON]'s criterion to predict the anisotropy with the material yield strength. By only implementing the plastic flow curve for the 0\({}^{\circ}\) orientation, it has been shown that it is possible to predict the mechanical response of the material in different orientations.