CN104794277A - Rubber block temperature field analogue simulation method considering cord thread - Google Patents
Rubber block temperature field analogue simulation method considering cord thread Download PDFInfo
- Publication number
- CN104794277A CN104794277A CN201510187662.9A CN201510187662A CN104794277A CN 104794277 A CN104794277 A CN 104794277A CN 201510187662 A CN201510187662 A CN 201510187662A CN 104794277 A CN104794277 A CN 104794277A
- Authority
- CN
- China
- Prior art keywords
- block rubber
- node
- cord
- energy density
- temperature field
- 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.)
- Granted
Links
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention provides a rubber block temperature field analogue simulation method considering a cord thread. The method is characterized by comprising the following steps that 1, a three-dimensional geometric model of a rubber block is built; 2, a mesh model is built; 3, a rubber block mesh model with the cord thread is built; 4, the unit strain energy density of the rubber block is extracted; 5, the strain energy density of rubber block joints is computed accurately; 6, the joint heat generation rate is computed; 7, temperature field analysis is carried out on the rubber block. The rubber block temperature field analogue simulation method considering the cord thread carries out temperature field computation and analysis on the rubber block model with the joints as the foundation, so that an accurate analogue simulation result is obtained, and a theoretical guidance is provided for actual application of the rubber block.
Description
Technical field
The present invention relates to analogue simulation field, particularly a kind of heat-dissipating of block rubber and the simulation method of heat transfer model.
Background technology
Block rubber parts meet the requirement of shaft coupling properties because of its damping capacity and elastic performance, become the ideal material of shaft coupling.But the good damping capacity of block rubber, makes its Conversion of Energy produced by dynamic torque at work be a large amount of heat power consumptions, if these heats can not be delivered in environment in time fully, then impact in various degree can be caused on block rubber parts.Light then component capabilities is declined, heavy then the time can be caused to damage.Therefore, the thermal property of block rubber parts has very large researching value.
The mechanism of block rubber parts heat-dissipating at work, only have comparatively general heat power consumption at present, namely kinetic transformation is the energy conversion equation of interior energy.Then few for the research refineing to block rubber parts heat-dissipating heat dissipation model operationally and Temperature Distribution.
Summary of the invention
The present invention be directed to that above-mentioned problem carries out, object is to provide that a kind of algorithm is simple, structure accurate block rubber temperature field simulation analytical approach.
The present invention for achieving the above object, have employed following technical scheme:
The invention provides a kind of simulation method taking into account the block rubber temperature field of cord, it is characterized in that, comprise the following steps: step one, set up the 3-D geometric model of block rubber; Step 2, sets up grid model, carries out stress and strain model, and by the mesh refinement 3 ~ 5 times in region larger for deflection in 3-D geometric model, obtain grid model to the 3-D geometric model that step one is set up; Step 3, sets up the block rubber grid model of band cord, is imported by grid model in ABAQUS software, rebar unit is adopted to generate cord, first cord is defined on shell, then the shell with cord is embedded in grid model, obtain the block rubber grid model being with cord; Step 4, extract the unit strain energy density of block rubber, under given operating mode, calculate the stress-strain field distribution situation of the block rubber grid model of band cord, obtain destination file, with Python programming, from destination file, extract unit strain energy density and corresponding unit and node serial number; Step 5, the strain energy density of accurate Calculation block rubber node, first the unit be connected with a node is found out, the number of the unit that statistics is connected, then by be connected unit strain energy density be added, finally by summation divided by unit number, obtain the strain energy density of node, adopt the strain energy density of identical other all nodes of method program calculation, computing formula is as follows:
in formula, W
ifor the strain energy density of i node, W
eifor the strain energy density of unit be connected with i node, n
eifor the unit sum be connected with i node; Step 6, computing node heat generation rate, the heat generation rate formula of block rubber is: q
i=2 π f β W
i, q in formula
ifor the heat generation rate of block rubber i node, f is loading frequency, and β is the dissipation factor of rubber, W
ifor the strain energy density of i node, calculate the heat generation rate of each node, then the block rubber grid model of band cord is imported ANSYS software, utilize APDL Programming with Pascal Language to apply heat generation rate to each node; Step 7, carries out temperature field analysis to block rubber, applies temperature boundaries condition, the block rubber grid model of ANSYS software to band cord is adopted to carry out temperature field analysis, wherein, according to the degree of mesh refinement in the calculated case determining step two of step 4, until the numerical convergence of step 4.
Further, according to the simulation method taking into account the block rubber temperature field of cord provided by the present invention, such feature can also be had, also comprise: step 8, Rationality Assessment is carried out to temperature field analysis result, if surface temperature and the environment temperature of the block rubber grid model of band cord differ greatly, then show to consider the impact of heat radiation on block rubber, revise the temperature boundaries condition in step 7, apply heat emissivity coefficient, temperature field analysis is re-started to the block rubber grid model of band cord.
The effect of invention and effect
According to the simulation method taking into account the block rubber temperature field of cord provided by the present invention, because utilize the block rubber grid model of ABAQUS software to band cord to carry out stress-strain analysis, be extracted unit strain energy density and node serial number, employing is shared algorithm equally and is calculated node strain energy density, finally the block rubber grid model of band cord and node strain energy density are imported the temperature field analysis that ANSYS carries out block rubber model, therefore unit strain energy density can be converted into node strain energy density by the method, based on node, Temperature calculating and analysis are carried out to block rubber model, thus obtain accurate simulation result, for the practical application of block rubber provides theoretical direction, the method has significant scientific research value and economic benefit.
Accompanying drawing explanation
Fig. 1 is the calculation flow chart of the simulation method in the block rubber temperature field taking into account cord;
Fig. 2 is the structural representation of coupling with rubber pad;
Fig. 3 is the grid model figure of coupling with rubber pad;
Fig. 4 is cord cross-sectional view;
Fig. 5 is cord ply grid chart;
Fig. 6 is 1/6th block rubber grid model figure of band cord;
Fig. 7 is the strain energy density figure of 1/6th block rubber grid models of band cord;
Fig. 8 is the temperature field analysis illustraton of model of 1/6th block rubber grid models of band cord;
Fig. 9 is the heat generation rate distribution plan of 1/6th block rubber grid models of band cord;
Figure 10 is the bulk temperature field pattern of 1/6th block rubber grid models of band cord;
Figure 11 is the hull-skin temperature field pattern of 1/6th block rubber grid models of band cord;
Figure 12 is the longitudinal profile thermo parameters method figure of 1/6th block rubber grid models of band cord
Figure 13 is the transverse section thermo parameters method figure of 1/6th block rubber grid models of band cord.
Embodiment
Below in conjunction with accompanying drawing, the simulation method taking into account the block rubber temperature field of cord provided by the present invention is elaborated.
< embodiment >
Fig. 1 is the calculation flow chart of the simulation method in the block rubber temperature field taking into account cord.
As shown in Figure 1, the simulation method taking into account the block rubber temperature field of cord comprises the following steps:
Step S-1, sets up the 3-D geometric model of block rubber.
Fig. 2 is the structural representation of coupling with rubber pad.
Pro/E software is utilized to set up the 3-D geometric model of block rubber as shown in Figure 2.
Step S-2, sets up grid model.
Fig. 3 is the grid model figure of coupling with rubber pad.
The 3-D geometric model that step S-1 sets up is imported to hypermesh software and carries out stress and strain model.Consider that block rubber, under high pulling torque effect, larger geometry deformation can occur, during grid division, answer the distribution of reasonable Arrangement grid, in the larger region of distortion by mesh refinement 3 ~ 5 times, with reduce grid stressed time relative deformation, ensure the accuracy of result.Obtain grid model as shown in Figure 3.
Step S-3, sets up the block rubber grid model of band cord.
Fig. 4 is cord cross-sectional view;
Fig. 5 is cord ply grid chart;
Fig. 6 is 1/6th block rubber grid model figure of band cord.
The grid model that step S-2 generates is preserved into " .inp " file, then imports in ABAQUS software, to rubber and metal parts definition material characteristic, add Cyclic Symmetry constraint, set up the coupling constraint on reference point and the coupling with rubber pad face of cylinder.Staff cultivation is applied to one of them Coupling point, another Coupling point is decontroled to the degree of freedom pivoted.By cord model as shown in Figure 4, import to grid division in ABAQUS, obtain grid chart as shown in Figure 5.Because coupling with rubber pad is symmetrical structure, so carry out sunykatuib analysis to its 1/6th model.When inserting cord in block rubber, first cord is defined on shell, then the shell with cord is embedded in grid model, obtain the block rubber grid model of band cord as shown in Figure 6.
Step S-4, extracts the unit strain energy density of block rubber.
Fig. 7 is the strain energy density figure of 1/6th block rubber grid models of band cord.
Under given operating mode, calculate the stress-strain field distribution situation of the block rubber grid model of band cord, obtain destination file and unit strain energy density figure as shown in Figure 7.With Python programming, from " .odb " file, extract unit strain energy density, and corresponding unit and node serial number.
Step S-5, the strain energy density of accurate Calculation block rubber node.
By the method computing node strain energy density of programming in MATLAB.For some nodes, first find out the unit be connected with this node, the number of unit that statistics is connected, then the strain energy density of the unit be connected be added, finally by summation divided by unit number, obtain the strain energy density of this node.
Computing formula is as follows:
In formula, W
ifor the strain energy density of i node, W
eifor the strain energy density of unit be connected with i node, n
eifor the unit sum be connected with i node.
Said method is adopted to calculate the strain energy density of other all nodes.
Step S-6, computing node heat generation rate.
The heat generation rate formula deriving block rubber according to the heat-dissipating principle of block rubber is:
q
i＝2πfβW
i，
In formula, q
ifor the heat generation rate of block rubber i node, f is loading frequency, and β is the dissipation factor of rubber, W
ifor the strain energy density of i node.
Calculate the heat generation rate of each node, then 1/6th block rubber grid models of band cord are imported ANSYS software, utilize APDL Programming with Pascal Language to apply heat generation rate to each node.
Step S-7, carries out temperature field analysis to block rubber.
Fig. 8 is the temperature field analysis illustraton of model of 1/6th block rubber grid models of band cord.
Definition material attribute in hypermesh, saves as " .cdb " file, imports in ANSYS, obtains temperature field analysis model as shown in Figure 8.
Block rubber parts have good damping capacity, suitable vibration moment of torsion can be born, but simultaneously because damping heat-dissipating makes its internal temperature raise, because carry-over moment itself does not produce power consumption, only have vibration moment of torsion can cause rubber elastic element heat-dissipating, so in realistic model, only temperature field analysis is carried out to the bonding flange on block rubber parts and both sides thereof, ignore the parts such as diaphragm unit, end plate and pressing plate.
Fig. 9 is the heat generation rate distribution plan of 1/6th block rubber grid models of band cord.
Temperature field analysis is carried out to block rubber, heat-dissipating and heat dissipation balancing be considered.Specifically comprise following two parts:
(1) heat-dissipating: be transformed by the viscoelasticity loss of block rubber, to node applying heat generation rate as shown in Figure 9;
(2) dispel the heat: inner at rubber and rubber and steel part surface of contact exists heat transfer, to rubber and steel part definition thermal conductivity; There is thermal convection in the surface of contact of rubber and air, to contact show to define convection transfer rate to rubber with air; Block rubber entirety produces heat radiation to external world, and because block rubber is surperficial and circumstance of temperature difference is little, thermal-radiating impact is very little, therefore ignores.
After applying above-mentioned temperature boundaries condition, ANSYS software is adopted to carry out temperature field analysis to block rubber temperature field analysis model as shown in Figure 8.
Step S-8, analysis result.
Figure 10 is the bulk temperature field pattern of 1/6th block rubber grid models of band cord;
Figure 11 is the hull-skin temperature field pattern of 1/6th block rubber grid models of band cord;
Figure 12 is the longitudinal profile thermo parameters method figure of 1/6th block rubber grid models of band cord
Figure 13 is the transverse section thermo parameters method figure of 1/6th block rubber grid models of band cord.
Temperature field analysis result, as shown in Figure 10 ~ 13, can learn block rubber thermo parameters method situation that is overall and surperficial and each section from analysis result.
According to thermo parameters method figure, sunykatuib analysis result is evaluated, if block rubber model surface temperature and environment temperature differ greatly, then show to consider the impact of heat radiation on block rubber.Now should revise the boundary condition, and the heat emissivity coefficient of definition block rubber, re-starts temperature field analysis to temperature field analysis model.
The effect of embodiment and effect
According to the simulation method taking into account the block rubber temperature field of cord of the present embodiment, because utilize the block rubber grid model of ABAQUS software to band cord to carry out stress-strain analysis, be extracted unit strain energy density and node serial number, employing is shared algorithm equally and is calculated node strain energy density, finally the block rubber grid model of band cord and node strain energy density are imported the temperature field analysis that ANSYS carries out block rubber model, therefore unit strain energy density can be converted into node strain energy density by the method, based on node, Temperature calculating and analysis are carried out to block rubber model, thus obtain accurate simulation result, for the practical application of block rubber provides theoretical direction, the method has significant scientific research value and economic benefit.
In addition, the method carries out analysis on its rationality to analysis result after temperature field analysis, can the error that causes of the condition of revising the boundary in time, ensure analysis result rationally and accurately.
The simulation method taking into account the block rubber temperature field of cord provided by the present invention, is not merely defined in the content described in above embodiment.These are only the present invention conceive under basic explanation, and according to any equivalent transformation that technical scheme of the present invention is done, all should protection scope of the present invention be belonged to.
In addition, in the present embodiment in step 2 by the mesh refinement 3 ~ 5 times in region larger for deflection, for other models, can first by distortion compared with large regions mesh refinement to a certain degree, observe the result of calculation of step 4, step 4 if can not restrain, then increases mesh refinement degree, until can obtain the result of calculation restrained.
Claims (2)
1. take into account the simulation method in the block rubber temperature field of cord, it is characterized in that, comprise the following steps:
Step one, sets up the 3-D geometric model of block rubber;
Step 2, sets up grid model, carries out stress and strain model, and by the mesh refinement in region larger for deflection in described 3-D geometric model, obtain grid model to the 3-D geometric model that described step one is set up;
Step 3, set up the block rubber grid model of band cord, described grid model is imported in ABAQUS software, rebar unit is adopted to generate cord, first described cord is defined on shell, again the described shell with described cord is embedded in described grid model, obtain the block rubber grid model being with cord;
Step 4, extract the unit strain energy density of block rubber, under given operating mode, calculate the stress-strain field distribution situation of the block rubber grid model of described band cord, obtain destination file, with Python programming, from described destination file, extract unit strain energy density and corresponding unit and node serial number;
Step 5, the strain energy density of accurate Calculation block rubber node, first the unit be connected with a node is found out, add up the number of described connected unit, then by the strain energy density of described connected unit be added, finally by summation divided by unit number, obtain the strain energy density of described node, adopt the strain energy density of identical other all nodes of method program calculation
Computing formula is as follows:
In formula, W
ifor the strain energy density of i node, W
eifor the strain energy density of unit be connected with i node, n
eifor the unit sum be connected with i node;
Step 6, computing node heat generation rate, the heat generation rate formula of block rubber is: q
i=2 π f β W
i,
Q in formula
ifor the heat generation rate of block rubber i node, f is loading frequency, and β is the dissipation factor of rubber, W
ifor the strain energy density of i node, calculate the heat generation rate of each described node, then the block rubber grid model of described band cord is imported ANSYS software, utilize APDL Programming with Pascal Language to apply described heat generation rate to each described node;
Step 7, carries out temperature field analysis to block rubber, applies temperature boundaries condition, adopts the block rubber grid model of ANSYS software to described band cord to carry out temperature field analysis,
Wherein, determine the degree of mesh refinement in described step 2 according to the calculated case of described step 4, until the numerical convergence of described step 4.
2. the simulation method taking into account the block rubber temperature field of cord according to claim 1, is characterized in that, also comprise:
Step 8, Rationality Assessment is carried out to temperature field analysis result, if surface temperature and the environment temperature of the block rubber grid model of described band cord differ greatly, then show to consider the impact of heat radiation on block rubber, revise the temperature boundaries condition in described step 7, apply heat emissivity coefficient, temperature field analysis is re-started to the block rubber grid model of described band cord.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510187662.9A CN104794277B (en) | 2015-04-20 | 2015-04-20 | The simulation method in the block rubber temperature field of meter and cord |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510187662.9A CN104794277B (en) | 2015-04-20 | 2015-04-20 | The simulation method in the block rubber temperature field of meter and cord |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104794277A true CN104794277A (en) | 2015-07-22 |
| CN104794277B CN104794277B (en) | 2017-09-12 |
Family
ID=53559068
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510187662.9A Active CN104794277B (en) | 2015-04-20 | 2015-04-20 | The simulation method in the block rubber temperature field of meter and cord |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104794277B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110955999A (en) * | 2019-11-29 | 2020-04-03 | 青岛森麒麟轮胎股份有限公司 | Tire three-dimensional temperature field simulation method |
| CN113435001A (en) * | 2021-05-10 | 2021-09-24 | 重庆大学 | Method for predicting meshing temperature field of plastic gear |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB261808A (en) * | 1925-06-02 | 1926-12-02 | Cecil Featherstone Hammond | Improvements in and connected with gaseous fuel burners |
| CN103880455A (en) * | 2014-02-12 | 2014-06-25 | 嘉兴皓特特种陶瓷有限公司 | Preparation method of steel cord enhanced graphite silicon carbide crucible |
-
2015
- 2015-04-20 CN CN201510187662.9A patent/CN104794277B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB261808A (en) * | 1925-06-02 | 1926-12-02 | Cecil Featherstone Hammond | Improvements in and connected with gaseous fuel burners |
| CN103880455A (en) * | 2014-02-12 | 2014-06-25 | 嘉兴皓特特种陶瓷有限公司 | Preparation method of steel cord enhanced graphite silicon carbide crucible |
Non-Patent Citations (2)
| Title |
|---|
| 张昊 等: "全钢子午线轮胎应力和生热研究", 《轮胎工业》 * |
| 王延顼 等: "采用计算机模拟预测滚动轮胎的温度分布状况", 《现代橡胶技术》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110955999A (en) * | 2019-11-29 | 2020-04-03 | 青岛森麒麟轮胎股份有限公司 | Tire three-dimensional temperature field simulation method |
| CN113435001A (en) * | 2021-05-10 | 2021-09-24 | 重庆大学 | Method for predicting meshing temperature field of plastic gear |
| CN113435001B (en) * | 2021-05-10 | 2022-09-06 | 重庆大学 | Method for predicting meshing temperature field of plastic gear |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104794277B (en) | 2017-09-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108984920B (en) | Direct fluid-solid coupling heat transfer analysis method for engine cooling water jacket | |
| CN103279627B (en) | A kind of heat-machinery based on finite element-wearing and tearing coupling analysis method for numerical simulation | |
| CN106626449B (en) | Consider the design method in the composite material V-type component autoclave molding tooling type face of curing deformation | |
| Dang et al. | Accelerated loosely-coupled CFD/CSD method for nonlinear static aeroelasticity analysis | |
| CN102855347A (en) | Finite-element analysis method for thermal mechanical strength of exhaust manifold | |
| CN112580236B (en) | Rapid analysis method for nonlinear dynamic response of thermal protection connection structure | |
| Peng et al. | Numerical investigation of the effects of structural geometric and material nonlinearities on limit-cycle oscillation of a cropped delta wing | |
| CN115618503B (en) | Method for simulating additive process and optimizing process of rudder wing structure | |
| Huang et al. | Coupling finite difference method with finite particle method for modeling viscous incompressible flows | |
| CN104239625A (en) | Corrective fluid motion equation linear iteration-based steady state solution method | |
| CN104794277A (en) | Rubber block temperature field analogue simulation method considering cord thread | |
| CN106980734A (en) | A kind of RPV sealing property computational methods for considering material and contact nonlinear | |
| Wu et al. | A novel dynamic isogeometric reanalysis method and its application in closed-loop optimization problems | |
| Meng et al. | Modelling of anisotropic beam for rotating composite wind turbine blade by using finite-difference time-domain (FDTD) method | |
| CN105404751A (en) | Method for realizing network unification based on heat-force-electromagnetic field | |
| Xue et al. | Efficient transient thermal analysis based on spectral element time domain method with curvilinear hexahedrons | |
| CN104699912A (en) | Strength calculation method of hub and spindle connecting bolt of wind generator set | |
| CN109657368A (en) | A kind of optimization method of air conditioner heat radiator | |
| CN104573206A (en) | Profile steel section thermal shape and size design method based on finite element thermal coupling | |
| CN113591416B (en) | Thermosetting coupling analysis method suitable for steam turbine She Pianliu | |
| Lu et al. | High order large eddy simulations for a transonic turbine blade using hybrid unstructured meshes | |
| Zhang et al. | A high-order flux reconstruction/correction procedure via reconstruction method for shock capturing with space-time extension time stepping and adaptive mesh refinement | |
| CN112364552A (en) | High-pressure cylinder dynamic thermal stress analysis method based on finite element | |
| Hu et al. | An o (n) complexity recursive algorithm for multi‐flexible‐body dynamics based on absolute nodal coordinate formulation | |
| Başkut et al. | Development of a closely coupled procedure for dynamic aeroelastic analyses |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |