Influence of fiber alignment on pseudoductility and microcracking in a cementitious carbon fiber composite material

This research examines the effect of fiber alignment on the performance of an exceptionally tough 3D-printable short carbon fiber reinforced cementitious composite material, the flexural strength of which can exceed 100 N/mm2. The material shows pseudoductility caused by strain-hardening and microcracking. An extrusion-based manufacturing process allows accurate control over the spatial alignment of the fibers’ orientation, since extrusion through a tight nozzle leads to nearly unidirectional alignment of the fibers with respect to the directional movement of the nozzle. Specimens were investigated using mechanical tests (flexural and tensile load), augmented by non-destructive methods such as X-ray 3D computed tomography and acoustic emission analysis to gain insight into the microstructure. Additionally, digital image correlation is used to visualize the microcracking process. X-ray CT confirms that about 70% of fibers show less than 10° deviation from the extrusion direction. Systematic variations of the fiber alignment with respect to the direction of tensile load show that carbon fibers enhance the flexural strength of the test specimens as long as their alignment angle does not deviate by more than 20° from the direction of the acting tensile stress. Acoustic emission analysis is capable of evaluating the spatiotemporal degradation behavior during loading and shows consistent results with the microstructural damage observed in CT scans. The strong connection of fiber alignment and flexural strength ties into a change from ductile to brittle failure caused by degradation on a microstructural level, as seen by complementary results acquired from the aforementioned methods of investigation.


Introduction
Unreinforced cementitious materials are well known for possessing high compressive strength while being weak under flexural and tensile loads. Such loads bring about a very sudden brittle failure, meaning a single macroscopic crack leading to instant destruction of the loaded specimen. Common technological solutions to mitigate poor behavior under tensile load and suddenness of the failure include the introduction of reinforcing steel bars or alternatively various kinds of fibers. Most reinforcement fibers are made of minerals, glass, polymer, steel or carbon [1][2][3][4][5]. Fiber reinforcement does not necessarily lead to improvement of tensile strength, as the fibers' mechanical properties might be unsuitable; they may be below the critical fiber length or the bond between fiber and matrix might not be strong enough to support full transition of loads from one material to another, which has a significant impact on the macroscopic response [6]. Such composites will show a slightly improved toughness, as the fibers are pulled out after cracks started separating the matrix [4,7]. Generally, such systems are labeled as 'Fiber Reinforced Concrete (FRC)' in the literature [8]. However, tailoring chemical compatibility and micromechanical parameters of cementitious materials towards fiber reinforcement allows for the creation of a material that no longer shows brittle single crack failure but rather a significant pseudoductility. This behavior shows itself in a strain-hardening response beyond purely elastic deformation capabilities. Such composites also show enhanced tensile and flexural strength. These materials have been dubbed 'Engineered Cementitious Composites (ECC)' or 'Strain-Hardening Cementitious Composites (SHCC)' in literature. [9,10]. The pseudoductile behavior is known to be attributed to a change of the microstructure. At stresses beyond the linear-elastic phase, small-scale cavities appear. Instead of combining into a macroscopic crack and critically failing instantly (brittle fracture) or a slow but continuous growth of the cavities leading to a decrease in stress (strain-softening), two effects combine. (a) The presence of the fibers stops the growth of the cavities and (b) the stiff fibers begin to transmit the stresses, which leads to an overall increase in the stress-strain response (strain-hardening). Thus, once the material undergoes strain-hardening, a large number of tightly spaced microcracks appear and those cracks are bridged by the fibers. Therefore, such materials are not only characterized by an increase in tensile and flexural strength but also by a considerably higher capacity for deformation. [11][12][13]. Typically these composites are made using polymer fibers, i.e. polyvinylalcohol (PVA) and high-density polyethylene (HDPE) fibers in particular [14]. Recent research has shown that addition of stiffer and stronger fibers, such as carbon fibers, coupled with the usage of an extrusion-based alignment process, a cementitious composite, which exhibits a flexural strength in excess of 100 N/mm 2 , can be created [15]. Such materials are highly interesting for applications in additive manufacturing ('3D printing') [16]. The extreme flexural strength makes the composite an enticing material to reduce the amount of continuous steel reinforcement necessary, as some tensile load can be handled by the cementitious composite itself. Opposed to reinforcing steel, carbon fibers show no tendency to corrode or degrade under atmospheric conditions, leading to a higher expected lifetime of the final structure. The low crack widths also lead to other desirable traits, such as high durability by hindered diffusion of aggressive compounds into the matrix and self-healing properties [17][18][19][20]. Reinforcement with carbon fibers grants the material a degree of electrical conductivity, which allows the material to be potentially used in heating elements [21] or in the field of structural health monitoring [22].
The aforementioned extrusion process leads to a nearly unidirectional alignment of the fibers. As the fibers are pushed through a tight nozzle, their alignment will follow the movement direction of the nozzle. The effectiveness of this principle has been shown in prior research on Portland cement [15,16], aluminate cement and mixed cementitious binders [23] and also with polymer resins [24][25][26]. The material has been successfully used to create 3D printed samples with high flexural strength [16,27].
While the first serious research in the field of additive manufacturing in the construction sector can be traced back to the 1990s [28,29], such methods had problems scaling up to structural levels and penetrating the market up until recently [30][31][32]. Especially the fact that reinforcing materials still have to be manually inserted into the structures and cannot be automatically incorporated into the printing process is a problem that has only recently been acknowledged by research [33]. As the process described in this paper is analogous to the 3D-printing method commonly described as 'Liquid Deposition', in which strands of a flowable material are deposited layer by layer to create a desired shape after hardening, this poses an elegant solution to this problem, as no external reinforcement is needed and a fully automatic process is realized. As fiber alignment follows printing direction, the production of structural members specifically designed to withstand complex loading scenarios by matching up printer movement to occurring principle tensile stress paths should be feasible.
One specific problem that occurs when using inherently anisotropic materials such as fibers is the fact that exact knowledge and control of fiber alignment plays a huge role in the interpretation of experimental results. As traditional mold-casting techniques do not allow for customization in this regard, little information between the direct correlation of fiber alignment and material behavior can be found in actual data. Applying non-destructive testing techniques is ambitious due to the high degree of heterogeneity and anisotropy [34]. This research focuses on the evaluation of fiber-matrix interaction by correlating macroscopic methods such as mechanical testing with highly localized non-destructive methods of analysis, such as digital image correlation (DIC), ultrasound (US), acoustic emission analysis (AEA) and X-ray 3D computed tomography (X-ray CT).

