There exist no large-series human data linking contact stress exposure to an articular joint's propensity for developing osteoarthritis because contact stress analysis for large numbers of subjects remains impractical. The speed and simplicity of discrete element analysis (DEA) for estimating contact stresses makes its application to this problem highly attractive, but to date DEA has been used to study only a small numbers of cases. This is because substantial issues regarding its use in population-wide studies have not been addressed. Chief among them are developing fast and robust methods for model derivation and the selection of boundary conditions, establishing accuracy of computed contact stresses, and including capabilities for modeling in-series structural elements (e.g., a meniscus). This article describes an implementation of DEA that makes it feasible to perform subject-specific modeling in articular joints in large population-based studies.
Keywords computational joint modeling; osteoarthritis; contact stressThe development of expeditious methods suitable for determining in vivo articular contact stress distributions is critical to understanding normal articular function and the mechanopathology of osteoarthritis (OA). Excessive functional loads habitually applied to at-risk joints (shoulder, wrist, hip, knee, ankle) can negatively impact the health of articular cartilage, eventually leading to degeneration (Buckwalter & Brown, 2004). Many computational models have been developed to determine contact stresses in a variety of joints, mostly using finite element analysis (FEA) . For the most part, however, these studies have involved lengthy and concerted efforts in well-controlled research settings, and have yielded articular contact stresses for only small numbers of cases. The high logistical cost of performing 3-D contact FEA on a subject-or patient-specific basis makes FEA of limited attraction for stress analysis in large epidemiologic studies (Englund et al., 2007;Felson & Nevitt, 2004;Lester, 2008). Advancement of expeditious methods for determining in vivo articular contact stress distributions will enhance the understanding of normal articular function, and of the mechano-pathology of OA. Chao and colleagues popularized the use of rigid body spring modeling (now commonly known as discrete element analysis, or DEA) in biomechanics, as a simple numerical framework to model articular contact (An et al., 1990;Genda et al., 2001;Iwasaki et al., 1998;Li et al., 1994;Volokh et al., 2007). The DEA formulation involves treating bones as rigid bodies, and the cartilage at a given joint as an array of compressive-only springs distributed over the articulating bone surfaces. This modeling approach has been used to study contact mechanics in various diarthrodial joints, most recently progressing to fully 3D subject-specific geometry (Elias & Cosgarea, 2007;Elias et al., 2004;Iwasaki et al., 1998). However, it has rarely been used to study large numbers of cases. To our knowledge, the largest 3-D ...