Modified precise time step integration method of structural dynamic analysis

The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method （e.g., an algorithm with an arbitrary order of accuracy） is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.

An Explicit-Implicit Mixed Staggered Asynchronous Step Integration Algorithm in Structural Dynamics

Many engineering applications need to analyse the system dynamics on the macro and micro level,which results in a larger computational effort.An explicit-implicit asynchronous step algorithm is introduced to solve the structural dynamics in multi-scale both the space domain and time domain.The discrete FEA model is partitioned into explicit and implicit parts using the nodal partition method.Multiple boundary node method is adopted to handle the interface coupled problem.In coupled region,the implicit Newmark coupled with an explicit predictor corrector Newmark whose predictive wave propagates into the implicit mesh.During the explicit subcycling process,the variables of boundary nodes are solved directly by dynamics equilibrium equation.The dissipation energy is dynamically determined in accordance with the energy balance checking.A cantilever beam and a building two numerical examples are proposed to verify that the method can greatly reduce the computing time while maintaining a high accuracy.

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.

Precise integration method without inverse matrix calculation for structural dynamic equations

The precise integration method proposed for linear time-invariant homogeneous dynamic systems can provide accurate numerical results that approach an exact solution at integration points. However, difficulties arise when the algorithm is used for non-homogeneous dynamic systems due to the inverse matrix calculation required. In this paper, the structural dynamic equalibrium equations are converted into a special form, the inverse matrix calculation is replaced by the Crout decomposition method to solve the dynamic equilibrium equations, and the precise integration method without the inverse matrix calculation is obtained. The new algorithm enhances the present precise integration method by improving both the computational accuracy and efficiency. Two numerical examples are given to demonstrate the validity and efficiency of the proposed algorithm.

ELEMENT－BY－ELEMENT MATRIX DECOMPOSITION ANDSTEP－BY－STEP INTEGRATION METHOD FOR TRANSIENTDYNAMIC PROBLEMS

In this paper a general matrix decomposition scheme as well as an element-by-clement relaxation algorithm combined with step-by -step integration method is presented for transient dynamic problems thus the finite element method can be fromforming global stiffness matrix global mass matrix as well as solyin large scale sparse equations Theory analysis and numerical results show that the presented matrix decomposition scheme is the optimal one The presented algoithm has else physicalmeaning and can be busily applied to finite element codes

A new simple method of implicit time integration for dynamic problems of engineering structures

This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditional and unconditional stability regions of the parameters. When the control parameters in the method are optimally taken in their regions, the accuracy may be improved to reach of order three. It is found that the new scheme can achieve lower numerical amplitude dissipation and period dispersion than some of the existing methods, e.g. the Newmark method and Zhai＇s approach, when the same time step size is used. The region of time step dependent on the parameters in the new scheme is explicitly obtained. Finally, some examples of dynamic problems are given to show the accuracy and efficiency of the proposed scheme applied in dynamic systems.

COORDINATION OF THE MOTION PARAMETERS WITHIN STEP-BY-STEP INTEGRATION FOR DYNAMIC EQUATION

A method is presented that coordinates the calculation of the displacement, velocity and acceleration of structures within the time-steps of different types of step-by-step integration. The dynamic equation is solved using an energy equation and the calculating data of the original method. The method presented is better than the original method in terms of calculating postulations and is in better conformity with the system's movement. Take the Wilson-θ method as an example. By using the coordination process, the calculation precision has been greatly im proved (reducing the errors by approximately 90% ), and the greater part of overshooting of the calculation result has been eliminated. The study suggests that the mal-coordination of the motion parameters within the time-step is the major factor that contributes to the result errors of step-by-step integration for the dynamic equation.

On time-step in structural seismic response analysis under ground displacement/acceleration

There are two models in use today to analyze structural responses when subjected to earthquake ground motions, the Displacement Input Model （DIM） and the Acceleration Input Model （AIM）. The time steps used in direct integration methods for these models are analyzed to examine the suitability of DIM. Numerical results are presented and show that the time-step for DIM is about the same as for AIM, and achieves the same accuracy. This is contrary to previous research that reported that there are several sources of numerical errors associated with the direct application of earthquake displacement loading, and a very small time step is required to define the displacement record and to integrate the dynamic equilibrium equation. It is shown in this paper that DIM is as accurate and suitable as, if not more than, AIM for analyzing the response of a structure to uniformly distributed and spatially varying ground motions.

