Analysis of statically indeterminate non-uniform bar problem in post elastic domain by an iterative variational method

The present study investigates the growth of elastic–plastic front of a statically indeterminate non-uniform bar in post-elastic regime. The solutions of statically indeterminate bar problems are critical in general, because they are not amenable to a ready analytical solution. A clamped axially loaded bar problem becomes indeterminate when the load is concentric, and it result in a singularity point in the domain. In the present bar problem more such singularity points arise when the bar is in post-elastic state, at higher magnitude of concentrated load and the other points come from the yield front location. The computational domain is divided into sub-domains based on the location of singularity points. The formulation is based on von-Mises yield criterion and for linear strain hardening type material behavior. The governing equation is derived through an extension of a variational method in elasto-plastic regime and solution is obtained by using Galerkin's approximation principle. The approximate solution further needs an iterative method to locate the growth in the yield front. The solution algorithm is implemented with the help of MATLAB^® computational simulation software and validation of the formulation is carried out successfully for some reduced problems. The effect of geometry parameters like aspect ratio, slenderness ratio and the type of taperness on the post-elastic performance of the bar is investigated and the relevant results are obtained in dimensionless form. The term bar used in this paper is in generic sense and hence the formulation is applicable for all one dimensional elements, e.g., rods, pipes, truss members, etc.     

Size-dependent forced vibration of non-uniform bi-directional functionally graded beams embedded in variable elastic environment carrying a moving harmonic mass

In this study, general non-uniform material-varying micro-beam models under a moving harmonic load/mass are investigated. Material variation is modeled by combining axial and thickness material grading models using exponential, linear, parabolic and sigmoidal functions. Beam is assumed to be resting on an elastic foundation and in this linear foundation model, foundation modulus is assumed to vary axially with respect to space variable in a non-linear manner ignoring the effect of mass density of foundation on the behavior of micro-beam. Cross-section variation through the length is formulated for both thickness and width variation. Governing equations for such comprehensive beam model is achieved using Hamilton's principle in conjunction with modified couple stress theory to add the scale-effects and solved by discussing explicit and implicit finite element methods with using various-steps and Wilson-theta method. Current methodology is verified using previous studies on simplified problems. A comprehensive parametric study is presented in order to indicate the influence of each design, material and fundamental terms on the forced vibration behavior of such structures under a moving harmonic/constant load/mass. It is shown that by appropriately choosing the material variation in bidirectional functionally graded beams dynamic vibration behavior of such structures could change significantly. Moreover, it is shown that varying cross-section, elastic foundation and type of harmonic moving mass can change the dynamic reaction of the general micro-beam model. From the influence of modified couple stress term on mechanical behavior of such structures it is concluded that this term has crucial effect in varying the dynamic deflections and it is important to acknowledge it in analyzing such structures.     

杆件轴向受迫振动的Galerkin有限元EEP法自适应求解

基于单元能量投影(element energy projection,EEP)法自适应分析在杆件静力问题以及离散系统运动方程组中所取得的成果,以直杆轴向受迫振动为例,研究并建立了一种在时间域和一维空间域同时实现自适应分析的方法.该方法在时间和空间两个维度都采用连续的Galerkin有限元法(finite element method,FEM)进行求解,根据半离散的思想,由空间有限元离散将模型问题的偏微分控制方程转化为离散系统运动方程组,对该方程组进行时域有限元自适应求解;然后再基于空间域超收敛计算的EEP解对空间域进行自适应,直至最终的时空网格下动位移解答的精度逐点均满足给定误差限要求.文中对其基本思想、关键技术和实施策略进行了阐述,并给出了包括地震波输入下的典型算例以展示该法有效可靠.     

Lattice Boltzmann simulation for high-speed compressible viscous flows with a boundary layer

In this study, the lattice Boltzmann method is employed for simulating high-speed compressible viscous flows with a boundary layer. The coupled double-distribution-function lattice Boltzmann method proposed by Li et al. (2007) is employed because of its good numerical stability and non-free-parameter feature. The non-uniform mesh construction near the wall boundary in fine grids is combined with an appropriate wall boundary treatment for the finite difference method in order to obtain accurate spatial resolution in the boundary layer problem. Three typical problems in high-speed viscous flows are solved in the lattice Boltzmann simulation, i.e., the compressible boundary layer problem, shock wave problem, and shock boundary layer interaction problem. In addition, in-depth comparisons are made with the non-oscillatory and non-free-parameter dissipation (NND) scheme and second order upwind scheme in the present lattice Boltzmann model. Our simulation results indicate the great potential of the lattice Boltzmann method for simulating high-speed compressible viscous flows with a boundary layer. Further research is needed (e.g., better numerical models and appropriate finite difference schemes) because the lattice Boltzmann method is still immature for high-speed compressible viscous flow applications.     

