THE INTEGRATION METHODS OF VACCO DYNAMICS EQUATION OFNONLINEAR NONHOLONOMIC SYSTEM

This paper presents the integration methods for vacco dynmmies equations of nonlinear nonholononic system,First.vacco dynamies equations are written in the canonical form and the field form.second the gradient methods the single-componentmethods and the field method are used to integrate the dynamics equations of the corresponding holonomic system respectively.And considering the restriction of nonholonomic construint to the initial conditions the solutions of Vacco dynamics cquations of nonlinear nonholonomic system are obtained.

OPTIMIZATION OF ADAPTIVE DIRECT METHOD FOR APPROXIMATE SOLUTION OF INTEGRAL EQUATIONS OF SEVERAL VARIABLES

This paper determines the exact error order on optimization of adaptive direct methods of approximate solution of the class of Fredholm integral equations of the second kind with kernel belonging to the anisotropic Sobolev classes, and also gives an optimal algorithm.

A PREDICT-CORRECT NUMERICAL INTEGRATION SCHEME FOR SOLVING NONLINEAR DYNAMIC EQUATIONS

A new numerical integration scheme incorporating a predict-correct algorithm for solving the nonlinear dynamic systems was proposed in this paper. A nonlinear dynamic system governed by the equation v = F（v, t） was transformed into the form as v = Hv ＋ f（v, t）. The nonlinear part f（v, t） was then expanded by Taylor series and only the first-order term retained in the polynomial. Utilizing the theory of linear differential equation and the precise time-integration method, an exact solution for linearizing equation was obtained. In order to find the solution of the original system, a third-order interpolation polynomial of v was used and an equivalent nonlinear ordinary differential equation was regenerated. With a predicted solution as an initial value and an iteration scheme, a corrected result was achieved. Since the error caused by linearization could be eliminated in the correction process, the accuracy of calculation was improved greatly. Three engineering scenarios were used to assess the accuracy and reliability of the proposed method and the results were satisfactory.

ON DIRECT METHOD OF SOLUTION FOR A CLASS OF SINGULAR INTEGRAL EQUATIONS

In this article, by introducing characteristic singular integral operator and associate singular integral equations （SIEs）, the authors discuss the direct method of solution for a class of singular integral equations with certain analytic inputs. They obtain both the conditions of solvability and the solutions in closed form. It is noteworthy that the method is different from the classical one that is due to Lu.

Adaptive Reduced Basis Methods Applied to Structural Dynamic Analysis

The reduced basis methods (RBM) have been demonstrated as a promising numerical technique for statics problems and are extended to structural dynamic problems in this paper. Direct step-by-step integration and mode superposition are the most widely used methods in the field of the finite element analysis of structural dynamic response and solid mechanics. Herein these two methods are both transformed into reduced forms according to the proposed reduced basis methods. To generate a reduced surrogate model with small size, a greedy algorithm is suggested to construct sample set and reduced basis space adaptively in a prescribed training parameter space. For mode superposition method, the reduced basis space comprises the truncated eigenvectors from generalized eigenvalue problem associated with selected sample parameters. The reduced generalized eigenvalue problem is obtained by the projection of original generalized eigenvalue problem onto the reduced basis space. In the situation of direct integration, the solutions of the original increment formulation corresponding to the sample set are extracted to construct the reduced basis space. The reduced increment formulation is formed by the same method as mode superposition method. Numerical example is given in Section 5 to validate the efficiency of the presented reduced basis methods for structural dynamic problems.

Towards a Unified Single Analysis Framework Embedded with Multiple Spatial and Time Discretized Methods for Linear Structural Dynamics

