无限长圆柱体中广义电磁热弹耦合问题研究

基于L-S广义热弹性理论,研究了处于磁场中无限长理想圆柱导体在边界受热冲击作用时的电磁热弹耦合问题.建立了广义电磁热弹耦合的有限元方程,为避免积分变换方法求解带来的精度丟失.采用将有限元方程直接在时间域求解的方法,得到了圆柱体中的温度、位移、应力、感应磁场和感应电场的分布规律,反映了热的波动性及电磁热弹的耦合效应.结果表明,将有限元方程直接在时间域求解,可以获得各物理量的准确分布.得到温度在热波波前处的阶跃,准确地反应热波的波动效应.     

Branching of an axially symmetric convective flow

The stability of an axisymmetric convective gas flow with respect to finite disturbances of the bottom temperature is studied numerically using a finite-difference method. The convective flow considered approximately describes the free convection developing in the atmosphere due to heating of the substrate surface. The temperature disturbance used increases the intensity of one of two possible flows and suppresses the other flow, with the opposite signs of the vortices. Using the methods of numerical experiment, the corresponding problem of branching of the solution is examined. It is shown that the transformation of one flow into the other far from the threshold of the onset of convection requires substantial disturbances of the temperature.     

有限覆盖径向点插值方法理论及其应用

数值流形方法能够统一地处理连续与非连续变形问题,有限覆盖技术是这种方法的核心。无网格方法前处理过程比较简单,径向点插值法是其中的一种计算格式。本文将有限覆盖技术与径向点插值方法相结合发展了有限覆盖径向点插值无网格方法,综合了数值流形方法与点插值方法的各自优点,能够有效地处理连续与非连续性问题,由此所构造的形函数具有Kronecker δ-函数属性,能够有效地处理位移边界条件。本文在阐述了这种方法基本原理的基础上,通过算例分析与数值计算论证了本文所建议方法的可靠性及其有效性。     

Time discretization methods in the computation of unsteady flow

对于不同非定常流动问题, 采用合适的时间离散方法,可有效提高数值精度和计算效率. 本文在总结传统时间离散方法的基础上,对近些年发展的非线性频域法、谐波平衡法、经典时间谱方法、时间谱元法、时间有限差分法等进行了系统地总结.根据离散形式的不同,将上述方法分为时域推进法、频域谐波法、时域配点法和混合方法4大类.首先简要介绍了各类方法的数学思想以及研究进展,并重点比较了(准)周期性非定常流动计算中各方法的精度、效率以及适用范围.然后, 对各种时间离散格式的特点进行总结,并就不同的非定常流动问题如何选择合适的时间离散方法给予了建议.最后, 对这些新型时间离散格式在工程中的应用进行了简要介绍,并对其发展方向进行展望.     

Thermal reliability analysis of a BLDC motor in a high-speed axial fan by the accelerated-life test and numerical methods

The thermal reliability of a closed-type BLDC motor for a high-speed fan is analyzed by an accelerated-life testing and numerical methods in this paper. Since a module and a motor part are integrated in a closed case, heat generated from a rotor in a motor and electronic components in the PCB module cannot be effectively removed to the outside. Therefore, the module will easily fail due to high temperature. The experiment for measuring the temperature and the surface heat flux of the electronic components is carried out to predict their surface temperature distributions and main heat sources. The accelerated-life test is accomplished to formulate the life equation depending on the environmental temperature. Moreover, the temperature of the PCB module is different from the environmental temperature since the heat generated from the motor cannot be effectively dissipated owing to the motor’s structure. Therefore a numerical method is used to predict the temperature of the PCB module, which is one of the life equation parameter, according to the environment. By numerically obtaining the maxima of the thermal stress and strain of the electronic components according to the operation environments with the temperature results, the fatigue cycle can be estimated.     

Non-Fourier heat conduction in an exponentially graded slab

The present article investigates one-dimensional non-Fourier heat conduction in a functionally graded material by using the differential transformation method. The studied geometry is a finite functionally graded slab, which is initially at a uniform temperature and suddenly experiences a temperature rise at one side, while the other side is kept insulated. A general non-Fourier heat transfer equation related to the functionally graded slab is derived. The problem is solved in the Laplace domain analytically, and the final results in the time domain are obtained by using numerical inversion of the Laplace transform. The obtained results are compared with the exact solution to verify the accuracy of the proposed method, which shows excellent agreement.     

Study on generalized magneto-thermoelastic problems by FEM in time domain

This paper presents an investigation of temperature, displacement, stress, and induced magnetic field in a half space perfectly-conductive plate. Finite element equations regarding generalized magneto-thermoelasticity problems with two relaxation times (i.e., the G-L theory) are derived using the principle of virtual work. For avoiding numerical complication involved in inverse Laplace and Fourier transformation and low precision thereof, the equations are solved directly in time-domain. As a numerical example, the derived equation is used to investigate the generalized magneto-thermoelastic behavior of a semi-infinite plate under magnetic field and subjecting to a thermal shock loading. The results demonstrate that FEM can faithfully predict the deformation of the plate and the induced magnetic field, and most importantly can reveal the sophisticated second sound effect of heat conduction in two-dimensional generalized thermoelastic solids, which is usually difficult to model by routine transformation methods. A peak can be observed in the distribution of stress and induced magnetic field at the heat wave front and the magnitude of the peak decreases with time, which can not be obtained by transformation methods. The new method can also be used to study generalized piezo-thermoelastic problems.The project supported by the National Natural Science Foundation of China (10132010 and 10472089)The English text was polished by Yunming Chen     

