有限元网格自动生成及自动调整技术

本文给出了一种任意平面域内三角形网格自动生成及自动调整算法，并编制了相应的程序。通过对带Ｕ形槽的三点变曲试件的分析，证明本文算法是可行的，并且具有稳定性能好、收敛速度快的优点。     

A two-scale second-moment turbulence closure based on weighted spectrum integration

A two-scale second-moment turbulence closure has been derived based on the weighted integration of the dynamic equation for the covariance spectrum. The goal is to close the Reynolds stress equations with two additional scalar equations that provide separately the scales of the spectral energy transfer and of the turbulence energy dissipation rate. Such a model should provide better prediction of nonequilibrium turbulent flows. The derivation consists of analytical integration of the wave-number-weighted covariance spectrum using a model of the spectral equations with an assumed simple representation of the shape of the energy spectrum. The resulting closure consists of a set of three tensorial equations, one for the Reynolds stress and two for length scale tensors, the latter representing the energy containing- and dissipative eddies respectively. The trace of the two tensor-scale equations leads to a set of two scalar scale parameters. In the equilibrium limit, the model reduces to the standard second-moment single-scale closure. The approach makes it also possible to derive the scale equations in a more systematic manner as compared with the common single-scale and other multi-scale models. The performance of the model in capturing the scale dynamics is illustrated by predictions of several generic homogeneous and inhomogeneous unsteady flows, demonstrating the expected response of the two scale equations. PACS 03.50.De, 04.20-q, 42.65-k     

The Secondary Bifurcation of an Aeroelastic Airfoil Motion: Effect of High Harmonics

The nonlinear dynamical response of a two-degree-of-freedom aeroelastic airfoil motion with cubic restoring forces is investigated. A secondary bifurcation after the primary Hopf (flutter) bifurcation is detected for a cubic hard spring in the pitch degree-of-freedom. Furthermore, there is a hysteresis in the secondary bifurcation: starting from different initial conditions the motion may jump from one limit cycle to another at different fluid flow velocities. A high-order harmonic balance method is employed to investigate the possible bifurcation branches. Furthermore, a numerical time simulation procedure is used to confirm the stable and unstable bifurcation branches.     

Second order effects in an elastic half-space acted upon by a non-uniform normal load

A closed form solution to the second order elasticity problem, when an isotropic compressible elastic half-space undergoes a deformation owing to a non-uniformly distributed normal load, is presented. The method of integral transform is employed and the case when loading is distributed, in accordance with Hertz's law, is discussed. The limiting solution for incompressible isotropic elastic material is also derived. Numerical calculations for the second order elastic material for the displacement and the normal stress in thez-direction are carried out. It is found that, in comparision to the linear elastic case, the displacement increases and the normal stress decreases in the second order elastic material

---

Sommario Si presenta una soluzione del secondo ordine in forma chiusa del problema di un semi-spazio elastico isotropo compressibile soggetto a un carico normale distribuito non uniformemente, ottenuta mediante il metodo delle trasformate integrali. Si discute il caso del carico distribuito secondo la legge di Hertz; si deriva anche la soluzione limite per materiale elastico isotropo incompressibile. Si ottengono numericamente alcuni risultati relativi a spostamenti e tensioni normali nella direzione z: si trova che, a paragone del caso elastico lineare, gli spostamenti aumentano e le tensioni diminuiscono nel materiale elastico del secondo ordine.

     

Resonance Vibrations and Dissipative Heating of Thin-Walled Laminated Elements Made of Physically Nonlinear Materials

The dynamic thermomechanical problem for thin-walled laminated elements is formulated based on the geometrically linear theory and Kirchhoff–Love hypotheses. A simplified model of vibrations and dissipative heating of structurally inhomogeneous inelastic bodies under harmonic loading is used. The mechanical properties of materials are described using strain-dependent complex moduli. A nonstationary vibration-heating problem is solved. The dissipative function, derived from the stationary solution, is used to specify internal heat sources. The amplitude–frequency characteristics and spatial distributions of the main field variables are studied for a sandwich beam subjected to forced vibrations     

Partial Uniqueness and Obstruction to Uniqueness in Inverse Problems for Anisotropic Elastic Media

