Finite element formulation for dynamics of planar flexible multi-beam system

In some previous geometric nonlinear finite element formulations, due to the use of axial displacement, the contribution of all the elements lying between the reference node of zero axial displacement and the element to the foreshortening effect should be taken into account. In this paper, a finite element formulation is proposed based on geometric nonlinear elastic theory and finite element technique. The coupling deformation terms of an arbitrary point only relate to the nodal coordinates of the element at which the point is located. Based on Hamilton principle, dynamic equations of elastic beams undergoing large overall motions are derived. To investigate the effect of coupling deformation terms on system dynamic characters and reduce the dynamic equations, a complete dynamic model and three reduced models of hub-beam are prospected. When the Cartesian deformation coordinates are adopted, the results indicate that the terms related to the coupling deformation in the inertia forces of dynamic equations have small effect on system dynamic behavior and may be neglected, whereas the terms related to coupling deformation in the elastic forces are important for system dynamic behavior and should be considered in dynamic equation. Numerical examples of the rotating beam and flexible beam system are carried out to demonstrate the accuracy and validity of this dynamic model. Furthermore, it is shown that a small number of finite elements are needed to obtain a stable solution using the present coupling finite element formulation.     

Geometric nonlinear effects on the planar dynamics of a pivoted flexible beam encountering a point-surface impact

Flexible-body modeling with geometric nonlinearities remains a hot topic of research by applications in multibody system dynamics undergoing large overall motions. However, the geometric nonlinear effects on the impact dynamics of flexible multibody systems have attracted significantly less attention. In this paper, a point-surface impact problem between a rigid ball and a pivoted flexible beam is investigated. The Hertzian contact law is used to describe the impact process, and the dynamic equations are formulated in the floating frame of reference using the assumed mode method. The two important geometric nonlinear effects of the flexible beam are taken into account, i.e., the longitudinal foreshortening effect due to the transverse deformation, and the stress stiffness effect due to the axial force. The simulation results show that good consistency can be obtained with the nonlinear finite element program ABAQUS/Explicit if proper geometric nonlinearities are included in the floating frame formulation. Specifically, only the foreshortening effect should be considered in a pure transverse impact for efficiency, while the stress stiffness effect should be further considered in an oblique case with much more computational effort. It also implies that the geometric nonlinear effects should be considered properly in the impact dynamic analysis of more general flexible multibody systems.     

高阶旋转体时域有限差分法

旋转体时域有限差分法(BOR FDTD)是模拟电磁波与旋转对称目标体相互作用的一种高效算法,但目前的BOR FDTD对时间和空间的差分精度均只有二阶.通过把电磁场各分量作泰勒展开,并结合麦克斯韦方程,提出了BOR FDTD的三种高阶算法.相对于目前的BOR FDTD,三种高阶算法均不增加存储量.通过模拟简谐波与介质球的相互作用并和理论值作比较,证实了高阶算法在提高计算精度方面的有效性.此外还分析了三种高阶算法的优、缺点,给出了关于高阶算法应用的建议.     

有限差分法的并行化计算实现

有限差分法是求解偏微分方程近似解的一种重要的数值方法。并行化计算可提高复杂计算问题的效率,二维场中拉普拉斯方程的差分格式非常适合并行化方法的计算。如何将串行部分并行化以提高大规模计算的效率,MPI（消息传递接口）是实现并行程序设计的标准之一。虚拟进程（MPI_PROC_NULL）是MPI中的假想进程,它的引用可简化MPI编程中的通信部分,引入虚拟进程编写代码,可实现有限差分方法的并行化计算。     

Combined analytical solution of overland flow and sediment transport

