复合材料弹性结构的高精度多尺度算法与数值模拟

１．引言 由于复合材料结构物理参数的非均匀各向异性以及细部几何构形的复杂性,在计算它的位移场、应力、应变场时,传统的有限元法因网格生成困难和计算规模太大而难以实现．７０年代初,Ｉ．Ｂａｂｕｓｋａ,Ｊ．Ｌ．Ｌｉｏｎｓ等人针对复合材料弹性结构提出了均匀化方法,见文［１］,数值实验表明,均匀化方法对描述复合材料弹性结构的有效材料常数及刚度性质是有效的,但它不能刻画应力和应变场的局部变化,因而作为复合材料强度理论的判断依据,均匀化理论显然是不够的．为此,Ｊ．Ｌ．Ｌｉｏｎｓ,Ｏ．Ａ．Ｏｌｅｉｎｉｋ等分别就…     

Large eddy simulation of turbulent incompressible flows by a three‐level finite element method

The variational multiscale method provides a methodical framework for large eddy simulation of turbulent flows. In this work, a particular implementation in the form of a three‐level finite element method separating large resolved, small resolved, and unresolved scales is proposed. Residual‐free bubbles are used for the numerical approximation of the small‐scale momentum equation. A stabilizing term is added, in order to take into account the effect of the small‐scale continuity equation. This implementation guarantees the stability of the method without further provisions and offers substantial computational savings on the small‐scale level. Furthermore, it is accounted for the unresolved scales by a specific dynamic modelling procedure. The method is tested for two different turbulent flow situations. Copyright © 2005 John Wiley & Sons, Ltd.     

Experimental determination of cohesive failure properties of a photodegradable copolymer

In this paper we present a methodology to measure the material traction-separation relation for a poly(ethylene carbon monoxide) copolymer, ECO. This material exhibits a ductile-to-brittle transition when subjected to ultraviolet irradiation and undergoes a change of failure mechanism, from shear yielding to crazing, in the process. Single edge notched tension specimens of ECO irradiated for 50 h were used to generate slow-speed stable crack growth, predominantly from material crazing. Full-field measurements of the in-plane deformation around the growing crack tip were performed using the optical technique of digital image correlation. A multicamera setup was used in which simultaneous measurement of both the far-field displacement and that directly surrounding the craze was performed. The far-field results were used to obtain a value of the energy release rate supplied to the crack tip region. The near-tip results were used to extract a material traction-separation law in the regime of steady-state crack growth. A softening traction-separation relation was measured. The area under the traction-separation curve was then compared with the simultaneous far-field measurements and the agreement was very good (within 6.5%) validating the experimental methodology used.     

Analysis of the heterogeneous multiscale method for parabolic homogenization problems

The heterogeneous multiscale method (HMM) is applied to various parabolic problems with multiscale coefficients. These problems can be either linear or nonlinear. Optimal estimates are proved for the error between the HMM solution and the homogenized solution.

     

Bridging the scales in nano engineering and science

This review article describes various multiscale approaches, development of which was spurred by the emergence of nanotechnology. The multiscale approaches are grouped into two main categories: information-passing and concurrent. In the concurrent multiscale methods both, the discrete and continuum scales are simultaneously resolved, whereas in the information-passing schemes, the discrete scale is modelled and its gross response is infused into the continuum scale. Most of the information-passing approaches provide sublinear computational complexity, (i.e., scales sublinearly with the cost of solving a fine scale problem), but the quantities of interest are limited to or defined only on the coarse scale. The issues of appropriate scale selection and uncertainty quantification are also reviewed.     

Multiscale methods for elliptic homogenization problems

In this article we study two families of multiscale methods for numerically solving elliptic homogenization problems. The recently developed multiscale finite element method [Hou and Wu, J Comp Phys 134 (1997), 169–189] captures the effect of microscales on macroscales through modification of finite element basis functions. Here we reformulate this method that captures the same effect through modification of bilinear forms in the finite element formulation. This new formulation is a general approach that can handle a large variety of differential problems and numerical methods. It can be easily extended to nonlinear problems and mixed finite element methods, for example. The latter extension is carried out in this article. The recently introduced heterogeneous multiscale method [Engquist and Engquist, Comm Math Sci 1 (2003), 87–132] is designed for efficient numerical solution of problems with multiscales and multiphysics. In the second part of this article, we study this method in mixed form (we call it the mixed heterogeneous multiscale method). We present a detailed analysis for stability and convergence of this new method. Estimates are obtained for the error between the homogenized and numerical multiscale solutions. Strategies for retrieving the microstructural information from the numerical solution are provided and analyzed. Relationship between the multiscale finite element and heterogeneous multiscale methods is discussed. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Simulation of plastic deformation in glassy polymers: Atomistic and mesoscale approaches

The mechanism of deformation in glasses is very different from that of crystals, even though their general behavior is very similar. In this study, we investigated the deformation of polycarbonate on the atomistic scale with molecular dynamics and on the continuum scale with a new simulation approach. The results indicated that high atomic/segmental mobility and low local density enabled the formation (nucleation) of highly deformed regions that grew to form plastic defects called plastic shear transformations. A continuum-scale simulation was performed with the concept of plastic shear transformations as the basic region of deformation. The continuum simulations were able to predict the primary and secondary creep behavior. The slope of the secondary creep depended on the interactions between the plastic shear transformations. © 2005 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 43: 994-1004, 2005     

粘弹性介质中波场重建方法

传统的波动方程波场重建基于完全弹性介质，不能获得满意的地震资料分辨率，本基于能描述大地吸收弹性介质中地震波传播的斯托克斯波动方程．提出了一种新的多尺度粘弹性波动方程波场重建的方法，根据地震波传播核函数的物理特性，研究了一种新的物理小波，提出了小波多尺度波场重建方法。达到对吸收信息的补偿，提高地震资料的分辨率。     

应用小波边缘检测技术分析荧光寿命显微像

研究了小波图象边缘检测技术在荧光寿命显微像分析中的应用，利用小波分析的多尺度特性进行了荧光寿命显微像的边缘提取，结果表明，这种方法非常有效。