讲义分析成果cma

1、Lesson content:Overview of Failure Criteria Plug-inWorkshop 11: Cooper Hood Failure Criteria Lesson 8: Failure Criteria Plug-in1 hoursOverview of Failure Criteria Plug-in (1/16)User-defined failure criteriaComposites Modeler allows users to extend the failure criteria available in AbaqusOverview of

2、Failure Criteria Plug-in (2/16)User interfaceOverview of Failure Criteria Plug-in (3/16)Set failure properties in Abaqus before running analysisOverview of Failure Criteria Plug-in (4/16)Output can be written to a text fileComposite Failure Results-Date : 2011-06-23 Time : 10:34:05Database : Cowling

3、_pressure_load.odbInstance : COWLING-1Step : Step-1 Frame: 1User : MikeSummary of Failure Analysis Results-Failure Criterion : TsaiWuFailure Basis : StressLoadcase : Subcase : Result : Stress TensorElement : 1217Section Point : 1Integration Point : 1Reserve Factors : 1.34200032E-004Material Values U

4、sed-Material Name : T300Xt Xc Yt Yc S12 S23 S31 F12 1.00000000E+007 1.00000000E+006 5.00000000E+004 5.00000000E+004 6.00000000E+005 0.00000000E+000 0.00000000E+000 5.00000000E-001 Summary of Element Failure Analysis Results-Element Sec.Point Int.Point Material Reserve Factors 1 1 1 1 7.64311384E-003

5、 2 45 1 1 2.61123385E-003 3 1 1 1 5.23332832E-003 4 1 1 1 3.23418365E-003Overview of Failure Criteria Plug-in (5/16)Output is stored in the current ODB fileOverview of Failure Criteria Plug-in (6/16)Plot of reserve factors: Tsai-WuOverview of Failure Criteria Plug-in (7/16)Adding your own failure cr

6、iteriaNeed to use C+Take and modify the template files provided. Add one function Run() to calculate the failure output results based on the inputs (stress/strain/allowables)Possibly add checkAllowable() to validate inputsOverview of Failure Criteria Plug-in (8/16)Interface functionsclass CSltAddinF

7、ailureLaminaPuck : public ISltAddinFailureLaminapublic:CSltAddinFailureLaminaPuck(void) ;virtual CSltAddinFailureLaminaPuck(void) ;/ Implementation of the ISltIUnknown interfacevirtual int AddRef( ) ;virtual int Release( ) ;virtual int QueryInterface( const char *pszID, void *ppInteface ) ;/ Get inp

8、ut requirements.virtual int GetNumAllowables(int &iNumAllowables) ;virtual int GetAllowableLabel(int iIndex, char *pszLabel, int iSize) ;/ Get output capabilities.virtual int GetNumResults(int &iNumResults) ;virtual int GetResultLabel(int iResultIndex, char *pszLabel, int iszSize) ;virtual int GetNu

9、mResultValues(int iResultIndex, int &iNumdValues, int &iNumiValues, int &iNumszValues) ;/ Set basis.virtual int SetBasis( EBasis eBasis) ;/ Input methods.virtual int SetAllowables(int iNumAllowables, double *pAllowables) ;virtual int SetLoadInDirection(ELoadDirection eDirection, double dLoad) ;virtu

10、al int CheckAllowables() ;/ Run.virtual int Run() ;/ Reinitialise.virtual int Clear() ;/ Implementation of the ISCFailtureLaminate interface/ Identification methods.virtual int GetName( char *pszName, int iSize ) ;virtual int GetProviderName(char *pszName, int iSize) ;virtual int GetProviderWebsite(

11、char *pszName, int iSize) ;virtual int GetVersion(int &iMajor, int &iMinor, int &iBuild) ;/ Get results.virtual int GetdResult(int iResultIndex, int idValue, double &dResultValue) ;virtual int GetiResult(int iResultIndex, int iiValue, int &iResultValue) ;virtual int GetszResult(int iResultIndex, int

12、 iszValue,char *pszResultValue, int iszSize) ;/ Message methods.virtual int GetNumMessages( int &iNumMessages ) ;virtual int GetMessageLength( int iMessageIndex, int &iMessageLength ) ;virtual int GetMessageString(int iMessageIndex, EMessageType &eMessageType, char *pszMessage, int iMessageLength) ;

13、Overview of Failure Criteria Plug-in (9/16)Example code: Puckint CSltAddinFailureLaminaPuck:Run()/ Check the allowable iReturn = CheckAllowables() ;if( iReturn != 0 ) return 1 ; double m_eps = 1.e-8; /small stress value to match for zero stress/retain these symbols to ensure accuracy of e

14、quationsdouble Xt, Xc, R22p, R22m, R12, R22A, p12p, p12m, p22p, p22m ;double sigma1, sigma2, tau21, tau21c, FF, temp ;/strength parametersXt = m_allowables0; /tensile fibre failure strengthXc = m_allowables1; pressive fibre failure strengthR22p = m_allowables2;/R perp+ = transverse tensile strengthR

15、22m = m_allowables3; /R perp- = transverse compressive strengthR12 = m_allowables4; /R perp parallel = shear strength/shape parametersp12p = m_allowables5;/p perp parallel +p12m = m_allowables6;/p perp parallel -p22p = m_allowables7;/p perp perp +/input stressessigma1=m_loadsISltAddinFailureLamina:e

16、XX; /stress along fibre directionsigma2=m_loadsISltAddinFailureLamina:eYY; /stress orthogonal to fibre directiontau21=m_loadsISltAddinFailureLamina:eXY; /shear stress in 1 direction on 2 faceCourtesy: Michael SutcliffeOverview of Failure Criteria Plug-in (10/16)Example code: Puck (contd)Check input

17、data, create temporary variables/Check whether all stresses less than eps - returns m_iresults(FAILURE_MODE_IX) = 0 / and m_dresults(FAILURE_INDEX_IX)=0if (abs(sigma1)=m_eps & abs(sigma2)=m_eps & abs(tau21)=0.)/Mode Am_dresultsFAILURE_MODE_IX=FAILURE_MODE_IFF_A_INT;temp = pow(tau21/R12),2)+pow(1-p12

18、p*R22p/R12),2)*pow(sigma2/R22p),2);m_dresultsFAILURE_INDEX_IX = sqrt(temp)+p12p*sigma2/R12;m_dresultsFAILURE_ANGLE_IX = 0.;m_dresultsRESERVE_FACTOR_IX = GetReserveFactorFromFailureIndex( m_dresultsFAILURE_INDEX_IX ) ;m_iresultsRESERVE_FACTOR_IX = FAILURE_MODE_IFF_A_INT ;else if (abs(tau21/sigma2)=0.

19、)/tensile axial stressFF = sigma1/Xt;if (FF m_dresultsFAILURE_INDEX_IX)m_dresultsFAILURE_MODE_IX = FAILURE_MODE_TENS_INT;m_dresultsFAILURE_INDEX_IX = FF;m_dresultsFAILURE_ANGLE_IX = 0.;m_dresultsRESERVE_FACTOR_IX = GetReserveFactorFromFailureIndex( m_dresultsFAILURE_INDEX_IX ) ;m_iresultsRESERVE_FACTOR_IX = FAILURE_MODE_TENS_INT ;else pressive axial stressFF =-sigma1/Xc;if (FF m_dresultsFAILURE_INDEX_IX)m_dresultsFAILURE_MODE_IX = P_INT;m_dresultsFAILURE_INDEX_IX