Computational modelling of building damage due to tunnelling

Underground constructions in soft ground may lead to settlement damage to existing buildings. In The Netherlands the situation is particularly complex, because of the combination of soft soil, fragile pile foundations and brittle, unreinforced masonry façades. The tunnelling design process in urban areas requires a reliable risk damage assessment. In the engineering practice the current preliminary damage assessment is based on the limiting tensile strain method (LTSM). Essentially this is an uncoupled analysis, in which the building is modelled as an elastic beam subject to imposed Greenfield settlements and the induced tensile strains are compared with a limit value for the material. The soil-structure interaction is included only as a ratio between the soil and the building stiffness. In this paper, a coupled approach is evaluated. The soil-structure interaction in terms of normal and shear behaviour is represented by interface elements and a cracking model for masonry is included. This project aims to improve the existing damage classification system for masonry buildings subjected to tunnel-induced settlement, in order to evaluate the necessity of strengthening techniques or mitigation measures.

Linearization Method and Linear Complexity

We focus on the relationship between the linearization method and linear complexity and show that the linearization method is another effective technique for calculating linear complexity. We analyze its effectiveness by comparing with the logic circuit method. We compare the relevant conditions and necessary computational cost with those of the Berlekamp-Massey algorithm and the Games-Chan algorithm. The significant property of a linearization method is that it needs no output sequence from a pseudo-random number generator (PRNG) because it calculates linear complexity using the algebraic expression of its algorithm. When a PRNG has n [bit] stages (registers or internal states), the necessary computational cost is smaller than O(2n). On the other hand, the Berlekamp-Massey algorithm needs O(N2) where N ( 2n) denotes period. Since existing methods calculate using the output sequence, an initial value of PRNG influences a resultant value of linear complexity. Therefore, a linear complexity is generally given as an estimate value. On the other hand, a linearization method calculates from an algorithm of PRNG, it can determine the lower bound of linear complexity.

encoding-based method for pattern matching

A new algorithm for compiling pattern matching is presented. Different from the traditional traversal-based approaches, it can represent a sequence of patterns as an integer by an encoding method and translate equations into case-expressions. The algorithm is simple to implement, and efficient for a kind of patterns, i.e. simple and dense patterns. This method can be complementary to traditional approaches.

U-transformation-finite element method in 3-dimensional analysis of orthotropic laminates

For an orthotropic laminate, an equivalent system with doubly cyclic periodicity is introduced. Then a 3-dimensional finite element model for the equivalent system is transformed into the unitary space, where the large finite element matrix equation is decoupled into some small matrix equations. Such a decoupling very efficiently reduces the computational effort. For an orthotropic laminate with four clamped edges, no exact elasticity solution is available, and the deflection values predicted by different methods have a considerable difference each other for a small length-to-thickness ratio. The present predictions are the largest because the present method is a full 3-dimensional finite element analysis without superfluous constraints. Illustrative numerical examples are presented to observe the distributions of stresses through the thickness of the laminates. (C) 2010 Elsevier Ltd. All rights reserved.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been developed to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

A global shallow water model using high order multi-moment constrained finite volume method and icosahedral grid

A novel accurate numerical model for shallow water equations on sphere have been developed by implementing the high order multi-moment constrained finite volume (MCV) method on the icosahedral geodesic grid. High order reconstructions are conducted cell-wisely by making use of the point values as the unknowns distributed within each triangular cell element. The time evolution equations to update the unknowns are derived from a set of constrained conditions for two types of moments, i.e. the point values on the cell boundary edges and the cell-integrated average. The numerical conservation is rigorously guaranteed. in the present model, all unknowns or computational variables are point values and no numerical quadrature is involved, which particularly benefits the computational accuracy and efficiency in handling the spherical geometry, such as coordinate transformation and curved surface. Numerical formulations of third and fourth order accuracy are presented in detail. The proposed numerical model has been validated by widely used benchmark tests and competitive results are obtained. The present numerical framework provides a promising and practical base for further development of atmospheric and oceanic general circulation models. (C) 2009 Elsevier Inc. All rights reserved.

