Generalized synchronization of spatiotemporal chaos in a weighted complex network

The nodes of the network are composed of the spatiotemporal chaos systems. The relations between the nodes are built through a weighted connection and the nonlinear terms of the chaos systems themselves are taken as coupling functions. The structure of the coupling functions between the connected nodes and the range of the control gain are obtained based on Lyapunov stability theory. It is proven that generalized chaos synchronization of the weight complex network can be realized even if the coupling strength between the nodes is adopted as any weight value. Subsequently, the catalytic reaction diffusion system which has spatiotemporal chaos behavior is taken as example, and simulation results show the effectiveness of the synchronization principle.     

自然单元法计算裂纹与材料边界问题

提出一种新的非凸边界上自然单元法形函数计算方法,通过边界结点限制点对间的邻点关系,对包括裂纹和材料边界在内的各种类型的非凸边界具有统一的处理原则,所得到的近似函数在边界结点间具有线性插值性.     

重分析边界处理与广义主从控制关系

本文讨论了在结构重分析中局部网格细分时，修改区域与未修改区域交界处的位移连续性问题．文中提出了广义主从控制的概念，把边界上的新增节点作为从节点，同时受多点控制，保证了结构的连续性，并对有限元结构中的常见单元给出了边界的位移控制关系，而且此关系已应用于DDJ程序系统，实现了算法及程序设计．     

Attack robustness of cascading model with node weight

Considering the node weight and the effect of the neighboring nodes, we introduce a method to define the initial load of a node in a network and propose a cascading model. Our aim is to explore how to allocate the initial load and select some nodes to be protected so as to maximize the network robustness against cascading failures. According to the distribution of the node weight and the normalized effect of the neighboring nodes, we focus on investigating different roles of high-load and low-load nodes and the correlation between some parameters in our model and the strongest robust level of a network against cascading failures. We obtain by the numerical simulations the optimal values of the parameters in our model at which the network can reach the strongest robust level against cascading failures. In addition, we find that the effect of two attacks strongly depends on the parameters of the node weight and the normalized effect of the neighboring nodes, i.e., the network robustness of attacking the low-load nodes has a positive correlation with the parameter of the node weight, while has a negative correlation with the parameter of the normalized reaction of the neighboring nodes. While the result of attacking the high-load nodes is almost on the contrary. Finally, we verify the numerical results by the theoretical analysis.     

Spatiotemporal chaos synchronization between uncertain complex networks with diverse structures

Spatiotemporal chaos synchronization between uncertain complex networks with diverse structures is investigated. The identification law of unknown parameters and the adaptive law of the configuration matrix element in state equations of network nodes are determined based on stability theory, and the conditions of realizing spatiotemporal chaos synchronization between uncertain complex networks with different structures are discussed and obtained. Further, the Fisher–Kolmogorov system with spatiotemporal chaotic behavior is taken as the nodes of drive and response networks to imitate the experiment. It is found that the synchronization performance between two networks is very stable.     

多尺度六面体单元网格连接方法研究

提出多结点六面体单元的结点形函数构建方法,解决了数字化分析时两个多尺度有限元网格在接触界面由于单元结点不匹配而导致的结点属性不能连续传递问题.首先将全局坐标下的六面体单元及其表面上多结点通过等参逆变换转成局部坐标下的规则六面体单元及其表面上多结点;在规则单元中,以每个结点为基点,分别沿三个正交的局部坐标方向在单元内寻找...     

非匹配结点的有限单元等参插值方法及其应用研究

针对两零件的异构网格单元结点在接触界面不能相互匹配导致结点属性不能连续过渡和传递的问题,提出非匹配结点的有限单元等参插值方法,通过构建非匹配结点的形函数和修正原结点的形函数,将结点属性值的影响范围限制在可控的局部区域,从而实现两异构网格结点属性在接触界面的连续过渡和传递。通过两个异构的四边形单元网格的结点属性在接触界面的过渡实例和啮合齿轮的接触分析应用,验证了该方法的正确性和有效性。     

一种适用于下承式钢桁结合桥计算的板梁组合单元

下承式钢桁结合桥具有建筑高度较低、行车时燥声、震动较小、刚度较大等优点,是高速铁路桥梁中比较理想的结构形式之一。为了能较合理和方便地分析下承式钢桁结合桥的力学行为,本文提出了一种考虑板梁共同作用和相对滑移的板梁组合单元。该单元以矩形板的4个结点的结点位移和板梁结合面3个点的相对滑移作为单元的基本自由度,板的位移模式通过对常规4结点矩形平板壳单元的形函数修正得到,梁的位移模式则根据板和梁之间的变形协调条件来确定。文末的算例表明,本文的解与试验结果和通用工程软件解较为吻合,说明了本文理论和方法的正确性。     

A binary-tree based model for rate-independent polycrystals

