An excel-based abaqus parametric modeling and analysis method for corrugated steel web composite channel beam

Abstract

The application discloses an Abaqus parameterized modeling and analysis method for a corrugated steel web combined channel girder based on Excel, and comprises the following steps: (1) a user determines the TGCW size, material and finite element analysis type required by the user according to project data; (2) a program running environment is configured, a third-party library function xlrd2 is installed, and the parameter information obtained in the step (1) is filled into an Excel input table of the corrugated steel web combined channel girder finite element analysis parameters; (3) TGCW parameter information in an Excel file is read through a Python script file; (4) TGCW related model information is defined through the Python script; (5) the Python script is run on an Abaqus platform, a TGCW finite element model is generated, and simulation analysis corresponding to the input parameters is carried out. The application can complete the modeling work of a complex model without relying on third-party modeling software, can realize an automatic modeling process, and is convenient for realizing parameterized analysis on the corrugated steel web combined channel girder.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of three-dimensional finite element modeling and numerical simulation of bridges, and in particular to an Abaqus parametric modeling and analysis method for a corrugated steel web composite channel girder based on Excel. BACKGROUND

[0002] At present, China is in a high-speed development stage of urban transportation and intercity high-speed railway projects, and some original structural forms cannot well adapt to the higher requirements of new projects on structural bearing capacity, stiffness and durability. Combining the corrugated steel web structure and the prestressed channel girder structure to form a new structural system, the corrugated steel web prestressed composite channel girder, will well inherit the advantages of each structure and be suitable for application in rail transit bridges, especially urban rail transit bridges and subway projects.

[0003] Nowadays, the analysis of bridge structures cannot be separated from finite element simulation technology. Abaqus is a powerful finite element engineering simulation software with a large element library and material library, which can analyze complex solid mechanics and structural mechanics systems and handle highly nonlinear problems. Its powerful analysis capability and reliability make it widely used in industrial research in various countries. Python is a powerful object-oriented programming language that allows scripting and rapid development on multiple platforms. The scripting interface provided by the secondary development environment of the Abaqus software (hereinafter referred to as the "Abaqus scripting interface") is a customized development based on the Python language. Users can execute commands through Python script files to control the pre-processing and post-processing processes of Abaqus.

[0004] Generally, for complex three-dimensional models, when users use Abaqus to perform pre-processing operations, they not only need to specify the size information of each part of the structure, but also need to disassemble and assemble it to meet the modeling needs of complex structures. This process requires a lot of time and high model reconstruction cost, and cannot meet the user's demand for parameter analysis of complex structures. In this case, it is particularly important to provide an efficient and convenient three-dimensional finite element parametric modeling method for engineering designers to help them quickly understand and master the performance of this structure. SUMMARY

[0005] The technical problem to be solved by the present application is to provide an Abaqus parametric modeling and analysis method for a corrugated steel web composite channel girder based on Excel, which can complete the modeling work of complex models without relying on third-party modeling software, realize an automatic modeling process, and facilitate parameterized analysis of the corrugated steel web composite channel girder.

[0006] To solve the above technical problems, the present application provides a kind of Abaqus parameterized modeling and analysis method based on Excel corrugated steel web combined channel beam, comprising the following steps:

[0007] (1) user determines TGCW size, material by project data, the type of finite element analysis required by user;Finite element analysis type includes frequency mode analysis, static load analysis, dynamic load analysis parameter information;

[0008] (2) program running environment configuration, install third-party library function xlrd2, fill in the parameter information obtained in step (1) into the finite element analysis parameter Excel input table of corrugated steel web combined channel beam;

[0009] (3) TGCW parameter information in Excel file is read by Python script file;

[0010] (4) TGCW related model information is defined by Python script, model information includes Part, Property, Assembly, Mesh, Step, Interaction, Load, Job;

[0011] (5) Abaqus platform is used to run Python script, and TGCW finite element model is generated and simulation analysis corresponding to input parameter is carried out.

