Search results

The purpose of this paper is to present the application of the boundary element method (BEM) in linear elastic fracture mechanics for analysis of fatigue crack propagation problems in mixed-mode (I+II) using a robust academic software named BemCracker2D and its graphical interface BemLab2D.

The methodology consists in calculating elastic stress by conventional BEM and to carry out an incremental analysis of the crack extension employing the dual BEM (DBEM). For each increment of the analysis, the stress intensity factors (SIFs) are computed by the J-Integral technique, the crack growth direction is evaluated by the maximum circumferential stress criterion and the crack growth rate is computed by a modified Paris equation, which takes into account an equivalent SIF to obtain the fracture Modes I and II. The numerical results are compared with the experimental and/or BEM values extracted from the open literature, aiming to demonstrate the accuracy and efficiency of the adopted methodology, as well as to validate the robustness of the programs.

The paper addresses the numerical simulation of fatigue crack growth. The main contribution of the paper is the introduction of a software for simulating two-dimensional fatigue crack growth problems in mixed-mode (I+II) via the DBEM. The software BemCracker2D coupled to the BemLab2D graphical user interface (GUI), for pre/post-processing, are very complete, efficient and versatile and its does make relevant contributions in the field of fracture mechanics.

The main contribution of the manuscript is the development of a GUI for pre/post-processing of 2D fracture mechanics problems, as well as the object oriented programming implementation. Finally, the main merit is of educational nature.

Hybrid methods, wherefore numerical and experimental data are used to calculate a critical parameter, have been used for several years with great success in Experimental Mechanics and, in particular, in fracture mechanics. The purpose of this paper is to report on the comparison of the strain field from numerical modelling forecasts against the experimental data obtained with the digital image correlation method under Mode II loading in fatigue testing. The numerical dual boundary element method has been established in the past as a very reliable method near singular regions where stresses tend to grow abruptly. The results obtained from the strain data near the crack tip were used in Williams expansion and agree fairly well with both the numerical results and the analytical solution proposed for pure Mode II testing.

The work presented in this note is experimental. The proposed methodology is of an hybrid experimental/numerical nature in that a numerical stress intensity factor calculation hinges upon a stress field obtained with an image method.

The obtained results are an important step towards the development of a practical tool for crack behaviour prediction in fatigue dominated events.

The results also stress the necessity of improving the experimental techniques to a point where the methods can be applied in real-life solicitations outside of laboratory premises.

Although several research teams around the globe are presently working in this field, the present research topic is original and the proposed methodology has been presented initially by the research team years ago.

A dual integral formulation for a cracked bar under torsion is derived, and a dual boundary element method is implemented. It is shown that as the thickness of the crack becomes thinner, the ill‐posedness for the linear algebraic matrix becomes more serious if the conventional BEM is used. Numerical experiments for solution instability due to ill‐posedness are shown. To deal with this difficulty, the hypersingular equation of the dual boundary integral formulation is employed to obtain an independent constraint equation for the boundary unknowns. For the sake of computational efficiency, the area integral for the torsion rigidity is transformed into two boundary integrals by using Green’s second identity and divergence theorem. Finally, the torsion rigidities for cracks with different lengths and orientations are solved by using the dual BEM, and the results compare well with the analytical solutions and FEM results.

The purpose of this paper is to develop a methodology for in situ stress intensity factor (SIF) determination that can be used for the analysis of cracked structures. The technique is based on digital image correlation (DIC) combined with an overdetermined algorithm.

The linear overdeterministic algorithm for calculating the SIF based on stress values around the crack tip is applied to a strain field obtained by DIC.

As long as the image quality is sufficiently high, a good accuracy can be obtained for the measured SIF. The crack tip can be automatically detected based on the same strain field. The use of the strain field instead of the displacement field, eliminates problems related to the rigid body motion of the analysed structure.

In future works, based on the applied techniques, the SIF of complex cracked plane stress structures can be accurately determined in real engineering applications.

The paper demonstrates application of known techniques, refined for other applications, also the use of stress field for SIF overdeterministic calculations.

This paper aims to provide SIF calculation method for engineering application.

In this paper, the stress intensity factors (SIFs) calculation method is applied to the anisotropic Ni-based single crystal film cooling holes (FCHs) structure.

Based on contour integral, the anisotropic SIFs analysis finite element method (FEM) in Ni-based single crystal is proposed. The applicability and mesh independence of the method is assessed by comparing the calculated SIFs using mode of plate with an edge crack. Anisotropic SIFs can be calculated with excellent accuracy using the finite element contour integral approach. Then, the effect of crystal orientation and FCHs interference on the anisotropic SIFs is clarified. The SIFs of FCH edge crack in the [011] orientated Ni-based single crystal increases faster than the other two orientations. And the SIF of horizontal interference FCHs edge crack is also larger than that of the inclined interference one.

The SIFs of the FCH edge crack in the turbine air-cooled blade are innovatively computed using the sub-model method. Both the Mode I and II SIFs of FCHs edge crack in blade increase with crack growing.

This paper aims to develop a new isogeometric boundary element formulation based on non-uniform rational basis splines (NURBS) curves for solving Reissner’s shear-deformable plates.

