A finite deformation analysis of polycrystal by considering the influence of grain boundary

By combining grain boundary (GB) and its influence zone, a micromechanic model for polycrystal is established for considering the influence of GB. By using the crystal plasticity theory and the finite element method for finite deformation, numerical simulation is carried out by the model. Calculated results display the microscopic characteristic of deformation fields of grains and are in qualitative agreement with experimental results.

Triple-region structure for turbulent-flow in a square duct - a finite-element approach

The Reynolds-averaged Navier-Stokes equations for describing the turbulent flow in a straight square duct are formulated with two different turbulence models. The governing equations are then expanded as a multi-deck structure in a plane perpendicular to the streamwise direction, with each deck characterized by its dominant physical forces as commonly carried out in analytical work using triple-deck expansion. The resulting equations are numerically integrated using higher polynomial (H-P) finite element technique for each cross-sectional plane to be followed by finite difference representation in the streamwise direction until a fully developed state is reached. The computed results using the two different turbulence models show fair agreement with each other, and concur with the vast body of available experimental data. There is also general agreement between our results and the recent numerical works anisotropic k-epsilon turbulence model.

A new formulated method of a quasi-compatible finite-element and its application in fracture-mechanics

Based on the local properties of a singular field, the displacement pattern of an isoparametric element is improved and a new formulated method of a quasi-compatible finite element is proposed in this paper. This method can be used to solve various engineering problems containing singular distribution, especially, the singular field existing at the tip of cracks. The singular quasi-compatible element (SQCE) is constructed. The characteristics of the elements are analysed from various angles and many examples of calculations are performed. The results show that this method has many significant advantages, by which, the numerical analysis of brittle fracture problems can be solved.

Experimental and finite element studies of a large arch dam

Forced vibration field tests and finite element studies have been conducted on Morrow Point (arch) Dam in order to investigate dynamic dam-water interaction and water compressibility. Design of the data acquisition system incorporates several special features to retrieve both amplitude and phase of the response in a low signal to noise environment. These features contributed to the success of the experimental program which, for the first time, produced field evidence of water compressibility; this effect seems to play a significant role only in the symmetric response of Morrow Point Dam in the frequency range examined. In the accompanying analysis, frequency response curves for measured accelerations and water pressures as well as their resonating shapes are compared to predictions from the current state-of-the-art finite element model for which water compressibility is both included and neglected. Calibration of the numerical model employs the antisymmetric response data since they are only slightly affected by water compressibility, and, after calibration, good agreement to the data is obtained whether or not water compressibility is included. In the effort to reproduce the symmetric response data, on which water compressibility has a significant influence, the calibrated model shows better correlation when water compressibility is included, but the agreement is still inadequate. Similar results occur using data obtained previously by others at a low water level. A successful isolation of the fundamental water resonance from the experimental data shows significantly different features from those of the numerical water model, indicating possible inaccuracy in the assumed geometry and/or boundary conditions for the reservoir. However, the investigation does suggest possible directions in which the numerical model can be improved.

Finite-Element Simulations of Earthquakes

This thesis discusses simulations of earthquake ground motions using prescribed ruptures and dynamic failure. Introducing sliding degrees of freedom led to an innovative technique for numerical modeling of earthquake sources. This technique allows efficient implementation of both prescribed ruptures and dynamic failure on an arbitrarily oriented fault surface. Off the fault surface the solution of the three-dimensional, dynamic elasticity equation uses well known finite-element techniques. We employ parallel processing to efficiently compute the ground motions in domains containing millions of degrees of freedom.

Using prescribed ruptures we study the sensitivity of long-period near-source ground motions to five earthquake source parameters for hypothetical events on a strike-slip fault (M_w 7.0 to 7.1) and a thrust fault (M_w 6.6 to 7.0). The directivity of the ruptures creates large displacement and velocity pulses in the ground motions in the forward direction. We found a good match between the severity of the shaking and the shape of the near-source factor from the 1997 Uniform Building Code for strike-slip faults and thrust faults with surface rupture. However, for blind thrust faults the peak displacement and velocities occur up-dip from the region with the peak near-source factor. We assert that a simple modification to the formulation of the near-source factor improves the match between the severity of the ground motion and the shape of the near-source factor.

For simulations with dynamic failure on a strike-slip fault or a thrust fault, we examine what constraints must be imposed on the coefficient of friction to produce realistic ruptures under the application of reasonable shear and normal stress distributions with depth. We found that variation of the coefficient of friction with the shear modulus and the depth produces realistic rupture behavior in both homogeneous and layered half-spaces. Furthermore, we observed a dependence of the rupture speed on the direction of propagation and fluctuations in the rupture speed and slip rate as the rupture encountered changes in the stress field. Including such behavior in prescribed ruptures would yield more realistic ground motions.

Energy stable and momentum conserving hybrid finite element method for the incompressible Navier-Stokes equations

