基于双尺度渐近分析的有限元算法

1．引言正如文山所说，由于复合材料和周期结构的材料系数ail（x）在局部区域内间断且跳跃性很大，加上区域内含有周期性洞穴或裂缝，且周期长度很小．一般而言，直接采用有限元方法进行数值模拟，其计算量大得惊人，甚至难以实现．文山针对这种特征，提出了一种可计算的双尺度渐近分析模式，本文在此基础上给出了相应的有限元算法，它包括：1．周期解在一个基本构造上的有限元计算；2．边界层的有限元计算．同时，给出了相应的误差分析．2．周期解的有限元计算首先考虑下列形式的边值问题；其中把，代E尸（on叫，iii（0关于E—（EI，ZZ…     

Convex central configurations for the n-body problem

We give a simple proof of a classical result of MacMillan and Bartky (Trans. Amer. Math. Soc. 34 (1932) 838) which states that, for any four positive masses and any assigned order, there is a convex planar central configuration. Moreover, we show that the central configurations we find correspond to local minima of the potential function with fixed moment of inertia. This allows us to show that there are at least six local minimum central configurations for the planar four-body problem. We also show that for any assigned order of five masses, there is at least one convex spatial central configuration of local minimum type. Our method also applies to some other cases.     

Three New Chromone Derivatives Produced by Phomopsis sp. HNY29‐2B from Acanthus ilicifolius Linn.

Three new chromone derivatives, phomochromenones A‐C ( 1–3 ), and one known chaetocyclinone B ( 4 ) were obtained from the cultures of Phomopsis sp. HNY29‐2B isolated from the mangrove Acanthus ilicifolius Linn., which was collected from the South China Sea. Their structures were determined by the analysis of ¹D NMR and ²D NMR as well as mass spectroscopic data. The absolute configurations of 1 and 2 were assigned by quantum chemical calculations of the electronic circular dichroism (ECD ) spectra. Compound 3 is the third example of alkaloids possessing the unique chromeno [3,2‐c ] pyridine nucleus. In the bioactivity assay, compound 4 showed cytotoxicity against human prostate cancer cell lines (PC ‐3 and DU145 cells) with the IC₅₀ values of 8.13 and 3.59 µmol/L, respectively.     

A DQ based approach to simulate the vibrations of buckled beams

In this paper, the nonlinear planar response of a hinged–hinged buckled beam to a primary-resonance excitation of its first vibration mode is computed by a new numerical scheme. The beam is subjected to an axial force beyond the critical load of the first buckling mode and to a transverse harmonic excitation. The nonlinear dynamical problem is solved by deducing directly the discretized equations governing the problem thanks to a new approach, here called DQ based approach, since it is based on the application of the quadrature rules of the DQM. As it will be shown, for the problem here considered, the minimum number of degrees of freedom to be retained to limit the numerical errors is four. Computer simulations of the dynamic behaviour of the discretized system are conducted by means of the IDQ method, a method proposed and recently generalized by the author. A sequence of supercritical period-doubling bifurcations leading to chaos, snapthrough motions and quasi-periodic motions can be observed, similarly to some cases existing in literature.     

On the numerical solution of tracked vehicle dynamic equations

In this investigation, the solution of the nonlinear dynamic equations of the multibody tracked vehicle systems are obtained using different procedures. In the first technique, which is based on the augmented formulation that employes the absolute Cartesian coordinates and Lagrange multipliers, the generalized coordinate partitioning of the constraint Jacobian matrix is used to determine the independent coordinates and the associated independent differential equations. An iterative Newton-Raphson algorithm is used to solve the nonlinear constraint equations for the dependent variables. The numerical problems encountered when one set of independent coordinates is used during the simulation of large scale tracked vehicle systems are demonstrated and their relationship to the track dynamics is discussed. The second approach employed in this investigation is the velocity transformation technique. One of the versions of this technique is discussed in this paper and the numerical problems that arise from the use of inconsistent system of kinematic equations are reported. In the velocity transformation technique, the tracked vehicle system is assumed to consist of two kinematically decoupled subsystems; the first subsystem consists of the chassis, the rollers, the sprocket and the idler, while the second subsystem consists of the track which is represented as a closed kinematic chain that consists of rigid links connected by revolute joints. It is demonstrated that the use of one set of recursive equations leads to numerical difficulties because of the change in the track configuration. Singular configurations can be avoided by repeated changes in the recursive equations. The sensitivity of the predictor-corrector multistep numerical integration schemes to the method of formulating the state equations is demonstrated. The numerical results presented in this investigation are obtained using a planner tracked vehicle model that consists of fifty four rigid bodies.     

