Calculation method and device for thermal self-stress of steel-concrete composite beams

Abstract

The present invention provides a method and device for calculating the temperature self-stress of a steel-concrete composite beam. The method comprises: performing plane grid division on the beam cross section of the steel-concrete composite beam, determining the cross section coordinate system, and obtaining the coordinates of each grid point; determining the temperature corresponding to each grid point; and calculating the temperature based on the coordinates and corresponding temperatures of all grid points and the elastic modulus E of concrete based on the condition that the internal force of the beam cross section is zero in a self-equilibrium state. c , steel elastic modulus E s , concrete linear expansion coefficient p c and the steel linear expansion coefficient p s , the axial strain ε of the beam section is calculated c , curvature around the y axis #imgabs0# and curvature around the z axis #imgabs1# according to the grid coordinates (y, z), t y,z 、E c 、E s 、p c 、p s , ε c , #imgabs2# and #imgabs3# calculate the temperature self-stress σ corresponding to the grid coordinate (y, z) t Through the present invention, the temperature self-stress of any steel-concrete composite beam structure under any temperature field can be accurately calculated.

Description

Technical Field

[0001] The present invention relates to the technical field of structural analysis and calculation, and in particular to a method and device for calculating the temperature self-stress of a steel-concrete composite beam. Background Art

[0002] Steel-concrete composite beams are widely used in bridge construction because they fully utilize the mechanical properties of steel and concrete materials, have the advantages of improving material utilization, saving materials and low cost.

[0003] Steel and concrete have different coefficients of thermal expansion. When the external temperature changes, the heat conduction rates of steel and concrete differ, resulting in a temperature difference between the steel beam and the concrete slab. Therefore, due to these two factors, thermal self-stress will be generated in the steel-concrete composite beam.

[0004] However, currently used methods for calculating thermal self-stress in composite beams assume that the temperature field exhibits a linear distribution of temperature differences along the height of the beam section, primarily encompassing the overall temperature, vertical temperature differences along the height of the beam section, and gradient temperatures. However, in reality, the temperature field distribution in composite beams also exhibits transverse and nonlinear temperature variations, such as transverse temperature differences and exponential temperature differences along the beam section, which cannot be accurately accounted for by traditional thermal self-stress calculation methods. For composite beams with a transverse slope or a non-horizontal top surface, the temperature curve changes perpendicular to the transverse slope. Traditional thermal effects calculations for composite beams are performed along the vertical height of the beam section, making it impossible to accurately calculate these effects. Manual equivalent substitution is required to convert the beam section into a flat-slope section in order to calculate the thermal self-stress of the composite beam. This cumbersome equivalent substitution process does not accurately calculate the thermal self-stress of the composite beam. Summary of the Invention

[0005] The main purpose of the present invention is to provide a method and device for calculating the temperature self-stress of steel-concrete composite beams, aiming to solve the technical problems of incompleteness and inaccuracy of the temperature self-stress calculation methods of steel-concrete composite beams in the prior art.

[0006] In a first aspect, the present invention provides a method for calculating the thermal self-stress of a steel-concrete composite beam, the method comprising:

[0007] Perform plane meshing on the cross section of the steel-concrete composite beam, determine the cross section coordinate system, and obtain the cross section coordinates of each grid point;

[0008] Determine the temperature corresponding to each grid point. The temperature corresponding to the beam section coordinate (y, z) is recorded as t y，z ;

[0009] Based on the condition that the internal force of the beam section is zero in the self-equilibrium state, according to the beam section coordinates and corresponding temperatures of all grid points, and the concrete elastic modulus E c , steel elastic modulus E s , concrete linear expansion coefficient p c and the steel linear expansion coefficient p s , the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis

[0010] According to the beam section coordinates (y, z), t y，z 、E c 、E s 、p c 、p s , ε c 、 and Calculate the temperature self-stress σ corresponding to the beam section coordinates (y, z) t .

[0011] Optionally, the beam cross-section internal force in the self-equilibrium state is zero, according to the beam cross-section coordinates and corresponding temperatures of all grid points, and the concrete elastic modulus E c , steel elastic modulus E s , concrete linear expansion coefficient p c and the steel linear expansion coefficient p s , the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The steps include:

[0012] Based on the first set of equations, the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The first set of equations is:

[0013]

[0014]

[0015]

[0016] Where (y, z) is the coordinate of the beam section, t y，z is the temperature corresponding to the beam section coordinate (y, z), E and p are the elastic modulus and linear expansion coefficient of the material corresponding to the area where the beam section coordinate (y, z) is located.