Overview and scope
This publication aims to shed light on the correlation between the general response of a carbon-fiber reinforced composite material during mechanical testing (e.g. ultimate strength and capacity for deformation under flexural and tensile load) and the mechanism of the resulting damage, in particular the relation between in-plane fiber alignment and the resulting microcracking behavior.
The feasibility of measuring fiber orientation and damage indicators (especially in-volume growth of microcracks) using high resolution X-ray CT coupled with advanced segmentation techniques on miniature bending beams is tested and the data is correlated to macroscopic mechanical behavior. We intend to show the reasons behind a transition from ductile to brittle behavior in the material as fiber alignment starts to deviate from the direction of the acting tensile load.
Alongside the CT, other powerful techniques are used to obtain a more complete picture of the fracture process. In particular DIC for in-situ observations of surface crack growth (using larger dogbone-shaped samples) [35] and acoustic emission analysis (using miniature bending beams) for further analysis of failure mechanism by means of frequency analysis are applied. Setup of miniature beams during 3-point bending tests and definitions of the coordinate system and angles used can be seen in Fig. 1.

Sample preparation of miniature beams
All test specimens were prepared as described in a previous publication of our working group [15], which describes a formulation tailored to allow for superior fiber dispersion and bond strength between fiber and matrix The cement used was a CEM I 52.5 R type supplied by Schwenk Zement KG, sourced from the cement plant Karlstadt. Carbon fibers were purchased form Teijin Ltd. being commercialized under the product name Tenax-J HT C261. The carbon fibers, which are uniformly chopped down to a length of 3 mm, have the following specifications: diameter: 7 lm, tensile strength: 4000 N/mm 2 , Young's modulus: 238 GPa. Prior to usage in the cement paste, the fibers were oxidatively heat treated in an open furnace at 425°C for 2 h to remove the sizing and to oxidize the carbon fiber's surface, improving wettability and fiber-matrix bonding [15,36,37]. Barium sulfate is added to increase the X-ray absorption of the cementitious matrix and thus to increase the electron density contrast with respect to other embedded materials, especially the carbon fibers [38]. All components used for sample preparation are given in Table 1.
Cement (CEM I 52.5 R, Schwenk Zement KG, Ulm, Germany), silica fume (EFACO, Egyptian Ferro-Alloys Company, Edfu, Egypt) and barium sulfate (Acros Organics, Geel, Belgium) are added into the mixing container and dry mixed by hand. Deionized water and superplasticizer (MasterGlenium ACE 430, BASF Construction Solutions GmbH, Trostberg, Germany) are premixed in a glass beaker. The liquid components are mixed into the binder using a Heidolph Hei-TORQUE Precision 400 overhead mixer (Heidolph Instruments GmbH & CO. KG, Schwabach, Germany) at 400 RPM for 90 s until all solids are dispersed homogenously. After scraping undispersed remains off the stirrer and container walls, the paste is mixed again for 90 s at 2000 RPM. After this step, the carbon fibers (Tenax-J HT C261, 3 mm in length, 7 lm in diameter, Teijin Ltd., Tokyo, Japan) are admixed into the paste and stirred at 70 RPM for 30 s. The fiber-reinforced paste is extruded through a disposable 10 mL syringe with a nozzle diameter of 2 mm (B. Braun Melsungen AG, Melsungen, Germany). Fibers were found to stay mostly intact throughout the process with only very little reduction in length visible. Some further information is included in the Supplementary Material to this paper. Admixing the fibers leads to a very noticeable increase in viscosity of the paste. Flow table tests using a Hägermann table (cone diameter of 100 mm) show a flow value of 181 mm for the paste without fibers, 152 mm for the paste containing 1 vol% of fibers and 117 mm for the paste containing 3 vol% of fibers after 15 shocks.
This process leads to an almost unidirectional fiber alignment along the movement path of the nozzle; the fiber orientation within representative volume elements of the test specimens are quantified by X-ray CT measurements. To obtain samples with a specific inplane fiber alignment, guide templates assisting the extrusion of the cement paste along predefined lines were printed and put under a transparent piece of polycarbonate plastic. The syringe was guided along the lines of the template, producing a slab measuring approximately 100 9 50 9 3 mm. The process can be seen in Fig. 2.
The slabs were allowed to harden for 24 h at 100% relative humidity by storing it within a sealed desiccator over water. After this period, it was stored under Fig. 1 Setup of sample and description of angles during testing and alignment analysis. a Coordinate system with a hypothetical angled object. The angle between x-and y-axis describes the inplane alignment (called h-angle or yaw). The angle coming from the positive z-axis and the object to be described is the out-ofplane alignment (called u-angle or pitch). b Orientation of bending beams within the coordinate system. Beams are aligned along their length with the x-axis and along their width with the y-axis. b Subsequently, fibers that lie in direction of the tensile force occurring during bending are labeled as having a h-angle of 0°. c u-angles describe the upward and downward tilt, with an angle of 90°being fully parallel to the x-y-plane Volumetric dosage of fibers is based on the calculated total volume of the fresh paste water for another 6 days. Following this, the sample was allowed to dry and stored within a sealed desiccator at 59% relative humidity for 21 days. The humidity level was adjusted by maintaining a reservoir of oversaturated sodium bromide solution within the desiccator. All tests were carried out after 28 days of storage. After 7 days of storage, the samples were manually lapped down to a thickness of approximately 2 mm using tungsten carbide powder. This procedure also made sure that remaining roughness of the sample surface was eliminated. For 3-point bending tests the slabs were then cut into miniature beam-shaped specimens of the dimensions 60 9 4 9 2 mm using a Buehler IsoMet low speed saw (Buehler, Lake Bluff, IL, USA).

