Generalized polynomial chaos expansion by reanalysis using static condensation based on substructuring

This paper presents a new computational method for forward uncertainty quantification(UQ)analyses on large-scale structural systems in the presence of arbitrary and dependent random inputs.The method consists of a generalized polynomial chaos expansion(GPCE)for statistical moment and reliability analyses associated with the stochastic output and a static reanalysis method to generate the input-output data set.In the reanalysis,we employ substructuring for a structure to isolate its local regions that vary due to random inputs.This allows for avoiding repeated computations of invariant substructures while generating the input-output data set.Combining substructuring with static condensation further improves the computational efficiency of the reanalysis without losing accuracy.Consequently,the GPCE with the static reanalysis method can achieve significant computational saving,thus mitigating the curse of dimensionality to some degree for UQ under high-dimensional inputs.The numerical results obtained from a simple structure indicate that the proposed method for UQ produces accurate solutions more efficiently than the GPCE using full finite element analyses(FEAs).We also demonstrate the efficiency and scalability of the proposed method by executing UQ for a large-scale wing-box structure under ten-dimensional(all-dependent)random inputs.

Grouping response method for equivalent static wind loads based on a modified LRC method

Wind loading is one of the most important loads for controlling the design of large-span roof structures. Equivalent static wind loads, which can generally aim at determining a specific response, are widely used by structural designers. A method for equivalent static wind loads applicable to multi-responses is proposed in this paper. A modified load- response-correlation （LRC） method corresponding to a particular peak response is presented, and the similarity algorithm implemented for the group response is described. The main idea of the algorithm is that two responses can be put into one group if the value of one response is close to that of the other response, when the structure is subjected to equivalent static wind loads aiming at the other response. Based on the modified LRC, the grouping response method is put forward to construct equivalent static wind loading. This technique can simultaneously reproduce peak responses for some grouped responses. To verify its computational accuracy, the method is applied to an actual large-span roof structure. Calculation results show that when the similarity of responses in the same group is high, equivalent static wind loads with high accuracy and reasonable magnitude of equivalent static wind distribution can be achieved.

Method of reverberation ray matrix for static analysis of planar framed structures composed of anisotropic Timoshenko beam members

Based on the method of reverberation ray matrix（MRRM）, a reverberation matrix for planar framed structures composed of anisotropic Timoshenko（T） beam members containing completely hinged joints is developed for static analysis of such structures.In the MRRM for dynamic analysis, amplitudes of arriving and departing waves for joints are chosen as unknown quantities. However, for the present case of static analysis, displacements and rotational angles at the ends of each beam member are directly considered as unknown quantities. The expressions for stiffness matrices for anisotropic beam members are developed. A corresponding reverberation matrix is derived analytically for exact and unified determination on the displacements and internal forces at both ends of each member and arbitrary cross sectional locations in the structure. Numerical examples are given and compared with the finite element method（FEM） results to validate the present model. The characteristic parameter analysis is performed to demonstrate accuracy of the present model with the T beam theory in contrast with errors in the usual model based on the Euler-Bernoulli（EB） beam theory. The resulting reverberation matrix can be used for exact calculation of anisotropic framed structures as well as for parameter analysis of geometrical and material properties of the framed structures.

Static and free vibration analyses of functionally graded porous variable-thickness plates using an edge-based smoothed finite element method

The main purpose of this paper is to present numerical results of static bending and free vibration of functionally graded porous(FGP) variable-thickness plates by using an edge-based smoothed finite element method(ES-FEM) associate with the mixed interpolation of tensorial components technique for the three-node triangular element(MITC3), so-called ES-MITC3. This ES-MITC3 element is performed to eliminate the shear locking problem and to enhance the accuracy of the existing MITC3 element. In the ES-MITC3 element, the stiffness matrices are obtained by using the strain smoothing technique over the smoothing domains formed by two adjacent MITC3 triangular elements sharing an edge. Materials of the plate are FGP with a power-law index(k) and maximum porosity distributions(U) in the forms of cosine functions. The influences of some geometric parameters, material properties on static bending, and natural frequency of the FGP variable-thickness plates are examined in detail.

An improved outer pipe method for expansive pressure measurement of static cracking agents

Static cracking agent(SCA)is actively investigated as an alternative to explosive blasting for rock breakage due to its immense expansion property.SCA can eliminate the negative effects of shock,noise and harmful gases encountered in explosive blasting processes.Accurate measurement and deep understanding of the expansive properties of SCAs are important in their industrial application.An improved outer pipe method(OPM),termed the upper end surface method(UESM),is proposed in this paper to overcome the shortcomings of the OPM in the expansive pressure measurement of SCAs.Numerical simulation is used to proof the concept and a mathematical model established to relate the internal pressure and the radial strains at different positions in the upper end surface method test equipment.The new equipment is calibrated using oil pressure and strain measurements.The calibrated equipment is then used to measure the expansion pressure of SCA at three different water contents to proof its potential.The differences in the measurements with OPM and UESM at three different moisture contents are less than 4%.The experimental results confirm the accuracy and applicability of the more user friendly and less expensive UESM in the measurement of the expansive pressures of SCAs.