We propose a novel computational framework that is capable of employing different time integration algorithms and different space discretized methods such as the Finite Element Method,particle methods,and other spatial methods on a single body sub-dividedintomultiple subdomains.This is in conjunctionwithimplementing thewell known Generalized Single Step Single Solve(GS4)family of algorithms which encompass the entire scope of Linear Multistep algorithms that have been developed over the past 50 years or so and are second order accurate into the Differential Algebraic Equation framework.In the current state of technology,the coupling of altogether different time integration algorithms has been limited to the same family of algorithms such as theNewmarkmethods and the coupling of different algorithms usually has resulted in reduced accuracy in one or more variables including the Lagrange multiplier.However,the robustness and versatility of the GS4 with its ability to accurately account for the numerical shifts in various time schemes it encompasses,overcomes such barriers and allows a wide variety of arbitrary implicit-implicit,implicit-explicit,and explicit-explicit pairing of the various time schemes while maintaining the second order accuracy in time for not only all primary variables such as displacement,velocity and acceleration but also the Lagrange multipliers used for coupling the subdomains.By selecting an appropriate spatialmethod and time scheme on the area with localized phenomena contrary to utilizing a single process on the entire body,the proposed work has the potential to better capture the physics of a given simulation.The method is validated by solving 2D problems for the linear second order systems with various combination of spatial methods and time schemes with great flexibility.The accuracy and efficacy of the present work have not yet been seen in the current field,and it has shown significant promise in its capabilities and effectiveness for general linear dynamics through numerical examples.

Flow Dynamics of a Spiral-groove Dry-gas Seal

The dry-gas seal has been widely used in different industries. With increased spin speed of the rotator shaft, turbulence occurs in the gas film between the stator and rotor seal faces. For the micro-scale flow in the gas film and grooves, turbulence can change the pressure distribution of the gas film. Hence, the seal performance is influenced. However, turbulence effects and methods for their evaluation are not considered in the existing industrial designs of dry-gas seal. The present paper numerically obtains the turbulent flow fields of a spiral-groove dry-gas seal to analyze turbulence effects on seal performance. The direct numerical simulation （DNS） and Reynolds-averaged Navier-Stokes （RANS） methods are utilized to predict the velocity field properties in the grooves and gas film. The key performance parameter, open force, is obtained by integrating the pressure distribution, and the obtained result is in good agreement with the experimental data of other researchers. Very large velocity gradients are found in the sealing gas film because of the geometrical effects of the grooves. Considering turbulence effects, the calculation results show that both the gas film pressure and open force decrease. The RANS method underestimates the performance, compared with the DNS. The solution of the conventional Reynolds lubrication equation without turbulence effects suffers from significant calculation errors and a small application scope. The present study helps elucidate the physical mechanism of the hydrodynamic effects of grooves for improving and optimizing the industrial design or seal face pattern of a dry-gas seal.

Generation of Dynamic Grids and Computation of Unsteady Transonic Flows around Assemblies

Algebraic methods and rapid deforming techniques are used to generate three-dimensional boundary-fitted dynamic grids for assemblies. The conservative full-potential equation is solved by a time-accurate approximate factorization algorithm and internal Newton iterations. An integral boundary layer method based on the dissipation integral is used to account for viscous effects. The computational results about unsteady transonic forces on wings, bodies and control surfaces are in agreement with experimental data.

Dynamic interaction of twin vertically overlapping lined tunnels in an elastic half space subjected to incident plane waves

The scattering of plane harmonic P and SV waves by a pair of vertically overlapping lined tunnels buried in an elastic half space is solved using a semi-analytic indirect boundary integration equation method. Then the effect of the distance between the two tunnels, the stiffness and density of the lining material, and the incident frequency on the seismic response of the tunnels is investigated. Numerical results demonstrate that the dynamic interaction between the twin tunnels cannot be ignored and the lower tunnel has a significant shielding effect on the upper tunnel for high-frequency incident waves, resulting in great decrease of the dynamic hoop stress in the upper tunnel; for the low-frequency incident waves, in contrast, the lower tunnel can lead to amplification effect on the upper tunnel. It also reveals that the frequency-spectrum characteristics of dynamic stress of the lower tunnel are significantly different from those of the upper tunnel. In addition, for incident P waves in low-frequency region, the soft lining tunnels have significant amplification effect on the surface displacement amplitude, which is slightly larger than that of the corresponding single tunnel.

Dynamic rupture process of the 1999 Chi-Chi,Taiwan,earthquake

