An approach based on the superelement theory is proposed, and it is applied to model the car's body-in-white as well as to dynamic simulation and optimization. This approach can improve the calculation speed and d...An approach based on the superelement theory is proposed, and it is applied to model the car's body-in-white as well as to dynamic simulation and optimization. This approach can improve the calculation speed and do the dynamic optimization among substructures respectively in the car's body design. To meet the car's design of harshness, a dynamic optimal design model, based on the mean square of vertical displacement response at two points of the car floor, is proposed. Satisfactory results are achieved.展开更多
To estimate the geometric characteristics, especially wet areas and section areas, of three dimensional numerical conceptual aircraft models, a method based on surface elements is proposed. On the premise that numerou...To estimate the geometric characteristics, especially wet areas and section areas, of three dimensional numerical conceptual aircraft models, a method based on surface elements is proposed. On the premise that numerous surface elements are generated to represent each component surface, a component wet area of the surface is estimated by adding up the areas of such elements that are not covered by any other component surfaces. The elements are also used to get the section polygons of such composite surfaces as the whole aircraft at a given body station, then a section area is approximated with the sum of trapezoidal areas between such sides of polygons that are not covered by any other component and a reference axis. Practical application to a computer aided aircraft conceptual design system shows that the methed is applicable to different kinds of conceptual aircraft models and its precision is satisfying to the conceptual design.展开更多
To minimize the reactive power of the converter of the control winding in the novel dual stator-winding induction generator based on the PWM converter, design features of the induction generator with a rectified load ...To minimize the reactive power of the converter of the control winding in the novel dual stator-winding induction generator based on the PWM converter, design features of the induction generator with a rectified load are proposed. The optimization method of excited capacitors to minimize the reactive power of the control winding at a variable speed is given. The calculation capacity of the machine with a diode bridge rectifier load is proposed. To achieve global searching, the integrated method with the improved real-coded genetic algorithm and the twodimensional finite element method (FEM) is introduced. Design results of the sample show that reactive power can be reduced by the method, and the converter capacity can be decreased to 1/3 of output rated power at the speed ratio of 1 : 3, thus reducing the volume and the mass of the inverter.展开更多
Flat metal strips are deformed progressively into profiles with certain geometry sections by a series of successive rotating rolls at the room temperature, and it is called the cold roll-forming. An effective method i...Flat metal strips are deformed progressively into profiles with certain geometry sections by a series of successive rotating rolls at the room temperature, and it is called the cold roll-forming. An effective method is given for roll design to avoid major defects via finite element method (FEM) simulation using a booting model. The simulation gives reasonable fit to the actual product in the two major roll-forming defects, i.e. , edge wave and springback. The redesigned rolls using the approach of multi-stand FEM simulation can effectively control these two defects.展开更多
An integrated system FSOP2D,including modules for the shape optimal modeling,structural analysis,sensitivity analysis,optimal method library and post- processing,is developed.By selecting fictitious loads as the desig...An integrated system FSOP2D,including modules for the shape optimal modeling,structural analysis,sensitivity analysis,optimal method library and post- processing,is developed.By selecting fictitious loads as the design variables that has a linear relationship with the grid point locations and using design sensitivity analysis of the domain method,it is easier to solve the velocity field.In the course of optimal iterations,mesh distortion is kept to a minimum,sensitivity derivatives of object function,stress constraints and displacement constraints are derived.Computation of sensitivity analysis is achieved in the system.Two engineering examples are used to prove the system's effectiveness,the optimal results can successfully be obtained by lesser number of iterations.展开更多
This paper discusses some aspects of finite element computation,such as the automatic generation of finite element ,refinement of mesh,process of node density, distribution of load,optimum design and the drawing o...This paper discusses some aspects of finite element computation,such as the automatic generation of finite element ,refinement of mesh,process of node density, distribution of load,optimum design and the drawing of stress contour, and describes the developing process of software for a planar 8 node element.展开更多
The effect of flow passage length in the die cavity and extrusion wheel velocity on the shape of aluminum sheath during the continuous extrusion sheathing process was analyzed by using finite element methods based on ...The effect of flow passage length in the die cavity and extrusion wheel velocity on the shape of aluminum sheath during the continuous extrusion sheathing process was analyzed by using finite element methods based on software DEFORM 3D and experimentally validated. The results show that by increasing the flow passage length, the velocity of metal at the cross-section of sheath tends toward uniformity, the values of the bending angles of sheath gradually approach the ideal value of zero and the cross-section exhibits a better shape. The extrusion wheel velocity has negligible effects on the bending shape and cross-section of the sheath product when a long flow passage is used.展开更多
Based on the engineering background of the Jiangxinzhou Bridge in Nanjing, issues related to the spatial main saddle of the self-anchored suspension bridge are studied. The refinement finite element model is establish...Based on the engineering background of the Jiangxinzhou Bridge in Nanjing, issues related to the spatial main saddle of the self-anchored suspension bridge are studied. The refinement finite element model is established by the secondary development technology based on the platform of the general finite element program, and a reasonable load pattern is used in its spatial structural analysis, by which its path of force transference and stress distribution are obtained. Matched with the spatial main cable, the tangency point correction method is also discussed. The results show that the lateral wall stress of the saddle groove is higher than the stress within the wall due to the role of lateral forces in the finished bridge state; the horizontal volume force of the main cable can generate a gradient distributed vertical extrusion pressure on the saddle clamping device and the main saddle body; the geometric nonlinear effect of the self- anchored suspension bridge cable system in the construction process is significant, which can be reflected in the spatial tangent point position of the main cable with the main saddle changes a lot from free cable to finished cable.展开更多
From a perspective of theoretical study, there are some faults in the models of the existing object-oriented programming languages. For example, C# does not support metaclasses, the primitive types of Java and C# are ...From a perspective of theoretical study, there are some faults in the models of the existing object-oriented programming languages. For example, C# does not support metaclasses, the primitive types of Java and C# are not objects, etc. So, this paper designs a programming language, Shrek, which integrates many language features and constructions in a compact and consistent model. The Shrek language is a class-based purely object-oriented language. It has a dynamical strong type system, and adopts a single-inheritance mechanism with Mixin as its complement. It has a consistent class instantiation and inheritance structure, and the ability of intercessive structural computational reflection, which enables it to support safe metaclass programming. It also supports multi-thread programming and automatic garbage collection, and enforces its expressive power by adopting a native method mechanism. The prototype system of the Shrek language is implemented and anticipated design goals are achieved.展开更多
In this work, power efficient butterfly unit based FFT architecture is presented. The butterfly unit is designed using floating-point fused arithmetic units. The fused arithmetic units include two-term dot product uni...In this work, power efficient butterfly unit based FFT architecture is presented. The butterfly unit is designed using floating-point fused arithmetic units. The fused arithmetic units include two-term dot product unit and add-subtract unit. In these arithmetic units, operations are performed over complex data values. A modified fused floating-point two-term dot product and an enhanced model for the Radix-4 FFT butterfly unit are proposed. The modified fused two-term dot product is designed using Radix-16 booth multiplier. Radix-16 booth multiplier will reduce the switching activities compared to Radix-8 booth multiplier in existing system and also will reduce the area required. The proposed architecture is implemented efficiently for Radix-4 decimation in time(DIT) FFT butterfly with the two floating-point fused arithmetic units. The proposed enhanced architecture is synthesized, implemented, placed and routed on a FPGA device using Xilinx ISE tool. It is observed that the Radix-4 DIT fused floating-point FFT butterfly requires 50.17% less space and 12.16% reduced power compared to the existing methods and the proposed enhanced model requires 49.82% less space on the FPGA device compared to the proposed design. Also, reduced power consumption is addressed by utilizing the reusability technique, which results in 11.42% of power reduction of the enhanced model compared to the proposed design.展开更多
The hydro-hammer sampler is a new type of sampler compared with traditional ones. An important part of this new offshore sampler is that the structure of the core cutter has a significant effect on penetration and cor...The hydro-hammer sampler is a new type of sampler compared with traditional ones. An important part of this new offshore sampler is that the structure of the core cutter has a significant effect on penetration and core recovery. In our experiments, a commercial finite element code with a capability of simulating large-strain frictional contact between two or more solid bodies is used to simulate the core cutter-soil interaction. The effects of the cutting edge shape, the diameter and the edge angle on penetration are analyzed by non-liner transient dynamic analysis using a finite element method (FEM). Simulation results show that the cutter shape clearly has an effect on the penetration and core recovery. In addition, the penetration of the sampler increases with an increase in the inside diameter of the cutter, but decreases with an increase in the cutting angle. Based on these analyses, an optimum structure of the core cutter is designed and tested in the north margin of the Dalian gulf. Experiment results show that the penetration rate is about 16.5 m/h in silty clay and 15.4 m/h in cohesive clay, while the recovery is 68% and 83.3% resoectively.展开更多
For the tunnel crossing active fault,the damage induced by fault movement is always serious.To solve such a problem,a detailed anti-faulting tunnel design process for Urumqi subway line 2 was introduced,and seven thre...For the tunnel crossing active fault,the damage induced by fault movement is always serious.To solve such a problem,a detailed anti-faulting tunnel design process for Urumqi subway line 2 was introduced,and seven three-dimensional elastic-plastic finite element models were established.The anti-faulting design process included three steps.First,the damage of tunnel lining from different locations of fault rupture surfaces was analyzed.Then,the analysis of the effect on tunnel buried depth was given.Finally,the effect of the disaster mitigation method on the flexible joint was verified and the location of the flexible joint was discussed.The results show that when the properties of surrounding rock at the tunnel bottom grows soft,the tunnel deformation curve is smoother and tunnel damage induced by fault movement is less serious.The vertical displacement change ratio of secondary linings along the tunnel axis may be the main factor to cause shear damage to the tunnel.The interface between the hanging wall and fracture zone is defined as the most adverse fault rupture surface.The tunnel damage was reduced with the decrease in the tunnel buried depth as more energy was dissipated by overburden soil and the differential uplift zone of soil became more diffuse.The method of the flexible joint can reduce the tunnel damage significantly and the disaster mitigation effect of different locations on the flexible joint is different.The tunnel damage is reduced by the greatest degree when the flexible joint is located on the fault rupture surface.展开更多
A microelectromechanical Digital to Analog Converter (DAC) based on Weighted Gap (WG) principle is described,which is analogous to the weighed resistor DAC in electronic circuits.To convert the input of binary vol...A microelectromechanical Digital to Analog Converter (DAC) based on Weighted Gap (WG) principle is described,which is analogous to the weighed resistor DAC in electronic circuits.To convert the input of binary voltage to the output of analog displacement,the gaps are proposed to be employed as a scale factor.A finite element method is used to simulate the performance of the DAC.To reduce the error,the structure design is optimized and the maximum error of 0 002μm is obtained.展开更多
The main aim of this paper is to have an accurate analysis on the famous Adini's element for the second order problems under to the anisotropic meshes. We firstly show that the interpolation of Adini's element satis...The main aim of this paper is to have an accurate analysis on the famous Adini's element for the second order problems under to the anisotropic meshes. We firstly show that the interpolation of Adini's element satisfy the anisotropic property. Then the optimal error estimate is obtained without the regularity assumption on the meshes.展开更多
A new set of binary sequences-Periodic Complementary Binary Sequence Pair (PCSP) is proposed. A new class of block design-Difference Family Pair (DFP) is also proposed.The relationship between PCSP and DFP, the proper...A new set of binary sequences-Periodic Complementary Binary Sequence Pair (PCSP) is proposed. A new class of block design-Difference Family Pair (DFP) is also proposed.The relationship between PCSP and DFP, the properties and existing conditions of PCSP and the recursive constructions for PCSP are given.展开更多
A 32 kbit OTP(one-time programmable)memory for MCUs(micro-controller units)used in remote controllers was designed.This OTP memory is used for program and data storage.It is required to apply 5.5V to BL(bit-line)and 1...A 32 kbit OTP(one-time programmable)memory for MCUs(micro-controller units)used in remote controllers was designed.This OTP memory is used for program and data storage.It is required to apply 5.5V to BL(bit-line)and 11V to WL(word-line)for a OTP cell of 0.35μm ETOX(EEPROM tunnel oxide)type by MagnaChip.We use 5V transistors on column data paths to reduce the area of column data paths since they require small areas.In addition,we secure device reliability by using HV(high-voltage)transistors in the WL driver.Furthermore,we change from a static logic to a dynamic logic used for the WL driver in the core circuit.Also,we optimize the WD(write data)switch circuit.Thus,we can implement them with a small-area design.In addition,we implement the address predecoder with a small-area logic circuit.The area of the designed 32 kbit OTP with 5V and HV devices is 674.725μm×258.75μm(=0.1745mm2)and is 56.3% smaller than that using 3.3V devices.展开更多
This paper proposes a theoretical method using finite element analysis(FEA) to calculate the plastic collapse loads of pressure vessels under internal pressure,and compares the analytical methods according to three cr...This paper proposes a theoretical method using finite element analysis(FEA) to calculate the plastic collapse loads of pressure vessels under internal pressure,and compares the analytical methods according to three criteria stated in the ASME Boiler Pressure Vessel Code. First,a finite element technique using the arc-length algorithm and the restart analysis is developed to conduct the plastic collapse analysis of vessels,which includes the material and geometry non-linear properties of vessels. Second,as the mechanical properties of vessels are assumed to be elastic-perfectly plastic,the limit load analysis is performed by em-ploying the Newton-Raphson algorithm,while the limit pressure of vessels is obtained by the twice-elastic-slope method and the tangent intersection method respectively to avoid excessive deformation. Finally,the elastic stress analysis under working pressure is conducted and the stress strength of vessels is checked by sorting the stress results. The results are compared with those obtained by experiments and other existing models. This work provides a reference for the selection of the failure criteria and the calculation of the plastic collapse load.展开更多
文摘An approach based on the superelement theory is proposed, and it is applied to model the car's body-in-white as well as to dynamic simulation and optimization. This approach can improve the calculation speed and do the dynamic optimization among substructures respectively in the car's body design. To meet the car's design of harshness, a dynamic optimal design model, based on the mean square of vertical displacement response at two points of the car floor, is proposed. Satisfactory results are achieved.
文摘To estimate the geometric characteristics, especially wet areas and section areas, of three dimensional numerical conceptual aircraft models, a method based on surface elements is proposed. On the premise that numerous surface elements are generated to represent each component surface, a component wet area of the surface is estimated by adding up the areas of such elements that are not covered by any other component surfaces. The elements are also used to get the section polygons of such composite surfaces as the whole aircraft at a given body station, then a section area is approximated with the sum of trapezoidal areas between such sides of polygons that are not covered by any other component and a reference axis. Practical application to a computer aided aircraft conceptual design system shows that the methed is applicable to different kinds of conceptual aircraft models and its precision is satisfying to the conceptual design.
文摘To minimize the reactive power of the converter of the control winding in the novel dual stator-winding induction generator based on the PWM converter, design features of the induction generator with a rectified load are proposed. The optimization method of excited capacitors to minimize the reactive power of the control winding at a variable speed is given. The calculation capacity of the machine with a diode bridge rectifier load is proposed. To achieve global searching, the integrated method with the improved real-coded genetic algorithm and the twodimensional finite element method (FEM) is introduced. Design results of the sample show that reactive power can be reduced by the method, and the converter capacity can be decreased to 1/3 of output rated power at the speed ratio of 1 : 3, thus reducing the volume and the mass of the inverter.
文摘Flat metal strips are deformed progressively into profiles with certain geometry sections by a series of successive rotating rolls at the room temperature, and it is called the cold roll-forming. An effective method is given for roll design to avoid major defects via finite element method (FEM) simulation using a booting model. The simulation gives reasonable fit to the actual product in the two major roll-forming defects, i.e. , edge wave and springback. The redesigned rolls using the approach of multi-stand FEM simulation can effectively control these two defects.
文摘An integrated system FSOP2D,including modules for the shape optimal modeling,structural analysis,sensitivity analysis,optimal method library and post- processing,is developed.By selecting fictitious loads as the design variables that has a linear relationship with the grid point locations and using design sensitivity analysis of the domain method,it is easier to solve the velocity field.In the course of optimal iterations,mesh distortion is kept to a minimum,sensitivity derivatives of object function,stress constraints and displacement constraints are derived.Computation of sensitivity analysis is achieved in the system.Two engineering examples are used to prove the system's effectiveness,the optimal results can successfully be obtained by lesser number of iterations.
文摘This paper discusses some aspects of finite element computation,such as the automatic generation of finite element ,refinement of mesh,process of node density, distribution of load,optimum design and the drawing of stress contour, and describes the developing process of software for a planar 8 node element.
基金Project (51175055) supported by the National Natural Science Foundation of ChinaProject (201102020) supported by the Natural Science Foundation of Liaoning Province, ChinaProject (200921085) supported by the Liaoning BaiQian Wan Talents Program, China
文摘The effect of flow passage length in the die cavity and extrusion wheel velocity on the shape of aluminum sheath during the continuous extrusion sheathing process was analyzed by using finite element methods based on software DEFORM 3D and experimentally validated. The results show that by increasing the flow passage length, the velocity of metal at the cross-section of sheath tends toward uniformity, the values of the bending angles of sheath gradually approach the ideal value of zero and the cross-section exhibits a better shape. The extrusion wheel velocity has negligible effects on the bending shape and cross-section of the sheath product when a long flow passage is used.
基金The National High Technology Research and Development Program of China(863 Program)(No.2006AA04Z416)the National Science Fund for Distinguished Young Scholars(No.50725828)
文摘Based on the engineering background of the Jiangxinzhou Bridge in Nanjing, issues related to the spatial main saddle of the self-anchored suspension bridge are studied. The refinement finite element model is established by the secondary development technology based on the platform of the general finite element program, and a reasonable load pattern is used in its spatial structural analysis, by which its path of force transference and stress distribution are obtained. Matched with the spatial main cable, the tangency point correction method is also discussed. The results show that the lateral wall stress of the saddle groove is higher than the stress within the wall due to the role of lateral forces in the finished bridge state; the horizontal volume force of the main cable can generate a gradient distributed vertical extrusion pressure on the saddle clamping device and the main saddle body; the geometric nonlinear effect of the self- anchored suspension bridge cable system in the construction process is significant, which can be reflected in the spatial tangent point position of the main cable with the main saddle changes a lot from free cable to finished cable.
基金The National Science Fund for Distinguished Young Scholars (No.60425206)the National Natural Science Foundation of China (No.60633010)the Natural Science Foundation of Jiangsu Province(No.BK2006094)
文摘From a perspective of theoretical study, there are some faults in the models of the existing object-oriented programming languages. For example, C# does not support metaclasses, the primitive types of Java and C# are not objects, etc. So, this paper designs a programming language, Shrek, which integrates many language features and constructions in a compact and consistent model. The Shrek language is a class-based purely object-oriented language. It has a dynamical strong type system, and adopts a single-inheritance mechanism with Mixin as its complement. It has a consistent class instantiation and inheritance structure, and the ability of intercessive structural computational reflection, which enables it to support safe metaclass programming. It also supports multi-thread programming and automatic garbage collection, and enforces its expressive power by adopting a native method mechanism. The prototype system of the Shrek language is implemented and anticipated design goals are achieved.
文摘In this work, power efficient butterfly unit based FFT architecture is presented. The butterfly unit is designed using floating-point fused arithmetic units. The fused arithmetic units include two-term dot product unit and add-subtract unit. In these arithmetic units, operations are performed over complex data values. A modified fused floating-point two-term dot product and an enhanced model for the Radix-4 FFT butterfly unit are proposed. The modified fused two-term dot product is designed using Radix-16 booth multiplier. Radix-16 booth multiplier will reduce the switching activities compared to Radix-8 booth multiplier in existing system and also will reduce the area required. The proposed architecture is implemented efficiently for Radix-4 decimation in time(DIT) FFT butterfly with the two floating-point fused arithmetic units. The proposed enhanced architecture is synthesized, implemented, placed and routed on a FPGA device using Xilinx ISE tool. It is observed that the Radix-4 DIT fused floating-point FFT butterfly requires 50.17% less space and 12.16% reduced power compared to the existing methods and the proposed enhanced model requires 49.82% less space on the FPGA device compared to the proposed design. Also, reduced power consumption is addressed by utilizing the reusability technique, which results in 11.42% of power reduction of the enhanced model compared to the proposed design.
基金Project 20002070005126 supported by the China Geological Survey
文摘The hydro-hammer sampler is a new type of sampler compared with traditional ones. An important part of this new offshore sampler is that the structure of the core cutter has a significant effect on penetration and core recovery. In our experiments, a commercial finite element code with a capability of simulating large-strain frictional contact between two or more solid bodies is used to simulate the core cutter-soil interaction. The effects of the cutting edge shape, the diameter and the edge angle on penetration are analyzed by non-liner transient dynamic analysis using a finite element method (FEM). Simulation results show that the cutter shape clearly has an effect on the penetration and core recovery. In addition, the penetration of the sampler increases with an increase in the inside diameter of the cutter, but decreases with an increase in the cutting angle. Based on these analyses, an optimum structure of the core cutter is designed and tested in the north margin of the Dalian gulf. Experiment results show that the penetration rate is about 16.5 m/h in silty clay and 15.4 m/h in cohesive clay, while the recovery is 68% and 83.3% resoectively.
基金The National Natural Science Foundation of China(No.41572276)the National Key Research and Development Program of China(No.2017YFC0805400).
文摘For the tunnel crossing active fault,the damage induced by fault movement is always serious.To solve such a problem,a detailed anti-faulting tunnel design process for Urumqi subway line 2 was introduced,and seven three-dimensional elastic-plastic finite element models were established.The anti-faulting design process included three steps.First,the damage of tunnel lining from different locations of fault rupture surfaces was analyzed.Then,the analysis of the effect on tunnel buried depth was given.Finally,the effect of the disaster mitigation method on the flexible joint was verified and the location of the flexible joint was discussed.The results show that when the properties of surrounding rock at the tunnel bottom grows soft,the tunnel deformation curve is smoother and tunnel damage induced by fault movement is less serious.The vertical displacement change ratio of secondary linings along the tunnel axis may be the main factor to cause shear damage to the tunnel.The interface between the hanging wall and fracture zone is defined as the most adverse fault rupture surface.The tunnel damage was reduced with the decrease in the tunnel buried depth as more energy was dissipated by overburden soil and the differential uplift zone of soil became more diffuse.The method of the flexible joint can reduce the tunnel damage significantly and the disaster mitigation effect of different locations on the flexible joint is different.The tunnel damage is reduced by the greatest degree when the flexible joint is located on the fault rupture surface.
文摘A microelectromechanical Digital to Analog Converter (DAC) based on Weighted Gap (WG) principle is described,which is analogous to the weighed resistor DAC in electronic circuits.To convert the input of binary voltage to the output of analog displacement,the gaps are proposed to be employed as a scale factor.A finite element method is used to simulate the performance of the DAC.To reduce the error,the structure design is optimized and the maximum error of 0 002μm is obtained.
基金the Henan Natural Science Foundation(072300410320)the Henan Education Department Foundational Study Foundation(200510460311)
文摘The main aim of this paper is to have an accurate analysis on the famous Adini's element for the second order problems under to the anisotropic meshes. We firstly show that the interpolation of Adini's element satisfy the anisotropic property. Then the optimal error estimate is obtained without the regularity assumption on the meshes.
基金Supported by National Natural Science Foundation of China (69972042),Natural Science Fund of Hebei Provice(599245)and Science Foundation of Yanshan University
文摘A new set of binary sequences-Periodic Complementary Binary Sequence Pair (PCSP) is proposed. A new class of block design-Difference Family Pair (DFP) is also proposed.The relationship between PCSP and DFP, the properties and existing conditions of PCSP and the recursive constructions for PCSP are given.
基金Project supported by the Second Stage of Brain Korea 21 Projects,Korea
文摘A 32 kbit OTP(one-time programmable)memory for MCUs(micro-controller units)used in remote controllers was designed.This OTP memory is used for program and data storage.It is required to apply 5.5V to BL(bit-line)and 11V to WL(word-line)for a OTP cell of 0.35μm ETOX(EEPROM tunnel oxide)type by MagnaChip.We use 5V transistors on column data paths to reduce the area of column data paths since they require small areas.In addition,we secure device reliability by using HV(high-voltage)transistors in the WL driver.Furthermore,we change from a static logic to a dynamic logic used for the WL driver in the core circuit.Also,we optimize the WD(write data)switch circuit.Thus,we can implement them with a small-area design.In addition,we implement the address predecoder with a small-area logic circuit.The area of the designed 32 kbit OTP with 5V and HV devices is 674.725μm×258.75μm(=0.1745mm2)and is 56.3% smaller than that using 3.3V devices.
基金Project (Nos. 2006BAK04A02-02 and 2006BAK02B02-08) supported by the National Key Technology R&D Program, China
文摘This paper proposes a theoretical method using finite element analysis(FEA) to calculate the plastic collapse loads of pressure vessels under internal pressure,and compares the analytical methods according to three criteria stated in the ASME Boiler Pressure Vessel Code. First,a finite element technique using the arc-length algorithm and the restart analysis is developed to conduct the plastic collapse analysis of vessels,which includes the material and geometry non-linear properties of vessels. Second,as the mechanical properties of vessels are assumed to be elastic-perfectly plastic,the limit load analysis is performed by em-ploying the Newton-Raphson algorithm,while the limit pressure of vessels is obtained by the twice-elastic-slope method and the tangent intersection method respectively to avoid excessive deformation. Finally,the elastic stress analysis under working pressure is conducted and the stress strength of vessels is checked by sorting the stress results. The results are compared with those obtained by experiments and other existing models. This work provides a reference for the selection of the failure criteria and the calculation of the plastic collapse load.