Thermoelastic simulations Answered


Does OnScale support thermoelastic simulations? In other words, I would like to simulate a laser pulse thermal excitation and then measure displacements or stresses generated on a metallic structure generated by thermal gradients. How to apply a temperature profile as excitation force?



  • Official comment

    Hi Piotr,


    Yes this is possible in OnScale. Here is an example that should help you get an idea of how to set this up.


    Best Regards,


    Comment actions Permalink
  • Hi Chloe,

    Thank you for sharing the file. It is very helpful. I have modified my code following your example and it seems to be OK however, it can not generate stresses from the thermal excitation. I have activated the thermal solver, and applied boundary conditions in a way as below:

      side xmin absr
      side xmax absr

      side ymin symm
      side ymax absr

      side zmin free

      defn zmxt tmpr func
      symb bxmin = 0.0049 * $coordFactor - $box
      symb bymin = -1e-10 * $coordFactor - $box
      symb bzmin = 0.001 * $coordFactor - $box
      symb #get { is js ks } clsnode $bxmin $bymin $bzmin
      symb bxmax = 0.0051 * $coordFactor + $box
      symb bymax = 0.0001 * $coordFactor + $box
      symb bzmax = 0.001 * $coordFactor + $box
      symb #get { ie je ke } clsnode $bxmax $bymax $bzmax
      node $is $ie $js $je $ks $ke alum void

    Could I send to you my code for a veryfication?



    Comment actions Permalink
  • Hi Piotr,


    Yes of course, send your files to us and I will take a look for you.


    Best Regards,


    Comment actions Permalink
  • Hi Chloe and Piotr,

    I try to do same thing - modelling thermo-elastic generation of US. I studied the linked example from Chloe and was basically able to adapt it to my needs. Except for one thing: If in this example I change the thermal source definition (line 165) from

    defn ymxt tmpr func


    defn ymxt flux func

    it seems to be no longer working - temperatures and displacements are all zero. However, using "flux" and setting "cupl off" produces a fine temperature field. 

    Is it supposed not to work using "flux"? If so, does this have physical reasons?


    Comment actions Permalink
  • Comment actions Permalink

Please sign in to leave a comment.

Didn't find what you were looking for?

New post