Methods of analysing ropes. The extension–torsion method

Two new approaches are used for calculating the stress–strain state of a rope and its stiffnesses. The first approach relies on the theory of fibrous composites and Saint-Venant's solution for a cylinder with helical anisotropy. The second approach is based on the solution by the finite element method of the three-dimensional problem of elasticity theory for a solid inhomogeneous cylinder formed by a finite number of elastic fibres arranged in helical lines and connected by a weak filler (in the sense that its Young's modulus is several orders of magnitude less than the Young's modulus of the fibre). The behaviour of the stiffness when the modulus of elasticity of the filler tends to zero is analysed, and the results of the limiting transition are discussed. The numerical results obtained are compared with calculations by other well-known applied theories.     

Size-dependent non-linear dynamics of curvilinear flexible beams in a temperature field

A mathematical model of the loss of dynamic stability of curvilinear size-dependent MEMS and NEMS elements embedded in a temperature field and subjected to large deflections was derived and studied. The fundamental governing dynamical equations of MEMS/NEMS members were yielded by Hamilton's principle. The investigations were based on combining the modified couple stress theory, the first-order approximation kinematic (Euler–Bernoulli) model, the von Kármán geometric non-linearity, and the Duhamel–Neumann law regarding the temperature input (the beam material is elastic, isotropic and there are no constraints imposed on the temperature distribution). The temperature field was defined by solving a heat transfer equation. The computational algorithm was based on the finite difference method and the Runge–Kutta method. The numerical methods were validated by estimating the temporal and spatial convergence and reliability of the obtained solution was validated with the Lyapunov exponents obtained by qualitatively different methods. A few case studies related to the loss of stability, the magnitude of the size-dependent parameter, the type and intensity of the temperature input, and the parameters of uniformly distributed transverse load were investigated.     

Comparison of different plasticity criteria for trabecular bone failure modelling

Different plasticity criteria applied to failure analysis of trabecular bone are compared. A cylindrical sample of bovine trabecular bone is mechanically tested in uniaxial compression/tension with 2% applied strain. Obtained response in compared to responses obtained using finite element model of trabecular bone inner structure subjected to the same loading conditions. FE model is reconstructed from micro–CT images. Elastic material properties at the level of trabecula are determined using nanoindentation. Compared plasticity criteria are based on these elastic material properties, i.e. on Young's modulus of elasticity from nanoindentation. The objective of the paper is to demonstrate the importance to reflect the anisotropic plasticity and to evaluate variation in obtained response when the anisotropy is neglected. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

正交各向异性薄板弯曲问题分裂模量有限元法

讨论了建立分裂模量有限元法的必要性,推导了正交各向异性薄板弯曲问题分裂模量变分原理的泛函,以此为基础建立了该问题的分裂模量有限元法。该模型的特点是其中含有一个被称为分裂因子的参数,通过算例说明:适当调整分裂因子的值,可以达到调整有限元模型的刚度、降低有限元刚度矩阵的谱条件数、克服常规有限元病态问题的目的,最后分析了克服病态问题的机理。     

不同模量弯压柱的解析解

依据不同模量弹性理论用流动坐标系及分段积分法导出复合荷载作用下弯压柱的解析解,建立了中性轴、应力、应变、位移的计算公式,并编制相应的有限元程序进行计算,与解析解进行误差对比．最后对不同模量计算结果与经典力学同模量计算结果进行分析对比,得出两种理论计算结果的差异,并提出对该类结构计算的合理建议．     

Exact relations for effective tensors of composites: Necessary conditions and sufficient conditions

Typically the elastic and electrical properties of composite materials are strongly microstructure dependent. So it comes as a nice surprise to come across exact formulae for effective moduli that are universally valid no matter what the microstructure. Such exact formulae provide useful benchmarks for testing numerical and actual experimental data and for evaluating the merit of various approximation schemes. They can also be regarded as fundamental invariances existing in a given physical context. Classic examples include Hill's formulae for the effective bulk modulus of a two‐phase mixture when the phases have equal shear moduli, Levin's formulae linking the effective thermal expansion coefficient and effective bulk modulus of two‐phase mixtures, and Dykhne's result for the effective conductivity of an isotropic two‐dimensional polycrystalline material. Here we present a systematic theory of exact relations embracing the known exact relations and establishing new ones. The search for exact relations is reduced to a search for matrix subspaces having a structure of special Jordan algebras. One of many new exact relations is for the effective shear modulus of a class of three‐dimensional polycrystalline materials. We present complete lists of exact relations for three‐dimensional thermoelectricity and for three‐dimensional thermopiezoelectric composites that include all exact relations for elasticity, thermoelasticity, and piezoelectricity as particular cases. © 2000 John Wiley & Sons, Inc.     

