Modeling Geometrically Nonlinear FG Plates: A Fast and Accurate Alternative to IGA Method Based on Deep Learning

Isogeometric analysis (IGA) is known to showadvanced features compared to traditional finite element approaches.Using IGA one may accurately obtain the geometrically nonlinear bending behavior of plates with functionalgrading (FG). However, the procedure is usually complex and often is time-consuming. We thus put forward adeep learning method to model the geometrically nonlinear bending behavior of FG plates, bypassing the complexIGA simulation process. A long bidirectional short-term memory (BLSTM) recurrent neural network is trainedusing the load and gradient index as inputs and the displacement responses as outputs. The nonlinear relationshipbetween the outputs and the inputs is constructed usingmachine learning so that the displacements can be directlyestimated by the deep learning network. To provide enough training data, we use S-FSDT Von-Karman IGA andobtain the displacement responses for different loads and gradient indexes. Results show that the recognition erroris low, and demonstrate the feasibility of deep learning technique as a fast and accurate alternative to IGA formodeling the geometrically nonlinear bending behavior of FG plates.

Geometrically Nonlinear Flutter Analysis Based on CFD/CSD Methods and Wind Tunnel Experimental Verification

This study presents a high-speed geometrically nonlinear flutter analysis calculation method based on the highprecision computational fluid dynamics/computational structural dynamics methods.In the proposed method,the aerodynamic simulation was conducted based on computational fluid dynamics,and the structural model was established using the nonlinear finite element model and tangential stiffness matrix.First,the equilibrium position was obtained using the nonlinear static aeroelastic iteration.Second,the structural modal under a steady aerodynamic load was extracted.Finally,the generalized displacement time curve was obtained by coupling the unsteady aerodynamics and linearized structure motion equations.Moreover,if the flutter is not at a critical state,the incoming flow dynamic pressure needs to be changed,and the above steps must be repeated until the vibration amplitude are equal.Furthermore,the high-speed geometrically nonlinear flutter of the wing-body assemblymodel with a high-aspect ratio was investigated,and the correctness of the method was verified using high-speed wind tunnel experiments.The results showed that the geometric nonlinearity of the large deformation of the wing caused in-plane bending to become a key factor in flutter characteristics and significantly decreased the dynamic pressure and frequency of the nonlinear flutter compared to those of the linear flutter.

New tangent stiffness matrix for geometrically nonlinear analysis of space frames

A three-dimensional beam element is derived based on the principle of stationary total potential energy for geometrically nonlinear analysis of space frames. A new tangent stiffness matrix, which allows for high order effects of element deformations, replaces the conventional incremental secant stiffness matrix. Two deformation stiffness matrices due to the variation of axial force and bending moments are included in the tangent stiffness. They are functions of element deformations and incorporate the coupling among axial, lateral and torsional deformations. A correction matrix is added to the tangent stiffness matrix to make displacement derivatives equivalent to the commutative rotational degrees of freedom. Numerical examples show that the proposed dement is accurate and efficient in predicting the nonlinear behavior, such as axial-torsional and flexural-torsional buckling, of space frames even when fewer elements are used to model a member.

基于实际预压载作业的自升式平台环境载荷图谱绘制新思路

The environmental load chart is an important technical support required for the jack-up drilling platform to facilitate its adaptation to different operating waters and ensure the safety of operation.This chart is a crucial part of the platform operation manual.The chart data are closely related to external factors such as water depth,wind,wave,and current conditions of the working water,as well as to the structural characteristics of the platform itself and the number of variable loads.This study examines the platform state under extreme wind,wave,and current conditions during preloading.In addition,this study focuses on the difference between the ultimate reaction force of the pile leg during preloading and the reaction force of the pile leg without considering any environmental load before preloading.Furthermore,the relationship between the difference and the new reaction force of the pile leg caused by the combination of different environmental conditions is established to facilitate the construction of a new form of environmental load chart.The newly formed chart is flexible and simple;thus,it can be used to evaluate the environmental adaptability of the platform in the target well location and provides the preloading target demand or variable load limit according to the given environmental constraints.Moreover,the platform can perform personalized preloading operations,thereby improving its capability to cope with complex geological conditions,such as reducing punch-through risks.This condition reduces the load on jacking system devices and increases its service life.