[0017] Optionally, the beam cross-section internal force in the self-equilibrium state is zero, according to the beam cross-section coordinates and corresponding temperatures of all grid points, and the concrete elastic modulus E c , steel elastic modulus E s , concrete linear expansion coefficient p c and the steel linear expansion coefficient p s , the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The steps include:

[0018] According to the beam section coordinates of all grid points, the equivalent area A of the beam section is obtained by Gaussian integration x , Equivalent static moment S around the y-axis y , Equivalent static moment S around the z axis z , equivalent moment of inertia around the y-axis I y , equivalent moment of inertia around the z-axis I z and product of inertia I yz ;

[0019] Based on the second set of equations, the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The second set of equations is:

[0020]

[0021]

[0022]

[0023] Where (y, z) is the coordinate of the beam section, t y，z is the temperature corresponding to the beam section coordinate (y, z), A c is the concrete area, A s is the concrete area, A x 、S y 、S z , I y , I z and I yz is the equivalent area of ​​the beam section, the equivalent static moment around the y-axis, the equivalent static moment around the z-axis, the equivalent moment of inertia around the y-axis, the equivalent moment of inertia around the z-axis and the product of inertia, E b is the elastic modulus of concrete E c and steel elastic modulus E s The ratio, p c is the linear expansion coefficient of concrete, p s is the linear expansion coefficient of steel.

[0024] Optionally, the origin of the cross-sectional coordinate system is the centroid of the composite beam cross section;

[0025] The above-mentioned condition based on the zero internal force of the beam section in the self-equilibrium state is based on the beam section coordinates and corresponding temperatures of all grid points, as well as the concrete elastic modulus E c , steel elastic modulus E s , concrete linear expansion coefficient p c and the steel linear expansion coefficient p s , the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The steps include:

[0026] According to the beam section coordinates of all grid points, the equivalent area A of the beam section is obtained by Gaussian integration x , equivalent moment of inertia around the y-axis I y , equivalent moment of inertia around the z-axis I z and product of inertia I yz ;

[0027] Based on the third set of equations, the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The third program group is:

[0028]

[0029]

[0030]

[0031] Where (y, z) is the coordinate of the beam section, t y，z is the temperature corresponding to the beam section coordinate (y, z), A c is the concrete area, A s is the concrete area, A x , I y , I z and I yz is the equivalent area of ​​the beam section, the equivalent static moment around the y-axis, the equivalent moment of inertia around the z-axis, and the product of inertia, E b is the elastic modulus of concrete E c and steel elastic modulus E s The ratio, p c is the linear expansion coefficient of concrete, p s is the linear expansion coefficient of steel.

[0032] Optionally, the beam cross-section coordinates (y, z), t y，z 、Ec 、E s 、p c 、p s , ε c 、 and Calculate the temperature self-stress σ corresponding to the beam section coordinates (y, z) t The steps include:

[0033] Based on the self-stress formula, the temperature self-stress σ corresponding to the beam section coordinates (y, z) is calculated t , the self-stress formula is:

[0034]

[0035] Where (y, z) is the coordinate of the beam section, t y，z is the temperature corresponding to the beam section coordinate (y, z), E and p are the elastic modulus and linear expansion coefficient of the material corresponding to the area where the beam section grid coordinate (y, z) is located, ε c 、 and are the axial strain, bending curvature around the y-axis, and bending curvature around the z-axis of the beam section.

[0036] In a second aspect, the present invention further provides a device for calculating the temperature self-stress of a steel-concrete composite beam, the device comprising:

[0037] The meshing module is used to perform plane meshing on the cross section of the steel-concrete composite beam, determine the cross section coordinate system, and obtain the mesh coordinates of each mesh point;

[0038] The temperature determination module is used to determine the temperature corresponding to each grid point. The temperature corresponding to the grid coordinate (y, z) is recorded as t y，z ;

[0039] The first calculation module is used to calculate the concrete elastic modulus E based on the grid coordinates and corresponding temperatures of all grid points under the condition that the internal force of the beam section is zero in the self-equilibrium state. c , steel elastic modulus E s , concrete linear expansion coefficient p c and the steel linear expansion coefficient p s , the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis

[0040] The second calculation module is used to calculate the grid coordinates (y, z), t y，z 、E c 、E s 、p c、p s , ε c 、 and Calculate the temperature self-stress σ corresponding to the grid coordinate (y, z) t .

[0041] Optionally, the first calculation module is used to:

[0042] Based on the first set of equations, the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The first set of equations is:

[0043]

[0044]

[0045]

[0046] Where (y, z) is the grid coordinate, t y，z is the temperature corresponding to the grid coordinate (y, z), E and p are the elastic modulus and linear expansion coefficient of the material corresponding to the area where the grid coordinate (y, z) is located.

