A method for predicting residual stress and deformation of T-joint welding

Abstract

The present invention provides a T-joint welding residual stress and deformation prediction method, comprising the following steps: S1: establishing a steel plate T-joint welding three-dimensional model in ABAQUS, assigning material properties, meshing the three-dimensional model, defining the welding finite element simulation trend, setting the welding step, cooling step and boundary conditions; S2: writing a double-wire MAG welding program language through the Dlux program in Visual Studio software; S3: submitting the coordinate information of the meshed T-joint welding three-dimensional model in ABAQUS, and the set material properties, welding step, cooling step and boundary conditions to the Visual Studio software, activating the Visual Studio software and then running the Dflux language controlled by the Inter Fortran program to perform welding finite element analysis. The T-joint welding residual stress and deformation prediction method provided by the present invention can more intuitively explore relevant data such as welding temperature field, post-welding residual stress, welding angle deformation, etc., and can greatly save production costs.

Description

Technical Field

[0001] The invention belongs to the technical field of welding, and in particular relates to a method for predicting residual stress and deformation of T-joint welding. Background Art

[0002] When T-joints are actually welded on site, only the stable state of the arc and the formation of the weld can be observed. The welding temperature and residual stress after welding cannot be measured. There are errors in the manual measurement of welding angle deformation, and the actual test is repeated many times. The existing T-joint welding has the following defects:

[0003] 1. During welding, there is no special equipment to measure the temperature field distribution of the weldment and the residual stress after welding. Although the welding angle deformation can be measured, there are human subjective errors. At present, only physical fracture, macroscopic and hardness tests can be performed on T-joints, and the data cannot meet the actual needs at all.

[0004] 2. Since the magnitude of stress generated during welding and the degree of angular deformation are not clear, the problems of angular deformation and fatigue failure after welding of T-joints have not been effectively solved.

[0005] 3. The degree of welding angle deformation can be reduced through actual field tests, but due to different plate thicknesses and different sizes of weld legs, the number of tests is large, the time spent is long, and the production cost is high.

[0006] Patent publication number CN 112257314 A discloses a residual stress prediction method for heat-resistant steel welding, including: step 1, establishing a geometric model; step 2, applying a double ellipsoid heat source to gradually calculate the temperature field; step 3, calculating the tissue field; step 4, adding the latent heat generated by the martensitic phase transformation as an endogenous heat source to the temperature field calculation; step 5, using the obtained temperature field as a predefined field, inputting the mechanical property parameters of the material, and calculating the stress field. During the calculation, the volume change, yield strength change and phase transformation plasticity caused by the martensitic phase transformation are calculated. The influence of the stress field.

[0007] Although this scheme can predict the residual stress of welding, it is not applicable to the prediction of welding stress of T-joints. Summary of the invention

[0008] In view of the deficiencies in the prior art, the present invention provides a method for predicting residual stress and deformation of T-joint welding.

[0009] The present invention provides a method for predicting residual stress and deformation of T-joint welding, comprising the following steps:

[0010] S1: Establish a three-dimensional model of steel plate T-joint welding in ABAQUS, assign material properties, mesh the three-dimensional model, define the welding finite element simulation direction, and set the welding step, cooling step and boundary conditions;

[0011] S2: Write the double-wire MAG welding program language through the Dlux program in Visual Studio software;

[0012] S3: Submit the coordinate information of the meshed three-dimensional model of the T-joint welding in the ABAQUS, as well as the set material properties, welding steps, cooling steps and boundary conditions to the Visual Studio software, activate the Visual Studio software and then run the Dflux language controlled by the Inter Fortran program to perform welding finite element analysis.

[0013] Preferably, the specific steps of writing the twin-wire MAG welding program language are as follows:

[0014] S21: setting welding parameters and double ellipsoid molten pool size, and obtaining a heat source equation according to the welding parameters and the double ellipsoid molten pool size, wherein the double ellipsoid molten pool size includes the front half zone size, the rear half zone size, the molten depth and the molten width, and the welding parameters include welding current, welding voltage, welding power coefficient, welding speed, welding starting point and welding direction;

[0015] S22: rotating the double ellipsoid molten pool by a certain angle so that the direction of the double ellipsoid molten pool is perpendicular to the weld surface;

[0016] S23: define the welding starting point position, obtain the movement trajectory equation according to the welding starting point position and the angle of rotation of the double ellipsoid molten pool in step S3, and obtain the heat source movement equation according to the heat source equation and the movement trajectory equation.