3-Point bending tests on miniature beams
The flexural strength of miniature beams was tested using a ZwickRoell zwickiLine Z5.0 universal testing machine with a 5 kN load cell attached. Testing was carried out in deformation controlled mode at 1 mm/ min. The force (F) measured by the machine is converted into flexural stress (f ct,fl ) using the following equation: where l is the span between the supports (set up at 50.0 ± 0.1 mm), w and h are the width and height of the sample, respectively. Dimensions of the sample were determined using a digital caliper prior to testing. Unless specified otherwise, strain (e) was calculated from deflection (D) at the center of the miniature beam according to the following equation: where h is the height of the sample and l is the span between the supports (set up at 50.0 ± 0.1 mm).

Visualization of crack patterns using digital image correlation
Since the small size of the miniature bending beams severely limits the usage of more sensitive deformation measurement tools such as strain gauges, larger test specimens are produced according to the process described in Sect. 2.2 using a 3D-concrete-printer developed at the Technical University of Munich. The tests are carried out under uniaxial tensile load using dogbone-shaped samples with a carbon fiber content of 1% by volume. They measure 450 mm in length, the cross section at the midpoint of the specimen measures 50 9 50 mm, the bone then expands towards both ends to a cross section of 100 9 50 mm (Fig. 3a). The specimens are tested with clamped ends in order to minimize eccentric loading caused by crack development and displacement of the center line (Fig. 3b). To achieve this, the dogbone is glued with duromer adhesive MC-DUR 1280 (MC-Bauchemie Müller GmbH & Co. KG, Bottrop, Germany) into bearing plates connected to the testing machine.
As the degradation process of the material beyond the linear-elastic region is strongly characterized by small and very localized deformation processes, strain gauges and DIC are used to provide a more accurate description of sample deformation. Especially DIC is a b Guiding a syringe along the path of an angled template, samples with specific alignment can be created. c Optical microscope image of aligned fibers in the crack of a test specimen after failure powerful tool to visualize the formation and growth of cracks across a whole specimen.
The dogbone-shaped specimens were tested using a SCHENK 400 universal testing machine with a 400 kN load cell attached. The testing program was deformation controlled at 0.18 mm/min. The maximum force (F) measured by the load cell is converted into tensile strength (f ct ) by normalizing force to the load bearing area (A) of the cross section according to: Deformation was measured using strain gauges FLAB-3-11 (Tokyo Measuring Instruments Laboratory Co., Ltd., Tokyo, Japan) and a GOM Aramis stereocamera (GOM GmbH, Braunschweig, Germany).
Additionally, flexural strength is tested on beamshaped specimens in a 3-point-bending setup. The dimensions are an elongated version of the standardized beam referenced in DIN EN 196-1. They measure 40 9 40 9 220 mm, the effective span is 180 mm and the load is applied at the midpoint. The standardized beam specimens were prepared with a carbon fiber content of 2% by volume, as the simpler geometry allowed for printing with a higher fiber content. All specimens were produced according to the mixture design given in Table 1.
Using photogrammetry combined with DIC, local strain development on the surface can be made visible and crack widths can be measured. Even cracks not visible to the naked eye can be visualized without further time consuming preparation (e.g. dye penetration) [39]. For this reason, DIC has been widely used to characterize fracture processes in construction materials [40,41]. The specimens are prepared by applying white paint on the surface and spraying a random dot pattern with a black graphite lacquer on top. During testing, the prepared area is photographed continuously. The first picture taken is used as reference in an unloaded stage to which the following pictures are correlated. To track displacements, the area is divided into small facets and due to the random dot pattern, each facet has its own fingerprint and can be found in every picture. The movement of the facets is used to calculate real displacements and strains.
2.5 X-ray 3D computed tomography (X-ray CT) measurements X-ray CT measurements were performed on miniature beam samples using a Phoenix nanotom m scanner (GE Inspection Technologies LP, Lewistown, PA, USA). One scan comprised of 2000 single images taken with an integration time of 2000 ms. The X-ray was created using an acceleration voltage of 70 kV and a beam current of 190 lA. The voxel size for all scans was 1.9 lm. To obtain as much information about the fracture process zone as possible, each sample was scanned four times at different positions along its length. Reconstruction was carried out for each partial scan using the software Phoenix datos|x. The reconstruction process was calibrated using the datos|x module 'agc' ('automatic geometry calibration') and a ring-artifact reducing filter and a Gaussian smoothing filter were applied. The reconstructed scan was saved into.raw format and imported into the metrology software Dragonfly Pro 4.1 (Object Research Systems (ORS) Inc., Montreal, Canada) for analysis. The resolution of the scan was scaled down by a factor of 2 from the original voxel size of 1.9 lm to 3.8 lm during the import process to decrease computational demand There, the reconstructed multiscans were aligned using the stitching functions of ORS Dragonfly in the overlapping regions.

