Why does OptiStruct result differ greatly from those of ANSYS?

13

Answers

  • Altair Forum User
    Altair Forum User
    Altair Employee
    edited November 2020

    Hi @Prakash Pagadala,

    I try to convert the elements to 2nd order. But the results is similar to the first order, and it also differ from Ansys.

    How to modify integration points? I cannot find the place to define it in hypermesh.

    <?xml version="1.0" encoding="UTF-8"?>2017-09-20_164131.thumb.png.8f0d7a45c8b606ec287035d51958d99c.png<?xml version="1.0" encoding="UTF-8"?>2017-09-20_164253.png.341545325df792d22fb3ac7f06893474.png

     

    Roy

    Hi @Prakash Pagadala,

    Could you give me some help on my initial doubts? I'm really confused about that. From Ansys, we know that it use a second order element SOLID186 with full integration. How can we do this in OptiStruct? I want to know how this different result coming from.

    Thank you in advance.

     

    Roy

    Unable to find an attachment - read this blog

  • Altair Forum User
    Altair Forum User
    Altair Employee
    edited September 2017

     

     

     

    I saw the displacement results are similar, just stress is different. I think it is because stress is not output at outer position.

    Enter panel analysis>control cards>GLOBAL OUTPUT REQUEST, slide down to activate 'STRESS' and select 'LOCATION' as 'ALL'

    Run anaylysis again and open Hyperview>Contour> select stress result and activate 'Use corner data'

    The max stress values are 133 and 143(XX) , I think they are comparable with ansys

     

    <?xml version="1.0" encoding="UTF-8"?>Capture.PNG.4680d25a044fe87ac374ac5ca39bd6ab.PNG

    Hi @tinh,

    Could you show me some message about 'corner data'? From the help manual of HyperView, I only found following message. However, I don't think I have understood the meaning.

    2017-09-22_121526.png.d6e9537250c0f843e3c4cc693f1abfab.png

     

    2017-09-22_121553.png.9b6b24eb2368ff5a11ba1b1c72b55967.png

     

    Roy

  • tinh
    tinh Altair Community Member
    edited September 2017

    Hi,

    when you do not use corner data => stress is value at element centroid => stress is constant in each element

    Capture.PNG.2994e8a4c7d0992b7b92701526af53f6.PNG

     

     

    when you use corner data => stress is interpolate from center to element corners => stress is varying in each element:

    Capture2.PNG.2fff904303a90da1beeda1598ce0d474.PNG

     

     

    So when you use corner data => you got stress value at outer face of solid (interpolate), normally higher than at centroid of solid element

    if you cover solid by membrane, you got more exact stress result. because stress at membrane elements does not need to be interpolated from centroid value

  • Q.Nguyen-Dai
    Q.Nguyen-Dai Altair Community Member
    edited September 2017

    Hi,

    ...

    when you use corner data => stress is interpolate from center to element corners => stress is varying in each element:

     

    I think stress will be interpolated from integration points of each element, not from 'center'.

     

  • tinh
    tinh Altair Community Member
    edited September 2017

    Yes, 

    sorry I am wrong.

    and also the value at centroid is came from IP values, isn't it?

    maybe solver output at center as default to save disk space

  • Altair Forum User
    Altair Forum User
    Altair Employee
    edited September 2017

    Hi Roy,

     

    You can use shell elements with minimal thickness to simulate membrane effect. 

     

    @tinh

     

    It is normal that stress varies within element. So I think the centroid stress value is from the IP value

  • Altair Forum User
    Altair Forum User
    Altair Employee
    edited September 2017

    Yes, 

    sorry I am wrong.

    and also the value at centroid is came from IP values, isn't it?

    maybe solver output at center as default to save disk space

    What if more than one integration points exist in one element? I think in this situation the stress in one element may be not constant. 

    In Ansys software, it also has nodal stress besides elemental stress. I don't know the stress in OptiStruct is elemental or nodal?

     

    Roy

  • Altair Forum User
    Altair Forum User
    Altair Employee
    edited September 2017

    Hi Roy,

     

    You can use shell elements with minimal thickness to simulate membrane effect. 

     

    @tinh

     

    It is normal that stress varies within element. So I think the centroid stress value is from the IP value

    Hi Prakash,

    However the shell nodes in optistruct have six degree of freedom which is not compatible with solid nodes only with 3 translational d.o.f. 

     

    Roy

  • Q.Nguyen-Dai
    Q.Nguyen-Dai Altair Community Member
    edited September 2017

    Yes, 

    sorry I am wrong.

    and also the value at centroid is came from IP values, isn't it?

    maybe solver output at center as default to save disk space

     

    I think 'centroid' value is just 'mean' value over Integration Points?

  • Altair Forum User
    Altair Forum User
    Altair Employee
    edited September 2017

     

    I think 'centroid' value is just 'mean' value over Integration Points?

    I also want to know this. However, I have not found any information about this. For Ansys, it seems that the mean element stress is the average value of its nodes stress. 

  • Altair Forum User
    Altair Forum User
    Altair Employee
    edited September 2017

    Hi Prakash,

    How to define second order full integration solid elements in OptiStruct? Ansys use the second order SOLID186 element and I want to check it in OptiStruct.

    In addition, could you explain what's the locking effect in solid?

     

    Best Wishes

    Roy

  • Altair Forum User
    Altair Forum User
    Altair Employee
    edited September 2017

    Hi Roy,

     

    You can use ISOP to define full integration

    <?xml version="1.0" encoding="UTF-8"?>****.PNG

  • Q.Nguyen-Dai
    Q.Nguyen-Dai Altair Community Member
    edited November 2020

    @Roy Duan: Here' are two tests with HEXA20+Skin (membrane) in SAMCEF & Optistruct:

     

     

    Unable to find an attachment - read this blog

  • Altair Forum User
    Altair Forum User
    Altair Employee
    edited November 2020

    Hi Roy,

     

    You can use ISOP to define full integration

    <?xml version="1.0" encoding="UTF-8"?>****.PNG

    Hi @Prakash Pagadala,

    The stress results are exactly the same no matter how the ISOP parameter in PSOLID card values. Therefore, I don’t think the value of ISOP can change the element integration points. What’s more, from the manual of PSOLID card, the ISOP parameter is a special integration schemes for elasto-plastic nonlinear quasi-static analysis. But my example is a linear static analysis. So how to change the element integration point formulation in OptiStruct for linear static analysis? I think this is a fundamental set up for finite element analysis.

    <?xml version="1.0" encoding="UTF-8"?>2017-09-26_203519.thumb.png.f0919e016bf7b25e6f419247f769642e.png

     

    The model and the pdf file shows detail result. 

     

    Roy

     

    Unable to find an attachment - read this blog

  • Altair Forum User
    Altair Forum User
    Altair Employee
    edited November 2020

    @Roy Duan: Here' are two tests with HEXA20+Skin (membrane) in SAMCEF & Optistruct:

     

     

    Hi @Q.Nguyen-Dai,

    It's strange I got a different result from yours. I have uploaded my model. The model came from your attachments, but I have change some parameter in skin pshell card because some errors. Could you show me how to set the right value to skin pshell card? Thank you.

    2017-09-26_204559.png.33cb41c9880cb072a63b51b981e00ff6.png2017-09-26_204433.png.a1a6000d533b05a858ad3e5cb3b5f4e6.png

    Unable to find an attachment - read this blog

  • Altair Forum User
    Altair Forum User
    Altair Employee
    edited November 2020

    Hi all,

    From some test I found that the OptiStruct stress result is similar to Workbench element mean stress result (please see the attachments). So I suspect the stress output by OptiStruct may be element mean stress. One element mean stress means the average stress of all its nodes. For example, the element mean stress of a first order hexa element is the average stress of its 8 nodes which attached to the element.

    <?xml version="1.0" encoding="UTF-8"?>2017-09-26_210703.thumb.png.a17f98e83105a28936c2cb7dfceb06cc.png

     

    1. So does anyone can tell me which stress result had been output by OptiStruct? 

    2. How to output element integration point stress in OptiStruct? I think maybe the integration point stress should be compared since it had highest precision. 

    3. How many default  integration points for different elements in OptiStruct? And how to change the integration points?

     

    Best Wishes

    Roy

    Unable to find an attachment - read this blog

  • Q.Nguyen-Dai
    Q.Nguyen-Dai Altair Community Member
    edited September 2017

    @Roy Duan: Yes, I forgot to tell you modify pshell card for membrane behavior only.

    Here's my used pshell card:

    pshell_membrane.png.157302857d46584139bdb6469c279738.png

     

  • Altair Forum User
    Altair Forum User
    Altair Employee
    edited September 2017

    @Roy Duan: Yes, I forgot to tell you modify pshell card for membrane behavior only.

    Here's my used pshell card:

    <?xml version="1.0" encoding="UTF-8"?>pshell_membrane.png.157302857d46584139bdb6469c279738.png

     

    Hi @Q.Nguyen-Dai,

    I have gotten the same result as you by setting the pshell card as you recommend. But how to verify this is a membrane element not a shell element? Thank you.

     

    Roy

  • tinh
    tinh Altair Community Member
    edited December 2017

    Hi all,

    I am again with this topic because we talked about 'corner data' of stress

    when I perform NLGEOM analysis, select GLOBAL_OUTPUT_REQUEST > stress > CORNER

    but I don't see available corner data by Hyperview.

    Is there any limited of this output request?

    this is my *.fem file 

    Unable to find an attachment - read this blog

  • Altair Forum User
    Altair Forum User
    Altair Employee
    edited December 2017

    Hi Tinh,

     

    Conner stress data is not available for NLGEOM analysis type. 

  • tinh
    tinh Altair Community Member
    edited December 2017

    Thanks Prakash

    Is there any specific reason ? I am supprise with this :huh:/emoticons/default_huh.png' title=':huh:' />

  • Altair Forum User
    Altair Forum User
    Altair Employee
    edited December 2017

    TInh,

     

    I need to check about this. But what I know is because RADIOSS does not support.

     

    Maybe you can try with NLSTAT with large displacement

  • tinh
    tinh Altair Community Member
    edited December 2017

     

    This's special technique we use often with SAMCEF's solver to catch 'surface stress' of solid.

    From solid mesh we extract the external faces and define them as a very small thickness (1.0E-5 mm) membrane.

     

    Hi Mr Q.Nguyen

    I have a doubt about above method although some documents also record that it may give more accurate result (altair's document said that, too).

    Is it because SAMCEF does not output stress at grid points?

    Optistruct can output stress at grid points (use card STRESS>LOCATION, or GPSTRESS)

    and I saw those values is nearer to calculation than shell membrane stress, when i do a simple example

     

    a round section bar, fixed one end and loaded on the other one.

    - in beam model, stress (XX) is 10.19 (~calculation by c*M/I)

    - in solid model, stress is 10.22

    - in membrane shell, stress is 9.7

     

    maybe, the method should not be applied in linear static analysis?

    Unable to find an attachment - read this blog

  • Q.Nguyen-Dai
    Q.Nguyen-Dai Altair Community Member
    edited December 2017

    Maybe Optistruct does an 'extrapolation' from integration points?

     

  • Altair Forum User
    Altair Forum User
    Altair Employee
    edited May 2018

    Hi Prakash,

    I tried modal analysis in Hyper mesh with Optistruct interface and in ANSYS workbench. I found very big difference in the natural frequencies. Can u help me through this??

  • Q.Nguyen-Dai
    Q.Nguyen-Dai Altair Community Member
    edited May 2018

    Hi Prakash,

    I tried modal analysis in Hyper mesh with Optistruct interface and in ANSYS workbench. I found very big difference in the natural frequencies. Can u help me through this??

     

    Optistruct is GOOD software.

    ANSYS is GOOD software too.

    So if you found 'very big difference', I'm sure you make 'very big mistake' :)/emoticons/default_smile.png' srcset='/emoticons/smile@2x.png 2x' title=':)' width='20' />

  • Altair Forum User
    Altair Forum User
    Altair Employee
    edited May 2018

    Hi @Rohit Ch

     

    Can you share result files from both solvers and also OptiStruct input file

  • Altair Forum User
    Altair Forum User
    Altair Employee
    edited May 2018

    Thanks Nguyen-Dai for your valuable information.

     

    I already know that something wrong input i have given. But i am unable to find the same. So i was just asking for the different possibilities of mistakes.

  • Altair Forum User
    Altair Forum User
    Altair Employee
    edited May 2018

    Hi @Rohit Ch

     

    Can you share result files from both solvers and also OptiStruct input file 

    Sorry I cant share the result file as it is confidential. But i can share values of natural frequencies.

    in Optistuct first natural frequency is 13.87 Hz

    in ANSYS wb first natural frequency is 594 Hz

     

    So can u plz tell me what are the possibilities of wrong input data??

  • Altair Forum User
    Altair Forum User
    Altair Employee
    edited May 2018

    in ANSYS wb first natural frequency is 594 Hz

     

    Did you validate the result?

     

    First, we need to investigate from which solver is the error is coming from.

     

    WIth just results values, it is difficult to come to a conclusion.