Implicit TVD schemes applied to gas-droplet detonation calculation

ＩｎｔｒｏｄｕｃｔｉｏｎＴＶＤ(ＴｏｔａｌＶａｒｉａｔｉｏｎＤｉｍｉｎｉｓｈｉｎｇ)ｓｃｈｅｍｅｓｈａｖｅｂｅｅｎｗｉｄｅｌｙａｐｐｌｉｅｄｆｏｒｎｕｍｅｒｉｃａｌｃａｌｃｕｌａｔｉｏｎｏｆｆｌｕｉｄｄｙｎａｍｉｃｐｒｏｂｌｅｍｓｓｉｎｃｅＡｍｉＨａｒｔｅｎｄｅｖｅｌｏｐｅｄｎｅｗｈｉｇｈ＿ｒｅｓｏｌｕｔｉｏｎＴＶＤｆｉｎｉｔｅｄｉｆｆｅｒｅｎｃｅｓｃｈｅｍｅｓｉｎｔｈｅｅａｒｌｙｅｉｇｈｔｉｅｓ.Ｂｕｔｗｅｈａｖｅｎｏｔｆｏｕｎｄｔｈｅｒｅｐｏｒｔｓａｂｏｕｔｔｈｅａｐｐｌｉｃａｔｉｏｎｉｎ2Ｄｎ…     

Nonlinear Dynamics Modeling of Mechanical Periodicity of End Diastolic Volume of Left Ventricle

ＮｏｍｅｎｃｌａｔｕｒｅＸｎ 1 ＴｈｅＥＤＶｏｆｔｈｅ (ｎ 1 )ｔｈｈｅａｒｔｂｅａｔＸｎ 　ＴｈｅＥＤＶｏｆｔｈｅｎｔｈｈｅａｒｔｂｅａｔＹｎ 　ＴｈｅｅｊｅｃｔｖｏｌｕｍｅｏｆｔｈｅｎｔｈｈｅａｒｔｂｅａｔＵｎ 　ＴｈｅｍａｓｓｏｆｔｈｅｂｌｏｏｄｔｈｒｏｕｇｈｖｅｉｎｔｏｈｅａｒｔＰｓ 　ＴｈｅｓｙｓｔｏｌｉｃｐｒｅｓｓｕｒｅｏｆｌｅｆｔｖｅｎｔｒｉｃｌｅＰｄ 　ＴｈｅｄｉａｓｔｏｌｉｃｐｒｅｓｓｕｒｅｏｆｌｅｆｔｖｅｎｔｒｉｃｌｅＰｖ 　ＴｈｅｐｒｅｓｓｕｒｅＰａ 　Ｔｈｅａｒｔｅｒｉａｌｐｒｅｓ…     

The conservation law and the radiation of general relativity

ＴＨＥＣＯＮＳＥＲＶＡＴＩＯＮＬＡＷＡＮＤＴＨＥＲＡＤＩＡＴＩＯＮＯＦＧＥＮＥＲＡＬＲＥＬＡＴＩＶＩＴＹＬａｍＹｕ－ｙｉｕ（林汝耀）（ＤｅｐａｒｔｍｅｎｔｏｆＡｐｐｌｉｅｄＭａｔｈｅｍａｔｉｃｓ）ＨｏｎｇＫｏｎｇＰｏｌｙｔｅｃｈｎｉｃ．ＨｏｎｇＫｏｎ...     

INVARIANT MANIFOLDS AND THEIR STABILITY IN A THREE－DIMENSIONAL MEASURE－PRESERVING MAPPING SYSTEMS

ＩＮＶＡＲＩＡＮＴＭＡＮＩＦＯＬＤＳＡＮＤＴＨＥＩＲＳＴＡＢＩＬＩＴＹＩＮＡＴＨＲＥＥ－ＤＩＭＥＮＳＩＯＮＡＬＭＥＡＳＵＲＥ－ＰＲＥＳＥＲＶＩＮＧＭＡＰＰＩＮＧＳＹＳＴＥＭＳＬｉｕｊｉｅ（刘杰）（ＲｅｃｅｉｖｅｄＤｅｃ．１．１９９４；Ｃｏｍｍｕｎｉｃ...     

Study on Exact Analytical Solutions for Two Systems of Nonlinear Evolution Equations

