Product integration methods for second-kind Abel integral equations

The construction and convergence of high-order product integration methods for the second-kind Abel equation are discussed and the results of De Hoog and Weiss are generalised. Backward difference methods are introduced, and numerical results are presented which verify the theoretical rates of convergence.     

NUMERICAL SIMULATION OF DEFORMATION BEHAVIOR OF 22MnB5 BORON STEEL AT ELEVATED TEMP ERATURES AND EXPERIMENTAL VERIFICATION

The effects of strain, strain rate and temperature on the mechanical behavior of 22MnB5 boron steel deformed isothermally under uniaxial tension tests and the experimental characterization of 22MnB5 boron steel in the austenitic region have been investigated. Based on the crystal plasticity theory and thermal kinematics, an improved integration model is presented. In this model, the elastic deformation gradient is the integration variable of the governing equation, which contains not only the elastic deformation but also the thermal effects. In the coupled thermo- mechanical process, this model can reveal the evolution of microstructures such as the rotation of a single crystal and the slip systems in each of them. The plastic behavior of the boron steel can be well described by the presented model.     

AN ADAPTIVE UNSTRUCTURED TRI-TREE ITERATIVE SOLVER FOR MIXED FINITE ELEMENT FORMULATION OF THE STOKES EQUATIONS

An iterative adaptive equation solver for solving the implicit Stokes equations simultaneously with tri-tree grid generation is developed. The tri-tree grid generator builds a hierarchical grid structure which is mapped to a finite element grid at each hierarchical level. For each hierarchical finite element grid the Stokes equations are solved. The approximate solution at each level is projected onto the next finer grid and used as a start vector for the iterative equation solver at the finer level. When the finest grid is reached, the equation solver is iterated until a tolerated solution is reached. In order to reduce the overall work, the element matrices are integrated analytically beforehand. The efficiency and behaviour of the present adaptive method are compared with those of the previously developed iterative equation solver which is preconditioned by incomplete LU factorization with coupled node fill-in. The efficiency of the incomplete coupled node fill-in preconditioner is shown to be largely dependent on the global node numbering. The preconditioner is therefore tested for the natural node ordering of the tri-tree grid generator and for different ways of sorting the nodes.     

An explicit formulation for multiscale modeling of bcc metals

Many materials for specialized applications exhibit a body-centered cubic structure; e.g., tantalum, vanadium, barium and chromium. In addition, the successful modeling of body-centered cubic (bcc) metals is a necessary step toward modeling of common structural materials such as iron. Implicit formulations for this class of materials exist [e.g., Stainier, L., Cuitiño, A., Ortiz, M., 2002. A micromechanical model of hardening, rate sensitivity, and thermal softening in bcc crystals. Journal of the Mechanics and Physics of Solids 50 (7), 1511–1545; Kuchnicki, S., Radovitzky, R., Cuitiño, A., Strachan, A., Ortiz, M., 2007. A pressure-dependent multiscale model for bcc metals], but are impractical to resolve large-scale dynamic deformation processes. In this article, we describe a procedure analogous to Kuchnicki et al. [Kuchnicki, S., Cuitiño, A., Radovitzky, R., 2006. Efficient and robust constitutive integrators for single-crystal plasticity modeling. International Journal of Plasticity 22 (10), 1988–2011]. wherein we construct an explicit formulation for the multiscale physics models. This update is based on the model of Kuchnicki et al. (in preparation) using a power law representation for the plastic slip rates. The existing implicit form of the model provides qualitative matching with experiments at quasi-static strain rates. The model is recast in an explicit form and applied first to a high quasi-static strain rate to verify that the two forms of the model return similar predictions for similar input parameters. The explicit model is also applied to several high strain rates, showing that it captures characteristic features observed in experimental tests of high-rate deformations, such as the drop in stress immediately after yield that is present in split Hopkinson pressure bar (SHPB) experiments. This test provides qualitative evidence that the model is suitable for high-strain-rate applications. The utility of the model is further demonstrated by a one-dimensional simulation of a SHPB test. Finally, a test case modeling pressure impact of a Tantalum plate using 600,000 elements is shown. The simulations show that the explicit model is capable of recovering the salient features of the experiments while integrating the constitutive update in a robust manner.     

关于边界元方法中柯西主值积分的探讨

本文对边界元方法中的各类积分根据其奇异性作分类，并对主值积分的收敛条件、变量替换等进行了讨论，又给出了变替换附加项显式。文中提供的主值积分配项消奇术在边界元方法中是有普遍意义的。     

车-桥系统耦合振动响应的简便计算

依据振动理论推导出了二自由度模型车辆与桥梁系统竖向耦合振动微分方程，采用模态分析的离散化方法，将复杂的偏微分方程问题转化为变系数常微分方程问题，并将微分方程数值积分的Runge-Kutta方法引入到该时变系统的振动响应计算中，使复杂的耦合响应问题得到简便的解决。通过算例验证了该方法的有效性和简便性。该方法只需要直接数值积分，具有公式简单，编程方便，计算速度快等优点，特别适合于工程实际问题的计算，并且不仅适用于匀速运动车辆，也适用于变速运动车辆。     

Correct Solutions of Plane Elastic Problems for a Half-Plane

Continuous and integrable solutions and one-to-one relationships between boundary forces and displacements are found through the direct integration of the differential equations of the plane elastic problem for a half-plane with boundary conditions for either forces or displacements or with mixed boundary conditions. The necessary equilibrium conditions for forces and the compatibility conditions for displacements that ensure the correctness of the solutions are formulated     

基于加性四元数误差方程的惯性/天文姿态组合算法

针对采用旋转四元数误差进行的组合导航误差建模中状态方程的非线性化问题,提出一种新的惯性/天文组合姿态组合算法,以姿态加性四元数误差和陀螺漂移为状态变量,推导系统线性化状态方程,并以天文导航和惯导姿态四元数之差为量测量,建立系统量测方程,然后利用卡尔曼滤波实现对该组合模式的信息融合,仿真分析表明,所设计的基于姿态四元数误差和陀螺漂移的组合模式能够有效估计系统状态误差,姿态误差0.02°左右,验证了其有效性,可避免较为复杂的非线性滤波器的使用,为工程实践提供了理论支持.     

求解二点边值问题打靶法的一种改进方法

本文讨论两点边值问题的数值解法,利用最小二乘法修正打靶法所需初始参数,将边值问题转化为相应的初值问题与左梯度控制方程合并,引入精细积分法进行求解。通过数值算例表明,该方法是一种求解两点边值问题的有效方法。