Geometric nonlinear dynamic analysis of curved beams using curved beam element

Instead of using the previous straight beam element to approximate the curved beam,in this paper,a curvilinear coordinate is employed to describe the deformations,and a new curved beam element is proposed to model the curved beam.Based on exact nonlinear strain-displacement relation,virtual work principle is used to derive dynamic equations for a rotating curved beam,with the effects of axial extensibility,shear deformation and rotary inertia taken into account.The constant matrices are solved numerically utilizing the Gauss quadrature integration method.Newmark and Newton-Raphson iteration methods are adopted to solve the differential equations of the rigid-flexible coupling system.The present results are compared with those obtained by commercial programs to validate the present finite method.In order to further illustrate the convergence and efficiency characteristics of the present modeling and computation formulation,comparison of the results of the present formulation with those of the ADAMS software are made.Furthermore,the present results obtained from linear formulation are compared with those from nonlinear formulation,and the special dynamic characteristics of the curved beam are concluded by comparison with those of the straight beam.     

Determining local values of the heat transfer coefficient on a fin surface

The paper discusses methods of establishing the course of a boiling curve that base on measurements of temperature distribution on a fin surface. Two cases of a local correlation between the heat transfer coefficient and the wall superheat are considered. The first is the case of a smooth copper fin for which the boiling curve is defined using the power-law formula. It is essential to determine, however, how sensitive this method is to measurement errors. In the other case, the local values of the heat transfer coefficient are approximated by means of a power-law polynomial. The method efficiency is studied using the example of a fin surface modified with a laser beam.     

Formulation of thermo-hydro-mechanical coupling behavior of unsaturated soils based on hybrid mixture theory

Thermo-Hydro-Mechanical （THM） coupling pro- cesses in unsaturated soils are very important in both theoretical researches and engineering applications. A coupled formulation based on hybrid mixture theory is derived to model the THM coupling behavior of unsaturated soils. The free-energy and dissipative functions for different phases are derived from Taylor＇s series expansions. Constitutive relations for THM coupled behaviors of unsaturated soils, which include deformation, entropy change, fluid flow, heat conduction, and dynamic compatibility conditions on the interfaces, are then established. The number of field equations is shown to be equal to the number of unknown variables; thus, a closure of this coupling problem is established. In addition to modifications of the physical conservation equations with coupling effect terms, the constitutive equations, which consider the coupling between elastoplastic deformation of the soil skeleton, fluid flow, and heat transfer, are also derived.     

Heat transfer from extended surfaces subject to variable heat transfer coefficient

The present article investigates the effect of locally variable heat transfer coefficient on the performance of extended surfaces (fins) subject to natural convection. Fins of different profiles have been investigated. The fin profiles presently considered are namely; straight and pin fin with rectangular (constant diameter), convex parabolic, triangular (conical) and concave parabolic profiles and radial fins with constant profile with different radius ratios. The local heat transfer coefficient was considered as function of the local temperature and has been obtained using the available correlations of natural convection for each pertinent extended surface considered. The performance of the fin has been expressed in terms of the fin efficiency. Comparisons between the present results for all fins considered and the results obtained for the corresponding fins subject to constant heat transfer coefficient along the fin are presented. Comparisons, i.e. showed an excellent agreement with the experimental results available in the literature. Results show that there is a considerable deviation between the fin efficiency calculated based on constant heat transfer coefficient and that calculated based on variable heat transfer coefficient and this deviation increases with the dimensionless parameter m.     

Data reduction for air-side performance of fin-and-tube heat exchangers

The present study focuses on the data reduction method to obtain the air-side performance of fin-and-tube heat exchangers. The data reduction methodology for air-side heat transfer coefficients in the literature is not based on a consistent approach. This paper recommends standard procedures for dry surface heat transfer in finned-tube heat exchangers having water on the tube-side. Inconsistencies addressed include the -NTU relationships, calculation of the tube-side heat transfer coefficient, calculation of fin efficiency, and whether entrance and exit loss should be included in the reduction of friction factors. Use of the recommended standardized methodology will provide more meaningful data for use in the development of correlations, or for performance comparison purposes.     

Effects of thermophoresis particle deposition and of the thermal conductivity in a porous plate with dissipative heat and mass transfer