RETHINKING TO FINITE DIFFERENCE TIME-STEP INTEGRATIONS

The numerical time step integrations of PDEs are mainly carried out by the finitedifference method to date. However,when the time step becomes longer, it causes theproblem of numerical instability,. The explicit integration schemes derived by the singlepoint precise integration method given in this paper are proved unconditionally stable.Comparisons between the schemes derived by the finite difference method and theschemes by the method employed in the present paper are made for diffusion andconvective-diffusion equations. Nunierical examples show the superiority of the singlepoint integration method.

Applying the Dynamic Two-Step Method to Forecast Remaining Oil Distribution of Lower Series ,Xiaermen Oilfield

The distribution of remaining oil is often described qualitatively. The remaining oil distributed in the whole reservoir is calculated according to the characteristics of the space distribution of the saturation of remaining oil. Logging data are required to accomplish this. However, many such projects cannot be completed. Since the old study of remaining oil distribution could not be quantified efficiently, the ＂dynamic two-step method＂ is presented. Firstly, the water cut of every flow unit in one well at one time is calculated according to the comprehensive water cut of a single well at one time. Secondly, the remaining oil saturation of the flow unit of the well at one time is calculated based on the water cut of the flow unit at a given time. The results show that ＂dynamic two-step method＂ has characteristics of simplicity and convenience, and is especially suitable for the study of remaining oil distribution at high water-cut stage. The distribution of remaining oil presented banding and potato form, remaining oil was relatively concentrated in faultage neighborhood and imperfect well netting position, and the net thickness of the place was great. This proposal can provide an effective way to forecast remaining oil distribution and enhance oil recovery, especially applied at the high water-cut stage.

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.

Precise integration methods based on the Chebyshev polynomial of the first kind

This paper introduces two new types of precise integration methods based on Chebyshev polynomial of the first kind for dynamic response analysis of structures, namely the integral formula method （IFM） and the homogenized initial system method （HISM）. In both methods, nonlinear variable loadings within time intervals are simulated using Chebyshev polynomials of the first kind before a direct integration is performed. Developed on the basis of the integral formula, the recurrence relationship of the integral computation suggested in this paper is combined with the Crout decomposed method to solve linear algebraic equations. In this way, the IFM based on Chebyshev polynomial of the first kind is constructed. Transforming the non-homogenous initial system to the homogeneous dynamic system, and developing a special scheme without dimensional expansion, the HISM based on Chebyshev polynomial of the first kind is able to avoid the matrix inversion operation. The accuracy of the time integration schemes is examined and compared with other commonly used schemes, and it is shown that a greater accuracy as well as less time consuming can be achieved. Two numerical examples are presented to demonstrate the applicability of these new methods.

NUMERICAL METHOD FOR THREE-DIMENSIONAL NONLINEAR CONVECTION-DOMINATED PROBLEM OF DYNAMICS OF FLUIDS IN POROUS MEDIA

For the three-dimensional convection-dominated problem of dynamics of fluids in porous media, the second order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward. Fractional steps techniques are needed to convert a multi-dimensional problem into a series of successive one-dimensional problems. Some techniques, such as calculus of variations, energy method, multiplicative commutation rule of difference operators, decomposition of high order difference operators, and the theory of prior estimates are adopted. Optimal order estimates are derived to determine the error in the second order approximate solution. These methods have already been applied to the numerical simulation of migration-accumulation of oil resources and predicting the consequences of seawater intrusion and protection projects.

Applying a Step Approach Method in Solving the Multi-Frequency Radiation From a Complex Obstacle