Time-domain finite element method to generalized diffusion-elasticity problems with the concentration-dependent elastic constants and the diffusivity

The investigations of mechanical-diffusion coupling are of great importance for the micro-electromechanical devices under non-uniform concentration environment, especially with the development of energy storage technology for a rapid charging system. In recent years there have been many experimental and theoretical studies show that the elastic constants and the diffusivity depend on the concentration of diffusing substances. In view of this, present work aims to study generalized diffusion-elasticity problems considering the concentration-dependent elastic constants and the diffusivity by time-domain finite element method. By using principle of virtual work, the obtained nonlinear finite element equations are solved directly in time domain to minimize precision losses in the application of integrated transformation method, and then the nonlinear solutions can be obtained. As numerical examples, the developed method is used to investigate the transient response of a thick circular plate subjected to the shock loading of the concentration. The results demonstrate that the developed method can faithfully predict the deformation of structure and most importantly the diffusive wave feature in both one-/two-dimensional solids whilst it is commonly difficult to model, especially for two-dimensional case, by using transform method. Parametric studies are performed to evaluate and discuss the effects of concentration-dependent elastic constants and diffusivity on the structural dynamic responses.     

旋转输液管动力稳定性理论分析北大核心CSCD

基于Lagrange原理和假设模态法建立了旋转输液管的动力学模型.通过降阶升维的方法求解系统的特征值问题,并分析了旋转输液管自由振动特性.得到了不同端部集中质量和转速下,系统特征值随流速升高的演变轨迹.揭示了临界流速随系统参数的变化规律.研究发现,内部流体的流动对旋转输液管动力学特性存在显著影响.在某些参数组合下,系统低阶模态能够形成不同形式的内共振关系.预示了旋转输液管模型蕴含丰富的动力学现象.     

Investigation on a rotating FGPM circular disk under a coupled hygrothermal field

In this paper, the distributions of the temperature, moisture, displacement and stress of a functionally graded piezoelectric material (FGPM) circular disk rotating around its axis at a constant angular velocity under a coupled hygrothermal field are presented by a numerical method. The material properties of the FGPM circular disk are assumed to vary along the radial coordinate exponentially. First, the coupled hygrothermal field along the radius of a rotating circular disk is achieved by solving the coupled hygrothermal equations, and then the dynamic equilibrium is solved by utilizing the finite difference method. Finally, numerical results show the effects of functionally graded index, inner radius, angular speed and hygrothermal index on the hygrothermal behaviors of the FGPM circular disk. The results can be useful for the optimal design of rotating FGPM circular disks under a coupled hygrothermal field.     

On the stability of inhomogeneously viscoelastic reinforced bars

There is investigated the stability of inhomogeneously ageing reinforced viscoelastic bars. It is assumed that the strains and stresses in the reinforcement are related by Hooke's law. The properties of the matrix material are described by equations of the theory of viscoelasticity of inhomogeneously ageing solids /1,2/. Under different boundary conditions for the ends of the bar and loading methods an expression is set up for the critical force in stability problems in an infinite time interval. The stability definition taken corresponds to the Liapunov stability definition for the motion of dynamical systems. Estimates of the critical time when the magnitude of the deflection of a viscoelastic bar reaches a given value are obtained for stability problems in a finite time interval. The formulation for the stability problem in a finite time interval starts from the definition of stability of motion of dynamical systems by taking its beginning from the Chetaev work. The dependence of the critical time on the inhomogeneity and the reinforcing parameter is investigated numerically. The stability of viscoelastic unreinforced bars was studied in /3,4/, A survey and bibliography of research associated with the stability problem for viscoelastic bars are available in /5–8/.     

A composite Chebyshev finite difference method for nonlinear optimal control problems