[0047] Optionally, the first calculation module is used to:

[0048] According to the grid coordinates of all grid points, the equivalent area A of the beam section is obtained by Gaussian integration x , Equivalent static moment S around the y-axis y , Equivalent static moment S around the z axis z , equivalent moment of inertia around the y-axis I y , equivalent moment of inertia around the z-axis I z and product of inertia I yz ;

[0049] Based on the second set of equations, the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The second set of equations is:

[0050]

[0051]

[0052]

[0053] Where (y, z) is the grid coordinate, t y，z is the temperature corresponding to the grid coordinate (y, z), A cis the concrete area, A s is the concrete area, A x 、S y 、S z , I y , I z and Iy z is the equivalent area of ​​the beam section, the equivalent static moment around the y-axis, the equivalent static moment around the z-axis, the equivalent moment of inertia around the y-axis, the equivalent moment of inertia around the z-axis and the product of inertia, E b is the elastic modulus of concrete E c and steel elastic modulus E s The ratio, p c is the linear expansion coefficient of concrete, p s is the linear expansion coefficient of steel.

[0054] Optionally, the origin of the cross-section coordinate system in the meshing module is the centroid of the composite beam cross section, and the first calculation module is used to:

[0055] According to the grid coordinates of all grid points, the equivalent area A of the beam section is obtained by Gaussian integration x , equivalent moment of inertia around the y-axis I y , equivalent moment of inertia around the z-axis I z and product of inertia I yz ;

[0056] Based on the third set of equations, the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The third program group is:

[0057]

[0058]

[0059]

[0060] Where (y, z) is the grid coordinate, t y，z is the temperature corresponding to the grid coordinate (y, z), A c is the concrete area, A s is the concrete area, A x , I y , I z and I yz is the equivalent area of ​​the beam section, the equivalent static moment around the y-axis, the equivalent moment of inertia around the z-axis, and the product of inertia, E b is the elastic modulus of concrete E c and steel elastic modulus E s The ratio, p c is the linear expansion coefficient of concrete, ps is the linear expansion coefficient of steel.

[0061] Optionally, the second calculation module is used to:

[0062] Based on the self-stress formula, the temperature self-stress σ corresponding to the grid coordinate (y, z) is calculated t , the self-stress formula is:

[0063]

[0064] Where (y, z) is the grid coordinate, t y，z is the temperature corresponding to the grid coordinate (y, z), E and p are the elastic modulus and linear expansion coefficient of the material corresponding to the area where the grid coordinate (y, z) is located, ε c 、 and are the axial strain, bending curvature around the y-axis, and bending curvature around the z-axis of the beam section.

[0065] In the present invention, the structure of the steel-concrete composite beam and the temperature field distribution of the beam cross section do not affect the applicability and accuracy of the calculation method. Through the present invention, the temperature self-stress of any steel-concrete composite beam structure under any temperature field can be accurately calculated. BRIEF DESCRIPTION OF THE DRAWINGS

[0066] Figure 1 Schematic diagram of a flow chart of a method for calculating the temperature self-stress of a steel-concrete composite beam in one embodiment of the present invention;

[0067] Figure 2 Schematic diagram of the cross section of a steel-concrete composite beam in one embodiment of the present invention;

[0068] Figure 3 for Figure 2 Meshing diagram of the beam section shown;

[0069] Figure 4 Schematic diagram of the hardware structure of a device for calculating the temperature self-stress of a steel-concrete composite beam in one embodiment of the present invention.

[0070] The purpose, features and advantages of the present invention will be further described with reference to the accompanying drawings and in conjunction with the embodiments. DETAILED DESCRIPTION

[0071] It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.

[0072] In a first aspect, an embodiment of the present invention provides a method for calculating the temperature self-stress of a steel-concrete composite beam.

[0073] Figure 1A schematic flow chart of a method for calculating the temperature self-stress of a steel-concrete composite beam in one embodiment of the present invention is shown.

[0074] Reference Figure 1 In one embodiment, a method for calculating the thermal self-stress of a steel-concrete composite beam includes the following steps:

[0075] S11. Perform plane meshing on the cross section of the steel-concrete composite beam, determine the cross section coordinate system, and obtain the grid coordinates of each grid point;

[0076] In this embodiment, the beam cross section is perpendicular to the axis of the steel-concrete composite beam; that is, the beam cross section is the cross section of the steel-concrete composite beam. After plane meshing, multiple grid points are determined on the beam cross section. After selecting a cross-sectional coordinate system as needed, the grid coordinates of each grid point are obtained. It will be appreciated that plane meshing is applicable to any composite beam cross section.