In this paper,a step approach method in the time domain is developed to calculate the radiated waves from an arbitrary obstacle pulsating with multiple frequencies.The computing scheme is based on the Boundary Integral Equation and derived in the time domain;thus,the time-harmonic Neumann boundary condition can be imposed.By the present method,the values of the initial conditions are set to zero,and the approach process is carried forward in a loop from the first time step to the last.At each time step,the radiated pressure on each element is updated.After several loops,the correct radiated pressures can be obtained.A sphere pulsating with a monopole frequency in an infinite acoustic domain is calculated first.This result is compared with the analytical solution,and both of them are in good agreement.Then,a complex-shaped radiator is taken as the studied case.The pulsating frequency of this case is multiple,and the waves propagate in half space.It is shown that the present method can treat multiple-frequency pulsation well,even when the radiator is a complex shape,and a robust convergence can be attained quickly.

ZnO micro-nano composite hydrophobic film prepared by the three-step method

The hydrophobicity of the lotus leaf is mainly due to its surface micro-nano composite structure. In order to mimic the lotus structure, ZnO micro-nano composite hydrophobic films were prepared via the three-step method. On thin buffer films of SiO2, which were first fabricated on glass substrates by the so,gel dip-coating method, a ZnO seed layer was deposited via RF magnetron sputtering. Then two different ZnO films, micro-nano and micro-only flowerlike structures, were grown by the hydrothermal method. The prepared films have different hydrophobic properties after surface modification. The structures of the obtained ZnO films were characterized using x-ray diffraction and field-emission scanning electron microscopy. A conclusion that a micro-nano composite structure is more beneficial to hydrophobicity than a micro-only structure was obtained through research into the effect of structure on hydrophobic properties.

基于STEP的三维CAD和PDM系统集成方法研究

研究了基于STEP标准的产品信息模型、CAD数据交换接口、产品数据交换接口,提出了基于STEP中性机制和产品信息模型的三维CAD和PDM系统集成框架,该框架有助于实现广义CAD功能CAD/CAE/CAPP/CAM/CAQ/PDM/ERP之间的系统集成,异地、异构系统在企业间的集成,能够为企业设计和管理信息化提供一体化解决方案。

基于STEP标准的BOM数据结构表达

针对BOM的研究对制造业信息化进程及效果的重要性,结合某动力股份有限公司CIMS应用示范工程的实际需要,研究如何通过STEP标准的部分实体来描述BOM数据结构。介绍了表达BOM数据结构所需引用的两个STEP集成资源构件,产品定义模式和产品结构模式,然后通过一个机械装配产品的范例来详细阐述了基于STEP标准描述和表达BOM数据结构的过程。

aecXML与STEP在建筑业信息集成中的应用

利用web实现建筑业信息集成是建筑业信息化的发展方向,而平台的不统一性和数据的异构性是信息集成化发展的瓶颈。综合运用aecXML和STEP技术实现信息的标准化传递,可以克服这个问题。文章介绍了aecXML和STEP的基本概念及其各自的优缺点,探讨了基于aecXML和STEP的建筑产品数据共享和交换结构模型,以及基于aecXML的建设项目信息共享与交换结构模型,并对两种技术的综合应用前景进行了展望。

Adaptive explicit Magnus numerical method for nonlinear dynamical systems

Based on the new explicit Magnus expansion developed for nonlinear equations defined on a matrix Lie group, an efficient numerical method is proposed for nonlinear dynamical systems. To improve computational efficiency, the integration step size can be adaptively controlled. Validity and effectiveness of the method are shown by application to several nonlinear dynamical systems including the Duffing system, the van der Pol system with strong stiffness, and the nonlinear Hamiltonian pendulum system.

中外合作办学人才培养新模式研究——以“混凝土结构”课程为例

随着国家对中外合作办学的支持,越来越多的高校致力于对中外合作办学模式的探讨。青岛理工大学土木工程学院以“混凝土结构”课程为例,针对中美合作办学培养新模式构建了经济全球化背景下源源不断培养国际人才的教学体系:中外教师先后对同一门课程进行线下课堂与远程授课紧密结合的“三位一体”教学,实行“分步走一体化”的流程,保证学生能够立足于国内外的混凝土建筑案例,夯实专业基础,充分拓宽国际视野。与传统中外合作办学模式相比,新型模式避免了赴外求学可能面临的种种风险,提高了学生对外文专业知识的了解程度。最终达到培养着眼于国际视野的高素质技术人才的目标,保证各种突发情况下国际化人才教育培养不停摆。