纵横载荷作用下功能梯度梁的弯曲和过屈曲

基于一阶非线性梁理论和物理中面概念,导出了纵横向载荷作用下功能梯度材料(FGM)梁非线性弯曲和过屈曲问题的控制方程,并获得了该问题的精确解;据此解研究了梯度材料性质、外载荷、横向剪切变形以及边界条件等因素对功能梯度材料梁非线性力学行为的影响,分析中假设功能梯度材料性质只沿梁厚度方向,并按成分含量的幂指数函数形式变化。结果表明:纵横载荷共同作用下,功能梯度梁的弯曲构形将有无限多个;随着梯度指数的增大,梁的变形减小,临界载荷升高;随着长高比的增大,横向剪切变形的影响减小。     

Alternating sign matrices with one −1 under vertical reflection

We define a bijection that transforms an alternating sign matrix A with one −1 into a pair (N,E) where N is a (so called) neutral alternating sign matrix (with one −1) and E is an integer. The bijection preserves the classical parameters of Mills, Robbins and Rumsey as well as three new parameters (including E). It translates vertical reflection of A into vertical reflection of N. A hidden symmetry allows the interchange of E with one of the remaining two new parameters. A second bijection transforms (N,E) into a configuration of lattice paths called “mixed configuration.”     

正四面体套中心构型的存在唯一性

Two cases of the nested configurations in R^3 consisting of two regular quadrilaterals are discussed. One case of them do not form central configuration, the other case can be central configuration. In the second case the existence and uniqueness of the central configuration are studied. If the configuration is a central configuration, then all masses of outside layer are equivalent, similar to the masses of inside layer. At the same time the following relation between r(the ratio of the sizes) and mass ratio b = m/m must be satisfied b=24(3的立方根)（3r^2 2r 3）^-3/2-8（1-r）|1-r|^-3-3(6r的立方根)/24(3的立方根)（3 r）（3r^2 2r 3）^-3/2-8r（1-r）|1-r|^-3-3(6r^-2的立方根)in which the masses at outside layer are not less than the masses at inside layer, and the solution of this kind of central configuration is unique for the given ratio (b) of masses.     

Constitutive modeling of the mechanics associated with crystallizable shape memory polymers

Shape memory polymers are novel materials that can be easily formed into complex shapes, retaining memory of their original shape even after undergoing large deformations. The temporary shape is stable and return to the original shape is triggered by a suitable mechanism such as heating. In this paper, we develop constitutive equations to model the mechanical behavior of crystallizable shape memory polymers. Crystallizable shape memory polymers are called crystallizable because the temporary shape is fixed by a crystalline phase, while return to the original shape is due to the melting of this crystalline phase. The modeling is done using a framework that was developed recently for studying crystallization in polymers ([28], [25], [27], [31]) and is based on the theory of multiple natural configurations. In this paper we formulate constitutive equations for the original amorphous phase and the semi-crystalline phase that is formed after the onset of crystallization. In addition we model the melting of the crystalline phase to capture the return of the polymer to its original shape. The model has been used to simulate a typical uni-axial cycle of deformation, the results of this simulation compare very well with experimental data. In addition to this we also simulate circular shear of a hollow cylinder and present results for different cases in this geometry. Received: January 5, 2005     

Charge-conserving electron density averaging for a set of nuclear configurations

The Löwdin-Inverse Löwdin Transformation (LILT) technique combined with the Diophantine Density Matrix Purification (DDMP) method is suggested for a charge-preserving density averaging approach for a range of nuclear configurations.