周期荷载下半透水边界饱和地基的固结

在Terzaghi一维固结假设的基础上，分析了半透水边界饱和土层在周期荷载作用下有效应力比的变化规律。通过基本方程，结合初边值条件，利用Laplace变换，得到了随时间任意变化荷载作用下的有效应力解。通过典型算例，结合编制的数值Laplace逆变换程序，对几种常见的荷载作用下的情况进行了讨论，研究了参数的影响。通常所研究的透水或不透水边界条件只是本文的特例。本文的结论可以用于指导工程实践。     

基于周期性边界条件的炭黑填充橡胶复合材料力学行为的预测

在分析炭黑填充橡胶复合材料的宏观与细观特征之间联系的基础上,提出了具有随机分布形态的代表性体积单元,推导并应用了周期性细观结构的边界约束条件,建立了三维多颗粒夹杂代表性体积单元的数值模型,对炭黑填充橡胶复合材料的宏观力学行为进行了模拟仿真。研究表明,该模型通过周期性边界条件的约束保证了宏观结构变形场和应力场的协调性;计算得到的炭黑填充橡胶复合材料的弹性模量明显高于未填充橡胶材料,并随着炭黑颗粒所占体积分数的增加而增大;该模型对复合材料有效弹性模量的预测结果与实验结果吻合较好,而且比Bergstrom三维模型的预测结果更好,证实了该模型能够用于炭黑颗粒增强橡胶基复合材料有效性能的模拟分析。     

Finite Element Multi-mode Approach to Thermal Postbuckling of Functionally Graded Plates

Postbuckling analysis of functionally graded ceramic-metal plates under temperature field is presented using finite element multi-mode method. The three-node triangular element based on the Mindlin plate theory is employed to account for the transverse shear strains, and the von-Karman nonlinear strain-displacement relation is utilized considering the geometric nonlinearity. The effective material properties are assumed to vary through the thickness direction according to the power law distribution of the volume fraction of constituents. The temperature distribution along the thickness is determined by one dimensional Fourier equations of heat conduction. The buckling mode shape solved from eigen-buckling analysis is adopted as the assumed mode function to reduce the degrees of freedom of nonlinear postbuckling equilibrium equations. The postbuckling response is obtained by solving the nonlinear equilibrium equations, and compared with the Newton- Raphson numerical results. The effects of boundary conditions, material gradient index and temperature distribution on postbuckling behavior are examined.     

On the Somigliana stress identity in elasticity

The paper presents and analytical investigation of the properties of the integral terms in the Somigliana stress identity which are associated with a perceived hypersingular nature of the integral equations in three dimensional elasticity. The nature of the integral equations is, in fact, found to be non-hypersingular, thereby permitting direct evaluation of the jump terms in the stress identity for a solution point taken, in the limit, to the surface of the body. The continuity requirements on the boundary conditions are found to be more liberal than previously reported. A weakly-singular form of the Somigliana identity is found that is easily used for BEM implementations that use Gaussian integrations. Demonstration of the boundary form of the Somigliana stress identity is given for a three dimensional elasticity problem.     

Matching boundary conditions for diatomic chains

We propose a class of efficient matching boundary conditions to suppress spurious reflection for multiscale computations of one dimensional diatomic chains. This provides the first local effective treatment of both acoustic and optical phonons. Adopting the extended zone scheme of the dispersion relation, we design a class of force boundary conditions by enforcing perfect absorption at certain selected wave numbers. Reflection suppression is improved by involving more neighboring atoms in the condition. The effectiveness of the proposed matching boundary conditions is demonstrated by reflection coefficient analysis, numerical tests, and comparisons with the time history treatment.     

Experimental study of dynamic buckling of plates under fluid–solid slamming

This paper describes the experimental investigation of the elastic–plastic dynamic buckling properties of rectangular plates under in-plane fluid–solid slamming. Based on the observation of elastic–plastic dynamic response characteristics of the plates, a dynamic buckling criterion and a dynamic yielding criterion are defined. The corresponding critical impulse are determined from the experimental results for each tested plate. The effect of different boundary conditions on the elastic–plastic dynamic buckling properties of plates is also examined. The results indicate that dynamic buckling always takes place elastically for the types of rectangular plates tested under fluid–solid slamming. The dynamic buckling modes of the plates are governed by the plate fundamental transverse free vibration mode. It is also found that boundary conditions strongly affect the dynamic buckling properties of plates subjected to fluid–solid slamming loads. Strengthening plate boundary constraint is a very effective way to enhance the plates’ ability to resist dynamic buckling.     