X-ray CT segmentation procedure
The main challenge during the segmentation procedure is the fact that the contrast between carbon fibers and cementitious matrix is rather low and the greyscale values of each material overlap instead of being separated into singular peaks in the histogram. Thus, a segmentation procedure, which is purely based on threshold values for each material, is not feasible. Instead, the scan is filtered through a local equalization algorithm, which emphasizes local deviation of greyscale values. As there is very little deviation within the fibers themselves, they will appear as black lines, within a bright (albeit noisy) matrix. The black portion of the locally equalized picture can then be segmented by manually defining a threshold value. Remaining artifacts and noise can be removed by performing a connectivity analysis and transforming connected voxels into discrete objects. By deleting objects that do not fit the geometric parameters of a fiber from the segmentation (e.g. spherical objects with an aspect ratio approaching 1 and objects which deviate too far from the fibers geometrical volume), a good fit of segmented voxels to the fiber material in the scan can be achieved. Segmentation strategies based on similar principles are often found in the literature when purely grey value-based segmentation fails. [42,43]. Microcracks were also segmented from images, which were passed through the local equalization filter. Since cracks are relatively discrete objects, they were not segmented using a global greyscale threshold, but rather by using a region-filling tool, which adds a connected region of the same greyscale value to the current segmentation. Objects that were falsely identified as belonging to a crack (e.g. fibers bridging through the crack) were removed by skeletonizing and then dilating the segmented region.

Acoustic emission analysis
Acoustic emission testing on small specimens made of various fiber reinforced materials was conducted by several working groups in the past. These experiments include a variety of different materials such as SiC reinforced ceramics [44], glass-fiber reinforced poylmers (GFRP) [45] and steel-fiber reinforced concrete [46]. The majority of publications consider carbonfiber reinforced polymers (CFRP) and attempt the determination of failure mechanics [47,48]. To achieve this, source localization frequency analysis and event counting methods are applied. Within this study, acoustic emission testing is used to monitor the damage evolution in miniature beam specimens with a carbon fiber content of 1 vol% under static three-point bending load [49]. The objective was the spatial and temporal resolution of crack formation and fiber failure. Four small-sized resonant acoustic emission (AE) sensors were attached to the specimen using hot adhesive. Due to the resonant frequency response of the used AE sensors and the requirement to perform a frequency analysis, waveforms were recorded by a Laser-Doppler-Vibrometer (LDV) in addition. The used LDV provided a quasilinear frequency response function between 30 kHz and 24 MHz and a resolution of 0.1 pm/Hz 0.5 . Load was applied using a Wolpert 10 kN load frame and a 500 N load cell. Acoustic waves, emerging from fiber breakage or matrix fracture due to mechanical load and propagating through the specimen were detected by the AE sensors and the LDV. The induced voltage signal was then transferred to an Elsys eight channel data acquisition system, where it was digitized at 20 MHz sampling rate. In Fig. 4 an overview of the experimental setup is provided. Data analysis was done with own code based on MATLABÒ.
In order to obtain the source coordinates of the acoustic emissions, the crucial step is the picking of the onset. Since there are several thousand acoustic emission events recorded for each experiment, automated onset picking methods must be applied. In this study, a picking method based on the Akaike Information Criterion (AIC) is used [50]. A first guess for the onset is determined by an adaptive threshold operating on the Hilbert-Transform of the signal, as described in [51]. Provided correct onset times at each sensor, it is possible to calculate event location x e , phase propagation velocity c and origin time t 0 by solving the equation system given in Eq. (4) in a leastsquare sense.
Ajb ¼ S xi denotes the sensor position in x direction and t ai is the arrival time at sensor i. The unknown parameter a 1;1 correspond to t 0 ,a 1;2 to x e = c , and a 1;2 equals to 1 = c Ã S x2 . Thus for every event we obtain the propagation velocity of the respective wave type. Together with the AIC-curvature value (introduced as 'DD value' in the literature [50]) the events can be classified into real transient signals and noise or reverberations. The DD value represents the second derivative of the AIC function in the vicinity of the global minimum. Furthermore, a quality estimation of the onset picking process can be derived. Figure 5 shows the relation between event velocity and DD value within a two-dimensional histogram. Events with a DD value above 380 and velocity ranging between 1000 and 5000 m/s, are highlighted by a blue background and are considered for further processing. Low velocity events result either from secondary arriving phases, reflections or picking errors. To ensure consistent localization, the velocity for three and four arrival times is calculated. Deviations between those two solutions larger than a 500 m/s are considered as false picks rather than secondary phases. The threshold is derived from stochastic input into the equation system while solving for corrupted synthetic arrival times. When considering the sampling frequency, it was found that cumulative picking errors of 0.5 ls or 10 samples lead to velocity errors of 1000 m/s. Therefore, 500 m/s is considered as reasonable to avoid false localizations. Following successful source localization, source mechanics can be examined. Research conducted by different authors as indicated in Table 2 (see also e.g. [52]) suggests that the frequency content of acoustic emissions are suitable for failure mode identification. They performed acoustic emission measurements on fiber reinforced materials during different loading regimes and derived relations between frequency and failure modes. Table 2 shows these frequency-to-failuremode relations from preceding research. The correlations described are supplemented with details on the used sensor equipment and the tested materials. The major influence on the found frequencies is the used sensor technology and the evaluation method.

Fiber alignment
Fiber alignment was measured using pristine miniature bending beams, which were later subjected to bending load. The damaged samples are then scanned again to analyze crack patterns. As X-ray CT scans are reconstructed into a 3-dimensional image, it is possible to analyze the alignment of both in-plane fiber angle (h-angle) and out-of-plane angle (u-angle). Figure 1 shows how the sample is oriented for alignment analysis and how description of those angles is handled. The in-plane alignment (h-angle) is given as angular deviation from the direction of the occurring tensile force during mechanical testing. Fibers, which are perfectly aligned to the tensile force are given an alignment angle of 0°.
The fibers are segmented according to the procedure laid out in Sect. 2.6. Once a satisfactory segmentation is achieved, the alignment angles can be directly exported from ORS Dragonfly using the 'connected components analysis' tool. Excerpts from the tomographic scans showing the overlap between segmented fibers and original scans are provided in the Supplementary Material to this paper. Overall, the  segmented voxels overlap very well with the fibers visible in the X-ray scan. Some oversegmentation can be seen in the form of segmented pores or non-fibrous grains, which might lead to underestimation of the overall alignment, i.e. the segmentation data may artificially suggest a slightly broadened alignment distribution or outliers that do not reflect the actual fibers. Figure 6 shows the distribution of measured fiber alignment in each template. Figure 6a shows the inplane alignment (h-angle), while Fig. 6b shows the out-of-plane alignment (u-angle). In both cases, the red area of the probability-density distributions represents the 70% quantile around the mode value, while the grey tail-ends represent the 15% highest and lowest values measured. In general, the distributions in Fig. 6a show good correlation between the expected orientation angle defined by the guide lines of the templates and the measured in-plane orientation values. Figure 6b shows that the out-of-plane alignment occurs mostly parallel to the x-y-plane. The fibers show a u-angle alignment of mostly 90°, i.e. they show very little upward or downward tilt. In both cases, the 70% quantile suggests that fiber orientation of both angles mainly deviates in the range of ± 10°, which is in good agreement with prior research on this system [15]. Deviations larger than 20°are very rare for both angles, even in samples with a somewhat broader distribution (e.g. the 40°sample).

