Implicit and explicit integration in the solution of the absolute nodal coordinate differential/algebraic equations

This investigation is concerned with the use of an implicit integration method with adjustable numerical damping properties in the simulation of flexible multibody systems. The flexible bodies in the system are modeled using the finite element absolute nodal coordinate formulation (ANCF), which can be used in the simulation of large deformations and rotations of flexible bodies. This formulation, when used with the general continuum mechanics theory, leads to displacement modes, such as Poisson modes, that couple the cross section deformations, and bending and extension of structural elements such as beams. While these modes can be significant in the case of large deformations, and they have no significant effect on the CPU time for very flexible bodies; in the case of thin and stiff structures, the ANCF coupled deformation modes can be associated with very high frequencies that can be a source of numerical problems when explicit integration methods are used. The implicit integration method used in this investigation is the Hilber–Hughes–Taylor method applied in the context of Index 3 differential-algebraic equations (HHT-I3). The results obtained using this integration method are compared with the results obtained using an explicit Adams-predictor-corrector method, which has no adjustable numerical damping. Numerical examples that include bodies with different degrees of flexibility are solved in order to examine the performance of the HHT-I3 implicit integration method when the finite element absolute nodal coordinate formulation is used. The results obtained in this study show that for very flexible structures there is no significant difference in accuracy and CPU time between the solutions obtained using the implicit and explicit integrators. As the stiffness increases, the effect of some ANCF coupled deformation modes becomes more significant, leading to a stiff system of equations. The resulting high frequencies are filtered out when the HHT-I3 integrator is used due to its numerical damping properties. The results of this study also show that the CPU time associated with the HHT-I3 integrator does not change significantly when the stiffness of the bodies increases, while in the case of the explicit Adams method the CPU time increases exponentially. The fundamental differences between the solution procedures used with the implicit and explicit integrations are also discussed in this paper.     

正装结构非线性应力应变累加规律研究

基于有限位移理论的正装结构非线性分析理论对解决超大跨度、柔性结构的非线性计算有重大意义.文中探讨了正装结构非线性的分析特点,研究了其应变场与应力场的Kirchhoff应力张量与Lagrange应变张量的适用性,提出了正装结构非线性分析中应力场与应交场的累加规律,导出了拖动坐标法的虚功增量方程,以此对杆系结构非线性分析常用的CR法和UL列式进行了精度比较分析.文中研究成果可为正装结构的非线性分析提供理论指导.     

利用凸优化进行多用户OFDM通信系统的简化设计

以有限发射功率、迫零均衡为限定条件,利用现代凸优化理论和算法进一步进行多用户通信系统收发器的简化设计。将半正定规划模型等效转化为二阶锥规划模型。二阶锥规划模型将更加简单而且计算复杂度获得很大的改善。借助优化软件SeDuMi的计算机仿真表明:设计的收发器比现存收发     

拉氏数值模拟中的人工粘性方法

较好的人工粘性需要满足较小的计算开销、不能去除真实具有的涡运动等条件.提出一种应用于拉氏数值模拟中基于Lew人工粘性,同时增加了限制器的人工粘性方法.可以有效减少数值模拟结果对网格的依赖;采用特征值限制器控制施加的人工粘性大小,通过限制器能够区分激波压缩和等熵压缩;方便应用在二维、三维,结构网格或者非结构网格上.     

三聚氰胺与鯡鱼精DNA相互作用研究

采用UV、荧光光谱、电化学和粘度法等方法研究了三聚氰胺(MM)与鲱鱼精DNA的作用机制, 并采用凝胶电泳考察了MM对pBR322质粒DNA迁移作用的影响. 当pH为7.0时, 动力学参数的计算结果表明: 加入DNA前后MM电荷转移系数α分别为0.27和0.26, 速率常数k_s分别为5.58 s^-1和5.65 s^－1; 热力学研究表明MM与DNA的结合常数K约为10⁵ mol^－1•L, MM与DNA可能是以沟槽式相结合.     

基于点插值的配点型无网格法解Helmholtz问题

基于点插值法的思想,用三角函数作为基函数在局部支持域内构造具有Kroneckerδ函数性、单位分解性、高阶连续性、再生性和紧支性的形函数.用配点法离散微分方程,得到了具有稀疏带状性的系数矩阵,用GMERS方法求解代数方程组,分别研究了Helmholtz问题的边界层问题和波传播问题.通过数值算例可以发现,给出的数值结果非常接近于精确解,且随着节点的增加,其精确度越来越高,具有良好的收敛性.     

丁二烯气相聚合的负载型稀土催化剂研究

合成了负载型稀土催化剂,并将其用于丁二烯气相聚合.研究表明,硅胶负载Nd(naph)₃-Al(i-Bu)₃-Al(i-Bu)₂Cl催化体系具有相当高的催化活性和立体定向性,其最佳组成:n(Al)/n(Nd)=40～60,n(Cl)/n(Nd)=3～7.添加适量单体丁二烯或异戊二烯,尤其是添加丁二烯,可成倍提高催化活性,随着聚合反应的进行,其速率在10min内迅速增加并达到峰值,随后降低,动力学行为属衰减型.所得聚丁二烯的分子量在几十万至百万,凝胶含量小于6%,顺-1,4结构含量达到98%左右.     

Multisymplectic Geometry for the Seismic Wave Equation

The multisymplectic geometry for the seismic wave equation is presented in this paper. The local energy conservation law, the local momentum evolution equations, and the multisymplectic form are derived directly from the variational principle. Based on the covariant Legendre transform, the multisymplectic Hamiltonian formulation is developed. Multisymplectic discretization and numerical experiments are also explored.     

A continuum model accounting for defect and mass densities in solids with inelastic material behaviour

In this paper the analysis of structures with inelastic material behaviour is considered taking into account the evolution of defects and changes in mass density. The underlying kinematical concept of an oriented continuum is general enough to describe the micro- and macrobehaviour of material bodies appropriately. Based on the logical and consistent variational arguments for a Lagrangian functional the dynamic balance laws, boundary and transversality conditions, all related to the evolution of defect density and mass changes, are derived for macro- and microstresses of deformational as well as of configurational type. The adopted procedure, which formally leaves the balance laws unaltered, leads to the additional balance law for changes in defect density and additional boundary conditions for the changes in mass and defect densities. Driving forces or affinities, associated with the evolution of defect and mass densities, and a generalization of the J-integral representing the thermodynamic forces on defects are obtained. A nonlocal constitutive model accounting for changes in the defect density is presented.