[0017] Preferably, the heat source movement equation satisfies the following conditions:

[0018] FLUX(1)=heat1*shape1;

[0019] FLUX(2)=heat2*shape2;

[0020] Where FLUX(1) is the heat flux of the first half, FLUX(2) is the heat flux of the second half, heat1 is the heat input source of the first half, heat2 is the heat input source of the second half, hape1 is the moving trajectory of the first half, and shape2 is the moving trajectory of the second half.

[0021] Preferably, the heat source equation satisfies the following conditions:

[0022] heat1=6.0*sqrt(3.0)*q / (a*b*c*PI*sqrt(PI))*f;

[0023] heat2=6.0*sqrt(3.0)*q / (aa*b*c*PI*sqrt(PI))*(2.0-f);

[0024] Wherein, q is the heat input size, a is the size of the front half of the heat source, b is the width of the heat source, c is the depth of the heat source, PI is the coordinate information of the welding three-dimensional model, and f is the heat proportion of the front hemisphere;

[0025] Preferably, q satisfies the following conditions:

[0026] q = u*i*effi;

[0027] Wherein, u is the welding voltage, i is the welding current, and effi is the welding efficiency.

[0028] Preferably, the movement trajectory equation satisfies the following conditions:

[0029] shape1=exp(-3.0*((x-x0)*cos(betahudu1)+(y-y0)*sin(betahudu1))*2 / a*2-3.0*( (y-y0)*cos(betahudu1)-(x-x0)*sin(betahudu1))*2 / b*2-3.0*((z-z0+d))*2 / c*2);

[0030] shape2=exp(-3.0*((x-x10)*cos(betahudu2)+(y-y10)*sin(betahudu2))*2 / aa*2-3.0* ((y-y10)*cos(betahudu2)-(x-x0)*sin(betahudu2))*2 / b*2-3.0*((z-z0+d))*2 / c*2);

[0031] Wherein, x, y, z are coordinate information, betahudu1 is the first half of the arc, and betahudu2 is the second half of the arc;

[0032] Preferably, betahudu1 and betahudu2 satisfy the following conditions:

[0033] betahudu1=beta1 / 180.0*PI;

[0034] betahudu2=beta2 / 180.0*PI;

[0035] Among them, beta1 and beta2 are the rotation angles of the double ellipsoid melt pool;

[0036] Preferably, the mesh division is a sweeping mesh division, the mesh density of the T-joint welding three-dimensional model in the weld area and the heat affected zone is greater than the mesh density in other areas, and the mesh is set as a temperature-displacement coupling unit reduced integration.

[0037] Preferably, the material properties include yield stress, expansion coefficient, conductivity and elastic modulus corresponding to different temperatures; the boundary conditions include setting the ambient temperature, convection heat transfer coefficient and radiation coefficient;

[0038] The welding step is to simulate the welding time used in the actual welding process; the cooling step is to simulate the cooling time used in the actual welding process.

[0039] The present invention also provides a computer-readable storage medium on which a computer program is stored. When the program is executed by a processor, the method for predicting residual stress and deformation of T-joint welding is implemented.

[0040] The T-joint welding residual stress and deformation prediction method provided by the present invention has the following beneficial effects: it can more intuitively explore relevant data such as welding temperature field, post-welding residual stress, welding angle deformation, etc., which can greatly save production costs.

[0041] The data of the temperature field, welding residual stress field and welding deformation at any position of the weldment during welding are supplemented with the time-varying data, and the relevant characteristic curves can be drawn, solving the problem of insufficient data in conventional physical property tests. For the welding angle deformation data of different plate thicknesses and different welding angles of T-joints, there is no need to spend a lot of manpower and material resources to conduct related tests, saving costs and improving work efficiency.

[0042] The relevant data of temperature field during welding, welding residual stress and welding angular deformation can be obtained. The analysis level of welded test plate is deepened and easy to understand, which reduces the amount of welding process testing and reduces the test cost, and has high practical value.

[0043] It provides a good reference standard for improving the quality of T-joint welds and optimizing welding parameters. BRIEF DESCRIPTION OF THE DRAWINGS

[0044] The above and other purposes, features and advantages of the present invention will become more clear by more specific description of the preferred embodiments of the present invention shown in the accompanying drawings. The same reference numerals indicate the same parts in all the accompanying drawings, and the accompanying drawings are not deliberately scaled to the actual size, and the focus is on illustrating the main purpose of the present invention.