ＩｎｔｒｏｄｕｃｔｉｏｎＤｕｒｉｎｇｔｈｅｃｏｕｒｓｅｏｆｓｔｕｄｙｉｎｇｔｈｅｗａｔｅｒｗａｖｅ,ｍａｎｙｃｏｍｐｌｅｔｅｌｙｉｎｔｅｇｒａｂｌｅｍｏｄｅｌｓｗｅｒｅｏｂｔａｉｎｅｄ ,ｓｕｃｈａｓＫｄＶｅｑｕａｔｉｏｎ ,ｍＫｄＶｅｑｕａｔｉｏｎ ,(2 1 )＿ｄｉｍｅｎｓｉｏｎａｌＫＰｅｑｕａｔｉｏｎ ,ｃｏｕｐｌｅｄＫｄＶｅｑｕａｔｉｏｎｓ,ｖａｒｉａｎｔＢｏｕｓｓｉｎｅｓｑｅｑｕａｔｉｏｎｓ ,ＷＫＢｅｑｕａｔｉｏｎｓｅｔｃ .[1- 13 ].Ｉｎｏｒｄｅｒｔｏｆｉｎｄｅｘｐｌｉｔｉｃｅｘａｃｔｓｏｌｕｔｉｏ…     

Iterative process for lipschitz local strictly pseudocontractive mappings

ＩＴＥＲＡＴＩＶＥＰＲＯＣＥＳＳＦＯＲＬＩＰＳＣＨＩＴＺＬＯＣＡＬＳＴＲＩＣＴＬＹＰＳＥＵＤＯＣＯＮＴＲＡＣＴＩＶＥＭＡＰＰＩＮＧＳＤｅｎｇＬｅｉ（邓磊）（ＣｈｏｎｇｑｉｎｇＴｅａｃｈｅｒ＇ｓＣｏｌｌｅｇｅ,Ｙｏｎｇｃｈｕａｎ,Ｃｈｏｎｇｑｉｎｇ）Ｄ...     

INVARIANT MANIFOLDS AND THEIR STABILITY IN A THREE－DIMENSIONAL MEASURE－PRESERVING MAPPING SYSTEMS

ＩＮＶＡＲＩＡＮＴＭＡＮＩＦＯＬＤＳＡＮＤＴＨＥＩＲＳＴＡＢＩＬＩＴＹＩＮＡＴＨＲＥＥ－ＤＩＭＥＮＳＩＯＮＡＬＭＥＡＳＵＲＥ－ＰＲＥＳＥＲＶＩＮＧＭＡＰＰＩＮＧＳＹＳＴＥＭＳＬｉｕＪｉｅ（刘杰）（ＲｅｃｅｉｖｅｄＤｅｃ．１．１９９４;Ｃｏｍｍｕｎｉｃ...     

[0,k_i]_1~m-FACTORIZATIONS ORTHOGONAL TO A SUBGRAPH

1　ＢａｓｉｃＤｅｆｉｎｉｔｉｏｎｓａｎｄＮｏｔａｔｉｏｎｓＷｅｏｎｌｙｄｅａｌｗｉｔｈｆｉｎｉｔｅｕｎｄｉｒｅｃｔｅｄｓｉｍｐｌｅｇｒａｐｈｓ.ＬｅｔＧｂｅａｇｒａｐｈｗｉｔｈｖｅｒｔｅｘｓｅｔＶ(Ｇ)ａｎｄｅｄｇｅｓｅｔＥ(Ｇ) .Ｆｏｒｅａｃｈｘ∈Ｖ(Ｇ) ,ｗｅｄｅｎｏｔｅｔｈｅｄｅｇｒｅｅｏｆｘｉｎＧｂｙｄＧ(ｘ) .Ｌｅｔｇ ,ｆｂｅｔｗｏｉｎｔｅｇｅｒ＿ｖａｌｕｅｄｆｕｎｃｔｉｏｎｓｄｅｆｉｎｅｄｏｎＶ(Ｇ)ｓｕｃｈｔｈａｔｇ(ｘ)≤ｆ(ｘ)ｆｏｒａｌｌｘ∈Ｖ(Ｇ) .Ａ (ｇ ,ｆ)＿ｆａｃｔｏｒｏｆＧｉｓａｓ…     

BASIC EQUATIONS, THEORY AND PRINCIPLES OF COMPUTATIONAL STOCK MARKET (Ⅱ)——BASIC PRINCIPLES

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｒｅａｒｅｍａｎｙｐｒｉｎｃｉｐｌｅｓｉｎｍｅｃｈａｎｉｃｓ,ｓｕｃｈａｓｔｈｅＳａｉｎｔ＿Ｖｅｎａｎｔ’ｓｐｒｉｎｃｉｐｌｅ,ｔｈｅＤ’Ａｌｅｍｂｅｒｔ’ｓｐｒｉｎｃｉｐｌｅ,ｖａｒｉａｔｉｏｎａｌｐｒｉｎｃｉｐｌｅｓ...     