[0012] Preferably, in step (2), after xlrd2 installation is completed, user needs to copy the corresponding library function file under site-packages folder in system Python to site-packages folder in Abaqus installation directory;And all attached programs need to be placed in the same folder, all attached programs include Excel file, and Python script file;When user inputs parameter information, the frame column without input mark in table should not be changed at will, and input should be carried out in specified area.

[0013] Preferably, in step (3), TGCW parameter information in Excel file is read by Python script file, and comprises the following steps:

[0014] (31) all objects in Abaqus module are imported to provide access to default model database;Symbol constant module is imported to call predefined constants in Abaqus;The third line of code imports all objects in caeModules module;RegionToolset is imported to call Region function;Library function in math is imported;Xlrd2 is imported to call third-party library function therein to carry out Excel reading work;

[0015] (32) Open the Excel file to read parameter data by the open_workbook() function; get the target worksheet data by name using the sheet_by_name() function; use the for loop function to parameterize and analyze multiple sets of data in the worksheet, and after the loop ends, use the release_resources() function to release the data and delete the temporary data.

[0016] Preferably, in step (4), the TGCW-related model information defined by the Python script includes the following steps:

[0017] (41) Initialize the name and path of the result report and perform initial calculation on the number of waves contained in the corrugated steel web CSW. Calculate the number of complete waves and the residual incomplete wavelength contained in the CSW according to the length of the corrugated steel web beam and the wavelength of the corrugated plate obtained from the input file. Determine the geometric parameters of the residual corrugated plate through the ratio of the residual wavelength to the complete wavelength;

[0018] (42) Initialize the model, and the model naming rule is 'TGCW_parameter analysis variable_iteration number'; if no iteration analysis is needed, the model name is TGCW_normal_1;

[0019] (43) Create component objects: concrete upper flange plate (Con_up), concrete bottom plate (Con_bottom), corrugated steel web steel beam top and bottom plate (CSW_top, CSW_bottom), corrugated steel web steel beam web and residual web (CSW_depart, CSW_depart2);

[0020] (44) Define model information such as materials and cross-sectional properties;

[0021] (45) Assemble the three-dimensional finite element model of the corrugated steel web composite beam bridge: concrete upper flange plate × 2, concrete bottom plate (Con_bottom), corrugated steel web steel beam × 2, ordinary steel reinforcement skeleton, prestressed steel reinforcement; wherein the corrugated steel web steel beam is composed of top and bottom plates, webs × n and residual webs; after creating instances of each component by calling the Instance() constructor function, continue to call the translate() function to move the component instances to the specified position;

[0022] (46) Create analysis steps according to user input parameters, apply boundary conditions and external loads;

[0023] (47) Use the Tie function to establish the constraint coupling relationship between the concrete top and bottom plates and the corrugated steel web steel beam;

[0024] (48) Divide the mesh, submit and run the analysis job.

[0025] Preferably, in step (5), the user can run the Python script on the Abaqus platform in multiple ways:

[0026] (51) Run the Python script file in the Abaqus Command command line window:

[0027] abaqus cae script=TGCW_establish.py

[0028] abaqus cae noGUI=TGCW_establish.py

[0029] If the user has no need to supplement the model, it is recommended to use the third command to realize automatic pre-processing operation;

[0030] (52) Run the script from the Abaqus / Cae File menu bar Run Script;

[0031] (53) Run the script from the command line interface, enter execfile('TGCW_establish.py') in the command line.

[0032] The beneficial effects of the present application are: (1) Without relying on third-party modeling software (Autocad, Catia, etc.), the modeling work of complex models can be completed, which not only reduces the learning cost of the user but also avoids errors that may occur due to incompatibility between different software during model import; (2) The automatic modeling process can be realized, which greatly saves the user's time, and only needs the user to input the related modeling parameters into the Excel table and run the Python script file to complete the modeling work; (3) It is convenient to realize parameterized analysis of the corrugated steel web composite channel beam; when the user needs to modify the parameters in batches, only needs to modify the related parameters of the Excel table in batches and run the Python script file to complete batch modeling. BRIEF DESCRIPTION OF DRAWINGS

[0033] Figure 1 The method flowchart of the present application.

[0034] Figure 2(a) is an elevation view of the corrugated steel web composite channel beam of the present application (unit: mm).