Numerical Exploration of Asymmetrical Impact Dynamics: Unveiling Nonlinearities in Collision Problems and Resilience of Reinforced Concrete Structures

This study provides a comprehensive analysis of collision and impact problems’ numerical solutions, focusing ongeometric, contact, and material nonlinearities, all essential in solving large deformation problems during a collision.The initial discussion revolves around the stress and strain of large deformation during a collision, followedby explanations of the fundamental finite element solution method for addressing such issues. The hourglassmode’s control methods, such as single-point reduced integration and contact-collision algorithms are detailedand implemented within the finite element framework. The paper further investigates the dynamic responseand failure modes of Reinforced Concrete (RC) members under asymmetrical impact using a 3D discrete modelin ABAQUS that treats steel bars and concrete connections as bond slips. The model’s validity was confirmedthrough comparisons with the node-sharing algorithm and system energy relations. Experimental parameterswere varied, including the rigid hammer’s mass and initial velocity, concrete strength, and longitudinal and stirrupreinforcement ratios. Findings indicated that increased hammer mass and velocity escalated RC member damage,while increased reinforcement ratios improved impact resistance. Contrarily, increased concrete strength did notsignificantly reduce lateral displacement when considering strain rate effects. The study also explores materialnonlinearity, examining different materials’ responses to collision-induced forces and stresses, demonstratedthrough an elastic rod impact case study. The paper proposes a damage criterion based on the residual axialload-bearing capacity for assessing damage under the asymmetrical impact, showing a correlation betweendamage degree hammer mass and initial velocity. The results, validated through comparison with theoreticaland analytical solutions, verify the ABAQUS program’s accuracy and reliability in analyzing impact problems,offering valuable insights into collision and impact problems’ nonlinearities and practical strategies for enhancingRC structures’ resilience under dynamic stress.

Numerical evaluation of buckling behavior in space structure considering geometrical parameters with joint rigidity

The buckling behavior of single layer space structure is very sensitive. The joint rigidity, moreover, is one of the main factors of stability which may determine the entire failure behavior. Thus, the reasonable stiffness of joint system, which is neither total pin assumption nor perfect fix condition, is very important to apply to the real single layer space one. Therefore, the purpose of this work was to investigate the buckling behavior of single layer space structure, using the development of the upgraded stiffness matrix for the joint rigidity. To derive tangential stiffness matrix, a displacement function was assumed using translational and rotational displacement at the node. The geometrical nonlinear analysis was simulated not only with perfect model but also with imperfect one. As a result, the one and two free nodal numerical models were investigated using derived stiffness matrix. It was figured out that the buckling load increases in proportion to joint rigidity with rise-span ratio. The stability of numerical model is very sensitive with the initial imperfection, responding of bifurcation in the structure.

Topology Optimization of Compliant Mechanisms with Geometrical Nonlinearities Using the Ground Structure Approach

The majority of topology optimization of compliant mechanisms uses linear finite element models to find the structure responses.Because the displacements of compliant mechanisms are intrinsically large,the topological design can not provide quantitatively accurate result.Thus,topological design of these mechanisms considering geometrical nonlinearities is essential.A new methodology for geometrical nonlinear topology optimization of compliant mechanisms under displacement loading is presented.Frame elements are chosen to represent the design domain because they are capable of capturing the bending modes.Geometrically nonlinear structural response is obtained by using the co-rotational total Lagrange finite element formulation,and the equilibrium is solved by using the incremental scheme combined with Newton-Raphson iteration.The multi-objective function is developed by the minimum strain energy and maximum geometric advantage to design the mechanism which meets both stiffness and flexibility requirements, respectively.The adjoint method and the direct differentiation method are applied to obtain the sensitivities of the objective functions. The method of moving asymptotes（MMA） is employed as optimizer.The numerical example is simulated to show that the optimal mechanism based on geometrically nonlinear formulation not only has more flexibility and stiffness than that based on linear formulation,but also has better stress distribution than the one.It is necessary to design compliant mechanisms using geometrically nonlinear topology optimization.Compared with linear formulation,the formulation for geometrically nonlinear topology optimization of compliant mechanisms can give the compliant mechanism that has better mechanical performance.A new method is provided for topological design of large displacement compliant mechanisms.

