Application of the modified iteration method to nonlinear vibration of corrugated circular plates

In this paper, the modified iteration method is successfully extended to investigate the nonlinear free vibration of corrugated circular plates with full corrugations. The analytical relation for the amplitude-frequency response of corrugated circular plates is obtained and discussions on the influences of geometrical parameters on vibration behaviours of corrugated circular plates are made. The present results are practically important in the design of elastic elements in precision instruments.     

Large deflection problem of thin orthotropic circular plate with variable thickness under uniform pressure

Basic equations for large deflection theory of thin orthotropic circular plate with variable thickness are derived in this paper. The modified iteration method is adopted to solve the large deflection problem of thin orthotropic circular plate with variable thickness under uniform pressure. If ε=0, then the solution derived from the result in this paper coincides completely with the result given by J. Nowinski (using perturbation method) for solving large deflection problem of thin orthotropic circular plate with constant thickness under uniform pressure.     

Large deplection of corrugated circular plate with a plane central region under the action of concentrated loads at the center

This paper solves the problem of large deflections for corrugated circular plates with a plane central region under the action of concentrated loads at the center by means of the non-linear bending theory for anisotropic circular plates. Using the modified iteration method, the characteristic relation of the plate is obtained. This formula may be applied directly to design of elastic elements of measuring instruments.     

LARGE DEFLECTION PROBLEM OF THIN ORTHOTROPIC CIRCULAR PLATE ON ELASTIC FOUNDATION WITH VARIABLE THICKNESS UNDER UNIFORM PRESSURE

ＬＡＲＧＥＤＥＦＬＥＣＴＩＯＮＰＲＯＢＬＥＭＯＦＴＨＩＮＯＲＴＨＯＴＲＯＰＩＣＣＩＲＣＵＬＡＲＰＬＡＴＥＯＮＥＬＡＳＴＩＣＦＯＵＮＤＡＴＩＯＮＷＩＴＨＶＡＲＩＡＢＬＥＴＨＩＣＫＮＥＳＳＵＮＤＥＲＵＮＩＦＯＲＭＰＲＥＳＳＵＲＥ（王嘉新）（刘杰）ＬＡＲＧ...     

Nonlinear stability of thin elastic circular shallow spherical shell under the action of uniform edge moment

In this paper, nonlinear stability of thin elastic circular shallow spherical shell under the.action of uniform edge moment is considered by the modified iteration method to obtain second and third approximations to decide the upper and lower critical loads. Results are plotted in curves for the engineering use and are compared with results of Hu Hai-chang’s. We also investigate the neighbour situation of the critical point, i.e. the double points of the upper and lower critical loads and denote the range of validity of the second approximation. In the end, we obtain the special case, the design formulas of rigidity and stress as well as the corresponding curves as v=0.3 of large deflection of circular plate under the same load. These results are also compared with Huang Tse-yen’s.     

Nonlinear analysis of a corrugated circular plate under combined lateral loading

In this paper, nonlinear bending of a corrugated circular plate with a plane central region under the combined action of uniformly distributed load and a concentrated load at the center has been investigated by using large deflection theories of isotropic and anisotropic circular plates. The quite accurate analytical solutions for rigidly as well as loosely clamped edge conditions have been obtained following the modified iteration method. Supported by the Science Fund of the Chinese Academy of Sciences.     

An iteration method for integral equations arising from axisymmetric loading problems

Let the concentrated forces and the centers of pressure with unknown density functions x(ξ) and y(ξ) respectively be distributed along the axis z outside the solid, then one can reduce an axismmetric loading problem of solids of revolution to two simultaneous Fredholm integral equations. An iteration method for solving such equations is duscussed. A lemma equivalent to E. Rakotch’s contractive mapping theorem and a theorem concerning the convergent proof of the iteration method are presented.     

Splitting method for two-dimensional phreatic flow equation

ＳＰＬＩＴＴＩＮＧＭＥＴＨＯＤＦＯＲＴＷＯ－ＤＩＭＥＮＳＩＯＮＡＬＰＨＲＥＡＴＩＣＦＬＯＷＥＱＵＡＴＩＯＮＸｕＳｈａｏｈｕｉ（徐绍辉）ＺｈｕＸｕｅｙｕ（朱学愚）ＺｈｕＧｕｏｒｏｎｇ（朱国荣）（ＲｅｃｅｉｖｅｄＭａｙ１１,１９９５;Ｃｏｍｍｕｎｉｃａｔ...     

‘Velocity’ Finite Element Method for Dynamic Problem