[0035] Figure 2(b) is a cross-sectional view of the corrugated steel web composite channel beam of the present application (unit: mm).

[0036] Figure 2(c) is a cross-sectional view of the corrugated steel web composite channel beam of the present application (unit: mm).

[0037] Figure 3The schematic diagram of the general reinforcement arrangement of the corrugated steel web combined channel girder of the application (unit: mm).

[0038] Figure 4 The schematic diagram of the Excel platform user input interface of the application.

[0039] Figure 5 The schematic diagram of the TGCW three-dimensional model and each component of the model of the application.

[0040] Figure 6 The schematic diagram of the CSW three-dimensional model and each component assembly of the model of the application.

[0041] Fig. 7(a) is a schematic diagram of the contact condition of the TGCW three-dimensional model of the application.

[0042] Fig. 7(b) is a schematic diagram of the boundary condition of the TGCW three-dimensional model of the application. DETAILED DESCRIPTION

[0043] A 25m standard span simply supported corrugated steel web combined channel girder is parameterized modeled, and on this basis, the fundamental frequency analysis and mode shape extraction are carried out.

[0044] The construction of the test girder TGCW in this example is shown in Fig. 2(a) and Fig. 2(b).

[0045] The diameter of the general reinforcement arranged in the test girder is 8mm, arranged in the concrete bottom plate and the upper flange plate. The transverse and vertical spacing of the longitudinal reinforcement in the concrete bottom plate is 150mm and 70mm respectively, and the transverse and vertical spacing of the longitudinal reinforcement in the concrete top plate is 120mm and 60mm respectively. The spacing of the stirrup along the length direction of the girder is 150mm. The specific construction of the general reinforcement is shown in Figure 3 .

[0046] The size of the corrugated steel web steel girder is basically the same, the steel girder is 26000mm long and 330mm high. Referring to the size of the commonly used corrugated steel web of the bridge, the wavelength of the corrugated steel web is 1200mm, the wave height is 200mm, and the steel plate thickness is 4mm, as shown in the figure. The upper flange plate of the steel girder is 40mm thick and 520mm wide, and is welded to the corrugated steel web as shown in Fig. 2(c).

[0047] The concrete used in the test is made on site in the laboratory, and the main materials involved include cement, river sand, gravel, and water reducing agent. The general reinforcement includes Q235 steel reinforcement with a nominal diameter of 8mm and 12mm, and the prestressed reinforcement uses 1x7 standard type-15.2-1860-Ⅱ-GB / T 5224-1995 steel strand.

[0048] The prestressed steel strand adopts 7-strand low-relaxation steel strand with a nominal diameter of 15.24 mm and a theoretical tensile strength of 1860 MPa. The ordinary steel bar adopts a steel bar with a nominal diameter of 8 mm and a theoretical yield strength of 235 MPa. The corrugated steel web steel beam includes a web, a flange plate and an end cross plate, all of which are made of weathering steel Q345q with a theoretical yield strength of 345 MPa. The test specimens for material tensile test are taken from the web, the flange plate, the end cross plate and the ordinary steel bar. The steel plates are made into standard test components, and the steel bars are cut to a certain length. The tensile strength and yield strength of them are tested respectively.

[0049] As Figure 1 shown, an Excel-based Abaqus parameterized modeling and analysis method for corrugated steel web composite channel beams, comprising the following steps:

[0050] Step one: obtain the size, material and other parameter information of TGCW by consulting project data or actual measurement; and determine the type of finite element analysis to be performed as self-vibration analysis.

[0051] Step two: configure the program running environment.

[0052] (1) download the attached program of the present application, and install the third-party library function xlrd2;

[0053] (2) after the installation of xlrd2 is completed, copy the corresponding library function file in the site-packages folder in the system Python to the site-packages folder in the Abaqus installation directory (for example, if Abaqus2019 is installed in the D disk, then the D:\SIMULIA\Abauqs\2019\Tools\SMApy\python2.7\Lib\site-packages folder is the target position for copying), and thus the corresponding plug-in installation is completed;

