熱座屈解析
Theory and Physics
What is Thermal Buckling?
Professor, can buckling occur due to heat? Even when no force is applied?
Good question. Thermal buckling is buckling caused by compressive stress generated by temperature changes. If a structure can expand freely, thermal stress does not occur, but if expansion is constrained, compressive stress is generated.
Is the bending of railroad rails on hot days also thermal buckling?
Exactly. Rails are constrained in axial expansion by the sleepers. The axial stress generated by a temperature rise $\Delta T$ is:
Here $\alpha$ is the coefficient of linear expansion. For steel ($\alpha \approx 12 \times 10^{-6}$ /°C, $E = 200$ GPa) with $\Delta T = 40$ °C, $\sigma_{th} = 96$ MPa. If this exceeds the buckling stress, the rail bends sideways.
96 MPa from a 40°C temperature difference... that's larger than I thought.
Yes. Thermal stress is proportional to the degree of constraint. Under full constraint, it is $E\alpha\Delta T$, but if expansion is partially allowed, the stress decreases. Correctly evaluating the degree of constraint in real structures is the first step in thermal buckling analysis.
Governing Equations for Thermal Buckling
Is the formulation of thermal buckling different from mechanical buckling?
The framework of eigenvalue buckling is the same. However, the reference load becomes the thermal load:
1. Perform a thermal stress analysis given a temperature distribution $\Delta T(x,y,z)$
2. Construct $[K_\sigma]$ from the obtained thermal stress
3. Solve the eigenvalue problem $([K_0] + \lambda [K_\sigma])\{\phi\} = \{0\}$
$\lambda$ is the buckling thermal load factor. $\lambda \cdot \Delta T$ is the critical temperature rise.
There are cases where the temperature distribution is not uniform, right?
In real structures, temperature distribution is usually non-uniform. For example, a beam during a fire has a hot bottom surface and a relatively cool top surface. If only one side of a plate is heated, a bending moment due to the temperature gradient also occurs.
There are two effects of temperature:
- Membrane stress (in-plane compression/tension) — The driving force for buckling
- Bending moment (temperature gradient through the thickness) — Induces additional deformation
Thermal Buckling of Plates
Is there a theoretical solution for thermal buckling of plates?
Critical temperature rise for a uniformly heated rectangular plate with four edges constrained:
The square of the ratio of plate thickness $t$ to width $b$ is included. Thinner plates buckle at lower temperature rises.
Moreover, Young's modulus $E$ is not in this formula! Under ideal constraint conditions, the critical temperature depends only on geometric parameters ($t/b$) and the coefficient of linear expansion $\alpha$, not on the material's stiffness.
Being independent of material is counterintuitive...
Under full constraint, thermal stress is $E\alpha\Delta T$, and buckling stress is $E \cdot f(t/b)$, so when taking the ratio, $E$ cancels out. However, real structures are not fully constrained, so the influence of $E$ appears.
Thermal Buckling Problems in Real Structures
Please give me examples of real-world problems where thermal buckling is an issue.
Typical examples:
| Structure | Situation | Characteristics |
|---|---|---|
| Railroad Rails | Summer temperature rise | Lateral buckling of continuously welded rails |
| Pipelines | Transport of high-temperature fluid | Lateral buckling (snaking) of subsea pipelines |
| Aircraft Skin Panels | Aerodynamic heating during supersonic flight | Buckling of thin panels |
| Space Structures | Temperature difference between sunlit/shaded areas | Solar panels, Antennas |
| Welded Structures | Local heating during welding | Welding distortion (angular distortion, buckling distortion) |
| Buildings during Fire | Heating of steel beams | Axial force buckling of beams with constrained ends |
Subsea pipeline buckling sounds interesting.
Thermal buckling of subsea pipelines is a critical issue in the oil and gas industry. When hot crude oil flows through, the pipe tries to expand but is constrained by seabed friction. When the temperature rises beyond a certain point, the pipe snakes laterally (lateral buckling). There is a design method called "designed lateral buckling" to intentionally control this.
Summary
Let me organize the theory of thermal buckling.
Key points:
- Thermal buckling is caused by compressive stress from "constrained thermal expansion" — No buckling if expansion is free
- $\sigma_{th} = E\alpha\Delta T$ — Thermal stress under full constraint
- Evaluation of the degree of constraint is key — Real structures are neither fully constrained nor completely free
- In FEM, a two-step process: thermal stress analysis → eigenvalue buckling analysis — Can handle non-uniform temperature distributions
- Temperature gradient generates both membrane stress and bending — One-sided heating is particularly dangerous
The scary part about thermal buckling is that "buckling can occur without applying force." Temperature management directly affects structural safety.
Exactly. Thermal buckling is easily overlooked in design because the "load" is invisible. Structures used in environments with large temperature changes should always consider the possibility of thermal buckling.
Thermal Buckling and the Bimetallic Strip
Buckling due to heat occurs when compressive thermal stress reaches the buckling load. The "bimetallic strip" (a two-layer metal plate with different thermal expansion coefficients) experimented by Johann Christoph Baugarten (Germany) in the 1780s is precisely the prototype of thermal buckling. Thermostats operating on this principle were industrialized in the 1900s and are still used for automatic cutoff in power and air conditioning systems.
Physical Meaning of Each Term
- Inertia term (mass term): $\rho \ddot{u}$, i.e., "mass × acceleration". Haven't you experienced your body being thrown forward during sudden braking? That "feeling of being carried away" is precisely the inertial force. Heavier objects are harder to set in motion and harder to stop once moving. Buildings shake during earthquakes because the ground moves suddenly while the building's mass "gets left behind". In static analysis, this term is set to zero, which assumes "forces are applied slowly so acceleration can be ignored". It absolutely cannot be omitted in impact loads or vibration problems.
- Stiffness term (elastic restoring force): $Ku$ or $\nabla \cdot \sigma$. When you pull a spring, you feel a "force trying to return", right? That is Hooke's law $F=kx$, the essence of the stiffness term. Now a question — an iron rod and a rubber band, which stretches more under the same pulling force? Obviously the rubber. This "resistance to stretching" is Young's modulus $E$, which determines stiffness. A common misconception: "high stiffness ≠ strong". Stiffness is "resistance to deformation", strength is "resistance to failure" — different concepts.
- External force term (load term): Body force $f_b$ (gravity, etc.) and surface force $f_s$ (pressure, contact force, etc.). Think of it this way — the weight of a truck on a bridge is a "force acting on the entire contents" (body force), the force of the tire pushing on the road surface is a "force acting only on the surface" (surface force). Wind pressure, water pressure, bolt tightening force... all are external forces. A typical mistake here: getting the load direction wrong. Intending "tension" but it becomes "compression" — sounds like a joke, but it actually happens when coordinate systems are rotated in 3D space.
- Damping term: Rayleigh damping $C\dot{u} = (\alpha M + \beta K)\dot{u}$. Try plucking a guitar string. Does the sound continue forever? No, it gradually fades. Because vibration energy is converted to heat by air resistance and internal friction in the string. Car shock absorbers work on the same principle — intentionally absorbing vibration energy to improve ride comfort. What if damping were zero? Buildings would continue shaking forever after an earthquake. Since that doesn't happen in reality, setting appropriate damping is important.
Assumptions and Applicability Limits
- Continuum assumption: Treats material as a continuous medium, ignoring microscopic heterogeneity
- Small deformation assumption (for linear analysis): Deformation is sufficiently small compared to initial dimensions, and stress-strain relationship is linear
- Isotropic material (unless otherwise specified): Material properties are independent of direction (anisotropic materials require separate tensor definitions)
- Quasi-static assumption (for static analysis): Ignores inertial/damping forces, considering only equilibrium between external and internal forces
- Non-applicable cases: Large deformation/large rotation problems require geometric nonlinearity. Nonlinear material behavior like plasticity or creep requires constitutive law extensions.
Dimensional Analysis and Unit Systems
| Variable | SI Unit | Notes / Conversion Memo |
|---|---|---|
| Displacement $u$ | m (meter) | When inputting in mm, unify loads/elastic modulus to MPa/N system |
| Stress $\sigma$ | Pa (Pascal) = N/m² | MPa = 10⁶ Pa. Be careful of unit system inconsistency when comparing with yield stress |
| Strain $\varepsilon$ | Dimensionless (m/m) | Note the distinction between engineering strain and logarithmic strain (for large deformation) |
| Elastic modulus $E$ | Pa | Steel: ~210 GPa, Aluminum: ~70 GPa. Note temperature dependence |
| Density $\rho$ | kg/m³ | In mm system: tonne/mm³ (= 10⁻⁹ tonne/mm³ for steel) |
| Force $F$ | N (Newton) | Unify to N in mm system, N in m system |
Numerical Methods and Implementation
Thermal Buckling Analysis by FEM
Please explain the procedure for analyzing thermal buckling with FEM.
It's a two-step (or three-step) procedure.
Step 1: Determination of Temperature Field
Determine the temperature distribution $T(x,y,z)$ via steady-state or transient heat conduction analysis. For simple uniform heating, temperature can be assigned manually.
Step 2: Thermal Stress Analysis
Input the obtained temperature field into structural analysis to calculate the stress distribution due to thermal expansion. Boundary conditions (constraint conditions) greatly influence the results at this stage.
Step 3: Eigenvalue Buckling Analysis
Construct the geometric stiffness matrix from the thermal stress in Step 2 and solve the eigenvalue problem. $\lambda$ is the critical thermal load factor.
What about cases combined with mechanical loads?
When mechanical and thermal loads act simultaneously, apply the mechanical load as a preload and evaluate buckling against the additional thermal load. If temperature and mechanical load vary independently, note that the "proportional load" assumption of eigenvalue buckling breaks down.
Nastran
```
SOL 105
CEND
SUBCASE 1
TEMPERATURE(LOAD) = 100
SPC = 1
METHOD = 10
```
Specify temperature distribution with the TEMPERATURE card. Define nodal temperatures with TEMP/TEMPD cards. Set the thermal expansion coefficient (A) in MAT1 correctly.
Abaqus
Related Topics
なった
詳しく
報告