Personal tools
 

History MSE: multi-case

In this example time history curves are used to calibrate directly experiment vs. simulation curves with two analysis cases.

 

Solution with LS-OPTui

Solution with LS-OPTui

Solvers

Solvers

Start LS-OPTui.

→ Open the command file com.mse_his.

Solvers Panel

  1. Select the Solvers panel.
  2. For Command select your LS-DYNA executable.
  3. For Input File browse the file foam2.k in your directory.
  4. Enter the Name of Analysis Case: Case2.
  5. Push the Add button to create a new analysis case.
solvers_case22.png

Sampling

Sampling

Sampling

  1. Select the Sampling panel.
  2. Select the analysis case Case2 from the list.
  3. For POINT SELECTION select Duplicate.
  4. For Case choose Case1.
sampling2.png

Histories

Histories

Extract History

  1. Select the Histories panel.
  2. From the possible histories select NODOUT.
  3. For analysis Case select Case2.
  4. Type in the node ID: 288.
  5. For Direction select Z Component.
  6. Type in Disp2 for the History Name.
  7. Push the Add button.
history_disp2.png
  1. From the possible histories select RCFORC.
  2. Type in for Interface ID: 1.
  3. For Component select Z slave force.
  4. Type in Force2 for the History Name.
  5. Push the Add button.
history_force2.png
  1. From the possible histories select Crossplot.
  2. For z(t) type in -Disp2 (negative Disp1).
  3. For F(t) type in Force2.
  4. Type in F2_vs_d2 as History Name.
  5. Push the Add button to create a new history.
history_f2-vs-d2.png
  1. From the possible histories select File.
  2. Push the Browse button to find the file Test2.txt from your main directory.
  3. Type in FILE2 as History Name.
  4. Push the Add button to create a new history.
 
history_test2.png

Responses

Responses

Define the Responses

  1. Select the Responses panel.
  2. From the possible response types select: MeanSqErr.
  3. For Target curve select FILE2.
  4. For Computed curve select F2_vs_d2.
  5. Type in the label MSE2 for Response Name.
  6. Push the Add button.
response_mse2.png

Objective

Objective

Objectives

  1. Select the Objective panel.
  2. Select MSE from Response as the objective function.
  3. Leave the default 1.0 for Weight.
  4. Select MSE2 from Response as the objective function.
  5. Leave the default 1.0 for Weight.
objective_mse2.png

Run

Run

Run the Optimization

  1. Select the Run panel.
  2. For Number of iterations enter 3.
  3. Push the Run button to start the optimization.
run2.png

Com-file

Com-file

The created command file may look like this:

"history based MSE"
$ Created on Fri Jan 11 14:58:45 2008
solvers 2
$
$ WARNING -- NO RESPONSES ARE DEFINED
$
histories 8
$
$ DESIGN VARIABLES
$
variables 2
 Variable 'YMod' 700000
  Lower bound variable 'YMod' 500000
  Upper bound variable 'YMod' 2e+06
 Variable 'Yield' 1500
  Lower bound variable 'Yield' 500
  Upper bound variable 'Yield' 2000

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$      OPTIMIZATION METHOD   
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
Optimization Method SRSM

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$      SOLVER "Case1"
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$ DEFINITION OF SOLVER "Case1"
$
 solver dyna960 'Case1'
  solver command "/home/prak1/LS-DYNA/ls971_s_7600_1224_ia32_redhat90"
  solver input file "/home/prak1/LS-OPT/beispiele/parameter_identification/history_based_multi/foam1.k"
  solver check output on 
  solver compress d3plot off 
$ ------ Pre-processor --------
$   NO PREPROCESSOR SPECIFIED
$ ------ Metamodeling ---------
  solver order linear
  solver experiment design dopt
$ ------ Job information ------
  solver concurrent jobs 1
$
$ WARNING - NO RESPONSES DEFINED FOR SOLVER "Case1"
$
$
$ HISTORIES FOR SOLVER "Case1"
$
 history 'Disp1' "BinoutHistory -res_type nodout  -cmp z_displacement -id 296"
 history 'Force1' "BinoutHistory -res_type RCForc -cmp z_force  -id 1  -side SLAVE"
$
$ HISTORY EXPRESSIONS FOR SOLVER "Case1"
$
 history 'F_vs_d' expression {Crossplot("-Disp1","Force1")}

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$      SOLVER "Case2"
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$ DEFINITION OF SOLVER "Case2"
$
 solver dyna960 'Case2'
  solver command "/home/prak1/LS-DYNA/ls971_s_7600_1224_ia32_redhat90"
  solver input file "/home/prak1/LS-OPT/beispiele/parameter_identification/history_based_multi/foam2.k"
  solver check output on 
  solver compress d3plot off 
$ ------ Pre-processor --------
$   NO PREPROCESSOR SPECIFIED
$ ------ Metamodeling ---------
  solver order linear
   solver experiment duplicate 'Case1'
$ ------ Job information ------
  solver concurrent jobs 1
$
$ WARNING - NO RESPONSES DEFINED FOR SOLVER "Case2"
$
$
$ HISTORIES FOR SOLVER "Case2"
$
 history 'Disp2' "BinoutHistory -res_type nodout  -cmp z_displacement -id 288"
 history 'Force2' "BinoutHistory -res_type RCForc -cmp z_force  -id 1  -side SLAVE"
$
$ HISTORY EXPRESSIONS FOR SOLVER "Case2"
$
 history 'F2_vs_d2' expression {Crossplot("-Disp2","Force2")}

$
$ HISTORIES FROM FILES
$
 history 'Test1' file "/home/prak1/LS-OPT/beispiele/parameter_identification/history_based_multi/Test1.txt"
 history 'Test2' file "/home/prak1/LS-OPT/beispiele/parameter_identification/history_based_multi/Test2.txt"
composites 2
$
$ COMPOSITE EXPRESSIONS
$
 composite 'MSE' {MeanSqErr(Test1,F_vs_d)}
 composite 'MSE2' {MeanSqErr(Test2,F2_vs_d2)}
$
$ OBJECTIVE FUNCTIONS
$
 objectives 2
 objective 'MSE' 1
 objective 'MSE2' 1
$
$ THERE ARE NO CONSTRAINTS!!!
$
 constraints 0
$
$ JOB INFO
$
 concurrent jobs 1
 iterate param design 0.01
 iterate param objective 0.01
 iterate param stoppingtype or
 iterate 3
STOP

Results

Results

Confidence Intervals

Confidence Intervals

90% Confidence Intervals

The confidence intervals for the two parameters YMod and Yield can be found in the lsopt_report file in the main directory.

→ YMod is rather insignificant and Yield is significant.

conf_interval2.png

Opt History

Opt History

Optimization History

  1. For Type of Plot select Opt History.
  2. From Entity to Monitor select Composites → MSE2.
  3. Click near the iteration number for information and curve comparison.
opthistory_mse22.png
  1. Shows the iteration number.
  2. View the entity values for the current iteration.
  3. Push the MeanSqErr curves button to view a comparison of the history in LS-PrePost.
opthistory_mse33.png
History Comparison between the cross plot F2_vs_d2 (green line) and the imported file (red line) in LS-PrePost.compare_history.png

Download

Download

The complete data set (input and command files) is available for download as tar.gz history_MSE_multi-case.tar.gz or zip-file history_MSE_multi-case.zip.