Numerical analysis of shear-band bifurcation

A numerical analysis of bifurcation in shear band pattern is presented to help understand the distributions of the velocity variation in the shear band. Comparison with the results of analytic method indicates that: (1) the critical strain is irrelevant to the relative width of the shear band; (2) the variations along the direction normal to the band have indeed the controlling effect whilst the effect of variations along the tangential direction is negligible. The project is supported by the National Natural Science Foundation of China.     

Strain analysis for extrusion of powder metals

A method is described for determining the strain distribution during axisymmetric extrusion of powder metals. Density variations in these materials are included in the analysis. The method examines these changes together with the laminar flow pattern to determine strain rates and, ultimately, strains.     

Design Sensitivity Analysis of Truss Structures with Elastoplastic Material*

ABSTRACT

A continuum-based design sensitivity analysis (DSA) method is presented for configuration (or layout) design of nonlinear structural systems with rate-independent elastoplastic material. Configuration design variables are characterized by shape and orientation changes of the structural component. A continuum-based shape DSA method that utilizes the material derivative of continuum mechanics is extended to account for effects of shape and orientation variations. The incremental analysis method, with updated Lagrangian formulation, is used to derive the design sensitivity for the nonlinear structural system.

To derive the design sensitivity, incremental energy and load forms are utilized. The first variations of energy and load forms and the static response with respect to configuration design variables are described using the material derivative. Direct differentiation is utilized to obtain the first variation of the performance measure explicitly in terms of variations of configuration design variables. With the consistent tangent stiffness matrix employed at the end of each load step to compute the design sensitivity, it is found that no iterations are necessary to compute design sensitivity. In addition, the linear design velocity is used to account for configuration design changes, with the velocity field being updated at each load step of the incremental analysis.     

“Non-linear normal modes” and the generalized Ritz method in the problems of vibrations of non-linear elastic continuous systems

In the study of natural vibrations of non-linear elastic systems it is shown that the mode shape of the vibration can vary with the amplitude as well as the frequency, and that the amplitude frequency relation is strongly affected by constraints imposed on the mode shape in an approximate solution. A method is developed which assumes the approximate solution in the form of a truncated series in which, instead of the set of coefficients, the set of functions of spatial variables is unknown and then determined by a procedure that can be regarded as a generalization of the Ritz method. The problem of variations of the normal mode shapes and of the associated natural frequencies with the amplitude is illustrated by two examples of beams with non-linear boundary conditions, and the amplitude-frequency relation is compared to that corresponding to the a priori assumed linear normal mode solution. Further possible consequences of the mode shape amplitude variations in forced, resonant motion of nonconservative systems are also indicated.     

环境温度影响下基于支持向量机与强化飞蛾扑火优化算法的结构稀疏损伤识别

结构处于自然环境中常会受到环境温度变化的影响,引起实测动力响应出现较大误差,进一步影响对结构健康状况的判定.另外,基于优化算法的损伤识别在反演损伤位置及量化损伤程度时,易出现局部最优解,且计算效率低下.针对以上难题,本文提出一种结合支持向量机与强化飞蛾扑火优化算法的损伤识别方法,用于对环境温度影响下的结构稀疏损伤进行识...     

磁头/磁盘气体润滑特性计算与分析

提出了1种计算超薄气膜润滑轴承压力分布的有限差分法,在此基础上对某磁头快速或缓慢偏离平衡位置造成的压力变化及气体轴承刚度进行分析.结果表明:随着气膜特征高度降低,初始纵翻倾角减小,气体轴承刚度增大;垂直于磁盘方向的微小位移不会对磁头平衡造成威胁,但磁头纵翻造成的翻转力矩变化较复杂,需要加以控制;磁头飞行姿态的突变将引起明显的挤压效应,挤压效应的强弱与初始飞行姿态有关.     

Optimal design of spiral casing tongue and wicket gate angle by decomposition method

A method of optimal design of Francis turbine tongue and wicket gate angle for given spiral casing is proposed. The potential flow in the doubly connected domain is decomposed into basic and circulation flows. The intensity of circulation is then calculated by the least-squares method minimizing the error function equal to the sum of squares of differences between given and calculated circumferential velocities in the outflow boundary nodes. It is shown that the error function has a sharp minimum, which qualifies the proposed method as well defined. For given numerical example, the variations in the outflow velocity angles are much smaller for optimal than for already used non-optimal design. A finite element method is used, with originally developed pre- and post-processor and frontal solver suited for personal computers.     

A generalized self-consistent method accounting for fiber section shape

A three-phase confocal elliptical cylinder model is proposed for fiber-reinforced composites, in terms of which a generalized self-consistent method is developed for fiber-reinforced composites accounting for variations in fiber section shapes and randomness in fiber section orientation. The reasonableness of the fiber distribution function in the present model is shown. The dilute, self-consistent, differential and Mori–Tanaka methods are also extended to consider randomness in fiber section orientation in a statistical sense. A full comparison is made between various micromechanics methods and with the Hashin and Shtrikman’s bounds. The present method provides convergent and reasonable results for a full range of variations in fiber section shapes (from circular fibers to ribbons), for a complete spectrum of the fiber volume fraction (from 0 to 1, and the latter limit shows the correct asymptotic behavior in the fully packed case) and for extreme types of the inclusion phases (from voids to rigid inclusions). A very different dependence of the five effective moduli on fiber section shapes is theoretically predicted, and it provides a reasonable explanation on the poor correlation between previous theory and experiment in the case of longitudinal shear modulus.     

Large-eddy simulation of turbulent channel flow at transcritical states

We present well-resolved large-eddy simulations (LES) of a channel flow solving the fully compressible Navier–Stokes equations in conservative form. An adaptive look-up table method is used for thermodynamic and transport properties. A physically consistent subgrid-scale turbulence model is incorporated, that is based on the Adaptive Local Deconvolution Method (ALDM) for implicit LES. The wall temperatures are set to enclose the pseudo-boiling temperature at a supercritical pressure, leading to strong property variations within the channel geometry. The hot wall at the top and the cold wall at the bottom produce asymmetric mean velocity and temperature profiles which result in different momentum and thermal boundary layer thicknesses. Different turbulent Prandtl number formulations and their components are discussed in context of strong property variations.     

Difference scheme and numerical simulation based on mixed finite element method for natural convection problem

ntroductionLetΩ R2 beaboundeddomain .Weconsiderthefollowingnon_stationarynaturalconvectionproblem :Problem (Ⅰ )　Findu =(u1,u2 ) ,p ,andTsuchthat,foranyt1>0 ,ut- μΔu +(u· )u + p=λjT　　 ((x ,y ,t) ∈Ω× (0 ,t1) ) ,divu =0　　　　　　　　　 ((x ,y,t) ∈Ω× (0 ,t1) ) ,Tt-ΔT +λu· T =0　　 ((x,y,t) ∈Ω× (0 ,t1) ) ,u =0 ,T =0　　　　　　 ((x,y,t)∈ Ω× (0 ,t1) ) ,u(x ,y ,0 ) =0 ,　T(x,y,0 ) =f(x,y)　　 ((x,y) ∈Ω) ,whereuisthefluidvelocityvectorfield ,pthepressurefield ,Tthet…     

Multiple-Camera Instrumentation of a Single Point Incremental Forming Process Pilot for Shape and 3D Displacement Measurements: Methodology and Results

A multiple-camera system (more than two cameras) has been developed to measure the shape variations and the 3D displacement field of a sheet metal part during a Single Point Incremental Forming (SPIF) operation. The modeling of the multiple-camera system and the calibration procedure to determine its parameters are described. The sequence of images taken during the forming operation is processed using a multiple-view Digital Image Correlation (DIC) method and the 3D reconstruction of the part shape is obtained using a Sparse Bundle Adjustment (SBA) method. Two experiments that demonstrate the potentiality of the method are described.     

Measurement of Sediment Retention in a Sandy Tidal Flat Based on Pressure Drop Model

Sediment can either play an important role in subsurface environments as a food source for bacteria or deteriorate the subsurface environments by its retention. Thus, understanding sediment retention is useful for designing the management of subsurface environments. The pressure drop model derived from the Kozeny–Carman model is experimentally verified by the seepage flow in sand beds. It was found that the water head in the sand bed under steady-state flow and variations of the water head corresponding to changes in the boundary water head could be reproduced by the pressure drop model. As the porosity of the sand bed is taken into account in the pressure drop model, the sediment retention can be predicted from variations of the porosity. Experimental results showed that the water head in the sand bed varied due to sediment retention. This ensured that variances in the porosity of the sand bed could be predicted, leading to the investigation onto sediment retention. A method based on the pressure drop model is proposed to measure temporal variations of the water head in a sandy tidal flat and river water head. From field experiments, the temporal variations of the water head in the tidal flat could be predicted when the porosity of the tidal flat was used. Conversely, it is expected that sediment retention in the tidal flat can be predicted based on the variations of the porosity, if the water head in the tidal flat is observed temporally.     

Sensitivity Resonance and Attractor Morphing Quantified by Sensitivity Vector Fields for Parameter Reconstruction

A novel method of parameter variation reconstruction for systems exhibiting chaotic dynamics is presented. The algorithm reconstructs variations of system parameters without the need for explicit system equations of motion, or knowledge of the nominal parameter values. The concept of a sensitivity vector field (SVF) is developed. This construct captures geometrical deformations to the dynamical attractor of the system in state space. These fields are collected by means of a proposed unique approach referred to as point cloud averaging (PCA). PCA is applied to discrete time series data from the system with nominal parameter values (healthy) and the system with changed parameters. Test variations are reconstructed from an optimal basis of SVF snapshots which is generated by means of proper orthogonal decomposition. The method is applied to two system models, a magneto-elastic oscillator (MEO) and an atomic force microscope (AFM). The method is shown to be highly accurate, and capable of identifying multiple simultaneous variations. The success of the method as applied to an AFM and a MEO indicates a potential for highly accurate readings by exploiting the geometric features of observed chaotic vibrations. An exciting new phenomenon referred to as sensitivity resonance was also observed, and some implications regarding its use in further improving algorithm performance are discussed.An earlier version of this work has been presented at the 20-th Biennial Conference on Mechanical Vibration and Noise, Long Beach, 2005.     

A Note on the Iterative Solution of Stress Problems for Plates and Shallow Shells

A technique based on a refined iterative theory and the numerical method of local variations is developed and used to determine the stress–strain state of transversely isotropic shallow shells and plates. All the components of the stress–strain state and boundary-layer effects are taken into account. The solutions are analyzed for accuracy and convergence.     

Sensitivity of probabilistic fracture mechanics analysis to variations in statistical parameters

Probability of failure (p_f) of a structure is usually calculated for a specified set of statistical parameters (mean, standard deviation, and probability distribution) that characterize random variables. This approach may not be efficient in cases where one would like to know the effect of variations in statistical parameters on the probability of failure. A method based on generating and analyzing randomly selected statistical parameters is proposed. The method consists of generating databases of mean and coefficient of variation (COV = mean/standard deviation) values of relevant fracture mechanics variables through a random process. The method was applied to surface cracks in a flat wide plate loaded under elastic conditions. Probability of failure was calculated for each database record using the first-order reliability method (FORM). Multiple linear regression analyses of the database records were performed with p_f as dependent variable and statistical parameters as independent variables. The predicted p_f values were in very good agreement with the directly calculated p_f values for the specified variations of statistical parameters (±10%, ±15%, and ±20%), except those for fracture toughness and tensile stress, where variations should be limited to ±10% and ±15% ranges.     

Computation of Liquid Hydrazine Depressurisation with a Mach-Uniform Staggered Scheme

A numerical method (pressure-correction method using a staggered grid) is coupled to a thermodynamic model for compressible liquid hydrazine. The method is applied to the venting of liquid hydrazine into space, during which the fluid undergoes a large pressure drop. Below the saturation pressure vaporisation occurs. This takes place near the outlet and induces variations of temperature, which may cause solidification and pipe clogging. In order to assess the risk of phase changes, numerical simulations of the venting line have been performed using a quasi one-dimensional approach. The numerical method can handle compressible flows of fluids with nonconvex equation of state at the low Mach numbers that occur during hydrazine venting. A numerical study of the liquid behaviour during strong depressurisation is performed. The method is validated using experimental data, and allows prediction of pressure evolution and vaporisation location along the pipe. This revised version was published online in July 2006 with corrections to the Cover Date.     

复合材料尖劈和接头端部奇性场的反平面问题研究

提出了一个基于位移的分析尖劈端部奇性位移场和应力场反平面问题的非协调元特征法．该方法与过去原有求解裂纹尖端近似场的有限元特征法有几点不同：(1)导出虚功原理的出发点为二维扇区的散度原理;(2)有限元的单元形式为非协调元;(3)尖劈端部邻域内的位移场假定没有采用奇异变换技术．运用该方法给出了求解正交各向异性复合材料尖劈端部附近奇性应力指数、奇性位移和应力角分布函数的算例．计算结果表明,该方法较原来的有限元特征法所用的单元少而且精度高．     

The suppression of a dynamic instability of an elastic body using feedback control

A theoretical and experimental study is made to determine the feasibility of controlling a thin cantilevered beam subject to a (nonconservative) follower force. A theoretical model is developed using the equations for a thin beam under initial stress and Galerkin's method. An experiment is constructed with the capability of using a variety of feedback loops to control a thin aluminum beam with a tip jet mounted parallel to the chord. A particular control system is chosen for study and an increase of follower force required to destabilize the beam of over 65 per cent is recorded. The theoretical results show good correlation with the experimentally determined stability boundaries and frequency variations with follower force.