Simulation Methods for Functional Materials

Functional materials embrace a broad area, ranging from functional information materials to special polymers, from special chemicals for printing to matenais used in making paints. Inasmuch as most tunctional materials are particulate, the present contribution is considered pertinent to the present FORUM.     

随机结构正交展开分析的Ritz动力聚缩法

针对随机结果正交展开理论计算上的弱点 ,本文在分析扩阶矩阵特性的基础上 ,于 Ritz模态向量子空间中对扩阶方程实现动力聚缩 ,大大提高了正交展开理论对实际工程问题的分析能力。分析实例表明 :即使结构参数具有很大变异性 (如δ =0 .4 )时 ,该算法依然能理想地与 Monte Carlo法模拟结果相吻合 ,计算时间则远远小于 Monte Carlo模拟法。同时 ,分析例证再一次强化了在结构动力分析中考虑结构参数随机性的必要性     

各向异性本构关系在板料成形数值模拟中的应用

对几种能表达面内各向异性的屈服准则Hill、Barlat-Lian、Barlat进行了比较。以弹性变形服从各向同性广义虎克定律的情况下，给出了基于张量算法推导的弹塑性本构关系的一般表达式，并由此导出了相应屈服准则的弹塑性本构关系的显式表达。借助ABAQUS软件本构模块用户子程序接口，分别实现了这些屈服准则在ABAQUS的嵌入。以模拟方形盒的拉延过程为例，分析了不同的屈服准则在板料成形过程数值模拟中的应用。模拟结果表明，基于弹塑性本构关系一般表达所列出的相应屈服准则的显式表达式是正确的；在采用壳元来模拟板料成形时，采用Barlat准则的模拟结果和采用Barlat-Lian准则的结果差别不大。     

延性材料平面冲击波基本实验数值模拟的若干问题

假设屈服强度正比于剪切模量,分别构建了四种强度模型,并用于数值模拟Ly12铝的平面冲击波试验的质点速度剖面以及钨合金的平面冲击波试验的锰铜计应力记录.研究发现,基于Steinberg等的屈服强度模型比较符合实验.还对平面冲击波试验的层裂现象进行了简化的数值模拟,其中,采用了该文作者所提出的一种基于空穴聚集的层裂模型,该模型的特点是分别引入了层裂过程的早期及后期应力松弛方程.     

L-leap-accelerating the stochastic simulation of chemically reacting systems

Presented here is an L-leap method for accelerating stochastic simulation of well-stirred chemically reacting systems,in which the number of reactions occurring in a reaction channel with the largest propensity function is calculated from the leap condition and the number of reactions occurring in the other reaction channels are generated by using binomial random variables during a leap.The L-leap method can better satisfy the leap condition.Numerical simulation results indicate that the L-leap method can obtain better performance than established methods.     

促进其线性频散特征另一种形式的Bousinesq方程

Ｂｏｕｓｉｎｅｓｑ方程能够用于模拟表面重力波传播过程中的折射、绕射、反射以及浅化，非线性作用等现象．用不同垂直积分方法所得到的二维Ｂｏｕｓｓｉｎｅｓｑ方程形式具有不同的线性频散特征．采用两个不同的水深层的水平速度变量组合，推导出一个新形式的Ｂｏｕｓｉｎｅｓｑ方程．通过对其参数的设置可得到精确的线性频散解Ｐａｄｅ近似４阶精度．其适用范围已由原来的浅水，向深水拓进．相速误差小于２％，其拓展适用范围可达到０８个波长水深．应用所得到的新型Ｂｏｕｓｉｎｅｓｑ方程，采用有限差分法，对经典工况进行了数值模拟，其计算结果表明，计算值与物模实验值吻合较好．这说明本文新形式的Ｂｏｕｓｓｉｎｅｓｑ方程对变水深非线性效应所产生的能量频散有着较为精确的描述     

A quasi-flow constitutive model with strain-rate dependence

In this paper, the proposed is a quasi-flow constitutive model with strain-rate sensitivity for elastic plastic large deformation. The model is based on the Quasi-flow Corner theory, and is suitable for the sheet metal forming process simulation with a variable punch machine velocity. Uniaxial tensile tests and deep-drawing tests of a circular blank with square punch are carried out and numerically simulated. The consistency between the experimental and the numerically simulated results shows the validity of the present new constitutive model. The project supported by the Scientific Foundation of National Outstanding Youth of China (10125208), the National Natural Science Foundation of China (19832020), and the National Education Committee of China     

自紧身管临界裂纹尺寸的概率断裂力学研究

基于概率断裂力学理论和Mome Caflo模拟方法，本文进行了自紧身管临界裂纹尺寸的可靠性研究。自紧身管内表面的疲劳裂纹考虑为半椭圆形式。裂纹尖端处的应力强度因子由内压和自紧残余应力共同产生。自紧残余应力采用了符合身管材料具有强化和包辛格效应性能推导的公式计算，它产生的应力强度因子通过权函数方法得到。根据断裂准则，可计算出自紧身管的临界裂纹尺寸。实例分析表明，对数正态分布为临界裂纹尺寸的最佳分布，同时给出了在各种置信度和可靠度下自紧身管的临界裂纹尺寸。     

Melnikov方法在输流管混沌运动研究中的应用

对基础简谐运动激励下两端固定输流管道的混沌运动进行了研究，推导出了系统的运动方程，确定了系统存在的平衡点及其稳定性，计算出了未扰系统的同宿轨道，并利用Melnikov方法得到了系统发生混沌运动时参数需满足的临界条件，同时还利用相平面图和：Poincare映射等方法对管道的混沌运动进行了数值模拟，通过比较发现，由Melnikov方法确定的临界参数值要稍小于数值模拟中首次观察到混沌运动时对应的临界值。     

Theoretical statistical solution and numerical simulation of heterogeneous brittle materials

The analytical stress-strain relation with heterogeneous parameters is derived for theheterogeneous brittle materials under a uniaxial extensional load, in which the distributions of theelastic modulus and the failure strength are assumed to be statistically independent. This theoreticalsolution gives an approximate estimate of the equivalent stress-strain relations for 3-D heterogeneousmaterials. In one-dimensional cases it may provide comparatively accurate results. The theoreticalsolution can help us to explain how the heterogeneity influences the mechanical behaviors, Further, anumerical approach is developed to model the non-linear behavior of three-dimensional heterogeneousbrittle materials. The lattice approach and statistical techniques are applied to simulate the initialheterogeneity of heterogeneous materials. The load increment in each loading stage is adaptivelydetermined so that the better approximation of the failure process can be realized. When the maximumtensile principal strain exceeds the failure strain, the elements are considered to be broken, which canbe carried out by replacing its Young‘s modulus with a very small value. A 3-D heterogeneous brittlematerial specimen is simulated during a full failure process. The numerical results are in good agreementwith the analytical solutions and experimental data.