解任意形状扁轴对称体Stokes流动的连续奇点分布法

本文以Sampson球形无穷级数作为基本奇点,采用分段等强度和分段二次抛物分布两种体内连续分布法解任意形状扁轴对称体的Stokes流动.通过扁球的无界绕流问题,对这两种方法的收敛性,精度和适用范围做了检验和比较.结果表明,在一定的范围内,无论是阻力系数或压力分布,它们的计算结果都和精确解符合得很好,而且,随着分布函数逼近程度的提高,其收敛性得到改善,适用范围也随之扩大.作为一般算例,分别用这两种方法解决了卡西尼扁卵形体的绕流问题,得到了一致的结果.最后,用分段二次连续分布法计算了具有一定生理意义的红细胞体的Stokes流动,首次得到了它的阻力系数和表面压力分布.     

厚球壳与实心球轴对称问题的一般解

本文试图从更一般的三维问题基本方程出发研究任意厚球壳与实心球的轴对称问题.对于受任意轴对称载荷的厚球壳和实心球体,文中运用加权残值法给出了以Legendre级数表示的一般解.     

在任意不定常温度场和任意法向动载荷联合作用下中心开孔圆底扁球壳的动力问题

本文得出了在任意不定常温度场和任意法向动载荷联合作用下中心开孔圆底扁球壳的动力问题的解析解.我们假设温度沿壳体厚度直线分布.在第一部分.我们研究了常用边界条件下的中心开孔圆底扁球壳的自由振动.作为例子,我们计算了一边缘夹紧的扁球壳的自然基频(m=0),所得结果与E.Reissner[1]的结果作了比较.频率方程的解法是钱伟长[2]提出来的.这将在附录3中介绍.在第二部分,我们研究了在任意谐温度场和任意谐法向动载荷联合作用下的中心开孔圆底扁球壳的强迫振动.在第三部分,我们研究了在任意不定常温度场和任意法向动载荷联合作用下的具有初始条件的上述壳体的强迫振动.在附录1和2中,我们讨论了如何用应力函数来表示位移边界条件和m=1情形的边界条件.     

简支梯形底扁球壳自由振动问题的准Green函数方法

以简支梯形底扁球壳的自由振动问题为例,详细阐明了准Green函数方法的思想.即利用问题的基本解和边界方程构造一个准Green函数,此函数满足了问题的齐次边界条件,采用Green公式,将简支梯形底扁球壳自由振动问题的振形控制微分方程化为两个耦合的第二类Fredholm积分方程.边界方程有多种选择,在选定一种边界方程的基础上,可以通过建立一个新的边界方程来表示问题的边界,以克服积分核的奇异性.最后由积分方程的离散化方程组有非平凡解的条件,求得固有频率.数值结果表明,该方法具有较高的精度.     

解任意形状非细长长轴对称体Stokes流动的一种新方法

本文提出离散奇点线分布法和连续奇点线分布法解决任意形状非细长长轴对称体的Stokes流动。取Sampson球形无穷级数为基本奇点。通过长球无界绕流问题检验了上述两种方法的收敛性及精度。与精确解比较表明,在一定的细长比下无论在阻力系数或压力分布上,这两种方法都具有良好的收敛性及高精度的计算结果,而且连续奇点线分布法和离散奇点线分布法相比较具有更好的收敛性能。最后,作为计算任意形状非细长长轴对称体的一个例子,利用这两种方法计算了卡西尼卵形体的阻力系数及压力分布值,得到了收敛的一致结果。     

球形扁壳超临界变形的步进求和计算

本文用步进求和法计算了球形扁壳第二类失稳问题,在球扁壳超临界变形计算上给出了优于一级近似结果,解决了该问题无法求二级近似解的困难.算例表明步进求和法收敛于二级近似解.     

中心受集中载荷的固定夹支边圆板和圆底扁球壳的卡门方程的精确解

由于卡门方程的非线性性和耦合性,使得寻求精确解的困难很大。迄今为止,除了少数未从数学上严格证明其收敛性的精确解外,大多数均采用近似方法求解。本文将卡门方程化为非线性奇异耦合的积分方程组,运用迭代法求得了连续函数序列。通过证明其一致收敛性,得到了中心受集中载荷作用的固定夹紧边界的圆板和圆底扁球壳的卡门方程的精确解的解析式及其收敛性证明。     

