Numerical study of nonlinear interaction between a crack and elastic waves under an oblique incidence

A Finite Element (FE) model is proposed to study the interaction between in-plane elastic waves and a crack of different orientations. The crack is modeled by an interface of unilateral contact with Coulombs friction. These contact laws are modified to take into account a pre-stress _σ0${\sigma }_0$ that closes the crack. Using the FE model, it is possible to obtain the contact stresses during wave propagation. These contact stresses provide a better understanding of the coupling between the normal and tangential behavior under oblique incidence, and explain the generation of higher harmonics. This new approach is used to analyze the evolution of the higher harmonics obtained as a function of the angle of incidence, and also as a function of the excitation level. The pre-stress condition is a governing parameter that directly changes the nonlinear phenomenon at work at the interface and therefore the harmonic generation. The diffracted fields obtained by the nonlinear and linear models are also compared.     

Acoustic wave propagation in inhomogeneous,layered waveguides based on modal expansions and hp-FEM

A new model is presented for harmonic wave propagation and scattering problems in non-uniform, stratified waveguides, governed by the Helmholtz equation. The method is based on a modal expansion, obtained by utilizing cross-section basis defined through the solution of vertical eigenvalue problems along the waveguide. The latter local basis is enhanced by including additional modes accounting for the effects of inhomogeneous boundaries and/or interfaces. The additional modes provide implicit summation of the slowly convergent part of the local-mode series, rendering the remaining part to be fast convergent, increasing the efficiency of the method, especially in long-range propagation applications. Using the enhanced representation, in conjunction with an energy-type variational principle, a coupled-mode system of equations is derived for the determination of the unknown modal-amplitude functions. In the case of multilayered environments, h$h$- and p$p$-FEM have been applied for the solution of both the local vertical eigenvalue problems and the resulting coupled mode system, exhibiting robustness and good rates of convergence. Numerical examples are presented in simple acoustic propagation problems, illustrating the role and significance of the additional mode(s) and the efficiency of the present model, that can be naturally extended to treat propagation and scattering problems in more complex 3D waveguides.     

Numerical study of contact forces for crack closure analysis

Plasticity induced crack closure (PICC) has been widely studied using numerical models. Different numerical parameters can be considered to quantify the opening level, namely one based on the analysis of contact stresses at minimum load. A modified version of this parameter is proposed here, based on nodal contact forces instead of contact stresses. The predictions were found to be similar to those obtained from the contact status of 2nd node behind crack tip. The PICC_contact parameter was also found to be very consistent and adequate for parametric studies of the influence of different physical parameters. The contributions to the opening stress intensity factor of different points along crack flank were found to strongly decrease with distance to crack tip. The cumulative K_open between the crack tip and a distance of 0.1 mm was found to vary from 30% to 100%, increasing with stress ratio, R. Finally, a K solution was developed for punctual forces applied on crack flank and compared with a literature solution for infinite plates. A good agreement was found for plane strain state but significant differences of about 10% were found for plane stress state.     

Acoustic radiation induced by bubble motion in compressible fluid

Based on the theory of compressible fluid, a three-dimension boundary element method is utilized to research the motion of bubble. The far-field noise radiation during the growth and contraction is calculated by the Kirchhoff formula and the Ffowcs Williams-Hawkings （FW-H） formula with a fixed radiation surface being arranged at the near-field of bubble as a new acoustic source. The results show that the amplitude of the sound pressure induced by non-spherical bubble is lower than that of spherical bubble in the contraction phase. The retardance effect is more obvious when the observer is farther away from the bubble. In the anaphase of contraction, the observer with the maximum amplitude of sound pressure moves up with the obvious jet. Larger buoyance parameters will generate lower sound pressure amplitudes in the anaphase, while larger intensive parameters will cause higher sound pressure amplitudes in the whole procedure of bubble motion.     

Temperature-dependence of acoustic fatigue life for thermal protection structures

Acoustic fatigue life evaluation is essential for thermal protection structures due to the severe thermo-acoustic load in service. A study on temperaturedependence of acoustic fatigue life for a C/SiC panel is presented in this paper. Effects of temperature on both the structural responses and the S-N curves are investigated. The Dirlik method is adopted to predict the fatigue life of a C/SiC panel at three different temperatures respectively. Significant differences are observed from the results of numerical simulations between the fatigue lives of the panel in the three cases. The temperature-dependence of acoustic fatigue life of a C/SiC panel is verified, and fatigue test of the material needs to be more attentively performed.     

Some basic contact problems in couple stress elasticity