A compact finite difference method for solving Burgers' equation

In this paper, a high‐order accurate compact finite difference method using the Hopf–Cole transformation is introduced for solving one‐dimensional Burgers' equation numerically. The stability and convergence analyses for the proposed method are given, and this method is shown to be unconditionally stable. To demonstrate efficiency, numerical results obtained by the proposed scheme are compared with the exact solutions and the results obtained by some other methods. The proposed method is second‐ and fourth‐order accurate in time and space, respectively. Copyright © 2009 John Wiley & Sons, Ltd.     

Convergence and exact solutions of spline finite strip method using unitary transformation approach

The spline finite strip method(PSM) is one of the most popular numerical methods for analyzing prismatic structures.Efficacy and convergence of the method have been demonstrated in previous studies by comparing only numerical results with analytical results of some benchmark problems.To date,no exact solutions of the method or its explicit forms of error terms have been derived to show its convergence analytically. As such,in this paper,the mathematical exact solutions of spline finite strips in the plat...     

On the use of some numerical damping methods of spurious oscillations in the case of elastic wave propagation

The use of finite element and finite difference methods of spatial and temporal discretization for solving structural dynamics problems gives rise to purely numerical errors. Among the many numerical methods used to damp out the spurious oscillations occurring in the high frequency domain, it is proposed here to analyse and compare the well-known Bulk Viscosity method, which modifies the stresses calculations and a method recently presented by Tchamwa and Wielgosz, which is based on a modification of an explicit time integration algorithm. The efficiency of both methods is evaluated in a 2-D axisymmetric compressive test.     

Variational principles of elastic-viscous dynamics in laplace transformation form,F. E. M. formulation and numerical method

The author gives variational principles of elastic-viscous dynamics in spectral resolving form^[1], it will be extended to Laplace transformation form in this paper, mixed variational principle of shell dynamics and variational principle of dynamics of elastic-viscous-porous media are concerned, for the latter, F. E. M. formulation has been worked out.Variational principles in Laplace transformation form have concise forms, for the sake of utilizing F. E. M. conveniently it is necessary to find values of preliminary time function at some instants, when values of Laplace transformation at some points are known, but there are no efficient methods till now. In this paper, a numerical method for finding discrete values of preliminary function is presented, from numerical example we see such a method is efficient.By combining both methods stated above, variational principles in Laplace transformation form and numerical method, a quite wide district of solid dynamic problems can be solved by ths aid of digital computers.     

具有非线性表面换热系数半圆柱薄壳在快速冷却过程中的热分析

用实验和数值计算相结合的方法，得到半圆柱壳体快速冷却过程中内外表面的非线性表面换热系数。在此基础上，用有限单元法对半圆柱壳体的热应力和变形进行了分析。在数值计算中，模拟钢的CCT图，计算了奥氏体、珠光体、贝氏体和马氏体的体积百分比，并将热物理性质和力学性能处理为相变体积百分比和温度的函数。所得结果表明，在半圆柱壳体快速冷却过程中的热应力和变形计算中，有必要考虑非线性表面换热系数、相变等非线性效应。     

Active control for a flexible beam with nonlinear hysteresis and time delay

A new active control method is presented to attenuate vibrations of a flexible beam with nonlinear hysteresis and time delay. The nonlinear and hysteretic behavior of the system is illustrated by the Bouc—Wen model. By specific transformation and augmentation of state parameters, we can convert the motion equation of the system with explicit time delay to the standard state space representation without any explicit time delay. Then the instantaneous optimal control method and Runge—Kutta method in fourth-order are applied to the controller design with time delay. Finally, in order to verify the effectivity of the time-delay controller proposed, numerical simulations are implemented. It is indicated by the simulation results that the control performance will deteriorate if neglect the time delay in process of the controller design and proposed time delay controller works well with both small and large time delay problems.     

An orthogonal self-stress matrix for efficient analysis of cyclically symmetric space truss structures via force method

The analysis of structures is normally carried out through displacement method while the force method is considered as an alternative approach for this purpose and used on occasion. The generation of compatibility conditions (the transpose of self-stress matrix) is one of the major and complicated parts of any structural analysis using force method. In this paper, an efficient method is proposed for producing orthogonal self-stress matrix related to space truss structures with cyclic symmetry. This is actually performed by eigen-decomposition of a special matrix having the same null basis as in equilibrium matrix. Then, the advantages of the obtained compatibility conditions are demonstrated with respect to different formulations such as standard force method, eigen force method and integrated force method. Finally, the efficiency of the presented method is comprehensively compared with three well-known numerical methods and tested on a set of practical examples. The results indicate clearly the significant superiority of the proposed approach in terms of both computational time and the accuracy of the results.     