After analysing the essential features of successive integration method taking displacement as variable by N. M. Newmark and E. L. Wilson et al, this paper presents a Velocity Element Method, taking velocity as variable for the solution of the initial value problem.A simplified scheme is offered for the non-damping system, and the stability is also discussed. Owing to the fact that this simplified scheme for non-damping and apparent static damping is explicit in form, it is unnecessary to solve the algebraic system of equations at every time interval, consequently the amount of computation is greatly reduced. For non-linear dynamic problems, this scheme may be used to obtain fairly good initial values for iteration.An extended form of elocity Element is presented for the arbitrary damping system. For the non-linear cases, the incremental Velocity iteration scheme is adopted and its convergence proved. Some discussions have been given on artificial damping and the effect of the parameter.Finally, the results of numerical calculatio of some typical problem are given in the appendix.     

REPRESENTATION THEOREM AND ONE-ITERATION THEOREM FOR FREDHOLM INTEGRAL EQUATION OF THE FIRST KIND Ax=y

In this paper,two theorems are presented.The representation theorem states:if the Fredholm integral equation of the first kind Ax=y,with bounded L_2 kernel,has a uniquesolution(?),then(?),whereThe one-iteration theorem states:(?)can be achieved in one iteration by(?)=x_0 g_0 A*(y-Ax_0)iff one of the following conditions is satisfied.     

Modified iterative method for augmented system