[0077] Figure 2 A schematic structural diagram of a cross-section of a steel-concrete composite beam in one embodiment of the present invention is shown; Figure 3 Shown Figure 2 Meshing diagram of the beam section shown.

[0078] Reference Figure 2 and Figure 3 Specifically, in a steel-concrete composite beam, a spatial rectangular coordinate system (xyz) is established with the axial direction as the x-axis, the transverse direction as the y-axis, and the height direction as the z-axis. Correspondingly, the cross-sectional coordinate system of the beam section is the plane rectangular coordinate system (yz).

[0079] As an optional implementation method, based on the geometric characteristics of the beam section of the steel-concrete composite beam, the finite element analysis software ABQUS is used to divide the beam section into plane meshes. The ABQUS inp file is exported through automesh to determine the section coordinate system and obtain the grid coordinates of each grid point.

[0080] S12. Determine the temperature corresponding to each grid point. The temperature corresponding to the grid coordinate (y, z) is recorded as t y,z ;

[0081] In this embodiment, there are no restrictions on the temperature field distribution form of the beam cross section; the temperature corresponding to the grid point can be determined under any temperature field. For example, the temperature change curve can change in the vertical direction, the horizontal direction, the direction perpendicular to the cross slope line, etc., and the change pattern can be linear or nonlinear (exponential distribution, gradient distribution, etc.). During use, the temperature of each grid point on the beam cross section is defined according to the preset temperature field distribution form. Therefore, this embodiment is applicable to any combination of beam cross sections under any temperature field.

[0082] Specific to Figure 2 、 Figure 3 In the embodiment shown, there is a transverse slope on the top of the steel-concrete composite beam. The temperature change curve changes in a direction perpendicular to the transverse slope line, and the change law is an exponential distribution, which is expressed as T(d)=T0·e -d , where T(d) is the temperature of the target point, T0 is the temperature of the top surface of the horizontal slope, and d is the vertical distance from the target point to the top surface of the horizontal slope. In this embodiment, the target point is a grid point.

[0083] S13. Based on the condition that the internal force of the beam section is zero in the self-equilibrium state, according to the grid coordinates and corresponding temperatures of all grid points, and the elastic modulus E of concrete c , steel elastic modulus E s , concrete linear expansion coefficient p c and the steel linear expansion coefficient p s , the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis

[0084] In this embodiment, the internal forces of the beam section include the axial force N along the x-axis, the bending moment M around the y-axis, and the y and the bending moment M about the z-axis z , N, M y and M z The grid coordinates and corresponding temperatures of all grid points, as well as the elastic modulus of concrete E c , steel elastic modulus E s , concrete linear expansion coefficient p c , steel linear expansion coefficient p s , axial strain ε of the beam section c , curvature around the y-axis and the curvature around the z-axis Integral representation. In N, M y and M z Under the condition that both are zero, we can solve the system of equations and get ε c 、 and For example, the elastic modulus of concrete is 34500 MPa and the linear expansion coefficient of concrete is 1.0e -5 , the elastic modulus of steel is 210000MPa, and the linear expansion coefficient of steel is 1.2e -5 .

[0085] Specifically, N, M y and M z The expression is:

[0086]

[0087]

[0088]

[0089] Where (y, z) is the grid coordinate, t y，z is the temperature corresponding to the grid coordinate (y, z), E and p are the elastic modulus and linear expansion coefficient of the material in the area where the grid coordinate (y, z) is located. If the grid coordinate (y, z) is in the concrete area, then E = E c 、p=p c , if the grid coordinate (y, z) is in the steel area, then E=E s 、p=p s .

[0090] In the first embodiment, step S13 specifically includes:

[0091] Based on the first set of equations, the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The first set of equations is:

[0092]

[0093]

[0094]

[0095] Where (y, z) is the grid coordinate, t y，z is the temperature corresponding to the grid coordinate (y, z), E and p are the elastic modulus and linear expansion coefficient of the material corresponding to the area where the grid coordinate (y, z) is located.

[0096] In this embodiment, the first set of equations is based on N, M y and M z The expression of N, M y and M z The condition of both being zero is established directly.

[0097] In the second implementation, step S13 specifically includes:

[0098] According to the grid coordinates of all grid points, the equivalent area A of the beam section is obtained by Gaussian integration x , Equivalent static moment S around the y-axis y , Equivalent static moment S around the z axis z , equivalent moment of inertia around the y-axis I y , equivalent moment of inertia around the z-axis I z and product of inertia I yz ;

[0099] Based on the second set of equations, the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The second set of equations is:

[0100]

[0101]

[0102]

[0103] Where (y, z) is the grid coordinate, t y，z is the temperature corresponding to the grid coordinate (y, z), A c is the concrete area, A s is the concrete area, A x 、S y 、S z , I y , I z and I yz is the equivalent area of ​​the beam section, the equivalent static moment around the y-axis, the equivalent static moment around the z-axis, the equivalent moment of inertia around the y-axis, the equivalent moment of inertia around the z-axis and the product of inertia, E b is the elastic modulus of concrete E c and steel elastic modulus E s The ratio, p c is the linear expansion coefficient of concrete, p s is the linear expansion coefficient of steel.

[0104] In this embodiment, based on the first set of equations, the intermediate set of equations is obtained by expressing them according to different materials:

[0105]

[0106]

[0107]

[0108] Where (y, z) is the grid coordinate, t y，z is the temperature corresponding to the grid coordinate (y, z), A c is the concrete area, A s For concrete area, E c is the elastic modulus of concrete, E s is the elastic modulus of steel, p c is the linear expansion coefficient of concrete, p s is the linear expansion coefficient of steel.

[0109] By classifying and summarizing the intermediate equations, we can obtain the second equations in this embodiment. First, we obtain A by Gaussian integral. x 、S y 、S z , I y , I z and I yz , and then use these values ​​as known quantities, together with the grid coordinates and corresponding temperatures of all grid points, as well as the elastic modulus and linear expansion coefficient of concrete and steel, to solve the second set of equations. This helps to simplify the solution process, reduce the amount of calculations, and improve calculation efficiency.

[0110] In the third embodiment, the origin of the cross-section coordinate system is the centroid of the composite beam cross section, and step S13 specifically includes:

[0111] According to the grid coordinates of all grid points, the equivalent area A of the beam section is obtained by Gaussian integration x , equivalent moment of inertia around the y-axis I y , equivalent moment of inertia around the z-axis I z and product of inertia I yz ;

[0112] Based on the third set of equations, the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The third set of equations is:

[0113]

[0114]

[0115]

[0116] Where (y, z) is the grid coordinate, t y，z is the temperature corresponding to the grid coordinate (y, z), A c is the concrete area, A s is the concrete area, A x , I y , I z and I yz is the equivalent area of ​​the beam section, the equivalent static moment around the y-axis, the equivalent moment of inertia around the z-axis, and the product of inertia, E b is the elastic modulus of concrete E c and steel elastic modulus E s The ratio, p c is the linear expansion coefficient of concrete, p s is the linear expansion coefficient of steel.

[0117] In this embodiment, when determining the cross-sectional coordinate system in step S11, the centroid of the composite beam cross section is used as the origin of the cross-sectional coordinate system. y 、S z are all 0, which helps to further simplify the solution process, reduce the amount of calculation and improve the calculation efficiency.

[0118] S14, according to the grid coordinates (y, z), t y，z 、E c 、E s 、p c 、p s , ε c 、 and Calculate the temperature self-stress σ corresponding to the grid coordinate (y, z) t .

[0119] In this embodiment, for the target grid point, the grid coordinates (y, z), t y，z 、E c 、E s 、p c 、p s , ε c 、 and The temperature self-stress σt at the target grid point can be calculated.

[0120] Specifically, based on the self-stress formula, the temperature self-stress σ corresponding to the grid coordinate (y, z) is calculated as t , the self-stress formula is:

[0121]

[0122] Where (y, z) is the grid coordinate, t y，z is the temperature corresponding to the grid coordinate (y, z), E and p are the elastic modulus and linear expansion coefficient of the material corresponding to the area where the grid coordinate (y, z) is located, ε c 、 and are the axial strain, bending curvature around the y-axis, and bending curvature around the z-axis of the beam section.

[0123] Therefore, in this embodiment, the structure of the steel-concrete composite beam and the temperature field distribution of the beam cross section will not affect the applicability and accuracy of the calculation method. Through this embodiment, the thermal self-stress of any steel-concrete composite beam structure under any temperature field can be accurately calculated.

[0124] In a second aspect, an embodiment of the present invention further provides a device for calculating the temperature self-stress of a steel-concrete composite beam.

[0125] Figure 4The figure shows a hardware structure diagram of a device for calculating the temperature self-stress of a steel-concrete composite beam in one embodiment of the present invention.

[0126] Reference Figure 4 In one embodiment, the device for calculating the temperature self-stress of a steel-concrete composite beam includes:

[0127] A meshing module 10 is used to perform plane meshing on the cross section of the steel-concrete composite beam, determine the cross section coordinate system, and obtain the mesh coordinates of each mesh point;

[0128] The temperature determination module 20 is used to determine the temperature corresponding to each grid point. The temperature corresponding to the grid coordinate (y, z) is recorded as t y，z ;

