随机桁架结构的非平稳随机动力响应分析

本文研究了随机桁架结构在非平稳随机激励下的动力响应问题。在利用随机因子法分析随机结构动力特性的基础上，给出了一种分析随机结构非平稳随机响应的新方法。从结构非平稳随机响应的表达式出发，同时考虑桁架结构的物理参数、几何尺寸的随机性，利用求解随机变量函数矩的方法和求解随机变量数字特征的代数综合法，导出了随机桁架结构在非平稳随机激励下位移响应均方值和应力响应均方值的均值、方差和变异系数的计算表达式。通过算例，分析了结构物理参数和几何尺寸的随机性对结构位移响应均方值和应力响应均方值随机变量随机性的影响。本文方法具有对随机结构进行一次动力分析便可求得动力响应的数字特征，且可以考察结构任一参数的随机性对结构非平稳随机响应分析结果的影响之优点。     

随机结构有限元分析的递推求解方法

提出了一种新的谱随机有限元分析方法——递推求解方法。该方法将随机结构的随机响应表示成非正交多项式展式,建立了和摄动法类似的一系列确定的递推方程,并通过确定性有限元方法对这些递推方程进行静力问题求解。算例表明,当随机量出现较大涨落时,计算结果相对于传统摄动法有不小的改进。     

Sensitivity and randomness in homogenization of periodic fiber-reinforced composites via the response function method

The main issue this paper addresses is the derivation and implementation of a general homogenization method, including the simultaneous determination of sensitivity gradients and probabilistic moments of the effective elasticity tensor. This is possible with an application of the perturbation method based on Taylor expansion and with the effective modules method. The computational procedure is implemented using plane strain analysis carried out with the finite element method (program MCCEFF) and the symbolic computations system MAPLE. The sensitivity gradients and probabilistic moments are commonly determined on the basis of partial derivatives for the homogenized elasticity tensor, calculated using the response function method with respect to some composite parameters. They are subjected separately to a normalization procedure (in deterministic analysis) and the relevant algebraic combinations (for the stochastic case). This enriched homogenization procedure is tested on a periodic fiber-reinforced two component composite, where the material parameters are taken as design variables and then, the input random quantities. The results of computational analysis are compared against the results of the central finite difference approach in the case of sensitivity gradients determination as well as the direct Monte-Carlo simulation approach. This numerical methodology may be further applied not only in the context of the homogenization method, but also to extend various discrete computational techniques, such as Boundary/Finite element and finite difference together with various meshless methods.     

Multiscale homogenization of n-component composites with semi-elliptical random interface defects

Effective elastic characteristics of periodic multicomponent composite materials with random interface defects are studied in the paper. The defects are assumed to be semi-elliptical and lying with major semi axes along the interfaces, where minor and major semi-axes as well as the defects number are given as input random variables. The homogenization approach has a multiscale character—some algebraic approximation is used first to calculate effective elastic parameters of the interphase including all defects located at the same interface. Equations for interphase random elastic parameters are obtained using MAPLE symbolic mathematics in conjunction with probabilistic generalized perturbation method. A different homogenization method is applied at the micro scale, where the cell problem is solved numerically using the Finite Element Method (FEM) program. Since the composites considered exhibit random variations of both elastic properties and the interface defects, the overall homogenized characteristics must be obtained as random quantities, which is realized on the micro scale by the Monte-Carlo simulation. The proposed interface defects model obeys the porosity effects resulting from the nature of some matrices in engineering composites as well as the interface cracks appearing as a result of composites ageing during static or fatigue fracture.     

Analysis of uncertainty effects due to microstructural disorder in cellular or porous materials