Application of essential boundary conditions in mesh‐free methods: a corrected collocation method

Collocation methods for applying essential boundary conditions are defined as those methods in which conditions are enforced exactly at a discrete set of boundary nodes. In mesh‐free methods, this is usually accomplished by replacing rows of the matrix equations which result from discretization of the weak form with equations which ensure the enforcement of boundary conditions. In this paper, an inconsistency in this method is pointed out, and a correction is derived. Numerical test are done on one‐ and two‐dimensional equations; it is shown that convergence rates decrease with the use of the invalid traditional collocation and are restored with the corrected method. Copyright © 2000 John Wiley & Sons, Ltd.     

Microstructure modelling and prediction of effective elastic properties of 3D multiphase and multilayer braided composite

Abstract

This paper is focused on the microstructure modelling and evaluation of effective elastic properties of three-dimensional multiphase and multilayer braided composite. Regarding the multiscale characteristics of the composite, the microstructure modelling is carried out sequentially from fibre to tow scale. The geometrical configuration of the microstructure is first analysed, and mathematical relations among different geometrical parameters are derived on each scale. Second, effective elastic properties are obtained based on the sequential homogenisation from fibre to tow scale. A strain energy based method is proposed to evaluate effective elastic properties with specific boundary conditions imposed on the microstructure. Numerical results obtained by the proposed method and the microstructure model show a good agreement with the results measured experimentally.     

Trefftz boundary elements—multi‐region formulations

This paper is concerned with an effective numerical implementation of the Trefftz boundary element method, for the analysis of two‐dimensional potential problems, defined in arbitrarily shaped domains. The domain is first discretized into multiple subdomains or regions. Each region is treated as a single domain, either finite or infinite, for which a complete set of solutions of the problem is known in the form of an expansion with unknown coefficients. Through the use of weighted residuals, this solution expansion is then forced to satisfy the boundary conditions of the actual domain of the problem, leading thus to a system of equations, from which the unknowns can be readily determined. When this basic procedure is adopted, in the analysis of multiple‐region problems, proper boundary integral equations must be used, along common region interfaces, in order to couple to each other the unknowns of the solution expansions relative to the neighbouring regions. These boundary integrals are obtained from weighted residuals of the coupling conditions which allow the implementation of any order of continuity of the potential field, across the interface boundary, between neighbouring regions. The technique used in the formulation of the region‐coupling conditions drives the performance of the Trefftz boundary element method. While both of the collocation and Galerkin techniques do not generate new unknowns in the problem, the technique of Galerkin presents an additional and unique feature: the size of the matrix of the final algebraic system of equations which is always square and symmetric, does not depend on the number of boundary elements used in the discretization of both the actual and region‐interface boundaries. This feature which is not shared by other numerical methods, allows the Galerkin technique of the Trefftz boundary element method to be effectively applied to problems with multiple regions, as a simple, economic and accurate solution technique. A very difficult example is analysed with this procedure. The accuracy and efficiency of the implementations described herein make the Trefftz boundary element method ideal for the study of potential problems in general arbitrarily‐shaped two‐dimensional domains. Copyright © 1999 John Wiley & Sons, Ltd.     

Three dimensional boundary formulations for nonlinear thermal shape sensitivities

Implicit differentiation of the discretized boundary integral equations governing the conduction of heat in three dimensional (3D) solid objects, subjected to nonlinear boundary conditions, and with temperature dependent material properties, is shown to generate an accurate and economical approach for the computation of shape sensitivities. The theoretical formulation for primary response (surface temperature and normal heat flux) sensitivities and secondary response (surface tangential heat flux components and internal temperature and heat flux components) sensitivities is given. Iterative strategies are described for the solution of the resulting sets of nonlinear equations and computational performances examined. Multi-zone analysis and zone condensation strategies are demonstrated to provide substantial computational economies in this process for models with either localized nonlinear boundary conditions or regions of geometric insensitivity to design variables. A series of nonlinear sensitivity example problems are presented that have closed form solutions. Sensitivities computed using the boundary formulation are shown to be in excellent agreement with these exact expressions.     