In this study, we preliminarily investigated the dynamic rupture process of the 1999 Chi-Chi, Taiwan, earthquake by using an extended boundary integral equation method, in which the effect of ground surface can be exactly included. Parameters for numerical modeling were carefully assigned based on previous studies. Numerical results indicated that, although many simplifications are assumed, such as the fault plane is planar and all heterogeneities are neglected, distribution of slip is still consistent roughly with the results of kinematic inversion, implying that for earthquakes in which ruptures run up directly to the ground surface, the dynamic processes are controlled by geometry of the fault to a great extent. By taking the common feature inferred by various kinematic inversion studies as a restriction, we found that the critical slip-weakening distance Dc should locate in a narrow region [60 cm, 70 cm], and supershear rupture might occur during this earthquake, if the initial shear stress before the mainshock is close to the local shear strength.

Conservation Law Classification and Integrability of Generalized Nonlinear Second-Order Equation

This paper is concerned with the generalized nonlinear second-order equation.By the direct construction method,all of the first-order multipliers of the equation are obtained,and the corresponding complete conservation laws(CLs) of such equations are provided.Furthermore,the integrability of the equation is considered in terms of the conservation laws.In addition,the relationship of multipliers and symmetries of the equations is investigated.

Analysis of the role of branching angle in the dynamic rupture process on a 3-D branching fault system

The fault branching phenomenon,which may heavily influence the patterns of rupture propagation in fault systems,is one of the geometric complexities of fault systems that is widely observed in nature.In this study,we investigate the effect of the branching angle on the rupture inclination and the interaction between branch planes in two-fork branching fault systems by numerical simulation and theoretical analysis based on Mohr’s circle.A friction law dependent on normal stress is used,and special attention is paid to studying how ruptures on the upper and lower branch planes affect the stress and rupture on each other separately.The results show that the two branch planes affect each other in different patterns and that the intensity of the effect changes with the branching angle.The rupture of the lower branch plane has a negative effect on the rupture of the upper branch plane in the case of a small branching angle but has almost no negative effect in the case of a large branching angle.The rupture of the upper branch plane,however,suppresses the rupture of the lower branch plane regardless of whether the branching angle is large or small.

Exact Solutions of the Harry-Dym Equation

The aim of this paper is to generate exact travelling wave solutions of the Harry-Dym equation through the methods of Adomian decomposition, He＇s variational iteration, direct integration, and power series. We show that the two later methods are more successful than the two former to obtain more solutions of the equation.

Classification of All Single Traveling Wave Solutions to (3 + 1)-Dimensional Breaking Soliton Equation

In order to get the exact traveling wave solutions to nonlinear partial differential equation, the complete discrimination system for polynomial and direct integral method are applied to the considered equation. All single traveling wave solutions to the equation can be obtained. As an example, we give the solutions to (3 + 1)-dimensional breaking soliton equation.

A two-step unconditionally stable explicit method with controllable numerical dissipations

A family of unconditionally stable direct integration algorithm with controllable numerical dissipations is proposed. The numerical properties of the new algorithms are controlled by three parameters α, β and γ. By the consistent and stability analysis, the proposed algorithms achieve the second-order accuracy and are unconditionally stable under the condition that α≥-0.5, β≤ 0.5 and γ≥-(1+α)/2. Compared with other unconditionally stable algorithms, such as Chang's algorithms and CR algorithm, the proposed algorithms are found to be superior in terms of the controllable numerical damping ratios. The unconditional stability and numerical damping ratios of the proposed algorithms are examined by three numerical examples. The results demonstrate that the proposed algorithms have a superior performance and can be used expediently in solving linear elastic dynamics problems.

High accuracy compact finite difference methods and their applications

Numerical simulation of complex flow fields with multi-scale structures is one of the most important and challenging branches of computational fluid dynamics. From linear analysis and numerical experiments it has been discovered that the higher-order accurate method can give reliable and efficient computational results, as well as better resolution of the complex flow fields with multi-scale structures. Compact finite difference schemes, which feature higher-order accuracy and spectral-like resolution with smaller stencils and easier application of boundary conditions, has attracted more and more interest and attention.