We consider the inverse problem of identifying the density and elastic moduli for three-dimensional anisotropic elastic bodies, given displacement and traction measurements made at their surface. These surface measurements are modelled by the dynamic Dirichlet-to-Neumann map on a finite time interval. For linear or nonlinear anisotropic hyperelastic bodies we show that the displacement-to-traction surface measurements do not change when the density and elasticity tensor in the interior are transformed tensorially by a change of coordinates fixing the surface of the body to first order. Our main tool, a new approach in inverse problems for elastic media, is the representation of the equations of motion in a covariant form (following Marsden and Hughes, 1983) that preserves the underlying physics.In the case of classical linear elastodynamics we then investigate how the type of anisotropy changes under coordinate transformations. That is, we analyze the orbits of general linear, anisotropic elasticity tensors under the action by pull-back of diffeomorphisms that fix the surface of the elastic body to first order, and derive a pointwise characterization of parts of the orbits under this action. For example, we show that the orbit of isotropic elastic media, at any point in the body, consists of some transversely isotropic and some orthotropic elastic media. We then derive the first uniqueness result in the inverse problem for anisotropic media using surface displacement-traction data: uniqueness of three elastic moduli for tensors in the orbit of isotropic elasticity tensors. ^†Partially supported by an MSRI Postdoctoral Fellowship. Research at MSRI is supported in part by NSF grant DMS-9850361. This work was conducted while the first author was a Gibbs Instructor at Yale University. ^‡Partially supported by an MSRI Postdoctoral Fellowship, and by NSF grant DMS-9801664 (9996350).     

多孔饱和半空间上刚体垂直振动的轴对称混合边值问题

研究圆柱形刚体在多孔饱和半空间上的垂直振动．首先应用Ｈａｎｋｅｌ变换求解多孔饱和固体的动力基本方程———Ｂｉｏｔ波动方程．然后按混合边值条件建立多孔饱和半空间上刚体垂直振动的对偶积分方程，用Ａｂｅｌ变换化对偶积分方程为第二类Ｆｒｅｄｈｏｌｍ积分方程．文末给出了多孔饱和半空间表面动力柔度系数的计算曲线．     

TC4钛合金自然破片的引燃机理分析

对两种热处理条件不同的ＴＣ４钛合金圆筒在内部爆轰加载下的自然破片形成过程作了研究，并对破片的宏观断裂特征和细观断裂机理进行了分析。研究结果表明，一些引燃战斗部采用钛合金壳体的根本原因在于利用壳体破裂时形成的高温绝热剪切带和钛元素本身在高温下的剧烈氧化性。     

Selective Excitation with Asymmetric Adiabatic Pulses for NMR Spectroscopy

Existing selective pulses are mainly constructed in the forms of classically shaped pulses, such as the Gaussian pulses, or generated by using numerical optimization methods. However, all of these pulses are highly sensitive to radiofrequency (RF) intensity variation, which means their performance is highly dependent on the accuracy and stability of the RF intensity. Even a slight RF intensity deviation can cause severe degradation in the excitation profile. To solve this problem, we propose a method for narrow selective excitation by sequential application of a pair of phase‐opposite asymmetric adiabatic pulses, all within two scans. By retaining the adiabatic character, the new method is highly robust to RF intensity variation. Moreover, it has flexible excitation bandwidth, ranging from line‐selective to narrow‐band‐selective pulses. The method is tested both in numerical simulations and solution‐state NMR experiments.     

Multicomponent Molecular Orbital–Climbing Image–Nudged Elastic Band Method to Analyze Chemical Reactions Including Nuclear Quantum Effect

To analyze the H/D isotope effects on hydrogen transfer reactions in XHCHCHCHY?XCHCHCHYH (X, Y=O, NH, or CH₂) including the nuclear quantum effect of proton and deuteron, we propose a multicomponent molecular orbital‐climbing image‐nudged elastic band (MC_MO–CI–NEB) method. We obtain not only transition state structures but also minimum‐energy paths (MEPs) on the MC_MO effective potential energy surface by using MC_MO–CI–NEB method. We find that nuclear quantum effect affects not only stationary‐point geometries but also MEPs and electronic structures in the reactions. We clearly demonstrate the importance of including nuclear quantum effects for H/D isotope effect on rate constants (k_H/k_D).