SmartAnalyze

class opstool.analysis.SmartAnalyze(analysis_type='Transient', **kargs)

The SmartAnalyze is a class to provide OpenSeesPy users a easier way to conduct analyses. Original Tcl version Author: Dr. Dong Hanlin (http://www.hanlindong.com/2019/opensees-converge/) Here’s the converted python version, with some modifications.

Parameters

analysis_type: str, default=”Transient”

Assign the analysis type, “Transient” or “Static”.

Other Parameters that control convergence

TEST RELATED:

testType: str, default=”EnergyIncr”

Identical to the testType in OpenSees test command. Choices see http://opensees.berkeley.edu/wiki/index.php/Test_Command.

testTol: float, default=1.0e-12

The initial test tolerance set to the OpenSees test command. If tryLooseTestTol is set to True, the test tolerance can be loosen.

testIterTimes: int, default=10

The initial number of test iteration times. If tryAddTestTimes is set to True, the number of test times can be enlarged.

testPrintFlag: int, default=0

The test print flag in OpenSees Test command. Choices see http://opensees.berkeley.edu/wiki/index.php/Test_Command.

tryAddTestTimes: bool, default=False

If True,the number of test times will be enlarged if the last test norm is smaller than normTol, the enlarged number is specified in testIterTimesMore. Otherwise, the number of test times will always be equal to testIterTimes.

normTol: float, default=1.e3

Only useful when tryAddTestTimes is True. If unconverge, the last norm of test will be compared to normTol. If the norm is smaller, the number of test times will be enlarged.

testIterTimesMore: int, default=50

Only useful when tryaddTestTimes is True. If unconverge and norm is ok, the test iteration times will be set to this number.

tryLooseTestTol: bool, default=False

If this is set to True, if unconverge at minimum step, the test tolerance will be loosen to the number specified by looseTestTolTo.the step will be set back.

looseTestTolTo: float, default=1.0

Only useful if tryLooseTestTol is True. If unconvergance at the min step, the test tolerance will be set to this value.

ALGORITHM RELATED:

tryAlterAlgoTypes: bool, default=False

If True, different algorithm types specified in algoTypes will be tried during unconvergance. If False, the first algorithm type specified in algoTypes will be used.

algoTypes: list[int], default=[10, 20, 30, 40, 50]

A list of flags of the algorithms to be used during unconvergance. The integer flag is documented in the following section. Only useful when tryAlterAlgoTypes is True. The first flag will be used by default. The algorithm command in the model will be ignored. If you need other algorithm, try a user-defined algorithm. See the following section.

STEP RELATED:

initialStep: float, default=None

Specifying the initial Step length to conduct analysis. If None, equal to dt.

relaxation: float, between 0 and 1, default=0.5

A factor that is multiplied by each time the step length is shortened.

minStep: float, default=1.e-6

The step tolerance when shortening the step length. If step length is smaller than minStep, special ways to converge the model will be used according to try- flags.

LOGGING RELATED:

printPer: int, default=10

Print to the console every several trials.

debugMode: bool, default=False

Print as much information as possible.

Tip

• 0: Linear

• 1: Linear -initial

• 2: Linear -factorOnce

• 10: Newton

• 11: Newton -initial

• 12: Newton -initialThenCurrent

• 20: NewtonLineSearch

• 21: NewtonLineSearch -type Bisection

• 22: NewtonLineSearch -type Secant

• 23: NewtonLineSearch -type RegulaFalsi

• 30: ModifiedNewton

• 31: ModifiedNewton -initial

• 40: KrylovNewton

• 41: KrylovNewton -iterate initial

• 42: KrylovNewton -increment initial

• 43: KrylovNewton -iterate initial -increment initial

• 44: KrylovNewton -maxDim 6

• 50: SecantNewton

• 51: SecantNewton -iterate initial

• 52: SecantNewton -increment initial

• 53: SecantNewton -iterate initial -increment initial

• 60: BFGS

• 70: Broyden

• 80: PeriodicNewton

• 90: User-defined0

Examples

Note

test() and algorithm() will run automatically in SmartAnalyze, but commands such as integrator() must be defined outside SmartAnalyze.

Example 1: Basic usage for Transient

>>> import opstool as opst
>>> ops.constraints('Transformation')
>>> ops.numberer('Plain')
>>> ops.system('BandGeneral')
>>> ops.integrator('Newmark', 0.5, 0.25)
>>> analysis = opst.SmartAnalyze(analysis_type="Transient")
>>> segs = analysis.transient_split(npts)
>>> for seg in segs:
>>>     analysis.TransientAnalyze(dt)

Example 2: Basic usage for Static

>>> import opstool as opst
>>> ops.constraints('Transformation')
>>> ops.numberer('Plain')
>>> ops.system('BandGeneral')
>>> protocol=[1, -1, 1, -1, 0]
>>> analysis = opst.SmartAnalyze(analysis_type="Static")
>>> segs = analysis.static_split(protocol, 0.01)
>>> for seg in segs:
>>>     analysis.StaticAnalyze(1, 1, seg)

Example 3: change control parameters

>>> analysis = opst.SmartAnalyze(analysis_type="Transient",
>>>                              algoTypes=[20, 30],
>>>                              printPer=20,
>>>                              tryAlterAlgoTypes=True,
>>>                             )

StaticAnalyze(node, dof, seg)

Single Step Static Analysis.

Parameters

nodeint

The node tag in the displacement control.

dofint

The dof in the displacement control.

segfloat

Each load step, i.e. each element returned by static_split.

Returns

None

TransientAnalyze(dt)

Single Step Transient Analysis.

Parameters

dtfloat

Time Step.

Returns

None

static_split(targets, maxStep=None)

Returns a sequence of substeps for static analysis, for use in outer analysis loops.

Parameters

targets: list

A list of target displacements, the first element must be positive.

maxStep: float, default=None

The maximum step length in the displacement control. If None, targets[1] - targets[0].

Returns

segs: list

A sequence of substeps for static analysis.

transient_split(npts)

Step Segmentation for Transient Analysis.

Parameters

nptsint

Total steps for transient analysis.

Returns

A list to loop.