[CIG-SHORT] Initial stress cell_data_fields settings and interpolation

[Brad Aagaard]

On 11/06/2015 11:06 AM, Chunyu Liu (cliu5) wrote:
> Hi all,
>
> Is there any way I can get initial stress data( not that resulting stress)  in .VTK file?
> I do not find any settings for the initial stress in “ cell_info_fields” and “cell-data_fields”.
>
> For example:
> cell_data_fields =[ stress, “stress_initial”]

This is an oversight on our part. It is not currently available. A 
workaround is to run a static simulation (single time step) without any 
deformation to get the "initial" stress.

> Q2: Interpolation problem
>
> 1.   I use linear interpolation.
> 2.   several initial stress lines. For example:
>
> // Columns are
>   25 // (1) x coordinates(m)
>   26 // (2) y coordinates(m)
>   27 // (3) z coordinates(m)
>   28 // (4) stress-xx (Pa)
>   29 // (5) stress-yy (Pa)
>   30 // (6) stress-zz (Pa)
>   31 // (7) stress-xy (Pa)
>   32 // (8) stress-yz (Pa)
>   33 // (9) stress-xz (Pa)
>   34 -222000.000000  -222000.000000  0.000000        -37264580.000000        -131340346.000000       0.000000        7981201.000000  0.000000        0.000000
>   35 -111000.000000  -222000.000000  0.000000        -36546640.000000        -129611870.000000       0.000000        8891624.000000  0.000000        0.000000
>   36 0.000000        -222000.000000  0.000000        -36008091.000000        -127917852.000000       0.000000        9781642.000000  0.000000        0.000000
>   37 111000.000000   -222000.000000  0.000000        -35531603.000000        -126206836.000000       0.000000        10653143.000000 0.000000        0.000000
>   38 222000.000000   -222000.000000  0.000000        -34957300.000000        -124402730.000000       0.000000        11503441.000000 0.000000        0.000000
>   39 -222000.000000  -111000.000000  0.000000        -22838144.000000        -116776288.000000       0.000000        11259071.000000 0.000000        0.000000
>   40 -111000.000000  -111000.000000  0.000000        -22551422.000000        -115173226.000000       0.000000        11932060.000000 0.000000        0.000000
>   41 0.000000        -111000.000000  0.000000        -22180910.000000        -113548630.000000       0.000000        12538850.000000 0.000000        0.000000
>   42 111000.000000   -111000.000000  0.000000        -21762040.000000        -111904430.000000       0.000000        13105260.000000 0.000000        0.000000
>   43 222000.000000   -111000.000000  0.000000        -21621862.000000        -110393807.000000       0.000000        13642334.000000 0.000000        0.000000
>   44 -222000.000000  0.000000        0.000000        -43329544.000000        -117636548.000000       0.000000        7273745.000000  0.000000        0.000000
>   45 -111000.000000  0.000000        0.000000        -42661628.000000        -115014364.000000       0.000000        8733876.000000  0.000000        0.000000
>   46 0.000000        0.000000        0.000000        -40696706.000000        -111989218.000000       0.000000        10180400.000000 0.000000        0.000000
>   47 111000.000000   0.000000        0.000000        -39202640.000000        -109388170.000000       0.000000        11592650.000000 0.000000        0.000000
>   48 222000.000000   0.000000        0.000000        -38123380.000000        -107149580.000000       0.000000        12943260.000000 0.000000        0.000000

Your points are not on a line. It looks like they are on a plane, so you 
should use area interpolation (dim = 2).

Regards,
Brad