A model of a rigid-plastic rate-independent polycrystalline aggregate wherein sub-aggregates are represented as the nodes of a binary tree is proposed. The lowest nodes of the binary tree represent grains. Higher binary tree nodes represent increasingly larger sub-aggregates of grains, culminating with the root of the tree, which represents the entire polycrystalline aggregate. Planar interfaces are assumed to separate the sub-aggregates represented by nodes in the binary tree. Equivalence between the governing equations of the model and a standard linear program is established. The objective function of the linear program is given by the plastic power associated with polycrystal deformation and the linear constraints are given by compatibility requirements between the sub-aggregates represented by sibling nodes in the binary tree. The deviatoric part of the Cauchy stress in each sub-aggregate is deduced as linear combinations of the Lagrange multipliers associated with the constraints. It is shown that the present model allows a generalization of Taylor’s principle to polycrystals. The proposed model is applied to simulate tensile, compressive, torsional, and plane-strain deformation of copper polycrystals. The predicted macroscopic response is in good agreement with published experimental data. The effect of the initial distribution of the planar interfaces separating the sub-aggregates represented by the binary tree on the predicted mechanical response in tension, compression and torsion is studied. Also, the role of constraints relaxation in simulations of plane strain compression is investigated in detail.     

利用高阶分区并行算法实现直接数值模拟

在基于MPI环境的分布式内存机群上,结合高阶WENO-RF格式的特点,实现了5阶WENO-RF格式的分区并行计算方法,计算精度不受分区和节点数量影响。使用该分区并行算法以三维可压缩时间发展混合层为例进行了直接数值模拟,验证了并行算法的准确性,表明机群并行运算可以显著扩展微机的计算能力,并行效率高,减少了计算的墙上时间,适合在小型高速局域网内进行大规模数值模拟计算。     

Extended truss theory with simplex constraints

This paper traces a way of generalization of the classical truss theory: in addition to the kinematic constraint expressing the distance between two nodes connected by a bar element, other similar constraints involving three and four nodes are introduced. Derived from energy principles, a general tangent stiffness formulation is given. Possible mechanical interpretations as well as problems of pre-stressing are also discussed.     

对流扩散方程的绝对稳定高阶中心差分格式

将作者提出的数值摄动算法改进为区分离散单元内上游和下游并分别对通量进行高精度重构的双重数值摄动算法,与原(单重)摄动算法相比,双重摄动算法既提高了格式精度又明显扩大了格式的稳定域范围.利用双重摄动算法,即分别利用上游和下游基点变量的摄动重构将高阶流体力学关系及迎风机制耦合进二阶中心格式之中,由此构建了对流扩散方程的对网格Reynolds数的任意值均稳定(绝对稳定)高精度(四阶和八阶精度)三基点中心TVD差分格式,通过解析分析以及3个算例计算证实了构建格式的优良性能;3个算例包括一维线性、非线性(Burgers方程)和二维变系数对流扩散方程.数值计算表明:构建的格式在粗网格下不振荡,构建格式在粗网格时的最大误差L_∞和均方误差L_2与二阶中心格式在细网格时的相应误差一致,对线性方程,构建格式在细网格下可达到L_2精度阶.     

Configurational forces in discrete elastic systems

The concept of configurational forces is applied to a simple, one-dimensional problem that is solved by finite elements. Both the exact solution and its finite-element approximation are provided in closed form. The total energy according to the approximate solution depends on the choice of the nodes. Any virtual shift of a node results in a virtual change of energy, which can be interpreted as the virtual work done by a configurational force acting on that node. It is shown that, in the case of equidistant nodes, the configurational forces acting on the interior nodes vanish. Also, the relation between the nodal configurational forces and the Eshelby stress resultant along the rod is investigated.     

Parallel ILU preconditioning,a priori pivoting and segregation of variables for iterative solution of the mixed finite element formulation of the Navier–Stokes equations

A parallel ILU preconditioning algorithm for the incompressible Navier–Stokes equations has been designed, implemented and tested. The computational mesh is divided into N subdomains which are processed in parallel in different processors. During ILU factorization, matrices and vectors associated with the nodes on the interface between the subdomains are communicated to the equation matrices to the adjacent subdomain. The bases for the parallel algorithm are an appropriate node ordering scheme and a segregation of velocity and pressure degrees of freedom. The inner nodes of the subdomain are numbered first and then the nodes on the interface between the subdomains. To avoid division by zero during the ILU factorization, the equations corresponding to the velocity degrees of freedom are assembled first in the global equation matrix, followed by the equations corresponding to the pressure degrees of freedom. Copyright © 2003 John Wiley & Sons, Ltd.     

Generalized synchronization of the fractional-order chaos in weighted complex dynamical networks with nonidentical nodes