An improved pseudo-static method for seismic resistant design of underground structures

This paper describes a commonly used pseudo-static method in seismic resistant design of the cross section of underground structures. Based on dynamic theory and the vibration characteristics of underground structures, the sources of errors when using this method are analyzed. The traditional seismic motion loading approach is replaced by a method in which a one-dimensional soil layer response stress is differentiated and then converted into seismic live loads. To validate the improved method, a comparison of analytical results is conducted for internal forces under earthquake shaking of a typical shallow embedded box-shaped subway station structure using four methods： the response displacement method, finite element response acceleration method, the finite element dynamic analysis method and the improved pseudo-static calculation method. It is shown that the improved finite element pseudo-static method proposed in this paper provides an effective tool for the seismic design of underground structures. The evaluation yields results close to those obtained by the finite element dynamic analysis method, and shows that the improved finite element pseudo-static method provides a higher degree of precision.

Static output feedback stabilization for second-order singular systems using model reduction methods

In this paper,the static output feedback stabilization for large-scale unstable second-order singular systems is investigated.First,the upper bound of all unstable eigenvalues of second-order singular systems is derived.Then,by using the argument principle,a computable stability criterion is proposed to check the stability of secondorder singular systems.Furthermore,by applying model reduction methods to original systems,a static output feedback design algorithm for stabilizing second-order singular systems is presented.A simulation example is provided to illustrate the effectiveness of the design algorithm.

Assessment of CH4 and C 2 surface emissions from Polesgo’s landfill(Ouagadougou, Burkina Faso) based on static chamber method

Methane (CH4) and carbon dioxide (C02) surface emissions from Polesgo's landfill (Ouagadougou, Burkina Faso) were measured using the static chamber technique in 2017 and 2018. The Polesgo's landfill was composed of four zones: Phase I, II, Phase III, and SP. The surface of Phase I was fully covered and its conditions are better for surface emission measurements. As results concerning the Phase I zone, the geospatial means flux rates of CH4 (657 mg m-2 h l in 2017 and 1210 mg m 2 h_, in 2018, respectively) are measured higher than the tolerable value reported in literature. The emitted CH4 or C 02 have permitted to locate higher surface emissions which are related to the cover state. The calculated gas collection efficiency (27.4% in 2017 and 23.0% in 2018) is low compared to those reported for landfills integrating landfill gas (LFG) extraction system. The carbon footprint calculations (24,966 tC02-eq 2017 and 40,025 tC02-eq in 2018, respectively) shown that Polesgo's landfill is a significant source of greenhouse gases (GHG) and its important potential for organic recovery can contribute to reduce the carbon footprint.

Three-Dimensional Static Analysis of Nanoplates and Graphene Sheets by Using Eringen’s Nonlocal Elasticity Theory and the Perturbation Method

A three-dimensional(3D)asymptotic theory is reformulated for the static analysis of simply-supported,isotropic and orthotropic single-layered nanoplates and graphene sheets(GSs),in which Eringen’s nonlocal elasticity theory is used to capture the small length scale effect on the static behaviors of these.The perturbation method is used to expand the 3D nonlocal elasticity problems as a series of two-dimensional(2D)nonlocal plate problems,the governing equations of which for various order problems retain the same differential operators as those of the nonlocal classical plate theory(CST),although with different nonhomogeneous terms.Expanding the primary field variables of each order as the double Fourier series functions in the in-plane directions,we can obtain the Navier solutions of the leading-order problem,and the higher-order modifications can then be determined in a hierarchic and consistent manner.Some benchmark solutions for the static analysis of isotropic and orthotropic nanoplates and GSs subjected to sinusoidally and uniformly distributed loads are given to demonstrate the performance of the 3D nonlocal asymptotic theory.

Static Analysis of Doubly-Curved Shell Structures of Smart Materials and Arbitrary Shape Subjected to General Loads Employing Higher Order Theories and Generalized Differential Quadrature Method

