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 functiona...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.展开更多
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 duri...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.展开更多
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...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.展开更多
This paper describes a method proposed for modeling large deflection of aircraft in nonlinear aeroelastic analysis by developing reduced order model(ROM).The method is applied for solving the static aeroelastic and ...This paper describes a method proposed for modeling large deflection of aircraft in nonlinear aeroelastic analysis by developing reduced order model(ROM).The method is applied for solving the static aeroelastic and static aeroelastic trim problems of flexible aircraft containing geometric nonlinearities;meanwhile,the non-planar effects of aerodynamics and follower force effect have been considered.ROMs are computational inexpensive mathematical representations compared to traditional nonlinear finite element method(FEM) especially in aeroelastic solutions.The approach for structure modeling presented here is on the basis of combined modal/finite element(MFE) method that characterizes the stiffness nonlinearities and we apply that structure modeling method as ROM to aeroelastic analysis.Moreover,the non-planar aerodynamic force is computed by the non-planar vortex lattice method(VLM).Structure and aerodynamics can be coupled with the surface spline method.The results show that both of the static aeroelastic analysis and trim analysis of aircraft based on structure ROM can achieve a good agreement compared to analysis based on the FEM and experimental result.展开更多
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...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.展开更多
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...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 c...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.展开更多
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 a...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.展开更多
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 stif...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.展开更多
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 difficul...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.展开更多
Instead of using the previous straight beam element to approximate the curved beam,in this paper,a curvilinear coordinate is employed to describe the deformations,and a new curved beam element is proposed to model the...Instead of using the previous straight beam element to approximate the curved beam,in this paper,a curvilinear coordinate is employed to describe the deformations,and a new curved beam element is proposed to model the curved beam.Based on exact nonlinear strain-displacement relation,virtual work principle is used to derive dynamic equations for a rotating curved beam,with the effects of axial extensibility,shear deformation and rotary inertia taken into account.The constant matrices are solved numerically utilizing the Gauss quadrature integration method.Newmark and Newton-Raphson iteration methods are adopted to solve the differential equations of the rigid-flexible coupling system.The present results are compared with those obtained by commercial programs to validate the present finite method.In order to further illustrate the convergence and efficiency characteristics of the present modeling and computation formulation,comparison of the results of the present formulation with those of the ADAMS software are made.Furthermore,the present results obtained from linear formulation are compared with those from nonlinear formulation,and the special dynamic characteristics of the curved beam are concluded by comparison with those of the straight beam.展开更多
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...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.展开更多
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 so...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.展开更多
In this paper, we present a strong-form framework for solving the boundary value problems with geometric nonlinearity, in which an incremental theory is developed for the problem based on the Newton-Raphson scheme. Co...In this paper, we present a strong-form framework for solving the boundary value problems with geometric nonlinearity, in which an incremental theory is developed for the problem based on the Newton-Raphson scheme. Conventionally, the finite ele- ment methods (FEMs) or weak-form based meshfree methods have often been adopted to solve geometric nonlinear problems. However, issues, such as the mesh dependency, the numerical integration, and the boundary imposition, make these approaches com- putationally inefficient. Recently, strong-form collocation methods have been called on to solve the boundary value problems. The feasibility of the collocation method with the nodal discretization such as the radial basis collocation method (RBCM) motivates the present study. Due to the limited application to the nonlinear analysis in a strong form, we formulate the equation of equilibrium, along with the boundary conditions, in an incremental-iterative sense using the RBCM. The efficacy of the proposed framework is numerically demonstrated with the solution of two benchmark problems involving the geometric nonlinearity. Compared with the conventional weak-form formulation, the pro- posed framework is advantageous as no quadrature rule is needed in constructing the governing equation, and no mesh limitation exists with the deformed geometry in the increment al-it erative process.展开更多
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 acousti...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.展开更多
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 inter...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.展开更多
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 c...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.展开更多
In order to analyze bellows effectively and practically, the finite_element_displacement_perturbation method (FEDPM) is proposed for the geometric nonlinear behaviors of shells of revolution subjected to pure bending ...In order to analyze bellows effectively and practically, the finite_element_displacement_perturbation method (FEDPM) is proposed for the geometric nonlinear behaviors of shells of revolution subjected to pure bending moments or lateral forces in one of their meridional planes. The formulations are mainly based upon the idea of perturba_ tion that the nodal displacement vector and the nodal force vector of each finite element are expanded by taking root_mean_square value of circumferential strains of the shells as a perturbation parameter. The load steps and the iteration times are not as arbitrary and unpredictable as in usual nonlinear analysis. Instead, there are certain relations between the load steps and the displacement increments, and no need of iteration for each load step. Besides, in the formulations, the shell is idealized into a series of conical frusta for the convenience of practice, Sander's nonlinear geometric equations of moderate small rotation are used, and the shell made of more than one material ply is also considered.展开更多
Based on the actual measured well depth, azimuth and oblique angles, a novel interpolation method to obtain the well axis is developed. The initial stress of drill string at the reference state consistent with well ax...Based on the actual measured well depth, azimuth and oblique angles, a novel interpolation method to obtain the well axis is developed. The initial stress of drill string at the reference state consistent with well axis can be obtained from the curvature and the tortuosity of well axis. By using the principle of virtual work, the formula to compute the equivalent load vector of the initial stress was derived. In the derivation,the natural (curvilinear) coordinate system was adopted since both the curvature and the tortuosity were generally not zero. A set of displacement functions fully reflecting the rigid body modes was used. Some basic concepts in the finite element analysis of drill string were clarified. It is hoped that the proposed method would offer a theoretical basis for handling the geometric nonlinear problem of the drill string in a 3-D larg edisplacement wellbore.展开更多
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 quad...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.展开更多
基金the National Natural Science Foundation of China(NSFC)under Grant Nos.12272124 and 11972146.
文摘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.
基金the authority of the National Natural Science Foundation of China(Grant Nos.52178168 and 51378427)for financing this research work and several ongoing research projects related to structural impact performance.
文摘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.
基金supported by National Science Foundation for Distinguished Young Scholars of China(Grant No.50825504)National Natural Science Foundation of China(Grant No.50775073)United Fund of Natural Science Foundation of China and Guangdong Province (Grant No.U0934004)
文摘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.
文摘This paper describes a method proposed for modeling large deflection of aircraft in nonlinear aeroelastic analysis by developing reduced order model(ROM).The method is applied for solving the static aeroelastic and static aeroelastic trim problems of flexible aircraft containing geometric nonlinearities;meanwhile,the non-planar effects of aerodynamics and follower force effect have been considered.ROMs are computational inexpensive mathematical representations compared to traditional nonlinear finite element method(FEM) especially in aeroelastic solutions.The approach for structure modeling presented here is on the basis of combined modal/finite element(MFE) method that characterizes the stiffness nonlinearities and we apply that structure modeling method as ROM to aeroelastic analysis.Moreover,the non-planar aerodynamic force is computed by the non-planar vortex lattice method(VLM).Structure and aerodynamics can be coupled with the surface spline method.The results show that both of the static aeroelastic analysis and trim analysis of aircraft based on structure ROM can achieve a good agreement compared to analysis based on the FEM and experimental result.
文摘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.
文摘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.
基金Supported by the National Natural Science Foundation of China Youth Fund under Grant No.51909148the Doctoral Research Initiation Fund of Shandong Jiaotong University under Grant No.BS2018001。
文摘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.
基金supported by the National Science Fund for Distinguished Young Scholars (No. 50725826).
文摘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.
基金Project(12 High-tech Urban C11) supported by High-tech Urban Development Program of Ministry of Land,Transport and Maritime Affairs,Korea
文摘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.
基金Project supported by the National 973 Program (No.2004CB719402), the National Natural Science Foundation of China (No. 10372030)the Open Research Projects supported by the Project Fund of the Hubei Province Key Lab of Mechanical Transmission & Manufacturing Engineering Wuhan University of Science & Technology (No.2003A16).
文摘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.
基金supported by the National Natural Science Foundation of China(10872126)Research Fund for the Doctoral Program of Higher Education of China(20100073110007)
文摘Instead of using the previous straight beam element to approximate the curved beam,in this paper,a curvilinear coordinate is employed to describe the deformations,and a new curved beam element is proposed to model the curved beam.Based on exact nonlinear strain-displacement relation,virtual work principle is used to derive dynamic equations for a rotating curved beam,with the effects of axial extensibility,shear deformation and rotary inertia taken into account.The constant matrices are solved numerically utilizing the Gauss quadrature integration method.Newmark and Newton-Raphson iteration methods are adopted to solve the differential equations of the rigid-flexible coupling system.The present results are compared with those obtained by commercial programs to validate the present finite method.In order to further illustrate the convergence and efficiency characteristics of the present modeling and computation formulation,comparison of the results of the present formulation with those of the ADAMS software are made.Furthermore,the present results obtained from linear formulation are compared with those from nonlinear formulation,and the special dynamic characteristics of the curved beam are concluded by comparison with those of the straight beam.
文摘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.
基金the National '973' Key Fundamental Research Projects of China(No.2003CB716207)the National '863' High-Tech Development Projects of China(No.2006AA04Z162)also the Australian Research Council(No.ARC-DP0666683).
文摘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.
基金Project supported by the Ministry of Science and Technology of Taiwan(No.MOST 104-2221-E-009-193)
文摘In this paper, we present a strong-form framework for solving the boundary value problems with geometric nonlinearity, in which an incremental theory is developed for the problem based on the Newton-Raphson scheme. Conventionally, the finite ele- ment methods (FEMs) or weak-form based meshfree methods have often been adopted to solve geometric nonlinear problems. However, issues, such as the mesh dependency, the numerical integration, and the boundary imposition, make these approaches com- putationally inefficient. Recently, strong-form collocation methods have been called on to solve the boundary value problems. The feasibility of the collocation method with the nodal discretization such as the radial basis collocation method (RBCM) motivates the present study. Due to the limited application to the nonlinear analysis in a strong form, we formulate the equation of equilibrium, along with the boundary conditions, in an incremental-iterative sense using the RBCM. The efficacy of the proposed framework is numerically demonstrated with the solution of two benchmark problems involving the geometric nonlinearity. Compared with the conventional weak-form formulation, the pro- posed framework is advantageous as no quadrature rule is needed in constructing the governing equation, and no mesh limitation exists with the deformed geometry in the increment al-it erative process.
基金supported by the National Natural Science Foundations of China(Nos.11872079,11572109)the Science and Technology Project of Hebei Education Department(No.QN2019135)Advanced Talents Incubation Program of the Hebei University(No.521000981285)。
文摘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.
基金supported by the National Natural Science Foundation of China (50725826)Specific Research on Cable-reinforced Membranes with Super Span and Complex Single-shell Structures of Expo Axis (08dz0580303)Shanghai Postdoctoral Fund (10R21416200)
文摘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.
基金supported by the Fundamental Research Funds for the Central Universities (SWJTU09CX012 and SWJTU11BR006)the Doctoral Fund for Youth Scholars of Ministry of Educationof China (No. 20110184120010)
文摘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.
文摘In order to analyze bellows effectively and practically, the finite_element_displacement_perturbation method (FEDPM) is proposed for the geometric nonlinear behaviors of shells of revolution subjected to pure bending moments or lateral forces in one of their meridional planes. The formulations are mainly based upon the idea of perturba_ tion that the nodal displacement vector and the nodal force vector of each finite element are expanded by taking root_mean_square value of circumferential strains of the shells as a perturbation parameter. The load steps and the iteration times are not as arbitrary and unpredictable as in usual nonlinear analysis. Instead, there are certain relations between the load steps and the displacement increments, and no need of iteration for each load step. Besides, in the formulations, the shell is idealized into a series of conical frusta for the convenience of practice, Sander's nonlinear geometric equations of moderate small rotation are used, and the shell made of more than one material ply is also considered.
文摘Based on the actual measured well depth, azimuth and oblique angles, a novel interpolation method to obtain the well axis is developed. The initial stress of drill string at the reference state consistent with well axis can be obtained from the curvature and the tortuosity of well axis. By using the principle of virtual work, the formula to compute the equivalent load vector of the initial stress was derived. In the derivation,the natural (curvilinear) coordinate system was adopted since both the curvature and the tortuosity were generally not zero. A set of displacement functions fully reflecting the rigid body modes was used. Some basic concepts in the finite element analysis of drill string were clarified. It is hoped that the proposed method would offer a theoretical basis for handling the geometric nonlinear problem of the drill string in a 3-D larg edisplacement wellbore.
文摘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.