The successive overrelaxation-like(SOR-like) method with the real parameters ω is considered for solving the augmented system. The new method is called the modified SOR-like(MSOR-like) method. The functional equation between the parameters and the eigenvalues of the iteration matrix of the MSOR-like method is given. Therefore, the necessary and sufficient condition for the convergence of the MSOR-like method is derived. The optimal iteration parameter ω of the MSOR-like method is derived. Finally, the proof of theorem and numerical computation based on a particular linear system are given, which clearly show that the MSOR-like method outperforms the SOR-like(Li, C. J., Li, B. J., and Evans, D. J. Optimum accelerated parameter for the GSOR method. Neural, Parallel Scientific Computations, 7(4), 453–462(1999)) and the modified symmetric SOR-like(MSSOR-like) methods(Wu, S. L., Huang, T. Z., and Zhao, X. L. A modified SSOR iterative method for augmented systems.     

Preferential Flow-Paths Detection for Heterogeneous Reservoirs Using a New Renormalization Technique

We have devised a renormalization scheme which allows very fast determination of preferential flow-paths and of up-scaled permeabilities of 2D heterogeneous porous media. In the case of 2D log-normal and isotropically distributed permeability-fields, the resulting equivalent permeabilities are very close to the geometric mean, which is in good agreement with a rigorous result of Matheron. It is also found to work well for geostatistically anisotropic media when comparing the resulting equivalent permeabilities with a direct solution of the finite-difference equations. The method works exactly as King's does, although the renormalization scheme was modified to obtain tensorial equivalent permeabilities using periodic boundary conditions for the pressure gradient. To obtain an estimation of the local fluxes, the basic idea is that if at each renormalization iteration all the intermediate renormalized permeabilities are stored in memory, we are able to compute -- ad reversum -- an approximation of the small-scale flux map under a given macroscopic pressure gradient. The method is very rapid as it involves a number of calculations that vary linearly with the number of elementary grid blocks. In this sense, the renormalization algorithm can be viewed as a rapid approximate pressure solver. The exact reference flow-rate map (for the finite-difference algorithm) was computed using a classical linear system inversion. It can be shown that the preferential flow paths are well detected by the approximate method, although errors may occur in the local flow direction.     

Exact and approximate solution of the regenerator equation for the case of high heat exchange and moderate heat capacity

If the inflow temperatures are constant the regenerator equation can be transformed into an integral equation. The latter is solved numerically by an iteration procedure. From the results, a simple approximate solution is found for the case of high heat exchange (30<<300) and moderate heat capacity (1.5<<3).     

非牛顿流体经波纹管流动的阻力特性

本文对Carreau流体经波纹管的蠕动流的阻力特性用差分方法作了数值求解。先施行坐标变换,提出一种保证周期性条件的混合迭代法,能较快得到问题的收敛解。流体的物质参数以及流动区域的几何参数对流动阻力特性的影响作了讨论,同时指出了直径均匀的毛管模型作为渗流模型的不足之处。     

均布载荷作用下具有光滑中心波纹膜片的非线性分析

采用轴对称旋转壳体的简化Reissner方程，研究了在均布载荷作用下具有光滑中心波纹膜片的非线性弯曲问题。应用格林函数方法，波纹膜片的非线性边值问题化为了非线性积分方程的求解。为了求解积分方程并防止发散，一个插值参数被引入到迭代格式中。计算表明，当载荷很小时，任何插值参数值均能保证迭代的收敛性，取插值参数值接近或等于1获得较快的收敛速度，而当载荷较大时，插值参数值不能取得过大。绘出了波纹膜片的特征曲线，得到的特征曲线可供设计参考。可以断言，当载荷不大时，特征曲线是近似线性的，随着载荷的增大，特征曲线开始向上弯曲，明显偏离线性。本文中提出的解决方法适应于任意轴向截面的波纹壳体。     

Viscous Film Flow down Corrugated Surfaces

A theoretical analysis of a downward viscous film flow on corrugated surfaces is reported. The study is based on Navier–Stokes equations (for one and twodimensional surfaces) and on an integral model (for a threedimensional surface with double corrugation). The calculations were carried out in a wide range of Reynolds numbers and geometric characteristics of the surface with due allowance for the surfacetension force. The shape of the free surface of the liquid film and other characteristics of the flow are calculated. It is shown that, in the case of a onedimensional surface, there exists a range of parameters where the flow is predominantly governed by surfacetension forces; this flow can be adequately treated with the integral approach. In this range of parameters, on the surface with double corrugation, the average quantities of the downward flow in wide corrugation valleys are determined by the finetexture geometry.     

Bending of corner-supported rectangular plate under concentrated load

In this paper the solution for the bending of corner-supported rectangular plate under concentrated load at any point (α/2, η) of the middle line of the plate is given by means of a conception called modified simply supported edges and the method of superposition. Some numerical example is presented. The solution obtained by this method checks very nicely with what was obtained by G.T. Shih^[3] by means of spline finite element method when η=d/2. This shows that this method of solution is satisfactory.     

基于能量变分法的钢底板波形钢腹板组合箱梁畸变效应分析

为合理分析钢底板波形钢腹板梯形箱梁的畸变效应,按各板件面内外抗弯刚度不变的原则将全截面等效为钢材,利用圣维南原理考虑顶底板对波形钢腹板的约束作用,修正畸变扇性坐标分布模式,基于能量变分法建立畸变控制微分方程。与已有文献及有限元进行对比分析,并研究腹板俯角和波形钢腹板厚度变化对畸变翘曲正应力的影响。结果表明,本文解析解与文献解及ANSYS解均吻合较好;基于圣维南原理修正后的扇形坐标分布模式更合理;利用本文等效方法亦可分析传统波形钢腹板组合箱梁的畸变效应;腹板俯角的设置有利于减小畸变翘曲正应力;波形钢腹板厚度变化对腹板与底板交接处的畸变翘曲正应力影响显著。     

Electrode layers in flows of an inviscid plasma with good conductivity

The structure of the electromagnetic electrode layers that are produced in flows across a magnetic field by a completely ionized and inviscid plasma with good conductivity and a high magnetic Reynolds number is examined in a linear approximation. Flow past a corrugated wall and flow in a plane channel of slowly varying cross section with segmented electrodes are taken as specific examples. The possibility is demonstrated of the formation of nondissipative electrode layers with thicknesses on the order of the Debye distance or electron Larmor radius and of dissipative layers with thicknesses on the order of the skin thickness, as calculated from the diffusion rate in a magnetic field [2].In plasma flow in a transverse magnetic field, near the walls, along with the gasdynamie boundary layers, which owe their formation to viscosity, thermal conductivity, etc. (because of the presence of electromagnetic fields, their structures may vary considerably from that of ordinary gasdynamic layers), proper electromagnetic boundary layers may also be produced. An example of such layers is the Debye layer in which the quasi-neutrality of the plasma is upset. No less important, in a number of cases, is the quasi-neutral electromagnetic boundary layer, in which there is an abrupt change in the frozen-in parameter k=B/p (B is the magnetic field and p is the density of the medium). This layer plays a special role when we must explicitly allow for the Hall effect and the related formation of a longitudinal electric field (in the direction of the veloeiryv of the medium). We will call this the magnetic layer. The magnetic boundary layer can be dissipative as well as noudissipative (see below). The dissipative magnetic layer has been examined in a number of papers: for an incompressible medium with a given motion law in [1], for a compressible medium with good conductivity in [2], and with poor conductivity in [3]. In the present paper, particular attention will be devoted to nondissipative magnetic boundary layers.