Search

The FLUKA Monte Carlo transport code is a well-known simulation tool in High Energy Physics. FLUKA is a dynamic tool in the sense that it is being continually updated and improved by the authors. Here we review the progresses achieved in the last year on the physics models. From the point of view of hadronic physics, most of the effort is still in the field of nucleus--nucleus interactions. The currently available version of FLUKA already includes the internal capability to simulate inelastic nuclear interactions beginning with lab kinetic energies of 100 MeV/A up the the highest accessible energies by means of the DPMJET-II.5 event generator to handle the interactions for greater than 5 GeV/A and rQMD for energies below that. The new developments concern, at high energy, the embedding of the DPMJET-III generator, which represent a major change with respect to the DPMJET-II structure. This will also allow to achieve a better consistency between the nucleus-nucleus section with the original FLUKA model for hadron-nucleus collisions. Work is also in progress to implement a third event generator model based on the Master Boltzmann Equation approach, in order to extend the energy capability from 100 MeV/A down to the threshold for these reactions. In addition to these extended physics capabilities, structural changes to the programs input and scoring capabilities are continually being upgraded. In particular we want to mention the upgrades in the geometry packages, now capable of reaching higher levels of abstraction. Work is also proceeding to provide direct import into ROOT of the FLUKA output files for analysis and to deploy a user-friendly GUI input interface.

A model which can be used to predict the two-dimensional nonlinear behavior of bridged cracks in orthotropic strips is presented. The results obtained using a singular integral equation formulation which incorporates the anisotropy rigorously show that, although the effects of anisotropy are significant, the nondimensional quantities employed by Cox and Marshall can generate nearly universal results (R-curves, for example) for different levels of relative anisotropy. The role of composite constituent properties in the behavior of bridged cracks is clarified.

This paper presents the preliminary results of an on-going study of the near-tip mechanics of mode-I cracking in brittle matrix composite laminates. A finite element model is developed within the context of two competing characteristic lengths present in the composite: the microstructural length (the thickness of the layers) and a macro-length (crack-length, uncracked ligament size, etc.). For various values of the parameters which describe the ratio of these lengths and the constituent properties, the stresses ahead of a crack perpendicular to the laminates are compared with those predicted by assuming the composite is homogeneous orthotropic. The results can be used to determine the conditions for which homogenization can provide a sufficiently accurate description of the stresses in the vicinity of the crack-tip.

Green's functions are derived for the plane elastostatics problem of a dislocation in a bimaterial strip. Using these fundamental solutions as kernels, various problems involving cracks in a bimaterial strip are analyzed using singular integral equations. For each problem considered, stress intensity factors are calculated for several combinations of the parameters which describe loading, geometry and material mismatch.

A solution is presented for the problem of a finite length crack branching off the interface between two bonded dissimilar isotropic materials. Results are presented in terms of the ratio of the energy release rate of a branched interface crack to the energy release rate of a straight interface crack with the same total length. It is found that this ratio reaches a maximum when the interface crack branches into the softer material. Longer branches tend to have smaller maximum energy release rate ratio angles indicating that all else being equal, a branch crack will tend to turn back parallel to the interface as it grows.

Green's functions are derived for the plane elastostatics problem of a dislocation in a bimaterial strip. Using these fundamental solutions as kernels, various problems involving cracks in a bimaterial strip are analyzed using singular integral equations. For each problem considered, stress intensity factors are calculated for several combinations of the parameters which describe loading, geometry and material mismatch.

The results are presented of a stress intensity factor analysis of semicircular surface cracks in the inner raceway of an engine bearing. The loading consists of a moving spherical Hertzian contact load and an axial stress due to rotation and shrink fit. A 3-D linear elastic Boundary Element Method code was developed to perform the stress analysis. The element library includes linear and quadratic isoparametric surface elements. Singular quarter point elements were employed to capture the square root displacement variation and the inverse square root stress singularity along the crack front. The program also possesses the capability to separate the whole domain into two subregions. This procedure enables one to solve nonsymmetric fracture mechanics problems without having to separate the crack surfaces a priori. A wide range of configuration parameters was investigated. The ratio of crack depth to bearing thickness was varied from one-sixtieth to one-fifth for several different locations of the Hertzian load. The stress intensity factors for several crack inclinations were also investigated. The results demonstrate the efficiency and accuracy of the Boundary Element Method. Moreover, the results can provide the basis for crack growth calculations and fatigue life prediction.

A solution is presented for the problem of a finite length crack branching off the interface between two bonded dissimilar isotropic materials. Results are presented in terms of the ratio of the energy release rate of a branched interface crack to the energy release rate of a straight interface crack with the same total length. It is found that this ratio reaches a maximum when the interface crack branches into the softer material. Longer branches tend to have smaller maximum energy release rate ratio angles indicating that all else being equal, a branch crack will tend to turn back parallel to the interface as it grows.

An algorithm is developed to construct the compliance matrix for a cracked solid in the integral-equation formulation of two-dimensional linear-elastic fracture mechanics. The integral equation is reduced to a system of algebraic equations for unknown values of the dislocation-density function at discrete points on the interval from -1 to 1, using the numerical procedure described by Gerasoulis (1982). Sample numerical results are presented, and it is suggested that the algorithm is especially useful in cases where iterative solutions are required; e.g., models of fiber-reinforced concrete, rocks, or ceramics where microcracking, fiber bridging, and other nonlinear effects are treated as nonlinear springs along the crack surfaces (Ballarini et al., 1984).

An algorithm is presented which can be used to develop compliance matrices for cracked bodies. The method relies on the numerical solution of singular integral equations with Cauchy-type kernels and provides an efficient and accurate procedure for relating applied loadings to crack opening displacements. The algorithm should be of interest to those performing repetitive calculations in the analysis of experimental results obtained from fracture specimens.