Computation of Viscous Uniform and Shear Flow over A Circular Cylinder by A Finite Element Method

The incompressible viscous uniform and shear flow past a circular cylinder is studied. The two-dimensional Navier-Stokes equations are solved by a finite element method. The governing equations are discretized by a weighted residual method in space. The stable three-step scheme is applied to the momentum equations in the time integration. The numerical model is firstly applied to the computation of the lid-driven cavity flow for its validation. The computed results agree well with the measured data and other numerical results. Then, it is used to simulate the viscous uniform and shear flow over a circular cylinder for Reynolds numbers from 100 to 1000. The transient time interval before the vortex shedding occurs is shortened considerably by introduction of artificial perturbation. The computed Strouhal number, drag and lift coefficients agree well with the experimental data. The computation shows that the finite element model can be successfully applied to the viscous flow problem.

A Finite Element Solution of Wave Forces on Submerged Horizontal Circular Cylinders

In this paper, a numerical model is established for estimating the wave forces on a submerged horizontal circular cylinder. For predicting the wave motion, a set of two dimensional Navier Stokes equations is solved numerically with a finite element method. In order to track the moving non linear wave surface boundary, the Navier Stokes equations are discretized in a moving mesh system. After each computational time step, the mesh is modified according to the changed wave surface boundary. In order to stabilize the numerical procedure, a three step finite element method is applied in the time integration. The water sloshing in a tank and wave propagation over a submerged bar are simulated for the first time to validate the present model. The computational results agree well with the analytical solution and the experimental data. Finally, the model is applied to the simulation of interaction between waves and a submerged horizontal circular cylinder. The effects of the KC number and the cylinder depth on the wave forces are studied.

THE ANALYSIS OF THIN WALLED COMPOSITE LAMINATED HELICOPTER ROTOR WITH HIERARCHICAL WARPING FUNCTIONS AND FINITE ELEMENT METHOD

In the present paper, a series of hierarchical warping functions is developed to analyze the static and dynamic problems of thin walled composite laminated helicopter rotors composed of several layers with single closed cell. This method is the development and extension of the traditional constrained warping theory of thin walled metallic beams, which had been proved very successful since 1940s. The warping distribution along the perimeter of each layer is expanded into a series of successively corrective warping functions with the traditional warping function caused by free torsion or free beading as the first term, and is assumed to be piecewise linear along the thickness direction of layers. The governing equations are derived based upon the variational principle of minimum potential energy for static analysis and Rayleigh Quotient for free vibration analysis. Then the hierarchical finite element method. is introduced to form a,. numerical algorithm. Both static and natural vibration problems of sample box beams axe analyzed with the present method to show the main mechanical behavior of the thin walled composite laminated helicopter rotor.

Finite Element Analysis of Fretting Behavior Between an Engine Cylinder Block and a Main Bearing Cap

The finite element analysis of the fretting behavior between a cylinder block and a main bearing cap is presented. The stresses, relative fretting slip, frettin g friction work parameter W and crack initiation location parameter Gon the fretting contact surface of the cylinder block are obtained and analyzed. It shows that the fretting fatigue problem of the cylinder block can be quantitat ively explained by WorG. The effects of pretightening force, friction factor and material combination of the cylinder block and the main bearing cap are studied. The computational results indicate that the fretting fatigue of the cylinder block can be allayed by increasing the elastic modulus of the cylinder block, but not by changing the other two factors.

Seismic analysis of cantilever earth retaining walls embedded in dry sand by simplified approaches and finite element method

In engineering practice simplified methods are essential to the seismic design of embedded earth retaining walls,as fullydynamic numerical analyses are costly,time-consuming and require specific expertise.Recently developed pseudostatic methods provide earth stresses and internal forces,even in those cases in which the strength of the soil surrounding the structure is not entirely mobilised.Semiempirical correlations or Newmark sliding block method provide an estimate of earthquake-induced permanent displacements.However,the use of these methods is hindered by uncertainties in the evaluation of a few input parameters,affecting the reliability of the methods.This study uses 1 D site response analyses and 2 D fully-dynamic finite element analyses to show that simplified methods can provide a reasonable estimate of the maximum bending moment and permanent displacements for stiff cantilever walls embedded in uniform sand,providing that a few input parameters are evaluated through semiempirical correlations and a simple 1 D site response analysis.

A Finite Element Analysis on MHD Free Convection Flow in Open Square Cavity Containing Heated Circular Cylinder