[0045] Figure 1A schematic diagram of a three-dimensional model of a T-joint welding provided in an embodiment of the present invention;

[0046] Figure 2 A temperature field distribution cloud diagram of two welding passes of a T-joint provided in an embodiment of the present invention;

[0047] Figure 3 A schematic diagram of a temperature variation curve of a characteristic point at the center of a weld provided by an embodiment of the present invention;

[0048] Figure 4 A cloud diagram of the equivalent residual stress distribution of two welds of a T-joint provided in an embodiment of the present invention;

[0049] Figure 5 A cloud diagram of longitudinal residual stress distribution of two welds of a T-joint provided in an embodiment of the present invention;

[0050] Figure 6 A schematic diagram of path 1 and path 2 in a two-pass welding simulation of a T-joint provided in an embodiment of the present invention;

[0051] Figure 7 A curve diagram of the transverse residual stress distribution along path 1 of two passes of welding of a T-joint provided in an embodiment of the present invention;

[0052] Figure 8 A graph showing the longitudinal residual stress distribution curve of two-pass welding of a T-joint along path 2 provided in an embodiment of the present invention;

[0053] Fig. 9 A distribution cloud diagram of welding deformation of two passes of a T-joint weld provided in an embodiment of the present invention (magnified 10 times);

[0054] Fig.10 A distribution curve diagram of angular deformation of two-pass welding of a T-joint along path 2 provided in an embodiment of the present invention;

[0055] Fig.11 A schematic diagram of the collaboration between ABAQUS, Visual Studio and Dlux programs provided in an embodiment of the present invention.

[0056] Fig.12 A flow chart of a method for predicting residual stress and deformation of T-joint welding provided by an embodiment of the present invention;

[0057] Fig.13 A schematic diagram of a specific flow chart of step S1 provided in an embodiment of the present invention;

[0058] Fig.14 A schematic diagram of a language flow for writing a double-wire MAG welding program provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0059] To facilitate understanding of the present invention, a more comprehensive description of the present invention will be given below with reference to the relevant drawings.

[0060] It should be noted that when an element is considered to be "connected" to another element, it may be directly connected to another element and integrated therewith, or there may be an intermediate element at the same time. The terms "mounted", "a part", "another part" and similar expressions used herein are for illustrative purposes only.

[0061] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as those generally understood by those skilled in the art. The terms used in the specification herein are only for the purpose of describing specific embodiments and are not intended to limit the present invention. The term "and / or" used herein includes any and all combinations of one or more related listed items.

[0062] refer to Figures 1 to 14 The embodiment of the present invention provides a method for predicting residual stress and deformation of T-joint welding. The method for predicting residual stress and deformation of T-joint welding provided by the embodiment of the present invention includes the following steps (refer to Fig.12 ):

[0063] S1: Establish a three-dimensional model of steel plate T-joint welding in ABAQUS, assign material properties, mesh the three-dimensional model, define the welding finite element simulation direction, and set the welding step, cooling step and boundary conditions.

[0064] S2: Write the double-wire MAG welding program language through the Dlux program in Visual Studio software. During the welding process, two heat sources are used simultaneously, two by two, and two by two parallel;

[0065] S3: Submit the coordinate information of the meshed three-dimensional model of the T-joint welding in ABAQUS, as well as the set material properties, welding steps, cooling steps and boundary conditions to the Visual Studio software, activate the Visual Studio software and run the Dflux language controlled by the Inter Fortran program to perform welding finite element analysis, refer to Fig.11 .

[0066] The T-joint welding residual stress and deformation prediction method provided in the embodiment of the present invention is submitted to the ABAQUS software, at which time the Visual Studio software is activated to run the Dflux language controlled by the Inter Fortran program, and then fed back to the ANAQUS software, thereby realizing the welding finite element analysis.

[0067] After the above steps, the finite element analysis animation of the welding temperature field, residual stress field, and angular deformation of the T-joint will be presented in the ABAQUS software after the calculation is completed. At the same time, the numerical curves of these three variables at different times and positions can also be obtained to prepare for welding analysis.

[0068] In a preferred embodiment, reference Fig.13 , step S1 is specifically as follows:

[0069] S11: Establish a three-dimensional model of steel plate T-joint welding in ABAQUS finite element analysis software. Figure 1 , the model size is consistent with the actual welding steel size to prevent the calculated results from being different from the actual results due to different sizes.