A hybrid method for the incompressible Navier-Stokes equations is presented. The method inherits the attractive stabilizing mechanism of upwinded discontinuous Galerkin methods when momentum advection becomes significant, equal-order interpolations can be used for the velocity and pressure fields, and mass can be conserved locally. Using continuous Lagrange multiplier spaces to enforce flux continuity across cell facets, the number of global degrees of freedom is the same as for a continuous Galerkin method on the same mesh. Different from our earlier investigations on the approach for the Navier-Stokes equations, the pressure field in this work is discontinuous across cell boundaries. It is shown that this leads to very good local mass conservation and, for an appropriate choice of finite element spaces, momentum conservation. Also, a new form of the momentum transport terms for the method is constructed such that global energy stability is guaranteed, even in the absence of a pointwise solenoidal velocity field. Mass conservation, momentum conservation, and global energy stability are proved for the time-continuous case and for a fully discrete scheme. The presented analysis results are supported by a range of numerical simulations. © 2012 Society for Industrial and Applied Mathematics.

Design of brushless doubly-fed machines based: On magnetic circuit modeling

This paper proposes a magnetic circuit model (MCM) for the design of a brushless doubly-fed machine (BDFM). The BDFM possesses advantages in terms of high reliability and reduced gearbox stages, and it requires a fractionally-rated power converter. This makes it suitable for utilization in offshore wind turbines. It is difficult for conventional design methods to calculate the flux in the stator because the two sets of stator windings, which have different pole number, form a complex flux pattern which is not easily determined using common analytical approaches. However, it is advantageous to predict the flux density in the teeth and air-gap at the initial design stage for sizing purposes without recourse finite element analysis. Therefore a magnetic circuit model is developed in this paper to calculate the flux density. A BDFM is used as a case study with FEA validation. © 1965-2012 IEEE.

Numerical modeling of deep soil mixing excavation support

The finite element method (FEM) is growing in popularity over the pressure diagram/hand calculation method for analysis of excavation systems in general and deep soil mixing excavations in particular. In this paper, a finite element analysis is used to study the behavior of a deep mixed excavation. Through the use of Plaxis (a FEM software program), the construction sequence is simulated by following the various construction phases allowing for deflections due to strut or anchor installation to be predicted. The numerical model used in this study simulates the soil cement columns as a continuous wall matching the bending stiffness of the actual wall. Input parameters based on laboratory tests and modeling assumptions are discussed. An example of the approach is illustrated using the Islais Creek Transport/Storage Project in San Francisco, California. Copyright ASCE 2006.

Modeling mode i fracture of bitumen films

The fracture behavior of thin films of bitumen in double cantilever beam (DCB) specimens was investigated over a wide range of temperature and loading rate conditions using finite-element analysis. The model includes a phenomenological model for the mechanical behavior of bitumen, implemented into a special-purpose finite-element user material subroutine, combined with a cohesive zone model (CZM) for simulating the fracture process. The finite-element model is validated against experimental results from laboratory tests of DCB specimens by comparing measured and predicted load-line deflection histories and fracture energy release rates. Computer simulation results agreed well with experimental data of DCB joints containing bitumen films in terms of peak stress, fracture toughness, and stress-strain history response. The predicted "normalized toughness," G=2h, was found to increase in a power-law manner with effective temperaturecompensated strain rate in the ductile region as previously observed experimentally. In the brittle regime, G=2h is virtually constant. The model successfully captured the ductile and brittle failure behavior of bitumen films in opening mode (tension) for stable crack growth conditions. © 2013 American Society of Civil Engineers.

Modeling and experiments on an isothermal fatigue test for solder joints

This paper investigates an isothermal fatigue test for solder joints developed at the NPL. The test specimen is a lap joint between two copper arms. During the test the displacement at the ends of the copper are controlled and the force measured. The modeling results in the paper show that the displacement across the solder joint is not equal to the displacement applied at the end of the specimen. This is due to deformation within the copper arms. A method is described to compensate for this difference. The strain distribution in the solder was determined by finite element analysis and compared to the distribution generated by a theoretical 'ideal' test which generates an almost pure shear mode in the solder. By using a damage-based constitutive law the shape of the crack generated in the specimen has been predicted for both the actual test and the ideal pure shear test. Results from the simulations are also compared with experimental data using SnAgCu solder.

Finite element idealization of a cold-formed steel portal frame

A simple linear beam idealization of a cold-formed steel portal frame is presented in which beam elements are used to idealize the column and rafter members, and rotational spring elements are used to represent the rotational flexibility of the joints. In addition, the beam idealization takes into account the finite connection length of the joints. Deflections predicted using the beam idealization are shown to be comparable to deflections obtained from both a linear finite element shell idealization and full-scale laboratory tests. Using the beam idealization, deflections under rafter load are divided into three components: Deflection due to flexure of the column and rafter members, deflection due to bolt-hole elongation, and deflection due to in-plane bracket deformation. Of these deflection components, the deflection due to bolt-hole elongation is the most significant and cannot, therefore, be ignored. Using the beam idealization, engineers can analyze and design cold-formed steel portal frames, including making appropriate allowances for connection effects, without the need to resort to expensive finite element shell analysis.