The effect of material combinations and relative crack size to the stress intensity factors at the crack tip of a bi-material bonded strip

Several types of singular stress fields may appear at the corner where an interface between two bonded materials intersects a traction-free edge depending on the material combinations. Since the failure of the multi-layer systems often originates at the free-edge corner, the analysis of the edge interface crack is the most fundamental to simulate crack extension. In this study, the stress intensity factors for an edge interfacial crack in a bi-material bonded strip subjected to longitudinal tensile stress are evaluated for various combinations of materials using the finite element method. Then, the stress intensity factors are calculated systematically with varying the relative crack sizes from shallow to very deep cracks. Finally, the variations of stress intensity factors of a bi-material bonded strip are discussed with varying the relative crack size and material combinations. This investigation may contribute to a better understanding of the initiation and propagation of the interfacial cracks.     

Determination of sensitivity coefficients of the elastic effective properties for periodic fiber-reinforced composites using the response function method

Derivation and implementation of the homogenization method including determination of sensitivity gradients of the effective elasticity tensor using combined numerical-analytical approach are addressed in this paper. This is possible thanks to an application of the numerical response function together with the effective moduli method known from classical homogenization theory. Computational procedure is implemented using 4-noded quadrilateral plane strain finite elements (program MCCEFF) and the symbolic computations system MAPLE. The sensitivity coefficients are determined on the basis of partial derivatives of the homogenized elasticity tensor calculated using the response function method with respect to all composite components’ elastic characteristics. They are further separately subjected to normalization procedure for a final comparison with each other. Such an enriched homogenization procedure is tested on the periodic fiber-reinforced two component composites; the results of computational analysis are compared to the results of the central finite difference approach applied before. Computational methodology proposed here may be further successively applied not only in the context of homogenization method but also to extend various discrete computational techniques like boundary/finite element, finite difference and volumes together with various meshless methods.     

金属玻璃切削加载的形变屈服准则与实验验证

以金属玻璃切削过程中的屈服形变为研究对象,引入源于岩土领域的M-C屈服准则来解决传统Treaca准则及Mises准则不能反映金属玻璃的应力敏感性问题。另外鉴于金属玻璃的温度敏感特性以及切削加工时较高的切削温升,将经过温度修正的改进型M-C屈服准则应用于金属玻璃的切削模型之中。切削力实验表明,基于传统屈服准则、不含温度项的M-C屈服准则以及经过温度修正的M-C屈服准则所建立的切削力模型中,后者的解析解与切削力实测值相比误差最小(平均误差8.92%),说明经过温度修正的M-C屈服准则可以较好地反映金属玻璃切削加载的切削力及材料形变过程,为后续金属玻璃切削机理的深入研究奠定了理论基础。     

A bridging transition technique for the combination of meshfree method with finite element method in 2D solids and structures

For certain continuum problems, it is desirable and beneficial to combine two different methods together in order to exploit their advantages while evading their disadvantages. In this paper, a bridging transition algorithm is developed for the combination of the meshfree method (MM) with the finite element method (FEM). In this coupled method, the MM is used in the sub-domain where the MM is required to obtain high accuracy, and the FEM is employed in other sub-domains where FEM is required to improve the computational efficiency. The MM domain and the FEM domain are connected by a transition (bridging) region. A modified variational formulation and the Lagrange multiplier method are used to ensure the compatibility of displacements and their gradients. To improve the computational efficiency and reduce the meshing cost in the transition region, regularly distributed transition particles, which are independent of either the meshfree nodes or the FE nodes, can be inserted into the transition region. The newly developed coupled method is applied to the stress analysis of 2D solids and structures in order to investigate its’ performance and study parameters. Numerical results show that the present coupled method is convergent, accurate and stable. The coupled method has a promising potential for practical applications, because it can take advantages of both the MM and FEM when overcome their shortcomings.     

Evaluating stress intensity factors due to weld residual stresses by the weight function and finite element methods

This paper presents a study on the application of the weight function and finite element methods to evaluate residual stress intensity factors in welded test samples. Three specimen geometries and various residual stress profiles were studied. Comparisons of the two different methods were made in terms of the accuracy, easiness to use, conditions and limitations. Calculated residual stress intensity factors by the two different methods are in general in good agreement for all the configurations studied. Computational issues involved in executing these methods are discussed. Some practical issues are also addressed, e.g. treatment of incomplete or limited residual stress measurements, influence of transverse residual stresses, and modelling residual stress in short-length specimens. The finite element method is validated by well-established weight functions and thus can be applied to complex geometries following the procedures recommended in this paper.     

An experimental and modeling study of the thermomechanical behavior of an ABS polymer structural component during an impact test

An experimental and modeling study of the thermomechanical behavior of an ABS polymer structural component during an impact test is presented. The structural component was a heel of a woman's shoe made of ABS polymer material reinforced or not by a pin. Kinematics and thermal full field measurement techniques were used to observe the material and structural component during preliminary experimental tensile and impact tests. With the kinematic fields it was possible to identify the stress–strain response, which takes the necking localization into account. Positive volume variations were also observed during these tensile tests, which were associated with the crazing damage mechanisms in this type of polymer. The thermal fields measured during these tests showed high temperature variations (a few K to 25 K) in the zone where strain was localized.     

An experimental and theoretical study on the prediction of forming limit diagrams using new BBC yield criteria and M–K analysis

In the present paper, a comprehensive study on the prediction of forming limit diagrams (FLDs) for an AA3003-O aluminium alloy is developed theoretically and experimentally. For obtaining the experimental FLDs, an out-of-plane formability test was performed based on the technique proposed by Ozturk and Lee [F. Ozturk, D. Lee, J. Mater. Process. Technol. 170 (2005) 247–253]. The classical Marciniak–Kuczynski (M–K) model and some new yield criteria are utilized to simulate the necking phenomenon and calculate the limit strains theoretically. The employed yield functions are: the BBC2000, BBC2002, and BBC2003 yield criteria proposed by Banabic et al. [D. Banabic, S.D. Comsa, T. Balan, in: Proceedings of the Cold Metal Forming 2000 Conference, Cluj-Napoca, 2000, p. 217; D. Banabic, T. Kuwabara, T. Balan, D.S. Comsa, D. Julean, Int. J. Mech. Sci. 45 (2003) 797–811; D. Banabic, H. Aretz, D.S. Comsa, L. Paraianu, Int. J. Plast. 21 (2005) 493–512]. To calibrate and determine each particular coefficients of performed yield functions an appropriate error-function is defined and minimized by a Newton algorithm. To compare the calculated yield stresses and r-values with experimental data a relative root mean square deviation method presented by Leacock [Alan G. Leacock, J. Mech. Phys. Solids 54 (2006) 425–444] is used. Work-hardening effects on the FLD are analyzed by using Swift and Voce hardening laws. The effect of yield surface on the prediction of numerical FLDs and the number of experimental anisotropy parameters on the accuracy of yield functions are also studied.     

An experimental and numerical study of the effects of length scale and strain state on the necking and fracture behaviours in sheet metals

Within sheet metal forming, crashworthiness analysis in the automotive industry and ship research on collision and grounding, modelling of the material failure/fracture, including the behaviour at large plastic deformations, is critical for accurate failure predictions. In order to validate existing failure models used in finite element (FE) simulations in terms of dependence on length scale and strain state, tests recorded with the optical strain measuring system ARAMIS have been conducted. With this system, the stress–strain behaviour of uniaxial tensile tests was examined locally, and from this information true stress–strain relations were calculated on different length scales across the necking region. Forming limit tests were conducted to study the multiaxial failure behaviour of the material in terms of necking and fracture. The failure criteria that were verified against the tests were chosen among those available in the FE software Abaqus and the Bressan–Williams–Hill (BWH) criterion proposed by Alsos et al, 2008. The experimental and numerical results from the tensile tests confirmed that Barba's relation is valid for handling stress–strain dependence on the length scale used for strain evaluation after necking. Also, the evolution of damage in the FE simulations was related to the processes ultimately leading to initiation and propagation of a macroscopic crack in the final phase of the tensile tests. Furthermore, numerical simulations using the BWH criterion for prediction of instability at the necking point showed good agreement with the forming limit test results. The effect of pre-straining in the forming limit tests and the FE simulations of them is discussed.     

Numerical and experimental determination of in-plane elastic properties of 2/2 twill weave fabric composites

This paper presents the in-plane elastic properties of 2/2 twill weave, T300 carbon/epoxy, woven fabric composite plates, obtained by both finite element analysis and experiments. A micromechanical, three-dimensional (3D) finite element model of the single layer unit cell of a 2/2 twill weave fabric composite is built, and a homogenization process is implemented. A unit cell is chosen such that it encloses the characteristic periodic repeat pattern in the fabric weave. Detailed geometry together with construction procedures for this new model are developed by using ANSYS Parametric Design Language (APDL). In this respect, the scope for altering the weave and yarn parameters is facilitated. Standard tensile and rail shear tests with modifications are performed for this kind of woven fabric composite. Elastic mechanical properties determined by experiments are presented, and the finite element model is verified. Satisfactory correlation between the predicted and experimental results are obtained.     

Analytical and finite element method design of quartz tuning fork resonators and experimental test of samples manufactured using photolithography 2: comprehensive analysis of resonance frequencies using Sezawa's approximations

Resonance frequency of a quartz tuning fork crystal for use in chips of code division multiple access, personal communication system, and a global system for mobile communication was comprehensively analyzed by an analytical method, Sezawa's approximations and the finite element method. A comparison was also made in a more detailed and comprehensive manner among resonance frequencies calculated by the Sezawa's approximations. From the finite element method analysis results, actual tuning fork crystals were fabricated using mass-production capable positive (subtractive) photolithography, selective etching and subsequent positive (subtractive) photoresist spray coating method. A target resonance frequency of was aimed at and a general scheme of commercially available 32.768 kHz tuning fork resonators was also followed in designing tuning fork geometry, tine electrode pattern and thickness. Comprehensive comparison was made among the modeled and experimentally measured resonance frequencies and the discrepancy explained and discussed. Finite element method analysis results quite closely agreed with the experimentally measured resonance frequencies (32.676-32.933 kHz) of the fabricated tuning fork samples measured at a vacuum level of 10⁻⁵ Torr. The difference between modeling and experimentally measured resonance frequency is attributed to the error in exactly manufacturing tuning fork tine width by photolithography. However, the tuning fork design using finite element method modeling must be modified comprehensively to optimize various design parameters affecting both the resonance frequency and other crystal parameters, most importantly crystal impedance (series resistance).     

The effect of microstructure and processing variables on the yield to ultimate tensile strength ratio in a Nb–Ti and a Nb–Ti–Mo line pipe steel

In a hot rolled Nb–Ti and a Nb–Ti–0.09%Mo micro-alloyed steel, the ratio of yield strength to tensile strength (YS/UTS) was found to be a function of the microstructure and cooling rate in those tests where no coiling simulation and no prior deformation. The coarse bainite or acicular ferrite, which was formed at high cooling rates, raised the YS/UTS ratio under these process conditions. With coiling simulation, the ratio was not sensitive to the cooling rate or the microstructure as coiling allows the recovery of dislocations, thereby decreasing the difference in dislocation density that had arisen between a low and a high cooling rate. Deformation with a 33% reduction below the nil-recrystallisation temperature (T_nr) prior to the transformation, led to a high YS/UTS ratio that ranged from 0.81 to 0.86. The prior deformation, therefore, had a stronger effect on the YS/UTS ratio than microstructural changes through cooling rate variations.     

Experimental and analytical studies on the prediction of forming limit diagrams

Metal forming processes are widely used in industrial productions, automobile bodies, food industries, oil refineries, and liquid and gas transmission systems. Analyzing these processes is very important to reduce wastes and optimize the processes. Study of some main factors such as physical and mechanical properties of material and its formability, die geometry, die material, lubrication and pressing speed has been the topic of many research projects. In this paper, forming limit diagrams (FLDs) for LC and ULC steels and the effect of different parameters like the work-hardening exponent, n, and the plastic strain ratio, r, on these diagrams have been evaluated and simulated using ABAQUS/Standard. In this case, Hill’s quadratic anisotropy function is assumed to be the yield function and the Atkins criterion is used as the failure criterion.     

Calculation of J-integrals using experimental and numerical data: Influences of ratio (a/W) and the 3D structure

Contrary to J-integral values calculated from the 2D numerical model, calculated J-integrals [1] from 3D specimen in the numerical and experimental cases are not very close with J-integral used in the literature and two distinct points are present. The first one is according to (a/W) and can be reduced, when this ratio is inferior to 0.2. The second is a structure problem and can be explain by local three-dimensional effects surrounding the crack tip. Two applications using polymer materials for large and minor deformations are experimented. A grid method is used to experimentally determine the in-plane displacement fields around a crack tip in a Single-Edge-Notch (SEN) tensile polyurethane and PMMA specimens. This indirect method composed of experimental in-plane displacement fields and of two theoretical formulations, allows the experimental J-integral to be determined and the results obtained by the numerical simulations to be confirmed.     

Numerical simulation of plasticity-induced fatigue crack closure with emphasis on the crack growth scheme: 2D and 3D analyses

In this paper a numerical simulation of plasticity-induced fatigue crack closure is performed using the finite element method. Emphasis is placed on the crack growth scheme usually adopted for modelling fatigue crack growth in crack closure problems. The number of load cycles between node releases usually reported in the literature has been, in general, one or two. The present work shows that increasing the number of load cycles between node releases has a strong effect on the opening stresses, particularly, under plane strain conditions and 3D fatigue cracks, in contrast plane stress shows little variation with increasing number of load cycles. This investigation also suggests that ratchetting may take place close to the crack tip in both plane strain and 3D crack problems. The problem of discontinuous crack closure under plane strain conditions, often reported in the literature, is also addressed.     

A note on the divergence problem arising in calculation of Green’s function for a 2D piezoelectric thin film strip containing straight line defects and free boundaries: Fourier approach for bodies of unrestricted anisotropy

In this paper we derive a set of novel formulas for computation of the Green’s function and the coupled electro-elastic fields in a 2D piezoelectric strip with free boundaries and containing a distribution of straight line defects. The strip is assumed to be of unrestricted anisotropy, but allowing piezoelectricity, and in this sense situation is more general than in the available literature where only cubic symmetry was investigated. We employ a set of already known analytic formulas for the Fourier amplitude of the Green’s function and the corresponding electro-elastic fields. The key novelty of this paper is solution for the divergence problem occurring during integration of the Fourier amplitude. This problem is caused by poles at k = 0 in various matrix components of the amplitude. From purely mathematical point of view such poles lead to quantities which do not tend to zero at infinity, and this situation is clearly unphysical. To resolve this issue it is demonstrated by means of rigorous analysis that when some additional physical conditions are imposed, physical fields exhibit regular behavior at infinity - the poles do not contribute. Nevertheless, they lead to irremovable numerical ∞ − ∞ uncertainties spreading over the whole domain of integration. This motivates us to compute exact formulas for all these poles to enable engineering calculations involving the system in question.