讨论多尺度模拟断裂

jianqiwei

多尺度，模拟断裂，这个应该是材料区的啊，为什么是数学呢？

zzmmxx

学习当中,希望大家有好的ANSYS文章,拿出来共享一下!!!

ydq23

多尺度断裂模拟在桥梁结构分析方面今年来也有应用，比如江苏的润扬长江公路大桥的钢箱梁结构分析就用到多尺度断裂模拟！！

feiyuzhen

学习一下

haows520

感觉粘性破裂好像也属于损伤，而脆性的多属于裂纹，微观是宏观的观测基础。

rock.li

多尺度灌水

ansys2007

一种模型势只能适用于一种原子,不能满足实际模拟的需要。

敦诚

找到一点关于损伤力学的介绍，是关于梯度依赖的，也算是多尺度的相关介绍：
混凝土类材料与结构的破坏失效是以材料内部裂纹萌生、扩展，并最终导致结构断裂破坏失去承载能力为特征的。固体力学的弹塑性与损伤耦合模型是目前较为常用的处理上述问题的理论工具。材料损伤演化的结果是非弹性变形的局部化，即：非弹性变形在变形后期随着外载增加而向物体的某一个局部集中，而物体的其它部分则处于弹性卸载状态。具体到混凝土类材料，这个过程为：原来成片的微裂纹区最后汇集成一条主要的宏观裂纹，形成结构断裂，宏观裂纹外的微裂纹随之闭合。
在数值计算方面，应用现有局部理论计算时，在载荷极点后区的数值结果通常严重依赖网格划分。因此，有必要采用非简单材料的非局部理论，以消除数值分析结果的网格依赖性，同时计入裂纹间的相互影响。另一方面，连续介质的非局部理论认为：物体内一个点的非弹性力学行为不仅取决于外载荷作用下该点的力学量的状态，而且与该点邻域内的点的力学量有关。该理论还认为：变形的一阶梯度不足以描述复杂材料的运动学及材料的力学行为。1909年，Cosserat首先提出了著名的Cosserat本构理论。其后，二十世纪六十年代，在高阶变形梯度等方面有了进一步的发展。限于求解手段，二十世纪八十年代以前，上述理论一直仅限于材料的理论研究，在固体力学结构分析中没有重大应用实例的报道。
二十世纪八十年代以来，非局部本构理论重新受到了重视。目前已有的理论模型除了上述的Cosserat模型外，主要的还有两大类：（1）“面积加权平均”的非局部非弹性模型。在这一类模型中，引入了材料的特征尺寸的概念，同时引入权函数来表示物质点间力学量相互影响的数学关系，从而确定了局部化区域的宽度，并有效地消除了数值结果的网格依赖性。（2）“梯度依赖”的非局部模型。这一类模型的理论及应用研究是目前国际上最为活跃的本构研究课题。由于数学上并没有严格的理论来规范本构行为中梯度依赖项的细节（阶数及形式等），目前的文献中有多种不同的梯度依赖形式。其中主要的有：（1）金属塑性研究的应变梯度理论。文[8]对国内外金属塑性梯度依赖模型进行了详细的综述， 文中则分别有详细的模型介绍和新进展报告。这一类模型主要包括高阶应变张量模型、偶应力理论等模型。以高阶应变模型为例，这一类模型的缺点是：由于3阶应变张量的引入，数值求解时需要引入特殊单元，计算量十分大，在国际上，未见有该类模型在混凝土的损伤塑性分析中应用的例子。（2）Aifantis的梯度依赖模型。Aifantis和他的合作者首先提出了一种应变软化塑性的梯度依赖模型：将等效塑性应变的Laplacian项和特征长度一起引入到塑性屈服条件中，从而使得物体内一点的弹塑性状态不仅取决于传统的等效应力和材料的屈服极限，还有赖于该点与相邻物质点之间的塑性应变场的Laplacian值和特征长度值的积。这一类模型数学表达简单，但是由于其梯度依赖项没有相应的热力学共轭量，因此在理论上是不完善的。可尽管其热力学理论不完善，但由于其数学表达简洁、实践上实用有效，Aifantis类的梯度依赖模型依然在混凝土的损伤塑性分析中得到了广泛应用，例如：DeBorst和他的合作者们将Aifantis的梯度依赖模型应用到混凝土标量损伤与塑性耦合的模型中，建立了梯度依赖的损伤塑性耦合模型，并进一步结合面积加权平均非局部理论的思想，提出了隐式的梯度依赖模型。
除了上述两类主要的梯度依赖塑性损伤本构模型之外，近年来一些研究者开始从热力学原理出发，提出混凝土的新的梯度依赖非局部本构模型：Nedjar将标量损伤的一阶梯度视作新的独立的内变量，推导了相应的虚功率原理的形式，给出了相应的计算格式。Liebe和Steinmann和Benallal一起，本着“能量耗散的非局部化”的理念，提出了“损伤一阶梯度依赖”的自由能表示式。潘一山和徐秉业等结合Aifantis的梯度依赖模型研究了岩石塑性屈服的尺寸效应问题。由于损伤和塑性是本质上不同的两类内变量，用塑性理论近似代替断裂损伤模型虽然有一定的效果，但是缺陷是明显的，目前这一做法在国际上愈来愈少见。关于特征长度值的定义和计算取值目前还没有普遍接受的定论。主要争议是：1)面积加权非局部模型定义特征长度为表示体元的尺寸，而梯度依赖模型则认为特征长度是非弹性行为相互影响的物质点间的最大距离。有少数人认为特征长度值是与混凝土中砂石颗粒尺寸相关的裂纹区宽度尺寸。另有人认为特征长度值只是梯度项对加载条件的影响系数，并且实际上在数值计算时有人也是这么做的；2)尽管多数人在计算中将特征长度值取为常量，但是这种做法并没有实验支持，而且近来Geers和deBorst等还提出了“时变型”的特征长度的理论模型，认为特征长度的值随着裂纹的扩展是变化的。
按照‘能量耗散非局部化’的原理，一个点上的非弹性过程的演化不仅与该点上的自由能有关，而且还受领域内物质点上的自由能影响。因此，本文认为梯度依赖损伤模型的特征长度应该是损伤演化能够相互影响的物质点间的最大距离，即损伤过程区的长度尺寸。同时认为特征长度值的平方（或二次函数）是‘梯度依赖项’对塑性损伤加载条件的影响系数。但是混凝土的力学行为受颗粒组分及试件尺寸等多种因素的影响，目前没有成熟的关于特征长度取值的理论模型。
**注，文【8】是指陈少华的那本书。

