Simulation and experimental studies of the repaired composite laminate

Firstly the compress experiment is undertaken to investigate the efficiency of repaired panels in this paper, and then modeling of the mechanical behavior of the repaired composite panel under compressive static load is conducted by using of the finite element method. The effect of geometric non‐linearity on the stress‐strain response is considered in the numeric analysis. Fatherly, the user material subroutine (UMAT) is integrated with the ABAQUS package with the geometric non‐linearity effect for studying the damage initiation and its progression in the composite structure, and quadrilateral, linear, thick shell elements (S8R) are adopted. Finally, the predicted strain distribution, damage evolution and strength of the laminate are compared with the test results.     

Scholte波波速和孔隙度关系的建立及实验验证

随着声学技术的快速发展,界面波在无损检测及环境参数反演领域的应用越来越广泛,但是目前理论与实践之间还存在一些差距,首先对基于Biot理论的流体饱和多孔介质声传播特性的理论进行了仿真分析,然后运用势函数结合多孔介质固体-流体边界条件建立了两半无限空间流-固界面Scholte波理论模型,最后对这一理论模型进行了实验验证。仿真结果表明:声波频率对各种波速几乎没有影响,但各种波速随孔隙度的变化较大;通过实验验证:Scholte波波速理论值与实际测量值之间的相对误差不超过0.24%。     

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.     

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.     

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 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.     

Two different small size testing technique for the creep residual life evaluation of an equiaxed Ni-based superalloy for GT blade and vane application

Two miniaturised creep testing techniques have been applied to a polycrystalline Ni-based superalloy employed for gas turbine blades and vanes. The test campaign was conducted on as-delivered material to assess both testing techniques and evaluate their effectiveness in determining creep residual life and creep residual properties of operated material. The first technique involved small size uniaxial creep testing, conducted on samples with a diameter lower than 3 mm. The diameter effect was deeply studied and a calibration was done to correlate results from standard creep test pieces with those from small diameter samples. The second technique involved the small punch creep test, and the assessment on virgin material enabled determination of the correlation factor, k_sp, to compare small punch test results with uniaxial creep tests. All the collected results were related to the coarse grain macrostructure that characterises this class of alloys and that increases the scatter of the collected experimental results.     

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.     

Investigation the effect of crystal orientation of nickel‐based single crystal superalloys on bending creep tests

The relationship between the three points bending creep test and the uni‐axial creep test on the single crystal superalloy was investigated by using crystal plasticity slip theory with a three‐dimensional (3D) finite element model. The purpose of the present work is to build the relationship between bending creep and conventional uni‐axial tensile creep in determining crystallographic creep parameters for face centered cubic (FCC) nickel‐based single crystal superalloys. To this aim, the bending creep performed on [001]‐, [011]‐, and [111]‐oriented nickel‐based single crystal superalloys were respectively investigated, and the data was compared with those obtained with uni‐axial tensile creep counterparts. It shows that the determination of crystallographic creep stress exponent is independent of crystallographic orientations, and the results agree reasonably well between bending creep test and uni‐axial tensile creep test. The findings may shed some light on understanding of the crystal structures and its time‐dependent deformation mechanisms with the bending creep method.     

Testing complete and compact radial basis functions for solution of eigenvalue problems using the boundary element method with direct integration

This paper analyses the performance of the main radial basis functions in the formulation of the Boundary Element Method (DIBEM). This is an alternative for solving problems modeled by non-adjoint differential operators, since it transforms domain integrals in boundary integrals using radial basis functions. The solution of eigenvalue problem was chosen to performance evaluation. Natural frequencies are calculated numerically using several radial functions and their accuracy is evaluated by comparison with the available analytical solutions and with the Finite Element Method as well. The standard radial basis functions have presented similar performance to compact radial functions, being even slightly superior.     

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.     

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).     

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.