[go: up one dir, main page]

CN108614913B - Calculation method and device for reinforcement ratio of concrete unit - Google Patents

Calculation method and device for reinforcement ratio of concrete unit Download PDF

Info

Publication number
CN108614913B
CN108614913B CN201810252692.7A CN201810252692A CN108614913B CN 108614913 B CN108614913 B CN 108614913B CN 201810252692 A CN201810252692 A CN 201810252692A CN 108614913 B CN108614913 B CN 108614913B
Authority
CN
China
Prior art keywords
unit
concrete
reinforcement
steel bar
units
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810252692.7A
Other languages
Chinese (zh)
Other versions
CN108614913A (en
Inventor
郭胜山
张爱静
李德玉
廖建新
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Institute of Water Resources and Hydropower Research
China Three Gorges Projects Development Co Ltd CTG
Original Assignee
China Institute of Water Resources and Hydropower Research
China Three Gorges Projects Development Co Ltd CTG
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Institute of Water Resources and Hydropower Research, China Three Gorges Projects Development Co Ltd CTG filed Critical China Institute of Water Resources and Hydropower Research
Priority to CN201810252692.7A priority Critical patent/CN108614913B/en
Publication of CN108614913A publication Critical patent/CN108614913A/en
Application granted granted Critical
Publication of CN108614913B publication Critical patent/CN108614913B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Architecture (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Reinforcement Elements For Buildings (AREA)

Abstract

The invention provides a method and a device for calculating the reinforcement ratio of a concrete unit, wherein the method comprises the following steps: the method comprises the steps of obtaining a concrete model and a steel bar model of a target building, adopting a preset geometric unit to divide the concrete model, obtaining a concrete finite element model, wherein the concrete finite element model comprises M concrete units, M is an integer larger than 0, adopting a preset linear unit to divide the steel bar model, obtaining the steel bar finite element model, the steel bar finite element model comprises N steel bar units, N is an integer larger than 0, and obtaining the reinforcement ratio of each concrete unit according to the position information of the steel bar units and the position information of the concrete units. The method is suitable for calculating the reinforcement ratio of the concrete units in more complex projects, and can provide data support for subsequent finite element calculation and analysis of the combined action of the reinforced concrete.

Description

混凝土单元配筋率的计算方法及装置Calculation method and device for reinforcement ratio of concrete elements

技术领域technical field

本发明涉及工程数值仿真领域,尤其涉及一种混凝土单元配筋率的计算方法及装置。The invention relates to the field of engineering numerical simulation, in particular to a method and a device for calculating the reinforcement ratio of a concrete unit.

背景技术Background technique

重力坝、拱坝等坝体是常见的大型水工建筑物,其中一些分布于我国的西部强震区,当发生地震时,坝体有开裂的风险,因此,需要通过坝面配筋来增加坝体的抗震性。一般地,需要计算工程中建筑物的混凝土的配筋率,通过混凝土的配筋率来反应工程的配筋工况,同时,通过控制混凝土的配筋率还可以控制工程中结构构件的破坏形态,避免超筋破坏和少筋破坏情况的发生。Dam bodies such as gravity dams and arch dams are common large-scale hydraulic structures, some of which are distributed in the strong earthquake area in the west of my country. When an earthquake occurs, the dam body is at risk of cracking. Therefore, it is necessary to increase the dam surface reinforcement. Seismic resistance of the dam body. Generally, it is necessary to calculate the reinforcement ratio of the concrete of the building in the project, and the reinforcement ratio of the concrete can be used to reflect the reinforcement conditions of the project. , to avoid the occurrence of super-strand damage and less-strength damage.

现有技术计算建筑物混凝土的配筋率的方法是通过钢筋体积除以混凝土体积得到混凝土的配筋率。现有技术的方法只适用于梁、柱、板等这类简单建筑物的配筋率的计算,对于工程中复杂的建筑物模型和复杂的配筋工况并不适用,并且通过现有技术获取的混凝土的配筋率无法为后续钢筋混凝土联合作用有限元分析提供数据支持。The method of calculating the reinforcement ratio of building concrete in the prior art is to obtain the reinforcement ratio of concrete by dividing the volume of reinforcement by the volume of concrete. The method in the prior art is only applicable to the calculation of the reinforcement ratio of such simple buildings as beams, columns, slabs, etc., and is not applicable to the complex building models and complex reinforcement conditions in the engineering, and through the existing technology The obtained reinforcement ratio of concrete cannot provide data support for the subsequent finite element analysis of reinforced concrete combined action.

发明内容SUMMARY OF THE INVENTION

针对现有技术中的缺陷,本发明提供一种混凝土单元配筋率的计算方法及装置。Aiming at the defects in the prior art, the present invention provides a method and device for calculating the reinforcement ratio of a concrete unit.

一方面,本发明提供一种混凝土单元配筋率的计算方法,所述方法包括:In one aspect, the present invention provides a method for calculating the reinforcement ratio of a concrete unit, the method comprising:

获取目标建筑物的混凝土模型和钢筋模型;Obtain the concrete model and reinforcement model of the target building;

采用预设立体几何单元对所述混凝土模型进行划分,获取混凝土有限元模型,所述混凝土有限元模型包括M个混凝土单元,M为大于0的整数;Divide the concrete model by using a preset three-dimensional geometric unit to obtain a concrete finite element model, where the concrete finite element model includes M concrete units, where M is an integer greater than 0;

采用预设线型单元对所述钢筋模型进行划分,获取钢筋有限元模型,所述钢筋有限元模型包括N个钢筋单元,N为大于0的整数;Use preset linear elements to divide the steel bar model, and obtain a steel bar finite element model, where the steel bar finite element model includes N steel bar elements, and N is an integer greater than 0;

根据所述钢筋单元的位置信息、所述混凝土单元的位置信息,获取各混凝土单元的配筋率。According to the position information of the reinforcing bar unit and the position information of the concrete unit, the reinforcement ratio of each concrete unit is acquired.

进一步地,根据所述钢筋单元的位置信息、所述混凝土单元的位置信息,获取各混凝土单元的配筋率之前,还包括:Further, before acquiring the reinforcement ratio of each concrete unit according to the location information of the reinforcing bar unit and the location information of the concrete unit, the method further includes:

根据所述混凝土单元的形函数、所述钢筋单元的位置信息、所述混凝土单元的位置信息,建立所述钢筋单元与所述混凝土单元的映射关系。According to the shape function of the concrete unit, the location information of the reinforcing bar unit, and the location information of the concrete unit, a mapping relationship between the reinforcing bar unit and the concrete unit is established.

进一步地,所述根据所述混凝土单元的形函数、所述钢筋单元的位置信息、所述混凝土单元的位置信息,建立所述钢筋单元与所述混凝土单元的映射关系,包括:Further, according to the shape function of the concrete unit, the location information of the reinforcing bar unit, and the location information of the concrete unit, establishing the mapping relationship between the reinforcing bar unit and the concrete unit includes:

根据所述混凝土单元的形函数、所述混凝土单元的位置信息,对所述混凝土单元等参变换,并采用牛顿迭代法计算,计算各所述钢筋单元端点的局部坐标;According to the shape function of the concrete unit and the position information of the concrete unit, isoparametric transformation of the concrete unit is performed, and the Newton iteration method is used for calculation to calculate the local coordinates of the end points of each of the reinforcing bar units;

将所述钢筋单元端点的局部坐标所在的混凝土单元,与所述钢筋单元建立映射关系。A mapping relationship is established between the concrete unit where the local coordinates of the end points of the reinforcing bar unit are located and the reinforcing bar unit.

进一步地,根据所述钢筋单元的位置信息、所述混凝土单元的位置信息,获取各混凝土单元的配筋率,包括:Further, according to the position information of the reinforcement unit and the position information of the concrete unit, the reinforcement ratio of each concrete unit is obtained, including:

根据所述钢筋单元与所述混凝土单元的映射关系,遍历各混凝土单元存在映射关系的钢筋单元;According to the mapping relationship between the reinforcing bar unit and the concrete unit, traverse the reinforcing bar units with the mapping relationship of each concrete unit;

根据各混凝土单元存在映射关系的钢筋单元中,各钢筋单元位于对应混凝土单元内的端点数,获取各混凝土单元的配筋率。The reinforcement ratio of each concrete element is obtained according to the number of end points of each reinforcement element located in the corresponding concrete element among the reinforcement elements in which each concrete element has a mapping relationship.

进一步地,根据各混凝土单元存在映射关系的钢筋单元中,各钢筋单元位于对应混凝土单元内的端点数,获取各混凝土单元的配筋率,包括:Further, obtain the reinforcement ratio of each concrete unit according to the number of end points of each reinforcement unit in the corresponding concrete unit in the reinforcement unit with the mapping relationship of each concrete unit, including:

计算各混凝土单元存在映射关系的钢筋单元在对应混凝土单元内的体积,其中,若钢筋单元位于对应混凝土单元内的端点数为1,则所述钢筋单元在对应混凝土单元内的体积为所述钢筋单元实际体积的1/2;若钢筋单元位于对应混凝土单元内的端点数为2,则所述钢筋单元在对应混凝土单元内的体积为所述钢筋单元实际体积;Calculate the volume of the reinforcement unit with the mapping relationship of each concrete unit in the corresponding concrete unit, wherein, if the number of end points of the reinforcement unit in the corresponding concrete unit is 1, then the volume of the reinforcement unit in the corresponding concrete unit is the reinforcement 1/2 of the actual volume of the unit; if the number of end points of the reinforcement unit in the corresponding concrete unit is 2, the volume of the reinforcement unit in the corresponding concrete unit is the actual volume of the reinforcement unit;

根据各混凝土单元存在映射关系的钢筋单元在对应混凝土单元内的体积,获取所述各混凝土单元的配筋率。The reinforcement ratio of each concrete unit is obtained according to the volume of the reinforcing bar unit with the mapping relationship of each concrete unit in the corresponding concrete unit.

另一方面,本发明还提供一种混凝土单元配筋率的计算装置,所述装置包括:In another aspect, the present invention also provides a device for calculating the reinforcement ratio of a concrete unit, the device comprising:

获取模块,用于获取目标建筑物的混凝土模型和钢筋模型;The acquisition module is used to acquire the concrete model and reinforcement model of the target building;

第一网格划分模块,用于采用预设立体几何单元对所述混凝土模型进行划分,获取混凝土有限元模型,所述混凝土有限元模型包括M个混凝土单元,M为大于0的整数;a first meshing module, configured to divide the concrete model by using preset three-dimensional geometric units to obtain a concrete finite element model, where the concrete finite element model includes M concrete units, where M is an integer greater than 0;

第二网格划分模块,用于采用预设线型单元对所述钢筋模型进行划分,获取钢筋有限元模型,所述钢筋有限元模型包括N个钢筋单元,N为大于0的整数;a second meshing module, configured to divide the reinforcing bar model by using preset line elements to obtain a reinforcing bar finite element model, where the reinforcing bar finite element model includes N reinforcing bar elements, and N is an integer greater than 0;

计算模块,用于根据所述钢筋单元的位置信息、所述混凝土单元的位置信息,获取各混凝土单元的配筋率。The calculation module is configured to obtain the reinforcement ratio of each concrete unit according to the position information of the reinforcement unit and the position information of the concrete unit.

进一步地,所述计算模块,还用于在根据所述钢筋单元的位置信息、所述混凝土单元的位置信息,获取各混凝土单元的配筋率之前,根据所述混凝土单元的形函数、所述钢筋单元的位置信息、所述混凝土单元的位置信息,建立所述钢筋单元与所述混凝土单元的映射关系。Further, the calculation module is further configured to, before obtaining the reinforcement ratio of each concrete unit according to the position information of the reinforcement unit and the position information of the concrete unit, according to the shape function of the concrete unit, the The location information of the reinforcement unit and the location information of the concrete unit establish a mapping relationship between the reinforcement unit and the concrete unit.

进一步地,所述计算模块,具体用于根据所述混凝土单元的形函数、所述混凝土单元的位置信息,对所述混凝土单元等参变换,并采用牛顿迭代法计算,计算各所述钢筋单元端点的局部坐标;Further, the calculation module is specifically used for isoparametric transformation of the concrete unit according to the shape function of the concrete unit and the position information of the concrete unit, and calculates by using the Newton iteration method to calculate each of the steel bar units. the local coordinates of the endpoint;

将所述钢筋单元端点的局部坐标所在的混凝土单元,与所述钢筋单元建立映射关系。A mapping relationship is established between the concrete unit where the local coordinates of the end points of the reinforcing bar unit are located and the reinforcing bar unit.

进一步地,所述计算模块,具体用于根据所述钢筋单元与所述混凝土单元的映射关系,遍历各混凝土单元存在映射关系的钢筋单元;Further, the calculation module is specifically configured to traverse the reinforcing bar units that have a mapping relationship in each concrete unit according to the mapping relationship between the reinforcing bar unit and the concrete unit;

根据各混凝土单元存在映射关系的钢筋单元中,各钢筋单元位于对应混凝土单元内的端点数,获取各混凝土单元的配筋率。The reinforcement ratio of each concrete element is obtained according to the number of end points of each reinforcement element located in the corresponding concrete element among the reinforcement elements in which each concrete element has a mapping relationship.

进一步地,所述计算模块,具体用于计算各混凝土单元存在映射关系的钢筋单元在对应混凝土单元内的体积,其中,若钢筋单元位于对应混凝土单元内的端点数为1,则所述钢筋单元在对应混凝土单元内的体积为所述钢筋单元实际体积的1/2;若钢筋单元位于对应混凝土单元内的端点数为2,则所述钢筋单元在对应混凝土单元内的体积为所述钢筋单元实际体积;Further, the calculation module is specifically used to calculate the volume of the reinforcing bar unit with the mapping relationship of each concrete unit in the corresponding concrete unit, wherein, if the number of end points of the reinforcing bar unit in the corresponding concrete unit is 1, then The volume in the corresponding concrete unit is 1/2 of the actual volume of the reinforcement unit; if the number of end points of the reinforcement unit in the corresponding concrete unit is 2, the volume of the reinforcement unit in the corresponding concrete unit is the reinforcement unit actual volume;

根据各混凝土单元存在映射关系的钢筋单元在对应混凝土单元内的体积,获取所述各混凝土单元的配筋率。The reinforcement ratio of each concrete unit is obtained according to the volume of the reinforcing bar unit with the mapping relationship of each concrete unit in the corresponding concrete unit.

本发明提供的混凝土单元配筋率的计算方法及装置中,通过获取目标建筑物的混凝土模型和钢筋模型,采用预设立体几何单元对所述混凝土模型进行划分,获取混凝土有限元模型,所述混凝土有限元模型包括M个混凝土单元,M为大于0的整数,采用预设线型单元对所述钢筋模型进行划分,获取钢筋有限元模型,所述钢筋有限元模型包括N个钢筋单元,N为大于0的整数,根据所述钢筋单元的位置信息、所述混凝土单元的位置信息,获取各混凝土单元的配筋率。即具体通过分别获取目标建筑物的混凝土有限元模型和钢筋有限元模型,并采用有限元等参变换形函数和牛顿迭代法进行计算,从而获取混凝土单元的配筋率,实现了能够适用于更多的复杂工程中混凝土单元配筋率的计算,并且获取的目标建筑物的混凝土单元的配筋率能够为后续钢筋混凝土联合作用有限元计算分析提供数据支持。In the method and device for calculating the reinforcement ratio of a concrete unit provided by the present invention, the concrete model and the reinforcement model of the target building are obtained, and the concrete model is divided by a preset three-dimensional geometric unit to obtain the concrete finite element model. The concrete finite element model includes M concrete units, where M is an integer greater than 0, and the steel bar model is divided by using preset line elements to obtain a steel bar finite element model, where the steel bar finite element model includes N steel bar units, N is an integer greater than 0, and the reinforcement ratio of each concrete unit is obtained according to the position information of the reinforcement unit and the position information of the concrete unit. That is, by obtaining the concrete finite element model and the steel finite element model of the target building respectively, and using the finite element isoparametric transformation function and the Newton iteration method for calculation, the reinforcement ratio of the concrete element is obtained, which can be applied to more It can calculate the reinforcement ratio of concrete elements in many complex projects, and the obtained reinforcement ratio of concrete elements of the target building can provide data support for the subsequent finite element calculation and analysis of reinforced concrete combined action.

附图说明Description of drawings

为了更清楚地说明本发明各实施例,下面将对各实施例描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。In order to explain the various embodiments of the present invention more clearly, the following briefly introduces the accompanying drawings used in the description of each embodiment. Obviously, the accompanying drawings in the following description are some embodiments of the present invention. For those of ordinary skill in the art, other drawings can also be obtained from these drawings without any creative effort.

图1是本发明提供的混凝土单元配筋率的计算方法实施例一的流程示意图;1 is a schematic flow chart of Embodiment 1 of a method for calculating the reinforcement ratio of a concrete unit provided by the present invention;

图2是混凝土有限元模型中局部坐标系示意图;Figure 2 is a schematic diagram of the local coordinate system in the concrete finite element model;

图3是本发明提供的混凝土单元配筋率的计算装置实施例一的结构示意图;3 is a schematic structural diagram of Embodiment 1 of the computing device for the reinforcement ratio of concrete units provided by the present invention;

图4是本发明提供的混凝土单元配筋率的计算装置实施例二的结构示意图。4 is a schematic structural diagram of Embodiment 2 of the device for calculating the reinforcement ratio of concrete units provided by the present invention.

具体实施方式Detailed ways

为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。In order to make the purposes, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments These are some embodiments of the present invention, but not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

图1为本发明提供的混凝土单元配筋率的计算方法实施例一的流程示意图,如图1所示,该方法包括:FIG. 1 is a schematic flowchart of Embodiment 1 of the method for calculating the reinforcement ratio of concrete units provided by the present invention. As shown in FIG. 1 , the method includes:

S101、获取目标建筑物的混凝土模型和钢筋模型。S101. Acquire a concrete model and a steel bar model of the target building.

根据实体目标建筑物在整体坐标系下建立分离式钢筋混凝土模型,也就是,根据实体目标建筑物在整体坐标系下分别建立混凝土模型和钢筋模型。并且,在该模型中假设混凝土与钢筋之间粘结良好,不会有相对滑移。The separate reinforced concrete model is established under the overall coordinate system according to the entity target building, that is, the concrete model and the steel bar model are respectively established under the overall coordinate system according to the entity target building. Also, in this model, it is assumed that the bond between the concrete and the rebar is good and there will be no relative slippage.

S102、采用预设立体几何单元对所述混凝土模型进行划分,获取混凝土有限元模型,所述混凝土有限元模型包括M个混凝土单元,M为大于0的整数。S102. Use a preset three-dimensional geometric unit to divide the concrete model to obtain a concrete finite element model, where the concrete finite element model includes M concrete units, where M is an integer greater than 0.

可选地,上述预设立体几何单元可为三维8节点六面体单元、三维6节点三棱柱单元、三维4节点四面体单元,需要说明的是,上述预设立体几何单元也可为其他立体几何单元,并不限定于上述三种立体几何单元,本发明对此不做限定。Optionally, the preset three-dimensional geometric unit may be a three-dimensional 8-node hexahedron unit, a three-dimensional 6-node triangular prism unit, and a three-dimensional 4-node tetrahedron unit. It should be noted that the above-mentioned preset three-dimensional geometric unit may also be other three-dimensional geometric units. , is not limited to the above three three-dimensional geometric units, which is not limited in the present invention.

在对上述混凝土模型进行划分,获取混凝土有限元模型时,预设立体几何单元的尺寸不宜过大,预设立体几何单元的尺寸越大,计算量虽然减小,但是得到的结果精确度低,预设立体几何单元的尺寸越小,得到的计算结果精确度更高,但是计算量较大,因此,针对不同的目标建筑物的混凝土模型可根据实际情况选取合适尺寸的立体几何单元进行划分,获取有限元模型。可选地一种实现方式中,预设立体几何单元各边长的尺寸可为1~2米。When dividing the above concrete model and obtaining the concrete finite element model, the size of the preset three-dimensional geometric element should not be too large. The smaller the size of the preset three-dimensional geometric unit, the higher the accuracy of the calculation result, but the larger the amount of calculation. Therefore, the concrete model for different target buildings can be divided into three-dimensional geometric units of suitable size according to the actual situation. Obtain a finite element model. Optionally, in an implementation manner, the length of each side of the preset three-dimensional geometric unit may be 1-2 meters.

进一步地,上述混凝土有限元模型包括M个混凝土单元,每个混凝土单元有多个节点。例如,采用边长为1米的正方体对混凝土模型进行划分,那么每个混凝土单元有8个节点。Further, the above-mentioned concrete finite element model includes M concrete elements, and each concrete element has a plurality of nodes. For example, if the concrete model is divided by a cube with a side length of 1 meter, then each concrete element has 8 nodes.

需要说明的是,在实际应用中,可以根据实际情况选择一种或多种立体几何单元对目标建筑物的混凝土模型进行划分,也就是说针对目标建筑物的混凝土模型的不同区域可以选择不同的立体几何单元进行划分,使获得的混凝土有限元模型更为精确,本发明对此不做限制。It should be noted that, in practical applications, one or more solid geometric units can be selected to divide the concrete model of the target building according to the actual situation, that is to say, different areas of the concrete model of the target building can be selected. The three-dimensional geometric unit is divided to make the obtained concrete finite element model more accurate, which is not limited in the present invention.

S103、采用预设线型单元对所述钢筋模型进行划分,获取钢筋有限元模型,所述钢筋有限元模型包括N个钢筋单元,N为大于0的整数。S103. Use preset line elements to divide the steel bar model to obtain a steel bar finite element model, where the steel bar finite element model includes N steel bar units, and N is an integer greater than 0.

由于钢筋是一种细长的材料,通常可忽略其横向抗剪强度,认为其仅承受顺长方向的轴力,因此,对钢筋模型采用线型单元进行划分获取钢筋有限元模型。例如:采用三维2节点线单元对钢筋模型进行划分。Since the steel bar is a slender material, its transverse shear strength can usually be ignored, and it is considered that it only bears the axial force in the longitudinal direction. Therefore, the steel bar model is divided by linear elements to obtain the steel bar finite element model. For example: use three-dimensional 2-node line elements to divide the reinforcement model.

其中,线型单元的尺寸越小,得到的计算结果精确度高,线型单元的尺寸越大,得到的结果精确度越低。为了精细模拟钢筋混凝土的联合作用效果,减小计算误差,可选地一种实现方式中,线型单元的尺寸可为混凝土立体几何单元尺寸的5%。例如,当采用边长为1米的正方体对混凝土模型进行划分,那么线型单元的尺寸为0.05米。Among them, the smaller the size of the linear element, the higher the accuracy of the obtained calculation result, and the larger the size of the linear element, the lower the accuracy of the obtained result. In order to accurately simulate the combined effect of reinforced concrete and reduce the calculation error, optionally in an implementation manner, the size of the linear unit may be 5% of the size of the concrete solid geometric unit. For example, when a concrete model is partitioned with a 1-meter cube, the linear element size is 0.05 meters.

进一步地,上述钢筋有限元模型包括有N个钢筋单元,每个钢筋单元包括有两个端点。并对上述N个钢筋单元的所有端点进行编号。Further, the above-mentioned finite element model of reinforcing bars includes N reinforcing bar units, and each reinforcing bar unit includes two end points. And number all the endpoints of the above N reinforcement elements.

S104、根据所述钢筋单元的位置信息、所述混凝土单元的位置信息,获取各混凝土单元的配筋率。S104. Acquire the reinforcement ratio of each concrete unit according to the location information of the reinforcing bar unit and the location information of the concrete unit.

本实施例,通过获取目标建筑物的混凝土模型和钢筋模型,采用预设立体几何单元对所述混凝土模型进行划分,获取混凝土有限元模型,所述混凝土有限元模型包括M个混凝土单元,M为大于0的整数,采用预设线型单元对所述钢筋模型进行划分,获取钢筋有限元模型,所述钢筋有限元模型包括N个钢筋单元,N为大于0的整数,根据所述钢筋单元的位置信息、所述混凝土单元的位置信息,获取各混凝土单元的配筋率。即具体通过分别获取目标建筑物的混凝土有限元模型和钢筋有限元模型,并采用有限元等参变换形函数和牛顿迭代法进行计算,从而获取混凝土单元的配筋率,实现了能够适用于更多的复杂工程中混凝土单元配筋率的计算,并且获取的目标建筑物的混凝土单元的配筋率能够为后续钢筋混凝土联合作用有限元计算分析提供数据支持。In this embodiment, by obtaining the concrete model and the steel bar model of the target building, the concrete model is divided by using preset three-dimensional geometric units, and the concrete finite element model is obtained. The concrete finite element model includes M concrete units, where M is If the integer is greater than 0, the steel bar model is divided by the preset line element, and the steel bar finite element model is obtained. The steel bar finite element model includes N steel bar units, and N is an integer greater than 0. The location information and the location information of the concrete unit are used to obtain the reinforcement ratio of each concrete unit. That is, by obtaining the concrete finite element model and the steel finite element model of the target building respectively, and using the finite element isoparametric transformation function and the Newton iteration method for calculation, the reinforcement ratio of the concrete element is obtained, which can be applied to more It can calculate the reinforcement ratio of concrete elements in many complex projects, and the obtained reinforcement ratio of concrete elements of the target building can provide data support for the subsequent finite element calculation and analysis of reinforced concrete combined action.

可选地一种实现方式中,上述实施例中根据所述钢筋单元的位置信息、所述混凝土单元的位置信息,获取各混凝土单元的配筋率,可以有以下过程:Optionally, in an implementation manner, in the above embodiment, the reinforcement ratio of each concrete unit is obtained according to the position information of the reinforcement unit and the position information of the concrete unit, and the following process may be performed:

首先,根据所述混凝土单元的形函数、所述钢筋单元的位置信息、所述混凝土单元的位置信息,建立所述钢筋单元与所述混凝土单元的映射关系。First, according to the shape function of the concrete unit, the location information of the reinforcing bar unit, and the location information of the concrete unit, a mapping relationship between the reinforcing bar unit and the concrete unit is established.

具体地,根据所述混凝土单元的形函数、所述混凝土单元的位置信息,对所述混凝土单元等参变换,并采用牛顿迭代法计算,计算各所述钢筋单元端点的局部坐标。Specifically, according to the shape function of the concrete unit and the position information of the concrete unit, isoparametric transformation of the concrete unit is performed, and the Newton iteration method is used for calculation to calculate the local coordinates of each end point of the steel bar unit.

将所述钢筋单元端点的局部坐标所在的混凝土单元,与所述钢筋单元建立映射关系。A mapping relationship is established between the concrete unit where the local coordinates of the end points of the reinforcing bar unit are located and the reinforcing bar unit.

这里需要说明的是,获取所述钢筋单元端点在某一混凝土单元内的局部坐标时,首先根据该混凝土单元建立局部坐标系。例如,采用三维八节点六面体单元对混凝土模型进行划分后,相应的混凝土单元示意图以及根据该混凝土单元建立的局部坐标系示意图如图2所示,当然也可以以该混凝土单元中的其他节点为原点建立局部坐标系,针对同一所述混凝土单元建立不同的局部坐标系并不影响计算结果。当采用其他类型的立体几何单元对混凝土模型进行划分时,相应的混凝土单元以及根据该混凝土单元建立的局部坐标系的方法与现有方法中类似,本发明在此不做详细描述。It should be noted here that, when obtaining the local coordinates of the end points of the reinforcing bar unit in a certain concrete unit, a local coordinate system is first established according to the concrete unit. For example, after the concrete model is divided by a three-dimensional eight-node hexahedron element, the schematic diagram of the corresponding concrete element and the schematic diagram of the local coordinate system established according to the concrete element are shown in Figure 2. Of course, other nodes in the concrete element can also be used as the origin. Establishing a local coordinate system and establishing different local coordinate systems for the same concrete unit does not affect the calculation results. When other types of solid geometric elements are used to divide the concrete model, the corresponding concrete elements and the method for establishing the local coordinate system based on the concrete elements are similar to the existing methods, and the present invention will not describe them in detail here.

在实际应用中,以混凝土单元的有限元等参变换为基础,根据钢筋单元的端点在整体坐标系下的坐标计算出其对应的局部坐标,进而判断钢筋单元的端点的局部坐标是否在该混凝土单元内。In practical applications, based on the finite element isoparametric transformation of the concrete element, the corresponding local coordinates are calculated according to the coordinates of the end points of the reinforcement element in the global coordinate system, and then it is judged whether the local coordinates of the end points of the reinforcement element are in the concrete within the unit.

局部坐标中几何形状规则的单元到整体坐标中几何形状扭曲的单元,坐标转换如下:From a geometrically regular cell in local coordinates to a geometrically distorted cell in global coordinates, the coordinate transformation is as follows:

Figure BDA0001608203480000071
Figure BDA0001608203480000071

其中,m为混凝土单元的节点个数,xi、yi、zi是某混凝土单元节点在整体坐标中的坐标值,ξ、η、ζ为局部坐标系,Ni(ξ,η,ζ)为用局部坐标表示的形函数,Ni(ξ,η,ζ)可以从有关有限元方法书籍中根据预设立体几何单元类型获取到。Among them, m is the number of nodes of the concrete element, xi , y i , zi are the coordinate values of a concrete element node in the global coordinate, ξ, η, ζ are the local coordinate system, N i (ξ,η,ζ ) is the shape function expressed in local coordinates, and Ni (ξ,η, ζ ) can be obtained from books on finite element methods according to the preset solid geometry element type.

根据公式(1),令:According to formula (1), let:

Figure BDA0001608203480000072
Figure BDA0001608203480000072

已知其中一个所述钢筋单元端点的整体坐标为

Figure BDA0001608203480000081
利用方程
Figure BDA0001608203480000082
求所述钢筋单元端点在该混凝土单元中的局部坐标
Figure BDA0001608203480000083
The global coordinates of one of the end points of the reinforcement element are known as
Figure BDA0001608203480000081
Use equations
Figure BDA0001608203480000082
Find the local coordinates of the end points of the reinforcement element in the concrete element
Figure BDA0001608203480000083

令:make:

Figure BDA0001608203480000084
Figure BDA0001608203480000084

求解:Solve:

Figure BDA0001608203480000085
Figure BDA0001608203480000085

根据牛顿迭代法可得:According to Newton's iteration method, we can get:

Figure BDA0001608203480000086
Figure BDA0001608203480000086

其中,公式(5)中

Figure BDA0001608203480000087
为第n+1步的迭代结果,
Figure BDA0001608203480000088
为第n步的迭代结果。Among them, in formula (5)
Figure BDA0001608203480000087
is the iterative result of the n+1th step,
Figure BDA0001608203480000088
is the iterative result of the nth step.

Figure BDA0001608203480000089
Figure BDA0001608203480000089

Figure BDA0001608203480000091
Figure BDA0001608203480000091

令:make:

Figure BDA0001608203480000092
Figure BDA0001608203480000092

其中,公式(8)中

Figure BDA0001608203480000093
称为雅克比Jacobi矩阵,
Figure BDA0001608203480000094
Figure BDA0001608203480000095
的逆矩阵。Among them, in formula (8)
Figure BDA0001608203480000093
is called the Jacobi matrix,
Figure BDA0001608203480000094
for
Figure BDA0001608203480000095
The inverse matrix of .

Figure BDA0001608203480000096
代入(5)、(6)、(8)中可得:Will
Figure BDA0001608203480000096
Substitute into (5), (6), (8) to get:

Figure BDA0001608203480000097
Figure BDA0001608203480000097

若满足|ξn+1n|≤ε,|ηn+1n|≤ε,|ζn+1n|≤ε,迭代结束,ε为设置的收敛容差,则整体坐标为

Figure BDA0001608203480000098
的所述钢筋单元端点在所述混凝土单元中的局部坐标为:If |ξ n+1n |≤ε, |η n+ 1 -η n |≤ε, |ζ n+1n |≤ε, the iteration ends, and ε is the set convergence tolerance, then The overall coordinates are
Figure BDA0001608203480000098
The local coordinates of the end points of the reinforcement element in the concrete element are:

Figure BDA0001608203480000099
Figure BDA0001608203480000099

进一步,根据所述钢筋单元端点在所述混凝土单元内的局部坐标与所述混凝土单元的位置信息进行判断。Further, the determination is made according to the local coordinates of the end points of the reinforcing bar unit in the concrete unit and the position information of the concrete unit.

若所述钢筋单元端点的局部坐标在所述混凝土单元内,则建立起所述钢筋单元端点号与所述混凝土单元号之间的映射关系。If the local coordinates of the end points of the reinforcing bar units are in the concrete unit, a mapping relationship between the end numbers of the reinforcing bar units and the concrete unit numbers is established.

若所述钢筋单元端点的局部坐标不在所述混凝土单元内,则遍历下一个混凝土单元,直至找到与所述钢筋单元端点对应的混凝土单元,并建立相应的映射关系。If the local coordinates of the end point of the reinforcing bar unit are not in the concrete unit, traverse the next concrete unit until a concrete unit corresponding to the end point of the reinforcing bar unit is found, and a corresponding mapping relationship is established.

需要说明的是,通过坐标位置判断所述钢筋单元端点的局部坐标是否在所述混凝土单元内的方法与现有技术类似,本发明在此不做详细描述。It should be noted that the method of judging whether the local coordinates of the end points of the reinforcing bar unit are in the concrete unit by the coordinate position is similar to the prior art, and the present invention will not describe it in detail here.

可选地一种实现方式中,若所述钢筋单元端点的局部坐标不在所述混凝土单元内,按照混凝土单元号由小到大的顺序进行遍历。Optionally, in an implementation manner, if the local coordinates of the end points of the reinforcing bar units are not in the concrete unit, the traversal is performed according to the order of the concrete unit numbers from small to large.

进一步地,根据所述钢筋单元与所述混凝土单元的映射关系,遍历各混凝土单元存在映射关系的钢筋单元。Further, according to the mapping relationship between the reinforcing bar unit and the concrete unit, traverse the reinforcing bar units with the mapping relationship of each concrete unit.

根据各混凝土单元存在映射关系的钢筋单元中,各钢筋单元位于对应混凝土单元内的端点数,获取各混凝土单元的配筋率。The reinforcement ratio of each concrete element is obtained according to the number of end points of each reinforcement element located in the corresponding concrete element among the reinforcement elements in which each concrete element has a mapping relationship.

具体地,针对某一所述混凝土单元,对所述钢筋单元进行遍历,根据上述建立起来的映射关系判断所述钢筋单元在所述混凝土单元内的钢筋单元端点个数,若钢筋单元位于对应混凝土单元内的端点数为1,则所述钢筋单元在对应混凝土单元内x、y、z方向上的体积为所述钢筋单元实际体积的1/2与钢筋单元方向矢量在x、y、z方向分量的乘积;若钢筋单元位于对应混凝土单元内的端点数为2,则所述钢筋单元在对应混凝土单元内x、y、z方向的体积为所述钢筋单元实际体积与钢筋单元方向矢量在x、y、z方向分量的乘积。当钢筋单元遍历结束后,分别累加该混凝土单元内所有钢筋单元x、y、z方向上的体积后除以所述混凝土的单元的体积,得到的计算结果即为该混凝土单元沿x、y、z方向上的配筋率ρx、ρy、ρzSpecifically, for a certain concrete unit, the reinforcement unit is traversed, and the number of end points of the reinforcement unit in the concrete unit is determined according to the mapping relationship established above. If the reinforcement unit is located in the corresponding concrete The number of endpoints in a unit is 1, then the volume of the reinforcing bar unit in the x, y, and z directions of the corresponding concrete unit is 1/2 of the actual volume of the reinforcing bar unit and the direction vector of the reinforcing bar unit is in the x, y, and z directions. If the number of end points of the reinforcement element in the corresponding concrete element is 2, the volume of the reinforcement element in the x, y, and z directions in the corresponding concrete element is the actual volume of the reinforcement element and the direction vector of the reinforcement element in x The product of the , y, and z direction components. When the traversal of the reinforcement element is completed, the volumes of all reinforcement elements in the concrete element in the x, y, and z directions are accumulated and divided by the volume of the concrete element. Reinforcement ratios ρ x , ρ y , ρ z in the z direction.

需要说明的是,在执行遍历时,若钢筋单元与混凝土单元不存在映射关系时,则钢筋单元在该混凝土单元内的体积为0。It should be noted that, when the traversal is performed, if there is no mapping relationship between the reinforcement element and the concrete element, the volume of the reinforcement element in the concrete element is 0.

需要说明的是,所述钢筋单元实际体积为所述钢筋单元横截面积与所述线型单元长度的乘积。It should be noted that the actual volume of the reinforcing bar unit is the product of the cross-sectional area of the reinforcing bar unit and the length of the linear unit.

进一步地,开始遍历所述混凝土有限元模型中的各所述混凝土单元,直至获取到所有所述混凝土单元的配筋率,结束对混凝土单元的遍历。Further, start to traverse each of the concrete units in the concrete finite element model until the reinforcement ratios of all the concrete units are obtained, and end the traversal of the concrete units.

可选地一种遍历方式为按照钢筋单元号与混凝土单元号由小到大的顺序进行遍历。Optionally, a traversal method is to traverse according to the order of the reinforcement unit numbers and the concrete unit numbers from small to large.

例如,混凝土有限元模型有M个混凝土单元,钢筋有限元模型有N个钢筋单元,那么针对第一个混凝土单元,开始遍历钢筋有限元模型中的N个钢筋单元。首先,根据映射关系判断第一个钢筋单元在上述第一个混凝土单元中的钢筋单元端点个数,获取第一个钢筋单元在第一个混凝土单元内x、y、z方向上的相应的钢筋单元体积,接着根据映射关系判断第二个钢筋单元在上述第一个混凝土单元中的钢筋单元端点个数,获取第二个钢筋单元在第一个混凝土单元内x、y、z方向上的相应的钢筋单元体积,以此类推,获取N个钢筋单元在第一个混凝土单元内x、y、z方向上的相应的钢筋单元体积。接着将钢筋有限元模型中N个钢筋单元在第一个混凝土单元内x、y、z方向上的相应的钢筋体积进行求和,再除以第一个混凝土单元的体积,获取第一个混凝土单元在x、y、z方向上的配筋率。For example, the concrete finite element model has M concrete elements, and the reinforcement finite element model has N reinforcement elements, then for the first concrete element, start to traverse the N reinforcement elements in the reinforcement finite element model. First, according to the mapping relationship, determine the number of end points of the first reinforcing bar element in the first concrete unit, and obtain the corresponding reinforcing bars of the first reinforcing bar unit in the x, y, and z directions of the first concrete unit. The unit volume, and then according to the mapping relationship, determine the number of end points of the second reinforcing bar unit in the first concrete unit, and obtain the corresponding x, y, and z directions of the second reinforcing bar unit in the first concrete unit. , and so on, to obtain the corresponding reinforcement unit volumes of N reinforcement elements in the x, y, and z directions of the first concrete element. Then, the corresponding reinforcement volumes of the N reinforcement elements in the reinforcement finite element model in the x, y, and z directions of the first concrete element are summed, and then divided by the volume of the first concrete element to obtain the first concrete element. The reinforcement ratio of the element in the x, y, and z directions.

进一步地,针对第二个混凝土单元,开始遍历所有钢筋有限元模型中的N个钢筋单元,通过上述同样的方法获取第二个混凝土单元在x、y、z方向上的配筋率。以此类推,遍历混凝土有限元模型中的M个混凝土单元,获取M个混凝土单元在x、y、z方向上的配筋率。Further, for the second concrete element, start to traverse all N reinforcement elements in the finite element model of reinforcement, and obtain the reinforcement ratio of the second concrete element in the x, y, and z directions by the same method as above. By analogy, the M concrete elements in the concrete finite element model are traversed, and the reinforcement ratios of the M concrete elements in the x, y, and z directions are obtained.

需要说明的是,本发明中也可以按照其他的顺序对钢筋单元和混凝土单元进行遍历,例如:可以按照钢筋单元号和混凝土单元号由大到小的顺序进行遍历,本发明对此不做限制。It should be noted that, in the present invention, the reinforcement unit and the concrete unit can also be traversed in other order. For example, the reinforcement unit number and the concrete unit number can be traversed in descending order, which is not limited in the present invention. .

上述根据所述钢筋单元的位置信息、所述混凝土单元的位置信息,获取各混凝土单元的配筋率,另一种实现方式,可以为:According to the position information of the reinforcing bar unit and the position information of the concrete unit, the reinforcement ratio of each concrete unit is obtained. Another implementation manner may be:

首先,根据所述混凝土单元的形函数、所述混凝土单元的位置信息,对所述混凝土单元等参变换,并采用牛顿迭代法计算,计算各所述钢筋单元端点的局部坐标;First, according to the shape function of the concrete element and the position information of the concrete element, the isoparametric transformation of the concrete element is carried out, and the Newton iteration method is used for calculation to calculate the local coordinates of the end points of each of the steel reinforcement elements;

根据各所述钢筋单元端点的局部坐标直接判断是否在所述混凝土单元内,若钢筋单元位于对应混凝土单元内的端点数为1,则所述钢筋单元在对应混凝土单元内x、y、z方向上的体积为所述钢筋单元实际体积的1/2与钢筋单元方向矢量在x、y、z方向分量的乘积;若钢筋单元位于对应混凝土单元内的端点数为2,则所述钢筋单元在对应混凝土单元内x、y、z方向的体积为所述钢筋单元实际体积与钢筋单元方向矢量在x、y、z方向分量的乘积。According to the local coordinates of the end points of each reinforcing bar unit, it is directly judged whether it is in the concrete unit. If the number of end points of the reinforcing bar unit in the corresponding concrete unit is 1, the reinforcing bar unit is in the x, y, and z directions of the corresponding concrete unit. The volume above is the product of 1/2 of the actual volume of the reinforcement element and the components of the direction vector of the reinforcement element in the x, y, and z directions; if the number of endpoints of the reinforcement element in the corresponding concrete element is 2, the reinforcement element is in the The volume corresponding to the x, y, and z directions in the concrete unit is the product of the actual volume of the reinforcing bar unit and the components of the direction vector of the reinforcing bar unit in the x, y, and z directions.

需要说明的是,计算各所述钢筋单元端点的局部坐标的方法与前述方法相同。It should be noted that the method for calculating the local coordinates of the end points of each of the reinforcing bar elements is the same as the aforementioned method.

本实施例通过获取目标建筑物的混凝土模型和钢筋模型,采用预设立体几何单元对所述混凝土模型进行划分,获取混凝土有限元模型,所述混凝土有限元模型包括M个混凝土单元,M为大于0的整数,采用预设线型单元对所述钢筋模型进行划分,获取钢筋有限元模型,所述钢筋有限元模型包括N个钢筋单元,N为大于0的整数,根据所述钢筋单元的位置信息、所述混凝土单元的位置信息,获取各混凝土单元的配筋率。即具体通过分别获取目标建筑物的混凝土有限元模型和钢筋有限元模型,并采用有限元等参变换形函数和牛顿迭代法进行计算,从而获取混凝土单元的配筋率,实现了能够适用于更多的复杂工程中混凝土单元配筋率的计算,并且获取的目标建筑物的混凝土单元的配筋率能够为后续钢筋混凝土联合作用有限元计算分析提供数据支持。In this embodiment, the concrete model and the steel bar model of the target building are obtained, and the concrete model is divided by the preset three-dimensional geometric elements, and the concrete finite element model is obtained. The concrete finite element model includes M concrete elements, and M is greater than or equal to an integer of 0, the steel bar model is divided by a preset line type element, and a steel bar finite element model is obtained. The steel bar finite element model includes N steel bar elements, and N is an integer greater than 0. According to the position of the steel bar element information, the location information of the concrete unit, and obtain the reinforcement ratio of each concrete unit. That is, by obtaining the concrete finite element model and the steel finite element model of the target building respectively, and using the finite element isoparametric transformation function and the Newton iteration method for calculation, the reinforcement ratio of the concrete element is obtained, which can be applied to more It can calculate the reinforcement ratio of concrete elements in many complex projects, and the obtained reinforcement ratio of concrete elements of the target building can provide data support for the subsequent finite element calculation and analysis of reinforced concrete combined action.

图3为本发明提供的混凝土单元配筋率的计算装置实施例一的结构示意图。本发明提供的混凝土单元配筋率的计算装置可以用于执行图1所示方法实施例的技术方案。如图3所示,本实施例提供的混凝土单元配筋率的计算装置包括:FIG. 3 is a schematic structural diagram of Embodiment 1 of the apparatus for calculating the reinforcement ratio of concrete units provided by the present invention. The device for calculating the reinforcement ratio of a concrete unit provided by the present invention can be used to implement the technical solution of the method embodiment shown in FIG. 1 . As shown in FIG. 3 , the calculation device for the reinforcement ratio of concrete units provided in this embodiment includes:

获取模块31,用于获取目标建筑物的混凝土模型和钢筋模型。The acquiring module 31 is used to acquire the concrete model and the steel reinforcement model of the target building.

第一网格划分模块32,用于采用预设立体几何单元对所述混凝土模型进行划分,获取混凝土有限元模型,所述混凝土有限元模型包括M个混凝土单元,M为大于0的整数。The first mesh dividing module 32 is configured to divide the concrete model by using preset three-dimensional geometric units to obtain a concrete finite element model, where the concrete finite element model includes M concrete units, where M is an integer greater than 0.

第二网格划分模块33,用于采用预设线型单元对所述钢筋模型进行划分,获取钢筋有限元模型,所述钢筋有限元模型包括N个钢筋单元,N为大于0的整数。The second meshing module 33 is configured to divide the reinforcing bar model by using preset line elements to obtain a reinforcing bar finite element model, where the reinforcing bar finite element model includes N reinforcing bar elements, and N is an integer greater than 0.

计算模块34,用于根据所述钢筋单元的位置信息、所述混凝土单元的位置信息,获取各混凝土单元的配筋率。The calculation module 34 is configured to obtain the reinforcement ratio of each concrete unit according to the position information of the reinforcement unit and the position information of the concrete unit.

本实施例提供的装置,即具体通过分别获取目标建筑物的混凝土有限元模型和钢筋有限元模型,并采用有限元等参变换形函数和牛顿迭代法进行计算,从而获取混凝土单元的配筋率,实现了能够适用于更多的复杂工程中混凝土单元配筋率的计算,并且获取的目标建筑物的混凝土单元的配筋率能够为后续钢筋混凝土联合作用有限元计算分析提供数据支持。The device provided in this embodiment specifically obtains the concrete finite element model and the reinforcement finite element model of the target building, and uses the finite element isoparametric transformation function and the Newton iteration method to calculate, so as to obtain the reinforcement ratio of the concrete element. , which can be applied to the calculation of the reinforcement ratio of concrete elements in more complex projects, and the obtained reinforcement ratio of concrete elements of the target building can provide data support for the subsequent finite element calculation and analysis of reinforced concrete combined action.

在本发明的一个实施例中,所述计算模块34,还用于在根据所述钢筋单元的位置信息、所述混凝土单元的位置信息,获取各混凝土单元的配筋率之前,根据所述混凝土单元的形函数、所述钢筋单元的位置信息、所述混凝土单元的位置信息,建立所述钢筋单元与所述混凝土单元的映射关系。In an embodiment of the present invention, the calculation module 34 is further configured to, before acquiring the reinforcement ratio of each concrete unit according to the position information of the reinforcement unit and the position information of the concrete unit, calculate the reinforcement ratio according to the concrete unit. The shape function of the unit, the location information of the reinforcing bar unit, and the location information of the concrete unit are used to establish the mapping relationship between the reinforcing bar unit and the concrete unit.

具体地,计算模块34,具体用于根据所述混凝土单元的形函数、所述混凝土单元的位置信息,对所述混凝土单元等参变换,并采用牛顿迭代法计算,计算各所述钢筋单元端点的局部坐标;Specifically, the calculation module 34 is specifically configured to perform isoparametric transformation of the concrete unit according to the shape function of the concrete unit and the position information of the concrete unit, and use the Newton iteration method for calculation to calculate the end points of each of the steel bar units. the local coordinates of ;

将所述钢筋单元端点的局部坐标所在的混凝土单元,与所述钢筋单元建立映射关系。A mapping relationship is established between the concrete unit where the local coordinates of the end points of the reinforcing bar unit are located and the reinforcing bar unit.

进一步地,计算模块34,具体用于根据所述钢筋单元与所述混凝土单元的映射关系,遍历各混凝土单元存在映射关系的钢筋单元;Further, the calculation module 34 is specifically configured to traverse the reinforcing bar units in which the mapping relationship exists in each concrete unit according to the mapping relationship between the reinforcing bar unit and the concrete unit;

根据各混凝土单元存在映射关系的钢筋单元中,各钢筋单元位于对应混凝土单元内的端点数,获取各混凝土单元的配筋率。The reinforcement ratio of each concrete element is obtained according to the number of end points of each reinforcement element located in the corresponding concrete element among the reinforcement elements in which each concrete element has a mapping relationship.

进一步地,计算模块34,具体用于计算各混凝土单元存在映射关系的钢筋单元在对应混凝土单元内的体积,其中,若钢筋单元位于对应混凝土单元内的端点数为1,则所述钢筋单元在对应混凝土单元内的体积为所述钢筋单元实际体积的1/2;若钢筋单元位于对应混凝土单元内的端点数为2,则所述钢筋单元在对应混凝土单元内的体积为所述钢筋单元实际体积;Further, the calculation module 34 is specifically used to calculate the volume of the reinforcing bar unit in the corresponding concrete unit that each concrete unit has a mapping relationship, wherein, if the number of endpoints of the reinforcing bar unit in the corresponding concrete unit is 1, then the reinforcing bar unit is in the corresponding concrete unit. The volume in the corresponding concrete unit is 1/2 of the actual volume of the reinforcement unit; if the number of endpoints of the reinforcement unit in the corresponding concrete unit is 2, the volume of the reinforcement unit in the corresponding concrete unit is the actual volume of the reinforcement unit. volume;

根据各混凝土单元存在映射关系的钢筋单元在对应混凝土单元内的体积,获取所述各混凝土单元的配筋率。The reinforcement ratio of each concrete unit is obtained according to the volume of the reinforcing bar unit with the mapping relationship of each concrete unit in the corresponding concrete unit.

图4为本发明提供的混凝土单元配筋率的计算装置实施例二的结构示意图。该装置可以包括:存储器401和处理器402。FIG. 4 is a schematic structural diagram of Embodiment 2 of the device for calculating the reinforcement ratio of a concrete unit provided by the present invention. The apparatus may include: a memory 401 and a processor 402 .

存储器401可以是独立的物理单元,与处理器402可以通过总线连接。存储器401、处理器402也可以集成在一起,通过硬件实现等。The memory 401 may be an independent physical unit, and may be connected to the processor 402 through a bus. The memory 401 and the processor 402 may also be integrated together, implemented by hardware, or the like.

存储器401用于存储实现以上方法实施例,处理器402调用该程序,执行以上装置执行的方法实施例的操作。The memory 401 is used to store and implement the above method embodiments, and the processor 402 invokes the program to execute the operations of the method embodiments performed by the above devices.

可选地,当上述实施例的方法中的部分或全部通过软件实现时,上述装置也可以只包括处理器。用于存储程序的存储器位于上述装置之外,处理器通过电路/电线与存储器连接,用于读取并执行存储器中存储的程序。Optionally, when part or all of the methods in the foregoing embodiments are implemented by software, the foregoing apparatus may also only include a processor. The memory for storing the program is located outside the above-mentioned device, and the processor is connected to the memory through a circuit/wire for reading and executing the program stored in the memory.

处理器可以是中央处理器(Central Processing Unit,CPU),网络处理器(Network Processor,NP)或者CPU和NP的组合。The processor may be a central processing unit (Central Processing Unit, CPU), a network processor (Network Processor, NP), or a combination of CPU and NP.

处理器还可以进一步包括硬件芯片。上述硬件芯片可以是专用集成电路(application-specific integrated circuit,ASIC),可编程逻辑器件(programmablelogic device,PLD)或其组合。上述PLD可以是复杂可编程逻辑器件(complexprogrammable logic device,CPLD),现场可编程逻辑门阵列(field-programmable gatearray,FPGA),通用阵列逻辑(generic array logic,GAL)或其任意组合。The processor may further include a hardware chip. The above hardware chip may be an application-specific integrated circuit (ASIC), a programmable logic device (PLD) or a combination thereof. The above-mentioned PLD may be a complex programmable logic device (CPLD), a field-programmable gate array (FPGA), a generic array logic (GAL) or any combination thereof.

存储器可以包括易失性存储器(volatile memory),例如随机存取存储器(random-access memory,RAM);存储器也可以包括非易失性存储器(non-volatilememory),例如快闪存储器(flash memory),硬盘(hard disk drive,HDD)或固态硬盘(solid-state drive,SSD);存储器还可以包括上述种类的存储器的组合。The memory may include volatile memory (volatile memory), such as random-access memory (random-access memory, RAM); the memory may also include non-volatile memory (non-volatile memory), such as flash memory (flash memory), A hard disk drive (HDD) or a solid-state drive (SSD); the memory may also include a combination of the above types of memory.

本领域普通技术人员可以理解:实现上述各方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成。前述的程序可以存储于一计算机可读取存储介质中。该程序在执行时,执行包括上述各方法实施例的步骤;而前述的存储介质包括:ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。Those of ordinary skill in the art can understand that all or part of the steps of implementing the above method embodiments may be completed by program instructions related to hardware. The aforementioned program can be stored in a computer-readable storage medium. When the program is executed, the steps including the above method embodiments are executed; and the foregoing storage medium includes: ROM, RAM, magnetic disk or optical disk and other media that can store program codes.

最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, but not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: The technical solutions described in the foregoing embodiments can still be modified, or some or all of the technical features thereof can be equivalently replaced; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the technical solutions of the embodiments of the present invention. scope.

Claims (4)

1. A method for calculating the reinforcement ratio of a concrete unit is characterized by comprising the following steps:
acquiring a concrete model and a steel bar model of a target building;
dividing the concrete model by adopting a preset solid geometry unit to obtain a concrete finite element model, wherein the concrete finite element model comprises M concrete units, and M is an integer greater than 0;
dividing the reinforcing steel bar model by adopting a preset linear unit to obtain a reinforcing steel bar finite element model, wherein the reinforcing steel bar finite element model comprises N reinforcing steel bar units, and N is an integer greater than 0;
acquiring the reinforcement ratio of each concrete unit according to the position information of the steel bar unit and the position information of the concrete unit;
the method for acquiring the reinforcement ratio of each concrete unit according to the position information of the steel bar unit and the position information of the concrete unit comprises the following steps:
according to the shape function of the concrete unit and the position information of the concrete unit, performing equal parameter transformation on the concrete unit, calculating by adopting a Newton iteration method, and calculating the local coordinate of each reinforcing steel bar unit end point;
establishing a mapping relation between the concrete unit where the local coordinates of the end points of the steel bar units are located and the steel bar units;
traversing the steel bar units with the mapping relation of the concrete units according to the mapping relation of the steel bar units and the concrete units;
acquiring the volume of the reinforcement unit with the mapping relation in the corresponding concrete unit according to the number of the end points of each reinforcement unit in the reinforcement unit with the mapping relation in the corresponding concrete unit;
and acquiring the reinforcement ratio of each concrete unit according to the volume of the reinforcement unit with the mapping relation in each concrete unit in the corresponding concrete unit.
2. The method of claim 1, wherein if the number of the end points of the reinforcement unit located in the corresponding concrete unit is 1, the volume of the reinforcement unit in the corresponding concrete unit is 1/2 of the actual volume of the reinforcement unit; and if the number of the end points of the steel bar unit in the corresponding concrete unit is 2, the volume of the steel bar unit in the corresponding concrete unit is the actual volume of the steel bar unit.
3. A device for calculating reinforcement ratio of a concrete unit, comprising:
the acquisition module is used for acquiring a concrete model and a steel bar model of a target building;
the first grid division module is used for dividing the concrete model by adopting a preset solid geometric unit to obtain a concrete finite element model, wherein the concrete finite element model comprises M concrete units, and M is an integer greater than 0;
the second mesh division module is used for dividing the steel bar model by adopting preset linear units to obtain a steel bar finite element model, wherein the steel bar finite element model comprises N steel bar units, and N is an integer greater than 0;
the calculation module is used for acquiring the reinforcement ratio of each concrete unit according to the position information of the steel bar unit and the position information of the concrete unit;
the calculation module is specifically configured to perform isoparametric transformation on the concrete unit according to a shape function of the concrete unit and position information of the concrete unit, calculate by using a newton iteration method, and calculate local coordinates of end points of each reinforcement unit;
establishing a mapping relation between the concrete unit where the local coordinates of the end points of the steel bar units are located and the steel bar units;
traversing the steel bar units with the mapping relation of the concrete units according to the mapping relation of the steel bar units and the concrete units;
acquiring the volume of the reinforcement unit with the mapping relation in the corresponding concrete unit according to the number of the end points of each reinforcement unit in the reinforcement unit with the mapping relation in the corresponding concrete unit;
and acquiring the reinforcement ratio of each concrete unit according to the volume of the reinforcement unit with the mapping relation in each concrete unit in the corresponding concrete unit.
4. The apparatus of claim 3, wherein if the number of the end points of the reinforcement unit located in the corresponding concrete unit is 1, the volume of the reinforcement unit in the corresponding concrete unit is 1/2 of the actual volume of the reinforcement unit; and if the number of the end points of the steel bar unit in the corresponding concrete unit is 2, the volume of the steel bar unit in the corresponding concrete unit is the actual volume of the steel bar unit.
CN201810252692.7A 2018-03-26 2018-03-26 Calculation method and device for reinforcement ratio of concrete unit Active CN108614913B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810252692.7A CN108614913B (en) 2018-03-26 2018-03-26 Calculation method and device for reinforcement ratio of concrete unit

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810252692.7A CN108614913B (en) 2018-03-26 2018-03-26 Calculation method and device for reinforcement ratio of concrete unit

Publications (2)

Publication Number Publication Date
CN108614913A CN108614913A (en) 2018-10-02
CN108614913B true CN108614913B (en) 2020-08-04

Family

ID=63658796

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810252692.7A Active CN108614913B (en) 2018-03-26 2018-03-26 Calculation method and device for reinforcement ratio of concrete unit

Country Status (1)

Country Link
CN (1) CN108614913B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109558629B (en) * 2018-10-18 2023-05-26 国核电力规划设计研究院有限公司 Method and device for generating reinforcement arrangement scheme
CN109537796B (en) * 2018-11-30 2020-12-15 陕西国防工业职业技术学院 A calculation method of reinforcement ratio of reinforced concrete slab based on strength reduction method
CN111444655B (en) * 2020-04-08 2021-03-16 中国水利水电科学研究院 A kind of static and dynamic reinforcement method of hydraulic building structure
CN111677283A (en) * 2020-07-30 2020-09-18 湖南中天建设集团股份有限公司 Method for assembling concrete structure
CN112883621B (en) * 2021-03-10 2022-09-02 陕西华山建设集团有限公司 Method for acquiring positions of steel bar joints in irregular plate column shear wall under finite element analysis
CN113821860B (en) * 2021-09-30 2024-04-12 广联达科技股份有限公司 Method, device, equipment and storage medium for determining break point on steel bar graphic element
CN120277789B (en) * 2025-06-04 2025-09-02 北京华科软科技有限公司 Intelligent analysis method for minimum reinforcement area of power transmission line pile foundation

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102286916B (en) * 2011-07-13 2014-10-29 东南大学 Method for determining time varying reliability of prestressed concrete box girder bridge
CN102651041B (en) * 2011-12-28 2016-12-28 同济大学 Prestressed cable basic engineering method in large-span space structure
US9857351B2 (en) * 2015-11-09 2018-01-02 University Of Dammam Concrete reinforcement assembly, method of installation, and method to determine cyclic load response
CN206440719U (en) * 2017-01-25 2017-08-25 厦门理工学院 The measuring system of Non-load Concrete internal stress in concrete structure

Also Published As

Publication number Publication date
CN108614913A (en) 2018-10-02

Similar Documents

Publication Publication Date Title
CN108614913B (en) Calculation method and device for reinforcement ratio of concrete unit
Yang et al. Optimal finite element model with response surface methodology for concrete structures based on Terrestrial Laser Scanning technology
CN116629079B (en) Method and device for constructing mixed finite element space and solving linear elastic mechanical problem
CN103838913B (en) The Finite Element of the curved bridge of curved box girder
CN108920786B (en) Section uncertainty analysis method based on Chebyshev polynomial fitting
Yang et al. A simple Galerkin meshless method, the fragile points method using point stiffness matrices, for 2D linear elastic problems in complex domains with crack and rupture propagation
CN105808884A (en) Method for predicting upper and lower bounds of bounded uncertain plane crack stress intensity factors based on fractal theory
CN112580212A (en) Iterative determination method, system and computer equipment for service life of concrete structure
CN106709219A (en) Area initial stress field inversion method and device under complicated geological conditions
JP5227384B2 (en) Simulation method using structured grid
CN111310349A (en) Data processing and analysis method suitable for continuous display of discrete element computing information
Tang et al. An efficient adaptive analysis procedure using the edge-based smoothed point interpolation method (ES-PIM) for 2D and 3D problems
CN111460715B (en) A Static Analysis Method of Natural Compensator Based on APDL Language
CN110543691B (en) Planar fracture grouting quasi-three-dimensional simulation method and system based on finite volume method
CN114444416B (en) Grid aggregation method and device for fluid dynamics simulation
CN115577597A (en) Simulation method, device, medium and equipment for solute hyperdiffusion in fracture channel
Tanbay et al. Numerical solution of the multigroup neutron diffusion equation by the meshless RBF collocation method
Izadian et al. The generalized finite difference method for solving elliptic equation on irregular mesh
CN113532308B (en) Ridge regression strain measurement method with initial value in digital image correlation
Yu et al. An h-adaptive numerical manifold method for solid mechanics problems
CN106909699B (en) Analysis method of thin plate vibration characteristics based on Galerkin strip transfer function
CN114492238A (en) Depth coupling method for one-dimensional and planar two-dimensional hydrodynamic models of rivers and lakes
CN114818523A (en) A FDM-DDA-Based Calculation Method for Unsaturated Transient Fluid-Structure Interaction
CN114781018A (en) Method for calculating tensile crack resistance reliability of asphalt concrete core wall
Gan et al. Superconvergent second derivative recovery technique and its application in a nonlocal damage mechanics model

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant