A global solution curve for a class of periodic problems,including the relativistic pendulum

Using continuation methods, we study the global solution structure of periodic solutions for a class of periodically forced equations, generalizing the case of relativistic pendulum. We obtain results on the existence and multiplicity of periodic solutions. Our approach is suitable for numerical computations, and in fact we present some numerically computed bifurcation diagrams illustrating our results.     

An inverse problem of Bilayer textile thickness determination in dynamic heat and moisture transfer

An inverse problem of bilayer textile thickness determination in dynamic heat and moisture transfer is presented satisfying the heat–moisture comfort level of human body. Heat and mass transfer law in bilayer textiles is displayed by proving the existence and uniqueness of solution to the coupled partial differential equations with initial-boundary value conditions. The finite difference method is employed to derive the numerical solution to partial differential equations. The regularized solution of the inverse problem is reformulated into solving function minimum problem through the Tikhonov regularization method. The golden section method is applied to solve the direct search problem and achieve the optimal solution to the inverse problem. Numerical algorithm and its numerical results provide theoretical explanation for textile materials research and development.     

On the multiple hyperbolic systems modelling phase transformation kinetics

We discuss Cahn’s time cone method modelling phase transformation kinetics. The model equation by the time cone method is an integral equation in the space-time region. First, we reduce it to a system of hyperbolic equations, and in the case of odd spatial dimensions, the reduced system is a multiple hyperbolic equation. Next, we propose a numerical method for such a hyperbolic system. By means of alternating direction implicit methods, numerical simulations for practical forward problems are implemented with satisfactory accuracy and efficiency. In particular, in the three dimensional case, our numerical method on the basis of reduced multiple hyperbolic equation is fast.     

On a general class of trigonometric functions and Fourier series†

We discuss a general class of trigonometric functions whose corresponding Fourier series can be used to calculate several interesting numerical series. Particular cases are presented.     

On some confidence intervals for estimating the mean of a skewed population

A number of methods are available in the literature to measure confidence intervals. Here, confidence intervals for estimating the population mean of a skewed distribution are considered. This note proposes two alternative confidence intervals, namely, Median t and Mad t, which are simple adjustments to the Student's t confidence interval. In order to compare the performance of these intervals, the following criteria are considered: (i) coverage probability; (ii) average width; and (iii) ratio of coverage to width. A simulation study has been undertaken to compare the performance of the intervals. The simulation study shows that for small sample size and moderate to highly skewed distributions, the proposed Median t performs the best in the sense of higher coverage, and the Mad t performs best in the sense of smaller confidence width. The proposed methods are very easy to calculate and are not overly computer-intensive, like Bootstrap confidence intervals. Some real-life examples have been considered that support the findings of the paper to some extent.     

Constructing Almost Symmetric Numerical Semigroups from Irreducible Numerical Semigroups

Every almost symmetric numerical semigroup can be constructed by removing some minimal generators from an irreducible numerical semigroup with its same Frobenius number.     

矩阵LU分解定理的简单证明及其数值实现

用简单的方法证明了矩阵LU分解定理,讨论了定理的推广以及定理相应的数值实现,并对《数值分析》课程教学方法改革进行了思考.     

Robin型离散Schwarz波形松弛算法的收敛性分析

Schwarz波形松弛(Schwarz waveform relaxation,SWR)是一种新型区域分解算法,是当今并行计算研究领域的焦点之一,但针对该算法的收敛性分析基本上都停留在时空连续层面.从实际计算角度看,分析离散SWR算法的收敛性更重要.本文考虑SWR研究领域中非常流行的Robin型人工边界条件,分析时空离散参数t和x、模型参数等因素对算法收敛速度的影响.Robin型人工边界条件中含有一个自由参数p,可以用来优化算法的收敛速度,但最优参数的选取却需要求解一个非常复杂的极小-极大问题.本文对该极小-极大问题进行深入分析,给出最优参数的计算方法.本文给出的数值实验结果表明所获最优参数具有以下优点:(1)相比连续情形下所获最优参数,利用离散情形下获得的参数可以进一步提高Robin型SWR算法在实际计算中的收敛速度,当固定t或x而令另一个趋于零时,利用离散情形下所获参数可以使算法的收敛速度具有鲁棒性(即收敛速度不随离散参数的减小而持续变慢).(2)相比连续情形下所获收敛速度估计,离散情形下获得的收敛速度估计可以更加准确地预测算法的实际收敛速度.     

CONSTITUTIVE MODELS AND HYPERGRAVITY PHYSICAL SIMULATION OF SOILS$^{\bf 1)}$

数值模拟和物理模拟是分析土体沉降和稳定性的主要手段. 本构模型作为描述土体应力应变关系的数学表达式, 是数值模拟的基础. 土体具有碎散性, 这一基本物理特性导致了其具有压硬性、摩擦性和剪胀性, 这是土的力学特性区别于金属的主要特征, 在土体的本构模型中必须反映这3个基本特性. 传统土力学将土体的变形和强度分离考虑, 分别采用弹性理论和基于刚塑性模型的极限平衡理论分析, 虽然应用广泛, 但由于不能全面地反映土的基本力学特性, 计算结果的精度常常难以满足定量分析的需要. 剑桥模型作为第一个全面反映压硬性、摩擦性和剪胀性的弹塑性本构模型, 实现了变形和强度的统一, 能较好地描述饱和正常固结黏土的应力应变关系, 被视为是现代土力学的开端; 统一硬化模型通过引入一个独特的硬化参数进一步发展了剑桥模型, 将适用范围扩大到超固结黏土. 作者认为, 未来岩土体本构模型研究的挑战是: 如何考虑岩土体在受力过程中土骨架相变与多场耦合, 以解决目前本构模型尚无法定量分析的能源、交通、环境和水利相关的重大岩土工程问题. 超重力物理模拟具有缩尺效应和缩时效应, 克服了常重力物理模拟中模型的应力水平低于原型的缺点, 特别适用于大尺度、长历时问题的模拟. 相较数值模拟, 超重力物理模拟的优势在于能够检验本构模型的合理性, 揭示本构模型无法描述的未知特性. 最后, 介绍了采用数值模拟和物理模拟联合分析大直径钢管桩水平受荷特性的工程案例.