[0, k i]1 m -Factorizations Orthogonal to a Subgraph

Let G be a graph,k₁,…,k_m be positive integers. If the edges of graph G can be decomposed into some edge disjoint. [0,k₁]-factor. F₁,…,[0,k_m]-factor F_m, then we can say F={F₁,…,F_m, is a [0,k_i]_m¹-factorization of G. If H is a subgraph with m edges in graph G and |E(H)∩E(F_i)|=1 for all 1≤i≤m, then we can call that F is orthogonal to H. It is proved that if G is a..[0,k₁+…+k_m-m+1]-graph, H is a subgraph with m edges in G, then graph G has a. [0,k_i]₁^m-factorization orthogonal to H.     

The strain energy criterion of mixed-mode brittle fracture

In this paper, we proposed a new cr tenon or mixea-moae brittle fracture, i.e., the strain energy criterion, which can be stated as (K_Ⅰ/K_Ⅰc)² +(K_Ⅰ/K_Ⅰc)²+(K_Ⅱ/K_Ⅱc)²=1. This criterion is shown to be in good agreement with known experimental data.In this paper, an experimental criterion:(K_Ⅰ/K_Ⅰc)^m+(K_Ⅱ/K_Ⅱc)ⁿ=1, 1≤_n^m≤2.is also proposed.     

Fixed point theorem of nonexpansive mappings in convex metric spaces

Let X be a convex metric space with the property that every decreasing sequence of nonenply dosed subsets of X with diameters tending to has menemptyintersection. This paper proved that if T is a mapping of a elosed conver nonempty subset K of X into itself satisfying the inequality:d(T_x,T_y)≤ad(x,t)+b{d(x,T_x)+d(y,T_y)}+c{d(x,T_x)+d(y,T_y)}for all x,y in K,where 0≤a<1,b≥0,c≥0,a+c≠0 and a+2b+3c≤1, then T has a unique fixed point in K.     

复杂凸多面体随机紧密堆积组构性能的数值研究:形状参数的影响

颗粒材料的宏观物理力学性能依赖于颗粒堆积体系的细观组构性能,研究颗粒堆积体系的组构性能有重要意义。然而,当前对颗粒堆积体系组构性能的研究集中于球、椭球和正则多面体等规则几何体,还未有对复杂凸多面体颗粒堆积体系组构性能的系统研究。本文基于旋转椭球面黄金螺旋网格构造了一组复杂凸多面体颗粒模型(Poly_κ-ngs),然后基于松弛算法获得了Poly_κ-ngs多面体的随机紧密堆积结构,最后研究了几何形状参数对Poly_κ-ngs多面体随机紧密堆积体系组构性能的影响。结果表明,长径比κ和顶点数量n_gs均对堆积体系的组构性能有影响,κ是主要影响因素。Poly_κ-ngs多面体随机紧密堆积结构中颗粒的位置分布均匀,长径比κ越接近1,顶点数量n_gs越大时,堆积结构表现出更强的位置长程有序性;颗粒方向分布不均匀,长径比κ越远离1,不均匀程度越高;最高堆积分数随长径比κ的增大先增大后减小,在κ=1时达到峰值;配位数分布服从高斯分布,平均配位数随形状参数的变化和堆积分数不同;面-面接触数量随长径比κ的增大先增大后减小,和堆积分数变化规律一致。本研究为复杂凸多面体颗粒的随机紧密堆积提供了数值模拟方案,得出的结论对含有凸多面体颗粒材料的设计和性能优化具有参考意义。     

Nonstationary pandom vibpation analysis of linear elastic structures with finite element method

At present, the finite element method is an efficient method for analyzing structural dynamic problems. When the physical quantities such as displacements and stresses are resolved in the spectra and the dynamic matrices are obtained in spectral resolving form, the relative equations cannot be solved by the vibration mode resolving method as usual. For solving such problems, a general method is put forward in this paper. The excitations considered with respect to nonstationary processes are as follows: P(t)={P_i(t)},P_i(t)=a_i(t)P_i(t), a_i(t) is a time function already known. We make Fourier transformation for the discretized equations obtained by finite element method, and by utilizing the behaviour of orthogonal increment of spectral quantities in random process^[1], some formulas of relations about the spectra of excitation and response are derived. The cross power spectral denisty matrices of responses can be found by these formulas, then the structrual safety analysis can be made. When a_i(t)=l (i= 1,2,…n), the. method stated in this paper will be reduced to that which is used in the special case of stationary process.     

On near-wall hot-wire measurements

A specially constructed hot-wire probe was used to obtain very near-wall velocity measurements in both a fully developed turbulent channel flow and flat plate boundary layer flow. The near-wall hot-wire probe, having been calibrated in a specially constructed laminar flow calibration rig, was used to measure the mean streamwise velocity profile, distributions of streamwise and spanwise intensities of turbulence and turbulence kinetic energy k in the viscous sublayer and beyond; these distributions compare very favorably with available DNS results obtained for channel flow. While low Reynolds number effects were clearly evident for the channel flow, these effects are much less distinct for the boundary layer flow. By assuming the dissipating range of eddy sizes to be statistically isotropic and the validity of Taylor's hypothesis, the dissipation rate ɛ _iso in the very near-wall viscous sublayer region and beyond was determined for both the channel and boundary layer flows. It was found that if the convective velocity U _c in Taylor's hypothesis was assumed to be equal to the mean velocity Uˉ at the point of measurement, the value of (ɛ⁺ _iso)₁ thus obtained agrees well with that of (ɛ ⁺)_DNS for y ⁺ ≥ 80 for channel flow; this suggests the validity of assuming U _c=Uˉ and local isotropy for large values of y ⁺. However, if U _c was assumed to be 10.6u _τ , the value of (ɛ⁺ _iso)₂ thus obtained was found to compare reasonably well with the distribution of (ɛ⁺ _iso)_DNS for y ⁺≤ 15. Received: 31 May 1999/Accepted: 20 December 1999     

海南岛清澜-八门湾潟湖潮汐水动力的模拟研究

基于DELFT3D模型研究了清澜-八门湾潟湖的水位分布和潮汐波内部结构，以及由于人类开垦引起的海岸线变化（以1962年、1985年和2008年为例）对潟湖水动力特性的影响。结果表明，清澜-八门湾潟湖潮汐是由多种分潮耦合而成的复杂驻波，其中K₁，O₁，M₂，S₂和M₄分潮的影响最大。由湾外向湾内传输，由于湾内红树林和浅滩引起的底部损耗增加，M₂，S₂，K₁和O₁分潮幅值减弱；M₄分潮幅值增强，表现出明显的浅水增幅效应；M₂和S₂分潮相位在文教河和文昌江领域表现出明显的干湿效应。不同年代海岸线的研究表明，1985年和2008年间，人类复垦导致潟湖及其潮汐汊道附近的红树林和滩涂区域严重破坏，海岸线缩减，引发了水位降低、纳潮量减少和潟湖潮汐汊道底摩擦弱化，从而削弱了干湿、潮呛和浅水效应。后果提示持续的人类复垦活动将会引发清澜-八门湾潮汐水动力环境的进一步恶化，可能导致未来发生更大的自然灾害。     

The analytical solution for Helmholtz boundary problem in non-horizontally stratified domains

There are N domains D_j(j=0,1,...,N-1) of different physical parameters in the whole space and their interfaces S_j,i+1 are non-horizontally smooth curved surfaces. The following boundary problem is called Hclinholiz boundary problem:V²H⁽¹⁾+K_jH⁽¹⁾=0 (j=0,1,…,N-1)(H⁽⁰⁾-H⁽¹⁾)=δ(S) (δ(S):generalized function)(H¹-Hⁱ⁺¹)=0 (j=0,1,…,N-2)The analytical solution of the above problem is given in this paper.     

核级石墨PFC2D模型的宏细观参数分析

颗粒流软件PFC能够很好地处理非连续介质力学的问题,模拟核级石墨的损伤和断裂过程具有独特的优势。在采用PFC^2D模拟核级石墨时,为了增加模拟的便捷性,实现了PFC^2D细观参数的快速标定。本文首先设计了四因素四水平的正交试验,分析了核级石墨的宏观参数与PFC^2D的细观参数关系;然后基于建立的核级石墨宏细观参数的关系,采用PFC^2D模拟石墨IG11的三点弯曲试验,并从位移场变形、荷载-位移曲线和断裂参数三个方面验证其宏细观参数关系的正确性。结果表明,弹性模量E随平节理模量E_c的增加线性增加;抗拉强度σ_t随平节理抗拉强度σ_b的增加线性增加;泊松比ν随平节理刚度比k_n/k_s的增加先增大后减小;数值模拟的位移场变化和荷载-裂缝口张开位移P-CMOD曲线与试验结果较吻合,此外,断裂参数的最大误差不超过3.57%。     

Experiments on the effect of acceleration on the drag of tapwater bubbles

Three-dimensional (3D) trajectories of spherical air bubbles passing through a converging part of a rectangular channel have been measured. Bubble diameters, d _b, were less than 1 mm and the Reynolds numbers, Re _b, for stagnant tapwater and for mean liquid velocity vˉ _L=0.25 m/s were in about same range. Received: 15 January 2001 / Accepted: 12 June 2001     

近源堆沉积土石体细观结构与抗剪强度指标关联分析

细观结构是认知土石体力学行为本质的关键科学问题。本文通过引入物理学和数学等方法，借助颗粒物质力学理论，从几何排列与接触力的空间分布来定量刻画水平固结与山前坡地堆积两种典型环境下土石体的细观结构特征，并建立其与抗剪强度指标的关联。研究表明，（1）两种环境的土石体在细观结构上存在较大差异。在几何排列上，水平固结环境下的土石体具有长程无序和短程有序的特点，坡地堆沉积环境下的土石体表现出了无序的无定形结构；在接触力与单位接触向量的空间分布上，两者较为相似，绝大多数接触力以小于均值接触力的形式存在，其概率密度曲线P（f）呈幂函数衰减；90%以上接触方位角集中在40°~160°和220°~340°范围内。（2）基于径向分布函数、接触力概率密度和单位接触向量分别定义细观结构的特征量K_a，K_s和K_o，发现三个特征量的增大对内摩擦角呈线性促进作用，对黏聚力呈非线性削弱作用。     

THE BÄCKLUND TRANSFORMATION AND NONLINEAR SUPERPOSITION FORMULA OF SOLUTIONS FOR THE LIOUVILLE’S EQUATION IN HIGHER DIMENSIONS

In this paper, we show that Backlund transformation derived by Leibbrandt et al. for the Liouville’s equation in three spatial dimensions, ▽²a=expa, ▽²=-∂_x²+∂_y²+∂_z² can be decomposed into several Backlund transformations for the same equation in two spatial dimensions, moreover, the superposition formula which is derived from this transformation is actually invalid, thus the discussions based on that formula is incorrect as well. We also considered some results about the Liouville’s equation in N spatial dimensions.     

An analysis of the effects of the prostatic hypertropy on he urinary flow in lower-urinary-tract

In this paper an analytic model corresponding to the collapsible tube for analysing the urinary flow in lower-urinary-tract is set up from physiologic background.By analysing the model it is found that the self-excited oscillations can both occur in the region of negative and positive slope of P_n-Q_n characteristic. So this paper extends the results of Conrad1’1, Griffiths121, Conrad, Cohen and McQueenI3! and others that the self-excited oscillations can only occur in the region of negative slope of P_n-Q_n characteristic. The effects of prostatic hypertropy on the flow parameters in lower-urinary-tract is discussed in detail by numerical calculations. The results show that it is possible to know the conditions of prostatic hypertropy according to the changes of bladder pressure, outlet urinary velocity and other parameters. From these results a theoretical method to detect and diagnose prostatic hypertropy is provided.     

Rheological aspects of dense lignite-water suspensions; time dependence, preshear and solids loading effects

An investigation of the rheological properties of dense lignite-water suspensions is reported here in order to evaluate the parameters that affect the rheology of these suspensions. Different types of particle size distributions were achieved via wet grinding with particle sizes down to 1 μm. Polyelectrolytes and surfactants were used to improve the stability and rheology of suspensions. The achieved solid volume fraction, φ, with acceptable rheological behavior was 0.45, while the ratio of φ/φ _m was close to 0.85; φ _m is the maximum solid volume fraction. Various types of flow tests were performed on lignite-water suspensions with various preshear times and stress levels. Typical shear rates varied in the range between 10^–2 to 10² s^–1. The low shear behavior was found to be quite different from the high shear behavior for the most concentrated suspensions, with a plateau value at low shear rates indicating the appearance of a yield stress. These plateau values, however, depend on shear history, which is responsible for the development of different structures in the sample. The viscosity curves corresponding to the ascending and descending parts of the flow curve were found to be different; these flow curve parts can be described either by well established models or by modified ones. The non-Newtonian time dependent behavior of the lignite-water suspensions is attributed to the high value of the ratio φ/φ _m , the polydispersity of the particle size distribution, and the non-spherical shape of the lignite particles.     

基于转捩/尺度适应模型与FW-H声学方程的气动噪声数值模拟研究

采用计算流体力学(CFD)同"声类比"相结合的方法进行噪声模拟,利用CFD数值模拟 MD30P30N多段翼型失速攻角附近流场以及其气动特性来校核近场精度,进一步通过结合可穿透数据面的FW-H声学方程进行气动噪声分析。为了准确捕捉近场流场信息,为噪声预测提供可靠的声源精度,本文基于k-ω剪应力输运(SST)湍流模型,建立了尺度适应分离流(SAS)模型,并采用γ-Re_θt转捩模型耦合k-ω SST湍流模型建立了边界层转捩数值模拟技术;充分利用尺度适应模型在边界层表现为雷诺平均(RANS)方法这一特点,将γ-Re_θt转捩模型与SAS模型结合,建立针对包含转捩、分离现象复杂流场的数值模拟技术。文中以RANS为控制方程,分别采用全湍流k-ω SST模型、γ-Re_θtSAS转捩/分离流模型对多段翼进行近场数值模拟计算,结合可穿透数据面的FW-H声学方程进行气动噪声预测,在分析其对流场及气动噪声影响的基础上,得出了几点结论。     

Havelock型格林函数振荡项数值积分的稳定性研究

Havelock型格林函数的传播项被积函数是高频振荡且奇异的复变函数,文献[4]引入变量代换获得了一种兼具积分效率和精度的积分方法, 本文研究了该方法的积分稳定性,发现该方法仍存在如下的积分困难: (1) θ=γ时复函数中分母为零引起的计算溢出;(2) θ=π/2是复函数在y及z方向偏导数的无穷间断点;(3) 场点与源点横坐标相同时伪奇异点变为真奇点。针对这些积分困难,采用极限公式计算θ=γ处复函数的值避免计算溢出;在保证积分精度的前提下采用截断方法略去θ=π/2邻近区域的积分消除无穷间断处的奇异;针对(3)采用分区法处理以避开原被积函数的高频振荡,并消除奇异性。伪奇异性存在的条件是场点必须在点源传播波的传播范围内,伪奇异点最多为2个。     

Analysis of surface crack on forward extruding bar during axisymmetric cup-bar combined extrusion process

In this paper, the kinematically admissible velocity field with surface crack on forward extruding bar is put forward during the axisymmetric cup-bar combined extrusion process, in accordance with the results of model experiments.On the basis of velocity field, the necessary condition for surface crack formation on the forward extruding bar is derived, with the help of upper bound theorem and the minimum energy principle. Meanwhile, the relationships between surface crack formation and combination of reduction in area for the part of forward and backward extursions (ε_b,ε_f) relative residual thickness of billet (T/R₀),frictional factor (m) or relative land length of ram and chamber (l_b|R₀,l_f|R₀) are calculated during the extrusion process. Therefore, whether the surface crack on forward exturding bar occurs can be predicted before extruding the lower-plasticity metals for axisymmetric cup-bar combined extrusion process.The analytical results agree very well with experimental results of aluminium alloy LY12 (ASTM 2024) and LC4 (ASTM 7075).     

BASIC EQUATIONS, THEORY AND PRINCIPLE OF COMPUTATIONAL STOCK MARKET (Ⅲ)-BASIC THEORIES

By basic equations, two basic theories are presented: 1. Theory of stock’s value v^*(t)=v^*(0) exp (ar₂^*t);2. Theory of conservation of stock’s energy. Let stock’s energy Φbe defined as a quadratic function of stock’s price v and its derivative v, Φ=Av²+Bv+Cv²+Dv, under the constraint of basic equation, the problem was reduced to a problem of constrained optimization along optimal path. Using Lagrange multiplier and Euler equation of variation method, it can be proved that Φkeeps conservation for any v,v. The application of these equations and theories on judgement and analysis of tendency of stock market are given, and the judgement is checked to be correct by the recorded tendency of Shenzhen and Shanghai stock markets.     

基于COMSOL的脆性材料相场断裂模型

相场法通过一系列微分方程描述材料断裂过程,避免了繁琐的裂纹面追踪,在模拟裂纹的萌生、扩展和分叉等方面具有优势。介绍了基于相场法的脆性材料断裂模型,给出了脆性材料断裂问题相场法控制方程的推导过程,提出了基于分步迭代法在COMSOL中实现脆性材料相场断裂模型的方法。再现了脆性材料单元模型和单边缺口平板受拉及受剪作用下的开裂过程,模拟的裂纹扩展路径与已有文献的结果相近,验证了程序的合理性。针对脆性材料相场断裂模型包含的诸多参数,采用Morris法对影响荷载-位移关系的脆性材料断裂模型参数进行了全局敏感性分析,结果表明,杨氏模量(E)、临界能量释放率(G_c)和位移增量(Δu_x)是影响模型荷载-位移关系输出结果的主要参数。基于COMSOL实现的相场断裂模型能够有效模拟脆性材料的裂纹萌生和扩展断裂过程,模型参数E,G_c和Δu_x对材料断裂性能的提升或模型参数反演效率的提高具有重要影响。     

柱塞泵多目标优化设计及CFD仿真分析

液压泵噪声是液压系统的主要噪声源,针对轴向柱塞泵的流致振动噪声,提出一种改善泵配流特性的设计方案。首先,根据柱塞泵的工作原理对柱塞腔压力特性和泵出口流量特性准确建模并求解。通过分析压力冲击和流量脉动对错配角（φ₀）的响应,得φ₀=4°为佳。利用一种多目标遗传算法（NSGA-Ⅱ）,以减小压力超调量和流量脉动率为目标,对三角槽结构进行了优化;并获得该多目标优化问题的Pareto最优解集,通过对最优解集的分析知,深度角θ₁=16°且宽度角θ₂=85°时较为理想。最后,为了验证模型的正确性,建立流体域计算流体动力学（CFD）模型,对比两种模型计算结果发现吻合较好,能够相互验证。利用CFD分析结果可视化的特点,从柱塞泵流场的角度,进一步分析了泵压力冲击以及流量脉动产生的原因。