半离散Jacobi-球面调和谱格式用于流体低马赫数流动

本文提出了一种用于解决流体在球内以低马赫数流动问题的半离散Jacobi-球面调和谱格式,并证明了它的稳定性和收敛性.所用技巧也可应用于球形区域上其它非线性问题.     

波纹壳的格林函数方法

应用轴对称旋转扁壳的基本方程,研究了在任意载荷作用下具有型面锥度的浅波纹壳的非线性弯曲问题·　采用格林函数方法,将扁壳的非线性微分方程组化为非线性积分方程组·　再使用展开法求出格林函数,即将格林函数展成特征函数的级数形式,积分方程就成为具有退化核的形式,从而容易得到非线性代数方程组·　应用牛顿法求解非线性代数方程组时,为了保证迭代的收敛性,选取位移作为控制参数,逐步增加位移,求得相应的载荷·　在算例中,研究了具有球面度的浅波纹壳的弹性特征·　结果表明,由于型面锥度的引入,特征曲线发生显著变化,随着荷载的增加,将出现类似扁球壳的总体失稳现象·　本文的解答符合实验结果·     

考虑横向剪切的对称层合圆柱正交异性扁球壳的非线性稳定问题

本文建立了计及横向剪切的复合材料对称层合圆柱正交异性扁球壳的大挠度理论。应用修正迭代法,研究了均布压力作用下对称层合圆柱正交异性固定边扁球壳的非线性稳定问题,得到了临界荷载的解析解。此解可直接应用到工程设计中。     

三维对流扩散方程的三种高精度分裂格式

在算子分裂法思想的基础上,将两种高精度的离散格式推广应用于三维对流扩散方程,同时对经典ADI格式的对流项做了改进,改进后的格式的对流项对空间具有4阶精度,而经典ADI格式对空间只有2阶精度,由此可见,提高了该格式的实用性．最后对两种典型的浓度场进行了数值模拟,将3种格式的计算结果与解析解以及其它传统差分格式的计算结果进行了对比,得出当Peclet数不大于5时,3种格式均获得了令人满意的数值结果,说明推广的这三种方法具有很高的准确性和可靠性．     

On the theory of dynamic thermoelasticity : PMM, vol. 42, no. 6, 1978, pp. 1093–1098

Use of the causality principle as radiation condition in dynamical problems of thermoelasticity is proposed. It follows from an analysis of the fundamental mathematical models describing the thermoelastic behavior of a continuous medium and used in the solution of specific problems, that some will yield physically unrealizable solutions. To eliminate the ambiguity in the solution which occurs, an approach is possible which has an explicit physical meaning and is based on the causality principle [1, 2]; it is required that the time source not yield a response earlier than the time of starting up of the source. Different kinds of radiation conditions of the Sommerfeld type are known in thermoelasticity problems [3 – 6].

To extract the unique solution in dynamical thermoelasticity problems, it is proposed in this paper to use the causality principle, which is equivalent to the requirement of analyticity of the solution in the upper half of the complex frequency plane; there are studied the analytic properties of the solutions of the fundamental boundary value problems for the models used most often for thermoelastic media, and there are made deductions about their physical realizability.     

Hamiltonian system-based new analytic free vibration solutions of cylindrical shell panels

This paper deals with the classical challenging free vibration problems of non-Lévy-type cylindrical shell panels, i.e., those without two opposite edges simply supported, by a Hamiltonian system-based symplectic superposition method. The governing equations of a vibrating cylindrical panel are formulated within the Hamiltonian system framework such that the symplectic eigen problems are constructed, which yield analytic solutions of two types of fundamental problems. By the equivalence between the superposition of the fundamental problems and the original problem, new analytic frequency and mode shape solutions of the panels with four different combinations of boundary conditions are derived. Comprehensive benchmark results are tabulated and plotted, which are useful for validation of other numerical/approximate methods. The primary advantage of the developed approach that no pre-determination of solution forms is needed enables one to pursue more analytic solutions of intractable shell problems.     

The method of fundamental solutions for linear diffusion-reaction equations

The paper presents an extension of the solution procedure based on the method of fundamental solutions proposed earlier in the literature for solving linear diffusion reaction equations in nonregular geometries in two and three dimensions. The solution procedure utilizes the fundamental solution to the problem along with boundary collocation to result in a grid-free numerical scheme. A new heuristic for source location is utilized along with orthogonal collocation in nonsmooth domains to improve the accuracy of the solution. The efficacy of the solution procedure is demonstrated for a variety of problems in nonregular simply and multiply connected geometries with nonuniform boundary conditions.     

