Latest ArticlesTo study the influence of steel truss web shear deformation on the deflection of steel truss web composite box beams, the beams were first decomposed into a laminated structure composed of top and bottom flanges and a steel truss. A steel truss web shear deformation angle function was introduced to establish an analytical model, and the flexural deformation of a simply supported beam was analyzed as an example. The effective stiffness of the cross-section was determined by combining Euler beam theory and the analytical solution, and the mid-span deflection was calculated using this effective stiffness. The flexural characteristics under different load conditions were analyzed and compared with the Euler beam theory. The influence of structural parameters such as steel truss web diameter, steel truss web wall thickness, and steel truss web inclination angle on the effective stiffness was also examined. The results show that considering steel truss web shear deformation provides an analytical solution closer to the finite element results, with a maximum error of 6.24%. Using the effective stiffness can effectively predict the mid-span deflection, with a maximum error of 3.64% compared with the analytical solution. Among the structural parameters affecting the effective stiffness, steel truss web wall thickness has the greatest influence, followed by steel truss web diameter and steel truss web inclination angle. Additionally, the effective stiffness is positively correlated with steel truss web diameter and steel truss wall thickness but negatively correlated with steel truss web inclination angle.
The calculation of the probability distribution of performance functions is a core issue in uncertainty quantification and reliability design, and the recently proposed equivalent expectation method (EEM) is an effective way to solve this problem. This paper proposes an improved EEM. Putting forward an empirical calculation formula for the standard deviation coefficient of auxiliary random variablesand obtaining a more accurate probability distribution of the auxiliary function. Meanwhile, aiming at the accuracy issue in calculating the probability distribution of theperformance function is proposed, the calculation formula for the PDF of the performance function is derived using only one auxiliary function. In the process of calculating the PDF, proposing an exact theoretical transformation of probability distribution from auxiliary functions to performance functions is proposed, resulting in a more accurate PDF of the performance function. Finally, the effectiveness and accuracy of the method are verified through three numerical examples. The results indicate that this method is suitable for computing the probability distribution of high-dimensional nonlinear or implicit performance functions.
Structural condition assessment is crucial for ensuring the safe services of structures, with structural damage detection (SDD) being a core component. In this paper, a novel SDD method is proposed based on the adaptive grasshopper algorithm and sparse regularization. It aims to tackle accuracy decline of SDD results and instability involving uncertainties and incomplete measurement, thereby achieving sparse-regularization-based structural condition assessment. Firstly, adaptive Lévy flight and elite opposition-based learning strategies are incorporated into the adaptive grasshopper algorithm to prevent the SDD process from falling into local optima and to enhance the stability of SDD results. Secondly, a modal parameter-based objective function with sparse regularization is formulated to increase the sparsity of SDD results, thereby improving SDD accuracy and robustness. The optimization results of competition-based evolutionary computation benchmark functions show that the adaptive grasshopper algorithm exhibits better global convergence and identification stability compared with its standard version. Numerical and experimental results for simply-supported beams indicate that the proposed method can ensure reliable SDD accuracy even in the case of incomplete measurements, and it possesses good noise robustness as well.
Small cross-sections of boundary elements easily induce the “internal tension” phenomenon of steel plate shear walls, making it difficult to fully utilize the seismic performance of buckling-restrained steel plate shear walls. The design of cross-sections of boundary elements is related to their internal force requirements, and analyzing the internal force requirements of boundary columns is meaningful. Based on the proposed buckling-restrained steel plate shear wall with multi-concrete panels (MBRSPSW), the analytical expressions for the internal force of the boundary column of the MBRSPSW were theoretically derived in this paper. Combined with the experimental research on the buckling-restrained steel plate shear wall horizontally assembled multi-concrete panels (H-MBRSPSW), its numerical model was established and verified. The internal force distributions of the boundary column obtained from numerical analysis and analytical calculation were further compared. The research results indicate that the inner steel plates in the MBRSPSW are divided into constrained regions and unconstrained regions. The axial force, shear force, and bending moment distributions of the boundary column from the analytical calculation results agree with those from the numerical analysis results, and the internal force calculation expressions of the boundary column are correct. The research results can provide a reference for the design of this category of steel plate shear walls.
Accurately simulating the impact characteristics of a dam break is of paramount significance for the prediction and mitigation of dam-break flow disasters. The B-spline material point method (BSMPM), as an improved algorithm of the material point method (MPM), effectively enhances computational accuracy and improves convergence. However, the BSMPM solves the governing equations based on a tensor grid rather than an Eulerian background grid. Moreover, its interpolation shape functions have a larger influence domain. Consequently, when solving problems involving fluid-structure coupling and contact, issues such as premature contact, difficulty in capturing contact interfaces, and challenges in calculating contact forces arise. Within the same tensor grid space of the BSMPM, accurate capture of contact interfaces is achieved based on the relative velocities and unit outward normals of the same nodes; employing the Greville Abscissa enables precise contact of contacting objects, thereby avoiding premature or spurious contact; through the Lagrange multiplier method, interface contact forces are accurately determined, thus constructing a high-precision contact algorithm for fluid-structure strongly coupled problems, facilitating research on the simulation of dam-break fluid impact with rigid and elastic obstacles, and enabling a comparison with existing experimental or simulated results. The results demonstrate that the simulated impact loads and structural deformation evolution patterns correspond well with existing experimental/simulated results. For rigid obstacles, the peak impact pressure exhibits concave parabolic growth and positive correlation exponential function growth with increasing water level and dam-break slope, respectively. For elastic obstacles, the peak impact pressure decreases exponentially with the increase in the height of the probing point. The feasibility and effectiveness of simulating dam-break flow impact problems using the BSMPM contact algorithm are validated, providing a new perspective for simulating dam-break flow impact problems.
Research on acoustic propagation in multiple fluids has important application values in naval architecture and ocean engineering, such as sound propagation in pipelines filled with water and air, and the detection of buried objects. There are two difficulties in solving such problems with the use of the classical finite element method: one is the serious numerical dispersion error in the finite element solutions under medium and high wave numbers; the other is the need to use refined mesh grids to discretize the fluids near the coupling interface. These difficulties lead to a large computational cost for the finite element method, and the manual intervention to generate refined grids. Compared with the finite element method, the weak-form meshfree method does not require traditional grids, and the dispersion error effect in its solution is much weaker, ensuring good computational accuracy and efficiency. However, the meshfree shape functions are usually discontinuous in the problem domain, resulting in the inability of the continuity condition of the acoustic particle velocity to be naturally satisfied on the interface. Therefore, this paper uses the penalty function method to reconstruct the continuity condition of the acoustic particle velocity on the interface, and proposes a Galerkin weak form suitable for meshfree methods for sound propagation in multiple fluids. Numerical analysis shows that the meshfree solutions is consistent with the reference solutions, and the computational accuracy and efficiency of the meshfree method can be higher than the finite element solutions.
Proppant transport in fractures is essentially a dense granular flow in a slot-shaped space. Applying the two-fluid method in numerical simulations of a field-scale particle flow is promising, but existing solid stress models cannot accurately describe the process of proppant accumulation. In this paper, the morphological change of a proppant pack under flow erosion was analyzed experimentally, and the important influence of cohesion on the change of the pack state was pointed out. Then, combined with the simulated results of a proppant transport and the results of a suspension apparent viscosity test, the influence of the particle radial distribution function on the solid kinetic pressure and the change trend of total solid pressure were analyzed, and the change rate of the solid friction pressure with the particle volume fraction was determined. Based on the granular matter theory and results of a direct shear test, the cohesion of the proppant pack was considered in the frictional viscosity model. The results show that the improved solid friction stress model can capture larger angles of the accumulation and settlement profiles, and correctly simulate the process of proppant accumulation.
This paper presents a Bezier triangle meshing method that considers both clipped and non-clipped forms for a single NURBS surface. The proposed method is applied to analyze isogeometric Kirchhoff-Love shell structures. The process begins by interpolating NURBS surfaces into Bezier surfaces. Subsequently, the topological relationship between the clipping curve and each parameter node is calculated within the parameter domain. A Bezier contour curve set is then generated in the parameter domain by selecting points along the clipping curve. Utilizing this contour curve set, a triangular mesh is generated in the parameter domain. Finally, the Bezier triangle mesh in the physical domain is created through a mapping method. The adaptability and robustness of the algorithm are verified through three models, and the mesh quality is assessed. The results demonstrate favorable overall mesh quality. Building upon this foundation, the paper illustrates the application of a rotation constraint between Kirchhoff-Love shell elements, using the penalty function method with Scordelis-Lo's Roof shell model as an example. The accuracy of Kirchhoff-Love shell elements based on Bezier triangles is subsequently validated.
A surrogate model is a new research idea direction of aerodynamic data generation. A traditional surrogate model relies on a large number of high-precision simulation sample points and their responses values to ensure the accuracy of the model. A multi-credibility surrogate model can reduce the computational cost while maintaining a certain accuracy by integrating multi-layer high and low credibility models, which is of great significance for reducing a missile development cycle. In this paper, the influence of different numbers of low-confidence sample points on the multi-confidence model and the optimal ratio of high-low-confidence sample points are studied, and a multi-confidence sampling method suitable for aerodynamics data is proposed. It is applied to the construction of three multi-credibility surrogate models in the prediction of missile aerodynamic data, among which the Co-Kriging model has the best comprehensive prediction effect. The recommended ratio of high and low confidence sample size is between 1: 4 and 1: 3.
A simple numerical implementation method is proposed for the Chaboche-type viscoplastic constitutive model coupled with Lemaitre anisotropic damage theory. Using the decoupled algorithm, the damage tensor is updated based on the forward difference format at the beginning of each incremental step. The damage tensor is considered as a constant in the discretization process of the constitutive equations. Based on the hypothesis of strain equivalence, the formulations containing only partial tensors are constructed in the effective deviatoric stress space, and the radial return process is simplified to solve a nonlinear scalar equation concerning the accumulated plastic strain increment. The numerical implementation method and the derivation of consistent tangent operator are provided based on the Voigt notation scheme. The comparison between the experimental data and the simulation results of isotropic scalar damage model under uniaxial and multiaxial stress states validates the effectiveness and high computational efficiency of this method. Numerical results under different time step sizes also indicate the good accuracy and stability.