A linear elastic Fuzzy Finite Element Method with two fuzzy input parameters

This work employs a fuzzy finite element method (FFEM) in order to characterize the non-statistical properties based on the possibility theory [2]. A linear elastic body with fuzzy input material parameters Young's modulus and Poisson's ratio is studied. To this end, the fuzzy parameters are represented by possibility distributions generated from sparse experimental data. A computational approach involving the α-level discretization technique [4] is used in order to calculate the possibility of the system response. Finally, our method is applied in a numerical example. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Parameter identification based on FE‐models by genetic algorithms

A realistic and reliable model is an important precondition for the simulation of revitalization tasks as well as for the estimation of properties of existing buildings. Within one theory the parameters of the model should be approximated best by gradually performed experiments and their analysis. Usually this kind of optimization problems leads into non‐convex non‐differentiable objective function spaces with high dimensions. Normally ore complex structures are modeled using finite element method. We present a method of identifying Young's modulus for a beam and a plate by using FE‐models and genetic optimization algorithms for parameter identification.     

Finite element based nonlinear dynamic response of elastically supported piezoelectric functionally graded beam subjected to moving load in thermal environment with random system properties

The second order statistics in terms of mean and standard deviation (SD) of normalized nonlinear transverse dynamic central deflection (NTDCD) response of un-damped elastically supported functionally graded materials (FGMs) beam with surface-bonded piezoelectric layers under the action of moving load are investigated in this paper. The random system properties such as Young's modulus, Poisson's ratio, density, thermal expansion coefficients, piezoelectric materials, volume fraction exponent and external loading are modeled as uncorrelated random variables. The basic formulation is based on higher order shear deformation theory (HSDT) with von-Karman nonlinear strain kinematics combined with Newton–Raphson technique through Newmark's time integrating scheme using finite element method (FEM). The non-uniform temperature distribution with temperature dependent material properties is taken into consideration for consideration of thermal loading. The one parameter Pasternak elastic foundation with Winkler cubic nonlinearity is considered as an elastic foundation. The stochastic based second order perturbation technique (SOPT) and direct Monte Carlo simulation (MCS) are adopted for the solution of nonlinear dynamic governing equation. The influences of volume fraction exponents, temperature increments, moving loads and velocity, nonlinearity, slenderness ratios, foundation parameters and external loadings with random system properties on the NTDCD are examined. The capability of present stochastic model in predicting the NTDCD statistics are compared by studying their convergence with the existing results those available in the literature.     

Estimation of material properties of cancellous bone using multiscale FEM

Cancellous bone is a spongy type of bone with voids filled by blood marrow. Without much loss of generality it can be modeled as a material with periodic microstructure where overall parameters can be calculated using homogenization methods. Here the multiscale finite element method is applied and the assumed representative volume element (RVE) is a cube with solid frame and fluid core. From the point of view of the finite element method the RVE is a combination of solid and shell elements. As the acoustic excitation is considered, a complex stiffness matrix and complex displacements appear in the solution of the problem. Calculation of overall properties is repeated for different geometries of the solid frame. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

基于变体积约束的阻尼材料微结构拓扑优化研究北大核心CSCD

阻尼复合结构的抑振性能取决于材料布局和阻尼材料特性.该文提出了一种变体积约束的阻尼材料微结构拓扑优化方法,旨在以最小的材料用量获得具有期望性能的阻尼材料微结构.基于均匀化方法,建立阻尼材料三维微结构有限元模型,得到阻尼材料的等效弹性矩阵.逆用Hashin-Shtrikman界限理论,估计对应于期望等效模量的阻尼材料体积分数限,并构建阻尼材料体积约束限的移动准则.将获得阻尼材料微结构期望性能的优化问题转化为体积约束下最大化等效模量的优化问题,建立阻尼材料微结构的拓扑优化模型.利用优化准则法更新设计变量,实现最小材料用量下的阻尼材料微结构最优拓扑设计.通过典型数值算例验证了该方法的可行性和有效性,并讨论了初始微构型、网格依赖性和弹性模量等对阻尼材料微结构的影响.     