[0070] S12: Define the material properties of the actual steel T-joint for the model drawn in step S11. Make the material properties assigned to the model consistent with the material properties of the actual steel plate. The material properties include the corresponding yield stress, expansion coefficient, conductivity and elastic modulus at different temperatures. In this embodiment, Q235 steel is taken as an example, and the material properties of Q235 steel are shown in Table 1.

[0071] Table 1

[0072]

[0073] S13: Mesh the model drawn in step S11. In meshing, the more meshes there are, the lower the speed of computer operation is. Therefore, the principle of meshing is to reduce the number of meshes under the condition of ensuring normal operation. Therefore, a transitional uniform meshing method is adopted, that is, a sweeping meshing method, that is, the meshes in the weld area and the heat-affected zone are relatively dense, and the meshes in the rest of the weldment are sparse after a mesh transition. The mesh density in the weld area and the heat-affected zone is greater than the mesh density in other areas; the mesh is set as a temperature-displacement coupling unit reduced integral, and the welding finite element simulation direction is defined. The temperature, stress, and deformation during the welding process can be obtained later.

[0074] S14: Set the welding analysis step of the model drawn in S11, which mainly includes welding step and cooling step. The welding step is the welding time used in simulating the actual welding process; the cooling step is the cooling time used in simulating the actual welding process.

[0075] S15: Set the boundary conditions of the model drawn in step S11. The boundary conditions include setting the ambient temperature, convection heat transfer coefficient and radiation coefficient. Since there is heat exchange between the weld and the base material, heat convection between the base material and the air during welding, and the actual ambient temperature on site have a certain influence on the test results. Therefore, this step is to simulate the heat dissipation of the welding test plate during the actual on-site welding process.

[0076] In a preferred embodiment, in step S2, the programming language of the twin-wire MAG welding mainly includes welding parameters, double ellipsoid molten pool size and welding angle.

[0077] refer to Fig.14 The specific steps for writing the double-wire MAG welding program language are as follows:

[0078] S21: setting welding parameters and double ellipsoid molten pool size, and obtaining a heat source equation according to the welding parameters and the double ellipsoid molten pool size.

[0079] In this embodiment, the size of the double ellipsoid molten pool includes the size of the front half area, the size of the rear half area, the molten depth and the molten width, so as to ensure that the size of the molten pool is consistent with the size of the molten pool formed in the actual welding process.

[0080] The welding parameters include welding current (determined according to the welding current in the actual welding process), welding voltage (determined according to the welding voltage in the actual welding process), welding power coefficient (different welding methods have different welding power coefficients, and its main function is to determine the welding power in the actual welding process), welding speed (determined according to the size of the model drawn in step S11 and the welding time of the welding step in step S14), welding starting point and welding direction (the same as the actual welding test).

[0081] The heat source equation satisfies the following conditions:

[0082] heat1=6.0*sqrt(3.0)*q / (a*b*c*PI*sqrt(PI))*f;

[0083] heat2=6.0*sqrt(3.0)*q / (aa*b*c*PI*sqrt(PI))*(2.0-f);

[0084] Among them, heat1 is the heat input source in the front half, heat2 is the heat input source in the back half, q is the heat input size, a is the size of the front half of the heat source, b is the width of the heat source, c is the depth of the heat source, PI is the coordinate information of the welding three-dimensional model, and f is the heat ratio of the front hemisphere.

[0085] q (heat input size) satisfies the following conditions:

[0086] q = u*i*effi;

[0087] Among them, u is the welding voltage, i is the welding current, and effi is the welding efficiency.

[0088] S22: Rotate the double ellipsoid molten pool by a certain angle so that the direction of the double ellipsoid molten pool is perpendicular to the weld surface. Since the model drawn in step S11 is a T-joint, and the direction of the molten pool in step S21 is perpendicular to the bottom plate of the T-joint at the weld, it is necessary to make the molten pool perpendicular to the weld direction, so the direction of the molten pool is rotated so that the molten pool is perpendicular to the weld.

[0089] S23: define the welding starting point position, obtain the movement trajectory equation according to the welding starting point position and the angle of rotation of the double ellipsoid molten pool in step S3, and obtain the heat source movement equation according to the heat source equation and the movement trajectory equation.

[0090] The heat source movement equation satisfies the following conditions:

[0091] FLUX(1)=heat1*shape1;

[0092] FLUX(2)=heat2*shape2;