Flexural strength
When referencing alignment angles during strength testing, the nomenclature for the in-plane alignment laid out in Fig. 1 is used. A fiber with the alignment angle of 0°is oriented in the same direction as the Fig. 6 Distributions of fiber alignment as laid out in Fig. 1. Distribution of a in-plane fiber alignment (h-angle) and b out-ofplane fiber alignment (u-angle). The filled in red area around mode value of the distribution represents a 70% quantile of wellaligned fibers, while the grey tail-ends represent the 15% lowest and highest values respectively tensile stress, a fiber with the orientation of 90°is oriented perpendicular to it. Results of the 3-point bending tests can be seen in Fig. 6. Figure 6a shows the stress-strain diagrams of samples containing 1 vol% of carbon fibers, Fig. 6b shows the stressstrain diagrams of samples containing 3 vol% of carbon fibers.
The relation between ultimate flexural strength and fiber alignment angle can be seen in Fig. 6c As can clearly be seen in the experiments, not only does ultimate flexural strength decrease as the fiber alignment angle shifts away from the direction of the tensile load, but strain at break (and with it the capacity for strain-hardening behavior) shows increasingly lower values as well. While samples at 0°alignment angle can show around 0.8% flexural strain before breaking, this value is already halved to values around 0.4% at angles of 20°(see Fig. 7a, b). Samples with an alignment angle of 40°show a strain at break of around 0.1% or less, which is similar to an unreinforced sample. In contrast to samples without fiber reinforcement, some post cracking behavior (most likely linked to fiber pull-out) can be seen.
The angular dependence of the stress-strain behavior can be quantified theoretically by taking into account that the stress-strain behavior is roughly piecewise linear (cf. Fig. 6a). Starting from an unloaded state, with increasing load, the materials behaves linear-elastically until microcracks begin to appear. When the strain is increased further, singlecrack propagation is substantially decreased and the stress-strain relationship behaves roughly linear again until just before failure. Assuming the fiber-volume concentration being small and the fibers being uniformly aligned, the stiffness tensor can be approximated by homogenization techniques or-under further assumptions-by the Mori-Tanaka technique [57] employing the Eshelby tensor [58], where the short fibers are approximated as prolate ellipsoids, cf. [59] for a recent discussion in the context of short-fiber composites. In the strain-hardening phase, numerical homogenization techniques based on representative volume elements containing fibers and microcracks can be used to compute an approximate stiffness tensor.
If the stiffness tensor for a 0°-orientation in Voigt notation is given by C, then the corresponding stiffness tensor describing the same material but rotated by the angle h in the x-y-plane can be computed as with transformation matrix T h and the transposed via simple (but technically involved) coordinatetransformation arguments and by exploiting symmetry properties [60,61]. Note that this transformation assumes a rotation of the entire material (i.e. representative volume element). This is likely a very good approximation in the undamaged case but less so in the case when microcracks have appeared as their orientation will depend on the direction of the applied stress, which is not rotated with the material. Furthermore, in the case of an ideal fiber reinforced material with fiber alignment along the x-axis (cf. Fig. 1b), the material can be approximated as transversely isotropic. The corresponding stiffness tensor C is described by only five independent material constants (instead of 21 independent constants in a fully anisotropic material): with longitudinal Young's modulus E L , transverse Young's modulus E T , m 12 Poisson's ratio for loading along the alignment axis x, m 23 Poisson's ratio in the plane perpendicular to the preferred direction x, inplane shear modulus G T and shear modulus in planes parallel to the preferred direction G L , and m 23 ¼ E T 2G T À 1 due to the isotropy assumption. This observation vastly reduces the measurement efforts when fully characterizing the undamaged reinforced concrete and possibly even the damaged material if highly symmetric crack surfaces can be considered a reasonably approximation.
3.3 Visualization of crack patterns on dogbone specimens using digital image correlation As described before, the plateau in the stress-strain diagram is caused by multi-microcracking over a large area of the specimen (cf. [9, 11,13]). In sharp contrast to other concrete structures, which show brittle failure and develop a single critical crack, the bridging action of the fibers leads to a state where the formation of a new crack is energetically preferred over the growth of an already existing crack [13]. DIC can serve as a powerful tool in the localized description of such phenomena. Figure 8 shows the behavior of a dogbone specimen under uniaxial tension. The first cracks appear after leaving the linear-elastic phase (point a in Fig. 8a, cracking pattern in image a in Fig. 8b), which then slowly expand and grow until failure, causing further strain in the dogbone while stress only increases minimally (points b-d in Fig. 8a, cracking patterns in the respective images in Fig. 8b). Once the strain-hardening effect has fully set in (points b-d in Fig. 8a), crack spacing is stable at values of around 10-15 mm. The stress-strain-relation during tensile testing can be described as bilinear with huge strain-atbreak values caused by a capacity for plastic deformation not present in unreinforced specimens. Unlike the behavior during flexural tests, the first microcracks appear when the stress is already very close to the ultimate strength. Still, there is a large amount of deformation capacity present in the material. These deformations result from elastic strain of the compound matrix and the growth of many microcracks, which further elongate the specimen. The average crack width measured from the image data is about 15 lm, while the maximum crack width shortly before failure is in the range of 30-80 lm.
Looking at a beam specimen with standardized geometry under 3-point bending load, the distance between two cracks stabilizes at values around 4 mm, as seen in Fig. 9b, thus cracks are closer in comparison to the dogbone-shaped samples tested under uniaxial tension. This is due to the fact that during tensile testing the load is distributed evenly across the cross section, whereas 3-point-bending leads to an area with concentrated peak stress directly under the loaded midpoint. Tensile tests are decidedly more sensitive to surface defects or pores, which lead to local stress maxima and subsequently, the first crack. In flexural tests, the load is a lot more concentrated under the point of load application; while macro pores still act as imperfections, the appearance of the first crack is locally much more constricted. Further increase in stress then forces more cracks to develop along the beam closer to the supports until one of the cracks grows to a fatal size and fiber rupture occurs [35]. The standardized beam geometry shows a lower ultimate flexural strength in comparison to the miniature beam specimens, hinting to either the influence of a sizeeffect or influences between the different manufacturing processes.