The problem of Magnetohydrodynamic (MHD) free convection heat transfer in a square open cavity containing a heated circular cylinder at the centre has been investigated in this work. As boundary conditions of the cavity, the left vertical wall is kept at a constant heat flux, bottom and top walls are kept at different high and low temperature respectively. The remaining side wall is open. Finite element analysis based on Galerkin weighted Residual approach is used to visualize the temperature distribution and fluid flow solving two-dimensional governing mass, momentum and energy equations for steady state, natural convection flow in presence of magnetic field in side an open square cavity. A uniformly heated circular cylinder is located at the centre of the cavity. The object of this study is to describe the effects of MHD on the thermal fields and flow in presence of such heated circular cylinder by visualization of graph. The investigations are conducted for different values of Rayleigh number (Ra) and Hartmann number (Ha). The results show that the temperature field and flow pattern are significantly dependent on the above mentioned parameters.

Impact of a Magnetic Dipole on Heat Transfer in Non-Conducting Magnetic Fluid Flow over a Stretching Cylinder

The thermal behavior of an electrically non-conducting magnetic liquid flowing over a stretching cylinder under the influence of a magnetic dipole is considered.The governing nonlinear differential equations are solved numerically using a finite element approach,which is properly validated through comparison with earlier results available in the literature.The results for the velocity and temperature fields are provided for different values of the Reynolds number,ferromagnetic response number,Prandtl number,and viscous dissipation parameter.The influence of some physical parameters on skin friction and heat transfer on the walls of the cylinder is also investigated.The applicability of this research to heat control in electronic devices is discussed to a certain extent.

Least Square Finite Element Method for Viscous Splitting of Unsteady Incompressible Navier–Stokes Equations

In order to solve unsteady incompressible Navier–Stokes(N–S) equations, a new stabilized finite element method,called the viscous-splitting least square FEM, is proposed. In the model, the N–S equations are split into diffusive and convective parts in each time step. The diffusive part is discretized by the backward difference method in time and discretized by the standard Galerkin method in space. The convective part is a first-order nonlinear equation.After the linearization of the nonlinear part by Newton’s method, the convective part is also discretized by the backward difference method in time and discretized by least square scheme in space. C0-type element can be used for interpolation of the velocity and pressure in the present model. Driven cavity flow and flow past a circular cylinder are conducted to validate the present model. Numerical results agree with previous numerical results, and the model has high accuracy and can be used to simulate problems with complex geometry.

Numerical simulation of flow past twin near-wall circular cylinders in tandem arrangement at low Reynolds number

Fluid flow past twin circular cylinders in a tandem arrangement placed near a plane wall was investigated by means of numerical simulations. The two-dimensional Navier-Stokes equations were solved with a three-step finite element method at a relatively low Reynolds number of Re -- 200 for various dimensionless ratios of 0.25 ≤ G/D ≤2.0 and 1.0 ≤ L/D ≤ 4.0, where D is the cylinder diameter, L is the center-to-center distance between the two cylinders, and G is the gap between the lowest surface of the twin cylinders and the plane wall. The influences of G/D and L/D on the hydrodynamic force coefficients, Strouhal numbers, and vortex shedding modes were examined. Three different vortex shedding modes of the near wake were identified according to the numerical results. It was found that the hydrodynamic force coefficients and vortex shedding modes are quite different with respect to various combinations of G/D and L/D. For very small values of G/D, the vortex shedding is completely suppressed, resulting in the root mean square （RMS） values of drag and lift coefficients of both cylinders and the Strouhal number for the downstream cylinder being almost zero. The mean drag coefficient of the upstream cylinder is larger than that of the downstream cylinder for the same combination of G/D and L/D. It is also observed that change in the vortex shedding modes leads to a significant increase in the RMS values of drag and lift coefficients.

Thermal Analysis of Turbine Blades with Thermal Barrier Coatings Using Virtual Wall Thickness Method

Avirtual wall thicknessmethod is developed to simulate the temperature field of turbine bladeswith thermal barrier coatings(TBCs),to simplify the modeling process and improve the calculation efficiency.The results show that the virtualwall thickness method can improve themesh quality by 20%,reduce the number ofmeshes by 76.7%and save the calculation time by 35.5%,compared with the traditional real wall thickness method.The average calculation error of the two methods is between 0.21%and 0.93%.Furthermore,the temperature at the blade leading edge is the highest and the average temperature of the blade pressure surface is higher than that of the suction surface under a certain service condition.The blade surface temperature presents a high temperature at both ends and a low temperature in themiddle height when the temperature of incoming gas is uniformand constant.The thermal insulation effect of TBCs is the worst near the air film hole,and the best at the blade leading edge.According to the calculated temperature field of the substrate-coating system,the highest thermal insulation temperature of the TC layer is 172.01 K,and the thermal insulation proportions of TC,TGO and BC are 93.55%,1.54%and 4.91%,respectively.

