This investigation evaluates, by the dislocation method, the dynamic stress intensity factors of cracked orthotropic half-plane and functionally graded material coating of a coating- substrate material due to the acti...This investigation evaluates, by the dislocation method, the dynamic stress intensity factors of cracked orthotropic half-plane and functionally graded material coating of a coating- substrate material due to the action of anti-plane traction on the crack surfaces. First, by using the complex Fourier transform, the dislocation problem can be solved and the stress fields are obtained with Cauchy singularity at the location of dislocation. The dislocation solution is utilized to derive integral equations for multiple interacting cracks in the orthotropic half-plane with functionally graded orthotropic coating. Several examples are solved and dynamic stress intensity factors are obtained.展开更多
The behaviors of an interface crack between dissimilar orthotropic elastic half-planes subjected to uniform tension was reworked by use of the Schmidt method.By use of the Fourier transform,the problem can be solved w...The behaviors of an interface crack between dissimilar orthotropic elastic half-planes subjected to uniform tension was reworked by use of the Schmidt method.By use of the Fourier transform,the problem can be solved with the help of two pairs of dual integral equations,of which the unknown variables are the jumps of the displacements across the crack surfaces.Numerical examples are provided for the stress intensity factors of the cracks.Contrary to the previous solution of the interface crack,it is found that the stress singularity of the present interface crack solution is of the same nature as that for the ordinary crack in homogeneous materials.When the materials from the two half planes are the same,an exact solution can be otained.展开更多
Vertically oriented carbon structures constructed from low-dimen-sional carbon materials are ideal frameworks for high-performance thermal inter-face materials(TIMs).However,improving the interfacial heat-transfer eff...Vertically oriented carbon structures constructed from low-dimen-sional carbon materials are ideal frameworks for high-performance thermal inter-face materials(TIMs).However,improving the interfacial heat-transfer efficiency of vertically oriented carbon structures is a challenging task.Herein,an orthotropic three-dimensional(3D)hybrid carbon network(VSCG)is fabricated by depositing vertically aligned carbon nanotubes(VACNTs)on the surface of a horizontally oriented graphene film(HOGF).The interfacial interaction between the VACNTs and HOGF is then optimized through an annealing strategy.After regulating the orientation structure of the VACNTs and filling the VSCG with polydimethylsi-loxane(PDMS),VSCG/PDMS composites with excellent 3D thermal conductive properties are obtained.The highest in-plane and through-plane thermal conduc-tivities of the composites are 113.61 and 24.37 W m^(-1)K^(-1),respectively.The high contact area of HOGF and good compressibility of VACNTs imbue the VSCG/PDMS composite with low thermal resistance.In addition,the interfacial heat-transfer efficiency of VSCG/PDMS composite in the TIM performance was improved by 71.3%compared to that of a state-of-the-art thermal pad.This new structural design can potentially realize high-performance TIMs that meet the need for high thermal conductivity and low contact thermal resistance in interfacial heat-transfer processes.展开更多
Black phosphorus nanotubes(BPNTs)may have good properties and potential applications.Determining thevibration property of BPNTs is essential for gaining insight into the mechanical behaviour of BPNTs and designingopti...Black phosphorus nanotubes(BPNTs)may have good properties and potential applications.Determining thevibration property of BPNTs is essential for gaining insight into the mechanical behaviour of BPNTs and designingoptimized nanodevices.In this paper,the mechanical behaviour and vibration property of BPNTs are studied viaorthotropic cylindrical shell model and molecular dynamics(MD)simulation.The vibration frequencies of twochiral BPNTs are analysed systematically.According to the results of MD calculations,it is revealed that thenatural frequencies of two BPNTs with approximately equal sizes are unequal at each order,and that the naturalfrequencies of armchair BPNTs are higher than those of zigzag BPNTs.In addition,an armchair BPNTs witha stable structure is considered as the object of research,and the vibration frequencies of BPNTs of differentsizes are analysed.When comparing the MD results,it is found that both the isotropic cylindrical shell modeland orthotropic cylindrical shell model can better predict the thermal vibration of the lower order modes of thelonger BPNTs better.However,for the vibration of shorter and thinner BPNTs,the prediction of the orthotropiccylindrical shell model is obviously superior to the isotropic shell model,thereby further proving the validity ofthe shell model that considers orthotropic for BPNTs.展开更多
To analyze the stress state of steel orthotropic deck pavement and provide reference for the design of the overlay, the inner stress state and strain distribution of surfacing under the load of the deformation of the ...To analyze the stress state of steel orthotropic deck pavement and provide reference for the design of the overlay, the inner stress state and strain distribution of surfacing under the load of the deformation of the whole bridge structure and tyre load are analyzed by the finite element method of submodeling. Influence of surfacing modulus on the strain state of the overlay is analyzed for the purpose of the optimal design of the overlay structure. Analysis results show that the deformation of the whole bridge structure has no evident influence on the stress state of the overlay. The key factor of the overlay design is the transverse tensile strain in the overlay above the upper edge of web plate of rib. The stress state of the overlay is influenced evidently by the modulus of rigidity transform overlay. And the stress state of the overlay can be optimized and lowered by increasing the modulus and thickness of rigidity transform overlay, The fatigue test has been done to evaluate the fatigue performance and modulus of different deck pavement materials such as epoxy asphalt, SBS modified asphalt, rosphalt asphalt which can provide reference for deck pavement structure design.展开更多
In this paper three important characteristics in piezoresistance for the orthotropic material are given and proved theoretically:(1) The piezoresistance on the principal axis of an orthotropic material is independent ...In this paper three important characteristics in piezoresistance for the orthotropic material are given and proved theoretically:(1) The piezoresistance on the principal axis of an orthotropic material is independent of shear strains/stresses, but correlated with the normal strains/stresses only;(2) On the principal axis of material, following relations between piezoconductivity and piezoresistivity exist η iikk =-(γ ii ) -2 ξ iikk =-(ρ ii ) 2ξ iikk λ iikk =-(γ ii ) -2 χ iikk =-(ρ ii ) 2χ iikk (3) A laminate composed of orthotropic laminae in different orientations is orthotropic for its average/effective properties.展开更多
The purpose of the present work is to study the buckling problem with plate/shell topology optimization of orthotropic material. A model of buckling topology optimization is established based on the independent, conti...The purpose of the present work is to study the buckling problem with plate/shell topology optimization of orthotropic material. A model of buckling topology optimization is established based on the independent, continuous, and mapping method, which considers structural mass as objective and buckling critical loads as constraints. Firstly, composite exponential function (CEF) and power function (PF) as filter functions are introduced to recognize the element mass, the element stiffness matrix, and the element geometric stiffness matrix. The filter functions of the orthotropic material stiffness are deduced. Then these filter functions are put into buckling topology optimization of a differential equation to analyze the design sensitivity. Furthermore, the buckling constraints are approximately expressed as explicit functions with respect to the design variables based on the first-order Taylor expansion. The objective function is standardized based on the second-order Taylor expansion. Therefore, the optimization model is translated into a quadratic program. Finally, the dual sequence quadratic programming (DSQP) algorithm and the global convergence method of moving asymptotes algorithm with two different filter functions (CEF and PF) are applied to solve the optimal model. Three numerical results show that DSQP&CEF has the best performance in the view of structural mass and discretion.展开更多
This paper analyses the bending of rectangular orthotropic plates on a Winkler elastic foundation.Appropriate definition of symplectic inner product and symplectic space formed by generalized displacements establish d...This paper analyses the bending of rectangular orthotropic plates on a Winkler elastic foundation.Appropriate definition of symplectic inner product and symplectic space formed by generalized displacements establish dual variables and dual equations in the symplectic space.The operator matrix of the equation set is proven to be a Hamilton operator matrix.Separation of variables and eigenfunction expansion creates a basis for analyzing the bending of rectangular orthotropic plates on Winkler elastic foundation and obtaining solutions for plates having any boundary condition.There is discussion of symplectic eigenvalue problems of orthotropic plates under two typical boundary conditions,with opposite sides simply supported and opposite sides clamped.Transcendental equations of eigenvalues and symplectic eigenvectors in analytical form given.Analytical solutions using two examples are presented to show the use of the new methods described in this paper.To verify the accuracy and convergence,a fully simply supported plate that is fully and simply supported under uniformly distributed load is used to compare the classical Navier method,the Levy method and the new method.Results show that the new technique has good accuracy and better convergence speed than other methods,especially in relation to internal forces.A fully clamped rectangular plate on Winkler foundation is solved to validate application of the new methods,with solutions compared to those produced by the Galerkin method.展开更多
When a body consists completely or even partly of viscoelastic materials, its response under static loading will be time-dependent. The adhesives used to glue together single plies in laminates usually exhibit a certa...When a body consists completely or even partly of viscoelastic materials, its response under static loading will be time-dependent. The adhesives used to glue together single plies in laminates usually exhibit a certain viscoelastic characteristic in a high temperature environment. In this paper, a laminated orthotropic rectangular plate with viscoelastic interfaces, described by the Kelvin-Voigt model, is considered. A power series expansion technique is adopted to approximate the time-variation of various field quantities. Results indicate that the response of the laminated plate with viscoelastic interfaces changes remarkably with time, and is much different from that of a plate with spring-like or viscous interfaces.展开更多
On the basis of the existing fundamental solutions ofdisplacements, further improvement is made, and then the generalfundamental solutions of both plane elastic and plane plasticproblems for ortho- tropic materials ar...On the basis of the existing fundamental solutions ofdisplacements, further improvement is made, and then the generalfundamental solutions of both plane elastic and plane plasticproblems for ortho- tropic materials are obtained. Two parametersbased on material constants a_1, a_2 are used to derive the rele-vant expressions in a real variable form. Additionally, an analyticalmethod of solving the singular integral for the internal stresses isintroduced, and the corresponding result are given. If a_1=a_2=1, allthe expres- sions obtained for orthotropy can be reduced to thecorresponding ones for isotropy. Because all these expres- sions andresults can be directly used for both isotropic problems andorthotropic problems, it is convenient to use them in engineeringwith the boundary element method (BEM).展开更多
A nonlocal continuum orthotropic plate model is proposed to study the vibration behavior of single-layer graphene sheets (SLGSs) using an analytical symplectic approach. A Hamiltonian system is established by introduc...A nonlocal continuum orthotropic plate model is proposed to study the vibration behavior of single-layer graphene sheets (SLGSs) using an analytical symplectic approach. A Hamiltonian system is established by introducing a total unknown vector consisting of the displacement amplitude, rotation angle, shear force, and bending moment. The high-order governing differential equation of the vibration of SLGSs is transformed into a set of ordinary differential equations in symplectic space. Exact solutions for free vibration are obtianed by the method of separation of variables without any trial shape functions and can be expanded in series of symplectic eigenfunctions. Analytical frequency equations are derived for all six possible boundary conditions. Vibration modes are expressed in terms of the symplectic eigenfunctions. In the numerical examples, comparison is presented to verify the accuracy of the proposed method. Comprehensive numerical examples for graphene sheets with Levy-type boundary conditions are given. A parametric study of the natural frequency is also included.展开更多
In this paper, the evaluation of stress intensity factor of plane crack problems for orthotropic plate of equal-parameter is investigated using a fractal two-level finite element method (F2LFEM). The general solutio...In this paper, the evaluation of stress intensity factor of plane crack problems for orthotropic plate of equal-parameter is investigated using a fractal two-level finite element method (F2LFEM). The general solution of an orthotropic crack problem is obtained by assimilating the problem with isotropic crack problem, and is employed as the global interpolation function in F2LFEM. In the neighborhood of crack tip of the crack plate, the fractal geometry concept is introduced to achieve the similar meshes having similarity ratio less than one and generate an infinitesimal mesh so that the relationship between the stiffness matrices of two adjacent layers is equal. A large number of degrees of freedom around the crack tip are transformed to a small set of generalized coordinates. Numerical examples show that this method is efficient and accurate in evaluating the stress intensity factor (SIF).展开更多
Orthotropic materials weakened by a doubly periodic array of cracks under far-field antiplane shear are investigated, where the fundamental cell contains four cracks of unequal size. By applying the mapping technique,...Orthotropic materials weakened by a doubly periodic array of cracks under far-field antiplane shear are investigated, where the fundamental cell contains four cracks of unequal size. By applying the mapping technique, the elliptical function theory and the theory of analytical function boundary value problems, a closed form solution of the whole-field stress is obtained. The exact formulae for the stress intensity factor at the crack tip and the effective antiplane shear modulus of the cracked orthotropic material are derived. A comparison with the finite element method shows the efficiency and accuracy of the present method. Several illustrative examples are provided, and an interesting phenomenon is observed, that is, the stress intensity factor and the dimensionless effective modulus are independent of the material property for a doubly periodic cracked isotropic material, but depend strongly on the material property for the doubly periodic cracked orthotropic material. Such a phenomenon for antiplane problems is similar to that for in-plane problems. The present solution can provide benchmark results for other numerical and approximate methods.展开更多
In this paper, double dissimilar orthotropic composite materials interfacial crack is studied by constructing new stress functions and employing the method of composite material complex. When the characteristic equati...In this paper, double dissimilar orthotropic composite materials interfacial crack is studied by constructing new stress functions and employing the method of composite material complex. When the characteristic equations' discriminants △1 〉 0 and △2 〉0, the theoretical formula of the stress field and the displacement field near the mode I interface crack tip are derived, indicating that there is no oscillation and interembedding between the interfaces of the crack.展开更多
Adopting the complex function approach, the paper studies the stress intensity factor in orthotropic bi-material interface cracks under mixed loads. With con- sideration of the boundary conditions, a new stress functi...Adopting the complex function approach, the paper studies the stress intensity factor in orthotropic bi-material interface cracks under mixed loads. With con- sideration of the boundary conditions, a new stress function is introduced to transform the problem of bi-material interface crack into a boundary value problem of partial dif- ferential equations. Two sets of non-homogeneous linear equations with 16 unknowns are constructed. By solving the equations, the expressions for the real bi-material elastic constant εt and the real stress singularity exponents λt are obtained with the bi-material engineering parameters satisfying certain conditions. By the uniqueness theorem of limit, undetermined coefficients are determined, and thus the bi-material stress intensity factor in mixed cracks is obtained. The bi-material stress intensity factor characterizes features of mixed cracks. When orthotropic bi-materials are of the same material, the degenerate solution to the stress intensity factor in mixed bi-material interface cracks is in complete agreement with the present classic conclusion. The relationship between the bi-material stress intensity factor and the ratio of bi-material shear modulus and the relationship be- tween the bi-material stress intensity factor and the ratio of bi-material Young's modulus are given in the numerical analysis.展开更多
This article deals with solutions of transient vibration of a rectangular viscoelastic orthotropic thin 2D plate for particular deformation models according to Flu¨gge and Timoshenko-Mindlin.The linear model,a ge...This article deals with solutions of transient vibration of a rectangular viscoelastic orthotropic thin 2D plate for particular deformation models according to Flu¨gge and Timoshenko-Mindlin.The linear model,a general standard viscoelastic body,of the rheologic properties of a viscoelastic material was applied.The time and coordinate curves of the basic quantities displacement,rotation,velocity,stress and deformation are compared.The results obtained by an approximate analytic method are compared with numerical results for 3D plate generated by FEM application and with experimental investigation.展开更多
In this paper, ihe probleins of nonlinear unsymmeirucal bending for cylindricallyorthotropic circular plale are sludied by using “ the method of two-variabie” ̄[1], and theuniformly valid asympiotic soluiions of Nth...In this paper, ihe probleins of nonlinear unsymmeirucal bending for cylindricallyorthotropic circular plale are sludied by using “ the method of two-variabie” ̄[1], and theuniformly valid asympiotic soluiions of Nth-order .lor ε_1 and Mth-order for ε_2 areobiained展开更多
A theoretical method for analyzing the axisymmetric plane strain elastodynamic problem of a non-homogeneous orthotropic hollow cylinder is developed. Firstly, a new dependent variable is introduced to rewrite the gove...A theoretical method for analyzing the axisymmetric plane strain elastodynamic problem of a non-homogeneous orthotropic hollow cylinder is developed. Firstly, a new dependent variable is introduced to rewrite the governing equation, the boundary conditions and the initial conditions. Secondly, a special function is introduced to transform the inhomogeneous boundary conditions to homogeneous ones. By virtue of the orthogonal expansion technique, the equation with respect to the time variable is derived, of which the solution can be obtained. The displacement solution is finally obtained, which can be degenerated in a rather straightforward way into the solution for a homogeneous orthotropic hollow cylinder and isotropic solid cylinder as well as that for a non-homogeneous isotropic hollow cylinder. Using the present method, integral transform can be avoided and it can be used for hollow cylinders with arbitrary thickness and subjected to arbitrary dynamic loads. Numerical results are presented for a non-homogeneous orthotropic hollow cylinder subjected to dynamic internal pressure.展开更多
The dynamic behavior of a rectangular crack in a three-dimensional (3D) orthotropic elastic medium is investigated under a harmonic stress wave based on the non-local theory. The two-dimensional (2D) Fourier trans...The dynamic behavior of a rectangular crack in a three-dimensional (3D) orthotropic elastic medium is investigated under a harmonic stress wave based on the non-local theory. The two-dimensional (2D) Fourier transform is applied, and the mixed- boundary value problems are converted into three pairs of dual integral equations with the unknown variables being the displacement jumps across the crack surfaces. The effects of the geometric shape of the rectangular crack, the circular frequency of the incident waves, and the lattice parameter of the orthotropic elastic medium on the dynamic stress field near the crack edges are analyzed. The present solution exhibits no stress singularity at the rectangular crack edges, and the dynamic stress field near the rectangular crack edges is finite.展开更多
The analysis presented here is to study the effect of non-homogeneity on thermally induced vibration of orthotropic visco-elastic rectangular plate of linearly varying thickness. Thermal vibrational behavior of non-ho...The analysis presented here is to study the effect of non-homogeneity on thermally induced vibration of orthotropic visco-elastic rectangular plate of linearly varying thickness. Thermal vibrational behavior of non-homogeneous rectangular plates of variable thickness having clamped boundary conditions on all the four edges is studied. For non–homogeneity of the plate material, density is assumed to vary linearly in one direction. Using the method of separation of variables, the governing differential equation is solved. An approximate but quite convenient frequency equation is derived by using Rayleigh-Ritz technique with a two-term deflection function. Time period and deflection at different points for the first two modes of vibration are calculated for various values of temperature gradients, non- homogeneity constant, taper constant and aspect ratio. Comparison studies have been carried out with non-homogeneous visco-elastic rectangular plate to establish the accuracy and versatility.展开更多
文摘This investigation evaluates, by the dislocation method, the dynamic stress intensity factors of cracked orthotropic half-plane and functionally graded material coating of a coating- substrate material due to the action of anti-plane traction on the crack surfaces. First, by using the complex Fourier transform, the dislocation problem can be solved and the stress fields are obtained with Cauchy singularity at the location of dislocation. The dislocation solution is utilized to derive integral equations for multiple interacting cracks in the orthotropic half-plane with functionally graded orthotropic coating. Several examples are solved and dynamic stress intensity factors are obtained.
文摘The behaviors of an interface crack between dissimilar orthotropic elastic half-planes subjected to uniform tension was reworked by use of the Schmidt method.By use of the Fourier transform,the problem can be solved with the help of two pairs of dual integral equations,of which the unknown variables are the jumps of the displacements across the crack surfaces.Numerical examples are provided for the stress intensity factors of the cracks.Contrary to the previous solution of the interface crack,it is found that the stress singularity of the present interface crack solution is of the same nature as that for the ordinary crack in homogeneous materials.When the materials from the two half planes are the same,an exact solution can be otained.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52130303,52327802,52303101,52173078,51973158)the China Postdoctoral Science Foundation(2023M732579)+2 种基金Young Elite Scientists Sponsorship Program by CAST(No.2022QNRC001)National Key R&D Program of China(No.2022YFB3805702)Joint Funds of Ministry of Education(8091B032218).
文摘Vertically oriented carbon structures constructed from low-dimen-sional carbon materials are ideal frameworks for high-performance thermal inter-face materials(TIMs).However,improving the interfacial heat-transfer efficiency of vertically oriented carbon structures is a challenging task.Herein,an orthotropic three-dimensional(3D)hybrid carbon network(VSCG)is fabricated by depositing vertically aligned carbon nanotubes(VACNTs)on the surface of a horizontally oriented graphene film(HOGF).The interfacial interaction between the VACNTs and HOGF is then optimized through an annealing strategy.After regulating the orientation structure of the VACNTs and filling the VSCG with polydimethylsi-loxane(PDMS),VSCG/PDMS composites with excellent 3D thermal conductive properties are obtained.The highest in-plane and through-plane thermal conduc-tivities of the composites are 113.61 and 24.37 W m^(-1)K^(-1),respectively.The high contact area of HOGF and good compressibility of VACNTs imbue the VSCG/PDMS composite with low thermal resistance.In addition,the interfacial heat-transfer efficiency of VSCG/PDMS composite in the TIM performance was improved by 71.3%compared to that of a state-of-the-art thermal pad.This new structural design can potentially realize high-performance TIMs that meet the need for high thermal conductivity and low contact thermal resistance in interfacial heat-transfer processes.
基金supported by the National Science Fund for Distin-guished Young Scholars(Grants No.11925205)the National Natural Science Foundation of China(Grant Nos.51921003 and U2341230).
文摘Black phosphorus nanotubes(BPNTs)may have good properties and potential applications.Determining thevibration property of BPNTs is essential for gaining insight into the mechanical behaviour of BPNTs and designingoptimized nanodevices.In this paper,the mechanical behaviour and vibration property of BPNTs are studied viaorthotropic cylindrical shell model and molecular dynamics(MD)simulation.The vibration frequencies of twochiral BPNTs are analysed systematically.According to the results of MD calculations,it is revealed that thenatural frequencies of two BPNTs with approximately equal sizes are unequal at each order,and that the naturalfrequencies of armchair BPNTs are higher than those of zigzag BPNTs.In addition,an armchair BPNTs witha stable structure is considered as the object of research,and the vibration frequencies of BPNTs of differentsizes are analysed.When comparing the MD results,it is found that both the isotropic cylindrical shell modeland orthotropic cylindrical shell model can better predict the thermal vibration of the lower order modes of thelonger BPNTs better.However,for the vibration of shorter and thinner BPNTs,the prediction of the orthotropiccylindrical shell model is obviously superior to the isotropic shell model,thereby further proving the validity ofthe shell model that considers orthotropic for BPNTs.
文摘To analyze the stress state of steel orthotropic deck pavement and provide reference for the design of the overlay, the inner stress state and strain distribution of surfacing under the load of the deformation of the whole bridge structure and tyre load are analyzed by the finite element method of submodeling. Influence of surfacing modulus on the strain state of the overlay is analyzed for the purpose of the optimal design of the overlay structure. Analysis results show that the deformation of the whole bridge structure has no evident influence on the stress state of the overlay. The key factor of the overlay design is the transverse tensile strain in the overlay above the upper edge of web plate of rib. The stress state of the overlay is influenced evidently by the modulus of rigidity transform overlay. And the stress state of the overlay can be optimized and lowered by increasing the modulus and thickness of rigidity transform overlay, The fatigue test has been done to evaluate the fatigue performance and modulus of different deck pavement materials such as epoxy asphalt, SBS modified asphalt, rosphalt asphalt which can provide reference for deck pavement structure design.
文摘In this paper three important characteristics in piezoresistance for the orthotropic material are given and proved theoretically:(1) The piezoresistance on the principal axis of an orthotropic material is independent of shear strains/stresses, but correlated with the normal strains/stresses only;(2) On the principal axis of material, following relations between piezoconductivity and piezoresistivity exist η iikk =-(γ ii ) -2 ξ iikk =-(ρ ii ) 2ξ iikk λ iikk =-(γ ii ) -2 χ iikk =-(ρ ii ) 2χ iikk (3) A laminate composed of orthotropic laminae in different orientations is orthotropic for its average/effective properties.
基金supported by the National Natural Science Foundation of China (Grants 11072009, 11172013)the Beijing Education Committee Development Project (Grant SQKM2016100 05001)the Beijing University of Technology Basic Research Fund (Grant 001000514313003)
文摘The purpose of the present work is to study the buckling problem with plate/shell topology optimization of orthotropic material. A model of buckling topology optimization is established based on the independent, continuous, and mapping method, which considers structural mass as objective and buckling critical loads as constraints. Firstly, composite exponential function (CEF) and power function (PF) as filter functions are introduced to recognize the element mass, the element stiffness matrix, and the element geometric stiffness matrix. The filter functions of the orthotropic material stiffness are deduced. Then these filter functions are put into buckling topology optimization of a differential equation to analyze the design sensitivity. Furthermore, the buckling constraints are approximately expressed as explicit functions with respect to the design variables based on the first-order Taylor expansion. The objective function is standardized based on the second-order Taylor expansion. Therefore, the optimization model is translated into a quadratic program. Finally, the dual sequence quadratic programming (DSQP) algorithm and the global convergence method of moving asymptotes algorithm with two different filter functions (CEF and PF) are applied to solve the optimal model. Three numerical results show that DSQP&CEF has the best performance in the view of structural mass and discretion.
基金supported by the National Natural Science Foundation of China (10772039 and 10632030)the National Basic Research Program of China (973 Program) (2010CB832704)
文摘This paper analyses the bending of rectangular orthotropic plates on a Winkler elastic foundation.Appropriate definition of symplectic inner product and symplectic space formed by generalized displacements establish dual variables and dual equations in the symplectic space.The operator matrix of the equation set is proven to be a Hamilton operator matrix.Separation of variables and eigenfunction expansion creates a basis for analyzing the bending of rectangular orthotropic plates on Winkler elastic foundation and obtaining solutions for plates having any boundary condition.There is discussion of symplectic eigenvalue problems of orthotropic plates under two typical boundary conditions,with opposite sides simply supported and opposite sides clamped.Transcendental equations of eigenvalues and symplectic eigenvectors in analytical form given.Analytical solutions using two examples are presented to show the use of the new methods described in this paper.To verify the accuracy and convergence,a fully simply supported plate that is fully and simply supported under uniformly distributed load is used to compare the classical Navier method,the Levy method and the new method.Results show that the new technique has good accuracy and better convergence speed than other methods,especially in relation to internal forces.A fully clamped rectangular plate on Winkler foundation is solved to validate application of the new methods,with solutions compared to those produced by the Galerkin method.
基金Project supported by the National Natural Science Foundation of China (No. 10432030) and NCET.
文摘When a body consists completely or even partly of viscoelastic materials, its response under static loading will be time-dependent. The adhesives used to glue together single plies in laminates usually exhibit a certain viscoelastic characteristic in a high temperature environment. In this paper, a laminated orthotropic rectangular plate with viscoelastic interfaces, described by the Kelvin-Voigt model, is considered. A power series expansion technique is adopted to approximate the time-variation of various field quantities. Results indicate that the response of the laminated plate with viscoelastic interfaces changes remarkably with time, and is much different from that of a plate with spring-like or viscous interfaces.
文摘On the basis of the existing fundamental solutions ofdisplacements, further improvement is made, and then the generalfundamental solutions of both plane elastic and plane plasticproblems for ortho- tropic materials are obtained. Two parametersbased on material constants a_1, a_2 are used to derive the rele-vant expressions in a real variable form. Additionally, an analyticalmethod of solving the singular integral for the internal stresses isintroduced, and the corresponding result are given. If a_1=a_2=1, allthe expres- sions obtained for orthotropy can be reduced to thecorresponding ones for isotropy. Because all these expres- sions andresults can be directly used for both isotropic problems andorthotropic problems, it is convenient to use them in engineeringwith the boundary element method (BEM).
基金support of the National Natural Science Foundation of China (Grant 11672054)the Research Grant Council of Hong Kong (11215415)the National Basic Research Program of China (973 Program) (Grant 2014CB046803)
文摘A nonlocal continuum orthotropic plate model is proposed to study the vibration behavior of single-layer graphene sheets (SLGSs) using an analytical symplectic approach. A Hamiltonian system is established by introducing a total unknown vector consisting of the displacement amplitude, rotation angle, shear force, and bending moment. The high-order governing differential equation of the vibration of SLGSs is transformed into a set of ordinary differential equations in symplectic space. Exact solutions for free vibration are obtianed by the method of separation of variables without any trial shape functions and can be expanded in series of symplectic eigenfunctions. Analytical frequency equations are derived for all six possible boundary conditions. Vibration modes are expressed in terms of the symplectic eigenfunctions. In the numerical examples, comparison is presented to verify the accuracy of the proposed method. Comprehensive numerical examples for graphene sheets with Levy-type boundary conditions are given. A parametric study of the natural frequency is also included.
文摘In this paper, the evaluation of stress intensity factor of plane crack problems for orthotropic plate of equal-parameter is investigated using a fractal two-level finite element method (F2LFEM). The general solution of an orthotropic crack problem is obtained by assimilating the problem with isotropic crack problem, and is employed as the global interpolation function in F2LFEM. In the neighborhood of crack tip of the crack plate, the fractal geometry concept is introduced to achieve the similar meshes having similarity ratio less than one and generate an infinitesimal mesh so that the relationship between the stiffness matrices of two adjacent layers is equal. A large number of degrees of freedom around the crack tip are transformed to a small set of generalized coordinates. Numerical examples show that this method is efficient and accurate in evaluating the stress intensity factor (SIF).
基金supported by the National Natural Science Foundation of China (No.10672008).
文摘Orthotropic materials weakened by a doubly periodic array of cracks under far-field antiplane shear are investigated, where the fundamental cell contains four cracks of unequal size. By applying the mapping technique, the elliptical function theory and the theory of analytical function boundary value problems, a closed form solution of the whole-field stress is obtained. The exact formulae for the stress intensity factor at the crack tip and the effective antiplane shear modulus of the cracked orthotropic material are derived. A comparison with the finite element method shows the efficiency and accuracy of the present method. Several illustrative examples are provided, and an interesting phenomenon is observed, that is, the stress intensity factor and the dimensionless effective modulus are independent of the material property for a doubly periodic cracked isotropic material, but depend strongly on the material property for the doubly periodic cracked orthotropic material. Such a phenomenon for antiplane problems is similar to that for in-plane problems. The present solution can provide benchmark results for other numerical and approximate methods.
基金the Natural Science Foundation of Shanxi Province(No.2007011008)
文摘In this paper, double dissimilar orthotropic composite materials interfacial crack is studied by constructing new stress functions and employing the method of composite material complex. When the characteristic equations' discriminants △1 〉 0 and △2 〉0, the theoretical formula of the stress field and the displacement field near the mode I interface crack tip are derived, indicating that there is no oscillation and interembedding between the interfaces of the crack.
基金supported by the National Key Basic Research Program of China(973 Program)(No.2009CB724201)the Science and Technology Major Project of the Ministry of Education of China(No.208022)+1 种基金the Postgraduate Scientific and Technological Innovation Project of Taiyuan University of Science and Technology(No.20125027)the Scientific Research Funds for Doctoral Students of Taiyuan University of Science and Technology(No.20122005)
文摘Adopting the complex function approach, the paper studies the stress intensity factor in orthotropic bi-material interface cracks under mixed loads. With con- sideration of the boundary conditions, a new stress function is introduced to transform the problem of bi-material interface crack into a boundary value problem of partial dif- ferential equations. Two sets of non-homogeneous linear equations with 16 unknowns are constructed. By solving the equations, the expressions for the real bi-material elastic constant εt and the real stress singularity exponents λt are obtained with the bi-material engineering parameters satisfying certain conditions. By the uniqueness theorem of limit, undetermined coefficients are determined, and thus the bi-material stress intensity factor in mixed cracks is obtained. The bi-material stress intensity factor characterizes features of mixed cracks. When orthotropic bi-materials are of the same material, the degenerate solution to the stress intensity factor in mixed bi-material interface cracks is in complete agreement with the present classic conclusion. The relationship between the bi-material stress intensity factor and the ratio of bi-material shear modulus and the relationship be- tween the bi-material stress intensity factor and the ratio of bi-material Young's modulus are given in the numerical analysis.
基金supported by the grant project GA CR 101/07/0946
文摘This article deals with solutions of transient vibration of a rectangular viscoelastic orthotropic thin 2D plate for particular deformation models according to Flu¨gge and Timoshenko-Mindlin.The linear model,a general standard viscoelastic body,of the rheologic properties of a viscoelastic material was applied.The time and coordinate curves of the basic quantities displacement,rotation,velocity,stress and deformation are compared.The results obtained by an approximate analytic method are compared with numerical results for 3D plate generated by FEM application and with experimental investigation.
文摘In this paper, ihe probleins of nonlinear unsymmeirucal bending for cylindricallyorthotropic circular plale are sludied by using “ the method of two-variabie” ̄[1], and theuniformly valid asympiotic soluiions of Nth-order .lor ε_1 and Mth-order for ε_2 areobiained
基金The project supported by the National Natural Science Foundation of China (10172075 and 10002016)
文摘A theoretical method for analyzing the axisymmetric plane strain elastodynamic problem of a non-homogeneous orthotropic hollow cylinder is developed. Firstly, a new dependent variable is introduced to rewrite the governing equation, the boundary conditions and the initial conditions. Secondly, a special function is introduced to transform the inhomogeneous boundary conditions to homogeneous ones. By virtue of the orthogonal expansion technique, the equation with respect to the time variable is derived, of which the solution can be obtained. The displacement solution is finally obtained, which can be degenerated in a rather straightforward way into the solution for a homogeneous orthotropic hollow cylinder and isotropic solid cylinder as well as that for a non-homogeneous isotropic hollow cylinder. Using the present method, integral transform can be avoided and it can be used for hollow cylinders with arbitrary thickness and subjected to arbitrary dynamic loads. Numerical results are presented for a non-homogeneous orthotropic hollow cylinder subjected to dynamic internal pressure.
基金Project supported by the National Natural Science Foundation of China(Nos.11272105 and 11572101)
文摘The dynamic behavior of a rectangular crack in a three-dimensional (3D) orthotropic elastic medium is investigated under a harmonic stress wave based on the non-local theory. The two-dimensional (2D) Fourier transform is applied, and the mixed- boundary value problems are converted into three pairs of dual integral equations with the unknown variables being the displacement jumps across the crack surfaces. The effects of the geometric shape of the rectangular crack, the circular frequency of the incident waves, and the lattice parameter of the orthotropic elastic medium on the dynamic stress field near the crack edges are analyzed. The present solution exhibits no stress singularity at the rectangular crack edges, and the dynamic stress field near the rectangular crack edges is finite.
文摘The analysis presented here is to study the effect of non-homogeneity on thermally induced vibration of orthotropic visco-elastic rectangular plate of linearly varying thickness. Thermal vibrational behavior of non-homogeneous rectangular plates of variable thickness having clamped boundary conditions on all the four edges is studied. For non–homogeneity of the plate material, density is assumed to vary linearly in one direction. Using the method of separation of variables, the governing differential equation is solved. An approximate but quite convenient frequency equation is derived by using Rayleigh-Ritz technique with a two-term deflection function. Time period and deflection at different points for the first two modes of vibration are calculated for various values of temperature gradients, non- homogeneity constant, taper constant and aspect ratio. Comparison studies have been carried out with non-homogeneous visco-elastic rectangular plate to establish the accuracy and versatility.