[0129] The first calculation module 30 is used to calculate the internal force of the beam section under the self-equilibrium state according to the grid coordinates and corresponding temperatures of all grid points and the elastic modulus E of concrete. c , steel elastic modulus E s , concrete linear expansion coefficient p c and the steel linear expansion coefficient p s , the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis

[0130] The second calculation module 40 is used to calculate the grid coordinates (y, z), t y，z 、E c 、E s 、p c 、p s , ε c 、 and Calculate the temperature self-stress σ corresponding to the grid coordinate (y, z) t .

[0131] Furthermore, in one embodiment, the first calculation module 30 is configured to:

[0132] Based on the first set of equations, the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The first set of equations is:

[0133]

[0134]

[0135]

[0136] Where (y, z) is the grid coordinate, t y，z is the temperature corresponding to the grid coordinate (y, z), E and p are the elastic modulus and linear expansion coefficient of the material corresponding to the area where the grid coordinate (y, z) is located.

[0137] Furthermore, in one embodiment, the first calculation module 30 is configured to:

[0138] According to the grid coordinates of all grid points, the equivalent area A of the beam section is obtained by Gaussian integration x , Equivalent static moment S around the y-axis y , Equivalent static moment S around the z axis z , equivalent moment of inertia around the y-axis I y , equivalent moment of inertia around the z-axis I z and product of inertia I yz ;

[0139] Based on the second set of equations, the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The second set of equations is:

[0140]

[0141]

[0142]

[0143] Where (y, z) is the grid coordinate, t y，z is the temperature corresponding to the grid coordinate (y, z), A c is the concrete area, A s is the concrete area, A x 、S y 、S z , I y , I z and I yz is the equivalent area of ​​the beam section, the equivalent static moment around the y-axis, the equivalent static moment around the z-axis, the equivalent moment of inertia around the y-axis, the equivalent moment of inertia around the z-axis and the product of inertia, E b is the elastic modulus of concrete E c and steel elastic modulus E s The ratio, p c is the linear expansion coefficient of concrete, p s is the linear expansion coefficient of steel.

[0144] Furthermore, in one embodiment, the origin of the cross-section coordinate system in the meshing module 10 is the centroid of the composite beam cross section, and the first calculation module 30 is used to:

[0145] According to the grid coordinates of all grid points, the equivalent area A of the beam section is obtained by Gaussian integration x , equivalent moment of inertia around the y-axis I y , equivalent moment of inertia around the z-axis I z and product of inertia I yz ;

[0146] Based on the third set of equations, the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The third set of equations is:

[0147]

[0148]

[0149]

[0150] Where (y, z) is the grid coordinate, t y，z is the temperature corresponding to the grid coordinate (y, z), A c is the concrete area, A s is the concrete area, A x , I y , I z and I yz is the equivalent area of ​​the beam section, the equivalent static moment around the y-axis, the equivalent moment of inertia around the z-axis, and the product of inertia, E b is the elastic modulus of concrete E c and steel elastic modulus E s The ratio, p c is the linear expansion coefficient of concrete, p s is the linear expansion coefficient of steel.

[0151] Furthermore, in one embodiment, the second calculation module 40 is configured to:

[0152] Based on the self-stress formula, the temperature self-stress σ corresponding to the grid coordinate (y, z) is calculated t , the self-stress formula is:

[0153]

[0154] Where (y, z) is the grid coordinate, t y，z is the temperature corresponding to the grid coordinate (y, z), E and p are the elastic modulus and linear expansion coefficient of the material corresponding to the area where the grid coordinate (y, z) is located, ε c 、 and are the axial strain, bending curvature around the y-axis, and bending curvature around the z-axis of the beam section.

[0155] Among them, the functional implementation of each module in the above-mentioned steel-concrete composite beam temperature self-stress calculation device corresponds to the various steps in the above-mentioned steel-concrete composite beam temperature self-stress calculation method embodiment, and its functions and implementation processes are no longer repeated here.

[0156] It should be noted that, in this document, the terms "comprises," "includes," or any other variations thereof are intended to encompass non-exclusive inclusion, such that a process, method, article, or system comprising a series of elements includes not only those elements but also other elements not explicitly listed, or elements inherent to such process, method, article, or system. In the absence of further limitations, an element defined by the phrase "comprising a ..." does not exclude the presence of other identical elements in the process, method, article, or system comprising the element.

[0157] The serial numbers of the above embodiments of the present invention are for description only and do not represent the advantages or disadvantages of the embodiments.

