CN103778344A - Double-scale calculation method for tube sheet strength - Google Patents
Double-scale calculation method for tube sheet strength Download PDFInfo
- Publication number
- CN103778344A CN103778344A CN201410050862.5A CN201410050862A CN103778344A CN 103778344 A CN103778344 A CN 103778344A CN 201410050862 A CN201410050862 A CN 201410050862A CN 103778344 A CN103778344 A CN 103778344A
- Authority
- CN
- China
- Prior art keywords
- tube
- tube sheet
- stress
- equivalent
- region
- 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.)
- Pending
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention provides a double-scale calculation method for tube sheet strength. Based on the homogenization theory, a tube bundle area is equivalent to a uniform and homogenous area, in an equivalent process, a tube sheet is regarded as a fiber enhanced material, a tube nozzle is regarded as fiber, and the attribute of the equivalent material is obtained according to the essential attributes of the tube sheet and a heat exchange tube material with the Mori-Tanaka method; the tube sheet after equivalence is simulated by finite element analysis software, and a stress moire pattern of the tube sheet equivalent area is obtained; a position with the largest stress is found, tube bundle restoration is performed in the position, and an obtained stress value is the maximum stress value of a tube distribution area in a real model; similarly, maximum stress positions of a non-tube distribution area and a junctions of the tube distribution area and the tube distribution area are determined in part of restoration models respectively; and finally, linear processing is performed on the maximum stress positions of the areas, and strength check is performed. The double-scale calculation method can reduce the hardware condition required by simulation analysis greatly, and shorten time required for simulation analysis of the tube sheet strength effectively.
Description
technical field:
The present invention relates to computing method technical field, refer in particular to a kind of two dimension calculation methods of tubesheet intensity.
background technology:
The sunykatuib analysis of prior art to tubesheet intensity, need to divide the grid of sufficient amount.A large amount of grids is by the process that causes the sunykatuib analysis of tubesheet intensity to be calculated, and calculated amount is excessive, needs to use high performance computer workstation to simulate, and result of calculation cannot be verified in common computer.For the tube sheet that has a large amount of mouths of pipe, number of nodes can reach up to ten million, only spends the time on calculating just to need a month even some months.
summary of the invention:
The object of this invention is to provide a kind of two dimension calculation methods of tubesheet intensity, solve grid in tubesheet intensity simulation process and divide the huge problem of calculated amount of too much bringing.
The object of the invention is to realize as follows: a kind of two dimension calculation methods of tubesheet intensity, this kind of method is equivalent to even same sex region based on homogenization theory by tube bank region, in the process of equivalence, tube sheet is considered as to fiber reinforcement shaped material, the mouth of pipe is considered as fiber, according to the base attribute of tube sheet and heat exchange tube material, by Mori-Tanaka method, obtain the attribute of equivalent material; Tube sheet after utilizing finite element analysis software to equivalence is simulated, and obtains the stress moire pattern in tube sheet equivalence region; Find the position of stress maximum, this position is stringing regional stress maximum in realistic model; Then restrain reduction in this position, the stress value obtaining is the maximum stress value in cloth area under control in realistic model, in like manner, determines respectively the position of non-tube distributing area, cloth area under control and non-tube distributing area intersection maximum stress in partial reduction model; Finally the maximum stress place in each region is carried out to linearization process, carry out strength check according to design standards.
The invention has the beneficial effects as follows, the method can reduce the required hardware condition of sunykatuib analysis greatly, effectively shortens the needed time of tubesheet intensity sunykatuib analysis.
embodiment:
A kind of two dimension calculation methods of tubesheet intensity, it is characterized in that: this kind of method is equivalent to even same sex region based on homogenization theory by tube bank region, in the process of equivalence, tube sheet is considered as to fiber reinforcement shaped material, the mouth of pipe is considered as fiber, according to the base attribute of tube sheet and heat exchange tube material, by Mori-Tanaka method, obtain the attribute of equivalent material; Tube sheet after utilizing finite element analysis software to equivalence is simulated, and obtains the stress moire pattern in tube sheet equivalence region; Find the position of stress maximum, this position is stringing regional stress maximum in realistic model; Then restrain reduction in this position, the stress value obtaining is the maximum stress value in cloth area under control in realistic model, in like manner, determines respectively the position of non-tube distributing area, cloth area under control and non-tube distributing area intersection maximum stress in partial reduction model; Finally the maximum stress place in each region is carried out to linearization process, carry out strength check according to design standards.
Embodiment one:
Tube sheet is equivalent in the process in even same sex region, first hollow pipe is equivalent to even solid pipe, then the region at pore place is equivalent to homogeneous area, finally tube sheet is equivalent to the region with even physical property.The compound substance that tube sheet is considered as, take tube sheet as matrix, heat exchanger tube is fiber.Heat exchanger tube is equivalent to even solid tubes, and take air as matrix, tube wall is fiber, utilizes mixed formulation to obtain the physical parameter of even solid tubes.The physical property of connecting pipe plate material, obtains the material property of the tube sheet after equivalence.
For the tube sheet that contains some mouths of pipe, for reducing the impact of boundary condition on equivalent region, the region of choosing the outer certain limit of tube sheet is that border applies region.
In fact pressure is applied to above the tube sheet of removing after pore, and the region of exerting pressure in simulation is different.Therefore according to area relationship, pressure is carried out to corresponding equivalence.In the process of dividing at grid, in order to improve mesh quality, take different yardsticks to divide to zones of different, most of region adopts hexahedral element, and small part is because geometric configuration reason is taked tetrahedron element, grid division.
According to the physical parameter of calculated equivalent back tube sheet, stress analysis is carried out in the region after utilizing finite element software to equivalence, obtains the stress envelope on tube sheet.
According to the stress distribution of equivalent rear region, can obtain the region that stress ratio is concentrated, these regions are the special region of being concerned about in tube sheet Stress calculation process, local reduction is carried out in these regions, true physical property according to these region material is carried out stress simulation analysis, thereby obtains the maximum stress value of tube sheet in realistic model.Along tube plate thickness direction, path is set thereon, carries out Stress Linearization processing, judge in conjunction with design standards whether intensity meets standard.
The related method of this invention can shorten greatly carries out the needed time of sunykatuib analysis to tubesheet intensity, utilize this kind of method can in common computer, carry out the tubesheet intensity analysis of labyrinth, the time that tubesheet intensity is carried out to analog computation can reach common simulation method tens of/mono-, even one of percentage.
Claims (1)
1. the two dimension calculation methods of tubesheet intensity, it is characterized in that: this kind of method is equivalent to even same sex region based on homogenization theory by tube bank region, in the process of equivalence, tube sheet is considered as to fiber reinforcement shaped material, the mouth of pipe is considered as fiber, according to the base attribute of tube sheet and heat exchange tube material, by Mori-Tanaka method, obtain the attribute of equivalent material; Tube sheet after utilizing finite element analysis software to equivalence is simulated, and obtains the stress moire pattern in tube sheet equivalence region; Find the position of stress maximum, this position is stringing regional stress maximum in realistic model; Then restrain reduction in this position, the stress value obtaining is the maximum stress value in cloth area under control in realistic model, in like manner, determines respectively the position of non-tube distributing area, cloth area under control and non-tube distributing area intersection maximum stress in partial reduction model; Finally the maximum stress place in each region is carried out to linearization process, carry out strength check according to design standards.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410050862.5A CN103778344A (en) | 2014-02-14 | 2014-02-14 | Double-scale calculation method for tube sheet strength |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410050862.5A CN103778344A (en) | 2014-02-14 | 2014-02-14 | Double-scale calculation method for tube sheet strength |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN103778344A true CN103778344A (en) | 2014-05-07 |
Family
ID=50570572
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410050862.5A Pending CN103778344A (en) | 2014-02-14 | 2014-02-14 | Double-scale calculation method for tube sheet strength |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103778344A (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102419121A (en) * | 2011-12-09 | 2012-04-18 | 中国船舶重工集团公司第七�三研究所 | Heat exchange unit of heat exchange device |
| US20130146810A1 (en) * | 2011-12-08 | 2013-06-13 | Basf Se | Process for Producing Water-Absorbing Polymer Fibres |
-
2014
- 2014-02-14 CN CN201410050862.5A patent/CN103778344A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130146810A1 (en) * | 2011-12-08 | 2013-06-13 | Basf Se | Process for Producing Water-Absorbing Polymer Fibres |
| CN102419121A (en) * | 2011-12-09 | 2012-04-18 | 中国船舶重工集团公司第七�三研究所 | Heat exchange unit of heat exchange device |
Non-Patent Citations (2)
| Title |
|---|
| 刘俊明 等: "高压厚管板的有限元分析计算", 《压力容器》 * |
| 崔鼎: "基于细观力学的管板式冷凝器结构分析", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102902847B (en) | Transmission grid single line diagram automatic design method based on combination primitives | |
| CN104281747A (en) | Fine mask screening process analysis method | |
| CN105307397B (en) | A kind of circuit board exposure method and device | |
| CN110889166B (en) | Lightweight design method for aviation bearing bracket based on thermal coupling constraint | |
| CN105489102B (en) | A kind of three-dimensional interactive training training system | |
| CN106650085A (en) | Finite element mesh-based method for generating hollow model by entity model | |
| CN105095086A (en) | Improved software test model and improved software test method | |
| CN107194036A (en) | A kind of large-span roof structure non-gaussian Numerical Simulation Methods of Wind Load | |
| CN107422376B (en) | A kind of non-stationary non-gaussian earthquake motion time history analogy method | |
| CN108595800A (en) | A kind of latticed shell structure modeling method | |
| CN109190317B (en) | Method for quickly and accurately creating grid BIM (building information modeling) model by using Dynamo | |
| CN104827575A (en) | Bidirectional migration method for two-dimensional random feeding of concrete aggregates | |
| CN104794309A (en) | Farming machine design method based on emulation technique | |
| CN103778344A (en) | Double-scale calculation method for tube sheet strength | |
| CN106650092A (en) | Design and construction collaboration platform system architecture | |
| CN109101671B (en) | Variable density and variable configuration three-dimensional lattice structure modeling method | |
| JP2017516130A5 (en) | ||
| CN110334450A (en) | A kind of multi-blocked structure grid generate in object plane erroneous projection restorative procedure | |
| JP2013206996A5 (en) | ||
| CN108280265A (en) | Distortion of the mesh interpolation method, system and medium based on novel radial basis function | |
| CN104732589A (en) | Rapid generation method of hybrid grid | |
| CN103077554B (en) | Based on the crop Interactive Design method and system of joint unit | |
| WO2023287362A3 (en) | Method and apparatus for managing transaction flow of virtual power plant, and computer device and storage medium thereof | |
| CN109101736A (en) | A kind of entire compensation method for high-tension plate | |
| CN104268337A (en) | Method and system for generating path of simulation support platform image connecting line |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20140507 |