In this paper, a composite Chebyshev finite difference method is introduced and is successfully employed for solving nonlinear optimal control problems. The proposed method is an extension of the Chebyshev finite difference scheme. This method can be regarded as a non-uniform finite difference scheme and is based on a hybrid of block-pulse functions and Chebyshev polynomials using the well-known Chebyshev–Gauss–Lobatto points. The convergence of the method is established. The nice properties of hybrid functions are then used to convert the nonlinear optimal control problem into a nonlinear mathematical programming one that can be solved efficiently by a globally convergent algorithm. The validity and applicability of the proposed method are demonstrated through some numerical examples. The method is simple, easy to implement and yields very accurate results.     

Alternative positional FEM applied to thermomechanical impact of truss structures

This paper presents a formulation to deal with dynamic thermomechanical problems by the finite element method. The proposed methodology is based on the minimum potential energy theorem written regarding nodal positions, not displacements, to solve the mechanical problem. The thermal problem is solved by a regular finite element method. Such formulation has the advantage of being simple and accurate. As a solution strategy, it has been used as a natural split of the thermomechanical problem, usually called isothermal split or isothermal staggered algorithm. Usual internal variables and the additive decomposition of the strain tensor have been adopted to model the plastic behavior. Four examples are presented to show the applicability of the technique. The results are compared with other authors’ numerical solutions and experimental results.     

Dynamic transport scheduling under multiple resource constraints

This paper presents a heuristic for the dynamic vehicle scheduling problem with multiple resource capacity constraints. In the envisaged application, an automated transport system using Automated Guided Vehicles, bottleneck resources are (1) vehicles, (2) docks for loading/unloading, (3) vehicle parking places, and (4) load storage space. This problem is hard, because interrelated activities (loading, transportation, unloading) at several geographical locations have to be scheduled under multiple resource constraints, where the bottleneck resource varies over time. Besides, the method should be suitable for real-time planning. We developed a dedicated serial scheduling method and analyzed its dynamic behavior using discrete event simulation. We found that our method is very well able to find good vehicle schedules satisfying all resource constraints. For comparison, we used a simple approach where we left out the resource constraints and extended the processing times by statistically estimated waiting times to account for finite capacities. We found that our newly designed method finds better schedules in terms of service levels.     

Nonlinear and bilinear models for chemical reactor control

The dynamic behavior of a continuously stirred tank reactor (CSTR) with an exothermic reversible reaction is studied. The balance equations of the reaction lead to a set of highly nonlinear differential equations. For system analysis and control synthesis the dynamic equation are rewritten as state space model. From this nonlinear model a bilinear model is derived. Then, two optimization problems are solved: The time optimal problem for the nonlinear model and the quadratic problem for the bilinear model. In case of the finite time bilinear-quadratic problem a modified Riccati approximation algorithm for a stabilizing feedback controller is presented.     

Semi-exact solution for thermo-mechanical analysis of functionally graded elastic-strain hardening rotating disks

In this paper, distributions of stress and strain components of rotating disks with non-uniform thickness and material properties subjected to thermo-elasto-plastic loading are obtained by semi-exact method of Liao’s homotopy analysis method (HAM) and finite element method (FEM). The materials are assumed to be elastic-linear strain hardening and isotropic. The analysis of rotating disk is based on Von Mises’ yield criterion. A two dimensional plane stress analysis is used. The distribution of temperature is assumed to have power forms with the hotter point located at the outer surface of the disk. A mathematical technique of transformation has been proposed to solve the homotopy equations which are originally hard to be handled. The domain of the solution has been substituted by a new domain through which the unknown variable has been taken out from the argument of the function. This makes the solution much easier. A numerical solution of the governing differential equations is also presented based on the Runge–Kutta’s method. The results of three methods are presented and compared which shows good agreements. This verifies the implementation of the HAM and demonstrates its applicability to provide accurate solution for a very complicated case of strongly high nonlinear differential equations with no exact solution. It is important to notice that compared with other methods, HAM needs significant more computation time and computer hardware requirements which limit its application for those problems that other methods can easily handle them.     

基于非均匀参数化的自由终端时间最优控制问题求解

针对自由终端时间最优控制问题,提出了一种基于非均匀控制向量参数化的数值解法.将控制时域离散化为不同长度的时间段,各时间段长度作为新的控制变量.通过引入标准化的时间变量,原问题转化为均匀参数化的固定终端时间最优控制问题.建立目标和约束函数的Hamilton函数,通过求解伴随方程获得目标和约束函数的梯度,采用序列二次规划(SQP)获得数值解.针对两个经典的化工过程自由终端时间最优控制问题进行仿真研究,验证了所提出算法的可行性和有效性.     