Acoustic emission analysis
The effect of different fiber orientations is reflected in the number of events as well as in the spatial and temporal distribution of the events, as can be observed in the normalized cumulative events over time for each tested miniature beam specimen (Fig. 10). The event count includes the time from loading start to ultimate failure of the beam, post-cracking behavior is not part of the analysis. The event number is normalized to the event count when the specimen fails. Miniature bending beams with an in-plane alignment of 0°show almost linear growth of event numbers until 80% percent of the total event count. Miniature beams with an in-plane alignment angle of 10°show the linear trend until 50%, reflecting the multiple phases of the stress-strain diagrams. The linear trend is followed by an exponential like growth, signaling impending failure. Higher deviations in alignment don't show linear trends at all, instead an exponential growth starts from the beginning, indicating sudden brittle failure.
The position of the final crack is measured by means of photogrammetry and shown by a red patch on the time position diagram. For miniature beam specimens with an alignment of 0°, the events are distributed lengthwise around a larger area of the specimen. In contrast, deviating angles show the events concentrated around the final crack position. This is especially visible before maximum load is reached as in Fig. 11a. Dye penetration tests show different fracture patterns for different fiber orientations. A comparison between the distribution of spatial events and crack pattern for a miniature beam with an in-plane-alignment angle of 0°is given in Fig. 11a. The spatial event distribution and the fracture positions match and show a bimodal distribution around the final crack. Event clusters are distributed between 27 and 37 mm, whereas results for miniature beams with an in-plane alignment angle of 40°deviate clearly (see Fig. 11b). The event distribution follows a unimodal distribution around the final fracture, which is reflected by the absence of visible fractures from dye penetration test. When disregarding single events, as shown in the histogram, the majority of all events are localized between 36 und 40.2 mm. The majority of all events is therefore within a range of two millimeters around the fatal crack. This distribution can either be attributed to uncertainty in localization or to the extent of fibers outside the crack and fiber pullout events.
Based on the temporal distribution, it can be stated that first event clusters are not necessarily attributed to the fatal fracture. Instead, as seen in the experiments and also described elsewhere in the literature, fractures  can develop later in the loading process [62]. This suggests that several flaws within the material cumulate to a major crack while other initial cracks and crack branches die out. This effect can be seen clearly in Fig. 11b, where the crack at position 1 produced the first localized events (2) in the experiment. Displacement signals captured by the LDV were decomposed into a frequency-time representation using a continuous wavelet transformation and a Gabor wavelet. Thus, the temporal occurrence of different frequencies becomes visible. Figure 12 shows characteristic events representing different frequency contents. Figure 12 provides characteristic signals, with increasing frequency content. The lowest frequencies shown in Fig. 12a can be linked to extended matrix cracking whereas the higher frequencies can be attributed to different mechanisms. Fiber fracture is generally attributed to high frequency content, and thus the signal in Fig. 12d is considered to be a fiber fracture. Figure 12b, c show intermediate frequencies which can be assigned to fiber pull out events. All signals exhibit large amplitudes at 85 kHz which probably stem from specimen eigenfrequencies. In order to distinguish between different fracture Fig. 11 a Spatial and temporal event distribution for a miniature beam with 0°fiber orientation. The region marked in red indicates the dimensions of the fatal fracture. b Spatial and temporal event distribution for miniature beams with 40°fiber orientation. Position 2 marks first localized events which correspond to the crack marked at position 1. Marker 3 is pointing to the fatal crack modes and eigenfrequencies, modal analysis measurements are required and will be considered in future studies.
During the experiment, the specimen will change its exact position and orientation relative to the fixed laser beam. Thus, the movement of the specimen leads to a change in the reflectivity of the laser and therefore a change in the Green's function. As a consequence, some signals are superimposed by heavy noise and cannot be analyzed. Therefore, a statistical approach for data analysis is not feasible and therefore this work is limited to a phenomenological representation of the results. As outcome of the frequency analysis the following three observations can be reported: • Most of the recorded events seem to include more than one fracture mechanism in the time frame of the transient signal and therefore, combined processes of matrix-and fiber fracture must be considered. Fig. 12 Selection of characteristic AE signals, which were recorded on a miniature beam with an alignment of 40°. The signals are arranged in ascending frequency content from (a-d).
a Example of a low frequency signal with maximum amplitudes at 25 kHz and is considered a matrix fracture event. b Significant amplitudes at 85 and 350 kHz, whereas 85 kHz is regarded to be a resonance within the specimen. c Signal with dominant frequency around 650 kHz can be attributed to fiber pullout. The highest frequency components recorded within the experiments are probably associated with fiber breakage are shown in (d) • Fiber breakage and pullout can be detected and can possibly be distinguished on frequency basis.
3.5 Damage patterns and crack spacing segmented using X-ray CT X-ray CT scans were carried out on miniature bending beams with a fiber content of 1 vol%. The results of the acoustic emission analysis indicate that the main fracture processes occur symmetrically around the loaded area at the midpoint of the miniature bending beams. Thus, an area of about 12 mm along the fracture zone is scanned during X-ray CT measurements and the damage is segmented. Usually, the fiber-reinforced miniature beams will not be completely split in two after failure and stay connected over a small area in the compressive zone. The process zone of samples that stay in one piece is scanned directly in the tomograph. For all following figures, the fatal crack is always visualized in red, while non-fatal microcracks are shown in teal. In the case of the miniature beam with 0°fiber alignment angle, the miniature beam split completely during testing, thus each half of the beam was scanned separately and the main crack is shown as an unsegmented large gap with only the microcracks being segmented and displayed in teal color.
Analysis of microstructural damage confirms the effects also seen in the AE tests. The strain-hardening response diminishes as seen in decreasing strain-atbreak values the more fiber alignment deviates from the direction of the occurring tensile stress. Also, transmission of stress onto the fibers becomes gradually worse and thus cracks cannot be effectively bridged anymore, leading to the sharp decrease in flexural strength and strain-at-break values seen in Fig. 7. As a consequence, fewer microcracks can be observed in the specimen. The sample with 0°a lignment angle (Fig. 13a) shows multiple, very closely spaced microcracks in its process zone. At a 10°deviation in alignment (Fig. 13b), the number of observable cracks is already more than halved. At 20°a ngular deviation (Fig. 13c), only 3 cracks remain visible. For the samples with alignment angles of 30°a nd 40° (Fig. 13d, e), only one crack could be observed within the process zone, showing that the fiber reinforcement is no longer capable of bearing the tensile load and explaining the brittle behavior of the specimens. A brief quantitative approach to measuring crack distances from the CT datasets can be found in the Supplementary Material to this paper.