[0093] Where FLUX(1) is the heat flux of the first half, FLUX(2) is the heat flux of the second half, heat1 is the heat input source of the first half, heat2 is the heat input source of the second half, hape1 is the moving trajectory of the first half, and shape2 is the moving trajectory of the second half.

[0094] The moving trajectory equation satisfies the following conditions:

[0095] shape1=exp(-3.0*((x-x0)*cos(betahudu1)+(y-y0)*sin(betahudu1))*2 / a*2-3.0*( (y-y0)*cos(betahudu1)-(x-x0)*sin(betahudu1))*2 / b*2-3.0*((z-z0+d))*2 / c*2);

[0096] shape2=exp(-3.0*((x-x10)*cos(betahudu2)+(y-y10)*sin(betahudu2))*2 / aa*2-3.0* ((y-y10)*cos(betahudu2)-(x-x0)*sin(betahudu2))*2 / b*2-3.0*((z-z0+d))*2 / c*2);

[0097] Among them, x, y, z are coordinate information, betahudu1 is the first half of the arc, and betahudu2 is the second half of the arc;

[0098] betahudu1 and betahudu2 meet the following conditions:

[0099] betahudu1=beta1 / 180.0*PI;

[0100] betahudu2=beta2 / 180.0*PI;

[0101] Among them, beta1 and beta2 are the rotation angles of the double ellipsoid melt pool.

[0102] The T-joint welding residual stress and deformation prediction method provided in this embodiment can finally be obtained in ABAQUS:

[0103] 1. Temperature field distribution cloud diagram of two welding passes of T-joint, reference Figure 2 , you can intuitively see the temperature changes of the weld and its plate during the welding process.

[0104] 2. Temperature change curve of characteristic point in the center of weld, reference Figure 3 , the temperature variation range of the weld center at different times and positions can be obtained.

[0105] 3. T-joint two-pass welding equivalent residual stress distribution cloud diagram, reference Figure 4 , the equivalent residual stress distribution of the entire welding test plate can be intuitively seen.

[0106] 4. Longitudinal residual stress distribution cloud diagram of two welds of T-joint, reference Figure 5 , we can see the stress distribution in the longitudinal direction of the weld.

[0107] 5.Reference Figure 6 , Schematic diagram of path 1 and path 2 in the two-pass welding simulation of the T-joint. The corresponding paths can be used to draw stress and deformation curves at different positions.

[0108] 6. Transverse residual stress distribution curve of two passes of T-joint welding along path 1, refer to Figure 7 , the magnitude of the transverse residual stress at different locations can be compared.

[0109] 7. The longitudinal residual stress distribution curve of the two-pass welding of the T-joint along the path 2, refer to Figure 8 , the magnitude of longitudinal residual stress at different positions can be compared.

[0110] 8. T-joint two-pass welding deformation distribution cloud diagram (magnified 10 times), reference Fig. 9 , the deformation of the weldment can be seen intuitively.

[0111] 9. Angular deformation distribution curve of two passes of T-joint welding along path 2, reference Fig.10 , compare the deformation size at different positions.

[0112] The present invention also provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the method for predicting residual stress and deformation of T-joint welding mentioned in any of the above-mentioned embodiments.

[0113] The present invention also provides an electronic device, which includes: one or more processors; a memory for storing one or more programs, and when the one or more programs are executed by one or more processors, the one or more processors implement the above-mentioned T-joint welding residual stress and deformation prediction method.

[0114] In the present application, unless otherwise clearly specified and limited, a first feature being “above” or “below” a second feature may mean that the first and second features are in direct contact, or the first and second features are in indirect contact through an intermediate medium. Moreover, a first feature being “above”, “above”, and “above” a second feature may mean that the first feature is directly above or obliquely above the second feature, or simply means that the first feature is higher in level than the second feature. A first feature being “below”, “below”, and “below” a second feature may mean that the first feature is directly below or obliquely below the second feature, or simply means that the first feature is lower in level than the second feature.

[0115] In the description of this specification, the description with reference to the terms "preferred embodiment", "further embodiment", "other embodiments" or "specific example" etc. means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, unless they are contradictory.

[0116] Although the embodiments of the present application have been shown and described above, it can be understood that the above embodiments are exemplary and cannot be understood as limitations on the present application. Ordinary technicians in this field can change, modify, replace and modify the above embodiments within the scope of the present application.

Claims

1. A method for predicting residual stress and deformation of T-joint welding, characterized in that: The steps include: S1: Establish a three-dimensional model of steel plate T-joint welding in ABAQUS, assign material properties, mesh the three-dimensional model, define the welding finite element simulation direction, and set the welding step, cooling step and boundary conditions; S2: Write the double-wire MAG welding program language through the Dlux program in Visual Studio software. The specific steps are as follows: S21: setting welding parameters and double ellipsoid molten pool size, and obtaining a heat source equation according to the welding parameters and the double ellipsoid molten pool size, wherein the double ellipsoid molten pool size includes the front half zone size, the rear half zone size, the molten depth and the molten width, and the welding parameters include welding current, welding voltage, welding power coefficient, welding speed, welding starting point and welding direction; S22: rotating the double ellipsoid molten pool by a certain angle so that the direction of the double ellipsoid molten pool is perpendicular to the weld surface; S23: Define the welding starting point position, obtain the movement trajectory equation according to the welding starting point position and the angle of rotation of the double ellipsoid molten pool in step S3, and obtain the heat source movement equation according to the heat source equation and the movement trajectory equation. The heat source movement equation satisfies the following conditions: FLUX(1)=heat1*shape1; FLUX(2)=heat2*shape2; Where FLUX(1) is the heat flux of the first half, FLUX(2) is the heat flux of the second half, heat1 is the heat input source of the first half, heat2 is the heat input source of the second half, shape1 is the moving trajectory of the first half, and shape2 is the moving trajectory of the second half; The heat source equation satisfies the following conditions: heat1=6.0*sqrt(3.0)*q / (a*b*c*PI*sqrt(PI))*f; heat2=6.0*sqrt(3.0)*q / (aa*b*c*PI*sqrt(PI))*(2.0-f); Wherein, q is the heat input size, a is the size of the front half of the heat source, b is the width of the heat source, c is the depth of the heat source, PI is the coordinate information of the welding three-dimensional model, and f is the heat proportion of the front hemisphere; S3: Submit the coordinate information of the meshed three-dimensional model of the T-joint welding in the ABAQUS, as well as the set material properties, welding steps, cooling steps and boundary conditions to the Visual Studio software, activate the Visual Studio software and then run the Dflux language controlled by the Inter Fortran program to perform welding finite element analysis.

2. The method for predicting residual stress and deformation of T-joint welding according to claim 1, characterized in that: The q satisfies the following conditions: q = u*i*effi; Wherein, u is the welding voltage, i is the welding current, and effi is the welding efficiency.

3. The method for predicting residual stress and deformation of T-joint welding according to claim 1, characterized in that: The moving trajectory equation satisfies the following conditions: shape1=exp(-3.0*((x-x0)*cos(betahudu1)+(y-y0)*sin(betahudu1))*2 / a*2-3.0*( (y-y0)*cos(betahudu1)-(x-x0)*sin(betahudu1))*2 / b*2-3.0*((z-z0+d))*2 / c*2); shape2=exp(-3.0*((x-x10)*cos(betahudu2)+(y-y10)*sin(betahudu2))*2 / aa*2-3.0* ((y-y10)*cos(betahudu2)-(x-x0)*sin(betahudu2))*2 / b*2-3.0*((z-z0+d))*2 / c*2); Wherein, x, y, z are coordinate information, betahudu1 is the first half of the arc, and betahudu2 is the second half of the arc; 4. The method for predicting residual stress and deformation of T-joint welding according to claim 3, characterized in that: The betahudu1 and betahudu2 meet the following conditions: betahudu1=beta1 / 180.0*PI; betahudu2=beta2 / 180.0*PI; Among them, beta1 and beta2 are the rotation angles of the double ellipsoid melt pool; 5. The method for predicting residual stress and deformation of T-joint welding according to claim 1, characterized in that: The mesh division is a sweeping mesh division, the mesh density in the weld area and the heat affected zone of the T-joint welding three-dimensional model is greater than the mesh density in other areas, and the mesh is set as a temperature-displacement coupling unit reduced integration.

6. The method for predicting residual stress and deformation of T-joint welding according to claim 1, characterized in that: The material properties include yield stress, expansion coefficient, conductivity and elastic modulus corresponding to different temperatures; the boundary conditions include setting ambient temperature, convection heat transfer coefficient and radiation coefficient; The welding step is to simulate the welding time used in the actual welding process; the cooling step is to simulate the cooling time used in the actual welding process.

7. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the program is executed by a processor, the method for predicting residual stress and deformation of T-joint welding as described in any one of claims 1 to 6 is implemented.

Images