Dynamic analysis of rotating double-tapered cantilever Timoshenko nano-beam using the nonlocal strain gradient theory

Based on the nonlocal strain gradient theory, the coupling nonlinear dynamic equations of a rotating double-tapered cantilever Timoshenko nano-beam are derived using the Hamilton principle. The equation of motion is discretized via the differential quadrature method. The effects of the angular velocity, nonlocal parameter, slenderness ratio, cross-section parameter, and taper ratios are examined and discussed. It is shown that taper ratios and cross-section parameter play a significant role in the vibration response of a rotating cantilever nano-beam. Further as rotational angular velocity increases, the taper ratios and cross-section parameter effect on the frequency response are increased for first modes of vibration.     

Positive solutions and dynamics of a finite difference reaction–diffusion system

We investigate the dynamics and methods of computation for some nonlinear finite difference systems that are the discretized equations of a time-dependent and a steady-state reaction–diffusion problem. The formulation of the discrete equations for the time-dependent problem is based on the implicit method for parabolic equations, and the computational algorithm is based on the method of monotone iterations using upper and lower solutions as the initial iterations. The monotone iterative method yields improved upper and lower bounds of the solution in each iteration, and the sequence of iterations converges monotonically to a solution for both the time-dependent and the steady-state problems. An important consequence of this method is that it leads to a bifurcation point that determines the dynamic behavior of the time-dependent problem in relation to the corresponding steady-state problem. This bifurcation point also determines whether the steady-state problem has one or two non-negative solutions, and is explicitly given in terms of the physical parameters of the system and the type of boundary conditions. Numerical results are presented for both the time-dependent and the steady-state problems under various boundary conditions, including a test problem with known analytical solution. These numerical results exhibit the predicted dynamic behavior of the time-dependent solution given by the theoretical analysis. Also discussed are the numerical stability of the computational algorithm and the convergence of the finite difference solution to the corresponding continuous solution of the reaction–diffusion problem. © 1993 John Wiley & Sons, Inc.     

圆柱壳在径向冲击载荷作用下的弹性脉冲屈曲

当圆柱壳承受径向脉冲载荷时，如果其径厚比大于一特定值，圆柱壳将产生弹性动力屈曲．本文根据有关实验结果，假定变形模态，采用Ｌａｇｒａｎｇｅ方法分析了有限长薄圆柱壳（ａ／ｈ＝４８０）在余弦冲击载荷作用下的弹性脉冲动力屈曲．导出了动力屈曲方程组，借助数值方法求解方程，并与有关计算结果进行了比较．     

Modelling and bifurcation analysis of internal cantilever beam system on a steadily rotating ring

Using Hamilton variation principle, a nonlinear dynamic model of the system with a finite deforming Rayleigh beam clamped radially to the interior of a rotating rigid ring, under the assumption that the constitutive relation of the beam is linearly elastic, is discussed. The bifurcation behavior of the simple system with the Euler-Bernoulli beam is also discussed. It is revealed that these two models have no influence on the critical bifurcation value and buckling solution in the steady state. Then we use the assumption model method to analyse the bifurcation behavior of the steadily rotating Euler-Bernoulli beam and get two different types of bifurcation behavior which physically exist. Finite element method and shooting method are used to verify the analytical results. The numerical results confirm our research conclusion. Project supported by the National Natural Science Foundation of China (Grant No. 19332022) and Space High Technology Foundation of China.     

Eigenvalue analysis of structures with interval parameters using the second-order Taylor series expansion and the DCA for QB

When the uncertainties in interval parameters are fairly large, the current analysis methods, which are usually based on the information of the first-order partial derivatives of eigenvalues, may not work well for the structural eigenvalue problem with interval parameters. To overcome this drawback, in this work, the structural eigenvalue problem with interval parameters is modeled as a series of QB (quadratic programming with box constrains) problems by taking advantage of the information of the second-order partial derivatives of eigenvalues. Then the series of QB problems would be solved by using the DCA (difference of convex functions algorithm) which is turn out be very effective for the QB problem. The specific examples, a concrete frame with sixty bars and a plate discretized with 300 finite elements, are given to show the effectiveness and feasibility of the proposed method compared with other methods.