Aim of the present study is an analysis of the effect of microstructural uncertainties on the scatter in the macroscopic material properties of highly porous materials consisting of metallic or other constituents. For the numerical analysis of the uncertainty effects, a probabilistic homogenization scheme is proposed. In contrast to direct Monte-Carlo approaches, the thermomechanical response of a limited number of pre-selected cases throughout the range of possible microstructures is analyzed. Their effective properties are determined by means of an energy based homogenization procedure. In a stochastic evaluation, the results of the individual computations are weighted with the probability of the occurrence of the underlying microstructures. As a result, the probability distributions for the effective properties are obtained. The basic uncertain microstructural properties considered in the investigation are the microstructural geometry and orientation, the local relative density and the local pore size distribution. In an application to an experimental data base from other sources, the approach proves to be accurate and numerically efficient compared to direct Monte-Carlo approaches. Parameter studies reveal that uncertainties in the local relative density are the most important factor leading to scatter in the macroscopic material properties of cellular materials.     

Stochastic thermal stresses in an FGM annular disc of variable thickness with spatially random heat transfer coefficients

The second-order statistics (i.e. mean and standard deviation) of the temperature and thermal stresses are evaluated in an axisymmetrically heated functionally graded annular disc of variable thickness with spatially random heat transfer coefficients (HTCs) on the major surfaces of the disc. This annular disc is assumed to have arbitrary variations in the HTCs and material composition along the radial direction only. The randomness in the HTCs is considered to be a random field. The stochastic temperature field is analysed by considering the annular disc to be a multilayered one with stepwise thickness variation, where each layer is assumed to have constant deterministic material properties and random HTCs. In order to evaluate the statistics, the Monte Carlo simulation method is applied to analytical solutions for the deterministic temperature and thermal stresses. The analytical solution for the thermal stresses is obtained through the use of a piecewise power function approximation for Young’s modulus. Numerical results demonstrate the effects of the magnitude of the HTC means, volume fraction distributions of the constitutive materials and thickness variation on the statistics of the temperature and thermal stresses.     

Micromechanical modeling of linear viscoelastic behavior of heterogeneous materials

Study of effective behavior of heterogeneous materials, starting from the properties of the microstructure, represents a critical step in the design and modeling of new materials. Within this framework, the aim of this work is to introduce a general internal variables approach for scale transition problem in linear viscoelastic case. A new integral formulation is established, based on the complete taking into account of field equations and differential constitutive laws of the heterogeneous problem, in which the effects of elasticity and viscosity interact in a representative volume element. Thanks to Green’s techniques applied to space convolution’s term, a new concentration relation is obtained. The step of homogenization is then carried out according to the self-consistent approximation. The results of the present model are illustrated and compared with those provided by Hashin’s and Rougier’s ones, considered as references, and by internal variables models such as those of Weng and translated fields.     

埋地管道随机振动的摄动分析

埋地管道的材料特性和沿线的土壤性质存在差异，对这种差异性采用随机参数描述有一定的合理性，因此在管道的随机振动分析中考虑结构参数的随机性是必要的，对于管道的抗震设计具有现实意义。对于在空间相关的地震地面随机激励下的埋地管道，将结构参数看作是随机变量，采用摄动分析法，推导了随机响应的相关函数和功率谱密度函数的解析表达式，大大方便了工程应用。针对某输油管道，给出了计算结果。     

Dielectric elastomer composites: A general closed-form solution in the small-deformation limit

A solution for the overall electromechanical response of two-phase dielectric elastomer composites with (random or periodic) particulate microstructures is derived in the classical limit of small deformations and moderate electric fields. In this limit, the overall electromechanical response is characterized by three effective tensors: a fourth-order tensor describing the elasticity of the material, a second-order tensor describing its permittivity, and a fourth-order tensor describing its electrostrictive response. Closed-form formulas are derived for these effective tensors directly in terms of the corresponding tensors describing the electromechanical response of the underlying matrix and the particles, and the one- and two-point correlation functions describing the microstructure. This is accomplished by specializing a new iterative homogenization theory in finite electroelastostatics (Lopez-Pamies, 2014) to the case of elastic dielectrics with even coupling between the mechanical and electric fields and, subsequently, carrying out the pertinent asymptotic analysis.Additionally, with the aim of gaining physical insight into the proposed solution and shedding light on recently reported experiments, specific results are examined and compared with an available analytical solution and with new full-field simulations for the special case of dielectric elastomers filled with isotropic distributions of spherical particles with various elastic dielectric properties, including stiff high-permittivity particles, liquid-like high-permittivity particles, and vacuous pores.     