The article proposes an Equivalent Single Layer(ESL)formulation for the linear static analysis of arbitrarily-shaped shell structures subjected to general surface loads and boundary conditions.A parametrization of the physical domain is provided by employing a set of curvilinear principal coordinates.The generalized blendingmethodology accounts for a distortion of the structure so that disparate geometries can be considered.Each layer of the stacking sequence has an arbitrary orientation and is modelled as a generally anisotropic continuum.In addition,re-entrant auxetic three-dimensional honeycomb cells with soft-core behaviour are considered in the model.The unknown variables are described employing a generalized displacement field and pre-determined through-the-thickness functions assessed in a unified formulation.Then,a weak assessment of the structural problem accounts for shape functions defined with an isogeometric approach starting fromthe computational grid.Ageneralizedmethodology has been proposed to define two-dimensional distributions of static surface loads.In the same way,boundary conditions with three-dimensional features are implemented along the shell edges employing linear springs.The fundamental relations are obtained from the stationary configuration of the total potential energy,and they are numerically tackled by employing the Generalized Differential Quadrature(GDQ)method,accounting for nonuniform computational grids.In the post-processing stage,an equilibrium-based recovery procedure allows the determination of the three-dimensional dispersion of the kinematic and static quantities.Some case studies have been presented,and a successful benchmark of different structural responses has been performed with respect to various refined theories.

General Expression of LRC Method in Estimation of Bridge’s Equivalent Static Wind Load

According to the relationship between load and response, the equivalent static wind load(ESWL) of a structure can be estimated by load-response correlation(LRC) method, which can be accurately used to estimate the background ESWL of a structure. The derivation of the classical expression of LRC formula is based on a specific command response at a critical position, and the ESWL distribution has only one form in this case. In this paper, a general expression of LRC formula is derived based on a specific command response at all positions. For the general expression, ESWLs can be expressed by load-response correlation coefficients, response-response correlation coefficients, RMS values of the fluctuating wind loads, and peak factor in the form of matrices. By comparing the expressions of LRC method, it was found that the classical expression was only one form of the general one. The general expression which introduces the response-response correlation coefficients provided more options for structural engineers to estimate ESWLs and offered further insights into the LRC method. Finally, a cable-stayed bridge, a rigid three span continuous girder bridge, and a suspension bridge were used to verify the correctness of the general expression of LRC method.

Analysis of the Results Obtained from the Application of the Two-Stage Method with Calculations of Some Statically Indeterminate Trusses

The paper presents results of calculations of forces in members of selected types of statically indeterminate trusses carriedout by application of the two-stage method of computations of such structural systems. The method makes possible to do the simple andapproximate calculations of the complex trusses in two stages, in each of which is calculated a statically determinate truss being anappropriate counterpart of the basic form of the statically indeterminate truss structure. Systems of the statically determinate trussesconsidered in the both stages are defined by cancelation of members, number of which is equal to the statically indeterminacy of thebasic truss. In the paper are presented outcomes obtained in the two-stage method applied for two different shapes of trusses and carriedout for various ways of removing of appropriate members from the basic trusses. The results are compared with outcomes gained due toapplication of suitable computer software for computation of the same types of trusses and for the same structural conditions.

A Modified Quasi-Static Method for CALM System Analysis

The quasi-static analysis method introduced by API RP 2P is well known and accepted as a very useful mooring analysis method. In the early design stage, this method is widely used for preliminary analysis and mooring parameter selection. However, the quasi-static method of API RP 2P is developed for single-floating-body condition, i. e., only one floating body is considered in the computation procedure. Difficulties arise when it is used for the analysis of a CALM system, which is comprised of two floating bodies (tanker and buoy). This paper presents an analysis procedure for a two-floating-body system based on the quasi-static procedure of API RP 2P with some modifications reflecting special characteristics of the CALM system. Finally, the analysis results of a CALM system are given to illustrate the use of this procedure.

Design Optimization for Generating a High Static Magnetic Field

This paper details the creation of a device capable of generating a powerful and consistent static magnetic field. This apparatus serves the purpose of quantifying the magnetostrictive strain found in materials like annealed cobalt ferrite and Terfenol-D, specifically those shaped as cylindrical rods. In our investigation, the use of static magnetic fields proves most advantageous. This choice is made to mitigate the generation of eddy currents, which would inevitably occur if the magnetic field intensity were varied. The fundamental idea behind this design involves employing a C-shaped iron core constructed from low-carbon mild steel. On this core, three coils are mounted, each capable of producing one-third of the required 9000 Oersted (Oe) magnetic field strength. The test specimen is situated within the “jaws” of the C-shaped core, thus completing the magnetic circuit. To manage the heat generated by each coil, a cooling system consisting of copper tubes is employed. These tubes facilitate the flow of air to dissipate the heat. To model and predict the magnetic field strength produced by the coils, finite element analysis (FEMM) software is utilized, and the results align closely with the anticipated outcomes. This design effectively generates a robust and unchanging magnetic field measuring a stable 9000 Oe in total. Consequently, this equipment finds utility in characterizing the magnetic properties of specific materials.