A fractional-order weighted complex network consists of a number of nodes, which are the fractional-order chaotic systems, and weighted connections between the nodes. In this paper, we investigate generalized chaotic synchronization of the general fractional-order weighted complex dynamical networks with nonidentical nodes. The well-studied integer-order complex networks are the special cases of the fractional-order ones. Based on the stability theory of linear fraction-order systems, the nonlinear controllers are designed to make the fractional-order complex dynamical networks with distinct nodes asymptotically synchronize onto any smooth goal dynamics. Numerical simulations are provided to verify the theoretical results. It is worth noting that the synchronization effect sensitively depends on both the fractional order ?? and the feedback gain k _i . Moreover, generalized synchronization of the fractional-order weighted networks can still be achieved effectively with the existence of noise perturbation.     

基于8节点平面单元的RC梁柱节点单元

为了分析RC框架结构的非线性滞回性能,基于平面8节点单元,本文提出了一个新的针对受循环荷载作用钢筋混凝土的梁柱节点单元．单元中梁与节点交界面和柱与节点交界面被划分成“节点截面”和“梁柱截面”,节点核心区的力学性能由8节点单元描述,而梁柱受力钢筋与节点核心区的粘结滑移由存在于“节点截面”和“梁柱截面”之间的8根弹簧控制,梁柱与节点之间的剪切由4根剪切弹簧表示．单元具有4个外节点和8个内节点,每个内节点具有2个自由度,每个外节点具有3个自由度,该3个自由度与普通梁单元一致,从而确保本单元能够同普通一维梁柱单元一起进行钢筋混凝土结构平面非线性分析．通过将内节点上的自由度依附到外节点上,单元在数值表现上具有4个节点和28个自由度．通过对比试验和模拟分析结果,验证了本模型适合于循环荷载作用下平面框架结构的非线性响应分析．     

NEW ALGORITHM OF COUPLING ELEMENT-FREE GALERKIN WITH FINITE ELEMENT METHOD

IntroductionMeshless methods, as a special numerical method, originated from1970s. Since thediffuse element method was proposed by Nayroleset al.[1]in1992, the meshless methodshave received wide attentions in the mechanics area, and have shown obvious adv…     

A 2-D lattice model for simulating the failure of paper

A new two-dimensional network model is proposed as a micromechanics model to simulate paper’s failure process due to sequentially breakages of fibers and/or bonds. Paper is approximated as a network composed of fibers any two of which link to each other by their intersecting point, namely so-called bond. Fibers distribute along three particular directions, leading to network’s macro-level isotropy. In the framework of finite element method, nodes correspond to fiber-to-fiber bonds, while elements are fiber segments between every two neighboring nodes and described by Timoshenko beam theory. Element breaks when its equivalent internal tensile stress reaches the tensile strength of fiber. Strength of nodes, i.e. fiber-to-fiber bonds is assumed to be dependant on shearing interaction between fibers, considering the dominant interaction is shearing in a plane problem. Numerical examples show the model’s capacity of reflecting basic failure characteristic in paper. Influences of fiber length and the ratio of fiber strength to bond strength are analyzed in detail.     

两套节点格林元嵌入式离散裂缝模型数值模拟方法

对于原始嵌入式离散裂缝模型(EDFM),在计算包含裂缝单元的基质网格内的压力分布时采用了线性分布假设,这导致了油藏开发早期对非稳态窜流量的计算精度不足.因此,本文提出了一种两套节点格林元法的EDFM数值模拟方法.两套节点格林元法的核心思想是将压力节点与流量节点区分开,一套压力节点设置在单元顶点,另一套流量节点设置在网格边的中点,满足局部物质守恒、具有二阶精度的同时,可适用于任意网格类型.本文将两套节点格林元法与EDFM耦合,采用了非稳态渗流控制方程的边界积分形式推导了基质网格与裂缝网格之间传质量的新格式,代替了线性分布假设以提高模拟精度;此外,修正后的EDFM能适应任意形态的基质网格剖分,拓展了原始EDFM仅适用于矩形基质网格、难以考虑复杂油藏边界的局限性.研究表明:通过对比商业模拟软件tNavigator?LGR模块与原始EDFM,验证了本文模型具有较高的早期计算精度;以复杂油藏边界-缝网-SRV分区模型为例,通过对比SFEM-COMSOL商业模拟软件,验证了本文模型处理复杂问题的适应性.本文研究可用于裂缝性油藏开发动态的精确模拟.     

AN h-TYPE ADAPTIVE FINITE ELEMENT

ＡＮｈ－ＴＹＰＥＡＤＡＰＴＩＶＥＦＩＮＩＴＥＥＬＥＭＥＮＴ￥ＸｕＸｉｎｇ（徐兴）ＬｉｎｇＤａｏｓｈｅｎｇ（凌道盛）ＤｕＱｉｎｇｈｕａ（杜庆华）ＤｉｎｇＨａｏｊｉａｎｇ（丁皓江）（ＺｈｅｊｉａｎｇＵｎｉｖｅｒｓｉｔｙ．Ｈａｎｇｚｈｏｕ３１００２７．Ｐ．...