Fully Nonlinear Simulations of Wave Resonance by An Array of Cylinders in Vertical Motions

The finite element method （FEM） is employed to analyze the resonant oscillations of the liquid confined within multiple or an array of floating bodies with fully nonlinear boundary conditions on the free surface and the body surface in two dimensions. The velocity potentials at each time step are obtained through the FEM with 8-node quadratic shape functions. The finite element linear system is solved by the conjugate gradient （CG） method with a symmetric successive overelaxlation （SSOR） preconditioner. The waves at the open boundary are absorbed by the combination of the damping zone method and the Sommerfeld-Orlanski equation. Numerical examples are given by an array of floating wedge- shaped cylinders and rectangular cylinders. Results are provided for heave motions including wave elevations, profiles and hydrodynamic forces. Comparisons are made in several cases with the results obtained from the second order solution in the time domain. It is found that the wave amplitude in the middle region of the array is larger than those in other places, and the hydrodynamic force on a cylinder increases with the cylinder closing to the middle of the array.

Short-Crested Waves Interaction with A Concentric Porous Cylinder System with Partially Porous Outer Cylinder

In this paper, based on the linear wave theory, the interaction of short-crested waves with a concentric dual cylindrical system with a partially porous outer cylinder is studied by using the scaled boundary finite element method （SBFEM）, which is a novel semi-analytical method with the advantages of combining the finite element method （FEM） with the boundary element method （BEM）. The whole solution domain is divided into one unbounded sub-domain and one bounded sub-domain by the exterior cylinder. By weakening the governing differential equation in the circumferential direction, the SBFEM equations for both domains can be solved analytically in the radial direction. Only the boundary on the circumference of the exterior porous cylinder is discretized with curved surface finite elements. Meanwhile, by introducing a variable porous-effect parameter G, non-homogeneous materials caused by the complex configuration of the exterior cylinder are modeled without additional efforts. Comparisons clearly demonstrate the excellent accuracy and computational efficiency associated with the present SBFEM. The effects of the wide range wave parameters and the structure configuration are examined. This parametric study will help determine the various hydrodynamic effects of the concentric porous cylindrical structure.

Initial residual stress experiment and simulation of thin-walled parts for layer removal method

Thin-walled parts have low stiffness characteristic. Initial residual stress of thin-walled blanks is an important influence factor on machining stability. The present work is to verify the feasibility of an initial residual stress measurement of layer removal method. According to initial residual stress experiment for casting ZL205 A aluminum alloy tapered thin-walled blank by a common method,namely hole-drilling method,three finite element models with initial residual stress are established to simulate the layer removal method in ABAQUS and ANSYS software. By analyzing the results of simulation and experiments,the cutting residual stress inlayer removal process has a significant effect on measurement results. Reducing cutting residual stress is helpful to improve accuracy of layer removal method.

Powell's optimal identification of material constants of thin-walled box girders based on Fibonacci series search method

A dynamic Bayesian error function of material constants of the structure is developed for thin-walled curve box girders. Combined with the automatic search scheme with an optimal step length for the one-dimensional Fibonacci series, Powell＇s optimization theory is used to perform the stochastic identification of material constants of the thin-walled curve box. Then, the steps in the parameter identification are presented. Powell＇s identification procedure for material constants of the thin-walled curve box is compiled, in which the mechanical analysis of the thin-walled curve box is completed based on the finite curve strip element （FCSE） method. Some classical examples show that Powell＇s identification is numerically stable and convergent, indicating that the present method and the compiled procedure are correct and reliable. During the parameter iterative processes, Powell＇s theory is irrelevant with the calculation of the FCSE partial differentiation, which proves the high computation efficiency of the studied methods. The stochastic performances of the system parameters and responses axe simultaneously considered in the dynamic Bayesian error function. The one-dimensional optimization problem of the optimal step length is solved by adopting the Fibonacci series search method without the need of determining the region, in which the optimized step length lies.

模块化防屈曲钢板墙抗震性能

考虑经济性和便利性,提出了模块化防屈曲钢板墙及其常用模数,然后基于已有的防屈曲钢板墙的简化计算模型,推导了在不同墙数量情况下,模块化防屈曲钢板墙充分发挥抗震性能时,周边梁的承载力和刚度需求,并提出了周边梁的设计方法.进而,对模块化防屈曲钢板墙和整体墙进行了对比验证试验及其有限元模拟分析.研究结果表明:模块化防屈曲钢板墙能够达到与整体钢板墙相同的抗震性能,且周边梁未发生明显转动和破坏,说明所提出的周边梁设计方法是合理的;在层间位移角为1/50时,防屈曲钢板墙上下周边梁均未进入塑性,且承载力和初始刚度等的有限元值与试验值误差小于6%,从而进一步验证了其抗震性能.

