Patent ID: 7451067

Claim:
A method for analysis of a cell structure for analyzing stress which occurs inside the cell structure due to a temperature distribution which occurs inside the cell structure or pressure applied from outside of the cell structure, the cell structure being in a shape of a tubular body including two end faces and a body face which connects the two end faces, a plurality of cells partitioned by wall sections being formed inside the tubular body in an axial direction of the tubular body, and the cell structure including a repeating structure, formed by the wall sections and the cells, the cell structure being a honeycomb structure, and the cells are parallel passages for fluids, the method comprising: a macroanalysis step which includes replacing the cell structure or a part of the cell structure with an anisotropic solid body having property values of equivalent rigidity characteristics, creating a finite element model of the anisotropic solid body based on the property values, applying an internal temperature distribution or an external pressure to the finite element model of the anisotropic solid body, calculating the stress, and selecting an aimed cell structure of which the stress should be calculated as a cell structure based on a stress distribution in the anisotropic solid body; and a microanalysis step which includes creating a finite element model of the aimed cell structure, and calculating the stress based on the finite element model of the aimed cell structure to determine a stress distribution in the aimed cell structure, wherein the rigidity characteristics of the anisotropic solid body in the macroanalysis step are expressed by the following numerical equation (1): ( σ ⁢ ⁢ x σ ⁢ ⁢ y σ ⁢ ⁢ z τ ⁢ ⁢ xy τ ⁢ ⁢ yz τ ⁢ ⁢ zx ) = ( K ⁢ ⁢ 11 K ⁢ ⁢ 12 K ⁢ ⁢ 13 0 0 0 K ⁢ ⁢ 21 K ⁢ ⁢ 22 K ⁢ ⁢ 23 0 0 0 K ⁢ ⁢ 31 K ⁢ ⁢ 32 K ⁢ ⁢ 33 0 0 0 0 0 0 K ⁢ ⁢ 44 0 0 0 0 0 0 K ⁢ ⁢ 55 0 0 0 0 0 0 K ⁢ ⁢ 66 ) ⁢ ( ɛ ⁢ ⁢ x ɛ ⁢ ⁢ y ɛ ⁢ ⁢ z γ ⁢ ⁢ xy γ ⁢ ⁢ yz γ ⁢ ⁢ zx ) ( 1 ) σx: X-axis direction normal stress, σy: Y-axis direction normal stress, σz: Z-axis direction (honeycomb passage direction) normal stress, τxy: Y-axis direction shear stress in a plane perpendicular to the X axis, τyz: Z-axis direction shear stress in a plane perpendicular to the Y axis, τzx: X-axis direction shear stress in a plane perpendicular to the Z axis, εx: X-axis direction tensile (or compression) strain, εy: Y-axis direction tensile (or compression) strain, εz: Z-axis direction tensile (or compression) strain, γxy: XY inplane shear strain, γyz: YZ inplane shear strain, γzx: ZX inplane shear strain, K11, K12, K13, K21, K22, K23, K31, K32, K33, K44, K55, and K66: moduli of elasticity.