The generalized displacements and tractions along the problem boundary are approximated as NURBS curves having the same rational B-spline basis functions used to describe the geometrical boundary of the problem. The source points positions are determined over the problem boundary by the well-known Greville abscissae definition. The singular integrals are accurately evaluated using the singularity subtraction technique.

Numerical examples are solved to demonstrate the validity and the accuracy of the developed formulation.

This formulation is considered to preserve the exact geometry of the problem and to reduce or cancel mesh generation time by using NURBS curves employed in computer aided designs as a tool for isogeometric analysis. The present formulation extends such curves to be implemented as a stress analysis tool.

The unstable dynamic propagation of multistage hydrofracturing fractures leads to uneven development of the fracture network and research on the mechanism controlling this phenomenon indicates that the stress shadow effects around the fractures are the main mechanism causing this behaviour. Further studies and simulations of the stress shadow effects are necessary to understand the controlling mechanism and evaluate the fracturing effect.

In the process of stress-dependent unstable dynamic propagation of fractures, there are both continuous stress fields and discontinuous fractures; therefore, in order to study the stress-dependent unstable dynamic propagation of multistage fracture networks, a series of continuum-discontinuum numerical methods and models are reviewed, including the well-developed extended finite element method, displacement discontinuity method, boundary element method and finite element-discrete element method.

The superposition of the surrounding stress field during fracture propagation causes different degrees of stress shadow effects between fractures and the main controlling factors of stress shadow effects are fracture initiation sequence, perforation cluster spacing and well spacing. The perforation cluster spacing varies with the initiation sequence, resulting in different stress shadow effects between fractures; for example, the smaller the perforation cluster spacing and well spacing are, the stronger the stress shadow effects are and the more seriously the fracture propagation inhibition arises. Moreover, as the spacing of perforation clusters and well spacing increases, the stress shadow effects decrease and the fracture propagation follows an almost straight pattern. In addition, the computed results of the dynamic distribution of stress-dependent unstable dynamic propagation of fractures under different stress fields are summarised.

A state-of-art review of stress shadow effects and continuum-discontinuum methods for stress-dependent unstable dynamic propagation of multiple hydraulic fractures are well summarized and analysed. This paper can provide a reference for those engaged in the research of unstable dynamic propagation of multiple hydraulic structures and have a comprehensive grasp of the research in this field.

Although organizational citizenship behaviour (OCB) is a concept associated with significant values within the modern workplace, many employees find it challenging to exhibit some necessary extra-role behaviours, such as helping co-workers complete their duties when a situation demands it. Currently, in the South African workspace, fostering OCB among employees is a concern to people practitioners. Specifically, extra-role types of behaviour are declining among professionals as 21st-century technologies promote remote-working policy, leaving employees to work robotically with computers and having no colleagues around to seek or render assistance with their duties. Moreover, professionals are overwhelmed with the timely and endless obligations received frequently and hardly have time and energy for extra-role behaviours. In addition, physical and psychological health-related concerns such as technology stress and career worries are among the contemporary issues affecting human resource (HR) management in this present time. This necessitates more scholarly actions in the niche of OCB, especially identifying and arresting its hindrances. However, a thorough review of the literature on OCB revealed a paucity of scientific reports in the areas of relationships between technostress, career concerns and OCB. Hence, the proposed chapter seeks to bridge the gap in the literature of OCB by theoretically exploring the relationships between technostress, career concerns and OCB in the professional services context in South Africa.

Based on the error analysis, the authors proposed a new kind of high accuracy boundary element method (BEM) (HABEM), and for the large-scale problems, the fast algorithm, such as adaptive cross approximation (ACA) with generalized minimal residual (GMRES) is introduced to develop the high performance BEM (HPBEM). It is found that for slender beams, the stress analysis using iterative solver GMRES will difficult to converge. For the analysis of slender beams and thin structures, to enhance the efficiency of GMRES solver becomes a key problem in the development of the HPBEM. The purpose of this paper is study on the preconditioning method to solve this convergence problem, and it is started from the 2D BE analysis of slender beams.

The conventional sparse approximate inverse (SAI) based on adjacent nodes is modified to that based on adjacent nodes along the boundary line. In addition, the authors proposed a dual node variable merging (DNVM) preprocessing for slender thin-plate beams. As benchmark problems, the pure bending of thin-plate beam and the local stress analysis (LSA) of real thin-plate cantilever beam are applied to verify the effect of these two preconditioning method.

For the LSA of real thin-plate cantilever beams, as GMRES (m) without preconditioning applied, it is difficult to converge provided the length to height ratio greater than 50. Even with the preconditioner SAI or DNVM, it is also difficult to obtain the converged results. For the slender real beams, the iteration of GMRES (m) with SAI or DNVM stopped at wrong deformation state, and the computation failed. By changing zero initial solution to the analytical displacement solution of conventional beam theory, GMRES (m) with SAI or DNVM will not be stopped at wrong deformation state, but the stress error is still difficult to converge. However, by GMRES (m) combined with both SAI and DNVM preconditioning, the computation efficiency enhanced significantly.

This paper presents two preconditioners: DNVM and a modified SAI based on adjacent nodes along the boundary line of slender thin-plate beam. In the LSA, by using GMRES (m) combined with both DNVM and SAI, the computation efficiency enhanced significantly. It provides a reference for the further development of the 3D HPBEM in the LSA of real beam, plate and shell structures.