Field fluctuations in multicomponent mixtures

Fluctuations of electrostatic and elastostatic fields in a random phase mixture may be characterized by mean values and square means of the fields in each component. Exact relations between the square means and the analytical properties of the effective moduli are established for isotropic mixtures. Moreover, a modified effective medium procedure for calculating the field fluctuations in mixtures with aggregate topology is proposed. Explicit results are given for mixtures of isotropic components and spherical grain shapes. Particularly strong fluctuations occur in strongly heterogeneous media near the percolation threshold.     

Numerical analysis of effective elastic properties of geomaterials containing voids using 3D random fields and finite elements

The purpose of the study is to investigate the influence of porosity and void size on effective elastic geotechnical engineering properties with a 3D model of random fields and finite element. The random field theory is used to generate models of geomaterials containing spatially random voids with controlled porosity and void size. A “tied freedom” analysis is developed to evaluate the effective Young’s modulus and Poisson’s ratio in an ideal block material of finite elements. To deliver a mean and standard deviation of the elastic parameters, this approach uses Monte-Carlo simulations and finite elements, where each simulation leads to an effective value of the property under investigation. The results are extended to investigate an influence of representative volume element (RVE). A comparison of the effective elastic stiffness of 2D and 3D models is also discussed.     

Periodic Homogenization and Material Symmetry in Linear Elasticity

Here homogenization theory is used to establish a connection between the symmetries of a periodic elastic structure associated with the microscopic properties of an elastic material and the material symmetries of the effective, macroscopic elasticity tensor. Previous results of this type exist but here more general symmetries on the microscale are considered. Using an explicit example, we show that it is possible for a material to be fully anisotropic on the microscale and yet the symmetry group on the macroscale can contain elements other than plus or minus the identity. Another example demonstrates that not all material symmetries of the macroscopic elastic tensor are generated by symmetries of the periodic elastic structure.     

Homogenization of Polycrystalline¶Aggregates

The effective elastic properties of a polycrystalline material depend on the single crystal elastic constants of the crystallites comprising the polycrystal and on the manner in which the crystallites are arranged. In this paper we apply the techniques of homogenization to put the problem of determining effective elastic constants in a precise mathematical framework that permits us to derive an expression for the effective elasticity tensor. We also study how the homogenized elasticity tensor changes as the probability characterizing the ensemble changes. Under the assumption that the field of orientations of the crystallographic axes of the crystallites is an independent random field, we show that our theory is compatible with the formulation used in texture analysis. In particular, we are able to prove that the physical assumption made by [10] in his study of weakly-textured polycrystals holds true. In addition, we establish some elementary bounds on the material constants that characterize the effective elasticity tensor of weakly-textured orthorhombic aggregates of cubic crystallites. Accepted: (June 15, 1999)     

Stochastic Simulations in Dynamic Soil–Structure Interaction

This paper deals with several numerical techniques that account for random excitations and random material parameters occurring in earthquake engineering. It focuses in sequence on the influence on the structural response of the variability of the incident field using filtering theory, on the soil variability by coupling Stochastic Finite Elements, integral operators and Monte-Carlo simulation, and finally on the influence on the site response of a random building distribution using periodic simulations. In each case, practical examples are given to help illustrate the numerical techniques and to underline the importance of randomness in earthquake engineering problems.     

Stochastic modeling of multi-filament yarns. I. Random properties within the cross-section and size effect