全机绕流Euler方程多重网格分区计算方法

全机三维复杂形状绕流数值求解只能采用分区求解的方法,本文采用可压缩Euler方程有限体积方法以及多重网格分区方法对流场进行分区计算。数值方法采用改进的van Leer迎风型矢通量分裂格式和MUSCL方法,基于有限体积方法和迎风型矢通量分裂方法,建立一套处理子区域内分界面的耦合条件。各个子区域之间采用显式耦合条件,区域内部采用隐式格式和局部时间步长等,以加快收敛速度。计算结果飞机表面压力分布等气动力特性与实验值进行了比较,二者基本吻合。计算结果表明采用分析“V”型多重网格方法,能提高计算效率,加快收敛速度达到接近一个量级。根据全机数值计算结果和可视化结果讨论了流场背风区域旋涡的形成过程。     

Implicit solution of harmonic balance equation system using the LU-SGS method and one-step Jacobi/Gauss-Seidel iteration

ABSTRACT

This paper presents efficient alternative numerical methods for an implicit solution of the harmonic balance equation system for analysing temporal periodic unsteady flows. The proposed method employs approximate factorisation to decouple the common residual term and the time spectral source term of a harmonic balance equation system when it is discretised implicitly. With this approximate factorisation, the complexity of implicit solution of the discrete system is greatly reduced. The common residual term can be dealt with using a lower-upper symmetric-Gauss-Seidel (LU-SGS) method and the time spectral source term is integrated using a Jacobi iteration (JI) or one step Gauss-Seidel (GS) iteration, leading to the LU-SGS/JI method or LU-SGS/GS method. The NASA stage 35 compressor and the 1.5 stage Aachen turbine were used to demonstrate the effectiveness of the proposed methods in stabilising solution and its advantages in comparison with the existing lower-upper symmetric-Gauss-Seidel/block Jacobi (LU-SGS/BJ) method. The LU-SGS/GS method and the LU-SGS/JI method are more robust than the LU-SGS/BJ method in stabilising solution. The LU-SGS/GS method also has faster and tighter convergence and lower memory consumption in comparison with the LU-SGS/BJ method.     

Calculation of coupled problem between temperature and phase transformation during gas quenching in high pressure

The gas quenching is a modern, effective processing technology. On the basis of nonlinear surface heat-transfer coefficient obtained by Cheng during the gas quenching, the coupled problem between temperature and phase transformation during gas quenching in high pressure was simulated by means of finite element method. In the numerical calculation, the thermal physical properties were treated as the functions of temperature and the volume fraction of phase constituents. In order to avoid effectual "oscillation" of the numerical solutions under smaller time step, the Norsette rational approximate method was used.     

Numerical study of flows of complex fluids between eccentric cylinders using transformation functions

In this paper, we investigate fluid flows between eccentric cylinders by means of two stream‐tube analyses. The first method considers a one‐to‐one global transformation function that allows the physical domain to be transformed into a mapped domain, used as computational domain, that involves concentric streamlines. The second approach uses local transformations and domain decomposition techniques to deal with mixed flow regimes. Both formulations are particularly adapted for handling time‐dependent constitutive equations, since particle‐tracking problems are avoided. Mass conservation is verified in both formulations and the relevant numerical procedure can be carried out using simple meshes built on the mapped streamlines. Fluids obeying anelastic and viscoelastic constitutive equations are considered in the calculations. The numerical results are consistent with those in the literature for the flow rates tested. Application of the method to the K‐BKZ memory‐integral constitutive equation highlights significant differences between the model predictions and those provided by more simple rheological models. Copyright © 2002 John Wiley & Sons, Ltd.     

A combined experimental and numerical method for the solution of generalized elasticity problems

This work presents the method for the investigation of three-dimensionally stressed bodies with arbitrary shape which are under the action of an outside system of arbitrary forces. The combined method is based on syntheses of photoelastic experimental methods (other experimental methods may also be used) and digital methods of discrete analysis. Experimental procedures are used for defining superfluous boundary conditions. The boundary-value problem with such boundary conditions is solved by numerical methods. This approach qualitatively changes the very essence of experimental methods and essentially widens their range. It reduces the amount of measurements required and, at the same time, allows one to obtain complete stress fields throughout a body in a short time. In comparison with numerical methods, the combined method increases the accuracy of problem solutions and, at the same time, reduces the time required for complete investigations.     

成形过程数值模拟的张量时空求解策略

成形过程的数值模拟和分析涉及几何、材料、接触摩擦等高度非线性的耦合作用，采用传统的增量算法会导致巨大计算量。本文建议采用张量时间函数的非增量时空算法，在整个时间和空间域上迭代求解。由于采用新的分离变量构思，以张量表达时间函数，可提高问题求解的速度和精度。但对问题的列式和数值求解方法提出了更高的要求。文中讨论其方法与数值实施。