[0158] The above are only preferred embodiments of the present invention and are not intended to limit the patent scope of the present invention. Any equivalent structure or equivalent process transformation made using the contents of the present invention description and drawings, or directly or indirectly applied in other related technical fields, are also included in the patent protection scope of the present invention.

Claims

1. A method for calculating the thermal self-stress of steel-concrete composite beams, characterized in that: The method for calculating the temperature self-stress of the steel-concrete composite beam includes: Perform plane meshing on the cross section of the steel-concrete composite beam, determine the cross section coordinate system, and obtain the grid coordinates of each grid point; Determine the temperature corresponding to each grid point. The temperature corresponding to the grid coordinate (y, z) is recorded as t y,z ; Based on the condition that the internal force of the beam section is zero in the self-equilibrium state, according to the grid coordinates and corresponding temperatures of all grid points, and the elastic modulus E of concrete c , steel elastic modulus E s , concrete linear expansion coefficient p c and the steel linear expansion coefficient p s , the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis According to the grid coordinates (y, z), t y,z 、E c 、E s 、p c 、p s , ε c 、 and Calculate the temperature self-stress σ corresponding to the grid coordinate (y, z) t .

2. The method for calculating the thermal self-stress of steel-concrete composite beams according to claim 1, wherein: The condition that the internal force of the beam section is zero in the self-equilibrium state is based on the grid coordinates of all grid points and the corresponding temperature, as well as the elastic modulus E of the concrete. c , steel elastic modulus E s , concrete linear expansion coefficient p c and the steel linear expansion coefficient p s , the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The steps include: Based on the first set of equations, the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The first set of equations is: Where (y, z) is the grid coordinate, t y,z is the temperature corresponding to the grid coordinate (y, z), E and p are the elastic modulus and linear expansion coefficient of the material corresponding to the area where the grid coordinate (y, z) is located.

3. The method for calculating the thermal self-stress of steel-concrete composite beams according to claim 1, wherein: The above-mentioned condition based on the zero internal force of the beam section in the self-equilibrium state is based on the coordinates of all grid points and the corresponding temperatures, as well as the elastic modulus E of the concrete. c , steel elastic modulus E s , concrete linear expansion coefficient p c and the steel linear expansion coefficient p s , the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The steps include: According to the coordinates of all grid points, the equivalent area A of the beam section is obtained by Gaussian integration x , Equivalent static moment S around the y-axis y , Equivalent static moment S around the z axis z , equivalent moment of inertia around the y-axis I y , equivalent moment of inertia around the z-axis I z and product of inertia I yz ; Based on the second set of equations, the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The second set of equations is: Where (y, z) is the grid coordinate, t y,z is the temperature corresponding to the grid coordinate (y, z), A c is the concrete area, A s is the concrete area, A x 、S y 、S z , I y , I z and I yz is the equivalent area of ​​the beam section, the equivalent static moment around the y-axis, the equivalent static moment around the z-axis, the equivalent moment of inertia around the y-axis, the equivalent moment of inertia around the z-axis and the product of inertia, E b is the elastic modulus of concrete E c and steel elastic modulus E s The ratio, p c is the linear expansion coefficient of concrete, p s is the linear expansion coefficient of steel.

4. The method for calculating the thermal self-stress of steel-concrete composite beams according to claim 1, wherein: The origin of the cross-sectional coordinate system is the centroid of the composite beam cross section; The condition that the internal force of the beam section is zero in the self-equilibrium state is based on the grid coordinates of all grid points and the corresponding temperature, as well as the elastic modulus E of the concrete. c , steel elastic modulus E s , concrete linear expansion coefficient p c and the steel linear expansion coefficient p s , the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The steps include: According to the grid coordinates of all grid points, the equivalent area A of the beam section is obtained by Gaussian integration x , equivalent moment of inertia around the y-axis I y , equivalent moment of inertia around the z-axis I z and product of inertia I yz ; Based on the third set of equations, the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The third program group is: Where (y, z) is the grid coordinate, t y,z is the temperature corresponding to the grid coordinate (y, z), A c is the concrete area, A s is the concrete area, A x , I y , I z and I yz is the equivalent area of ​​the beam section, the equivalent static moment around the y-axis, the equivalent moment of inertia around the z-axis, and the product of inertia, E b is the elastic modulus of concrete E c and steel elastic modulus E s The ratio, p c is the linear expansion coefficient of concrete, p s is the linear expansion coefficient of steel.