Network simulation method(NSM) is used to solve the laminar heat and mass transfer of an electricallyconducting,heat generating/absorbing fluid past a perforated horizontal surface in the presence of viscous and Joule heating problem. The governing partial differential equations are non-dimensionalized and transformed into a system of nonlinear ordinary differential similarity equations,in a single independent variable,η. The resulting coupled,nonlinear equations are solved under appropriate transformed boundary conditions. Computations are performed for a wide range of the governing flow parameters,viz Prandtl number,thermophoretic coeffcient(a function of Knudsen number),thermal conductivity parameter,wall transpiration parameter and Schmidt number. The numerical details are discussed with relevant applications. The present problem finds applications in optical fiber fabrication,aerosol filter precipitators,particle deposition on hydronautical blades,semiconductor wafer design,thermo-electronics and problems including nuclear reactor safety.     

Contact-impact formulation for multi-body systems using component mode synthesis

The efficiency and accuracy are two most concerned issues in the modeling and simulation of multi-body systems involving contact and impact. This paper proposed a formulation based on the component mode synthesis method for planar contact problems of flexible multi-body systems. A flexible body is divided into two parts： a contact zone and an un-contact zone. For the un-contact zone, by using the fixed-interface substructure method as reference, a few low-order modal coordinates are used to replace the nodal coordinates of the nodes, and meanwhile, the nodal coordinates of the local impact region are kept unchanged, therefore the total degrees of freedom （DOFs） are greatly cut down and the computational cost of the simulation is significantly reduced. By using additional constraint method, the impact constraint equations and kinematic constraint equations are derived, and the Lagrange equations of the first kind of flexible multi-body system are obtained. The impact of an elastic beam with a fixed half disk is simulated to verify the efficiency and accuracy of this method.     

Natural convection heat transfer from a thin rectangular fin with a line source at the base — a finite difference solution

This paper presents an analysis of the problem of a thin fin of finite thermal conductivity, with an isothermal line source at the base, dissipating heat to the surrounding air by natural convection. The horizontal surface to which the fin is attached is adiabatic so that heat is dissipated only through the fin. The temperature and velocity distributions in the field, the temperature profile in the fin, local Nusselt numbers along the fin and the average heat transfer coefficient of the fin are obtained by solving the governing equations in the field and the heat transfer equation in the fin simultaneously, using an explicit unsteady Finite Difference formulation leading to the steady state result. Numerical experiments are performed to study the influence of parameters namely the fin height, temperature of the heating source and the fin material on the average heat transfer coefficient. Comparison is made with fins of infinite thermal conductivity and the vertical isothermal flat plate.     

Plane Deformations in Incompressible Nonlinear Elasticity

The substantially general class of plane deformation fields, whose only restriction requires that the angular deformation not vary radially, is considered in the context of isotropic incompressible nonlinear elasticity. Analysis to determine the types of deformations possible, that is, solutions of the governing systems of nonlinear partial differential equations and constraint of incompressibility, is developed in general. The Mooney-Rivlin material model is then considered as an example and all possible solutions to the equations of equilibrium are determined. One of these is interpreted in the context of nonradially symmetric cavitation, i.e., deformation of an intact cylinder to one with a double-cylindrical cavity. Results for general incompressible hyperelastic materials are then discussed. The novel approach taken here requires the derivation and use of a material formulation of the governing equations; the traditional approach employing a spatial formulation in which the governing equations hold on an unknown region of space is not conducive to the study of deformation fields containing more than one independent variable. The derivation of the cylindrical polar coordinate form of the equilibrium equations for the nominal stress tensor (material formulation) for a general hyperelastic solid and a fully arbitrary cylindrical deformation field is also given. This revised version was published online in August 2006 with corrections to the Cover Date.     

Effect of variable ambient temperature on fin efficiency in a two-dimensional flow passage

The fin efficiency in a heat exchanger element that is a simplification of one row in a tube-and-fin heat exchanger was theoretically examined within wide ranges of the affecting variables: the conventional fin efficiency and the isothermal effectiveness of the heat exchanger. These variables are suggested for use also in the further studies. An analytical solution can be found for the case of a constant heat transfer coefficient. The ambient temperature variation alone decreases the fin efficiency less than 4%. The local heat transfer coefficient obtained from the numerical fluid flow simulations is strongly affected by the fin properties because the thermal boundary conditions for the fluid flow changes. On a poorly conducting fin surface the heat transfer coefficient in front of the fin base is much larger than on an isothermal fin because the heat flux is increasing in the flow direction. At low fin efficiencies this compensates for the decrease in fin efficiency due to ambient temperature variation.     