双重孔隙流体饱和介质弹性波散射二维IBIEM模拟

基于平面波势函数,采用间接边界积分方程法(indirect boundary integral equation method, IBIEM)研究了双孔隙流体饱和介质中弹性波入射下二维孔洞的散射特性。推导得到了双孔隙介质中全空间二维线源动力格林函数,并给出了各散射波的位移场和应力场。在数值精度验证的基础上,以双孔隙二维饱和全空间中孔洞为例,解决了平面P、SV波入射下的地震波散射问题。数值结果表明:双重孔隙介质中的位移幅值、环向应力幅值、孔隙压力变化规律与不同入射波形,入射频率,孔隙率和边界排水条件密切相关,位移幅值在低频(无量纲频率η≤2)入射时出现峰值。环向应力幅值与干土条件相比更为复杂,基质孔压与裂缝孔压的存在增大了双重孔隙饱和介质的能量效应,总体震动趋势大于干土条件,环向应力放大可达62%。

近断层脉冲型随机地震动模拟方法及结构动力可靠度研究

以集集脉冲型地震动记录为数据集建立了一种改进的参数化近断层脉冲型随机地震动模型.首先,采用MP03模型和人工蜂群算法从实际地震动记录提取低频脉冲部分,建立随机低频脉冲模型;采用谱表示-随机函数法生成高频部分,二者叠加生成近断层脉冲型随机地震动.然后,采用蒙特卡罗模拟(MCS)法和直接概率积分法(DPIM)分别研究非线性单自由度体系和高层框架结构在不同类型地震动作用下的动力可靠度.结果表明:DPIM与MCS法计算结果相比基本吻合,但计算效率得到较大提升,可以很好地用于非线性结构的时变可靠度研究;脉冲型地震动与非脉冲型地震动相比能够对结构造成更大的变形,极大降低了结构的可靠度,且近断层距越小影响越显著,低频脉冲是使结构可靠度降低的主要因素,而高频部分对结构可靠度的影响较小.

考虑碳-绿证耦合机制的电-气-热综合能源系统分布式鲁棒低碳经济调度

为了保护综合能源系统中不同运营主体的信息隐私,权衡系统运行低碳性和经济性,解决天然气网络鲁棒模型难以求解的问题,提出了考虑碳-绿证耦合机制的电-气-热综合能源系统分布式鲁棒低碳经济调度方法。考虑天然气网络和热力网络的动态特性、碳-绿证耦合机制,构建综合能源系统动态低碳经济调度模型;为了保护各网络运营商的数据隐私,根据能量耦合关系对综合能源系统解耦,建立电-气-热网络的分布式协同优化模型;在此基础上,考虑风电出力和多能负荷的不确定性,提出基于一致性交替方向乘子法的分布式鲁棒优化框架,并利用二阶锥对偶理论与交替优化方法实现含二阶锥约束和二元变量的鲁棒子问题求解。以修改的IEEE 39节点电力网络、比利时20节点天然气网络和15节点热力网络为算例进行仿真分析,验证所提方法能够在源荷不确定性条件下实现综合能源系统各网络分散自治运行,同时兼顾系统运行的低碳性和经济性。

综合能源系统建模与仿真:综述、思考与展望

建模和仿真将真实系统及其物理性质抽象化为适当的数学形式,并采用数学方法求解以模拟真实系统行为或性能。为综合能源系统在各时间尺度的物理过程建立恰当的数学模型并进行仿真验证,是分析系统运行规律、优化运行策略以及评估运行性能的基础。文中围绕电气热综合能源系统建模和仿真问题,首先介绍了电、气、热子系统和耦合设备的机理模型,总结了建模的研究进展与技术难点。进一步,从仿真框架、算法以及改进策略3个方面介绍了综合能源系统仿真技术所取得的研究进展。特别地,总结了仿真算法中数值法、半解析法、解析法的原理、应用现状及其难点。最后,围绕现存技术难点,从建模和仿真算法2个方面对未来研究进行了展望。