Experimental and Numerical Analysis of SMC Compression Molding in Conﬁned Regions—A Comparison of Simulation Approaches

: The compression molding process of sheet molding compound (SMC) is an economical manufacturing process for lightweight parts. However, molding defects, such as ﬁber matrix separation, and ﬁber re-orientation, may develop during the molding process in conﬁned regions, such as ribs and bosses. Hence, the mechanical properties of the composite depend on the local ﬁber architecture. Consequently, this work compares the predictive capabilities of tensor-based and directly modeled process simulation approaches regarding compression force, ﬁber volume content and ﬁber orientation on the example of honeycomb structures molded from SMC. The results are validated by micro-computed tomography and thermal gravimetric analysis. The ﬁber orientation in the honeycomb varies between individual samples because a sheet molding compound is macroscopically heterogeneous and thus the ﬁber architecture is strongly inﬂuenced by random events. Tensor-based ﬁber orientation models can not reliably predict ﬁber volume content and ﬁber orientation in the part’s thickness direction if there is a lack of scale separation. Therefore, directly modeled process simulations should be preferred in cases in which ﬁber length and mold dimensions prohibit scale separation. The prediction of ﬁber volume content is a difﬁcult task and no simulation can predict the severity of ﬁber matrix separation precisely in all cases.


Introduction
The transition from conventional vehicles to battery electric vehicles creates new challenges for the structural design process. The new crash management system (CMS) has to overcome the missing substructure of the combustion engine with regard to crash safety and the increase in mass caused by the battery system by means of improved support structures.
One approach is the local reinforcement of metallic crash elements with fiber reinforced polymers (FRP), which offer high mass-specific stiffness and strength. The energy absorption capacity of the crash elements is an important factor here. For the specific energy absorption, honeycomb structures show a high potential. As crash elements, they are already being used by various car manufacturers and are mainly made of discontinuous fiber-reinforced polymers. An important advantage of discontinuous FRP are their low material and manufacturing costs, as well as the high freedom of design. The excellent mechanical properties of FRP are directly related to the fiber length, fiber volume fraction and the fiber orientation in the component. An increase in these properties up to a certain level while maintaining the complexity and design freedom of the manufacturing process would thus have a positive effect on the crashworthiness of the hybrid metal-FRP-CMS [1,2].
One possibility is the use of sheet molding compound (SMC), which has become a very attractive material in recent years, especially for exterior parts in the automotive industry. SMC consists of a thermosetting resin reinforced with chopped fibers and in most cases mineral fillers. SMC components are usually manufactured in a compression molding process, which allows the processing of long-fiber-reinforced polymers. Compared with injection-molded parts, compression-molded parts are superior in terms of increased strength and toughness due to the higher fiber length.
During the compression molding of SMC, the initial fiber architecture is changed and is influenced by process parameters, e.g., the initial charge position or the manufacturing process conditions. At first, research on local changes in microstructure was mainly focused on fiber length and orientation rather than changes in FVC, which is commonly known as fiber matrix separation (FMS). It develops if fibers are not transported homogeneously with the matrix during mold filling. As a result, fiber content becomes non-uniform, even when averaged over macroscopic volumes. This leads to a lower than nominal FVC along the flow path and at the flow front in particular [3]. In this process, FMS results from a balance of forces between suspended fibers, fibers and mold surface, the elastic deformation of the fibers, and hydrodynamic forces acting as a result of fluid flow [4][5][6]. Basic experiments have shown that FMS is also formed in flat, round or rectangular compression-molded stacks [7,8] and depends, for instance, on the fiber-aspect ratio, the matrix viscosity [9] and the closing speed during production [10]. Usually, FMS is most pronounced in the top sections of ribs [11][12][13], which can be coupled with an increased fiber density in the rib base [14].
This fiber architecture has a direct influence on the mechanical behavior of the molded part. The effects of process parameters on the structural-mechanical component properties must therefore be mapped with a functional virtual process chain in order to create enhanced confidence in the class of discontinuous long-fiber-reinforced composites, as is already standard in metal processing with commercial software linking tools. Such a coupling has already been considered in short-fiber thermoplastic injection molding with a validation using a demonstrator and is also available for commercial use [15]. For the compression molding process of long-fiber-reinforced thermoplastics, Buck et al. [16] showed the effect of coupling process and structure simulation. Görthofer et al. [17] presented an approach for an integrated virtual process chain for SMC compression molding.
However, the base of any virtual process chain is a reliable and efficient process simulation, which is able to provide high-quality result data to map the material-specific characteristics of the specific material class. A brief overview of the state of research in fiber reorientation modeling is described below. Basically, a distinction is made between macroscopic models, which describe the fiber composite as a single-phase medium, and microscopic or mesoscopic approaches, which differentiate between fiber and matrix. Macroscopic models are used to predict the fiber orientation and pressure distribution in components and are already integrated in commercial software environments. Fiber orientation prediction started in 1922 with G. B. Jeffery's approach to Einstein's publication on the behavior of spherical particles in a liquid [18]. Jeffery extended Einstein's work [19] for ellipsoidal particles and described their rotation in a dilute Newtonian flow. The velocity field of the fluid is not influenced by the particle and the contact between particle and fluid is assumed to be perfect. This approach was extended by Folgar and Tucker, who developed a phenomenological model to predict the orientation distribution of fibers in concentrated suspensions [20]. Their derived model is regarded as the initial model for the prediction of fiber orientation in process simulation. Advani and Tucker [21] proposed the use of statistical moments of the fiber orientation distribution function (called fiber orientation tensors and denoted A here) to efficiently describe the orientation state in numerical models. Some experimental shortcomings, especially regarding the limitation of fiber length, have led to extended models by Wang et al. [22], Phelps et al. [23] and Tseng et al. [24]. These tensor-based models provide good results for FRP with fiber lengths smaller than the characteristic length of the component [25], i.e., if scale separation applies. This term refers to the assumption that the fiber architecture is a local state of material without any significant extension at the part scale. It may be violated if the fiber length is in the order of magnitude of geometrical features.
All macroscopic modeling approaches published to date are phenomenological and cannot detail the specific mechanics during fiber-fiber interaction. In highly concentrated fiber suspensions, however, fiber interactions are the dominant factor for modeling fiber re-orientation and distribution. For this reason, such methods seem to be of limited suitability for use in a virtual process chain for long-fiber-reinforced plastics with confined regions.
With simulation methods on the microscopic or mesoscopic level, fiber movements can be simulated for discretized fibers. In contrast to the phenomenological macroscopic fiber orientation models, modeling approaches on the microscopic scale allow a more precise approximation of the physical behavior [22,[26][27][28][29][30]. At the micro level, different modeling approaches can be used. They can be roughly divided into particle-based methods, in which the polymer matrix and the fibers are treated as particles, or element-based methods, in which fibers are treated as particles and the matrix as a continuous medium. In contrast to particle-based methods, which are by nature two-way coupled, element-based simulations are often solved with one-way coupling. The feedback from the fiber motion on the fluid is computationally complex, but can be integrated in a general or refined way [31,32].
In this work, two direct element-based simulation methods with different coupling mechanisms between fiber and matrix and two macroscopic models are investigated for the process simulation of the SMC compression molding process. Furthermore, they are tested for their predictive accuracy using a real structure manufactured with different process conditions, based on fiber orientation, distribution, fiber matrix separation, compression forces and short-shot flow front.
An important factor for this comparison is the reproducibility of fiber orientation and distribution under real conditions, especially in complex three-dimensional structures. To the authors' knowledge, previous works did not address this problem or only in geometrically simple structures, where the fiber length and component dimensions are at least one order of magnitude apart. Instead, this study examines the filling of a complex honeycomb structure with an SMC whose fiber length is in the same order of magnitude as the part's geometric features. The process simulations are compared and validated by their prediction of the compression force, fiber volume content and fiber orientation inside a complex ribbed honeycomb structure.

Materials and Specimen Production
The material processed is an SMC named HUP EJ 43529 produced by Polynt Composites Germany GmbH, Germany. It contains 54 wt % glass fiber bundles, 36.8 wt % vinyl ester resin and 9.2 wt % CaCO 3 (measured by thermal gravimetric analysis (TGA) of SMC plates with 95 % mold coverage). Therefore, the nominal fiber volume content φ n is approximately 36.6 %. Glass fibers are 50 mm long, have a diameter of 20 µm and are assembled in bundles comprising about 560 individual fibers.
The honeycomb structure is a demonstrator for a complex rib structure that will be molded into a metal profile which serves as an automotive CMS. The mold cavity features two honeycomb structures with different wall thicknesses and draft angles (see Figure 1a). The different geometries are chosen to determine the influence of rib width on the development of fiber matrix separation (FMS) and molding defects, such as pores. For this study, only the results of the thick honeycomb walls are considered. The influence of the geometry (i.e., different wall thicknesses) on fiber architecture is not a focus of this work.
The thick honeycomb consists of walls with 2.5 mm thickness at the top (see Figure 1c) and a draft angle of 2°. The base plate has an area of 310 mm by 310 mm and is 2 mm thick.    The compression molding is conducted on a double-column press of type RSP 160.42 SZ by Röcher GmbH & Co. KG Maschinenbau, Germany. To investigate the influence of different initial charge positions, two process configurations named Configuration A and Configuration B are specified. The initial charge position of Configuration A is underneath the honeycombs, whereas the initial charge is positioned opposite to the honeycombs in Configuration B (see Figure 1a). The process parameters are summarized in Table A1.
For the investigation, only one of three honeycomb cells, a coupon with six walls, is used. A coupon consisting of these six walls is preferred because the individual walls can be shown side by side, thus giving a better overview of the 3D fiber architecture. A pre-trial showed that this specific cell reliably contained the defects, which are the focus of this work, and it is filled regularly to obtain equivalent samples. Figure 1b shows the coupon and the representation of its walls side by side. The dotted red lines indicate the positions where the walls are cut into smaller specimens for TGA.
The fiber matrix separation index (FMSI) is introduced as a means to quantify the overall peculiarity of fiber matrix separation in the honeycomb. It is the average of the normalized difference between the experimentally found FVC in a specimen φ i and the nominal FVC φ n . The volume weighting ν i takes the different volumes of the specimens into account. Thus, the FMSI is calculated by where ν S denotes the sum over all weighting factors ν i . The specimen partition presented in Figure 1b is a trade-off between sufficient spatial resolution and reasonable workload. Three repetitions are conducted for each process configuration.

Matrix Rheology
The dynamic matrix viscosity is characterized by small-amplitude oscillatory shear (SAOS) in a plate-plate geometry with a diameter of 25 mm. The rheological measurements are conducted on a Haake Mars III by Thermo Fisher Scientific Inc., United States. The matrix material is stored for three weeks at room temperature in a closed container before the rheological measurements. This time is similar to the maturation time of the SMC before delivery and subsequent compression molding.

SMC In-Mold Rheology and Friction at the Tool Surface
The reliable material parameter identification is one of the most important steps for the process simulation, especially considering the prediction of press forces. For this reason, a near-process characterization was carried out employing a press rheometer on a Dieffenbacher COMPRESS PLUS DCP-G 3600/3200 AS hydraulic press by Dieffenbacher Holding GmbH & Co. KG., Germany. A press rheometer is a simple rectangular plate tool, which is equipped with several pressure sensors. The pressure sensors are used both to measure the local pressure and to identify the position of the current flow front. Such a tool was initially developed by Castro and Griffith [36] and improved by Dumont et al. [37][38][39], as well as Hohberg et al. [40][41][42].
For the press rheometer trials, a rectangular cavity of 800 mm × 450 mm and 9 pressure sensors were used at the Fraunhofer Institute for Chemical Technology (ICT), Pfinztal, Germany. The molding trials were performed with the standard temperatures given by the data sheet. The initial charge coverage varied between 37.5% (300 mm × 450 mm) and

Micro-Computed Tomography
Micro-computed tomography (µCT) is a non-destructive testing method which allows the analysis of a sample's internal material structure by reconstructing a 3D object from X-ray images. The resolution of the 3D object, the so-called voxel, depends on the distance between specimen and X-ray source, or rather the size of specimen and detector. For fiber orientation analysis of individual fibers, voxel size needs to be smaller than the fiber diameter. In contrast, the orientation analysis of fiber bundles only requires a resolution that is smaller than the fiber bundle diameter, as fibers tend to remain assembled in bundles during mold filling [43]. Table A3 gives an overview about the parameters used for the µCT of coupons.

Numerical Methods
All numerical models assume negligible curing [44][45][46] and a planar isotropic initial fiber orientation state, because fiber bundles fall irregularly onto the conveyor belt during SMC prepreg production. Besides the characterizations performed within this work, the SMC viscosity was also characterized by Toray Engineering Co., Ltd., Japan via press rheometry and its viscosity is modeled via a Cross-Andrade model [47]. During the rheological characterization of the SMC by Toray, there was no differentiation between the contributions of viscous forces and friction forces. Thus, the viscosity measurement by Toray includes the friction between SMC and press rheometer. This can be interpreted as treating the friction and viscosity in a combined way, but it does not allow for a separation of effects due to friction and the material itself, which can lead to incorrect results in cases that are different to the characterization conditions. Still, the characterization by Toray was used because it fits the material model used in 3D TIMON. The resulting material data are used as input for MOLDEX3D and 3D TIMON (see Appendix Table A2).

Macroscopic Process Simulation
Within this investigation, two commercial macroscopic process simulation tools are used: MOLDEX3D by CoreTech System Co., Ltd. and MOLDFLOW by Autodesk, Inc. In both tools, the gradients of the velocity field are used to calculate the fiber orientation evolution with one-phase macroscopic models, assuming that fibers and matrix have the same velocity and, therefore, that no fiber matrix separation can occur. In the models applied in this work, the fiber orientation does not affect the flow, as no anisotropic viscosity models are used. Thus, the fiber orientation calculation could be considered a post-processing step.

MOLDEX3D 2020
The part geometry and compression zone are modeled with triangular prisms with an edge length of 1.5 mm. The part is modeled by 17 elements over its thickness, while the compression zone is discretized by 50 elements in the compression direction. Overall, this results in 2,322,404 and 4,953,872 mesh elements in the part and compression zone, respectively. The compression profile of the press is prescribed as input for MOLDEX3D with a maximum compression force of 1500 kN. The friction coefficient for the contact between SMC and mold is set to 0. This seems appropriate, considering that the viscosity characterized by Toray already includes the friction between SMC and press rheometer. The iARD-RPR model and IBOF-closure are used for the fiber orientation analysis [48,49]. The coefficients for the iARD-RPR model are C I = 0.1, C M = 1 and α = 0.9, following [50].

MOLDFLOW 2021
For the simulation, the full geometry is used, since MOLDFLOW offers no symmetry boundary condition. The mold cavity and therefore the compression zone is modeled with tetrahedrons with an edge length of 3 mm and at least 26 elements over the part thickness, which results in 2,144,015 tetrahedrons elements for the cavity. For the press displacement, a closing speed profile as well as maximum compression force of 1500 kN were given as input, which is identical to the press profile applied in MOLDEX3D. The SMC was characterized for this study, as described in Section 3.3 with the resulting rheological parameters and friction parameters given in Table A5. For the fiber orientation analysis, the Reduced Strain Closure (RSC) Model is used together with the Orthotropic 3 closure [51,52]. The coefficient for the RSC model is C i = 0.1 and the RSC factor is 0.9 according to the MOLDEX3D simulation settings.

Direct Simulation Models at Microscale and Mesoscale
The term direct simulation models is used in this work to classify methods that directly model the motion of individual fiber particles in a suspension. Here, fibers are approximated by rods or spheroids and flexible fibers are represented by connected rods, spheres or spheroids. The fibers might interact through hydrodynamic, elastic, lubrication and friction forces with the surrounding fluid and each other [22,[26][27][28][29]. The effect of fibers on the fluid field can be incorporated, for example, following the work of [31]. One method employed in this work is based on the process simulation software 3D TIMON COMPOSITEPRESS™ by Toray Engineering Co., Ltd., which does not account for interactions between fibers. Another method is based on the approach by [32], which takes fiber-fiber interactions and the effect of fibers on the fluid field into account. The resulting bundle positions of the direct simulation models are exported to VTK and mapped to an evaluation mesh using a ParaView filter (the filter is available at https://github.com/nilsmeyerkit/paraview_map_ lines, accessed on 14 February 2022) to compute fiber orientation tensors and fiber volume fractions from the line elements. The evaluation mesh exactly represents the cells that are used as coupons during the thermal gravimetric analysis and enables the comparison between experimental results and the macroscopic models. The Direct Fiber Simulation (DFS), which is a feature of 3D TIMON [53,54], is used to simulate the motion of long fibers during SMC compression molding. The fiber movement is calculated in a post-processing step after a classical macroscopic process simulation with 3D TIMON COMPOSITEPRESS™ and fiber nodes move with the pre-computed velocity field [55]. This results in a one-way coupling, which means that fiber movement is dictated by the flow field but the flow field is not influenced by the presence or orientation of fibers [56]. Anisotropic flow behavior due to fiber orientation is neglected and fiber-fiber interactions are not considered [14]. Thus, the need to list neighboring fibers and to calculate interaction forces between neighboring fiber segments vanishes. Thereby, the computational effort is reduced immensely.
Fibers are modeled as rods that are connected by nodes which function as joints. During the simulation, joint nodes are moved independently from each other with the polymer flow, whereby a new temporary node position is calculated. Here, x n i is the position of node i at time step n, u n i the velocity of the surrounding fluid around node i at time step n and dt the time increment. This inevitably leads to an unrealistic elongation of fibers during the filling of complex geometries. Therefore, the DFS-algorithm continuously checks the length of each rod in each time step. Once an elongation threshold is exceeded, the temporary node positions are adjusted with a regularization scheme (x n+1 i → x n+1 i ) to maintain a constant fiber length [55]. Contrary to the other methods employed in this work, only the mold half containing the thick honeycomb is modeled to reduce the calculation run-time. This is deemed acceptable because the problem is quasi-symmetric and minor geometrical differences between thick and thin honeycomb hardly have an influence on the resulting flow field. Furthermore, this is favorable because in this study fibers are generated per unit volume of the initial charge. Thus, the number of fibers simulated and therefore calculation run-time were reduced. The part geometry and compression zone were modeled with a 3D voxel mesh with an edge length of 1.2 mm in x-, 1.2 mm in yand 0.7 mm in z-direction, respectively. Additionally, the mesh of the cavity around the part was refined with a 3D voxel mesh with an edge length of 1.2 mm in x-, 1.2 mm in yand 0.2 mm in z-direction. Overall, this results in 1,066,411 mesh elements for the cavity (see also Table A6 for computational parameters). The fibers are modeled with a length of 50 mm and 20 segments per fiber. A total of 145,130 fibers was generated with 2,902,600 fiber segments.
Contrary to the deformation mechanism reported by Barone and Caulk [57,58], the SMC/BMC-module of 3D TIMON assumes a no-slip boundary condition at the mold surface [47]. This results in high shear stresses and flow kinematics that do not match experimental results. Using the original compression force in the force-controlled compression phase would result in unrealistic long compression times. Hence, the maximum compression force was set to 10 5 kN, asserting that the simulation follows the velocity profile and fits the compression time of the compression molding trials (see Section 5.3). This approach appears reasonable because the compression force has little influence on the development of the fluid velocity field, except for the velocity profile in the thickness direction of the SMC stack. The compression speed profile can only be modeled by a series of constant velocities. Therefore, the compression profile of the actual press data can only be approximated and it should be mentioned that jumps or discontinuities in the velocity profile are nonphysical.

Direct Bundle Simulation
Direct Bundle Simulation (DBS) is an approach to directly model SMC processing at mesoscale and is implemented with multiple custom subroutines in Simulia ABAQUS [32]. DBS uses the observation that fiber bundles typically stay intact during the compression molding of SMC [59][60][61]. Thus, the DBS approach treats hundreds of fibers in a bundle as a common instance and therefore drastically reduces the amount of suspended objects in the direct simulation. This enables the accurate prediction of bundle architectures on the component level [32].
The matrix motion is computed based on mass balance and momentum balance, which are solved with a Coupled Eulerian-Lagrangian framework in ABAQUS by operator splitting [62]. This technique enables the large deformation of the polymer and contact with Lagrangian rigid body mold surfaces. The matrix is modeled as a Newtonian isotropic fluid and no heat transfer is considered in this example. However, the incorporation of coupling with the bundle motion renders an overall anisotropic flow behavior. Due to the explicit solving scheme, a weakly compressible equation of state is utilized to model quasi-incompressible behavior. The interaction between both tool surfaces and SMC matrix is modeled with a hard normal contact and a hydrodynamic friction according to with the relative velocity at the mold surface v s , the friction coefficient λ, power-law coefficient m and reference velocity v 0 . It is motivated by a thin lubrication layer of heated matrix at the mold surface [57]. The upper and lower mold surfaces are extracted from the CAD model of the part, placed at an initial distance of 12 mm and are meshed with rigid body elements of 2 mm average edge length. The lower mold remains fixed, while the upper tool is either controlled by a virtual press controller mimicking the press profile (Configuration B), or subjected to a constant force (Configuration A). Configuration A is simulated with a constant compression force, because the press comes effectively to a stop after reaching the stack and proceeds in a force-controlled mode in the trials with this configuration. The cavity is meshed with 1 mm × 1 mm × 1 mm Eulerian elements in the honeycomb area and 2.5 mm × 2.5 mm × 1 mm elements in the planar sections (see also Table A6 for computational parameters). Fiber bundles are represented by one-dimensional truss elements with 2.5 mm length and thus behave like threads. They may bend, but not present any mechanical resistance to bending. The total amount of truss elements is prescribed by the fiber volume fraction assuming 560 fibers per bundle with 20 µm diameter each and results in an equivalent bundle diameter of 0.24 mm. The exact amount varies slightly between individual initial bundle architectures due to the generation procedure. Furthermore, bundle elements are idealized as frictionless contact with other bundles, as well as the mold surface. The interaction between fiber bundles and the fluid phase is implemented with numerous subroutines which compute a body force on bundles based on the velocity relative to the surrounding fluid field ∆v. This hydrodynamic force on a bundle segment j is computed as with viscosity η, bundle radius R and an orthogonal direction q. The factors k d and k l are drag and lift coefficients that depend on the segment aspect ratio and orientation of the bundle [32]. The weighted hydrodynamic forces are used to apply an opposed body force on the fluid phase that ensures a force equilibrium and effectively generates a generally anisotropic flow. Thus, this approach allows the evaluation of fiber matrix separation, the computation of bundle curvature and an accurate prediction of the fiber architecture. In comparison with 3D TIMON, it accounts for fiber-fiber interactions and the effect of fibers on the flow. However, short-range lubrication between fiber bundles, as in [63], is not considered here.
The results of DBS are influenced by the fiber bundle architecture in the initial SMC stack. Thus, the simulation is run with three different initial fiber bundle architectures in each configuration to estimate the scatter due to different random realizations in the initial stack.

Matrix Rheology
The viscosity curve of the matrix shows that the matrix exhibits pronounced shear thinning at higher shear rates (see Figure 2). This behavior is typical for filled polymer systems, polymer melts and polymer solutions [64]. On the one hand, the number of entanglements between neighboring oligomers decreases with increasing shear rate. On the other hand, connections between oligomers and filler particles, as well as connections among filler particles themselves are lost. During frequency sweeps in the nonlinear regime, slip occurs between the plate geometry and the VE matrix at higher shear rates due to the high viscosity of the filled matrix. Hence, no reliable data of the shear thinning behavior are recorded. However, the deformation rates during compression molding of the specific case considered in this work are below 10 s −1 . Subsequently, the VE-matrix is modeled as a Newtonian fluid without temperature dependence in the Direct Bundle Simulation. According to Figure 2, the zero shear viscosity of the matrix is approximately 25,000 Pa s which is used as the input parameter for DBS.

SMC In-Mold Rheology and Friction at the Tool Surface
With these experimental data from Section 3.3, the parameterization for the lubrication layer and its friction parameters is carried out as described by Hohberg et al. [40][41][42]. Therefore, an analytical stress equation, considering the rheological as well as the friction stresses, is fitted to the different experimental molding trials pressure measurements. Here, the data are fitted to the Cross-Andrade model for the effective viscosity with the shear rateγ, zero shear viscosity η 0 , power-law coefficient n, critical shear stress τ * , temperature T, and data-fitting coefficients B and T b . The slip velocity, respectively, the friction, is modeled as with τ w the wall shear stress, the critical wall shear stress τ c and the slip coefficient α w . The resulting fitted values for all these model parameters are given in Table A5. Equivalently, the parameters of Equation (3) are determined. An exemplary fitting result for three representative pressure sensors (S1 at 147 mm from tool edge, S2 at 249 mm and S3 at 451 mm) for an 50% initial charge coverage and a closing speed of 2.5 mm s −1 is given in Figure 3 to illustrate the fitting quality.

Compression Profile
The press data of parts molded with Configuration A show that the compression force varies significantly between successive moldings after reaching the maximum compression force (see Figure 4). Thereby, the force drops by up to 300 kN. The velocity profile shows that the press cannot maintain the programmed closing speed of 10 mm s −1 upon contact with the SMC charge. Interestingly, the compression force does not reach its maximum value of 1500 kN at this point. Here, the upper mold almost comes to a complete stop at approximately 8 s, which might be caused by the grip that the insert with its rib structure exerts on the stack. Then, the compression speed increases again for a short moment and reaches up to 15 mm s −1 . Afterwards, the velocity drops to roughly 3 to 5 mm s −1 and decreases slowly from there on until the mold is closed completely.
The press data of parts molded with Configuration B (see Figure 4) show similar characteristics. The difference in compression force between successive moldings after reaching the maximum compression force is as large as 400 kN. Similarly to the press data of Configuration A, the velocity profile of Configuration B shows that the press cannot maintain the programmed closing speed of 10 mm s −1 upon contact with the SMC charge. Again, the compression force does not reach its maximum value of 1500 kN at this point. Here, the upper mold does not stop completely. Instead, the velocity drops to roughly 3 to 5 mm s −1 . Then, the compression speed increases again for a short moment and reaches up to 15 mm s −1 . Afterwards, the velocity drops to roughly 1 to 2 mm s −1 and decreases slowly from there on until the mold is closed completely.

Fiber Orientation
The reproducibility of fiber orientation states is assessed by scanning the back wall of three identical trials that are molded with process Configuration A. The walls are extracted from the identical location. Figure 5a-c show the three walls side by side. The same volume of each wall is analyzed with a set of evaluation meshes. The evaluation meshes are 2 mm thick in y-direction and 18 mm long in x-and z-direction, which are the directions in the plane of the walls. Each evaluation mesh has the same total size but the number of cells per evaluation mesh is different. Thus, only the cell size and number of cells varies in the in-plane directions (see Figure 5d,e). Cells located at the same position are compared regarding the fiber orientation components A xx , A yy and A zz . The standard deviation of respective components is plotted as a function of the number of cells in Figure 6.   A way to assess the reproducibility of a quantity is to compare its standard deviation (SD) to the standard deviation of a continuous uniform distribution (CUD). If the two standard deviations are approximately the same, this suggests that the variable depends on random processes of equal probability rather than systematic processes. The standard deviation of a CUD is with random variable X that is uniformly distributed on an interval between lower limit a and upper limit b [65]. The lower and upper limit for the out-of-plane orientation A yy are set to 0 and 0.15, respectively. This appears reasonable because the fiber bundles are significantly longer than the wall thickness. Subsequently, the lower and upper limit for the in-plane orientations A xx and A zz are set to 0 and 0.85, respectively. The standard deviation of a random variable X uniformly distributed yields 0.245 for the in-plane and 0.043 for the out-of-plane fiber orientation components, respectively. The averages of the standard deviations of the in-plane orientations A xx and A zz are roughly 46.9% and 38.3% of the standard deviation of a CUD. Here, the standard deviations scatter strongly. This trend increases for an increased number of cells per evaluation mesh. There is no correlation between cell position and standard deviation. Thus, the standard deviation of individual cells seems to be random as well. The average of the standard deviation of the out-of-plane orientation A yy is approximately 64.1% of the standard deviation of a CUD. In conclusion, there are strong indications that the fiber orientation states in the walls are determined by random events rather than by systematic processes.  The interpretation of these results is that due to a lack of scale separation, each molding represents only a single realization of random processes. Thus, the different initial microstructure in each SMC stack results in varying fiber orientation states in the honeycomb. Of course, the comparison of fiber orientation states via fiber orientation tensors is complicated once FMS leaves only few fiber bundles in the volume under consideration. In conclusion, the evaluation of fiber orientation states and the validation of process simulations via an orientation tensor is questionable if it depends only on a few entities whose existence and orientation are determined by random processes. According to the reasoning above, it is likely that the microstructure of each SMC stack is unique due to the macroscopically heterogeneous fiber bundle architecture of SMC and the multitude of fiber bundles in the stack. Considering the interactions between neighboring fiber bundles, which are determined by random processes due to the lack of scale separation, it may be concluded that the microstructure of each SMC part is unique as well. Following this line of reasoning, fiber matrix separation and fiber orientation would still feature significant scatter, even with an arbitrary number of repetitions.

Fiber Volume Content
This section shows the results of µCT and TGA of the thick coupons molded from an SMC with a nominal FVC of 36.6%. The coupon walls are ordered from I on the left to VI on the right. Subsequently, the simultaneous consideration of µCT and TGA enables the holistic analysis of fiber matrix separation and defects.

Configuration A
In Configuration A, the initial charge is positioned underneath the honeycombs. Thereby, the initial SMC stack overlaps walls I to III. Thus, fiber bundles are transported into the ribbed structures after a short flow path.
The honeycomb shows FMS only in the top sections of some walls with small areas with few or no fiber bundles (see Figure 7). Hereby, FMS develops in different positions with varying severity. The variance in severity and the location of FMS highlights the statistical nature of phenomena such as FMS, that develop in the absence of scale separation. The fibers stay assembled in bundles during compression molding; therefore, individual fiber bundles can be identified in the µCT scan. This matches the findings of other investigators [32,59,60] and confirms the assumption of the bundle approximation used for the Direct Bundle Simulation. Breakage of individual fibers or entire fiber bundles is not observed, which matches the findings of Trauth [66] and Motaghi [67]. The fiber bundles exhibit a strong curvature that results from complex interactions during mold filling. There are small cracks visible in matrix-rich sections. Only few pores can be seen in the junction of two adjacent walls. Three samples from the same process configuration, but different specimens, are analyzed via thermal gravimetric analysis and the resulting average FVC is illustrated in Figure 8. In general, the honeycomb contains FVC that are similar to the nominal 36.6 %. The FMSI equals 0.06 for the results shown in Figure 8. Here, the standard deviation is small compared with the average FVC of the individual sections. The small differences in FVC of adjacent sections could be an indication for locally occurring FMS, which is also visible in µCT scans of coupons molded with Configuration A (Figure 7). The thermal gravimetric analysis confirms that there is slight FMS in walls IV to VI, which do not overlap the area that is covered by the initial charge. The top sections of these ribs is filled later than the rest of the honeycomb. The matrix is likely heated during flow through the honeycomb structure and thereby its viscosity is lowered. Hence, the hydrodynamic force acting on the fiber bundles would decrease, which may facilitate the development of FMS.

Configuration B
In Configuration B, the initial charge is positioned opposite to the honeycombs. Thus, fiber bundles must flow over at least half the length of the base plate to reach the honeycomb structures.
Similar to Configuration A, a complex fiber bundle architecture is observed. Although the bundles have to flow a considerable distance and through complex geometries, fibers remain assembled in bundles during molding. The honeycomb shows pronounced FMS mainly in the downstream walls (see walls I to III in Figure 9). The severity of this FMS and the fiber orientation states vary notably between individual coupons. Again, this indicates that the development of FMS is affected by random processes.  Figure 10, which shows the TGA results for Configuration B, illustrates that the average FVC in the upstream walls is larger than the nominal FVC, while the average FVC in the downstream walls is lower than nominal. The occurrence of FMS in downstream walls most likely arises from fiber bundle jamming in the upstream walls. These walls are closest to the position of the initial charge in this configuration. Therefore, the SMC has to flow through these walls first to reach the downstream walls. The fiber bundles may jam along the flow path, thus hindering further fiber bundle transport and creating matrix-rich sections downstream. The FMSI equals 0.19 for the coupons shown in Figure 10 and is therefore more pronounced in the honeycomb molded with Configuration B than in those molded with Configuration A.  It is also important to note that walls IV and VI both deviate from the direction of flow by the same angle but have significantly different FVC in the top sections. Thereby, wall VI is connected to two other walls, whereas wall IV is only adjacent to wall III. This suggests that neighboring walls limit or enable long fiber transport, depending on their number, position and orientation. Thus, it might not only be the geometry of the rib under consideration that matters, but also the ribs, walls and screw domes attached to it.
TGA and µCT indicate a large scatter of FVC between consecutive moldings in a wall, if the process configuration facilitates the development of FMS in this wall. Similar to the varying fiber orientation states in Section 5.4, this scatter might originate from a stronger dependence of the development of FMS on the initial bundle architecture in each stack. Thus, FVC in coupons molded with Configuration A are more reproducible than in coupons molded with Configuration B. A higher number of repetitions might be helpful to reliably determine the FVC in these configurations, but the question remains if sufficient reproducibility can be achieved even with an arbitrary number of repetitions.

Fiber Bundle Architecture
For a better understanding of the resulting fiber architecture, whole components are pyrolyzed to provide an overview of the fiber architecture in the honeycomb as well as in the base plate. The residue of the complete part after pyrolysis shows that fibers stay assembled in bundles during compression molding, which is in agreement with the µCT results shown above.
In Configuration A (see Figure 11a), the fibers to the left of the honeycombs are mainly oriented in the flow direction, which is likely the result of differently directed forces at the fiber ends. One end of these fiber bundles is pulled inside the honeycomb, while the other end is drawn in the flow direction. This results in the straightening and orientation of fiber bundles in the flow direction. Additionally, the flow front lags behind around the honeycombs because a lot of material flows into the honeycombs themselves. Hence, large velocity gradients emerge at the flow front, which in turn contribute to the reorientation of fiber bundles in the flow direction.
The residue of a honeycomb molded with Configuration B after pyrolysis shows that the peculiarity of FMS increases along the flow path (see Figure 11b). In the upper right cell, FMS is most pronounced and is visible with the bare eye. The dark gray clumps are mineral fillers which do not burn at the temperatures used for incineration. A weld line forms between the right edge of the base plate and the honeycomb from the material that flows around the honeycomb geometry. This shows that weld lines may form as the result of a complex rib geometry and not only due to the presence of multiple SMC stacks.

Numerical Results
The following sections discuss the press data, the fiber volume content and fiber orientation in the thick honeycomb of Configurations A and B predicted by MOLDEX3D, MOLDFLOW, 3D TIMON and the Direct Bundle Simulation. The differences between experimental and numerical results are examined and explained in the context of the respective modeling approaches. Based on these results, the different process simulation methods are assessed.

Configuration A
The simulation of Configuration A with MOLDEX3D predicts that the maximum compression force of 1500 kN is reached almost immediately upon contact with the SMC stack (see Figure 12). This point is marked by a sudden drop of the compression velocity to roughly 0 mm s −1 . The compression force remains constant for the rest of the molding while the compression velocity shoots up to 10 mm s −1 . After a sudden drop, the compression velocity increases slowly up to 10 mm s −1 again and then drops back to 0 mm s −1 when the mold is closed almost completely. Afterwards, the compression force reaches 5 mm s −1 for a short time, during which the mold is closed completely.  The process simulation results from MOLDFLOW show for Configuration A a slight underestimation of the press force until the maximum press force is reached. After the maximum press force is reached, some problems with the press controller occur and the correct press force is not reached exactly. For the predicted closing speed profile, it can be observed that it matches the experimental data well. The combination of the underestimation of the press forces and the underestimation of the velocity peak point could be explained by the too-high compressibility of the used material, which has not been specifically characterized for this material. In comparison with the other simulation results, MOLDFLOW has the longest total compression time of approx. 16 s. The gap prediction is therefore better compared with the other tools, but nevertheless overestimates the total compression time due to some problems with the press controller.
The compression profile of the half geometry molded with Configuration A shows that the compression force predicted by 3D TIMON (15,670 KN) is about a factor of 10.4 higher than the force applied by the press on the full geometry during compression molding (1500 kN). The reason for the discrepancy between experimental and numerical results is the no-slip boundary condition on the mold surface. Usually, the outer layers of SMC are heated by the mold and lubricate the surface with a thin layer of low-viscosity matrix. This alters the flow kinematics and decreases the shear stress in the outer layers of the SMC stack. Thus, the pressure needed to facilitate flow is reduced. Here, the no-slip boundary condition induces high shear stresses. Thereby, the pressure necessary for flow is increased. The jumps in the force profile reflect the discontinuities in the velocity profile, which can only be modeled as a series of constant velocities. This would require an infinitely large acceleration and force acting between each step of the compression profile and is a non-physical representation of the press profile.
In the Direct Bundle Simulation, the compression molding in Configuration A is modeled as the force-controlled closing of the mold with 1500 kN due to the rapid stop of the physical press upon reaching the SMC stack. Here, the compression speed increases rapidly and decreases from there on until the mold is closed. The predicted compression time is significantly shorter than the actual compression time during compression molding due to the high initial closing speed. Employing the force boundary condition skips the initial pre-compression of the stack and saves computational time, but does not account for the rather small volume change in the stack due to the force application at this stage of the real process.

Configuration B
The simulation of Configuration B with MOLDEX3D predicts that the maximum compression force of 1500 kN is reached almost immediately upon contact with the SMC stack (see Figure 13). Similar to the simulation of Configuration A, the compression velocity drops and increases several times before the mold is closed completely.
The process simulation results from MOLDFLOW show for Configuration B a similar behavior to that of Configuration A. The press force is slightly underestimated, until the maximum press force is reached. Nevertheless, the velocity and gap prediction matches the experiments quite well.
The compression profile of the half geometry molded with Configuration B shows that the compression force predicted by 3D TIMON (15,270 KN) is about 10.2 times larger than the force needed by the press for the full geometry during compression molding (1500 kN). This discrepancy and the jumps in the force profile can be explained by the same arguments as for Configuration A.
The Direct Bundle Simulation predicts that the compression force increases rather slowly upon contact with the SMC stack. Here, the velocity profile follows the actual compression velocity data of the press until it switches to force control. The velocity profile correctly shows a sudden decrease in compression velocity once the maximum force of 1500 kN is reached. The predicted gap to closed state matches the press data well. The compression time predicted by the DBS is slightly shorter than the actual compression time.

Configuration A
MOLDEX3D predicts a uniform distribution of fibers in the honeycomb. Here, the FVC is equal to 44.8% (see Figure 14). Thus, the FMSI Moldex3D,A equals 0.22, which is about four times the experimental FMSI (FMSI TGA,A = 0.06). Kuhn found similar results during his investigation of FMS with MOLDEX3D [14]. Here, rib base and rib top contained just a slightly different FVC, which also differed only marginally from the nominal FVC. The use of the FMSI might seem counter-intuitive because the FVC is basically the same, not only in the coupon but in the rest of the geometry as well. Nevertheless, the predicted FVC still deviates from the nominal FVC of 36.6% and the FMSI is employed to quantify the deviation from the nominal FVC.
MOLDFLOW was not considered within the investigations of the fiber volume content, since no model for fiber matrix separation is available. Even if it was not considered here, such a model could be implemented via the Solver API.
The DFS predicts an extreme FMS with low FVC in the top sections and an overall low FVC in walls II and V, which are oriented normal to the flow direction. The FMSI DFS,A is roughly 0.75 which is about a factor of 12.5 larger than what is found by TGA (FMSI TGA,A = 0.06). Here, the result clearly shows the weak point of disregarding fiber-fiber interaction in a directly modeled process simulation. Because fibers do not interact with one another, any number of them may be located at the exact same position. Thus, it is possible to obtain FVC that are higher than the theoretical maximum of a square packing (78.5%) and hexagonal packing (90.7%).
The DBS predicts a pronounced FMS in Configuration A. Here, FMSI DBS,A is approximately 0.38, which is about 6.3 times larger than what is found by TGA (FMSI TGA,A = 0.06). The FVC in the top sections of walls IV, V and VI, which do not overlap the area, where the initial charge is placed, is lower than the nominal 36.6%. This is qualitatively consistent with the results of TGA (see Figure 8). There are large differences between the FVC in directly adjacent walls, for example, the center section of wall III. Here, the FVC is 67.00 ± 5.84 % and therefore larger than the experimentally observed FVC in the rib base or any adjacent section. The FVC decreases from base to top which is expected during fiber bundle jamming.

Configuration B
Similar to Configuration A, MOLDEX3D predicts a uniform FVC in Configuration B which is equal to 44.7 % (see Figure 15). The FMSI Moldex3D,B equals 0.22 and is therefore similar to the FSMI which results from TGA (FMSI TGA,B = 0.19).
The DFS predicts an extreme FMS with low FVC in wall II which is oriented normal to the flow direction. Fibers get stuck in the base of walls I and III, thereby lowering the FVC in wall II. Interestingly, the top section of wall I contains a higher FVC than its base or any of the surrounding sections in wall VI. The FVC in the base of walls I and III is higher than the maximum of a hexagonal packing (φ f,hex,max = 90.7%). The FMSI DFS,B is about 0.60, which is about 3.15 times larger than what is found by TGA (FMSI TGA,B = 0.19).
The DBS predicts a pronounced FMS in Configuration B with FMSI DBS,B = 0.34, which is about 1.8 times larger than what is found by TGA (FMSI TGA,B = 0.19). Here, the walls that are located up the flow path, i.e., walls IV, V and VI, contain a higher than nominal FVC. This is qualitatively consistent with the results of µCT and TGA (see Figures 9 and 10).
However, µCT and TGA show that fiber bundle jamming increases the FVC in the rib base and up the flow path, whereas the FVC in the walls, which are located down the flow path, is decreased. Contrary to the results of µCT and TGA, the DBS predicts that wall II contains a low FVC and that wall I contains a high FVC. Therefore, the DBS predicts that the lowest FVC is found in the wall positioned closest to the end of the flow path. The FVC in the center of wall VI is 60%. In contrast, the maximum FVC found by TGA in one coupon of Configuration B is 51.2%. Therefore, the effect of fiber bundle jamming inside the ribbed structure is overestimated.

Fiber Orientation
The discrete orientation tensors of each cell of the evaluation mesh are rotated by an orthogonal orientation tensor R according tõ such that the coordinate system aligns with the wall orientation. The evaluation of the normal componentÃ nn is shown in Figure 16. This component is of specific interest, because orientation normal to the walls is confined by the molds. It seems unreasonable that a significant portion of fiber bundles is aligned normal to the rib wall given their length of 50 mm and rib width of 2.5 mm. This line of reasoning is confirmed by µCT fiber orientation analysis, which averages toÃ nn = 0.127 ± 0.066 in three samples. As discussed in Section 5.4, a detailed spatially resolved evaluation of fiber orientation tensors at this scale is questionable due to the scatter between individual realizations, but the overall average supports the statement that the out-of-plane component is small. However, both macroscopic models indicate a relevant component in the normal direction. This is a consequence of transporting the initial planar orientation state up to the ribs without sufficient reorientation and no consideration of the confinement by mold walls. The direct simulation models do account for interaction with the mold walls and thus the ability of fibers to orient normal to the ribs is limited. Hence, both direct models predict a reasonable normal component that is non-zero due to fiber curvature and bundles that protrude to different regions of the adjacent geometry. Configuration B leads to an equivalent conclusion and is therefore omitted here. Other orientation components are not compared, because the large scatter does not provide an unambiguous ground truth (see Section 5.4).

Short Shots
In addition to the comparison of press forces and fiber orientations, the flow front during the mold filling of Configuration A is also investigated by so-called short shots. For this purpose, spacers of 1.6 mm and 4.8 mm were inserted between the upper and lower mold, resulting in an incomplete filling. The resulting short shots were then measured three-dimensionally with an ATOS system by GOM to obtain a digital surface model. The digitized short shots with 1.6 mm and 4.8 mm spacers are compared with the predicted flow front at equivalent processing stages by means of a pictorial superimposition of the filling patterns (see Table 1). The short shot of the 4.8 mm experiment shows a uniform filling of the honeycomb walls. In particular, the initial stack's occupancy area (see Figure 1a) is evenly pressed into the honeycombs. A comparison with the simulation results shows that the uniform filling of the honeycomb walls at the rear could only be predicted by the MOLDFLOW (4.8 mm spacer) and the DBS (4.8 mm spacer). The result of the MOLDFLOW calculation show a slight underestimation of the flow front progress, especially in the middle area, while the DBS tends to slightly overestimate the flow front. In general, however, the global filling pattern at this point is most accurately predicted by the DBS. The MOLDEX3D (4.8 mm spacer) and 3D TIMON (4.8 mm spacer) calculations both tend to have a gradient in the calculated flow front from the back to the front wall of the honeycomb. This is particularly evident in the MOLDEX3D result, where the rear area of the honeycomb is already completely filled, while the front area still shows no filling. A similar behavior, although not quite as pronounced, can also be seen in the results of the 3D TIMON simulation. Thus, these two methods do not represent the progress of the flow front sufficiently well at this point.
The short shot of the 1.6 mm experiment shows that the honeycomb geometry is already completely filled, except for one area in the front. This wall has an imperfection in the filling pattern, which can be explained by a fiber or scale separation effect. The fiber bundles, which are significantly longer than the filling structure, interact with each other and accumulate in the lower area of the honeycomb base.
This effect is matched very well by DBS (1.6 mm spacer) as both the location and the expression of the defect match the experiment most accurately at this point. 3D TIMON (1.6 mm spacer) also shows good agreement with the location of the resulting defect, but significantly underestimates it compared with the experiment. MOLDFLOW and MOLDEX3D (1.6 mm spacer) also show incomplete filling of the outer honeycomb wall. However, for both models the defect is predicted on the top of the outer honeycomb wall. In summary, this configuration was predicted very well by all investigated simulation methods within the respective model accuracy.

Compression Profile
The compression profile predicted by MOLDEX3D compares well to the press data of Configuration A; however, there is a noticeable difference between the predicted and actual compression time in Configuration B. It is likely that the lack of differentiation between friction and viscous contributions during the rheological characterization of the SMC leads to a higher determined viscosity. This in turn increases the compression force or time needed for mold filling. Thus, it would be useful to characterize the friction and viscous contributions separately during press rheology.
The press force is quite well predicted by MOLDFLOW in both cases. The press force is slightly underestimated in the beginning, but fits quite well with the start of the compression cycle and the switching point. Furthermore, the predicted gap after the switching point is predicted well, only the problems with the press force controller lead to some misconceptions with the press force and thus the gap height.
As already mentioned in Section 6.1, the compression force predicted by 3D TIMON at the end of compression molding is several times larger than the maximum force. The velocity profile of 3D TIMON roughly approximates the actual compression velocity. Accordingly, the predicted gap to closed state matches the actual press data well. Thus, due to some modeling approaches, which can not be modified in the current version, 3D TIMON does not describe the SMC compression molding material behavior in a suitable way. Therefore, it would be helpful to implement boundary conditions, such as a wall friction model, that model its flow kinematics more accurately.
In Configuration A, the Direct Bundle Simulation overestimates the initial closing speed peak, because the simulation does not account for the pre-compaction of the material prior to the application of the constant compression force. Additionally, the applied compression force is higher than during the trials, since the press did not adhere to the prescribed maximum force during experiments. This results in a significant underestimation of the total compression time. In Configuration B, the compression force predictions of the Direct Bundle Simulation compare well to the press data. The total compression time is less underestimated than in Configuration A.

Fiber Volume Content
As shown in Section 5.5, the variation in the position of the initial charge induces FMS in various sections of the honeycomb. The uncertainty measure of FVC is larger, where severe FMS develops. The mechanism behind FMS is most likely fiber bundle jamming in walls along the flow path and the base of the walls. This matches the findings of Kuhn [14] who showed by experiment and simulation that fiber jamming is the driving mechanism behind FMS in ribbed structures made from long-fiber-reinforced polymers.
The macroscopic simulation tool MOLDEX3D predicts the FVC to be basically constant throughout the honeycomb. There is no distinctive distribution of fiber bundles inside the coupon as a result of mold filling. Thus, the fiber bundle jamming and the following development of matrix-rich sections is not replicated. Therefore, the statistical description of fiber-matrix interaction and fiber-mold interaction in macroscopically heterogeneous fiber suspensions might not be useful to predict FMS. The model used for the calculation of fiber volume content is based on a particle migration model introduced by Morris and Boulay [68]. It is a suspension balance equation which describes the flow of rigid spherical particles. It is questionable how such a model could describe the complex interactions between individual fibers or fiber bundles as well as between fiber bundles and the mold. Tseng et al. [69] validated the suspension balance model by comparing their prediction with the experimental work of Vélez-García et al. who investigated the FMS in a center-gated disk molded by injection molding [70,71]. However, the material used was a short-fiber-reinforced polymer and the geometry had no curvature, ribs or bosses. Therefore, the model used by MOLDEX3D might not be suitable for the simulation of fiber matrix separation during the compression molding of ribbed structures molded from long-fiber-reinforced polymers.
In contrast to MOLDEX3D, the Direct Fiber Simulation predicts a distinctive distribution of fiber bundles inside the coupon. Here, the fibers accumulate in the rib base and the FVC predicted by DFS are larger than the maximum FVC of a hexagonal packing. The key to answering the question of why there are high FVC in some walls is the DFS algorithm. As already stated, the DFS algorithm checks whether fibers are stretched above a certain threshold. Once the threshold is exceeded, the positions of the fiber nodes are recalculated. This is similar to a setback and fiber motion is impeded. Thus, fibers become stuck or are slowed down in flow situations in which the setback is triggered continually. One example of such a situation is complex ribbed structures, such as the honeycomb. Here, the setback is triggered because the fibers are pulled inside the honeycomb and at the same time in the flow direction.
The DBS shows a distinctive distribution of fiber bundles and the FVC predicted by the DBS is generally larger than what is found by TGA. The discrepancy between the experimental and predicted FVC in Configuration A during the pull-out of bundles might occur due to the modeling of the matrix as a Newtonian fluid, the viscosity of which does not depend on temperature. The DBS predicts in Configuration B that the fiber bundle jamming results in low fiber volume contents in wall II, which is positioned closest to the end of the flow path. However, the results of µCT and TGA show that wall I contains significantly less fiber bundles than wall II, although this measurement is accompanied by a large standard deviation.

Fiber Orientation
As already mentioned in Section 5.4, the fiber orientation states in honeycomb structures compression-molded from SMC feature significant scatter. Still, a statement can be made about the meaningfulness of the predicted fiber orientation components. Since the fiber bundles are 50 mm long and the wall only about 2.5 mm thick at the top, the fiber orientation component in the thickness direction should be relatively small. Thus, the fiber orientation should be dominated by the in-plane fiber orientation components.
MOLDEX3D and MOLDFLOW predict normal orientation components that are significantly larger than what is predicted by the directly modeled process simulations. This highlights that tensor-based fiber orientation models are not suited to simulate the orientation kinetics in fiber suspensions where there is a lack of scale separation. This could be addressed in the future by the incorporation of constraints in tensor-based models, for example, by following the work of [72]. In the DFS and DBS, fiber re-orientation depends on the modeled physics and not on the choice of fiber orientation model, closure approximation and interaction coefficient. Hence, the use of directly modeled process simulations is advantageous for the simulation of fiber orientation states in SMC components with confined geometric regions.

Short Shots
The results from the short shots in the 4.8 mm case show the influence of the wall friction used in the simulation model. The results of models, which use either no friction or a no-slip boundary condition, show a linear progression of the filling from the back wall of the honeycomb to the front. The models using a modified wall friction model (see Equations (3) and (6)) to represent the lubricating film between the material and the wall during manufacture show a more uniform filling of the honeycomb cavity, which provides a more realistic representation of the experimental results. In general, the DBS shows the best match of the filling pattern for the two investigated short shots. This underlines the necessity to accurately model the basic physical effects, such as fiber interactions or friction effects, when calculating the flow of highly filled long-fiber-reinforced polymers.

Outlay and Benefit
The characterization of thermo-rheological material parameters for the process simulation of long-fiber-reinforced materials is an elaborate and time-consuming process. Therefore, it is advantageous to employ simulation methods that produce acceptable results with a minimal number of material parameters. It is necessary to characterize the SMC viscosity for MOLDEX3D, MOLDFLOW and the DFS. In contrast, the Direct Bundle Simulation needs only basic information about the initial FVC, fiber bundle geometry and matrix viscosity to run. However, there usually is no information about the rheological properties of the matrix of commercially available SMC materials. Furthermore, samples of the matrix material might not be provided by the SMC prepreg producer.
Unfortunately, the simulation run-times cannot be compared rigorously, because the computations are performed on different computers employing varying element types and spatial discretizations. Nonetheless, in general, the calculation run-times of directly modeled process simulations are usually longer than those of tensor-based models because of the additional calculation steps during the fiber simulation. Here, the DBS requires more computational resources than the DFS model due to the additional consideration of fiber interactions. Hence, tensor-based statistical models deliver results in a shorter time than directly modeled process simulations.

Conclusions
The predictive capabilities of different simulation approaches for the compression molding process of sheet molding compounds are evaluated and compared with experimental molding trials of a honeycomb structure with different process configurations. The assessment considers the compression profile as well as the fiber volume content and fiber orientation inside a confined region.
The development and severity of fiber matrix separation during the compression molding of sheet molding compounds is investigated experimentally by thermal gravimetric analysis and micro-computed tomography. Long flow paths facilitate fiber bundle jamming and the development of large matrix-rich sections. The fiber matrix separation index proved useful for the quantification and comparison of FMS in different configurations, as well as for the comparison between experimental and numerical results. Due to the macroscopically heterogeneous structure of sheet molding compounds, the location and severity of FMS and fiber orientation vary. Particularly in confined regions, the prediction quality of tensor-based models is considerably reduced, since fiber-fiber interactions, fiber-matrix interactions and fiber-mold interactions cannot be modeled explicitly.
Macroscopic process simulation tools are useful for a quick estimation of the mold filling pattern of ribbed structures. Combined with adequate material data, the compression profile can be predicted acceptably. Thus, macroscopic process simulations are valuable during the design phase. However, if the user is interested in the prediction of the fiber architecture in confined regions, a directly modeled process simulation should be preferred. They enable the simulation of local changes in fiber volume content. Additionally, directly modeled process simulation considers the constraint that is imposed on long fibers in confined regions by the mold. This leads to a more realistic out-of-plane fiber orientation component. However, the direct simulation of fibers requires more computational resources and calculation run-time than macroscopic process simulations. Acknowledgments: We would like to thank TORAY ENGINEERING CO., LTD. and CORETECH SYSTEM CO., LTD. for providing their software tools for research and testing purposes. We want to thank Apis mellifera for the groundbreaking work on efficient space filling with regular polygons.

Conflicts of Interest:
The authors declare no conflict of interest.   Figure A1. Average of the viscosity η (solid lines) of three press rheometry samples vs. shear rateγ of the SMC at T = 25°C and T = 70°C. The dashed lines indicate the viscosity function described by the Cross-Andrade model and the fit parameters shown in Table A2.