Calibration of Maxwell Material
Example posted by: Dr. Dimitrios G. Lignos (McGill University)
This example demonstrates how to conduct a calibration of a viscous damper using the maxwell model.
The files needed to analyze this structure in OpenSees are included here:
- The main file: File:Maxwell Calibrator.tcl
Supporting procedure files
- File:SquareSsection.tcl – displays a square fiber section
- File:ViscousD.tcl – contains the displacement loading history of the damper in units of mm
All files are available in a compressed format here: Calibration_Maxwell_example.zip
The rest of this example describes the model and shows the analysis results.
The units of the model are mm, N, and seconds.
The basic geometry of the viscous damper is defined by input variables for the length L=5000mm, and area A = 12000mm2. Two nodes are used for the geometry of the damper.
A square section is used to define the area A of the damper (Adamper = 12000.0mm2).
Force-Based Beam-Column elements are used to link the two nodes that define the geometry of the viscous damper with n=5 sections of integration and the damper section defined previously.
Node 1 is pinned. A roller is assigned to node 2
To model the viscous damper the Maxwell is used. The input parameters that are selected for the damper example are as follows: Axial Stiffness K = 500000.0 N/mm, Viscous Coefficient C=100000.0 N(s/mm)0.30 and exponent a=0.30.
The viscous damper is loaded with a displacement history that is called Damper.txt in the axial loading direction. A linear pattern is selected and the node of application is Node 2.
The recorders used in this example include:
- The localForce recorder to track the axial force of the viscous damper
- The stressStrain recorder to track the stress and strain history of the viscous damper in section 3.
In order to calibrate the viscous damper a multiple support excitation option is selected with application of displacement as the imposed motion. This is necessary since an integration step dt should be specified for integration of strain histories for each time step.
Comparison of OpenSees Model with MATLAB based Maxwell model
The force - displacement relationship from the maxwell damper are shown in Figure 2. In order to obtain the displacement of the damper we multiply the strain in section 3 of the damper times the length L=5000mm. A comparison of the response with a MATLAB based script is also shown in the same figure. Results are identical.