微观结构对复合材料弹性有效性能的影响

在渐近均匀化方法的基础上,用ANSYS参数设计语言建立了周期性边界条件,用ANSYS有限元程序对单胞进行求解,得到了复合材料的有效性能。分析了不同微观结构对材料有效性能的影响,并与实验和其它理论结果进行比较。得到了不同方向的方形纤维对于材料的有效模量和有效泊松比的影响。     

On stability and reflection‐transmission analysis of the bipenalty method in contact‐impact problems: A one‐dimensional,homogeneous case study

The stability and reflection‐transmission properties of the bipenalty method are studied in application to explicit finite element analysis of one‐dimensional contact‐impact problems. It is known that the standard penalty method, where an additional stiffness term corresponding to contact boundary conditions is applied, attacks the stability limit of finite element model. Generally, the critical time step size rapidly decreases with increasing penalty stiffness. Recent comprehensive studies have shown that the so‐called bipenalty technique, using mass penalty together with standard stiffness penalty, preserves the critical time step size associated to contact‐free bodies. In this paper, the influence of the penalty ratio (ratio of stiffness and mass penalty parameters) on stability and reflection‐transmission properties in one‐dimensional contact‐impact problems using the same material and mesh size for both domains is studied. The paper closes with numerical examples, which demonstrate the stability and reflection‐transmission behavior of the bipenalty method in one‐dimensional contact‐impact and wave propagation problems of homogeneous materials.     

Reflection of plane waves at the initially stressed surface of a fiber-reinforced thermoelastic half space with temperature dependent properties

In this paper, a model of two dimensional problem of generalized thermoelasticity for a fiber-reinforced anisotropic elastic medium under the effect of temperature dependent properties is established. Reflection phenomena of plane waves in an initially stressed thermoelastic medium is studied in the context of two theories proposed by Lord–Shulman and Green–Lindsay. Using proper boundary conditions, the amplitude ratios and energy ratios for various reflected waves are presented. The phase speeds, reflection coefficients and energy ratios are computed numerically with the help of MATLAB programming and are depicted graphically to show the effect of initial stress and temperature dependent properties. It is found that there is no dissipation of energy at the boundary surface during reflection. A comparison between the two theories is also depicted in the present investigation.

     

Skeletal reduction of boundary value problems

Boundary value problems posed over thin solids are amenable to a dimensional reduction in that one or more spatial variables may be eliminated from the governing equation, resulting in significant computational gains with minimal loss in accuracy. Extant dimensional reduction techniques rely on representing the solid as a hypothetical mid‐surface plus a possibly varying thickness. Such techniques are however hard to automate since the mid‐surface is often ill‐defined and ambiguous. We propose here a skeletal representation based dimensional reduction of boundary value problems. The proposed technique has the computational advantages of mid‐surface reduction, but can be easily automated. A systematic methodology is presented for reducing boundary value problems to lower‐dimensional problems over the skeleton of a solid. The theoretical properties of the proposed method are derived, and supported by representative numerical experiments. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical studies on the effective shear modulus of particle reinforced composites with an inhomogeneous inter-phase

Numerical studies on the effective shear modulus of particle reinforced composites with an inhomogeneous inter-phase are performed. The influences of many parameters to the equivalent shear modulus of composites are carefully analyzed, including the inter-phase thickness, variation of interfacial properties, boundary conditions and volume fraction of particles, etc. Numerical results show that the Poisson ratio can be assumed as a constant across the whole inter-phase zone in the computation. The form of property variation across the inter-phase also greatly affects the effective shear modulus of composite. Numerical results predicted by the rigid boundary conditions are remarkably higher than those by the free boundary conditions and the exact solutions. The reasonability and exactness of the available models for predicting the effective shear modulus of composites are accessed by the numerical results in the present work.     

Lubrication theory for micropolar fluids and its application to a journal bearing