Effect of induced magnetic field on natural convection in vertical concentric annuli

In the present paper,we have considered the steady fully developed laminar natural convective flow in open ended vertical concentric annuli in the presence of a radial magnetic field.The induced magnetic field produced by the motion of an electrically conducting fluid is taken into account.The transport equations concerned with the considered model are first recast in the non-dimensional form and then unified analytical solutions for the velocity,induced magnetic field and temperature field are obtained for the cases of isothermal and constant heat flux on the inner cylinder of concentric annuli.The effects of the various physical parameters appearing into the model are demonstrated through graphs and tables.It is found that the magnitude of maximum value of the fluid velocity as well as induced magnetic field is greater in the case of isothermal condition compared with the constant heat flux case when the gap between the cylinders is less or equal to 1.70 times the radius of inner cylinder,while reverse trend occurs when the gap between the cylinders is greater than 1.71 times the radius of inner cylinder.These fields are almost the same when the gap between the cylinders is equal to 1.71 times the radius of inner cylinder for both the cases.It is also found that as the Hartmann number increases,there is a flattening tendency for both the velocity and the induced magnetic field.The influence of the induced magnetic field is to increase the velocity profiles.     

Geometric nonlinear formulation for thermal-rigid-flexible coupling system

This paper develops geometric nonlinear hybrid formulation for flexible multibody system with large deformation considering thermal efect. Diferent from the conventional formulation, the heat flux is the function of the rotational angle and the elastic deformation, therefore, the coupling among the temperature, the large overall motion and the elastic deformation should be taken into account. Firstly,based on nonlinear strain–displacement relationship, variational dynamic equations and heat conduction equations for a flexible beam are derived by using virtual work approach,and then, Lagrange dynamics equations and heat conduction equations of the first kind of the flexible multibody system are obtained by leading into the vectors of Lagrange multiplier associated with kinematic and temperature constraint equations. This formulation is used to simulate the thermal included hub-beam system. Comparison of the response between the coupled system and the uncoupled system has revealed the thermal chattering phenomenon. Then, the key parameters for stability, including the moment of inertia of the central body, the incident angle, the damping ratio and the response time ratio, are analyzed. This formulation is also used to simulate a three-link system applied with heat flux. Comparison of the results obtained by the proposed formulation with those obtained by the approximate nonlinear model and the linear model shows the significance of considering all the nonlinear terms in the strain in case of large deformation. At last, applicability of the approximate nonlinear model and the linear model are clarified in detail.     

Boiling of HFE-7100 on straight pin fin

This paper, for the first time, experimentally elucidated pin fin boiling characteristics of saturated and subcooled HFE-7100, a CFC-substitute, under atmospheric pressure. Fin base temperature and heat flux data are, for the first time, measured along with the fin tip temperature. For a given liquid/fin combination there exist upper steady-state (USS) branch, lower steady-state (LSS) branch, and a large, unstable regime located in between the two branches. Zones with different stability characteristics are mapped according to boiling on fins with different aspect ratios. Liquid subcooling enhances heat transfer performance and boiling stability. Calculation of boiling curve qualitatively describes the experimental results.     

A velocity transformation method for the nonlinear dynamic simulation of railroad vehicle systems

The solution of the constrained multibody system equations of motion using the generalized coordinate partitioning method requires the identification of the dependent and independent coordinates. Using this approach, only the independent accelerations are integrated forward in time in order to determine the independent coordinates and velocities. Dependent coordinates are determined by solving the nonlinear constraint equations at the position level. If the constraint equations are highly nonlinear, numerical difficulties can be encountered or more Newton–Raphson iterations may be required in order to achieve convergence for the dependent variables. In this paper, a velocity transformation method is proposed for railroad vehicle systems in order to deal with the nonlinearity of the constraint equations when the vehicles negotiate curved tracks. In this formulation, two different sets of coordinates are simultaneously used. The first set is the absolute Cartesian coordinates which are widely used in general multibody system computer formulations. These coordinates lead to a simple form of the equations of motion which has a sparse matrix structure. The second set is the trajectory coordinates which are widely used in specialized railroad vehicle system formulations. The trajectory coordinates can be used to obtain simple formulations of the specified motion trajectory constraint equations in the case of railroad vehicle systems. While the equations of motion are formulated in terms of the absolute Cartesian coordinates, the trajectory accelerations are the ones which are integrated forward in time. The problems associated with the higher degree of differentiability required when the trajectory coordinates are used are discussed. Numerical examples are presented in order to examine the performance of the hybrid coordinate formulation proposed in this paper in the analysis of multibody railroad vehicle systems.     