5. The method for calculating the thermal self-stress of steel-concrete composite beams according to claim 1, wherein: According to the beam section coordinates (y, z), t y,z 、E c 、E s 、p c 、p s , ε c 、 and Calculate the temperature self-stress σ corresponding to the beam section coordinates (y, z) t The steps include: Based on the self-stress formula, the temperature self-stress σ corresponding to the beam section coordinates (y, z) is calculated t , the self-stress formula is: Where (y, z) is the coordinate of the beam section, t y,z is the temperature corresponding to the beam section coordinate (y, z), E and p are the elastic modulus and linear expansion coefficient of the material in the area where the beam section coordinate (y, z) is located, ε c 、 and are the axial strain, bending curvature around the y-axis, and bending curvature around the z-axis of the beam section.

6. A device for calculating the temperature self-stress of steel-concrete composite beams, characterized in that: The device for calculating the temperature self-stress of the steel-concrete composite beam comprises: The meshing module is used to perform plane meshing on the cross section of the steel-concrete composite beam, determine the cross-sectional coordinate system, and obtain the cross-sectional coordinates of each grid point; The temperature determination module is used to determine the temperature corresponding to each grid point. The temperature corresponding to the beam section coordinate (y, z) is recorded as t y,z ; The first calculation module is used to calculate the beam cross-section coordinates and corresponding temperatures of all grid points based on the condition that the internal force of the beam cross-section is zero in the self-equilibrium state, as well as the concrete elastic modulus E c , steel elastic modulus E s , concrete linear expansion coefficient p c and the steel linear expansion coefficient p s , the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The second calculation module is used to calculate the beam section coordinates (y, z), t y,z 、E c 、E s 、p c 、p s , ε c 、 and Calculate the temperature self-stress σ corresponding to the beam section coordinates (y, z) t .

7. The device for calculating the temperature self-stress of steel-concrete composite beams according to claim 6, characterized in that: The first calculation module is used for: Based on the first set of equations, the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The first set of equations is: Where (y, z) is the coordinate of the beam section, t y,z is the temperature corresponding to the beam section coordinate (y, z), E and p are the elastic modulus and linear expansion coefficient of the material corresponding to the area where the beam section coordinate (y, z) is located.

8. The device for calculating the temperature self-stress of steel-concrete composite beams according to claim 6, characterized in that: The first calculation module is used for: According to the beam section coordinates of all grid points, the equivalent area A of the beam section is obtained by Gaussian integration x , Equivalent static moment S around the y-axis y , Equivalent static moment S around the z axis z , equivalent moment of inertia around the y-axis I y , equivalent moment of inertia around the z-axis I z and product of inertia I yz ; Based on the second set of equations, the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The second set of equations is: Where (y, z) is the coordinate of the beam section, t y,z is the temperature corresponding to the beam section coordinate (y, z), A c is the concrete area, A s is the concrete area, A x 、S y 、S z , I y , I z and I yz is the equivalent area of ​​the beam section, the equivalent static moment around the y-axis, the equivalent static moment around the z-axis, the equivalent moment of inertia around the y-axis, the equivalent moment of inertia around the z-axis and the product of inertia, E b is the elastic modulus of concrete E c and steel elastic modulus E s The ratio, p c is the linear expansion coefficient of concrete, p s is the linear expansion coefficient of steel.

9. The device for calculating the temperature self-stress of steel-concrete composite beams according to claim 6, characterized in that: The origin of the cross-section coordinate system in the meshing module is the centroid of the composite beam cross section. The first calculation module is used to: According to the beam section coordinates of all grid points, the equivalent area A of the beam section is obtained by Gaussian integration x , equivalent moment of inertia around the y-axis I y , equivalent moment of inertia around the z-axis I z and product of inertia I yz ; Based on the third set of equations, the axial strain ε of the beam section is calculated c , curvature around the y-axis and the curvature around the z-axis The third program group is: Where (y, z) is the coordinate of the beam section, t y,z is the temperature corresponding to the beam section coordinate (y, z), A c is the concrete area, A s is the concrete area, A x , I y , I z and I yz is the equivalent area of ​​the beam section, the equivalent static moment around the y-axis, the equivalent moment of inertia around the z-axis, and the product of inertia, E b is the elastic modulus of concrete E c and steel elastic modulus E s The ratio, p c is the linear expansion coefficient of concrete, p s is the linear expansion coefficient of steel.

10. The device for calculating thermal self-stress of steel-concrete composite beams according to claim 6, characterized in that: The second calculation module is used for: Based on the self-stress formula, the temperature self-stress σ corresponding to the beam section coordinates (y, z) is calculated t , the self-stress formula is: Where (y, z) is the coordinate of the beam section, t y,z is the temperature corresponding to the beam section coordinate (y, z), E and p are the elastic modulus and linear expansion coefficient of the material in the area where the beam section coordinate (y, z) is located, ε c 、 and are the axial strain, bending curvature around the y-axis, and bending curvature around the z-axis of the beam section.

Images