GEOMETRICALLY NONLINEAR FINITE ELEMENT MODEL OF SPATIAL THIN-WALLED BEAMS WITH GENERAL OPEN CROSS SECTION

Based on the theory of Timoshenko and thin-walled beams, a new finite element model of spatial thin-walled beams with general open cross sections is presented in the paper, in which several factors are included such as lateral shear deformation, warp generated by nonuni- form torsion and second-order shear stress, coupling of flexure and torsion, and large displacement with small strain. With an additional internal node in the element, the element stiffness matrix is deduced by incremental virtual work in updated Lagrangian （UL） formulation. Numerical examples demonstrate that the presented model well describes the geometrically nonlinear property of spatial thin-walled beams.

A MESHLESS LOCAL PETROV-GALERKIN METHOD FOR GEOMETRICALLY NONLINEAR PROBLEMS

Nonlinear formulations of the meshless local Petrov-Galerkin （MLPG） method are presented for geometrically nonlinear problems. The method requires no mesh in computation and therefore avoids mesh distortion difficulties in the large deformation analysis. The essential boundary conditions in the present formulation axe imposed by a penalty method. An incremental and iterative solution procedure is used to solve geometrically nonlinear problems. Several examples are presented to demonstrate the effectiveness of the method in geometrically nonlinear problems analysis. Numerical results show that the MLPG method is an effective one and that the values of the unknown variable are quite accurate.

Geometrical nonlinear stability analyses of cable-truss domes

The nonlinear finite element method is used to analyze the geometrical nonlinear stability of cable truss domes with different cable distributions. The results indicate that the critical load increases evidently when cables, especially diagonal cables, are distributed in the structure. The critical loads of the structure at different rise span ratios are also discussed in this paper. It was shown that the effect of the tensional cable is more evident at small rise span ratio. The buckling of the structure is characterized by a global collapse at small rise span ratio; that the torsional buckling of the radial truss occurs at big rise span ratio; and that at proper rise span ratio, the global collapse and the lateral buckling of the truss occur nearly simultaneously.

TOPOLOGY SYNTHESIS OF GEOMETRICALLY NONLINEAR COMPLIANT MECHANISMS USING MESHLESS METHODS

This paper presents a new method for topology optimization of geometrical nonlinear compliant mechanisms using the element-free Galerkin method （EFGM）. The EFGM is employed as an alternative scheme to numerically solve the state equations by fully taking advantage of its capability in dealing with large displacement problems. In the meshless method, the imposition of essential boundary conditions is also addressed. The popularly studied solid isotropic material with the penalization （SIMP） scheme is used to represent the nonlinear dependence between material properties and regularized discrete densities. The output displacement is regarded as the objective function and the adjoint method is applied to finding the sensitivity of the design functions. As a result, the optimization of compliant mechanisms is mathematically established as a nonlinear programming problem, to which the method of moving asymptotes （MMA） belonging to the sequential convex programming can be applied. The availability of the present method is finally demonstrated with several widely investigated numerical examples.

Geometrically Nonlinear Random Responses of Stiffened Plates Under Acoustic Pressure