A New High Order Accurate Shock Capture Method with Wave Booster

A new kind of shock capturing method is developed. Before applying the high order accurate traditional scheme which is called as base scheme in this paper the fluid parameters are preconditioned in order to control the group velocity. The newly constructed scheme is high order accurate, simple, has high resolution of the shock, and less computer time consumed.

Computational study on the addition of HCN to methanimine catalyzed by formamidine and formamide

In this paper, we have carried out a theoretical study on the addition of HCN to methanimine with formamidine or formamide using second-order Moller-Plesset perturbation (MP2) method with 6-31 + G(d,p) basis sets. At MP2 level. a high-energy, intermediate has been located for each pathway. The addition of HCN to methanimine with formamidine has the lowest free energy barrier according to the calculations at MP2 level. (C) 2004 Elsevier B.V. All rights reserved.

a requirement traceability refinement method based on relevance feedback

Knowledge Systems Institute Graduate School

Investigation of Chemical Bond Characteristics, Thermal Expansion Coefficients and Bulk Moduli of alpha-R2MoO6 and R2Mo2O7 (R = rare earths) by Using a Dielectric Chemical Bond Method

Theoretical researches are performed on the alpha-R2MoO6 (R = Y, Gd, Tb Dy, Ho, Er, Tm and Yb) and pyrochlore-type R2Mo2O7 (R = Y, Nd, Sm, Gd, Tb and Dy) rare earth molybdates by using chemical bond theory of dielectric description. The chemical bonding characteristics and their relationship with thermal expansion property and compressibility are explored. The calculated values of linear thermal expansion coefficient (LTEC) and bulk modulus agree well with the available experimental values. The calculations reveal that the LTECs and the bulk moduli do have linear relationship with the ionic radii of the lanthanides: the LTEC decreases from 6.80 to 6.62 10(-6)/K and the bulk modulus increases from 141 to 154 GPa when R goes in the order Gd, Tb Dy, Ho, Er, Tm, and Yb in the alpha-R2MoO6 series; while in the R2Mo2O7 series, the LTEC ranges from 6.80 to 6.61 10(-6)/K and the bulk modulus ranges from 147 to 163 GPa when R varies in the order Nd, Sm, Gd, Tb and Dy.

Impact of difference accuracy on computational properties of vertical grids for a nonhydrostatic model

Starting from nonhydrostatic Boussinesq approximation equations, a general method is introduced to deduce the dispersion relationships. A comparative investigation is performed on inertia-gravity wave with horizontal lengths of 100, 10 and 1 km. These are examined using the second-order central difference scheme and the fourth-order compact difference scheme on vertical grids that are currently available from the perspectives of frequency, horizontal and vertical component of group velocity. These findings are compared to analytical solutions. The obtained results suggest that whether for the second-order central difference scheme or for the fourth-order compact difference scheme, Charny-Phillips and Lorenz ( L) grids are suitable for studying waves at the above-mentioned horizontal scales; the Lorenz time-staggered and Charny-Phillips time staggered (CPTS) grids are applicable only to the horizontal scales of less than 10 km, and N grid ( unstaggered grid) is unsuitable for simulating waves at any horizontal scale. Furthermore, by using fourth-order compact difference scheme with higher difference precision, the errors of frequency and group velocity in horizontal and vertical directions produced on all vertical grids in describing the waves with horizontal lengths of 1, 10 and 100 km cannot inevitably be decreased. So in developing a numerical model, the higher-order finite difference scheme, like fourth-order compact difference scheme, should be avoided as much as possible, typically on L and CPTS grids, since it will not only take many efforts to design program but also make the calculated group velocity in horizontal and vertical directions even worse in accuracy.