Determination of Stackelberg–Nash equilibria using a sensitivity based approach

A sensitivity based approach is presented to determine Nash solution(s) in multiobjective problems modeled as a non-cooperative game. The proposed approach provides an approximation to the rational reaction set (RRS) for each player. An intersection of these sets yields the Nash solution for the game. An alternate approach for generating the RRS based on design of experiments (DOE) combined with response surface methodology (RSM) is also explored. The two approaches for generating the RRS are compared on three example problems to find Nash and Stackelberg solutions. For the examples presented, it is seen that the proposed sensitivity based approach (i) requires less computational effort than a RSM-DOE approach, (ii) is less prone to numerical errors than the RSM-DOE approach, (iii) has the ability to find multiple Nash solutions when the Nash solution is not a singleton, (iv) is able to approximate nonlinear RRS, and (v) on one example problem, found a Nash solution better than the one reported in the literature.     

Two-dimensional steady-state general solution for isotropic thermoelastic materials with applications. I: General solutions and fundamental solutions

Four steady-state general solutions are derived in this paper for the two-dimensional equation of isotropic thermoelastic materials. Using the differential operator theory, three general solutions can be derived and expressed in terms of one function, which satisfies a six-order partial differential equation. By virtue of the Almansi’s theorem, the three general solutions can be transferred to three general solutions which are expressed in terms of two harmonic functions, respectively. At last, a integrate general solution expressed in three harmonic functions is obtained by superposing the obtained two general solutions. The proposed general solution is very simple in form and can be used easily in certain boundary problems. As two examples, the fundamental solutions for both a line heat source in the interior of infinite plane and a line heat source on the surface of semi-infinite plane are presented by virtue of the obtained general solutions.     

A branch-and-price approach for harvest scheduling subject to maximum area restrictions

Recently, research on exact methods has been undertaken to solve forest management problems subject to constraints on the maximum clearcut area by using the area restriction model approach. Three main basic integer programming models for these problems have been discussed in the literature: the so-called cluster, path and bucket formulations. Solving these models via branch-and-bound, where all variables and constraints are used a priori, is adequately suited for real problems of a small to medium size, but is not appropriate for larger problems. In this paper, we describe a branch-and-price approach for the cluster model, and we show that this formulation dominates the bucket model, by completing the results of the dominance relationships between the bounds of the three models. Branch-and-price was tested on real and hypothetical forests ranging from 45 to 2945 stands and temporal horizons ranging from three to twelve periods were employed. Results show that the solutions obtained by the proposed approach stood within 1% of the optimal solution and were achieved in a short computation time. It was found that branch-and-bound was unable to produce solutions for most forests from 850 stands with either eleven or an average number of stands per clearcut greater or equal than eight.     

The Eshelby tensor

The problem of the existence of a tensor that is inverse to the well-known Eshelby tensor, which connects the free homogeneous and hindered strains of an ellipsoidal elastic inclusion undergoing transformation, is investigated. It is shown that this tensor exists for inclusions in the form of oblate and prolate spheroids in isotropic elastic space. Certain applications are considered, in particular problems of determining the stresses in ellipsoidal rigid and rigid plastic inclusions.     

A data-driven neural network approach to optimal asset allocation for target based defined contribution pension plans

A data-driven Neural Network (NN) optimization framework is proposed to determine optimal asset allocation during the accumulation phase of a defined contribution pension scheme. In contrast to parametric model based solutions computed by a partial differential equation approach, the proposed computational framework can scale to high dimensional multi-asset problems. More importantly, the proposed approach can determine the optimal NN control directly from market returns, without assuming a particular parametric model for the return process. We validate the proposed NN learning solution by comparing the NN control to the optimal control determined by solution of the Hamilton–Jacobi–Bellman (HJB) equation. The HJB equation solution is based on a double exponential jump model calibrated to the historical market data. The NN control achieves nearly optimal performance. An alternative data-driven approach (without the need of a parametric model) is based on using the historic bootstrap resampling data sets. Robustness is checked by training with a blocksize different from the test data. In both two and three asset cases, we compare performance of the NN controls directly learned from the market return sample paths and demonstrate that they always significantly outperform constant proportion strategies.