[0054] Note: all programs (Excel files, Python script files, etc.) attached to the present application need to be placed in the same folder.

[0055] (3) fill the parameter information obtained in step one into the finite element analysis parameter Excel input table of the corrugated steel web composite channel beam, as shown in Figure 4 .

[0056] Step three: read the parameter information in the Excel file through the Python script file.

[0057] Further, the step three specifically includes the following steps.

[0058] (1) import the Abaqus related modules, and the code field is located in the 1-6 lines of the script file;

[0059] from abaqus import*

[0060] from abaqusConstants import*

[0061] from caeModules import*

[0062] import regionToolset

[0063] import math

[0064] import xlrd2

[0065] The first line of code imports all objects in the Abaqus module to provide access to the default model database;

[0066] The second line of code imports the symbol constant module to call the predefined constants in Abaqus;

[0067] The third line of code imports all objects in the caeModules module;

[0068] The fourth line of code imports regionToolset to call the Region function;

[0069] The fifth line of code imports the library function in math;

[0070] The sixth line of code imports xlrd2 to call the third-party library function therein for Excel reading work.

[0071] (2) Obtain parameter information in the Excel file;

[0072]

[0073]

[0074] The open_workbook() function opens the Excel file to read data; the third line of code obtains the target worksheet data by name; the fifth line of code and the following code block use the for loop function to parameterize analysis of multiple sets of data in the worksheet; after the loop ends, the release_resources() function is used to release the data and delete the temporary data.

[0075] Further, the step four specifically includes the following steps;

[0076] (1) Initialize variables and perform some preliminary calculations;

[0077] The system initializes variables for the name and path of the results report and performs initial calculations for the number of waves contained in the corrugated steel web (CSW). Based on the CSW beam length and corrugated plate wavelength obtained from the input file, the system calculates the number of complete waves and residual incomplete wavelengths contained in the CSW. The geometric parameters of the residual corrugated plate are determined by the ratio of the residual wavelength to the complete wavelength.

[0078] (2) Initialize the model;

[0079] modelname='TGCW_'+sheet_name.name+'-'+repr(iteration_count)

[0080] tgcwmodel=mdb.Model(name=modelname)

[0081] The model naming convention is 'TGCW_parameter analysis variable_iteration number'; if iterative analysis is not required, the model name is TGCW_normal_1; the iteration number is appended to the model name, so that the Abaqus / CAE model tree will contain all created models at the end of all simulation runs.

[0082] (3) Draw a 2D sketch, create the component, and add an extrusion feature to the sketch;

[0083] First, the ConstrainedSketch constructor is called to create a sketch object and draw the cross-sectional sketch; second, a part object is created and the BaseSolidExtrude constructor is called to add the extrusion feature and specify the extrusion length.

[0084] Perform the above two steps on the component objects to be created in the model: concrete upper flange plate (Con_up), concrete bottom plate (Con_bottom), corrugated steel web steel beam top and bottom plates (CSW_top, CSW_bottom), corrugated steel web steel beam web and residual web plate (CSW_depart, CSW_depart2).

[0085] (4) Define model information such as material and section properties;

[0086] Call the Material function to create a material and define its properties: elastic modulus, density, Poisson's ratio, etc.

[0087] tgcwmodel.Material(name='Q235')

[0088] tgcwmodel.materials['Q235'].Density(table=((CSW_density,),))

[0089] tgcwmodel.materials['Q235'].Elastic(table=((CSW_E, CSW_p),))

[0090] …

[0091] Create section and assign section property to the part.

[0092] tgcwmodel.HomogeneousSolidSection(name='Sec_Con_top', material='C40',

[0093] thickness=None)

[0094] p = tgcwmodel.parts['Con_up']

[0095] region=(p.cells,)

[0096] p.SectionAssignment(region=region, sectionName='Sec_Con_top', offset=0.0,

[0097] offsetType=MIDDLE_SURFACE, offsetField="",

[0098] thicknessAssignment=FROM_SECTION)

[0099] …

[0100] (5) Model assembly;

