| Matrix or Vector | Option | Shape Functions[1] | Integration Points |
|---|---|---|---|
| Stiffness Matrix | Longitudinal | Equation 11–15 | None |
| Torsional | Equation 11–18 | None | |
| Stress Stiffening Matrix | Longitudinal | Equation 11–7 and Equation 11–8 | None |
The user explicitly defines the force-deflection curve for COMBIN39 by the input of discrete points of force versus deflection. Up to 250 points on the curve can be defined, and are entered as real constants. The input curve must pass through the origin and must lie within the unshaded regions, if KEYOPT(1) = 1.
The input deflections must be given in ascending order, with the minimum change of deflection of:
 | (13–69) |
where:
| ui = input deflections (input as D1, D2, ... D250 on R or RMORE commands) |
|
| umax = most positive input deflection |
| umin = most negative input deflection |
During the stiffness pass of a given iteration, COMBIN39 will use the results of the previous iteration to determine which segment of the input force-deflection curve is active. The stiffness matrix and load vector of the element are then:
 | (13–70) |
 | (13–71) |
where:
| Ktg = slope of active segment from previous iteration (output as SLOPE) |
| F1 = force in element from previous iteration (output as FORCE) |
If KEYOPT(4) > 0, Equation 13–70 and Equation 13–71 are expanded to 2 or 3 dimensions.
During the stress pass, the deflections of the current equilibrium iteration will be examined to see whether a different segment of the force-deflection curve should be used in the next equilibrium iteration.
If KEYOPT(2) = 0 and if no force-deflection points are input for deflection less than zero, the points in the first quadrant are reflected through the origin (Figure 13.4: Input Force-Deflection Curve Reflected Through Origin).
If KEYOPT(2) = 1, there will be no stiffness for the deflection less than zero (Figure 13.5: Force-Deflection Curve with KEYOPT(2) = 1).
If KEYOPT(1) = 0, COMBIN39 is conservative. This means that regardless of the number of loading reversals, the element will remain on the originally defined force-deflection curve, and no energy loss will occur in the element. This also means that the solution is not path-dependent. If, however, KEYOPT(1) = 1, the element is nonconservative. With this option, energy losses can occur in the element, so that the solution is path-dependent. The resulting behavior is illustrated in Figure 13.6: Nonconservative Unloading (KEYOPT(1) = 1).
When a load reversal occurs, the element will follow a new force-deflection line passing through the point of reversal and with slope equal to the slope of the original curve on that side of the origin (0+ or 0-). If the reversal does not continue past the force = 0 line, reloading will follow the straight line back to the original curve (Figure 13.7: No Origin Shift on Reversed Loading (KEYOPT(1) = 1)).
If the reversal continues past the force = 0 line, a type of origin shift occurs, and reloading will follow a curve that has been shifted a distance uorig (output as UORIG) (Figure 13.8: Origin Shift on Reversed Loading (KEYOPT(1) = 1)).
A special option (KEYOPT(2) = 2) is included to model crushing behavior. With this option, the element will follow the defined tensile curve if it has never been loaded in compression. Otherwise, it will follow a reflection through the origin of the defined compressive curve (Figure 13.9: Crush Option (KEYOPT(2) = 2)).