The paper presents a thorough study of sources of randomness/disorder in multi-filament yarns that are relevant for their performance in textile reinforced concrete. For the analysis of bundles with infinite number of fibers we use a continuous analytical model with refined kinematic hypothesis. In parallel, we apply a more general discrete numerical bundle model covering the cases of finite number of fibers and spatial variations of material properties. The considered distributions of material properties include variations of parameters (local stiffness and strength) both from filament to filament and over the length of each filament. In the present paper, we study the influence of variations in filament properties across the bundle. In the companion paper (Vořechovský, M., Chudoba, R., 2006. Stochastic modeling of multi-filament yarns. II. Random properties over the length and size effect. International Journal of Solids and Structures), the variations of material properties over the length are considered. Both papers provide the basis for correct interpretation of the data obtained from the tensile test on multi-filament yarn with varied specimen length. As the actual goal of the present work is to derive yarn characteristics relevant for crack bridges in cementitious composites we systematically assess the significance of the studied sources of disorder in the context of extremely short bundle lengths.     

Consolidation Statistics Investigation via Thin Layer Method Analysis

In this paper, the Thin Layer Method (TLM) is adapted for solving one-dimensional primary consolidation problems. It is also combined with a stochastic formulation integrating Monte Carlo simulations to investigate primary consolidation of a random heterogeneous soil profile. This latter is modeled as a set of superposed layers extending horizontally to infinity, and having random properties. Spatial variability of soil properties is considered in the vertical direction only. Soil properties of interest are elastic modulus and soil permeability, modeled herein as spatially random fields. Lognormal distribution is chosen because it is suitable for strictly non-negative random variables, and enables analyzing the large variability of such properties. The statistics regarding final settlement and its corresponding time are investigated by performing a parametric study, which integrates the influence of variation coefficient of both elastic modulus and soil permeability, and vertical correlation length. Obtained results indicate that heterogeneity significantly influences primary consolidation of the soil profile, generating a quite different way of soil grain rearrangement and water pressure dissipation in comparison to the homogeneous case, and causing a delay in the consolidation process.     

含随机参数的非线性黏弹性有限元计算

进行了粗粒土与结构接触面单调和循环加载试验，基于宏细观测量结果, 扩展了 损伤概念以 描述该类接触面在受载过程中的物态演化, 及由于物态演化导致的力学特性从初始状态到最终 稳定状态的连续变化过程. 揭示了接触面损伤的细观物理基础主要是接触面内土的颗粒破碎 和剪切压密这两种物态演化；指出接触面的剪胀体应变可以划分为可逆性和不可逆性剪胀体 应变两部分，其中不可逆性剪胀体应变可作为接触面损伤发展的宏观量度，因此其归一化 形式可作为一种损伤因子的定义；提出了建立粗粒土与结构接触面一种损伤本构关系的基本思路.     

平面弹性力学充要的随机边界元

将平面弹性力学确定性的充分必要的边界积分方程推广到含材料常数随机的不确定问题中去，给出了位移的均值以及偏差的充分必要的边界积分方程。数值计算结果表明，和确定性的积分方程一样，习用的随机边界积分方程在退化尺度附近，无论是均值还是偏差都存在巨大的误差，而充要的随机边界积分方程则始终保持良好的精度     

Virtual improvement of ice cream properties by computational homogenization of microstructures

Optimal shape design of microstructured materials has recently attracted a great deal of attention in materials science. The shape and the topology of the microstructure have a significant impact on the macroscopic properties. This paper presents different computational models of random microstructures, to virtually improve the physical properties of ice cream. Several sensory properties of this heterogeneous material issued from food industry are directly controlled by the elastic and thermal conducting ones. The material effective elastic and thermal conducting properties are obtained through direct large scale numerical simulations. The different formulations address the problem of finding the shape of the representative microstructural element for random heterogeneous media that increase the elastic moduli and thermal conductivity compared to existing products. The computational models are established using finite element method and images of virtual microstructures. In this paper we propose a new model of microstructures. This model is constructed with hexagonal prismatic rods and plates with volume fractions around 0.7 for the hard phase represented by hexagons of ice. A comparison between three two-phase elastic heterogeneous microstructures models is drawn. This illustrates the concept of design of microstructures using computational homogenization tools.