18.04.04 · organismal-bio / musculoskeletal

Skeletal mechanics: bone remodeling, joint biomechanics, and the lever principles

stub3 tiersLean: nonepending prereqs

Anchor (Master): Nigg, B. M. & Herzog, W. — Biomechanics of the Musculo-skeletal System, 3rd ed. (2007)

Intuition Beginner

Bones are living tissue that constantly breaks down and rebuilds. Osteoblasts build new bone; osteoclasts dissolve old bone. This remodeling adapts bone to the loads you place on it (Wolff's law).

Joints connect bones. Some barely move (skull sutures); others — like your knee — move freely, cushioned by cartilage and lubricated by synovial fluid.

Your body works as a system of levers. Bones are rigid bars, joints are fulcrums, and muscles pull on bones to create rotation. Just like a seesaw, the arrangement of fulcrum, effort, and load determines whether you gain force or speed. Most joints sacrifice force for large, fast movements.

Visual Beginner

The three lever classes differ in the arrangement of fulcrum, effort, and load. Most joints in the human body are third-class levers, which sacrifice force for range of motion and speed.

Worked example Beginner

Your biceps inserts about 5 cm from the elbow joint. You hold a 2 kg weight in your hand, 30 cm from the joint. How much force must the biceps produce?

Step 1. The lever principle says effort times its distance equals load times its distance:

Step 2. The load force is .

Step 3. Solve for muscle force:

The biceps must produce six times the load force. This is the price of a third-class lever: the muscle works harder but the hand moves six times farther and faster than the muscle's attachment point.

Check your understanding Beginner

Formal definition Intermediate+

Bone remodeling

Bone is a composite material of hydroxyapatite crystals (a calcium phosphate mineral giving stiffness) embedded in a collagen matrix (giving toughness). Cortical (compact) bone forms the outer shell of long bones; trabecular (spongy) bone fills the ends and vertebrae. Both types are living tissue undergoing continuous renewal.

The basic multicellular unit (BMU) is the functional unit of bone remodeling. A BMU consists of a team of osteoclasts at the leading edge that resorb a tunnel of old bone, followed by osteoblasts that fill the tunnel with new bone. A single BMU progresses through bone at about 20-40 m/day and has a lifespan of about 6-12 months. At any given time, roughly one million BMUs are active throughout the adult skeleton. The remodeling cycle has three phases:

  1. Activation: pre-osteoclasts are recruited to a bone surface in response to microdamage or hormonal signals (parathyroid hormone, PTH, stimulates osteoclast activity; calcitonin inhibits it).
  2. Resorption: osteoclasts adhere to the bone surface, seal off a compartment, and secrete hydrochloric acid and proteolytic enzymes (cathepsin K, tartrate-resistant acid phosphatase) that dissolve mineral and digest collagen. Resorption takes about 1-3 weeks.
  3. Formation: osteoblasts deposit osteoid (unmineralised collagen matrix), which subsequently mineralises over 3-6 months. Osteoblasts that become trapped in the matrix differentiate into osteocytes, which act as mechanosensors.

Wolff's law (1892) states that bone adapts its internal architecture and external shape to the mechanical loads placed upon it. The trabecular struts in the femoral head, for example, align along the principal stress trajectories, forming a pattern that engineers would recognise as an optimised lattice. The cellular mechanism involves osteocytes in the lacuno-canalicular network sensing fluid flow caused by bone strain; these signals modulate osteoblast and osteoclast activity via RANKL/OPG signalling and sclerostin secretion.

The stress-strain relationship for bone is approximately linear (Hookean) in the physiological range. Cortical bone has a Young's modulus of about 15-20 GPa, an ultimate compressive strength of about 130-190 MPa, and an ultimate tensile strength of about 100-140 MPa. Bone fails in tension at lower strains (0.025) than in compression (0.045) because crack propagation proceeds more readily under tensile loading. The collagen matrix provides post-yield toughness by bridging cracks and dissipating energy.

Joint classification

Joints are classified by the tissue that connects the articulating bones:

  • Fibrous joints: bones connected by dense fibrous connective tissue. Examples: cranial sutures (synarthroses, essentially immobile), syndesmoses (interosseous membrane between radius and ulna, allowing slight movement).
  • Cartilaginous joints: bones connected by cartilage. Synchondroses (hyaline cartilage, e.g., costochondral joints, epiphyseal plates before fusion) allow growth but minimal movement. Symphyses (fibrocartilage, e.g., intervertebral discs, pubic symphysis) allow limited movement and absorb compressive loads.
  • Synovial joints: the most mobile and most complex class. The articulating surfaces are covered with hyaline articular cartilage, enclosed within a joint capsule lined by a synovial membrane that secretes synovial fluid. Synovial fluid is a dialysate of plasma enriched with hyaluronic acid, providing both lubrication and nutrient transport to the avascular cartilage.

Synovial joint subtypes are classified by the shape of the articular surfaces and the degrees of freedom:

  • Plane joints: flat surfaces, gliding motion (intercarpal joints).
  • Hinge joints: uniaxial, flexion/extension (elbow, interphalangeal).
  • Pivot joints: uniaxial, rotation (atlantoaxial, proximal radioulnar).
  • Condyloid joints: biaxial, flexion/extension + abduction/adduction (metacarpophalangeal).
  • Saddle joints: biaxial with concave-convex surfaces (carpometacarpal of the thumb).
  • Ball-and-socket joints: triaxial, all rotations (glenohumeral, hip).

Lever mechanics

The musculoskeletal system operates as a system of levers. A lever consists of a rigid bar (bone) that rotates around a fulcrum (joint). Two forces act: the effort (muscle force) and the load (resistance, typically gravity acting on a body segment or an external weight). The moment arm of a force is the perpendicular distance from the fulcrum to the line of action of that force. The torque (moment) produced by a force about the fulcrum is:

where is the force magnitude and is the moment arm. For rotational equilibrium, the sum of clockwise torques equals the sum of counterclockwise torques.

The three lever classes are distinguished by the relative positions of fulcrum (F), effort (E), and load (L):

  • Class 1 (F between E and L): mechanical advantage can be greater than, equal to, or less than one depending on the relative moment arms. Example: the atlanto-occipital joint — the fulcrum is the joint, the effort comes from posterior neck muscles pulling down on the occiput, and the load is the weight of the head anterior to the joint. The head rocks on the vertebral column like a seesaw.
  • Class 2 (L between F and E): mechanical advantage is always greater than one (the effort arm is longer than the load arm). Example: ankle plantarflexion during standing — the fulcrum is the metatarsal heads, the load is body weight acting through the tibia (ankle joint), and the effort is the calf muscles pulling via the Achilles tendon posterior to the ankle. The calf produces less force than the body weight because it acts through a longer lever arm.
  • Class 3 (E between F and L): mechanical advantage is always less than one. Example: elbow flexion — the biceps inserts only a few centimetres distal to the elbow joint (fulcrum), while the hand (load) is 30-35 cm away. The muscle must produce forces far exceeding the external load, but the hand moves through a much larger arc than the muscle's insertion point, providing speed and range of motion.

The mechanical advantage of a lever is the ratio of load to effort, which equals the ratio of the effort moment arm to the load moment arm:

Most muscles in the human body operate at a mechanical disadvantage (MA < 1) because they attach close to the joint. This arrangement trades force for displacement and speed, which is functionally useful: a small contraction of the biceps produces a large movement of the hand.

Key mechanism Intermediate+

Torque and moment arms

The turning effect of a force about a joint depends on both the magnitude of the force and its moment arm. The moment arm is not the anatomical distance from the joint to the muscle's insertion; it is the perpendicular distance from the joint centre to the line of action of the muscle force. As a joint moves through its range of motion, the muscle's line of action changes orientation, so the moment arm varies with joint angle. This variation is one reason muscles produce different torques at different joint positions even at constant activation.

For a single muscle acting across a joint, the torque is:

where is the muscle force (itself a function of length, velocity, and activation — see 18.04.01) and is the moment arm as a function of joint angle . The joint angle at which the moment arm is maximal typically corresponds to the position where the muscle's line of action is most nearly perpendicular to the bone.

When multiple muscles act across the same joint, the net torque is the vector sum of individual muscle torques:

Co-contraction of antagonistic muscles (e.g., biceps and triceps at the elbow) produces torques in opposite directions; the net torque determines whether the joint accelerates in flexion or extension. Co-contraction sacrifices net torque but increases joint stiffness and stability.

Joint reaction forces

Newton's third law requires that the joint surfaces bear the reaction forces transmitted through the articulation. For a simple two-dimensional analysis of elbow flexion holding a weight, the joint reaction force at the elbow can be estimated from static equilibrium:

The biceps force acts upward with moment arm . The weight of the forearm acts downward at the forearm's centre of mass, and the external load acts downward at the hand. For equilibrium:

The joint reaction force (directed laterally and distally along the humerus) is approximately:

in the vertical direction (simplified). In reality, the joint reaction force is a vector with both magnitude and direction, and during activities of daily living the joint reaction force at the hip can reach 2-3 times body weight and at the knee 3-4 times body weight during stair climbing.

Mechanical advantage and the musculoskeletal trade-off

The musculoskeletal system faces a fundamental design constraint: muscles that attach close to the joint (short moment arms) must produce large forces but gain speed and range of motion; muscles that attach far from the joint (long moment arms) gain mechanical advantage but sacrifice excursion. The biceps, with a moment arm of about 4-5 cm at the elbow, must produce roughly 7-8 times the load force when the hand is at full forearm length (about 35 cm). The masseter, by contrast, has a moment arm ratio closer to 1:1, giving it high bite force but limited jaw opening speed.

This trade-off is not a design flaw; it reflects the geometric constraint that the product of force and displacement (work) is conserved. A muscle working at a mechanical disadvantage of 7:1 produces 7 times the force but its insertion moves only 1/7 as far as the hand. The muscle's working range (a few centimetres of contraction) is amplified into a large arc of motion at the hand.

Exercises Intermediate+

Bone mechanics and joint biomechanics Master

Finite element analysis of bone

The mechanical behaviour of bone cannot be fully captured by simple beam theory because bone is anisotropic (different stiffness in different directions), heterogeneous (cortical versus trabecular), and viscoelastic (stiffness depends on loading rate). Finite element analysis (FEA) divides the bone geometry into a mesh of small elements, assigns material properties to each element, applies boundary conditions (loads and constraints), and solves the equations of elasticity numerically to obtain the stress and strain distribution throughout the bone.

Subject-specific FEA models are constructed from quantitative computed tomography (QCT) scans, which provide both geometry and bone density at each voxel. The elastic modulus at each element is computed from the CT Hounsfield units using empirical density-modulus relationships (e.g., for trabecular bone in MPa, where is apparent density in g/cm). These models predict fracture risk in osteoporotic patients with higher specificity than bone mineral density alone, because they account for geometry (bone size, cortical thickness) and loading direction — factors that DXA scans average out.

The validation of FEA against cadaveric fracture testing shows prediction errors of 10-25% for femoral fracture load, depending on model complexity and the accuracy of boundary conditions. The main sources of error are: (i) simplified material models that do not capture post-yield behaviour and damage accumulation; (ii) uncertain muscle and joint contact forces during the loading event (a fall, for example, involves dynamic impact forces that are difficult to reconstruct); (iii) patient-specific variation in bone quality (collagen cross-linking, microdamage accumulation) that is not captured by CT density.

Stress fractures

Stress fractures occur when repetitive sub-threshold loading accumulates microdamage faster than remodeling can repair it. The pathogenesis is a race between microcrack propagation and BMU-mediated repair. Normal bone contains microcracks that accumulate with age and are removed by targeted remodeling (osteocyte apoptosis near a crack triggers local BMU activation). When the loading rate exceeds the remodeling capacity — as in military recruits marching with heavy packs, distance runners ramping mileage too quickly, or ballet dancers training on hard surfaces — microcracks coalesce into a macroscopic fracture.

The most common sites are the tibial diaphysis (medial cortex, compression side), the femoral neck (superior cortex, tension side), the metatarsals (especially the second and third), and the pars interarticularis of the lumbar vertebrae (spondylolysis in adolescent athletes). The mechanical basis differs by site: tibial stress fractures result from repetitive bending loads during running, femoral neck fractures from cyclic compression and shear during weight-bearing, and pars defects from repetitive hyperextension loading in gymnastics and football.

Early stress fractures (bone scan or MRI positive, X-ray negative) heal with activity modification in 6-8 weeks. Advanced fractures with a visible fracture line on X-ray require longer protection. The key clinical insight is that stress fractures are a remodeling disease, not a single-event trauma: they develop over weeks of excessive loading and heal over weeks of relative rest, and the preventive strategy is to control the rate of increase in mechanical loading (the "10% rule" for runners: increase weekly mileage by no more than 10%).

Osteoporosis: DEXA scans and T-scores

Osteoporosis is defined clinically by the World Health Organization as a bone mineral density (BMD) T-score of -2.5 or below at the lumbar spine, total hip, or femoral neck, as measured by dual-energy X-ray absorptiometry (DEXA). The T-score compares the patient's BMD to the mean BMD of a young adult reference population:

T-scores between -1.0 and -2.5 indicate osteopenia; T-scores above -1.0 are normal. The Z-score compares the patient to age-matched controls and is used for younger patients.

The relationship between BMD and fracture risk is approximately exponential: each standard deviation decrease in BMD roughly doubles the fracture risk. However, BMD alone explains only about 60-70% of bone strength; the remainder is determined by bone geometry (cortical thickness, cross-sectional area), microarchitecture (trabecular connectivity, cortical porosity), and material properties (mineralisation, collagen quality). This is why two patients with the same T-score can have very different fracture risks.

The FRAX tool (WHO Fracture Risk Assessment Tool) integrates BMD with clinical risk factors (age, sex, body mass index, previous fracture, parental hip fracture, glucocorticoid use, rheumatoid arthritis, secondary osteoporosis, alcohol, smoking) to estimate 10-year fracture probability. Treatment thresholds vary by country but are generally set at a 10-year hip fracture probability of 3% or major osteoporotic fracture probability of 20%.

Pharmacological treatment targets different stages of the remodeling cycle. Antiresorptives (bisphosphonates, denosumab, selective oestrogen receptor modulators) reduce osteoclast activity or survival, slowing bone loss. Anabolic agents (teriparatide, abaloparatide, romosozumab) stimulate osteoblast activity, building new bone. Bisphosphonates bind to hydroxyapatite and are ingested by osteoclasts, where they inhibit farnesyl pyrophosphate synthase in the mevalonate pathway, preventing the prenylation of small GTPases required for osteoclast function and survival. Denosumab is a monoclonal antibody against RANKL, preventing osteoclast differentiation. Romosozumab is an anti-sclerostin antibody that removes the inhibition of Wnt signalling in osteoblasts, promoting bone formation while also reducing resorption.

Osteoarthritis: cartilage degradation

Osteoarthritis (OA) is the progressive degradation of articular cartilage accompanied by subchondral bone sclerosis, osteophyte formation, synovial inflammation, and joint space narrowing. It is a whole-organ disease of the synovial joint, not merely "wear and tear" of cartilage.

Articular cartilage is avascular, aneural, and alymphatic tissue composed of a sparse population of chondrocytes (about 2% of tissue volume) embedded in an extracellular matrix of type II collagen fibrils and aggrecan (a large proteoglycan that attracts water via its dense negative charge from chondroitin and keratan sulphate chains). The swelling pressure of the hydrated aggrecan, restrained by the collagen fibril network, gives cartilage its compressive stiffness and ability to distribute loads across the joint surface.

In OA, an imbalance between matrix synthesis and degradation shifts the net balance toward loss. Matrix metalloproteinases (MMP-1, MMP-3, MMP-13) and aggrecanases (ADAMTS-4, ADAMTS-5) cleave collagen and aggrecan respectively. The triggers include mechanical overload (obesity, joint malalignment, meniscal or ligament injury), inflammatory mediators (IL-1beta, TNF-alpha, IL-6 from inflamed synovium), and age-related changes in chondrocyte function (reduced synthetic capacity, senescence-associated secretory phenotype).

The biomechanics of OA progression involve a positive feedback loop: cartilage loss reduces the contact area and increases peak contact stress on the remaining cartilage, accelerating further degradation. Subchondral bone stiffens in response to altered loading (a Wolff's law response to the changed stress distribution), which further increases the stress transmitted to the overlying cartilage. Joint malalignment (varus or valgus at the knee) concentrates load on one compartment, creating a focal stress peak that drives compartment-specific OA.

Joint replacement (arthroplasty) replaces the degenerated articulating surfaces with metal and polyethylene (or ceramic) components. The biomechanical design of joint replacements aims to restore normal joint kinematics, minimise polyethylene wear particle generation (which causes osteolysis and implant loosening), and distribute loads to avoid stress shielding of the surrounding bone. Stress shielding occurs when the stiff implant carries load that the bone would normally bear, causing bone resorption (Wolff's law in reverse) and weakening the bone-implant interface.

Ligament mechanics

Ligaments connect bone to bone and resist excessive joint displacement. Their mechanical behaviour is characterised by a nonlinear stress-strain curve with four regions:

  1. Toe region (0-3% strain): collagen fibres straighten from their relaxed crimped state. The stiffness is low and the curve is concave upward.
  2. Linear region (3-8% strain): fibres are taut and bear load in proportion to strain. The slope is the tangent modulus (typically 300-500 MPa for the human anterior cruciate ligament).
  3. Microfailure region (8-12% strain): individual fibre bundles begin to fail progressively. The curve departs from linearity.
  4. Macrofailure (12-15% strain): complete rupture.

Ligaments exhibit viscoelastic behaviour: their stiffness and failure load depend on the rate of loading. At high loading rates (simulating an injury mechanism), the ligament is stiffer and absorbs more energy before failure. This rate-dependence is attributed to the fluid flow within the ground substance and the time-dependent recruitment of collagen fibres.

Creep is the progressive elongation of a ligament under constant load. When a ligament is subjected to a sustained sub-failure load, it elongates over time as collagen fibres rearrange and the ground substance flows. This has clinical relevance: prolonged sitting in a slouched posture creep-elongates the posterior spinal ligaments, reducing their ability to stabilise the spine and contributing to low back pain.

Stress relaxation is the decline in tension over time when a ligament is held at constant elongation. The ligament initially bears high tension but the force decays to about 60-80% of the initial value over 1-2 hours. Both creep and stress relaxation reflect the viscoelastic nature of the collagen-ground substance composite.

Tendon properties

Tendons connect muscle to bone and transmit muscle force to the skeleton. They are composed predominantly of type I collagen (about 85% dry weight) arranged in parallel fascicles. Tendons are stiffer than ligaments (modulus about 1-2 GPa) because their collagen fibres are more highly aligned.

Tendon has a lower strain at failure than muscle (8-10% versus 15-20% for passive muscle), which means that under extreme stretch the tendon is not necessarily the weakest link — the muscle-tendon junction or the bone-tendon insertion (enthesis) may fail first. The enthesis is a functionally graded material that transitions from tendon (collagen-dominated) to fibrocartilage to calcified fibrocartilage to bone, distributing stress across a compliant-to-stiff gradient and avoiding a sharp material discontinuity.

Tendon responds to loading with the same Wolff's-law logic as bone, but more slowly. Chronic exercise increases tendon cross-sectional area and stiffness; immobilisation decreases both. The metabolic activity of tendon is low (avascular in the mid-substance, relying on diffusion from the paratenon or synovial fluid), so adaptation takes months rather than weeks.

Bone healing phases

Fracture healing proceeds through four overlapping phases:

  1. Inflammatory phase (0-7 days): haematoma formation at the fracture site, invasion by inflammatory cells (neutrophils, macrophages), release of cytokines (BMPs, TGF-beta, PDGF, FGF) that recruit mesenchymal stem cells. The haematoma provides the initial scaffold and signalling environment.
  2. Soft callus phase (1-3 weeks): mesenchymal stem cells differentiate into chondrocytes and fibroblasts, producing a cartilaginous (soft) callus that bridges the fracture gap. The callus is initially too weak to bear weight but provides some stability.
  3. Hard callus phase (3-12 weeks): endochondral ossification converts the cartilaginous callus to woven bone. Osteoblasts deposit bone directly (intramembranous ossification) in well-vascularised regions and indirectly via a cartilage intermediate (endochondral ossification) in less vascularised regions. The hard callus is visible on X-ray as a bulge of new bone around the fracture.
  4. Remodeling phase (months to years): woven bone is gradually replaced by lamellar bone through BMU-mediated remodeling. The callus reshapes to approximate the original bone geometry (though a small deformity may persist if the fracture was not well reduced). The final result is bone that is mechanically comparable to the original, though the remodelled cortex may retain a slightly different architecture.

Secondary healing (with callus formation, as described above) requires some motion at the fracture site. Rigid internal fixation (plates and screws) suppresses callus formation and forces primary healing, in which osteoclasts tunnel directly across the fracture line followed by osteoblasts reconstructing Haversian systems. Primary healing is slower for the patient to bear weight but produces a more anatomical reduction.

Connections Master

Skeletal mechanics connects cellular bone biology to whole-body biomechanics. Wolff's law is a feedback loop: mechanical loading generates strain in bone tissue, which osteocytes sense and translate into remodeling signals that adjust bone mass and architecture. The same principle applies at the joint level: cartilage and subchondral bone adapt to habitual loading patterns, and the failure of these adaptations underlies osteoarthritis.

The lever analysis of joints links this unit directly to muscle physiology (18.04.01). Muscle force production (the Hill equation, the length-tension relationship) becomes the input to the torque calculation at each joint. The muscle's operating point on its force-length-velocity surface determines the torque it can produce at a given joint angle and angular velocity, and the lever geometry (moment arm as a function of joint angle) determines how that muscle force translates into joint torque and ultimately into movement.

The joint reaction forces computed from lever analysis are the loads that drive bone remodeling (Wolff's law) and cartilage loading (osteoarthritis risk). The hip joint reaction force of 2-3 times body weight during gait and 4-5 times body weight during running is the mechanical stimulus that maintains (or fails to maintain) femoral neck bone density. Understanding these forces is essential for orthopaedic implant design, fracture risk prediction, and rehabilitation planning.

The clinical dimension — osteoporosis, osteoarthritis, stress fractures, joint replacement — demonstrates the consequences when the musculoskeletal system's adaptive capacity is overwhelmed. Osteoporosis is the failure of bone to maintain mass; osteoarthritis is the failure of cartilage to maintain matrix; stress fractures are the failure of remodeling to keep pace with microdamage. Each has a mechanical aetiology rooted in the principles covered here.

Historical notes Master

Julius Wolff (1836-1902) published Das Gesetz der Transformation der Knochen in 1892, articulating the principle that bone architecture reflects its mechanical loading history. Wolff's observations were qualitative, based on dissecting room specimens and the patterns of trabecular alignment in the femoral head. The cellular mechanisms — osteocytes as mechanosensors, the RANKL/OPG axis, sclerostin — were discovered more than a century later. Wilhelm Roux (1850-1924) independently proposed a similar idea (the "functional adaptation" of tissues) at roughly the same time, and the combined principle is sometimes called the Wolff-Roux law.

The lever analysis of the musculoskeletal system has its origins in Giovanni Alfonso Borelli (1608-1679), whose De Motu Animalium (1680) applied Archimedes' lever principles to animal movement. Borelli calculated the forces required for standing, walking, and flying using static equilibrium, recognising that muscles operate at mechanical disadvantage and must therefore produce forces far exceeding the external loads. His work founded the discipline of biomechanics.

The mathematical theory of elasticity, developed by Cauchy, Navier, and Lam'e in the early 19th century, provided the continuum-mechanics foundation for bone stress analysis. The finite element method was developed by engineers in the 1950s-60s (Turner, Clough, Martin, and Topp; Zienkiewicz) for aerospace structures and was first applied to bone by Rybicki and Simonen in the 1970s. Subject-specific FEA from CT data became feasible in the 1990s with advances in computing power and medical imaging.

DEXA was developed in the 1980s (Mazess, Wahner) as a refinement of single-photon and dual-photon absorptiometry, providing precise, low-dose measurements of areal bone mineral density. The WHO T-score definition of osteoporosis (1994) standardised diagnosis worldwide. The recognition that BMD alone is an imperfect predictor of fracture risk has driven the development of trabecular bone score (TBS), hip structure analysis (HSA), and FEA-based strength estimates as supplements to DXA.

The modern understanding of osteoarthritis as a whole-joint disease, rather than simple cartilage wear, emerged from the work of Dieppe, Felson, and others in the 1990s-2000s, integrating biomechanical, inflammatory, and genetic factors. The discovery of aggrecanases (ADAMTS-4, ADAMTS-5) by Sandy and Verscharen in the late 1990s identified a specific enzymatic pathway for cartilage degradation, opening therapeutic targets that are now in clinical trials.

Bibliography Master

  1. Wolff, J. — Das Gesetz der Transformation der Knochen. August Hirschwald, Berlin (1892).

  2. Borelli, G. A. — De Motu Animalium. Ex Typographia Angeli Bernabò, Rome (1680).

  3. Nigg, B. M. & Herzog, W. — Biomechanics of the Musculo-skeletal System, 3rd ed. Wiley (2007).

  4. Sherwood, L. — Human Physiology, 9th ed. Cengage (2016).

  5. Silverthorn, D. U. — Human Physiology: An Integrated Approach, 8th ed. Pearson (2019).

  6. Currey, J. D. — Bones: Structure and Mechanics, 2nd ed. Princeton University Press (2006).

  7. Felson, D. T. — Osteoarthritis as a disease of mechanics. Osteoarthritis Cartilage 21, 10-15 (2013).

  8. Kanis, J. A. on behalf of WHO — Assessment of fracture risk and its application to screening for postmenopausal osteoporosis. WHO Technical Report Series 843 (1994).

  9. Rybicki, E. F. & Simonen, F. A. — Mechanics of oblong bone. In Advances in Bioengineering, ASME (1973).

  10. Martin, R. B., Burr, D. B. & Sharkey, N. A. — Skeletal Tissue Mechanics. Springer (1998).