A THREE-FLUID MODEL OF THE SAND-DRIVEN FLOW

ＡＴＨＲＥＥ－ＦＬＵＩＤＭＯＤＥＬＯＦＴＨＥＳＡＮＤ－ＤＲＩＶＥＮＦＬＯＷ￥(刘大有，董飞)ＬｉｕＤａｙｏｕ；ＤｏｎｇＦｅｉ（ＩｎｓｔｉｔｕｔｅｏｆＭｅｃｈａｎｉｃｓ，ＡｃａｄｅｍｉａＳｉｎｉｃａ，Ｂｅｉｊｉｎｇ１０００８０，Ｐ．Ｒ．Ｃｈｉｎａ）Ａｂｓ...     

一种空间缆索结构静力分析的解析元法

将空间缆索结构简化为具有拉伸刚度的质点系,给出了缆索结构空间解析元法的基本方程和求解方法,单元间的作用力与坐标变化的关系可以用解析法得到,对所得到的反映结构特性的质点系方程组进行力的平衡迭代,求解方程组.采用自动的动态可变步长的迭代方法,能够提高计算效率,保证收敛.这种方法既考虑了几何非线性,又适用于材料非线性的计算,比有限元法优越之处还在于,它不用求解线性方程组,所以适用范围广,允许求解多自由度的几何可变体系,而有限元法在求解此类问题时经常不收敛.     

Constrained parameter-splitting perturbation method for the improved solutions to the nonlinear vibrations of Euler–Bernoulli cantilevers

In this paper, a new method named constrained parameter-splitting perturbation method for improving the solutions obtained from the parameter-splitting perturbation method is proposed for solving the problems in some extremal cases, such as the strongly nonlinear vibration of an Euler–Bernoulli cantilever. The proposed method takes the advantages of both the perturbation method and the harmonic balance method. The idea is that the solution obtained by the parameter-splitting perturbation method is substituted into the equation of motion and then the accumulative error of the equation is minimized for determining the unknown splitting parameters under the constraints constructed under the frame of harmonic balance method. The forced vibration of an oscillator with cubic geometric nonlinearity and inertia nonlinearity and the forced vibration of a planar microcantilever beam with a lumped tip mass are studied as examples to reveal the efficacy of the proposed method. The inspection of the steady-state response including its stability is conducted by means of comparing the frequency-response curves obtained by the proposed method with those obtained by the numerical continuation method and harmonic balance method, respectively, to show the efficacy and the advantages of the proposed method. Meanwhile, the nonlinear ordering effect on the solutions of the proposed method is also studied by comparing the results obtained by using different nonlinear orderings in the systems. In the last, we found through convergence examinations that it is necessary to have corrections to the erroneous solution which are obtained by harmonic balance method and Floquet theory in stability analysis.

     

基于剩余推力法的地震滑坡永久位移研究

基于剩余推力法思想 ,结合Newmark有限滑动位移法 ,考虑了由于动力作用造成的孔隙水压力变化 ,对一种最为常见的边坡灾害—滑坡 ,提出了一种简便的估算地震动力永久位移的方法。对一实例用该法和快速拉格朗日元 (FLAC3D)进行对比计算 ,结果表明两者的结果基本接近 ,前者要保守一些。这就使得应用剩余推力法这一常规方法对滑坡进行真正意义上的动力时程分析成为可能.     

一种从离散模拟到连续介质弹性模拟的过渡方法

提出了一种从离散分子动力学模拟(MD)到连续介质弹性有限元计算分析(FEA)的过渡方法, 简称MD-FEA方法. 首先通过MD计算获得晶体材料原子的移动位置, 然后根据晶体结构的周期性特征构造连续介质假设下的有限单元变形模型, 进一步结合材料的力学行为本构关系获得应变和应力场. 为了检验MD-FEA方法的有效性, 将该方法应用于详细分析Al-Ni软硬组合两相材料纳米柱体的拉伸变形问题和基底材料为Al球形压头材料为金刚石的纳米压痕问题. 采用MD-FEA方法获得了上述两种问题的应力?应变场, 并将计算结果分别与传统MD方法中通过变形梯度计算的原子应变以及原子的位力应力进行了比较, 详细讨论了用MD-FEA方法计算的应力?应变场与传统MD原子应变和位力应力的区别, 并对MD-FEA方法的有效性及其相较于传统MD方法所具有的优势进行了探讨. 结论显示, MD-FEA方法与传统MD方法在应力?应变变化平缓的区域得到的结果接近, 但在变化剧烈的区域以及材料的表/界面区域, MD-FEA方法能够得到更加精确的结果. 同时, MD-FEA方法避免了传统MD方法中, 需要人为选取截断半径以及加权函数所导致的误差. 另外, 当应变较大时, MD-FEA方法计算的小应变与传统MD方法计算的格林应变存在一定差异, 因此, MD-FEA方法更适合应变较小的情形.      