An algorithm integrating reduced order model(ROM),equivalent linearization(EL),and finite element method(FEM)is proposed to carry out geometrically nonlinear random vibration analysis of stiffened plates under acoustic pressure loading.Based on large deflection finite element formulation,the nonlinear equations of motion of stiffened plates are obtained.To reduce the computation,a reduced order model of the structures is established.Then the EL technique is incorporated into FE software NASTRAN by the direct matrix abstraction program(DMAP).For the stiffened plates,a finite element model of beam and plate assembly is established,in which the nodes of beam elements are shared with shell elements,and the offset and section properties of the beam are set.The presented method can capture the root-mean-square(RMS) of the stress responses of shell and beam elements of stiffened plates,and analyze the stress distribution of the stiffened surface and the unstiffened surface,respectively.Finally,the statistical dynamic response results obtained by linear and EL methods are compared.It is shown that the proposed method can be used to analyze the geometrically nonlinear random responses of stiffened plates.The geometric nonlinearity plays an important role in the vibration response of stiffened plates,particularly at high acoustic pressure loading.

A new beam element for analyzing geometrical and physical nonlinearity

Based on Timoshenko＇s beam theory and Vlasov＇s thin-walled member theory, a new model of spatial thin-walled beam element is developed for analyzing geometrical and physical nonlinearity, which incorporates an interior node and independent interpolations of bending angles and warp and takes diversified factors into consideration, such as traverse shear deformation, torsional shear deformation and their coupling, coupling of flexure and torsion, and the second shear stress. The geometrical nonlinear strain is formulated in updated Lagarange （UL） and the corresponding stiffness matrix is derived. The perfectly plastic model is used to account for physical nonlinearity, and the yield rule of von Mises and incremental relationship of Prandtle-Reuss are adopted. Elastoplastic stiffness matrix is obtained by numerical integration based on the finite segment method, and a finite element program is compiled. Numerical examples manifest that the proposed model is accurate and feasible in the analysis of thin-walled structures.

Reliability analysis of circular concrete-filled steel tube with material and geometrical nonlinearity

Most of the previous research on concrete-filled steel tube is restricted to a deterministic approach. To gain clear insight into the random properties of circular concrete-filled steel tube, reliability analysis is carried out in the present study. To obtain the Structural nonlinear response and ultimate resistance capacity, material and geometrical nonlinear analysis of circular concrete-filled steel tube is performed with a three-dimensional degenerated beam ele- ment. Then we investigate the reliability of concrete-filled steel tube using the first-order reliability method combined with nonlinear finite element analysis. The influences of such parameters as material strength, slenderness, initial geo- metrical imperfection, etc. on the reliability of circular concrete-filled steel tube column are studied. It can be con- cluded that inevitable random fluctuation of those parameters has significant influence on structural reliability, and that stochastic or reliability methods can provide a more rational and subjective evaluation on the safety of CFT structures than a deterministic approach.

Geometrically Nonlinear Analysis of Structures Using Various Higher Order Solution Methods: A Comparative Analysis for Large Deformation

The suitability of six higher order root solvers is examined for solving the nonlinear equilibrium equations in large deformation analysis of structures.The applied methods have a better convergence rate than the quadratic Newton-Raphson method.These six methods do not require higher order derivatives to achieve a higher convergence rate.Six algorithms are developed to use the higher order methods in place of the Newton-Raphson method to solve the nonlinear equilibrium equations in geometrically nonlinear analysis of structures.The higher order methods are applied to both continuum and discrete problems(spherical shell and dome truss).The computational cost and the sensitivity of the higher order solution methods and the Newton-Raphson method with respect to the load increment size are comparatively investigated.The numerical results reveal that the higher order methods require a lower number of iterations that the Newton-Raphson method to converge.It is also shown that these methods are less sensitive to the variation of the load increment size.As it is indicated in numerical results,the average residual reduces in a lower number of iterations by the application of the higher order methods in the nonlinear analysis of structures.

Large deflection response-based geometrical nonlinearity of nanocomposite structures reinforced with carbon nanotubes