[0101] As Figure 5 shown, the three-dimensional finite element model of the corrugated steel web composite beam bridge includes: concrete upper flange plate x 2, concrete bottom plate (Con_bottom), corrugated steel web steel beam x 2, ordinary steel reinforcement framework, and prestressed steel. The corrugated steel web steel beam is composed of top and bottom plates, web x n, and residual web, as shown in Figure 6 .

[0102] After the code calls the Instance constructor to create an instance of each part, it continues to call the translate function to move the part instance to the specified location.

[0103] a1 = tgcwmodel.rootAssembly

[0104] p = tgcwmodel.parts['Con_bottom']

[0105] a1.Instance(name='Con_bottom-1', part=p, dependent=ON)

[0106] a1.translate(instanceList=('Con_bottom-1',), vector=(0.0, Con_bottom_h / 2, 0.0))

[0107] …

[0108] (6) According to the user input parameters to create analysis steps, such as boundary conditions and external loads, as shown in Figures 7(a) and 7(b).

[0109] According to the user input parameters in the Excel table to create static or frequency analysis steps;

[0110]

[0111] By defining the end face where the end is located, and creating constraints at the beam end;

[0112] …

[0113] tgcwmodel.DisplacementBC(name='BC-1', createStepName='Frequency',

[0114] region=region, u1=0.0, u2=0.0, u3=0.0, ur1=UNSET, ur2=0.0, ur3=0.0,

[0115] amplitude=UNSET, fixed=OFF, distributionType=UNIFORM, fieldName=",

[0116] localCsys=None)

[0117] …

[0118] (7) Use the Tie function to establish the constraint coupling relationship between the concrete top and bottom plates and the corrugated steel web steel beam, as shown in Figures 7(a) and 7(b).

[0119] a = tgcwmodel.rootAssembly

[0120] region1 = a.instances['CSW_minus'].surfaces['top']

[0121] region2 = a.instances['Con_up_minus'].surfaces['bottom']

[0122] tgcwmodel.Tie(name='CSW_Con_top_minus', master=region1, slave=region2,

[0123] positionToleranceMethod=COMPUTED, adjust=ON, tieRotations=ON, thickness=ON)

[0124] …

[0125] (8) Divide the mesh, submit and run the analysis job;

[0126] Use the regions parameter to set the mesh seed for the entire part instance, specify the seed size, unit type; Call the generateMesh function to divide the mesh.

[0127] p = tgcwmodel.parts['Con_bottom']

[0128] cells = p.cells.findAt(((0.0001,-0.0001,Con_bottom_l / 2),))

[0129] pickedRegions = (cells,)

[0130] p.setElementType(regions=pickedRegions, elemTypes=(elemType1, elemType2,

[0131] elemType3))

[0132] p.seedPart(size=Con_seed_size, deviationFactor=0.1, minSizeFactor=0.1)

[0133] p.generateMesh()

[0134] …

[0135] (9) Create, and submit the running analysis job;

[0136] jobname = modelname + '_frequency'

[0137] mdb.Job(name = jobname, model = modelname, description = "",

[0138] type = ANALYSIS, atTime = None, waitMinutes = 0, waitHours = 0, queue = None,

[0139] memory = 90, memoryUnits = PERCENTAGE, getMemoryFromAnalysis = True,

[0140] explicitPrecision = SINGLE, nodalOutputPrecision = SINGLE, echoPrint = OFF,

[0141] modelPrint = OFF, contactPrint = OFF, historyPrint = OFF, userSubroutine = "",

[0142] scratch = "", resultsFormat = ODB, multiprocessingMode = DEFAULT, numCpus = 1,

[0143] numGPUs = 0)

[0144] mdb.jobs[jobname].submit(consistencyChecking = OFF)

[0145] …

[0146] Step five: generate the TGCW finite element model and perform simulation analysis corresponding to the input parameters; the user can run the Python script in multiple ways on the Abaqus platform (see the summary of the invention); if the user has no need to supplement the model, it is recommended to run the command 'abaqus cae noGUI = TGCW_establish.py' in the Abaqus Command command line window to realize automatic pre-processing operation.

Claims