不可压黏性流问题的无网格法分析

不可压缩黏性流问题一般采用Navier-Stokes方程来描述，基于加权残值法，推导了问题的无网格伽辽金法（EFGM）离散Navier-Stokes方程，在时间域上采用分步方法计算，速度和压力由相互独立的方程以解耦的形式求解，并采用同阶移动最小二乘近似，在每一时间步中，对压力解和速度解采用了Newton-Raphson迭代法进行修正，最后将所得到的方法应用到剪切驱动空腔流问题中，验证了方法的有效性，且解的精度高、稳定性好。     

样条配点法分析结构动力响应的无条件稳定计算格式

本文以三次B样条函数作为位移时域函数,应用配点法,在任一时间步长中列出两个时刻的结构物运动方程式残值为零的条件并使结果满足了稳定条件得到一个计算结构物动力响应无条件稳定的计算格式。 这个计算格式于参数θ符合0.15<θ<0.05时是绝对无条件稳定的,研究证明,本文所提出的格式的精度根据三项指标:振幅衰减率 AD,周期伸长率 PE 及算法阻尼比都比较Wilson-θ法,Newmark法及Houbolt法的为优,例题计算也证明此点,即本文方法精度较上述三法为高,本文格式简单,准确,工作量不多,可以作为计算结构动力响应有效方法之一。     

捷联惯导系统最简多位置解析对准改进算法

最简多位置解析对准作为测量IMU偏置的简单方法,降低了传统多位置解析对准的工程复杂度,摆脱了伺服平台的限制,但是最简多位置解析对准目前使用的是几何方法,几何方法计算复杂,精度受先验信息的影响。针对这一问题,提出了最简多位置对准解析算法,通过解析算法可得到IMU偏置的解析解。通常两位置即可求得IMU偏置,特殊情况需要三位置,某些特殊位置可以计算特定轴向的偏置。解析法计算简单,不受先验信息的影响,并通过实例仿真证实该算法的有效性,计算误差主要由一阶近似误差和测量噪声引起。     

一种基于半参数回归的加速度计误差模型辨识方法

为了解决加速度计离心机试验中系统误差和未建模误差对加速度计模型辨识的影响,将统计学中的半参数回归方法引入到加速度计的模型辨识中,建立了加速度计的半参数回归模型,提出了一种基于最小二乘-半参数回归模型(LS-SPRM)的估计方法,该估计方法利用最小二乘法估计加速度计的误差模型系数,利用半参数回归方法估计加速度计测试中的系统误差,并通过检验残差是否为白噪声作为判断是否有系统误差的条件。在半参数模型的估计中,采用二阶段估计方法,利用三次样条函数进行非参数部分的估计,并讨论了光滑参数的选取方法。仿真试验结果表明,采用该方法能够较好地补偿由于系统误差和未建模误差带来的影响,使加速度计模型辨识的标准差较普通最小二乘法减小45%左右,估计的残差也减小了近一倍。     

大跨屋盖的完全三分量等效静力风荷载Refined three components method for determining equivalent static wind loads of large span roof

惯性风荷载法物理意义不明确,传统三分量叠加法忽略背景‐共振耦合效应,本文提出一种完全三分量的等效静力风荷载方法。利用完全三分量方法研究了开孔屋盖的等效静力风荷载,研究结果表明,阵风荷载因子法除在等效节点处有较好的精确性外,其他非等效节点处的等效静力风荷载均偏离正常分布,不适宜求解大跨屋盖等效静力风荷载;完全三分量法具有明显物理意义,其背景部分可采用全模态的L RC法实现,比截断模态的惯性风荷载法更加精确;忽略背景‐共振耦合效应的传统三分量法低估等效静力风荷载10％左右,完全三分量法可考虑背景‐共振耦合效应,相应于传统三分量法更加精确。     

扩散方程单内点精细积分法与差分法比较研究

一维扩散方程初值问题可以用全域或子域精细积分求解。子域积分可以采用不同数量的内点，单内点是其最简单的情况。当单内点精细积分中的传递函数即指数函数用其泰勒展开式的一阶近似来替代时，精细积分转化为差分方程。本文研究了这一对应关系。各种常见差分格式均找到了对应的单点精细积分格式，并在单点精细积分一般公式中得到了统一表达形式