## 1 Introduction
It opposes the assumption of the weaker behavior of weld lines. The previously mentioned studies stand as proof of the complexity of 3D printed parts behavior and reveal many issues to model such materials. Numerous studies investigate the effect of raster orientation on the mechanical behavior of 3D printed parts through experimentation [13; 14; 15; 16; 17; 18; 19]. To investigate the orthotropic elastic properties of 3D printed parts, [PERSON] et al. [20] developed computational models based on microstructures of 3D printed ABS using micromechanics of a representative volume element (RVE). In a similar fashion, [PERSON] et al. [21] employed asymptotic homogenization as a tool for predicting the mechanical properties of fused filament fabricated thermoplastics with line infill and grid infill structures. Both [PERSON] et al. [22] and [PERSON] et al. [2] used the extended finite element method (XFEM) with cohesive segments to capture the dependency of fracture behavior on build/raster orientations. The classical laminate theory was also utilized to model orthotropic mechanical properties of 3D printed parts [19; 23]. Although research has illuminated the effect of raster orientation through different modeling methods and experimentation, no study to date has studied complex deposition trajectories. Previous studies have almost exclusively focused on modeling unidirectional or bidirectional raster orientations. No numerical study to date has examined complex deposition trajectories. In this paper, the authors intend to propose a simplified finite element method that surmounts the weld lines and the air-gap problems by approaching the behavior in a global way. The model considers the material as homogenous and reproduces the built trajectories using local references assignment in the mesh element [24]. A transverse isotropic law is therefore defined at each integration point to model the behavior in the elastic region. The material constants were determined from a tensile specimen fabricated using the same parameters as the pre-structured material. Then, the Hill criterion coupled with a plastic flow model was identified to go until the fracture limit. The finite element software ABAQUS has been used to implement the identified elastoplastic behavior. Nevertheless, local phenomena such as filaments buckling and air gap repartition are not modeled. The related problems will be discussed later on where this numerical tool will be tested on the compact tension CT specimen [4] and the Single Edge Notched Bend SENB specimen defined in [7] where both the optimized deposition and the classical deposition trajectories will be modeled and compared. The structure of the paper is as follows. In Section 2, the modeling strategy adopted for the 3D printed material is presented. First, the local references assignment technique is explained. Then the pre-structured material behavior definition is presented. Next, the experimental setups and results are displayed. In Section 4, the implementation of the material's behavior in ABAQUS is detailed. In Section 5, the model will be used to predict the tensile tests. Afterward, it will be applied to the CT and SENB specimens. Here, its capability of implementing complex deposition trajectories will be checked and the limitations of the proposed method will be resumed in Section 6.