敦诚

以前参与过混凝土损伤梯度依赖理论方面的课题，但是当时不知道这就是多尺度问题，个人观点认为，由于损伤力研究的是亚微观现象，无法跟踪微裂纹这种微观现象，所以对这一发面心有怀疑。但是看到了一些文献，也做了一些数值计算，发现一些梯度依赖效应真实地存在在问题中，同时利用一些相关理论可以解决这些问题，才减少了对这一理论的猜测。

piteqiu

对这些问题没有研究过。对断裂的认识局限于做过几次试样拉伸断裂方面的实验。不过用到的理论都比较浅显

tt06601

第一次接触到多尺度（multi_scale）是在学习有限元重合网格法（s_version Fem）时，因有限元重合网格法思想中有全局和局部模型的概念，而局部可以相对全局在建模以及网格划分时有不同数量级的尺寸，后来老师说这种不是真正的多尺度，非要有个多少数量级，或者简单的说就是从宏观和微观的不同尺度研究的一种概念，我自己理解，这种从宏微观尺寸考虑的多尺度是否只是一种侠义的多尺度，广义的多尺度是否包含时间，位移等多种物理量上的一种多尺度，呵呵，还望各位多多指教！

tt06601

本人最近在做一个裂纹扩展的模拟，对于断裂的模拟，可能不是想象中那么简单，可以想到的困难，包括断裂机理如何，裂纹扩展的过程如何，以及模拟方法上，如何模拟裂纹，如何模拟裂纹扩展或轨迹，又涉及到不同方法下更新有限元网格时的效率问题等。

SUNSEAHE

来学习的。

ma

微尺度断裂问题，用分子动力学模拟会更加可靠，这在纳米力学领域是一个可靠的方向。

acamus

通过宏细观尺度的变化才能建立材料的宏观特性与细观结构之间的关系，就是通常所说的细观方法。还是由一定的难度

Lee2002hu

原帖由 wxgreat 于 2008-5-7 13:54 发表
对分子动力学模拟比较感兴趣，想把它应用到生物学领域，但是不知从何着手啊

直接使用gromacs就行，现在这个应用很广泛了。/。。
http://www.gromacs.org/
GROMACS is a versatile package to perform molecular dynamics, i.e. simulate the Newtonian equations of motion for systems with hundreds to millions of particles.
It is primarily designed for biochemical molecules like proteins and lipids that have a lot of complicated bonded interactions, but since GROMACS is extremely fast at calculating the nonbonded interactions (that usually dominate simulations) many groups are also using it for research on non-biological systems, e.g. polymers.
GROMACS supports all the usual algorithms you expect from a modern molecular dynamics implementation, (check the online reference or manual for details), but there are also quite a few features that make it stand out from the competition:

• GROMACS provides extremely high performance compared to all other programs. A lot of algorithmic optimizations have been introduced in the code; we have for instance extracted the calculation of the virial from the innermost loops over pairwise interactions, and we use our own software routines to calculate the inverse square root. The innermost loops are generated automatically in either C or Fortran at compile time, with optimizations adopted to your architecture. Assembly loops using SSE and 3DNow! multimedia instructions are provided for i386 processors, separate versions using all x86-64 registers are used on Opteron x86-64 and Xeon EM64t machines. This results in exceptional performance on inexpensive PC workstations, and for Pentium IV/Opteron processors there are also SSE2 double precision assembly loops. There are new manually tuned assembly loops for ia64 (both single and double precision), and last but certainly not least we have written Altivec assembly loops both for Linux and Mac OS X. Gromacs is normally 3-10 timesfaster than any other program; check the article in Journal of Molecular Modeling (reference can be found under resources) for a comparison benchmark.
• GROMACS is user-friendly, with topologies and parameter files written in clear text format. There is a lot of consistency checking, and clear error messages are issued when something is wrong. Since the C preprocessor is used, you can have conditional parts in your topologies and include other files. You can even compress most files and GROMACS will automatically pipe them through gzip upon reading.
• There is no scripting language - all programs use a simple interface with command line options for input and output files. You can always get help on the options by using the -h option, or use the extensive manuals provided free of charge in electronic or paper format. There is also an integrated graphical user interface available for all programs.
• As the simulation is proceeding, GROMACS will continuously tell you how far it has come, and what time and date it expects to be finished.
• Both run input files and trajectories are independent of hardware endian and can thus be read by any version GROMACS, even if it was compiled using a different floating-point precision. All files from GROMACS 2.0 can further be used in the new version 3!
• GROMACS can write coordinates using lossy compression, which provides a very compact way of storing trajectory data. The accuracy can be selected by the user.
• GROMACS comes with a large selection of flexible tools for trajectory analysis - you won't have to write any code to perform routine analyses. The output is further provided in the form of finished Xmgr/Grace graphs, with axis labels, legends, etc. already in place!
• A basic trajectory viewer that only requires standard X libraries is included, and several external visualization tools can read the GROMACS file formats.
• GROMACS can be run in parallel, using standard MPI communication.
• GROMACS contains several state-of-the-art algorithms that make it possible to extend the time steps is simulations significantly, and thereby further enhance performance without sacrificing accuracy or detail.
• The package includes a fully automated topology builder for proteins, even multimeric structures. Building blocks are available for the 20 standard aminoacid residues as well as some modified ones, the 4 nucleotide and 4 deoxinucleotide resides, several sugars and lipids, and some special groups like hemes and several small molecules.
• There is ongoing development to extend GROMACS with interfaces both to Quantum Chemistry and Bioinformatics/databases.
• GROMACS is Free Software, available under the GNU General Public License.

Lee2002hu

现在最常用的是MD+FEM方法。国内有本很好的书可以借鉴----

http://www.wl.cn/2355766
纳米力学与材料 - - 理论多尺度方法和应用(影印版)/微纳技术著作丛书合成与分析纳米物质特性的能力取得了革命性的进展，广泛应用于生物医学、机械、电子、精密材料以及军事工程等领域。纳米力学是研究和描述单个原子、系统和结构在各种载荷条件下响应的机械行为特性的学科，它的发展促进了该技术的进步。尤其是多尺度建模方法，它可以使此领域的工程师更好的理解微纳米材料。
本书由该领域内资深专家撰写，对纳米力学和材料的基本概念进行了介绍，侧重于多重尺度建模方法和技术的研究。本书内容包括：分子力学基础，微粒系统、晶格机械及现代多尺度建模理论等。
本书是一本相关领域电子工程师、材料科研工作者开发微纳米材料应用的全面的指南，同时也可作为微纳米力学和微纳米技术专业研究生的参考书。

huzhizi

这俺不太懂呀，路过

aresaran

原帖由 Lee2002hu 于 2008-5-12 18:34 发表
现在最常用的是MD+FEM方法。国内有本很好的书可以借鉴----

http://www.wl.cn/2355766
纳米力学与材料 - - 理论多尺度方法和应用(影印版)/微纳技术著作丛书合成与分析纳米物质特性的能力取得了革命性的进展，广泛 ...

不推荐将此方法太过宣传，因为DFM+MD+FEM的整体idea虽然合理，但是其中还有很多问题没有解决，

比如量级力学计算density function方法是完全没有参照的理论，是存在很多假设的理论，很多计算量子力学从根上就假设分子核不动，研究电子，所以可以说谁能算出来个结果，就可以声称自己是对的，完全没有比对，这个学科有很大的研究空间，可以发n多paper，但是绝对不能算作成熟；

Molecular dynamics是分子化学家提出来的基于牛顿运动力学的理论，不考虑电子，只研究分子核的运动问题，是强非线性问题；MD与FEM连接如何控制这些非线性响应，两种domain结合有很强的奇异性，等等，都是还没有成熟的问题，虽然MD＋FEM已经做了很久，可是本质问题应该是没有完全解决。

在这个行业，谁做的都有可能正确，也可能不正确，自己审视。

[ 本帖最后由 aresaran 于 2008-5-13 00:31 编辑 ]

seaseastar

不懂。不过普及一下科普。。。
谢谢大家。