A thermodynamic modeling of rubber was proposed for the purpose of evaluating long term thermal aging stability of rubber seals. The finite element analysis code that was developed based on the proposed formulations, apply the hyperelastic constitutive equation to mechanical property, un-stationary heat conduction equation to thermal property and the chemical kinetics theory to chemical reaction. Approximately 20 numerical parameters that required for the finite element analysis code, were determined by material test results and numerical simulations were conducted for a rubber seal in supposed kinematic geometries ; therefore, configurations of chemical reaction caused by thermal environments and softening phenomena of mechanical stiffness were evaluated easier than former methods.