Conclusion
Additions of chopped carbon fibers leads to a strainhardening cementitious material, the mechanical performance of which is intrinsically linked to its microcracking behavior. Alignment of carbon fibers during the extrusion process, by which test specimens are prepared, leads to an altered fracture process. Microcracks bridged by fibers demand higher energy to propagate further, rendering formation of new cracks energetically more favorable [63]. Thus, the material no longer shows sudden single-crack brittle failure, but rather pseudoductility and multiple cracking. Transmission of tensile load from the cementitious matrix to the fiber leads to significantly increased tensile and flexural strength.
Since acoustic emission analysis (AEA) is directly related to the cracking behaviour, appropriate measuring techniques have been applied based on former experimental investigations studying fiber reinforced concrete [34]. As one of only a few non-destructive testing techniques, AEA is able to investigate the spatiotemporal development of deteriorations like fiber breaking and matrix cracks. Signal-based techniques have been applied including a 1D location based on linear least squares and an automatic picking algorithm (AIC). A modified technique for the discrimination between noise and AE signals was used. To perform higher-order signal investigations, a LDV was additionally used for evaluations of the signal frequency content including wavelet transforms [64]. The application of the LDV-technique enabled for a recording of the true ground motion (and not normalized amplitudes only) originating from the emitted AE waves. In further studies, with adapted specimen geometries, this will be used to derive the true seismic moment [65] of the cracks and therefore the real size of the crack and the released energy can be calculated. This will in future be performed along with applications of ultrasound techniques including coda wave interferometry (CWI) [66]. CWI was already applied but the number of data was yet too little to be processed in a reliable way. Another promising approach tested elsewhere at the time being is the use of an etalon-based optical microphone for the detection of acoustic waves in a non-contact way [67].
Extrusion of the fibrous paste through a tight nozzle allows for excellent control of fiber alignment; CT measurement show a number of 70% of well-aligned fibers, which is close to numbers reported elsewhere on this system. Microcracking during loading can be seen in stress-strain curves as a strain-hardening phase, the process is very well observable by means of DIC and acoustic emission. X-ray CT scans taken after failure confirm the processes seen in the aforementioned methods. The strength of the material and the amount of microcracks formed correlates strongly with alignment of fibers in direction of occuring tensile forces. As fiber alignment deviates, both pseudoductility and strength decrease, indicating that transmission of tensile load from matrix to fiber gets worse. The reduced transmission of loads leads to a vastly reduced capability for microcracking in the material and thus a change from ductile to brittle material behavior. Effects of the reinforcement can be observed up until an angular deviation of 20°, at 30°h ardly any effects and at 40°no effects of an effective reinforcement can be seen anymore.
Acknowledgements The researchers extend their gratitude to everyone involved in the priority programme SPP 2020, as the working atmosphere between the groups is shaped by a sense of fruitful collaboration. The Chair of Solid State and Materials Chemistry thanks Schwenk Zement KG for kindly supplying the cement used in this research. Funding Open Access funding enabled and organized by Projekt DEAL. This research was funded by the German Research Foundation DFG as part of the 'Priority Programme SPP 2020: Cyclic deterioration of High-Performance Concrete in an experimental-virtual lab' (Grant Numbers VO 829/13-1, FI 1720/7-1, GR1664/13-1, PE 1464/6-1).

Fig. 13
Cracking patterns in miniature bending beams with a fiber content of 1 vol%. In-plane alignment angles in the samples are a 0°, b 10°, c 20°, d 30°and e 40°as defined in Fig. 1. The figures show a side view on the samples, with about 1 mm of the sample digitally removed to high-light crack progression. The fatal crack is shown in red, non-fatal microcracks are shown in teal color. The 0°sample separated completely during testing, the fatal crack is shown as large gap