A Two Phase Mixture Model based on Bone Observation

An optimization method is developed to describe the mechanical behaviour of the human cancellous bone. The method is based on a mixture theory. A careful observation of the behaviour of the bone material leads to the hypothesis that the bone density is controlled by the principal stress trajectories (Wolff's law). The basic idea of the developed method is the coupling of a scalar value via an eigenvalue problem to the principal stress trajectories. On the one hand this theory will permit a prediction of the reaction of the biological bone structure after the implantation of a prosthesis, on the other hand it may be useful in engineering optimization problems. An analytical example shows its efficiency. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Extended multiscale FEM for design of beams and frames with complex topology

A numerical method for design of beams and frames with complex topology is proposed. The method is based on extended multi-scale finite element method where beam finite elements are used on coarse scale and continuum elements on fine scale. A procedure for calculation of multi-scale base functions, up-scaling and downscaling techniques is proposed by using a modified version of window method that is used in computational homogenization. Coarse scale finite element is embedded into a frame of a material that is representing surrounding structure in a sense of mechanical properties. Results show that this method can capture displacements, shear deformations and local stress-strain gradients with significantly reduced computational time and memory comparing to full scale continuum model. Moreover, this method includes a special hybrid finite elements for precise modelling of structural joints. Hence, the proposed method has a potential application in large scale 2D and 3D structural analysis of non-standard beams and frames where spatial interaction between structural elements is important.     

Bivariate spline method for numerical solution of steady state Navier–Stokes equations over polygons in stream function formulation

We use the bivariate spline finite elements to numerically solve the steady state Navier–Stokes equations. The bivariate spline finite element space we use in this article is the space of splines of smoothness r and degree 3r over triangulated quadrangulations. The stream function formulation for the steady state Navier–Stokes equations is employed. Galerkin's method is applied to the resulting nonlinear fourth‐order equation, and Newton's iterative method is then used to solve the resulting nonlinear system. We show the existence and uniqueness of the weak solution in H²(Ω) of the nonlinear fourth‐order problem and give an estimate of how fast the numerical solution converges to the weak solution. The Galerkin method with C¹ cubic splines is implemented in MATLAB. Our numerical experiments show that the method is effective and efficient. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 147–183, 2000     

The material parameter design and finite element simulation of the quadrilateral thermal cloak device

In this paper, we first design a coordinate transformation and derive the anisotropic material parameters of the quadrilateral thermal cloak according to the transformation thermodynamics principle. Then, since the derived parameters are inherently anisotropic, we eliminate its anisotropy by considering the effective medium theory and use a layered structure of metamaterials composed of only two isotropic materials to design the cloak device. Finally, we simulate the performance of a perfect and layered thermal cloak by the finite element method. To the best of our knowledge, this is the first work to design and simulate the performance of this quadrilateral thermal cloak by the finite element method(FEM).     

Rigid‐plastic hybrid element analyses of the plane strain upsetting

A rigid‐plastic hybrid element method (HEM) for simulation of metal forming is developed. This method is a mixed approach of the rigid‐plastic domain‐BEM and the rigid‐plastic FEM based on the theory of compressible plasticity. Because the compatibilities of not only velocity but also velocity's derivative between the adjoining boundary elements and finite elements can be met, the velocities and the derivatives of the velocity can be calculated with the same precision for the rigid‐plastic HEM. Then, it is considered that the rigid‐plastic HEM is a more precise method in formulation than the conventional rigid‐plastic FEMs for which the compatibilities of velocity's derivative cannot be met. The plane strain upsetting processes with two friction factors are analyzed by the rigid‐plastic HEM in this article. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 726–737, 2002; Published online in Wiley InterScience (www.interscience.wiley.com); DOI 10.1002/num.10031.     

An adaptive approach for modeling a fiber-matrix composite with the FE2 method

In materials with a complicated microstructre [1], the macroscopic material behaviour is unknown. In this work a Fiber-Matrix composite is considered with elasto-plastic fibers. A homogenization of the microscale leads to the macroscopic material properties. In the present work, this is realized in the frame of a FE² formulation. It combines two nested finite element simulations. On the macroscale, the boundary value problem is modelled by finite elements, at each integration point a second finite element simulation on the microscale is employed to calculate the stress response and the material tangent modulus. One huge disadvantage of the approach is the high computational effort. Certainly, an accompanying homogenization is not necessary if the material behaves linear elastic. This motivates the present approach to deal with an adaptive scheme. An indicator, which makes use of the different boundary conditions (BC) of the BVP on microscale, is suggested. The homogenization with the Dirichlet BC overestimates the material tangent modulus whereas the Neumann BC underestimates the modulus [2]. The idea for an adaptive modeling is to use both of the BCs during the loading process of the macrostructure. Starting initially with the Neumann BC leads to an overestimation of the displacement response and thus the strain state of the boundary value problem on the macroscale. An accompanying homogenization is performed after the strain reaches a limit strain. Dirichlet BCs are employed for the accompanying homogenization. Some numerical examples demonstrate the capability of the presented method. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Fractal dimension of a random invariant set

In recent years many deterministic parabolic equations have been shown to possess global attractors which, despite being subsets of an infinite-dimensional phase space, are finite-dimensional objects. Debussche showed how to generalize the deterministic theory to show that the random attractors of the corresponding stochastic equations have finite Hausdorff dimension. However, to deduce a parametrization of a ‘finite-dimensional’ set by a finite number of coordinates a bound on the fractal (upper box-counting) dimension is required. There are non-trivial problems in extending Debussche's techniques to this case, which can be overcome by careful use of the Poincaré recurrence theorem. We prove that under the same conditions as in Debussche's paper and an additional concavity assumption, the fractal dimension enjoys the same bound as the Hausdorff dimension. We apply our theorem to the 2d Navier–Stokes equations with additive noise, and give two results that allow different long-time states to be distinguished by a finite number of observations.     

Exponential Finite Elements for a Phase Field Fracture Model

Sharp interface material models can be related to phase field models by introducing an order parameter, whose value is assigned to the different phases of a material. The elastic material law is coupled to the evolution equation of the order parameter and cracking is addressed as a phase transition problem instead of a moving boundary value problem. A regularization parameter ϵ controls the width of the diffuse cracks represented by the order parameter and the underlying sharp interface model can be recovered from the phase field model by the limit ϵ → 0. However, in numerical simulations using standard finite elements with linear shape functions, the minimum value of ϵ is restricted by the grid size and therefore the discretization of the crack field requires extensive mesh refinement for small values of ϵ. In this work, we construct special 2d shape functions which take into account the exponential character of the crack field and its dependence on the parameter ϵ. Especially in simulations with small values of ϵ and a rather coarse mesh, the elements with exponential shape functions perform significantly better than standard linear elements. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Effective bulk moduli of materials containing stochastically distributed nano-inhomogeneities with surface stresses

In the present contribution, a mathematical model for the investigation of the effective properties of a material with randomly distributed nano-particles is proposed. The surface effect is introduced via Gurtin-Murdoch equations describing properties of the matrix/nano-particle interface. They are added to the system of stochastic differential equations formulated within the framework of linear elasticity. The homogenization problem is reduced to finding a statistically averaged solution of the system of stochastic differential equations. These equations are based on the fundamental equations of linear elasticity, which are coupled with surface/interface elasticity accounting for the presence of surface tension. Using Green's function this system is transformed to a system of statistically non-linear integral equations. It is solved by the method of conditional moments. Closed-form expressions are derived for the effective moduli of a composite consisting of a matrix with randomly distributed spherical inhomogeneities. The radius of the nano-particles is included in the expression for the bulk moduli. As numerical examples, nano-porous aluminum and nano-porous gold are investigated assuming that only the influence of the interface effects on the effective bulk modulus is of interest. The dependence of the normalized bulk moduli of nano-porous aluminum on the pore volume fraction (for certain radii of nano-pores) are compared to and discussed in the context of other theoretical predictions. The effective Young's modulus of nano-porous gold as a function of pore radius (for fixed void volume fraction) and the normalized Young's modulus vs. the pore volume fraction for different pore radii are analyzed. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)