Orientation-based Continuum Damage Models for Rocks

A general formulation of the Helmholtz free energy used in thermodynamics of damage process of rocks is derived within a multi-scale framework. Such a physically-based thermodynamic state potential has a hybrid, discrete/continuum, nature in the sense tha

Molecular/Cluster Statistical Thermodynamics Methods To Simulate Quasi-Static Deformations At Finite Temperature

The rapid evolution of nanotechnology appeals for the understanding of global response of nanoscale systems based on atomic interactions, hence necessitates novel, sophisticated, and physically based approaches to bridge the gaps between various length and time scales. In this paper, we propose a group of statistical thermodynamics methods for the simulations of nanoscale systems under quasi-static loading at finite temperature, that is, molecular statistical thermodynamics (MST) method, cluster statistical thermodynamics (CST) method, and the hybrid molecular/cluster statistical thermodynamics (HMCST) method. These methods, by treating atoms as oscillators and particles simultaneously, as well as clusters, comprise different spatial and temporal scales in a unified framework. One appealing feature of these methods is their "seamlessness" or consistency in the same underlying atomistic model in all regions consisting of atoms and clusters, and hence can avoid the ghost force in the simulation. On the other hand, compared with conventional MD simulations, their high computational efficiency appears very attractive, as manifested by the simulations of uniaxial compression and nanoindenation. (C) 2008 Elsevier Ltd. All rights reserved.

simulation and interaction of fluid dynamics

In the fluid simulation, the fluids and their surroundings may greatly change properties such as shape and temperature simultaneously, and different surroundings would characterize different interactions, which would change the shape and motion of the fluids in different ways. On the other hand, interactions among fluid mixtures of different kinds would generate more comprehensive behavior. To investigate the interaction behavior in physically based simulation of fluids, it is of importance to build physically correct models to represent the varying interactions between fluids and the environments, as well as interactions among the mixtures. In this paper, we will make a simple review of the interactions, and focus on those most interesting to us, and model them with various physical solutions. In particular, more detail will be given on the simulation of miscible and immiscible binary mixtures. In some of the methods, it is advantageous to be taken with the graphics processing unit (GPU) to achieve real-time computation for middle-scale simulation.

CLIGEN非降水要素在黄土塬区的适用性评估

许多基于物理机制的水文和作物模型需要日序列气象数据来驱动,CLIGEN是为WEPP等模型产生气候输入文件的天气发生器,可以产生10个日序列气象变量来满足这种需要,但是其在中国的适用性需要进行评估。研究的目标是利用黄土高原陕西长武1957~2001年的气象数据评估CLIGEN产生非降水要素(最高温度、最低温度、露点温度、太阳辐射和风速)的能力。结果表明,CLIGEN对最高温度、最低温度和露点温度的模拟效果较好,对太阳辐射和极端气候事件的模拟效果较差,对风速的模拟效果最差。相关性检验表明CLIGEN很好地保持了气象要素的季节性,这对模拟农业生产是非常重要的;但是没有保留气象要素逐日的自相关和互相关性,进而导致产生的温度变化不符合连续渐变的规律。

基于物理的植物叶子形态变化过程仿真造型

文中提供了一种对植物叶子造型的新方法,可以非常真实地展现叶子在干枯、老化过程中几何形态的变化.基于植物学和物理原理,文中引入了双层结构模型表达叶子的力学结构,很好地模拟了不同种类植物叶子形状的多样性.文中首先分析了叶子形状变化的原因,即在叶子枯萎过程中,叶肉和叶脉由于各自不同的组织结构,导致了收缩比例不同,从而在物理模拟中,通过建立关于叶子基本结构的双层质点-弹簧模型,并对上下两层的不同参数的合理设置,很好地表现了叶肉和叶脉的不同力学特性.双层模型的相互作用,决定了叶子最终变形的方式和效果,由此可以得到非常接近真实树叶的各种叶子形态.

可重构星球探测机器人的OpenGL仿真实验平台研究

采用基于Visual C++和OpenGL的建模和运动仿真方法,对可重构星球探测机器人系统的三维运动仿真实验平台进行了研究,建立了一个多机器人系统的仿真实验平台。开发的实验平台可用于探索和验证机器人系统的工作原理、工作空间、多机器人协调算法、重构方法、系统集成技术等。在该平台上进行了机器人的运动学仿真和协调运动研究,验证了该仿真平台的有效性和机器人系统体系结构的合理性。