The problem considered is that of a steady laminar flow of an incompressible micropolar fluid in a journal bearing. A two dimensional flow field is considered and order of magnitude arguments are made, which reduce the governing balance equations to a system of coupled, ordinary differential equations. The equations are solved subject to appropriate boundary conditions and the bearing characteristics obtained. The effect of microstructure is elaborated through various graphs. The prominent feature of a micropolar fluid is an increased effective viscosity especially in thin films.     

Coupled heat conduction and thermal stress analyses in particulate composites

This study introduces two micromechanical modeling approaches to analyze spatial variations of temperatures, stresses and displacements in particulate composites during transient heat conduction. In the first approach, a simple micromechanical model based on a first order homogenization scheme is adopted to obtain effective mechanical and thermal properties, i.e., coefficient of linear thermal expansion, thermal conductivity, and elastic constants, of a particulate composite. These effective properties are evaluated at each material (integration) point in three dimensional (3D) finite element (FE) models that represent homogenized composite media. The second approach treats a heterogeneous composite explicitly. Heterogeneous composites that consist of solid spherical particles randomly distributed in homogeneous matrix are generated using 3D continuum elements in an FE framework. For each volume fraction (VF) of particles, the FE models of heterogeneous composites with different particle sizes and arrangements are generated such that these models represent realistic volume elements “cut out” from a particulate composite. An extended definition of a RVE for heterogeneous composite is introduced, i.e., the number of heterogeneities in a fixed volume that yield the same expected effective response for the quantity of interest when subjected to similar loading and boundary conditions. Thermal and mechanical properties of both particle and matrix constituents are temperature dependent. The effects of particle distributions and sizes on the variations of temperature, stress and displacement fields are examined. The predictions of field variables from the homogenized micromechanical model are compared with those of the heterogeneous composites. Both displacement and temperature fields are found to be in good agreement. The micromechanical model that provides homogenized responses gives average values of the field variables. Thus, it cannot capture the discontinuities of the thermal stresses at the particle-matrix interface regions and local variations of the field variables within particle and matrix regions.     

Essential boundary condition enforcement in meshless methods: boundary flux collocation method

Element‐free Galerkin (EFG) methods are based on a moving least‐squares (MLS) approximation, which has the property that shape functions do not satisfy the Kronecker delta function at nodal locations, and for this reason imposition of essential boundary conditions is difficult. In this paper, the relationship between corrected collocation and Lagrange multiplier method is revealed, and a new strategy that is accurate and very simple for enforcement of essential boundary conditions is presented. The accuracy and implementation of this new technique is illustrated for one‐dimensional elasticity and two‐dimensional potential field problems. Copyright © 2001 John Wiley & Sons, Ltd.     

Three-dimensional free vibration analysis of functionally graded annular sector plates with general boundary conditions

The main purpose of this paper is to investigate free vibration behaviors of functionally graded sector plates with general boundary conditions in the context of three-dimensional theory of elasticity. Generally, the material properties of functionally graded sector plates are assumed to vary continuously and smoothly in thickness direction. However, the changes in the material properties may occur in the other directions, such as radial direction. Therefore, two types of functionally graded annular sector plates are considered in the paper. In this work, both the Voigt model and Mori-Tanaka scheme are adopted to evaluate the effective material properties. Each of displacements of annular sector plate, regardless of boundary conditions, is expressed as modified Fourier series which consists of three-dimensional Fourier cosine series plus several auxiliary functions introduced to overcome the discontinuity problems of the displacement and its derivatives at edges. To ensure the validity and accuracy of the method, numerous examples for isotropic and functionally graded sector plates with various boundary conditions are presented. Furthermore, new results for functionally graded sector plates with elastic restraints are given. The effects of the material profiles and boundary conditions on the free vibration of the functionally sector plates are also studied.     

冷冻外科中一维相变问题的数值解

用变时间步长法求解了冷冻外科手术时模拟生物组织中的一维相变问题。计算方法经与文献比较,简便可行。运用本文的计算方法计算并讨论了不同边界条件下,相变区的温度场分布和无量纲相界面降温速率随Ｓｔｅ数的变化,计算分析表明,相变区的温度基本呈线性分布,无量纲相界面降温速率随Ｓｔｅ数增大而降低,所得结论可对冷刀设计和冷冻外科手术提供帮助。