PARALLEL ADAPTIVELY MODIFIED CHARACTERISTIC BASIS FUNCTION METHOD BASED ON STATIC LOAD BALANCE

Characteristic Basis Function Method (CBFM) is a novel approach for analyzing the ElectroMagnetic (EM) scattering from electrically large objects. Based on dividing the studied object into small blocks, the CBFM is suitable for parallel computing. In this paper, a static load balance parallel method is presented by combining Message Passing Interface (MPI) with Adaptively Modified CBFM (AMCBFM). In this method, the object geometry is partitioned into distinct blocks, and the serial number of blocks is sent to related nodes according to a certain rule. Every node only needs to calculate the information on local blocks. The obtained results confirm the accuracy and efficiency of the proposed method in speeding up solving large electrical scale problems.

落石被动柔性防护网非线性静力设计方法及验证

被动柔性防护网是一种常见的落石灾害防治措施。针对当前落石被动柔性防护网尚无适用的工程化设计方法的现状,以非线性计算过程中结构内能与设计防护能级相等作为加载完成的判别准则,建立了落石被动柔性防护网的非线性静力设计方法。基于LS-DYNA对250~5000 kJ共8种常用能级的被动柔性防护网开展非线性动力分析,明确了冲击力的分布范围及相对大小,形成了适用于不同能级落石被动柔性防护网非线性静力设计的荷载模式。为了有效模拟被动柔性防护网特殊的受力特征,开发了由索单元与受拉弹簧单元组成的非线性耗能单元用于耗能器模拟,开发了内部由多段索单元组成的滑移索单元用于滑移边界以及网环单元模拟。基于对网片冲击试验的反演和3500 kJ被动柔性防护网的设计及足尺冲击试验检验,验证了所提出的落石被动柔性防护网非线性静力设计方法的有效性。

基于离散元法的面粉颗粒接触参数标定试验

为设计和优化面粉输送设备,应用离散元法对面粉进行准确地工程建模和分析,需要对其接触参数进行必要的标定。该研究依据颗粒缩放理论,用“Hertz-Mindlin with Johnson-Kendall-Roberts”接触模型表征面粉颗粒间黏性的影响,提出了一种基于静/动态休止角的接触参数标定方法。运用正交试验方法,对接触参数的敏感性和方差分析,表明面粉颗粒间的滚动摩擦系数、面粉颗粒与不锈钢表面间的静摩擦系数、表面能对静态休止角的影响极显著(P<0.01),并且多组接触参数都可以模拟出与试验相同的静态休止角。进一步研究表明,面粉颗粒与不锈钢表面间的静摩擦系数的合理取值范围为0.2~0.4。通过2种填充率、4种转速下基于动态休止角的参数标定,将其中与试验最为吻合的一组参数作为标定结果,其值如下:面粉颗粒之间恢复系数为0.6、面粉颗粒之间静摩擦系数为0.2、面粉颗粒之间滚动摩擦系数为0.1、面粉颗粒与不锈钢容器表面之间恢复系数为0.6、面粉颗粒与不锈钢容器表面之间静摩擦系数为0.6、面粉颗粒与不锈钢容器表面之间滚动摩擦系数为0.5、表面能为0.12 J/m^(2)。使用该组参数对矩形容器中物料自由坍塌试验进行仿真,其结果与试验结果相符,验证了该标定方法的有效性。该研究提出的标定方法简单、易执行,对粉料输送设备的设计及优化具有一定的工程应用价值。

考虑弹流润滑效应的三点接触球轴承刚度特性分析

为获得润滑状态下三点接触球轴承更为准确的刚度特性,应考虑弹流润滑效应对轴承刚度的影响。文中基于拟静力学模型考虑高速离心力和陀螺力矩效应,根据给定轴承的结构参数和工况,计算滚动体与内外圈的法向接触载荷和各部件的运动速度。将拟静力学模型的计算结果和润滑介质参数代入弹流润滑模型,求解出滚动体与内外圈之间的压力分布和油膜厚度分布。进一步研究了转速、轴向载荷和润滑油的初始黏度对油膜压力和最小油膜厚度的影响。基于弹流润滑理论分析了转速和轴向载荷对轴承接触刚度、油膜刚度及综合刚度的影响。结果表明:转速的提高会大幅增加润滑油膜的整体厚度;润滑油初始黏度的增大会增加油膜厚度;随着轴承转速的提高,轴承的整体轴向刚度和轴向油膜刚度减小;随着轴向载荷的增大,轴承轴向刚度和轴向油膜刚度增大,且差值变化不大。