地震波模拟：非规则网格离散的数值解法、吸收边界及其实际应用

This dissertation presents a series of irregular-grid based numerical technique for modeling seismic wave propagation in heterogeneous media. The study involves the generation of the irregular numerical mesh corresponding to the irregular grid scheme, the discretized version of motion equations under the unstructured mesh, and irregular-grid absorbing boundary conditions. The resulting numerical technique has been used in generating the synthetic data sets on the realistic complex geologic models that can examine the migration schemes. The motion equation discretization and modeling are based on Grid Method. The key idea is to use the integral equilibrium principle to replace the operator at each grid in Finite Difference scheme and variational formulation in Finite Element Method. The irregular grids of complex geologic model is generated by the Paving Method, which allow varying grid spacing according to meshing constraints. The grids have great quality at domain boundaries and contain equal quantities of nodes at interfaces, which avoids the interpolation of parameters and variables. The irregular grid absorbing boundary conditions is developed by extending the Perfectly Matched Layer method to the rotated local coordinates. The splitted PML equations of the first-order system is derived by using integral equilibrium principle. The proposed scheme can build PML boundary of arbitrary geometry in the computational domain, avoiding the special treatment at corners in a standard PML method and saving considerable memory and computation cost. The numerical implementation demonstrates the desired qualities of irregular grid based modeling technique. In particular, (1) smaller memory requirements and computational time are needed by changing the grid spacing according to local velocity; (2) Arbitrary surfaces and interface topographies are described accurately, thus removing the artificial reflection resulting from the stair approximation of the curved or dipping interfaces; (3) computational domain is significantly reduced by flexibly building the curved artificial boundaries using the irregular-grid absorbing boundary conditions. The proposed irregular grid approach is apply to reverse time migration as the extrapolation algorithm. It can discretize the smoothed velocity model by irregular grid of variable scale, which contributes to reduce the computation cost. The topography. It can also handle data set of arbitrary topography and no field correction is needed.

A simple method to simulate shrinkage-induced cracking in cement-based composites by lattice-type modeling

A new numerical procedure is proposed to investigate cracking behaviors induced by mismatch between the matrix phase and aggregates due to matrix shrinkage in cement-based composites. This kind of failure processes is simplified in this investigation as a purely spontaneous mechanical problem, therefore, one main difficulty during simulating the phenomenon lies that no explicit external load serves as the drive to propel development of this physical process. As a result, it is different from classical mechanical problems and seems hard to be solved by using directly the classical finite element method (FEM), a typical kind of "load -> medium -> response" procedures. As a solution, the actual mismatch deformation field is decomposed into two virtual fields, both of which can be obtained by the classical FEM. Then the actual response is obtained by adding together the two virtual displacement fields based on the principle of superposition. Then, critical elements are detected successively by the event-by-event technique. The micro-structure of composites is implemented by employing the generalized beam (GB) lattice model. Numerical examples are given to show the effectiveness of the method, and detailed discussions are conducted on influences of material properties.

AMPS modeling of light J-V characteristics of a-Si based solar cells

AMPS (Analysis of microelectronic and photonic structures) mode,which was developed by Pennsylvania State University, has been used to module the light J-V characteristics of a-Si solar cells with a structure of TCO/p-a-SiC:H/i-a-Si:H/n-a-Si:H/ metal. The effects of valence band offset and contact barriers at p/i and TOC/p, n/metal interfaces on the light J-V characteristics have been examined. The modeling has qualitatively categorized and explained the non-ideal J-V behaviors (rollover, crossover, Voc shift,and rollunder) observed in a-Si based solar cells.

stochastic process algebra based software process simulation modeling

Chinese Acad Sci, ISCAS Lab Internet Software Technologies

Modeling and mesoscopic damage constitutive relation of brittle short-fiber-reinforced composites

Aimed at brittle composites reinforced by randomly distributed short-fibers with a relatively large aspect ratio, stiffness modulus and strength, a mesoscopic material model was proposed. Based on the statistical description, damage mechanisms, damage-induced anisotropy, damage rate effect and stress redistribution, the constitutive relation were derived. By taking glass fiber reinforced polypropylene polymers as an example, the effect of initial orientation distribution of fibers, damage-induced anisotropy, and damage-rate effect on macro-behaviors of composites were quantitatively analyzed. The theoretical predictions compared favorably with the experimental results.