Viscous falling film instability around a vertical moving cylinder

The present work discusses both the linear and nonlinear stability conditions of a viscous falling film down the outer surface of a solid vertical cylinder which moves in the direction of its axis with a constant velocity.After studying the linear conditions,a generalized nonlinear kinematic model is then derived to present the physical system.Applying the boundary conditions,analytical solutions are obtained using the long-wave perturbation method.In the first step,the normal mode method is used to characterize the linear behaviors.In the second step,the nonlinear film flow model is solved by using the method of multiple scales,to obtain Ginzburg-Landau equation.The influence of some physical parameters is discussed in both linear and nonlinear steps of the problem,and the results are displayed in many plots showing the stability criteria in various parameter planes.     

Numerical analysis of laminar natural convection in enclosures with fins attached to an active wall

Vertical enclosures with conducting fins attached to the cold wall were considered. Side walls were kept at constant but different temperatures, while horizontal top and bottom walls were insulated. A conjugate formulation was used for the mathematical formulation of the problem, and a computer program based on the control volume approach and the SIMPLE algorithm was developed. Computations were performed to investigate the effects of the fin configuration and Rayleigh number on the flow structure and heat transfer. It was observed that the heat transfer rate through an enclosure can be controlled by attaching fins to the wall(s) of the enclosure. At low Rayleigh numbers (conduction controlled regime), the heat transfer rate increases with the increasing number of fins and the fin length. However, at higher Rayleigh numbers (convection dominant regimes), the heat transfer rate can be decreased or increased by properly choosing the number of fins and the fin lengths. Received on 07 April 1997     

On the limitations of linear beams for the problems of moving mass-beam interaction using a meshfree method

This paper deals with the capabilities of linear and nonlinear beam theories in predicting the dynamic response of an elastically supported thin beam traversed by a moving mass. To this end, the discrete equations of motion are developed based on Lagrange’s equations via reproducing kernel particle method (RKPM). For a particular case of a simply supported beam, Galerkin method is also employed to verify the results obtained by RKPM, and a reasonably good agreement is achieved. Variations of the maximum dynamic deflection and bending moment associated with the linear and nonlinear beam theories are investigated in terms of moving mass weight and velocity for various beam boundary conditions. It is demonstrated that for majority of the moving mass velocities, the differences between the results of linear and nonlinear analyses become remarkable as the moving mass weight increases, particularly for high levels of moving mass velocity. Except for the cantilever beam, the nonlinear beam theory predicts higher possibility of moving mass separation from the base beam compared to the linear one. Furthermore, the accuracy levels of the linear beam theory are determined for thin beams under large deflections and small rotations as a function of moving mass weight and velocity in various boundary conditions.     

Research of air-cushion isolation effects on high arch dam reservoir

A three-dimensional (3D) finite element model of air-cushion isolated arch dam is presented with the nonlinear gas-liquid-solid multi-field dynamic coupling effect taken into account.In this model,the displacement formulation in Lagrange method,pressure formulation in Euler method,nonlinear contact model based on Coulomb friction law are applied to the air-cushion,reservoir and contraction joint domain,respectively.The dynamic response of Jinping I arch dam with a height of 305 m is analyzed using the seismic records of the Wenchuan Earthquake in 2008.Numerical results show that the air-cushion isolation reduces significantly the hydrodynamic pressure as well as the opening width for the contraction joints of high arch dam.     

Analysis of axial strain in one-dimensional loading by different models

Different phenomenological equations based on plasticity, primary creep （as a viscoplastic mechanism）, secondary creep （as another viscoplastic mechanism） and different combinations of these equations are presented and used to describe the material inelastic deformation in uniaxial test. Agreement of the models with experimental results and with the theoretical concepts and physical realities is the criterion of choosing the most appropriate formulation for uniaxial test. A model is thus proposed in which plastic deformation, primary creep and secondary creep contribute to the inelastic deformation. However, it is believed that the hardening parameter is composed of plastic and primary creep parts. Accordingly, the axial plastic strain in a uniaxial test may no longer be considered as the hardening parameter. Therefore, a proportionality concept is proposed to calculate the plastic contribution of deformation.