1. An Excel-based Abaqus parametric modeling and analysis method for corrugated steel web composite channel beam, characterized in that, Comprise the following steps: (1) The user determines the TGCW size, material, and the type of finite element analysis required by the user from the project data; the finite element analysis type includes frequency mode analysis, static load analysis, and dynamic load analysis parameter information; (2) The program running environment is configured, the third-party library function xlrd2 is installed, and the parameter information obtained in step (1) is filled into the corrugated steel web combined channel beam finite element analysis parameter Excel input table; (3) The TGCW parameter information in the Excel file is read through the Python script file; Specifically comprising the following steps: (31) Import all objects in the Abaqus module to provide access to the default model database; import the symbol constant module for calling predefined constants in Abaqus; the third line of code imports all objects in the caeModules module; import regionToolset for calling Region functions; import library functions in math; Import xlrd2 to call third-party library functions to read Excel; (32) Open the Excel file to read parameter data through the open_workbook function; use the sheet_by_name function to obtain target worksheet data by name; use the for loop function to perform parameterized analysis on multiple sets of data in the worksheet, and after the loop ends, use the release_resources function to release the data and delete the temporary data; (4) Define TGCW related model information through the Python script, including Part, Property, Assembly, Mesh, Step, Interaction, Load, and Job; specifically comprising the following steps: (41) Initialize the name and path of the result report execution variable, and perform initial calculation on the number of waves contained in the corrugated steel web CSW; calculate the number of complete waves and the residual incomplete wavelength of the corrugated steel web CSW according to the length of the corrugated steel web beam and the wavelength of the corrugated plate obtained from the input file, and determine the geometric parameters of the residual corrugated plate through the ratio of the residual wavelength to the complete wavelength; (42) Initialize the model, the model naming rule is 'TGCW_parameter analysis variable_iteration number'; if no iteration analysis is required, the model name is TGCW_normal_1; (43) Create component objects: concrete upper flange plate Con_up, concrete bottom plate Con_bottom, corrugated steel web steel beam top and bottom plate CSW_top, CSW_bottom, corrugated steel web steel beam web, and residual web CSW_depart, CSW_depart2; (44) Define material, cross-section attribute model information; (45) Assemble the three-dimensional finite element model of the corrugated steel web composite beam bridge: concrete top flange plate × 2, concrete bottom plate Con_bottom, corrugated steel web steel beam × 2, ordinary steel reinforcement framework, prestressed steel; wherein the corrugated steel web steel beam is composed of top and bottom plates, web × n and residual web; after calling the Instance constructor to create an instance of each component, continue to call the translate function to move the component instance to the specified location; (46) Create an analysis step according to the user input parameters, apply boundary conditions and external loads; (47) Use the Tie function to establish the constraint coupling relationship between the concrete top and bottom plates and the corrugated steel web steel beam; (48) Divide the mesh, submit and run the analysis job; (5) Run the Python script using the Abaqus platform, generate the TGCW finite element model and perform simulation analysis corresponding to the input parameters.

2. The Excel-based Abaqus parametric modeling and analysis method for corrugated steel web composite box-girder as claimed in claim 1, wherein, In step (2), after the xlrd2 installation is completed, the user needs to copy the corresponding library function file in the site-packages folder in the system Python to the site-packages folder in the Abaqus installation directory; and all attached programs, including Excel files and Python script files, are placed in the same folder; When the user inputs parameter information, the box column without input flag in the table cannot be changed arbitrarily, and input must be performed in the specified area.

3. The Excel-based Abaqus parametric modeling and analysis method for corrugated web composite box-girder as claimed in claim 1, wherein, In step (5), the user runs the Python script on the Abaqus platform in multiple ways: (51) Run the Python script file in the Abaqus Command command line window: abaqus cae script=TGCW_establish.py abaqus cae noGUI=TGCW_establish.py If the user has no need for supplementary operations on the model, use abaqus cae noGUI=TGCW_establish.py to realize automatic pre-processing operations; (52) Run the script from the Abaqus / Cae File menu bar Run Script; (53) Run the script from the command line interface, enter execfile(‘TGCW_establish.py’) in the command line.

Images