Microsoft Word - numero_26_art_10 R. Citarella et alii, Frattura ed Integrità Strutturale, 26 (2013) 92-103; DOI: 10.3221/IGF-ESIS.26.10 92 Coupled FEM-DBEM method to assess crack growth in magnet system of Wendelstein 7-X R. Citarella, M. Lepore Department of Mechanical Engineering, Politecnico di Milano, Milano (Italy) stefano.foletti@polimi.it, antonietta.loconte@polimi.it, stefano.salgarollo@mail.polimi.it, federico.bassi@mail.polimi.it Joris Fellinger, Victor Bykov, Felix Schauer Max Planck Institute for Plasma Physics, EURATOM Association, Wendelsteinstr. 1, 17491 Greifswald, Germany ABSTRACT. The fivefold symmetric modular stellarator Wendelstein 7-X (W7-X) is currently under construction in Greifswald, Germany. The superconducting coils of the magnet system are bolted onto a central support ring and interconnected with five so-called lateral support elements (LSEs) per half module. After welding of the LSE hollow boxes to the coil cases, cracks were found in the vicinity of the welds that could potentially limit the allowed number N of electromagnetic (EM) load cycles of the machine. In response to the appearance of first cracks during assembly, the Stress Intensity Factors (SIFs) were calculated and corresponding crack growth rates of theoretical semi-circular cracks of measured sizes in potentially critical position and orientation were predicted using Paris’ law, whose parameters were calibrated in fatigue tests at cryogenic temperature. In this paper the Dual Boundary Element Method (DBEM) is applied in a coupled FEM-DBEM approach to analyze the propagation of multiple cracks with different shapes. For this purpose, the crack path is assessed with the Minimum Strain Energy density criterion and SIFs are calculated by the J- integral approach. The Finite Element Method (FEM) is adopted to model, using the commercial codes Ansys or Abaqus;, the overall component whereas the submodel analysis, in the volume surrounding the cracked area, is performed by FEM (“FEM-FEM approach”) or alternatively by DBEM (“FEM-DBEM approach”). The “FEM-FEM approach” considers a FEM submodel, that is extracted from the FEM global model; the latter provide the boundary conditions for the submodel. Such approach is affected by some restrictions in the crack propagation phase, whereas, with the “FEM-DBEM approach”, the crack propagation simulation is straightforward. In this case the submodel is created in a DBEM environment with boundary conditions provided by the global FEM analysis; then the crack is introduced and a crack propagation analysis has been performed to evaluate the effects of the crack shape and of the presence of nearby cracks on the allowed number of EM load cycles. SOMMARIO. I nuclei superconduttori del sistema magnetico di Wendelstein 7-X (W7-X) sono collegati tramite bulloni su un anello di supporto centrale e connessi tra loro tramite cinque cosiddetti elementi di supporto laterale (LSE) per ogni metà modulo. A valle della saldatura tra gli elementi scatolari LSE e i nuclei superconduttori è stata evidenziata la presenza di cricche in prossimità delle saldature che potenzialmente potevano limitare il numero consentito N di cicli elettromagnetici da applicare alla macchina. In conseguenza dell’apparizione delle prime cricche durante la fase di assemblaggio, si è proceduto al calcolo degli Stress Intensity Factors (SIFs) e, tramite la formula di Paris (calibrata con test di fatica a temperatura criogenica), al calcolo delle corrispondenti velocità di accrescimento, per cricche di forma iniziale semicircolare e dimensioni come da evidenza sperimentale, collocate in posizione ed orientazione potenzialmente critiche. In un precedente http://dx.medra.org/10.3221/IGF-ESIS.26.10&auth=true http://www.gruppofrattura.it R. Citarella et alii, Frattura ed Integrità Strutturale, 26 (2013) 92-103; DOI: 10.3221/IGF-ESIS.26.10 93 lavoro si è ipotizzata una forma semicircolare per una singola cricca con limiti di accrescimento tali da non incorrere in situazioni di potenziale coalescenza tra cricche diverse, anche se tali ipotesi restrittive determinavano una notevole diminuzione del numero di cicli ammissibile. In questo lavoro si cerca di rilassare le limitazioni imposte dalle precedenti assunzioni, con la considerazione di cricche multiple interagenti e di forma iniziale anche semiellittica. A tale scopo si fa riferimento ad una procedura accoppiata FEM-DBEM (Dual Boundary Element Method) nella quale il percorso di propagazione è definito sulla base del criterio della Minima Densità di Energia Deformazionale mentre i SIF sono calcolati a partire dal J-integral.Si è proceduto alla modellazione del componente globale in ambiente FEM (Ansys o Abaqus) mentre l’analisi di un sottomodello (estratto dal precedente) e costituito dal volume circostante la cricca, è stata realizzata sia in ambiente FEM che DBEM, in modo da disporre di un confronto tra i SIF calcolati con due procedure numeriche differenti. Con l’approccio ibrido FEM-DBEM è stato possibile valutare gli effetti della geometria della cricca iniziale e della presenza di cricche adiacenti sul numero di cicli di carico ammissibili sulla macchina. KEYWORDS. FEM; DBEM; Superconductive coils; Crack propagation; Electro-magnetic cycles. INTRODUCTION he fivefold symmetric modular stellarator Wendelstein 7-X (W7-X) is currently under construction in Greifswald, Germany [1]. Each of the five modules are build up of two flip symmetric half modules (ten in total) which consist of 5 non-planar (NP) and 2 planar superconducting coils each to be operated at 4K (Fig. 1). The NP coil cases are made of cast stainless steel (SS) EN 1.3960. They are bolted onto a SS central support ring and welded together at the outboard side of the torus via the LSE’s consisting of 100 – 150 mm long hollow tubes of 30 – 35 mm thick forged SS EN 1.4429 (Fig. 1). After welding the LSEs (weld depth 15-30 mm), many surface cracks substantially larger than 8 mm (typical acceptance limit of EN 23277 [2]) were found with dye penetration tests at the accessible surfaces, particularly at the coil side of the weld within the cast steel, oriented parallel to the weld seam. Figure 1: Fivefold symmetric FE model of the magnet system of Wendelstein 7-X (left) and detail of a half module where the LSE's are shown (right). In a first paper [3], the averaged SIF along the crack front of these cracks were predicted using an analytical method and finite element method (FEM) models which were benchmarked against each other and against a single run with a Dual Boundary Element Method (DBEM) model. Assuming proportionality between the SIF and square root of the crack size ( K a  ) and an arbitrary maximum crack size, it was possible to derive the number of load cycles until failure from the SIF of the initial crack. The maximum radial crack size increase was limited to 10 mm to avoid coalescence of the T http://dx.medra.org/10.3221/IGF-ESIS.26.10&auth=true http://www.gruppofrattura.it R. Citarella et alii, Frattura ed Integrità Strutturale, 26 (2013) 92-103; DOI: 10.3221/IGF-ESIS.26.10 94 observed crack with adjacent cracks and to ensure that the crack remains small compared to the cross section, i.e. to ensure the proportionality rule K a  . Based on Paris' law parameters derived from 10 fatigue crack growth rate tests (FCGR) carried out at cryogenic temperature at Karlsruhe Institute of Technology (KIT), the predicted SIF could be related to cyclic crack growth [3]. In the current paper, the assessment is extended, for the most critical crack, using the DBEM [4-11]. This most critical crack was located in LSE-D05 in the cast steel (Fig. 1). In the DBEM method, the stress state and SIFs along the crack front are updated at each step of crack advance; as a consequence, deviation from the proportionality rule is automatically simulated. Moreover, an adjacent crack was included in the model which allows for the simulation of crack coalescence. So, no maximum crack size needs to be defined and crack growth can be continued until the critical SIF is reached, i.e. the SIF where unstable crack growth starts. In addition, the method calculates the crack growth rate along the crack front rather than assuming a pure radial growth as done in [3] with a pure FEM approach. So, it enables the prediction of crack growth of non-circular cracks, differently from what imposed by the FEM sub-modelling approach, where only semicircular crack fronts are modelled [3]. The depth of the crack is fundamentally unknown but from repair experience it was found that the crack depth was typically smaller than half the crack length. To verify the influence of the assumed crack shape, a semi-elliptical crack has been compared with the original semi-circular crack. In the next section, the DBEM and FEM models are presented and the results of the FEM and DBEM submodels are compared. In the DBEM submodel also the effect of an elliptical crack, adjacent to the main crack is evaluated. Finally the conclusions are given. MODELS AND RESULTS ince the mechanical behaviour of the magnet system is of crucial importance for the operation of W7-X, two independent global FEM models of the magnet system have been developed, respectively using the commercial codes Ansys and Abaqus, which are successfully benchmarked against each other. The displacements and generalised sectional forces and moments typically deviate less than 10% between both global models when loaded with a magnetic field on the plasma axis of 3 T (Tesla). An FEM submodel analysis around the most critical crack has been made in Abaqus: a crack is introduced in such submodel (loaded on the cutting boundaries with displacements from the Abaqus global model) and the SIF’s are calculated along the front. Two different DBEM submodels are considered: one extracted from Ansys (loaded on the cutting boundaries with displacements from the Ansys global model) and another one extracted from the Abaqus submodel, with a further close- up around the cracked area (again the boundary conditions on the cutting face are taken for the Abaqus submodel). FEM model and submodel Preparing the FEM submodel, it was found that it suffices to include in the submodel only the LSE without the adjacent coils. So only LSE D05 was modelled with the semi-circular crack of 14 mm diameter at the outside of the cast steel (Fig. 2). The crack is modelled as a seam in the mesh with elements following the contour of the crack. At the same time, it was found that the results are sensitive to the weld modelling. Notably, the weld does not penetrate through the entire thickness of the forged steel tube, so part of the cross section remains un-welded (see right picture of Fig. 2). If the un- welded area is at the outside of the cross section, it shields the crack from stresses perpendicular to the crack. Vice versa, the perpendicular stresses around the crack are increased when the un-welded area is inside of the cross section. For the critical crack under consideration, the stresses perpendicular to the crack are increased. Since the mesh is adapted to the crack shape and size, the model does not allow an automatic crack growth. In Fig. 3 the Von Mises stresses for the cast steel around the crack, under an EM load equal to 3T, are presented for both the Ansys and Abaqus models (and submodel): it is to be pointed out that the LSE geometry is less accurate in the Ansys global model than in the Abaqus global model so that a rigorous comparison between the stress scenarios is prevented. DBEM submodel The DBEM sub-modelling is aimed at analyzing the propagation of cracks with elliptical or semicircular shapes and, eventually, affected by nearby cracks. For this purpose, the Dual Boundary Element Method (DBEM) is applied in a coupled FEM-DBEM approach [12-15] with the crack path assessed by the Minimum Strain Energy density criterion and the Stress Intensity Factors (SIFs) calculated by the J-integral approach [4-6]. The Finite Element Method (FEM) is S http://dx.medra.org/10.3221/IGF-ESIS.26.10&auth=true http://www.gruppofrattura.it R. Citarella et alii, Frattura ed Integrità Strutturale, 26 (2013) 92-103; DOI: 10.3221/IGF-ESIS.26.10 95 adopted to model the overall component, whereas the submodel analysis in the volume surrounding the cracked area is performed by DBEM. The boundary conditions, applied on the DBEM local model, are the displacements calculated from the ANSYS global model (Fig. 4), with related material properties listed in Tab. 1, or from the Abaqus submodel (Fig. 5) with related material properties listed in Tab. 2. Figure 2: Abaqus submodel with highlight of the critical crack considered in LSE D05 (left) and cut through the weld to show the penetration depth of the weld and the cold contact zones (right). Figures 3: Comparison of Von Mises equivalent stress in the LSE as extracted from: Ansys global model with no crack introduction (left); Abaqus global model, with no crack introduction (centre) and with crack introduction (right). Figure 4: DBEM submodel extracted from Ansys global model, with highlight of imported displacement boundary conditions on the cutting surfaces and related subdivision in zones. In Fig. 4 the DBEM submodel is shown, as obtained by a Boolean operation of subtraction of a spherical domain from the overall model; the sphere center coordinates are x=5.56659 m, y=2.95385 m, z=0.476757 m and the radius is equal to 0.15 m (with such value the submodel is judged sufficiently bigger than the cracked area); in the same figure the submodel splitting in 5 zones is highlighted (zone 4 and 5 have orthotropic properties). http://dx.medra.org/10.3221/IGF-ESIS.26.10&auth=true http://www.gruppofrattura.it R. Citarella et alii, Frattura ed Integrità Strutturale, 26 (2013) 92-103; DOI: 10.3221/IGF-ESIS.26.10 96 Zone 1 EN 1.3960 Zone 2 EN 1.4429 Zone 3 Sand embedding Zone 4 Insulation Zone 5 Winding pack Young's modulus [GPa] 158 197 23.5 Ex=20.8 Ey= Ez =31.8 Gxy=Gxz=6.7 Gyz=7.20 Ex= Ey=26.5 Ez =41.2 Gxy=5.7 Gxz= Gyz=10.4 Poison's ratio [-] 0.3 0.3 0.32 yz=0.27 xzxy=0.21 yz=xz0.2 xy=0.26 Thermal expansion [µm/m] 10.2 10.38 12.1 x=20.4 y=z=7.61 x=y=14 z=14.7 Table 1: Material properties at 4 K temperature. Figure 5: DBEM submodel extracted from Abaqus submodel, with highlight of imported displacement boundary conditions on the cutting surfaces and related subdivision in zones (the welding affected area is shown in red). Zones 1, 2, 4 - EN 1.4429 Zone 3 (cracked) - EN 1.3960 Young's modulus [GPa] 197 158 Poison's ratio [-] 0.3 0.3 Thermal expansion [µm/m] 10.38 10.2 Table 2: Material properties at 4 K temperature. A DBEM thermal-stress analysis [9] is performed on the submodel extracted from Ansys (Fig. 4) in order to allow for the coil cooling from an ambient temperature equal to 293 K to the superconductive temperature equal to 4 K, and subsequent electromagnetic loading from 0 to 3 T. In Fig. 6 the DBEM submodels, with maximum principal stresses (MPa) are shown, in the configuration before and after the crack insertion (the initial crack is semicircular with a radius equal to 7 mm). A DBEM stress analysis (this time the thermal part is considered with negligible effects, namely the stress coming from cooling at 4 K are neglected) is performed on the submodel extracted from Abaqus (Fig. 5) with an electromagnetic loading equal to 3 T. In Fig. 7 the DBEM submodels, with maximum principal stresses (MPa) are shown for comparison with the results illustrated in Fig. 6: some limited differences pop up as expected due to the huge difference in submodel and global model mesh refinement, for not mentioning the geometry differences; there is also a small contribution to the difference coming from the different loading process (in one case there is no allowance for cooling effects). (4) (3)(1)(2) http://dx.medra.org/10.3221/IGF-ESIS.26.10&auth=true http://www.gruppofrattura.it R. Citarella et alii, Frattura ed Integrità Strutturale, 26 (2013) 92-103; DOI: 10.3221/IGF-ESIS.26.10 97 Figure 6: DBEM submodel before and after initial crack insertion, with highlight of Maximum Principal stresses (MPa). A, B, and C are the points obtained by the intersection between the crack front and the ellipse semiaxis. Figure 7: DBEM submodel (extracted from Abaqus submodel), before and after crack insertion, with highlight of Maximum Principal stresses (MPa). RESULTS FEM-FEM vs. FEM-DBEM approaches n Fig. 8 it is possible to compare the J-integral values calculated along the initial crack front with the 3 different approaches available: the FEM-DBEM approach provides two sets of results depending upon the software (Ansys or Abaqus) used for the global modelling, whereas the remaining set of results is related to the crack assessment done on the Abaqus submodel. As expected there is a non negligible discrepancy when considering the results coming from the Ansys global model but it is easy to find the reasons in the rough mesh and in the lack of a detailed geometric reproduction of the LSE; on the contrary a satisfactory correspondence is obtained between the results obtained by FEM and DBEM submodels when considering the Abaqus global model, apart from a non negligible discrepancy localised at the break through points that are more affected by extrapolation effects and mesh inaccuracy. With reference to the latter, with reference to DBEM analysis, when the J-path is built on a triangular element there is a loss of accuracy in the J calculation, but this can be overcome by further refining the mesh on the crack so as to minimise the impact of the last unavoidable triangular elements at break through point (see in Fig. 7 the close up of the mesh along the crack front). DBEM crack propagation analysis The simulation of crack propagation, in this work, is only performed with reference to the DBEM submodel extracted from the Ansys global model. The crack propagation law adopted is the Paris law (Eq.1) calibrated under cryogenic conditions (T=4K) [3]: mKC dN da  (1) with C=7.95·10-10 (with da/dN measured in mm/cycle and K in MPa√m) and m=3.23; the plane strain fracture toughness is IcK =106 MPa√m. The fatigue cycle considered is based on the coil cooling, from ambient temperature to 4 I http://dx.medra.org/10.3221/IGF-ESIS.26.10&auth=true http://www.gruppofrattura.it R. Citarella et alii, Frattura ed Integrità Strutturale, 26 (2013) 92-103; DOI: 10.3221/IGF-ESIS.26.10 98 K, and subsequent application of the electromagnetic load corresponding to a magnetic field cyclically varying from 0 to 3 T. After 6 steps of crack propagation, the crack front break through the LSE superior surface as shown in Fig. 9, where again the maximum principal stresses (MPa) and the final crack shape are shown. Figure 8: J-integral values calculated along the initial crack front with the 3 different available approaches. Figure 9: DBEM submodel, in the final cracked configuration, with highlight of Maximum Principal stresses (MPa). The mode I SIFs (model II and mode III SIFs are much lower) along the crack front for each step of crack advance are shown in Fig. 10: it is possible to see the sudden increase of KI, in the final step, along that part of the crack front that breaks through the upper LSE surface (see Fig. 9). One of the aims of this work was to check the influence on the fatigue life of a different initial crack front geometry, so it was also considered an alternative case with an initial semi-elliptical crack with the same value of a=7 mm and with a depth reduced from c=7 to c=5 mm. The fatigue cycles for each crack advance for the two aforementioned initial configurations are shown in Fig. 11: it is clear the impact of the initial crack geometry. The simulation is ended when the crack breaks through the upper LSE surface, because from now on the SIFs will rapidly increase approaching the fracture toughness. http://dx.medra.org/10.3221/IGF-ESIS.26.10&auth=true http://www.gruppofrattura.it R. Citarella et alii, Frattura ed Integrità Strutturale, 26 (2013) 92-103; DOI: 10.3221/IGF-ESIS.26.10 99 Figure 10: mode I SIFs (MPa√m) along the crack front for each step of crack advance. Figure 11: Semi-axis length (m) vs. fatigue cycles for an initial semicircular or elliptical crack. The second aim of this work was to analyse the reciprocal influence of adjacent cracks and consequently a multiple crack propagation analysis has been performed. The two initial cracks considered are respectively semicircular, with radius equal to 7 mm (the N.1 in Fig. 12), and semi- elliptical (N. 2 in Fig. 12) with semi-axis a and c (depth size) respectively equal to 7 mm and 5 mm; the distance between the centres of the two aforementioned cracks is equal to 28 mm; the corresponding stress scenario is illustrated in Fig. 12. The stress scenario and the crack configuration after four steps of crack propagation are shown in Fig. 13: it is possible to see that from this stage on the crack coalescence can be considered likely to happen. http://dx.medra.org/10.3221/IGF-ESIS.26.10&auth=true http://www.gruppofrattura.it R. Citarella et alii, Frattura ed Integrità Strutturale, 26 (2013) 92-103; DOI: 10.3221/IGF-ESIS.26.10 100 Figure 12: DBEM submodel, in the initial multiple cracked configuration, with highlight of Maximum Principal stresses (MPa). Figure 13: DBEM submodel, in the final multiple cracked configuration, with highlight of Maximum Principal stresses (MPa). Differently for the previous case, this time the SIFs are calculated by the crack opening displacement method (COD), [4] because there were difficulties in enforcing on the secondary crack a mesh along the crack front made of mainly quadrilateral elements and the J-integral accuracy, as previously said, is affected if the J-path is built on a triangular element. The calculated SIFs are shown for the two cracks at each crack advance in Figs. 14-16. Figure 14: mode I SIFs (MPa√m) along the crack front for each step of crack advance for cracks N. 1 and 2. It is to be pointed out the formula used (Eq. 2) for calculating the equivalent SIF used in the Paris formula is: )2)(( 22 IIIIIIeq KKKK  (2) http://dx.medra.org/10.3221/IGF-ESIS.26.10&auth=true http://www.gruppofrattura.it R. Citarella et alii, Frattura ed Integrità Strutturale, 26 (2013) 92-103; DOI: 10.3221/IGF-ESIS.26.10 101 Figure 15: mode II SIFs (MPa√m) along the crack front for each step of crack advance for cracks N. 1 and 2. Figure 16: mode III SIFs (MPa√m) along the crack front for each step of crack advance for cracks N. 1 and 2. Figure 17: Crack size vs. fatigue cycles. http://dx.medra.org/10.3221/IGF-ESIS.26.10&auth=true http://www.gruppofrattura.it R. Citarella et alii, Frattura ed Integrità Strutturale, 26 (2013) 92-103; DOI: 10.3221/IGF-ESIS.26.10 102 The fatigue cycles for each crack propagation step for the two aforementioned cracks are shown in Fig. 17: comparing the double crack results with the single crack results (Fig. 11) it is possible to observe a clear detrimental effect of the second crack on the main semicircular crack (e.g. in case of multiple cracks, the main crack depth reaches a size of about 2 cm after 4.5·104 cycles rather than after 7·104 cycles as in the single crack case). Again the simulation is ended when the main crack breaks through the upper LSE surface, because from now on the SIFs will rapidly increase towards the fracture toughness. CONCLUSIONS n this paper some partial results of a benchmarking activity between two different approaches for the crack growth assessment in a component of the magnet system of Wendelstein 7-X are presented: both approaches are based on the sub-modelling of the critical cracked domain, as extracted from an FEM overall model, but in one case the submodel is still FEM based whereas in the second case the submodel is DBEM based. The latter approach realises a synergic coupled approach between the two numerical procedures, providing a more flexible tool for analysing complex crack scenarios. The adopted codes are ANSYS and ABAQUS for the FEM modelling and BEASY for the DBEM modelling. The SIFs calculated with the FEM-FEM approach and with the FEM-DBEM approach can be successfully compared on the initial crack front. The effect of the initial crack shape and the presence of an adjacent crack have been evaluated: assuming a non-circular crack of 5 mm depth and 14 mm width, instead of the initial circular crack with 7 mm depth, leads to a reduction of crack growth and increases the number of cycles up to reaching KIc from 7·10⁴ to 8.5·10⁴. On the other side, an adjacent semi- circular crack of equal size of 7 mm radius at initially 28 mm centre to centre spacing accelerates crack growth and reduces the number of cycles to 4.5·10⁴. Therefore, when assessing individual cracks neglecting eventual adjacent cracks, it is necessary to prevent coalescence of the cracks by limiting crack growth to half the crack spacing on top of limiting Keff to KIc. REFERENCES [1] Bykov, V., Schauer, F., Egorov, K., Tereshchenko, A., van Eeten, P., Dübner, A., Sochor, M., Zacharias, D., Dudek, A., Chen, W., Czarkowski, P., Sonnerup, L., Fellinger, J., Hathiramani, D., Ye, M.Y., Dänner, W., W7-X team, Structural analysis of W7-X: Overview, Fusion Engineering and Design, 84 (2009) 215–219. [2] EN ISO 23277, CEN, Brussels 2010 [3] Fellinger, J., Bykov, V., Schauer, F., Assessment of cracks in lateral supports of the magnet system of Wendelstein 7- X, Fus. Eng. Des. (2012), http://dx.doi.org/10.1016/j.fusengdes.2012.11.021 (in press) [4] Citarella, R., Cricrì, G., Cricrì, Comparison of DBEM and FEM Crack Path Predictions in a notched Shaft under Torsion, Engineering Fracture Mechanics, 77 (2010) 1730-1749. [5] Citarella, R., Cricrì, G., Lepore, M., Perrella, M., DBEM and FEM Analysis of an Extrusion Press Fatigue Failure, in: A. Öchsner, L.F.M. da Silva, H. Altenbach (Eds.), Advanced Structured Materials with Complex Behaviour - Advanced Structured Materials, Springer-Verlag, Berlin, 3(2) (2010) 181-191. [6] Citarella, R., Perrella, M., Multiple surface crack propagation: numerical simulations and experimental tests, Fatigue and Fracture of Engineering Material and Structures, 28 (2005) 135-148. [7] Citarella, R., Cricrì, G., Armentani, E., Multiple crack propagation with Dual Boundary Element Method in stiffened and reinforced full scale aeronautic panels, in: Aliabadi M.H., Wen P. (Eds), Key Engineering Materials, Trans Tech Publications, Switzerland, doi:10.4028/www.scientific.net/KEM.560.129560, (2013) 129-155. [8] Citarella, R., Buchholz, F.-G., Comparison of crack growth simulation by DBEM and FEM for SEN-specimens undergoing torsion or bending loading, Engineering Fracture Mechanics, 75 (2008) 489–509. [9] Citarella, R., MSD Crack propagation on a repaired aeronautic panel by DBEM, Advances in Engineering Software, doi:10.1016/j.advengsoft.2011.02.014, 42(10) (2011) 887-901. [10] Citarella, R., Non Linear MSD crack growth by DBEM for a riveted aeronautic reinforcement, Advances in Engineering Software, doi:10.1016/j.advengsoft.2008.04.007, 40 (4), 2009. 253–259. [11] Armentani, E., Citarella, R., Sepe, R., FML Full Scale Aeronautic Panel Under Multiaxial Fatigue: Experimental Test and DBEM Simulation, EFM Special Issue on 'Multiaxial Fracture' doi:10.1016/j.engfracmech.2011.02.020, 78(8) (2011) 1717-1728. I http://dx.medra.org/10.3221/IGF-ESIS.26.10&auth=true http://www.gruppofrattura.it R. Citarella et alii, Frattura ed Integrità Strutturale, 26 (2013) 92-103; DOI: 10.3221/IGF-ESIS.26.10 103 [12] Citarella, R., Cricrì, G., Lepore, M., Perrella, M., Thermo-Mechanical Crack Propagation in Aircraft Engine Vane by Coupled FEM-DBEM Approach, Advances in Engineering Software, doi:10.1016/j.advengsoft.2013.07.006, (2014) 57-69. [13] Citarella, R., Cricrì, G., Lepore, M., Perrella, M., Assessment of Crack Growth from a Cold Worked Hole by Coupled FEM-DBEM Approach, Key Engineering Materials, © (2014) Trans Tech Publications, Switzerland, doi:10.4028/www.scientific.net/KEM.577-578.669, (2014) 577-578, 669-672. [14] Citarella, R., Carlone, P., Lepore, M., Palazzo, G., FEM-DBEM Investigation of the Influence of Process Parameters on Crack Growth in Aluminium Friction Stir Welded Butt Joints, Key Engineering Materials, , doi:10.4028/www.scientific.net/KEM.554-557.2118, 554-557 (2013) 2118-2126. [15] Citarella, R., Cricrì, G., A two-parameter model for crack growth simulation by combined FEM-DBEM approach, Advances in Engineering Software, doi:10.1016/j.advengsoft.2008.05.001., 40 (5) (2009) 363–377. http://dx.medra.org/10.3221/IGF-ESIS.26.10&auth=true http://www.gruppofrattura.it