Indentation tests have long been a standard method for material characterization due to the fact that they provide an easy, inexpensive, non-destructive and objective method of evaluating basic properties from small volumes of materials. As the contact scales in such experiments reduce progressively (micro to nano-scales) the internal material lengths become important and their effect upon the macroscopic response cannot be ignored. In the present study, we derive general solutions for three basic two-dimensional (2D) plane-strain contact problems within the framework of the generalized continuum theory of couple-stress elasticity. This theory introduces characteristic material lengths in order to describe the pertinent scale effects that emerge from the underlying microstructure and has proved to be very effective for modeling microstructured materials. By using this theory, we initially study the problem of the indentation of a deformable elastic half-plane by a flat punch, then by a cylindrical indentor, and finally by a shallow wedge indentor. Our approach is based on singular integral equations which have resulted from a treatment of the mixed boundary value problems via integral transforms and generalized functions. The results show significant departure from the predictions of classical elasticity revealing that it is inadequate to analyze indentation problems in microstructured materials employing only classical contact mechanics.     

Development of a machine integrated strain-based contact force sensor for pad foot soil compactors

An investigation was undertaken to explore the use of measurable pad strains on a non-vibratory pad foot roller to provide real time continuous evidence of compaction and contact force. Individual pads were instrumented with strain gages in a pattern chosen based on pad finite element analysis (FEA). Different pad–soil contact stress distributions were modeled to simulate a range of soil conditions. The FEA revealed that the contact stress distribution has a significant influence on the observed pad strain field, suggesting soil specific interpretation of pad strains in order to determine contact force. Results from uniaxial laboratory testing of pad loading on dry sand verified the FEA, i.e., experimental strains generally matched within 15% of FEA strains. The contact stress distribution was measured using tactile pressure sensors and found to be moderately parabolic. A soil specific empirical calibration factor relating vertical sidewall strains to contact force was determined. Field testing was performed on the dry sand with multiple instrumented pads installed on a Caterpillar CP56 roller. Pad strain magnitudes increased up to 250% during compaction from repeated passes of the roller. Using the empirical calibration factor, the estimated contact force was shown to increase with compaction, represented by the independently-measured soil unit weight.     

Parameterization and modelling of large off-road tyres for ride analyses: Part 2 – Parameterization and validation of tyre models

Every mathematical model used in a simulation is an idealization and simplification of reality. Vehicle dynamic simulations that go beyond the fundamental investigations require complex multi-body simulation models. The tyre–road interaction presents one of the biggest challenges in creating an accurate vehicle model. Many tyre models have been proposed and developed but proper validation studies are less accessible. These models were mostly developed and validated for passenger car tyres for application on relatively smooth roads. The improvement of ride comfort, safety and structural integrity of large off-road vehicles, over rough terrain, has become more significant in the development process of heavy vehicles. This paper investigates whether existing tyre models can be used to accurately describe the vertical behaviour of large off road tyres while driving over uneven terrain. [1] Presented an extensive set of experimentally determined parameterization and validation data for a large off-road tyre. Both laboratory and field test are performed for various loads, inflation pressures and terrain inputs. The parameterization process of four tyre models or contact models are discussed in detail. The parameterized models are then validated against test results on various hard but rough off-road terrain and the results are discussed.     

Toward a reduction of mesh imprinting

This paper is the initial investigation into a new Lagrangian cell‐centered hydrodynamic scheme that is motivated by the desire for an algorithm that resists mesh imprinting and has reduced complexity. Key attributes of the new approach include multidimensional construction, the use of flux‐corrected transport (FCT) to achieve second order accuracy, automatic determination of the mesh motion through vertex fluxes, and vorticity control. Toward this end, vorticity preserving Lax–Wendroff (VPLW) type schemes for the two‐dimensional acoustic equations were analyzed and then implemented in a FCT algorithm. Here, mesh imprinting takes the form of anisotropic dispersion relationships. If the stencil for the LW methods is limited to nine points, four free parameters exist. Two parameters were fixed by insisting that no spurious vorticity be created. Dispersion analysis was used to understand how the remaining two parameters could be chosen to increase isotropy. This led to new VPLW schemes that suffer less mesh imprinting than the rotated Richtmyer method. A multidimensional, vorticity preserving FCT implementation was then sought using the most promising VPLW scheme to address the problem of spurious extrema. A well‐behaved first order scheme and a new flux limiter were devised in the process. The flux limiter is unique in that it acts on temporal changes and does not place a priori bounds on the solution. Numerical results have demonstrated that the vorticity preserving FCT scheme has comparable performance to an unsplit MUSCL‐H algorithm at high Courant numbers but with reduced mesh imprinting and superior symmetry preservation. Copyright © 2014 John Wiley & Sons, Ltd.     

一种高阶非线性钢筋混凝土平板单元

在文献［１］高阶钢筋混凝土膜元的基础上，提出了一种具有１４个节点的高阶钢筋混凝土平板单元．钢筋混凝土材料模型仍采用Ｖｅｃｃｈｉｏ的抹平旋转裂缝模型，几何非线性仍采用所谓总体Ｌａｇｒａｎｇｅ列式法，非线性方程采用割线刚度位移增量迭代法．数值算例表明本文的方法是可靠的，高阶平板单元虽然列式复杂，但与低阶元相比，其计算量要少且精度要高