某在建商业建筑满足承载力要求的消能减震加固设计与分析

本文以新疆某商业楼加固工程为背景,针对该建筑工程施工时一、二层梁柱纵向钢筋等级选择与设计要求不符,造成多数框架梁、柱配筋不满足抗震承载力要求的问题,提出了传统加固技术与消能减震技术相结合的加固方案。基于PKPM和ETABS软件,通过大量数据分析对加固方案进行了优化和调整,同时采用PERFORM-3D软件分析了罕遇地震作用下结构的时程响应。得到结论:(1)减震加固后,结构满足“小震不坏、大震不倒”的设防目标;(2)通过合理布置黏滞阻尼墙,降低了地震作用,大幅度减小了传统加固的用量。

FINITE ELEMENT ANALYSIS OF LAMINAR FLOW AND HEAT TRANSFER IN A BUNDLE OF CYLINDERS

Two-dimensional Navier-Stokes equations and energy equation governing incompressible laminar flow past a bundle of cylinders were numerically solved by using the finite element method. The velocity correction method was used for time advancement, and spatial discretization was carried out with the Galerkin weighted residual method. Viscous flows past the cylinder banks arranged in in-line cylinder bundles and staggered cylinder bundles, coupled with heat transfer, were investigated for pitch-diameter ratios of 1.5 and 2.0 and the Reynolds numbers from 50 to 500. Flow structures and heat transfer behavior were discussed. The results obtained agree well with available numerical data.

Finite element analysis of second order wave radiation by a group of cylinders in the time domain

A finite element based numerical method is employed to analyze the wave radiation by multiple or a group of cylinders in the time domain. The nonlinear free surface and body surface boundary conditions are satisfied based on the perturbation method up to the second order. The first- and second-order velocity potential problems at each time step are solved through a Finite Element Method （FEM）. The matrix equation of the FEM is solved through iteration and the initial solution is obtained from the result at the previous time step. The three-dimensional （3-D） mesh required is generated based on a two-dimensional （2-D） hybrid mesh on a horizontal plane and its extension in the vertical direction. The hybrid mesh is generated by combining an unstructured grid away from cylinders and two structured grids near the cylinder and the artificial boundary. The fluid velocity on the free surface and the cylinder surface are calculated by using a differential method. Results for various configurations including the cases of two cylinders and four cylinders and a group of eighteen cylinders are obtained to show the joint influences of cylinders on the first- and second- order waves and forces, including the effects of spacing ratios and wave frequency on the second order waves and the mean force, in particular.

仿香蒲结构的侧向抗冲击性能研究

自然界中的香蒲茎秆能较好地抵御外部和自身荷载性能,其微观结构由大量直径不相同的细胞孔构成,基于其特点提出横截面的多尺度模型,并设计一种具有良好侧向抗冲击性能的仿香蒲结构。使用有限元方法研究仿香蒲结构的侧向抗冲击性能,对其相关参数进行分析,建立同质量的仿竹子、仿双节竹、仿仙人掌、仿马尾草、仿蜘蛛网和仿螳螂虾共6种仿生模型,并与仿香蒲结构模型进行对比。研究结果表明:横向胞壁间距、胞孔壁厚度、表皮厚度、胞孔直径这4个参数对仿香蒲结构的性能有显著影响,在较小的横向胞壁间距、中等的胞孔壁厚度、较大的表皮厚度以及由外到内先变大后变小的胞孔直径下能得到较高的比吸能;新型仿香蒲结构的比吸能和冲击力均明显比其他仿生结构的高,说明仿香蒲结构具有较强的吸能能力和较高的吸能效率。

单臂圆筒旋转超声电机优化设计

为了解决现有圆筒型超声电机中存在的定子结构复杂及振动耦合问题,设计了一种单臂圆筒旋转超声电机。所设计电机仅使用一组压电陶瓷片即可直接实现圆筒中两列驻波的激励,而非借助纵振陶瓷组和弯振陶瓷组的配合。介绍了单臂圆筒旋转超声电机的工作原理,确定了电机定子的结构参数,并基于有限元对其进行了仿真,在驱动齿所选质点处得到了椭圆的振动轨迹。仿真结果表明该优化设计能够简化电机结构,避免振型干涉,提高电机的机械输出能力。