This paper deals with the nonlinear large deflection analysis of functionally graded carbon nanotube-reinforced composite(FG-CNTRC)plates and panels using a finite element method.Based on the first-order shear deformation theory(FSDT),the proposed model takes into account the transverse shear deformations and incorporates the geometrical nonlinearity type.A C^0 isoparametric finite shell element is developed for the present nonlinear model with the description of large displacements and finite rotations.By adopting the extended rule of mixture,the effective material properties of FG-CNTRCs are approximated with the introduction of some efficiency parameters.Four carbon nanotube(CNT)distributions,labeled uniformly distributed(UD)-CNT,FG-VCNT,FG-O-CNT,and FG-X-CNT,are considered.The solution procedure is carried out via the Newton-Raphson incremental technique.Various numerical applications in both isotropic and CNTRC composite cases are performed to trace the potential of the present model.The effects of the CNT distributions,their volume fractions,and the geometrical characteristics on the nonlinear deflection responses of FG-CNTRC structures are highlighted via a detailed parametric study.

SIX-VARIABLE GEOMETRICAL NONLINEAR LAMINATED THEORY FOR LARGE DEFORMATION

A six-variable geometrical nonlinear shear deformation laminatedtheory is presented by which normal stress and strain distributioncan be calculated. by considering some affective factors that wereneglected under the finite deformation condition, an improved vonKarman geometrical non- linear deformation-strain relation is usedfor large deformation analysis. After analyzing the bending problemof laminated plates, and comparing it with 3-D elasticity solutionsand J.N.Reddy five-vari- able simple higher-order shear deformationlaminated theory, we can conclude that a satisfactory cal- culationprecision has been achieved, which shows that is especially suitablefor the calculation in the condition of large deformation and thelaminated thick plate analysis.

GEOMETRICALLY NONLINEAR ANALYSIS OF MINDLIN PLATEUSING THE INCOMPATIBLE BENDING ELEMENTSWITH INTERNAL SHEAR STRAIN

An approach of the incompatible elements with additional internal shear strain is,in the presem paper,suggested and applied to geometrically nonlinear analysis of Mi-ndlin plate bending problem.It provides a quite covenient way to avoid the whear loc-king troubles.An energy consistency condition for this kind of C°elements is offered.The nonlinear element formulations and some numerical results are presented.

2D Minimum Compliance Topology Optimization Based on a Region Partitioning Strategy

This paper presents an extended sequential element rejection and admission(SERA)topology optimizationmethod with a region partitioning strategy.Based on the partitioning of a design domain into solid regions and weak regions,the proposed optimizationmethod sequentially implements finite element analysis(FEA)in these regions.After standard FEA in the solid regions,the boundary displacement of the weak regions is constrained using the numerical solution of the solid regions as Dirichlet boundary conditions.This treatment can alleviate the negative effect of the material interpolation model of the topology optimization method in the weak regions,such as the condition number of the structural global stiffness matrix.For optimization,in which the forward problem requires nonlinear structural analysis,a linear solver can be applied in weak regions to avoid numerical singularities caused by the over-deformedmesh.To enhance the robustness of the proposedmethod,the nonmanifold point and island are identified and handled separately.The performance of the proposed method is verified by three 2D minimum compliance examples.

Optimal Design of Compliant Trailing Edge for Shape Changing

Adaptive wings have long used smooth morphing technique of compliant leading and trailing edge to improve their aerodynamic characteristics. This paper introduces a systematic approach to design compliant structures to carry out required shape changes under distributed pressure loads. In order to minimize the deviation of the deformed shape from the target shape, this method uses MATLAB and ANSYS to optimize the distributed compliant mechanisms by way of the ground approach and genetic algorithm （GA） to remove the elements possessive of very low stresses. In the optimization process, many factors should be considered such as airloads, input displacements, and geometric nonlinearities. Direct search method is used to locally optimize the dimension and input displacement after the GA optimization. The resultant structure could make its shape change from 0 to 9.3 degrees. The experimental data of the model confirms the feasibility of this approach.