With reference to the kinematic wave theory coupled with the hypothesis of constant linear velocity for the rating curve, rising limb analytical solutions have been calculated for overland flow, over an Hortonian-infiltrating surface, and sediment discharge. These analytical solutions are certainly easier to use than the numerical integration of the basic equations and they may be used to obtain an initial evaluation of the parameters of more complex models generally devised for complicated cases.Notation a exponent of the Horton law [T^–1] - b exponent of the rill erosion equation - B inter-rill erosion coefficient [ML^–m–2T ^m–1] - c sediment concentration [ML^–3] - c _o reference sediment concentration [ML^–3] - E _I inter-rill erosion [ML^–2T^–1] - E _R rill erosion [ML^–2T^–1] - f _c final infiltration rate of the soil [LT^–1] - f _o initial infiltration rate of the soil [LT^–1] - h flow depth [L] - h _o reference flow depth [L] - i infiltration rate [LT^–1] - k rill erosion coefficient [ML^–1–b T^–1] - K integration constant - L() Laplace transformation - m exponent of the inter-rill erosion equation - n Manning's coefficient [L^–1/3T] - p rainfall intensity [LT^–1] - q water discharge per unit width [L²T^–1] - q _s sediment discharge per unit width [ML^–1T^–1] - t time [T] - t _p ponding time [T] - x distance along the flow direction [L] Greek Letters coefficient of the stage-discharge equation [L^2–T^–1] - exponent of the stage-discharge equation - rill erosion coefficient [L^–1]     

黑磷纳米盘-层等离激元系统中的各向异性可调多阶强耦合

黑磷支持各向异性的表面等离激元,可用于设计具备更多功能的原理性器件。用时域有限差分法数值模拟了中红外到远红外波段基于黑磷的层-盘-层系统中不同等离激元模式之间的杂交行为。通过动态调节黑磷中的载流子浓度,可以实现两个晶格方向上强耦合现象的产生与控制。对不同模式间的耦合进行分析并计算,得到吸收光谱中的拉比分裂能最高可达42.9 meV。此外,还计算了偏振角度对各向异性强耦合的影响,其最高可以实现6个吸收频带。该模型可为构建未来在中远红外波段工作的基于二维材料紧凑型各向异性等离激元器件提供基础。     

Recent Advances on 3D Printing Technique for Thermal‐Related Applications

Advances in ink formulation and printing techniques make producing material systems with new and versatile characteristics and functionalities possible. Additive manufacturing or 3D printing enables fabricating complex structures at a faster production rate using different types of materials for various applications. Recently, 3D printing methods are being studied for thermal‐related applications. In this paper, the authors review recent progress of materials and printing techniques for thermal application devices using composite materials.

     

Min-degree constrained minimum spanning tree problem: New formulation via Miller–Tucker–Zemlin constraints

Given an undirected network with positive edge costs and a positive integer $d>2$, the minimum-degree constrained minimum spanning tree problem is the problem of finding a spanning tree with minimum total cost such that each non-leaf node in the tree has a degree of at least $d$. This problem is new to the literature while the related problem with upper bound constraints on degrees is well studied. Mixed-integer programs proposed for either type of problem is composed, in general, of a tree-defining part and a degree-enforcing part. In our formulation of the minimum-degree constrained minimum spanning tree problem, the tree-defining part is based on the Miller–Tucker–Zemlin constraints while the only earlier paper available in the literature on this problem uses single and multi-commodity flow-based formulations that are well studied for the case of upper degree constraints. We propose a new set of constraints for the degree-enforcing part that lead to significantly better solution times than earlier approaches when used in conjunction with Miller–Tucker–Zemlin constraints.     

The stochastic -hub center problem with service-level constraints

Small package delivery companies offer services where packages are guaranteed to be delivered within a given time-frame. With variability in travel time, the configuration on the hub-and-spoke delivery network is vital in ensuring a high probability of meeting the service-level guarantee. We present the stochastic p-hub center problem with chance constraints, which we use to model the service-level guarantees. We discuss analytical results, propose solution heuristics, and present the results from computational experiments.     

A two-asset stochastic model for long-term portfolio selection

In mean–variance (M–V) analysis, an investor with a holding period [0,T] operates in a two-dimensional space—one is the mean and the other is the variance. At time 0, he/she evaluates alternative portfolios based on their means and variances, and holds a combination of the market portfolio (e.g., an index fund) and the risk-free asset to maximize his/her expected utility at time T. In our continuous-time model, we operate in a three-dimensional space—the first is the spot rate, the second is the expected return on the risky asset (e.g., an index fund), and the third is time. At various times over [0,T], we determine, for each combination of the spot rate and expected return, the optimum fractions invested in the risky and risk-free assets to maximize our expected utility at time T. Hence, unlike those static M–V models, our dynamic model allows investors to trade at any time in response to changes in the market conditions and the length of their holding period. Our results show that (1) the optimum fraction y*(t) in the risky asset increases as the expected return increases but decreases as the spot rate increases; (2) y*(t) decreases as the holding period shortens; and (3) y*(t) decreases as the risk aversion parameter-γ is larger.