Home About us Contact
|
Home About us Contact | ||
|
|
||
Material Parameters (material + parameter)
Selected AbstractsFatigue crack initiation life estimation in a steel welded joint by the use of a two-scale damage modelFATIGUE & FRACTURE OF ENGINEERING MATERIALS AND STRUCTURES, Issue 5 2009N. LAUTROU ABSTRACT This work deals with the fatigue behaviour of S355NL steel welded joints classically used in naval structures. The approach suggested here, in order to estimate the fatigue crack initiation life, can be split into two stages. First, stabilized stress,strain cycles are obtained in all points of the welded joint by a finite element analysis, taking constant or variable amplitude loadings into account. This calculation takes account of: base metal elastic,plastic behaviour, variable yield stress based on hardness measurements in various zones of the weld, local geometry at the weld toe and residual stresses if any. Second, if a fast elastic shakedown occurs, a two-scale damage model based on Lemaitre et al.'s work is used as a post-processor in order to estimate the fatigue crack initiation life. Material parameters for this model were identified from two Wöhler curves established for base metal. As a validation, four-point bending fatigue tests were carried out on welded specimens supplied by ,DCNS company'. Two load ratios were considered: 0.1 and 0.3. Residual stress measurements by X-ray diffraction completed this analysis. Comparisons between experimental and calculated fatigue lives are promising for the considered loadings. An exploitation of this method is planned for another welding process. [source] Inverse estimation of material properties for sheet metalsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 2 2004K. M. Zhao The main objective of this paper is to estimate the material properties for sheet metals subjected to loading and reverse loading by using an inverse method. Cyclic three-point bending tests are conducted. Bending moments are computed from the measured data, namely, punch stroke, punch load, bending strain and bending angle. Bending moments are also calculated based on the selected constitutive model. Normal anisotropy and non-linear isotropic/kinematic hardening are considered. Material parameters are estimated by minimizing the difference between these two bending moments. Modified Levenberg,Marquardt method is used in the optimization procedure. Stress,strain curves are generated with the material parameters found in this way. Copyright © 2004 John Wiley & Sons, Ltd. [source] Numerical modelling of hydro-mechanical behaviour of collapsible soilsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 12 2008L. M. Costa Abstract This paper presents a numerical simulation of the response of unsaturated soil under different loads, with special attention devoted to the phenomenon of pore collapse during wetting. A coupled hydro-mechanical formulation, implemented in the computational code CODE_BRIGHT, is first presented. Then the model is applied to assess material parameter from laboratory tests and to simulate, as boundary value problems, new and appealing in situ tests designed in the Federal University of Pernambuco. The soil considered is a collapsible soil located in the semi-arid region of the northeast of Brazil. Two situations are analysed: (i) the soil is initially at its natural water content and then flooded at some determined stress level, (ii) the soil is previously flooded and then loaded. Comparison between output results of numerical simulations and experimental data shows a very good agreement, which validate to some extent the proposed experimental procedure and the model formulation. Copyright © 2007 John Wiley & Sons, Ltd. [source] Novel broadband asymptotic waveform evaluation method and its applications in parameter extractionINTERNATIONAL JOURNAL OF RF AND MICROWAVE COMPUTER-AIDED ENGINEERING, Issue 3 2007Guo-Bing Han Abstract A new technique to estimate the effective bandwidth of the asymptotic waveform evaluation (AWE) method is proposed, which guarantees the reliable application of AWE in a certain region. Furthermore, based on the effective bandwidth estimation, we develop a novel broadband AWE method, which is an extension and improvement of the conventional AWE method. The new method selects the multiple expansion points automatically by a simple compare algorithm and then implements the AWE process on every expansion point, thereby expanding the variable bandwidth of the AWE greatly. The variable here can be frequency, structure dimension, material parameter, etc. In the end, two examples of parameter extraction of microwave structures are given, which demonstrate that the new method has high efficiency and good accuracy in a broad band. © 2007 Wiley Periodicals, Inc. Int J RF and Microwave CAE, 2007. [source] Ein Trocknungskoeffizient für BaustoffeBAUPHYSIK, Issue 3 2009Gregor A. Scheffler Dr.-Ing. Berechnungsverfahren; Feuchte Wärme; Versuche Abstract Ein wesentliches Element der hygrothermischen Charakterisierung von Baustoffen ist der Trocknungsversuch. Im Gegensatz zu anderen Feuchtetransportexperimenten wie dem Diffusionsversuch oder dem Wasseraufnahmeexperiment ist es bislang nicht möglich, aus der Trocknung einen einfachen Kennwert abzuleiten. In vielen Fällen, beispielsweise in der Interaktion von Forschung und Industrie, aber auch beim praktischen Vergleich bzw. der Auswahl geeigneter Baustoffe wäre ein solcher Kennwert jedoch wünschenswert. Im vorliegenden Artikel wird zunächst die Bedeutung des Trocknungsversuches für die hygrische Charakterisierung von Baustoffen herausgestellt, aus der sich das Bestreben ableitet, das Trocknungsverhalten zu standardisieren und einen Einzahlen-Materialkennwert zu definieren. Nach einer die verschiedenen Einflussfaktoren der Trocknung differenzierenden Einleitung werden bestehende Ansätze für die Standardisierung des Trocknungsverlaufes bzw. die Ableitung eines Trocknungskoeffizienten vorgestellt. Die einhergehenden Probleme werden diskutiert und weitere Möglichkeiten evaluiert. Ein einfacher Trocknungskoeffizient, der sich aus dem Trocknungsverlauf ableiten lässt, wird definiert. Die Korrelation dieses Koeffizienten mit dem Wasseraufnahmekoeffizienten und dem Dampfdiffusionswiderstand wird analysiert. Sein zusätzlicher Informationsgehalt wird in diesem Zusammenhang kritisch hinterfragt. Im Ergebnis steht die Definition des Trocknungskoeffizienten als ein neuer, unabhängiger Materialkennwert, der die Feuchtetransporteigenschaften im Übergang zwischen hygroskopischem und gesättigtem Transport beschreibt. Mit diesem Kennwert ist es möglich, Baustoffe einfach und schnell hinsichtlich ihres Trocknungsverhaltens zu unterscheiden und zu beurteilen, was insbesondere bei feuchtesensitiven Materialien von Bedeutung ist. A drying coefficient for building materials. The drying experiment is an important element of the hygrothermal characterisation of building materials. Contrary to other moisture transport experiments as the vapour diffusion and the water absorption test, it is until now not possible to derive a simple coefficient for the drying. However, in many cases such a coefficient would be highly appreciated, e.g. in interaction of industry and research or for the distinction and selection of suitable building materials throughout design and practise. This article first highlights the importance of drying experiments for hygrothermal characterisation of building materials on which the attempt is based to standardize the drying experiment as well as to derive a single number material coefficient. The drying itself is briefly reviewed and existing approaches are discussed. On this basis, possible definitions are evaluated. Finally, a drying coefficient is defined which can be determined based on measured drying data. The correlation of this coefficient with the water absorption and the vapour diffusion coefficient is analyzed and its additional information content is critically challenged. As result, a drying coefficient has been derived and defined as a new and independent material parameter. It contains information about the moisture transport properties throughout the wide range of moisture contents from hygroscopic up to saturation. With this new and valuable coefficient, it is now possible to distinguish and select building materials quickly and easily by means of their drying behaviour. This is particularly important for moisture sensitive materials. [source] Reliability Aspects of Microsystems for Automotive Applications,ADVANCED ENGINEERING MATERIALS, Issue 4 2009Roland Müller-Fiedler Abstract The implementation of microsystems in automotive applications is certainly one of the driving forces for the success of MEMS as an industrial technology on mass production level. In many cases, automotive systems based on microsensors are critical to safety. Consequently, microsystems have to assure an accurate, reliable, and failsafe operation during the entire lifetime of the vehicle. Since, the car represents a harsh environment for electronic or mechanical systems and components, reliability issues of MEMS have moved more and more into the focus of research and development. In particular, reliability aspects related to packaging and assembly have become a key issue in lifetime investigations. The packaging of microsystems comprises a variety of materials and material combinations, that directly affect the stability of MEMS components. Therefore, sophisticated characterization methods are needed to extract the reliability-relevant material parameters. This paper gives an introduction into investigations of the stability assessment of glass frit bonding as well as new bonding technologies based on metallic sealing of MEMS devices. [source] Determination of the Gurson,Tvergaard damage model parameters for simulating small punch testsFATIGUE & FRACTURE OF ENGINEERING MATERIALS AND STRUCTURES, Issue 11 2010I. I. CUESTA ABSTRACT The objective of the final small punch test (SPT) is to determine the fracture properties of materials, such as fracture toughness, when not enough material is available for the conduct of conventional fracture tests. The damage model developed by Gurson, and subsequently modified by Tvergaard and Needleman (GTN), allows for the numerical simulation of the elastic-plastic behaviour until fracture. This model is based on several constitutive material parameters that must be calibrated if the model is to be properly applied. In this paper, we develop a consistent methodology for the identification of the GTN damage parameters based on the adjustment of the load-displacement curve obtained in the SPTs. The methodology presented is applicable to simulating other different SPTs with different thicknesses and test temperatures. Also, the three-dimensional modelling developed will be useful in the future for analysing the possible anisotropy exhibited by some materials. The next step in the simulation will be to determine its validity in other stress fields with different triaxiality ratios, like the one present in CT specimens, the ultimate goal being to allow for the estimation of the material fracture toughness. [source] Determining the life cycle of bolts using a local approach and the Dang Van criterionFATIGUE & FRACTURE OF ENGINEERING MATERIALS AND STRUCTURES, Issue 8 2006Y. FARES ABSTRACT The fatigue behaviour of bolts under axial load has always been considered from the component point of view for which fatigue limit is usually taken equal to 50 MPa, and few results are available to designers for limited lifetimes. Here, we take up this problem from a material point of view using a local approach. For each case of fatigue testing, using finite-element (FE) model of the bolt, we determine the stabilized local stress at the root of the first thread in contact with the nut. To characterize bolt behaviour with these numerical results, we use Dang Van multiaxial fatigue criterion for which we extend application to the medium fatigue life. These results can be correlated with the experimental numbers of cycles to failure to determine material parameters of the generalized criterion. Using statistical Gauss method, we can make lifetime predictions for any level of risk of failure. In addition, we propose an analytical model to rapidly determine the local stress condition from nominal loading data (mean stress and alternating stress). This model dispenses us from a new modelling if the bolt is stressed in the same manner as the bolts used for behaviour characterization. Using this model and the generalized criterion, it is extremely easy to make lifetime predictions whatever the risk considered. [source] Cohesive-zone modelling of the deformation and fracture of spot-welded jointsFATIGUE & FRACTURE OF ENGINEERING MATERIALS AND STRUCTURES, Issue 10 2005M. N. CAVALLI ABSTRACT The deformation and failure of spot-welded joints have been successfully modelled using a cohesive-zone model for fracture. This has been accomplished by implementing a user-defined, three-dimensional, cohesive-zone element within a commercial finite-element package. The model requires two material parameters for each mode of deformation. Results show that the material parameters from this type of approach are transferable for identical spot welds in different geometries where a single parameter (such as maximum stress) is not. The approach has been demonstrated using a model system consisting of spot-welded joints made from 5754 aluminium sheets. The techniques for determining the cohesive fracture parameters for both nugget fracture and nugget pullout are described in this paper. It has been demonstrated that once the appropriate cohesive parameters for a weld are determined, quantitative predictions can be developed for the strengths, deformations and failure mechanisms of different geometries with nominally identical welds. [source] Approximation methods for reliability-based design optimization problemsGAMM - MITTEILUNGEN, Issue 2 2007Irfan Kaymaz Abstract Deterministic optimum designs are obtained without considering of uncertainties related to the problem parameters such as material parameters (yield stress, allowable stresses, moment capacities, etc.), external loadings, manufacturing errors, tolerances, cost functions, which could lead to unreliable designs, therefore several methods have been developed to treat uncertainties in engineering analysis and, more recently, to carry out design optimization with the additional requirement of reliability, which referred to as reliability-based design optimization. In this paper, two most common approaches for reliability-based design optimization are reviewed, one of which is reliability-index based approach and the other performancemeasure approach. Although both approaches can be used to evaluate the probabilistic constraint, their use can be prohibitive when the associated function evaluation required by the probabilistic constraint is expensive, especially for real engineering problems. Therefore, an adaptive response surface method is proposed by which the probabilistic constraint is replaced with a simple polynomial function, thus the computational time can be reduced significantly as presented in the example given in this paper. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim) [source] Towards non-linear inversion for characterization of time-lapse phenomena through numerical modellingGEOPHYSICAL PROSPECTING, Issue 4 2003A. Abubakar ABSTRACT We compare two geophysical survey measurements of the same type made at different times in order to characterize the change in the geological medium during the elapsed time. The aim of this study is to develop a strategy using a full non-linear inversion algorithm as the interpretation tool. In this way, not only the location and the form of the changes are recovered, but also the changes in the material parameters of the geological medium can be estimated. In order to solve this fully non-linear problem, the so-called ,multiplicative regularized contrast source inversion' (MR-CSI) method is employed. The unique property of this iterative method is that it does not solve the forward problem at each iterative step. This makes it possible to use the non-linear inversion algorithm for large-scale computation problems. The numerical results show that by taking into account the non-linear nature of the problem, interpretation of the time-lapse data can be significantly improved, compared with that obtained using linear inversion. [source] Cover Picture: Colloidal Synthesis of Hollow Cobalt Sulfide Nanocrystals (Adv. Funct.ADVANCED FUNCTIONAL MATERIALS, Issue 11 2006Mater. Abstract Hollow nanocrystals have been synthesized through a mechanism analogous to the Kirkendall Effect. When a cobalt nanocrystal reacts with sulfur in solution, the outward diffusion of cobalt atoms is faster than the inward diffusion of sulfur atoms through the sulfide shell. The dominating outward diffusion of cobalt cations produces vacancies that can condense into a single void in the center of the nanocrystal at high temperatures. This process provides a general route to the synthesis of hollow nanostructures of a large number of compounds and is described in the Full Paper by A.,P. Alivisatos and co-workers on p.,1389. Formation of cobalt sulfide hollow nanocrystals through a mechanism similar to the Kirkendall Effect has been investigated in detail. It is found that performing the reaction at >,120,°C leads to fast formation of a single void inside each shell, whereas at room temperature multiple voids are formed within each shell, which can be attributed to strongly temperature-dependent diffusivities for vacancies. The void formation process is dominated by outward diffusion of cobalt cations; still, the occurrence of significant inward transport of sulfur anions can be inferred as the final voids are smaller in diameter than the original cobalt nanocrystals. Comparison of volume distributions for initial and final nanostructures indicates excess apparent volume in shells, implying significant porosity and/or a defective structure. Indirect evidence for fracture of shells during growth at lower temperatures was observed in shell-size statistics and transmission electron microscopy images of as-grown shells. An idealized model of the diffusional process imposes two minimal requirements on material parameters for shell growth to be obtainable within a specific synthetic system. [source] Probabilistic yielding and cyclic behavior of geomaterialsINTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 15 2010Kallol Sett Abstract In this paper, the novel concept of probabilistic yielding is used for 1-D cyclic simulation of the constitutive behavior of geomaterials. Fokker,Planck,Kolmogorov equation-based probabilistic elastic,plastic constitutive framework is applied for obtaining the complete probabilistic (probability density function) material response. Both perfectly plastic and hardening-type material models are considered. It is shown that when uncertainties in material parameters are taken into consideration, even the simple, elastic-perfectly plastic model captures some of the important features of geomaterial behavior, for example, modulus reduction with cyclic strain, which, deterministically, is only possible with more advanced constitutive models. Furthermore, it is also shown that the use of isotropic and kinematic hardening rules does not significantly improve the probabilistic material response. Copyright © 2010 John Wiley & Sons, Ltd. [source] Micromechanical viscoelasto-plastic models and finite element implementation for rate-independent and rate-dependent permanent deformation of stone-based materialsINTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 13 2010Qingli Dai Abstract This paper presents parallel and serial viscoelasto-plastic models to simulate the rate-independent and the rate-dependent permanent deformation of stone-based materials, respectively. The generalized Maxwell viscoelastic and Chaboche's plastic models were employed to formulate the proposed parallel and serial viscoelasto-plastic constitutive laws. The finite element (FE) implementation of the parallel model used a displacement-based incremental formulation for the viscoelastic part and an elastic predictor,plastic corrector scheme for the elastoplastic component. The FE framework of the serial viscoelasto-plastic model employed a viscoelastic predictor,plastic corrector algorithm. The stone-based materials are consisted of irregular aggregates, matrix and air voids. This study used asphalt mixtures as an example. A digital sample was generated with imaging analysis from an optically scanned surface image of an asphalt mixture specimen. The modeling scheme employed continuum elements to mesh the effective matrix, and rigid bodies for aggregates. The ABAQUS user material subroutines defined with the proposed viscoelasto-plastic matrix models were employed. The micromechanical FE simulations were conducted on the digital mixture sample with the viscoelasto-plastic matrix models. The simulation results showed that the serial viscoelasto-plastic matrix model generated more permanent deformation than the parallel one by using the identical material parameters and displacement loadings. The effect of loading rates on the material viscoelastic and viscoelasto-plastic mixture behaviors was investigated. Permanent deformations under cyclic loadings were determined with FE simulations. The comparison studies showed that the simulation results correctly predicted the rate-independent and rate-dependent viscoelasto-plastic constitutive properties of the proposed matrix models. Overall, these studies indicated that the developed micromechanical FE models have the abilities to predict the global viscoelasto-plastic behaviors of the stone-based materials. Copyright © 2009 John Wiley & Sons, Ltd. [source] Micromechanical analysis of failure propagation in frictional granular materialsINTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 15 2009Antoinette Tordesillas Abstract The extent to which the evolution of instabilities and failure across multiple length scales can be reproduced with the aid of a bifurcation analysis is examined. We adopt an elastoplastic micropolar constitutive model, recently developed for dense cohesionless granular materials within the framework of thermomicromechanics. The internal variables and their evolution laws are conceived from a direct consideration of the dissipative mechanism of force chain buckling. The resulting constitutive law is cast entirely in terms of the particle scale properties. It thus presents a unique opportunity to test the potential of micromechanical continuum formulations to reproduce key stages in the deformation history: the development of material instabilities and failure following an initially homogeneous deformation. Progression of failure, initiating from frictional sliding and rolling at contacts, followed by the buckling of force chains, through to macroscopic strain softening and shear banding, is reproduced. Bifurcation point, marking the onset of shear banding, occurred shortly after the peak stress ratio. A wide range of material parameters was examined to show the effect of particle scale properties on the progression of failure. Model predictions on the thickness and angle of inclination of the shear band and the structural evolution inside the band, namely the latitudinal distribution of particle rotations and the angular distributions of contacts and the normal contact forces, are consistent with observations from numerical simulations and experiments. Copyright © 2009 John Wiley & Sons, Ltd. [source] Vertical dynamic response of a rigid foundation embedded in a poroelastic soil layerINTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 11 2009Y. Q. Cai Abstract A simplified analytical method is presented for the vertical dynamic analysis of a rigid, massive, cylindrical foundation embedded in a poroelastic soil layer. The foundation is subjected to a time-harmonic vertical loading and is perfectly bonded to the surrounding soil in the vertical direction. The soil underlying the foundation base is represented by a single-layered poroelastic soil based on rigid bedrock while the soil at the side of the foundation is modeled as an independent poroelastic layer composed of a series of infinitesimally thin layers. The behavior of the soil is governed by Biot's poroelastodynamic theory and its governing equations are solved by the use of Hankel integral transform. The contact surface between the foundation base and the soil is smooth and fully permeable. The dynamic interaction problem is solved following standard numerical procedures. The accuracy of the present solution is verified by comparisons with the well-known solutions obtained from other approaches for both the elastodynamic interaction problem and poroelastodynamic interaction problem. Numerical results for the vertical dynamic impedance and response factor of the foundation are presented to demonstrate the influence of nondimensional frequency of excitation, soil layer thickness, poroelastic material parameters, depth ratio and mass ratio on the dynamic response of a rigid foundation embedded in a poroelastic soil layer. Copyright © 2008 John Wiley & Sons, Ltd. [source] Finite elements modelling of the long-term behaviour of a full-scale flexible pavement with the shakedown theoryINTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 1 2009Cyrille Chazallon Abstract Rutting, due to permanent deformations of unbound materials, is one of the principal damage modes of low-traffic pavements. Flexible pavement design methods remain empirical; they do not take into account the inelastic behaviour of pavement materials and do not predict the rutting under cyclic loading. A simplified method, based on the concept of the shakedown theory developed by Zarka for metallic structures under cyclic loadings, has been used to estimate the permanent deformations of unbound granular materials subjected to traffic loading. Based on repeated load triaxial tests, a general procedure has been developed for the determination of the material parameters of the constitutive model. Finally, the results of a finite elements modelling of the long-term behaviour of a flexible pavement with the simplified method are presented and compared with the results of a full-scale flexible pavement experiment performed by Laboratoire Central des Ponts et Chaussées. Finally, the calculation of the rut depth evolution with time is carried out. Copyright © 2008 John Wiley & Sons, Ltd. [source] Frost heave modelling using porosity rate functionINTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 8 2006Radoslaw L. Michalowski Abstract Frost-susceptible soils are characterized by their sensitivity to freezing that is manifested in heaving of the ground surface. While significant contributions to explaining the nature of frost heave in soils were published in late 1920s, modelling efforts did not start until decades later. Several models describing the heaving process have been developed in the past, but none of them has been generally accepted as a tool in engineering applications. The approach explored in this paper is based on the concept of the porosity rate function dependent on two primary material parameters: the maximum rate, and the temperature at which the maximum rate occurs. The porosity rate is indicative of ice growth, and this growth is also dependent on the temperature gradient and the stress state in the freezing soil. The advantage of this approach over earlier models stems from a formulation consistent with continuum mechanics that makes it possible to generalize the model to arbitrary three-dimensional processes, and use the standard numerical techniques in solving boundary value problems. The physical premise for the model is discussed first, and the development of the constitutive model is outlined. The model is implemented in a 2-D finite element code, and the porosity rate function is calibrated and validated. Effectiveness of the model is then illustrated in an example of freezing of a vertical cut in frost-susceptible soil. Copyright © 2006 John Wiley & Sons, Ltd. [source] Parameter sensitivity in finite element analysis with constitutive models of the rate typeINTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 2 2006Wolfgang Fellin Abstract Non-linear soil,structure interactions are usually analysed using an incremental finite element approach. There a constitutive subroutine provides for each element the stress increment for a given strain increment. In geotechnical calculations, uncertainties in material parameters and initial conditions are abundant. Sensitivity analysis can be a first step to account for such uncertainties. Sensitivities of the system response with respect to material parameters and initial conditions can be calculated by differentiating the whole numerical scheme. It turns out that the essential information from the constitutive subroutine are the derivatives of the stress increment with respect to the strain increment, as well as the derivatives with respect to material parameters and all state variables involved in the problem. We propose a method to compute these quantities numerically for any constitutive model that can be written in rate form and for any suitable integrator of such a model. We further present a concise way to supply the output of the sensitivity analysis to the designing engineer. Our theoretical investigations are illustrated with element tests and with a typical geotechnical application. Copyright © 2005 John Wiley & Sons, Ltd. [source] Inverse estimation of material properties for sheet metalsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 2 2004K. M. Zhao The main objective of this paper is to estimate the material properties for sheet metals subjected to loading and reverse loading by using an inverse method. Cyclic three-point bending tests are conducted. Bending moments are computed from the measured data, namely, punch stroke, punch load, bending strain and bending angle. Bending moments are also calculated based on the selected constitutive model. Normal anisotropy and non-linear isotropic/kinematic hardening are considered. Material parameters are estimated by minimizing the difference between these two bending moments. Modified Levenberg,Marquardt method is used in the optimization procedure. Stress,strain curves are generated with the material parameters found in this way. Copyright © 2004 John Wiley & Sons, Ltd. [source] Frictional granular mechanics: A variational approachINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 10 2010R. Holtzman Abstract The mechanical properties of a cohesionless granular material are evaluated from grain-scale simulations. Intergranular interactions, including friction and sliding, are modeled by a set of contact rules based on the theories of Hertz, Mindlin, and Deresiewicz. A computer-generated, three-dimensional, irregular pack of spherical grains is loaded by incremental displacement of its boundaries. Deformation is described by a sequence of static equilibrium configurations of the pack. A variational approach is employed to find the equilibrium configurations by minimizing the total work against the intergranular loads. Effective elastic moduli are evaluated from the intergranular forces and the deformation of the pack. Good agreement between the computed and measured moduli, achieved with no adjustment of material parameters, establishes the physical soundness of the proposed model. Copyright © 2009 John Wiley & Sons, Ltd. [source] Discrete element method for modelling solid and particulate materialsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 4 2007Federico A. Tavarez Abstract The discrete element method (DEM) is developed in this study as a general and robust technique for unified two-dimensional modelling of the mechanical behaviour of solid and particulate materials, including the transition from solid phase to particulate phase. Inter-element parameters (contact stiffnesses and failure criteria) are theoretically established as functions of element size and commonly accepted material parameters including Young's modulus, Poisson's ratio, ultimate tensile strength, and fracture toughness. A main feature of such an approach is that it promises to provide convergence with refinement of a DEM discretization. Regarding contact failure, an energy criterion based on the material's ultimate tensile strength and fracture toughness is developed to limit the maximum contact forces and inter-element relative displacement. This paper also addresses the issue of numerical stability in DEM computations and provides a theoretical method for the determination of a stable time-step. The method developed herein is validated by modelling several test problems having analytic solutions and results show that indeed convergence is obtained. Moreover, a very good agreement with the theoretical results is obtained in both elastic behaviour and fracture. An example application of the method to high-speed penetration of a concrete beam is also given. Copyright © 2006 John Wiley & Sons, Ltd. [source] A bounding surface plasticity model for cyclic loading of granular soilsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 14 2005N. Khalili Abstract A constitutive model for describing the stress,strain behaviour of granular soils subjected to cyclic loading is presented. The model is formulated using bounding surface theory within a critical state framework. A single set of material parameters is introduced for the complete characterization of the constitutive model. The shape of the bounding surface is based on experimental observations of undrained stress paths for loose samples. A mapping rule which passes through stress reversal points is introduced to depict the stress,strain behaviour during unloading and reloading. The effect of particle crushing is considered through a modified critical state line. Essential features of the model are validated using several experimental data from the literature. Both drained and undrained loading conditions are considered. The characteristic features of behaviour in granular soils subjected to cyclic loading are captured. Copyright © 2005 John Wiley & Sons, Ltd. [source] A conservative integral for bimaterial notches subjected to thermal stressesINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 6 2004Leslie Banks-Sills Abstract In this investigation, a conservative integral based on the Betti reciprocal principle is developed to obtain stress intensity factors for a bimaterial notch in which the body is subjected to a thermal load. The bonded materials are linear elastic, isotropic and homogeneous. According to the linear theory of elasticity, stresses in the neighbourhood of the notch tip are generally singular as a result of the mismatch of the elastic constants. Eigenvalues and eigenfunctions depend upon the mechanical properties and wedge angles. They may be real, complex or power-logarithmic. Real and complex eigenvalues are considered in this study. The stress intensity factor represents the amplitude of the stress singularity and depends upon material properties, geometry and load or temperature. Because of the highly singular behaviour of one of the integrals that is part of the conservative integral, the former is carried out by a hybrid analytical/numerical scheme. The finite element method is employed to obtain displacements caused by the temperature distribution in the body. The conservative integral is applied to several problems appearing in the literature. Both good agreement between those results and the ones obtained here, as well as path stability for all problems is attained. A wide range of material parameters is also studied. Copyright © 2004 John Wiley & Sons, Ltd. [source] Numerical study of consistency of rate constitutive equations with elasticity at finite deformationINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 9 2002Ruocheng Lin Abstract The present work is concerned with the numerical study of the elasticity consistency of the spatial rate equations using the conventional Oldroyd, Truesdell, Cotter,Rivlin, Jaumann and Green,Naghdi rates and the three novel co-rotational ,E - and ,¯L -based, logarithmic rates, and of the rotated material rate equation describing the relationship between the material time derivatives of the rotated Kirchhoff stress and material logarithmic strain. To this end, three integration procedures for updating stress are presented. The stress responses of several typical deformation processes are simulated. According to the numerical results we know that among the spatial rate equations only the logarithmic rate equation is consistent with elasticity under constant material parameters. Integrating the other spatial rate equations will provide path-dependent stress response. These numerical conclusions support the arguments in H. Xiao et al. (Acta Mechanica 1999; 138:31,50). The reasons leading to elasticity inconsistency of spatial rate equations are analysed. If the material parameters are assumed to be strain-dependent, the logarithmic rate equation loses also its elasticity-consistent property. The numerical results prove also that the spatial logarithmic and rotated material rate equations are equivalent to each other. Copyright © 2002 John Wiley & Sons, Ltd. [source] Numerical simulation of bubble and droplet deformation by a level set approach with surface tension in three dimensionsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 9 2010Roberto Croce Abstract In this paper we present a three-dimensional Navier,Stokes solver for incompressible two-phase flow problems with surface tension and apply the proposed scheme to the simulation of bubble and droplet deformation. One of the main concerns of this study is the impact of surface tension and its discretization on the overall convergence behavior and conservation properties. Our approach employs a standard finite difference/finite volume discretization on uniform Cartesian staggered grids and uses Chorin's projection approach. The free surface between the two fluid phases is tracked with a level set (LS) technique. Here, the interface conditions are implicitly incorporated into the momentum equations by the continuum surface force method. Surface tension is evaluated using a smoothed delta function and a third-order interpolation. The problem of mass conservation for the two phases is treated by a reinitialization of the LS function employing a regularized signum function and a global fixed point iteration. All convective terms are discretized by a WENO scheme of fifth order. Altogether, our approach exhibits a second-order convergence away from the free surface. The discretization of surface tension requires a smoothing scheme near the free surface, which leads to a first-order convergence in the smoothing region. We discuss the details of the proposed numerical scheme and present the results of several numerical experiments concerning mass conservation, convergence of curvature, and the application of our solver to the simulation of two rising bubble problems, one with small and one with large jumps in material parameters, and the simulation of a droplet deformation due to a shear flow in three space dimensions. Furthermore, we compare our three-dimensional results with those of quasi-two-dimensional and two-dimensional simulations. This comparison clearly shows the need for full three-dimensional simulations of droplet and bubble deformation to capture the correct physical behavior. Copyright © 2009 John Wiley & Sons, Ltd. [source] Numerical simulation of a single bubble by compressible two-phase fluidsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 6 2010Siegfried Müller Abstract The present work deals with the numerical investigation of a collapsing bubble in a liquid,gas fluid, which is modeled as a single compressible medium. The medium is characterized by the stiffened gas law using different material parameters for the two phases. For the discretization of the stiffened gas model, the approach of Saurel and Abgrall is employed where the flow equations, here the Euler equations, for the conserved quantities are approximated by a finite volume scheme, and an upwind discretization is used for the non-conservative transport equations of the pressure law coefficients. The original first-order discretization is extended to higher order applying second-order ENO reconstruction to the primitive variables. The derivation of the non-conservative upwind discretization for the phase indicator, here the gas fraction, is presented for arbitrary unstructured grids. The efficiency of the numerical scheme is significantly improved by employing local grid adaptation. For this purpose, multiscale-based grid adaptation is used in combination with a multilevel time stepping strategy to avoid small time steps for coarse cells. The resulting numerical scheme is then applied to the numerical investigation of the 2-D axisymmetric collapse of a gas bubble in a free flow field and near to a rigid wall. The numerical investigation predicts physical features such as bubble collapse, bubble splitting and the formation of a liquid jet that can be observed in experiments with laser-induced cavitation bubbles. Opposite to the experiments, the computations reveal insight to the state inside the bubble clearly indicating that these features are caused by the acceleration of the gas due to shock wave focusing and reflection as well as wave interaction processes. While incompressible models have been used to provide useful predictions on the change of the bubble shape of a collapsing bubble near a solid boundary, we wish to study the effects of shock wave emissions into the ambient liquid on the bubble collapse, a phenomenon that may not be captured using an incompressible fluid model. Copyright © 2009 John Wiley & Sons, Ltd. [source] Effects of microcompounding process parameters on the properties of ABS/polyamide-6 blends based nanocompositesJOURNAL OF APPLIED POLYMER SCIENCE, Issue 5 2008Güralp Özkoç Abstract Melt intercalation method was applied to produce acrylonitrile-butadiene-styrene/polyamide-6 (ABS/PA6) blends based organoclay nanocomposites using a conical twin-screw microcompounder. The blend was compatibilized using a maleated olefinic copolymer. The effects of microcompounding conditions such as screw speed, screw rotation-mode (co- or counter-), and material parameters such as blend composition and clay loading level on the morphology of the blends, dispersibility of nanoparticles, and mechanical properties were investigated. Furthermore, corotating screws were modified to achieve elongational flow which is efficient for obtaining dispersive mixing. The morphology was examined by SEM analysis after preferential extraction of the minor phase. Subsequently, the SEM micrographs were quantitatively analyzed using image analyzer software. The morphology of the blends indicated that processing with counter-rotation at a given screw speed yielded coarser morphology than that of processed with corotation. X-ray diffraction analysis showed that highest level of exfoliation is observed with increasing PA6 content, at 200 rpm of screw speed and in corotation mode. Also, the effects of screw speed, screw rotation mode, and screw modification were discussed in terms of XRD responses of the nanocomposites. The aspect ratio of the clay particles which were measured by performing image analysis on TEM micrographs exhibited a variation with processing conditions and they are in accordance with the modulus of the nanocomposites. © 2007 Wiley Periodicals, Inc. J Appl Polym Sci, 2008 [source] Experimental and model determination of human intervertebral disc osmoviscoelasticityJOURNAL OF ORTHOPAEDIC RESEARCH, Issue 8 2008Y. Schroeder Abstract Finite element (FE) models have become an important tool to study load distribution in the healthy and degenerated disc. However, model predictions require accurate constitutive laws and material properties. As the mechanical properties of the intervertebral disc are regulated by its biochemical composition and fiber-reinforced structure, the relationship between the constitutive behavior of the tissue and its composition requires careful consideration. While numerous studies have investigated the annulus fibrosus compressive and tensile properties, specific conditions required to determine model parameters for the osmoviscoelastic model are unavailable. Therefore, the objectives of this study were (1) to complement the existing material testing in the literature with confined compression and tensile tests on human annulus fibrosus and (2) to use these data, together with existing nucleus pulposus compression data to tune a composition-based, osmoviscoelastic material constitutive law. The osmoviscoelastic material constitutive law and the experimental data were used to describe the fiber and nonfiber properties of the human disc. The compressive material properties of normal disc tissue were Gm,=,1.23 MPa, M,=,1.57, and ,,=,1.964,×,10,16 m4/Ns; the tensile fiber material parameters were E0,=,77.0 MPa; E,,=,500 MPa, and ,,=,1.8,×,103 MPa,s. The goodness of fit ranged from 0.88 to 0.96 for the four experimental conditions evaluated. The constitutive law emphasized the interdependency of the strong swelling ability of the tissue and the viscoelastic nature of the collagen fibers. This is especially important for numerical models to further study the load sharing behavior with regard to disc degeneration and regeneration. © 2008 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 26:1141,1146, 2008 [source] Coarsening of Faceted CrystalsJOURNAL OF THE AMERICAN CERAMIC SOCIETY, Issue 3 2002Gregory S. Rohrer The influence of the nucleation energy barrier on the capillary-driven coarsening of faceted crystals that exchange material by diffusion is quantified. Our calculations are based on the assumption that the transport of material between particles must happen in series with the nucleation of partial layers on flat facets. Using a numerical model based on this idea, we simulate the time evolution of distributions of crystals that are made up of perfect faceted crystals (without step-producing defects), crystals containing step-producing defects, and mixtures of the two types. We find that the coarsening of a distribution containing only perfect faceted crystals is arrested at a size where the nucleation energy barrier becomes prohibitive. This critical size ranges from a few nanometers to several hundred nanometers, depending on material parameters and experimental conditions. When a small fraction of the crystals have step-producing defects (for these crystals the nucleation energy barrier vanishes), they can grow to large sizes at the expense of the perfect crystals and a bimodal grain size distribution is created. Based on these results, we hypothesize that when abnormal coarsening is observed in nature, it results from the presence of a small number of crystals with step-producing defects. [source] | |||||||||||||||||||||||||||||||||||||