A multiscale work-analysis approach for geotechnical structures

This paper presents a second-order work analysis in application to geotechnical problems by using a novel effective multiscale approach. To abandon complicated equations involved in conventional phenomenological models, this multiscale approach employs a micromechanically-based formulation, in which only four parameters are involved. The multiscale approach makes it possible a coupling of the finite element method (FEM) and the micromechanically-based model. The FEM is used to solve the boundary value problem (BVP) while the micromechanically-based model is utilized at the Gauss point of the FEM. Then, the multiscale approach is used to simulate a three-dimensional triaxial test and a plain-strain footing. On the basis of the simulations, material instabilities are analyzed at both mesoscale and global scale. The second-order work criterion is then used to analyze the numerical results. It opens a road to interpret and understand the micromechanisms hiding behind the occurrence of failure in geotechnical issues.     

Plane of elastic non-resisting tension material under foundation structures

In the frame of 2D-static problems one approaches the problem of elastic-NRT (not-resisting tension) semi-plane loaded on its limit line. This problem is intended to model the stress situation induced in the soil by a foundation structure. The solution, in terms of activated stress field, is searched for in the class of stress fields satisfying equilibrium and admissibility conditions, by applying an energy approach. Copyright © 2004 John Wiley & Sons, Ltd.     

显性语法知识传授中的发现原则

尽管语法在交际教学法鼎盛时期一度遭到冷遇，但现在又重新被英语语言教学所接纳，作为一种新的语法教学方式，发现原则综合了传统归纳发和演绎法的长处，概述了发现原则的理论依据，指出发现原则旨在鼓励学习者积极参与而非被动接受，符合当前认知主义对二语学习的认识。     

无规取向、无吸收非球形粒子光散射的T-矩阵方法

较详细地研究了无规取向、无吸收椭球体粒子的T－矩阵收敛问题。首先，简要概括了nmax的3种收敛方案和它们的基本特性。然后，应用1993年提出的数学收敛方法（M－方法）和1998年提出和物理收敛方法（P－方法）研究收敛问题。结果表明椭球粒子收敛精度对粒子的尺度参数，纵横比以及椭球体的种类（例如，长／扁椭球）有很强的依赖性。当粒子的尺度参数不太大时，甚至在极端纵横比的条件下，P－收敛方案优于M－收敛方案。     

关于非球形粒子光散射的T-矩阵数值计算方法的研究

首先介绍了无规取向轴对称非球形粒子光散射的T－矩阵数值方法。然后，进一步研究了T－矩阵计算的收敛问题，提出了一种新的关于无规取向无吸收非线形散射物体的物理收敛方法，我们研究了椭球粒子的收敛问题，结果表明我们的数值方法和T－矩阵计算程序是有效的， 收敛精度与粒子的尺度和形状有很强的依赖性，在一定的条件下，我们的物理收敛速度优于NASA Mishchenko的数学收敛速序。     

△T球体异常的三元解释法

在△Ｔ异常解释课题中，一个新的角度参量─—视磁化特征角的拟定，最大限度地简化了球体磁场表达式。从而才有可能通过特殊解法，给出了球体磁场的极值与拐点公式。作者据此建立了四种不受正常场制约、仅需异常的局部持征曲线即可全面解释三度异常的解释方法。屏障分解、叠加判别、三元迭代的设计，为解释三度体叠加异常给出了简单而易行的途径。     

Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

Tidal inundation mapping under enhanced land subsidence in Semarang,Central Java Indonesia

Tidal inundation by high tide under enhanced land subsidence is a damaging phenomenon and a major threat to the Semarang urban area in Indonesia. It impacts on economic activities, as well as the cost of an emergency program and causes interruption of pubic services, danger of infectious diseases and injury to human lives. This study examines a spatial analysis tool on the GIS-raster system for the tidal inundation mapping based on the subsidence-benchmark data and modified detail digital elevation model. Neighborhood operation and iteration model as a spatial analysis tool have been applied in order to calculate the encroachment of the tidal inundation on the coastal area. The resulting map shows that the tidal flood spreads to the lowland area and causes the inundation of coastal settlement, infrastructure, as well as productive agricultural land, i.e., the fish-pond area. The monitoring of the vulnerable area due to the tidal inundation under the scenario of extended land subsidence plays an important role in long-term coastal zone management in Semarang.     

Internal generation of waves: Delta source function method and source term addition method

In this study, we investigate two internal wave generation methods in numerical modeling of time-dependent equations for water wave propagation, i.e., delta source function method and source term addition method, the latter of which has been called the line source method in literatures. We derive delta source functions for the Boussinesq-type equations and extended mild-slope equations. By applying the fractional step splitting method, we show that the delta source function method is equivalent to the source term addition method employing the energy velocity. This suggests that the energy velocity should be used rather than the phase velocity for the transport of incident wave energy in the source term addition method. Finally, the performance of the delta source function method is verified by accurately generating nonlinear cnoidal waves as well as linear waves for horizontally one-dimensional cases.     

Freely floating-body simulation by a 2D fully nonlinear numerical wave tank

Nonlinear interactions between large waves and freely floating bodies are investigated by a 2D fully nonlinear numerical wave tank (NWT). The fully nonlinear 2D NWT is developed based on the potential theory, MEL/material-node time-marching approach, and boundary element method (BEM). A robust and stable 4th-order Runge–Kutta fully updated time-integration scheme is used with regriding (every time step) and smoothing (every five steps). A special φ_n-η type numerical beach on the free surface is developed to minimize wave reflection from end-wall and wave maker. The acceleration-potential formulation and direct mode-decomposition method are used for calculating the time derivative of velocity potential. The indirect mode-decomposition method is also independently developed for cross-checking. The present fully nonlinear simulations for a 2D freely floating barge are compared with the corresponding linear results, Nojiri and Murayama’s (Trans. West-Jpn. Soc. Nav. Archit. 51 (1975)) experimental results, and Tanizawa and Minami’s (Abstract for the 6th Symposium on Nonlinear and Free-surface Flow, 1998) fully nonlinear simulation results. It is shown that the fully nonlinear results converge to the corresponding linear results as incident wave heights decrease. A noticeable discrepancy between linear and fully nonlinear simulations is observed near the resonance area, where the second and third harmonic sway forces are even bigger than the first harmonic component causing highly nonlinear features in sway time series. The surprisingly large second harmonic heave forces in short waves are also successfully reproduced. The fully updated time-marching scheme is found to be much more robust than the frozen-coefficient method in fully nonlinear simulations with floating bodies. To compare the role of free-surface and body-surface nonlinearities, the body-nonlinear-only case with linearized free-surface condition was separately developed and simulated.