Dissertations / Theses on the topic 'Wing stress analysis'

The main aim of the this diploma thesis is the wing design of the UAV aircraft, including the appropriate material choice, calculation of the wing load and also strength analysis. Other goals include the design of the location and volume of fuel tanks, as well as the design of wingspan reduction after landing.

The amphibious fullmetall plane Mermaid is made by company Czech Aircraft Works in Kunovice. Subject of this diploma thesis was involve changes of structure, which pass through during development, to calculation. Within this published work was process: Load calculation of wing Stress analysis of wing Load calculation of aileron Load calculation of flap Load calculation of horizontail tail

In terms of future power generation in UK and Germany, offshore wind is the next big player with 40GW and 32GW capacity planned for installation in both countries respectively by 2030. The latest Round 3 of sites owned by the Crown Estate explore deeper water depths of up to 78m in the Irish Sea. Foundations for offshore wind structures consume around 25% of the total project cost therefore the design of support structures is the subject of this thesis. The current state-of-the-art support structure options available for offshore wind turbines have been outlined in this thesis with an evaluation of the preliminary design of monopile and jacket solutions. This assessment resulted in further studies into the loading acting on a monopile foundation along with research into the fatigue design of multiplanar tubular joints for jacket structures. Mathematical modelling of linear and nonlinear waves combined with the Morison equation was completed to check the effects of breaking waves on a monopile foundation. Results indicated that measured forces were up to a factor of 2.5 times greater than calculated values, which suggests that loads could be under predicted if the effects of breaking are not considered. The theoretical maximum wave height before breaking was then linked to wind speed and a comparison of overturning moments from the two loads was made. Wave loads dominated at water depths of around 30m for lower wind speeds but this depth decreased to around 12m as wind speeds approached cut-out of 25m/s. For deeper water depths and larger capacity turbines, jackets are the preferred design solution. Joint design in FLS is the critical aspect of jacket design with castings often required to provide adequate capacity. A review of stress concentration factors (SCF) for tubular joints indicated that the coded approach, which uses SCF equations for uniplanar joints, could be missing the multiplanar effects. Finite element (FE) modelling of multiplanar tubular joints was completed using ANSYS Workbench to examine the effects of loading in out-of plane braces. Carry-over of stress from the loaded brace of the joint to unloaded neighbouring braces was observed which implies the importance of modelling joints as multiplanar geometries. A parameter study in ANSYS Workbench covering 1806 different geometrical configurations and loads was carried out with a regression of the data to give new sets of SCF equations for multiplanar tubular joints. SCFs from these equations were improved compared to Efthymiou but difficulties were encountered when superimposing the output (including Efthymiou). Further work on the superposition of individual load cases was therefore recommended for future work.

Dissertação (mestrado) - Universidade Federal de Santa Catarina, Centro Tecnológico, Programa de Pós-Graduação em Engenharia Mecânica, Florianópolis, 2016.
Made available in DSpace on 2017-05-23T04:22:08Z (GMT). No. of bitstreams: 1 345248.pdf: 4540254 bytes, checksum: b61753d1576a6ca395e102d19e3c04e0 (MD5) Previous issue date: 2016
Abstract : Research institutes and industry have recently initiated the development of wind turbines for energy production at sites with lower wind speeds. The relatively high electric power production of this type of wind turbine is related to the growing of the swept area. Thus, in order to increase the swept area, the length of the new blades must be greater than the ones of traditional blades for horizontal axis wind turbines. The expansion of the rotors diameter of large wind turbines has some implications, such as the increase of gravity loads and new transportation logistics challenges. To overcome these challenges, advances in the blade technology have led to more efficient structural and aerodynamic design and optimized material usage. In this context, analytical models have to accurately predict aerodynamic loads acting on the blades and calculate the structural stresses developed during the turbine operation. Since composite materials have high strength/weight ratio compared to other structural materials and are flexible with respect to the fabrication process, they are more appropriate for blade applications. This dissertation presents an improvement in the aerodynamic model that takes into account an iterative procedure for determining Reynolds number instead of depending on experimental data for the process of airfoil selection along the blade span that maximize the power extracted from the wind. The multilayer shell theories and the relationship between the aerodynamic bending and torsional moments acting on the blade are presented. The in-plane normal and shear stresses on the thin-walled multilayer blade are determined by using the Shear Flow Theory. A case study is conducted on a 20 MW wind turbine developed in the Energy Research Centre of the Netherlands-ECN. The normal and shear stresses are calculated and the Tsai-Wu criterion is applied for the strength evaluation of the blade made of glass/epoxy. Results obtained with two stacking sequences are presented.

Institutos de pesquisa e a indústria iniciaram recentemente o desenvolvimento de turbinas eólicas para a produção de energia em locais com baixas velocidades de vento. A produção de energia elétrica relativamente alta deste tipo de turbina está relacionada com o crescimento da área varrida pelo rotor. Assim, de forma a aumentar a área varrida, o comprimento das novas pás deve ser maior que o das pás tradicionais de turbinas eólicas de eixo horizontal. O aumento do diâmetro dos rotores das turbinas eólicas tem algumas implicações, tais como o aumento das cargas gravitacionais e desafios no seu transporte. Para superá-los, avanços na tecnologia de pás têm levado a projetos estruturais e aerodinâmicos mais eficientes. Neste contexto, modelos analíticos têm de prever com precisão as cargas aerodinâmicas atuando sobre as pás, além de calcular as tensões estruturais desenvolvidas durante sua operação. Uma vez que materiais compostos têm uma elevada relação resistência/peso em comparação com outros materiais estruturais, além de serem flexíveis no que diz respeito ao processo de fabricação, estes materiais têm sido os mais investigados para aplicações em pás de turbina eólica. Esta dissertação apresenta uma melhoria no modelo aerodinâmico que leva em conta um processo iterativo para a determinação do número de Reynolds, em vez de depender de dados experimentais para o processo de escolha do aerofólio ao longo da envergadura da pá que maximiza a potência extraída do vento. As teorias de cascas multilaminares e a relação entre os momentos aerodinâmicos fletores e de torção que atuam sobre a pá são apresentados. As tensões normais e de cisalhamento atuantes no plano da pá multilaminar de parede fina são determinadas utilizando a Teoria do Fluxo de Cisalhamento. Um estudo de caso é conduzido em uma turbina eólica de 20 MW desenvolvida no Centro de Pesquisa em Energia dos Países Baixos (ECN). As tensões normais e de cisalhamento são calculadas e o critério de Tsai-Wu é aplicado para a avaliação da resistência da pá feita em vidro/epóxi. Os resultados obtidos com duas sequências de empilhamento de lâminas são apresentados.

The main purpose of this thesis is to develop a beam element in order to model the wing of a drone, made of composite materials. The proposed model consists of the framework for the structural design and analysis of long slender beam like structures, e.g., wings, wind turbine blades, and helicopter rotor blades, etc. The main feature consists of the addition of the coupling between axial and bending with torsional effects that may arise when using composite materials and the coupling stemming from the inhomogeneity in cross-sections of any arbitrary geometry. This type of modeling approach allows for an accurate yet computationally inexpensive representation of a general class of beam-like structures. The framework for beam analysis consists of main two parts, cross-sectional analysis of the beam sections and then using this section analysis to build up the finite element model. The cross-sectional analysis is performed in order to predict the structural properties for composite sections, which are used for the beam model. The thesis consists of the model to validate the convergence of the element size required for the cross-sectional analysis. This follows by the validation of the shell models of constant cross-section to assess the performance of the beam elements, including coupling terms. This framework also has the capability of calculating the strains and displacements at various points of the cross-section. Natural frequencies and mode shapes are compared for different cases of increasing complexity with those available in the papers. Then, the framework is used to analyze the wing of a drone and compare the results to a model developed in NASTRAN.
Master of Science

Bed connected support structures such as monopiles are expected to be impractical for water depths greater than 30 m and so there is increasing interest in alternative structure concepts to enable cost-effective deployment of wind and tidal stream turbines. Floating, moored platforms supporting multiple rotors are being considered for this purpose. This thesis investigates the dynamic response of such floating structures, taking into account the coupling between loading due to both turbulent flow and waves and the dynamic response of the system. The performance and loading of a single rotor in steady and quasi-steady flows are quantified with a Blade Element Momentum Theory (BEMT) code. This model is validated for steady flow against published data for two 0.8 m diameter rotors (Bahaj, Batten, et al., 2007; Galloway et al., 2011) and a 0.27 m diameter rotor (Whelan and Stallard, 2011). Time-averaged coefficients of thrust and power measured by experiment in steady turbulent flow were in agreement with BEMT predictions over a range of angular speeds. The standard deviation of force on the rotor is comparable to that on a porous grid for comparable turbulence characteristics. Drag and added mass coefficients are determined for a porous disc forced to oscillate normal to the rotor plane in quiescent flow and in the streamwise axis in turbulent flow. Added mass is negligible for the Keulegan Carpenter number range considered ( less than 1). The drag coefficient in turbulent flow was found to decay exponentially with number, to 2±10% for values greater than 0.5. These coefficients were found to be in good agreement with those for a rotor in the same turbulent flow with disc drag coefficient within 12.5% for less than 0.65. An extreme-value analysis is applied to the measured time-varying thrust due to turbulent flow and turbulent flow with waves to obtain forces with 1%, 0.1% and 0.01% probability of exceedance during operating conditions. The 1% exceedance force in turbulent flow with turbulence intensity of 12% is around 40% greater than the mean thrust. The peak force in turbulent flow with opposing waves was predicted to within 6% by superposition of the extreme force due to turbulence only with a drag force based on the relative wave-induced velocity at hub-height estimated by linear wave theory and with drag coefficient of 2.0. Response of a floating structure in surge and pitch is studied due to both wave- forcing on the platform defined by the linear diffraction code WAMIT and due to loading of the operating turbine defined by a thrust coefficient and drag coefficient. Platform response can either increase or decrease the loading on the rotor and this was dependant on the hydrodynamic characteristics of the support platform. A reduction of the force on the rotor is attained when the phase difference between the wave force on the support and the surface elevation is close to ± and when the damping of the support is increased. For a typical support and for a wave condition with phase difference close to , the 1% rotor forces were reduced by 8% when compared to the force obtained with a rotor supported on a stiff tower.

Structural dynamics is at the center of wind turbine tower design - excessive vibrations can be caused by a wide range of environmental and mechanical sources and can lead to reduced component life due to fatigue, noise, and impaired public perception of system integrity. Furthermore, periodic turbulent wind conditions can cause system resonance resulting in significantly increased structural loads. Structural vibration issues may become exacerbated in small wind applications where the analytical and experimental resources for system verification and optimization are scarce. This study combines several structural analysis techniques and packages them into a novel and integrated form that can be readily used by the small wind community/designer to gain insight into tower/rotor dynamic interaction, system modal characteristics, and to optimize the design for reduced tower loads and cost. The finite element method is used to model the tower structure and can accommodate various configurations including fixed monopole towers, guy-wire supported towers, and gin-pole and strut supported towers. The turbine rotor is modeled using the Equivalent Hinge-Offset blade model and coupled to the tower structure through the use of Lagrange’s Equations. Standard IEC Aeroelastic load cases are evaluated and transient solutions developed using the Modal Superposition Method and Runge-Kutta 4th order numerical integration. Validation is performed through comparisons to theoretical closed form solutions, physical laboratory test results, and peer studies. Finally a case study is performed by using the tool to simulate the Cal Poly Wind Power Research Center Wind Turbine and Tower System. Included in the case study is an optimization for hypothetical guy-wire placement to minimize tower stresses and maximize the tower’s natural frequency.

Photovoltaic (PV) solar panels and trackers represent one of the most common renewable energy technology which converts sunlight radiation into electrical energy. The solar trackers specifically are more complex structures because they involve mechanical devices, a supporting slender structure, and photovoltaic modules mounted and positioned on top of the supporting structure. Solar trackers are mounted on mobile supports or racks, in order to enable the rotation and tilt of the PV which thus maintains their optimum exposure to the incident sunlight. Solar trackers support structures should be designed for wind resistance during the operation and at stow position for its life span and this became a concern considering the new tendency of installing the solar trackers on the rooftop of low-rise or medium-rise buildings. The current research focused on performing site measurements of the wind-induced displacement for a dual-axis solar tracking system installed on the roof of the Mann Parking building of the University of Ottawa, for different azimuth, elevations. The supporting structure of the solar tracker was instrumented with 16 strain gauges and the strains developed in the metal truss members were measured during the months February 2015 and March 2015. The tracker was rotated and tilted at different angles through the duration of the experiment and the strains observed on each structural element were recorded. In order to estimate deflections of the supporting structure for wind speeds higher than the ones measured, a finite element (FE) model of the solar tracker was created and static analysis was performed for different inclinations using the SAP 2000 structural software. The experimental results were in agreement with the FE simulation results as the stresses obtained ranged between 1.02 × 107 Pa and 7.88 × 107 Pa. Lower attack angles between 45° and 60° were found to have significant effect on the elements of the solar tracker irrespective of the wind load magnitude. Operational attack angles between 65° and 75° were found to be safer positions as obtained displacements and stress analysis result showed that the supporting structure of the solar tracker was stable for wind speeds between 0 m/s and 33m/s in Ottawa region

Includes bibliographical references.
The coastal upwelling regions along northwest Africa (the Canary system), southwest Africa (the Benguela system), North America (the California system), and South America (the Peru-Chile system) were studied and compared on a seasonal timescale. A 17 -year NOAA Pathfinder Sea Surface Temperature (SST) dataset with a spatial resolution of 9km was used to describe the large-scale temporal and spatial variability of upwelling within the four regions. An upwelling index (OSST) was derived in order to add to the patterns of upwelling variability described from the sea surface temperature. The upwelling index was also used to describe the similarities and differences between the four upwelling regions. A 10-year ERS wind stress dataset with a spatial resolution of 10 x 10 was used to derive the offshore Ekman Transport in each of the regions. The offshore Ekman Transport was used to supplement the description of SST's in each system. Principal Component Analysis was used to investigate the variance structure of the anomalies of the sea surface temperature in each of the four regions. The results of the principal component analysis are interpreted in terms of the underlying physical dynamics.

To understand and to predict the weather and climate, numerical models are important tools and it is crucial that the controlling processes are described correctly. Since 70% of the global surface is covered with water the description how the ocean and atmosphere communicates has a considerable impact. The ocean–atmosphere exchange occurs through transport of momentum (friction) and heat, governed by turbulent eddies. The sea surface is also an important source of turbulence in both directions. The scales of the turbulent eddies cannot be resolved in ocean and climate models. Therefore, the turbulent exchanges have to be related to mean variables, such as wind speed and temperature differences. By using measurements, new methods to describe the air–sea exchange during two specific processes were developed. These processes are the so-called UVCN-regime (Unstable Very Close to Neutral stratification) and swell, i.e. waves which are not produced by the local wind. These processes were included in an ocean model and in a regional atmospheric climate model and the impact was investigated. The UVCN-regime enhances the heat transport significantly during the autumn and winter months in the ocean model. This results in a shallower well-mixed surface layer in the ocean. Wind-following swell reduces the surface friction, which is very important for the atmosphere. Some secondary effects in the climate model are reduced low-level cloud cover and reduced precipitation by more than 10% over sea areas. Locally and for short periods the impact is large. It is important to include the UVCN-regime and the swell impact in models, to make simulations more reliable.

Dans les expériences de Milgram sur la soumission à l’autorité (Milgram, 1963, 1965, 1974), les participants se sont vus ordonner d’administrer une série de chocs électriques d’intensité croissante à un autre participant (en réalité compère de l’expérimentateur) au nom d’une étude sur les effets de la punition sur l’apprentissage. Les résultats montrent que 62.5% des participants ont été jusqu’à infliger plusieurs chocs potentiellement mortels (condition standard ; Milgram, 1974). Ces résultats ont suscité un fort intérêt et sont toujours largement cités pour expliquer certains comportements destructeurs comme les actes de torture et de barbarie. Mais les travaux de Milgram ont également provoqué une forte controverse éthique et toute possibilité de réplication a été rapidement proscrite. Dans ce contexte, peu d’études expérimentales ont été réalisées et la question des mécanismes responsables de l’obéissance destructrice (OD) demeure sans réponse. La recherche récente a pu relancer l’étude expérimentale de l’OD par l’usage d’environnements immersifs. Ainsi, une récente étude IRMf reposant sur l’utilisation d’une version virtuelle du paradigme de Milgram montre que l’observation de la douleur de la victime dans ce contexte provoque un état de détresse personnelle chez les participants (i.e., réaction émotionnelle aversive centrée sur soi). Ce résultat suggère que l’OD pourrait être en partie la conséquence d’un défaut de régulation de la détresse provoquée par les mécanismes de résonance empathique. En nous appuyant sur la recherche récente en neurosciences sociales, nous avons fait l’hypothèse que la vulnérabilité au stress pourrait faciliter l’OD via l’exercice d’un contrôle inhibiteur sur la résonance empathique responsable d’une diminution de l’aversion pour l’atteinte à autrui. Nous avons réalisé six expériences visant (i) à examiner l’influence du tonus vagal (biomarqueur de la vulnérabilité au stress) sur l’autoritarisme de droite (prédicteur classique de l’OD) et sur l’OD, (ii) à manipuler expérimentalement la capacité des participants à exercer un contrôle inhibiteur durant la procédure d’obéissance, (iii) à explorer la relation entre ondes thêta (biomarqueur du contrôle inhibiteur) et OD, (iv) à examiner la relation entre OD et activité hémodynamique au niveau du cortex préfrontal ventromédian (incluant le cortex orbitofrontal) et du cortex préfrontal dorsolatéral, régions cérébrales fortement impliquées dans l’empathie et la cognition morale. L’obéissance a été mesurée à l’aide de l’« Immersive Video Milgram Obedience Experiment ». Dans ce qu’ils ont d’essentiel, nos résultats montrent : (i) qu’un moindre tonus vagal prédit l’autoritarisme de droite et l’OD, et que les participants obéissants ont exercé un effort cognitif couplé à une diminution du stress physiologique durant la procédure d’obéissance (études 2 et 3), (ii) que l’affaiblissement expérimental du contrôle inhibiteur via l’induction d’une « fatigue mentale » favorise la désobéissance et supprime l’influence de l’autoritarisme de droite sur l’OD (étude 4), (iii) qu’une augmentation de la puissance des ondes thêta prédit l’OD (étude 5), (iv) qu’une augmentation de l’oxy-hémoglobine au niveau du cortex préfrontal ventromédian droit prédit une moindre obéissance (étude 6). Dans leur ensemble, ces résultats supportent l’hypothèse voulant que les personnes présentant une plus grande vulnérabilité au stress exercent un contrôle inhibiteur sur leur résonance empathique dans un effort pour diminuer leur détresse, et que ce contrôle inhibiteur a pour conséquence une diminution des réponses émotionnelles aversives à l’atteinte à autrui et ainsi une augmentation de l’OD
In the Milgram's obedience experiments (Milgram, 1963, 1965, 1974), naive participants were ordered to administer increasingly severe electric shocks on a “learner” (a confederate) after being told that they were participating in an experiment on the effects of punishment on learning. Results revealed that 62.5% of the participants were willing to administer allegedly lethal electric shocks when ordered to do so (standard condition; Milgram, 1974). The Milgram's findings are still often cited when explaining destructive behaviors such as torture. The Milgram’s obedience studies have also been a target of ethical criticism and replication has been discouraged. In such a context, a very few experimental studies has been conducted since the Milgram’s experiments and the mechanisms responsible for destructive obedience remain unknown. Recent research reopens the door to direct empirical study of destructive obedience through the employment of immersive environments. A recent fMRI study showed that pain-related affective sharing in a virtual version of the Milgram paradigm elicited an aversive, self-oriented state of personal distress. This result suggests that low self-regulatory control of the shared affect evoked by the victim’s pain could be responsible for destructive obedience. Based on recent social neuroscience research, we hypothesized that stress vulnerability may facilitate destructive obedience through a mechanism of inhibitory control over empathic resonance responsible for decreased harm aversion. We conducted six studies aiming (i) to explore the influence of cardiac vagal tone (a biomarker of stress vulnerability) on right-wing authoritarianism (RWA, a classic predictor of destructive obedience) and on destructive obedience, (ii) to induce a self-regulatory fatigue in order to manipulate the participants’ abilities for inhibitory control during the obedience procedure, (iii) to explore the relation between theta oscillations (a biomarker of inhibitory control) and destructive obedience, (iv) to examine the relation between destructive obedience and hemodynamic response in the ventromedial prefrontal cortex (including the orbitofrontal cortex) and the dorsolateral prefrontal cortex, two brain areas highly involved in empathy and moral cognition. Obedience was measured using the “Immersive Video Milgram Obedience Experiment”. All in all, our results showed: (i) that lower vagal tone predicted higher RWA and destructive obedience, and that obedient participants exerted a cognitive effort associated to decreased physiological arousal (studies 2 and 3), (ii) that self-regulatory fatigue reduced destructive obedience and suppressed the influence of RWA, (iii) that increased theta power predicted destructive obedience (study 4), (iv) that increased oxygenated-hemoglobin in the right ventromedial prefrontal cortex predicts disobedience. On the whole, these results support the hypothesis that individuals with high in stress vulnerability exert an inhibitory control over their empathic resonance in an attempt to reduce their own distress, and that such a mechanism is responsible for decreased harm aversion and then destructive obedience

Atmospheric models are strongly dependent on the turbulent exchange of momentum, sensible heat and moisture (latent heat) at the surface. Oceans cover about 70% of the Earth’s surface and understanding the processes that control air-sea exchange is of great importance in order to predict weather and climate. In the atmosphere, for instance, hurricane development, cyclone intensity and track depend on these processes. Ocean waves constitute an obvious example of air-sea interaction and can cause the air-flow over sea to depend on surface conditions in uniquely different ways compared to boundary layers over land. When waves are generated by wind they are called wind sea or growing sea, and when they leave their generation area or propagate faster than the generating wind they are called swell. The air-sea exchange is mediated by turbulent eddies occurring on many different scales. Field measurements and high-resolution turbulence resolving numerical simulations have here been used to study these processes. The standard method to measure turbulent fluxes is the eddy covariance method. A spatial separation is often used between instruments when measuring scalar flux; this causes an error which was investigated for the first time over sea. The error is typically smaller over ocean than over land, possibly indicating changes in turbulence structure over sea. Established and extended analysis methods to determine the dominant scales of momentum transfer was used to interpret how reduced drag and sometimes net upward momentum flux can persist in the boundary layer indirectly affected by swell. A changed turbulence structure with increased turbulence length scales and more effective mixing was found for swell. A study, using a coupled wave-atmosphere regional climate model, gave a first indication on what impact wave mixing have on atmosphere and wave parameters. Near surface wind speed and wind gradients was affected especially for shallow boundary layers, which typically increased in height from the introduced wave-mixing. A large impact may be expected in regions of the world with predominant swell. The impact of swell waves on air-sea exchange and mixing should be taken into account to develop more reliable coupled Earth system models.

O σχεδιασμός και η ανάπτυξη μιας σύγχρονης αεροναυπηγικής κατασκευής περιλαμβάνει ως επιμέρους φάσεις (μεταξύ άλλων) τον αρχικό και τον προκαταρκτικό σχεδιασμό. Οι φάσεις αυτές έχουν ιδιαίτερη σημασία διότι εκεί δίνεται η αρχική μορφή και οι διαστάσεις της κατασκευής. Είναι γεγονός ότι η συμβατική σχεδίαση των βασικών δομικών στοιχείων των αεροσκαφών έχει φτάσει σε πολύ υψηλό επίπεδο βελτιστοποίησης που επιδέχεται πλέον μόνο μικρά περιθώρια περαιτέρω βελτίωσης. Οι σύγχρονες όμως απαιτήσεις των ελαφρών κατασκευών, όπως δραστική μείωση του βάρους, αύξηση του ωφέλιμου φορτίου κτλ. ωθεί τις αεροναυπηγικές βιομηχανίες στη δημιουργία δομών που ξεφεύγουν από τις παραδοσιακές (μη-συμβατικές δομές). Παράλληλα με τα παραπάνω γίνεται προσπάθεια για μερική αντικατάσταση μεταλλικών υλικών από σύνθετα υλικά στις πρωτεύουσες δομές αεροναυπηγικών κατασκευών. Για να σχεδιαστούν και να εξελιχθούν μη-συμβατικές αεροναυπηγικές δομές χωρίς να καταφύγει κάποιος σε εκτενείς πειραματικές δοκιμές, η σύγχρονη τάση είναι η ανάπτυξη και ο συνδυασμός προτύπων συμπεριφοράς στη λογική της εξομοίωσης των πειραματικών δοκιμών. Η εξομοίωση αυτή επιτυγχάνεται με τη βοήθεια ηλεκτρονικών υπολογιστών και κατάλληλων μεθόδων βασισμένων στη θεωρία των πινάκων (Πεπερασμένα Στοιχεία, Συνοριακά Στοιχεία κλπ.). Στη φάση του αρχικού και προκαταρκτικού σχεδιασμού η εφαρμογή των μεθοδολογιών προσομοίωσης δεν είναι πάντοτε εύκολη και απλή, λόγω των πολλαπλών αλλαγών στη γεωμετρία, το υλικό και τις κατασκευαστικές λεπτομέρειες που πραγματοποιούνται στη δομή κατά την επαναληπτική διαδικασία του σχεδιασμού. Επομένως, η αποκλειστική χρήση αριθμητικών μεθόδων ανάλυσης καθίσταται αναποτελεσματική από άποψη χρονικών απαιτήσεων, αν δεν συνοδεύεται από αναλυτικές ή ημιαναλυτικές προσεγγίσεις επιμέρους προβλημάτων του σχεδιασμού. Βασικό μέρος του προκαταρκτικού σχεδιασμού μιας πτέρυγας μη συμβατικής δομής αποτελεί η αποφυγή της αστοχίας του άνω τμήματός της, διότι οι λεπτότοιχες ενισχυμένες με δοκούς πλάκες που χρησιμοποιούνται στην κατασκευή υφίστανται λυγισμό λόγω των θλιπτικών φορτίσεων που κυρίως παραλαμβάνουν. Η διαστασιολόγηση των σύνθετων πλακών που φέρουν δοκούς ενίσχυσης στις κατασκευές αυτές απαιτούν συνήθως πλήθος επαναληπτικών υπολογισμών για διαφορετικές γεωμετρίες, φορτίσεις, οριακές συνθήκες κλπ. Η εξέταση της κάθε περίπτωσης μεμονωμένα με τη χρήση αριθμητικών μεθόδων καθιστά την επίλυση ολόκληρης της κατασκευής εξαιρετικά χρονοβόρα. Για το λόγο αυτό, στη φάση της αρχικής θεωρητικής μελέτης και της αρχικής διαστασιολόγησης η χρησιμοποίηση αναλυτικών μεθόδων για την εύρεση του κρίσιμου φορτίου λυγισμού πλακών με δοκούς ενίσχυσης οδηγεί στην εξοικονόμηση υπολογιστικού κόστους. Επομένως, η ανάπτυξη αναλυτικών ή ημιαναλυτικών μεθόδων προσδιορισμού των φορτίων λυγισμού ενισχυμένων με δοκούς συνθέτων πλακών και κελυφών θεωρείται πολύ σημαντική. Για τον σκοπό αυτό, στο πλαίσιο αυτής της διατριβής, αναπτύσσονται αναλυτικές και ημιαναλυτικές λύσεις για το λυγισμό πολύστρωτων πλακών ενισχυμένων με ενισχυτικές διαμήκεις δοκούς, οι οποίες ενσωματώνονται σαν κριτήρια στη μέθοδο διαστασιολόγησης της δομής. Η μεθοδολογία συμπληρώνεται με πλήθος άλλων κατάλληλων κριτηρίων για τον έλεγχο της αντοχής των δομικών στοιχείων της πτέρυγας καθώς και με κριτήρια για την επαναδιαστασιολόγηση των στοιχείων κατά την επαναληπτική διαδικασία της βελτιστοποίησης. Με τη μεθοδολογία που αναπτύσσεται διερευνούνται διατάξεις δομής πτερύγων από σύνθετα υλικά με πολυάριθμες κύριες δοκούς. Πιο συγκεκριμένα, αναπτύσσονται αναλυτικές/ημιαναλυτικές λύσεις ολικού και τοπικού λυγισμού πλακών που φέρουν δοκούς ενίσχυσης. Όσον αφορά τον ολικό λυγισμό αναπτύσσεται μια μεθοδολογία που βασίζεται στη μαθηματική μετατροπή μιας πλάκας που φέρει δοκούς ενίσχυσης σε μια ισοδύναμη ομογενή πλάκα. Η αναπτυχθείσα μεθοδολογία ομογενοποίησης των ενισχυμένων πλακών εμφανίζει σημαντικά πλεονεκτήματα σε σύγκριση με τις αντίστοιχες ήδη υπάρχουσες. Παράλληλα, η ενεργειακή μέθοδος Rayleigh-Ritz εφαρμόζεται για τη λύση προβλημάτων λυγισμού μερικώς ανισότροπων πλακών με ενισχυτικές δοκούς από σύνθετα υλικά, λαμβάνοντας διακριτά υπόψη τις ενισχυτικές δοκούς. Όσον αφορά το πρόβλημα του τοπικού λυγισμού, αναπτύσσεται μια νέα μεθοδολογία για την εύρεση των κρίσιμων φορτίων τοπικού λυγισμού λεπτότοιχης πλάκας με χρήση ενεργειακών μεθόδων. Το μαθηματικό μοντέλο που χρησιμοποιείται για την περίπτωση του τοπικού λυγισμού της επικάλυψης είναι η απομόνωση του τμήματος της επικάλυψης μεταξύ δυο ενισχυτικών δοκών και η αντικατάσταση της δυσκαμψίας της υπόλοιπης πλάκας με ελατήρια μεταβλητής δυσκαμψίας. Η μεθοδολογία αυτή επεκτείνεται και στον προσδιορισμό της μεταλυγισμικής συμπεριφοράς μιας πλάκας ενισχυμένης με διαμήκεις δοκούς. Οι παραπάνω μεθοδολογίες υπολογισμού του κρίσιμου φορτίου λυγισμού που αναπτύσσονται, στα πλαίσια αυτής της διατριβής, εφαρμόζονται στη διαστασιολόγηση πτέρυγας μη συμβατικής δομής από σύνθετα υλικά με πολυάριθμες κύριες δοκούς, σε αντίθεση με τις συμβατικές πτέρυγες (με δύο κύριες δοκούς). Η ανάλυση τάσεων της πτέρυγας πραγματοποιείται με τη βοήθεια της μεθόδου των πεπερασμένων στοιχείων. Η τελική διαστασιολόγηση επιτυγχάνεται με επαναληπτική διαδικασία βελτιστοποίησης βασισμένη σε αναλυτικές και ημιαναλυτικές σχέσεις. Με τον τρόπο αυτό, συγκρίνεται λεπτομερώς η συμβατική δομή πτέρυγας με 2 κύριες δοκούς και οι αντίστοιχες πτέρυγες με 4, 5 και 6 κύριες δοκούς. Για την περαιτέρω βελτιστοποίηση της συμπεριφοράς της πτέρυγας, διερευνάται η επίδραση που έχει η αλλαγή των μηχανικών ιδιοτήτων του υλικού και των επιτρεπόμενων ορίων παραμόρφωσης στη δυνατότητα ελαχιστοποίησης της μάζας της πτέρυγας. Υπολογίστηκε ότι κάτω από συγκεκριμένες συνθήκες η χρήση της μη συμβατικής πτέρυγας μπορεί να οδηγήσει σε μείωση μάζας μέχρι και 12%.
The design and development of a modern aerospace structure consists of many design stages. The most important stages are the conceptual and the preliminary where the initial sizing of the structure is obtained. It is known that the conventional design of the aircraft’s main structural members has reached a high optimization level, where margins for further improvement are small. The current demands of the lightweight structures such as weight reduction, payload increase etc. have led the aerospace industries develop unconventional structures and partially substitute the metallic materials of the primary structures with composites. The current trend of designing and evolving unconventional aerospace structures, without performing extended experimental tests, leads to the development of behavior models. The simulation of the experimental tests (through the behavior models) is achieved using high performance computers and numerical methods (Finite Element Method, Boundary Element Method etc). To apply simulation methods during the conceptual and preliminary stage is not an easy task. Most of the difficulties are the numerous geometrical, material parameters and the structural details that alter during the iterative process of the design. So, the exclusive usage of numerical analysis methods becomes very time consuming, if it is not accompanied by analytical or semi analytical methods of the sub-problems of the design. Part of the preliminary design of an unconventional wing structure is to prevent upper skin from failure. The stiffened panels that comprise the upper skin of the wing suffer from buckling due to the applied compressive loads. The sizing of the composite stiffened panels usually requires numerous of iterative calculations for various geometries, loading and boundary conditions etc. The examination of each case separately, with the use of numerical methods, results to time consuming analyses of the entire structure. Therefore, the development of appropriate analytical or semi analytical methods for estimating stiffened panels’ critical buckling load is of great importance. For this purpose, in the present thesis, analytical and semi analytical methodologies are developed for estimating the critical buckling load of stiffened panels. The developed methodologies are incorporated as design criteria in the sizing routine of the entire structure. The sizing routine comprises additional sizing criteria for checking the strength of wing’s structural members at each phase of the iterative process. Applying the developed sizing routine in various wing configurations made of composite materials, multispar wing designs are studied. Specifically, analytical and semi analytical methods for global and local buckling problems of stiffened panels are developed. The methodology of global buckling problems is based on the mathematical conversion of a stiffened panel to an equivalent homogeneous panel. The developed method of homogenization of stiffened panels appears to have significant advantages over the already existed homogenization methods. Additionally, the energy method Rayleigh-Ritz is applied for solving global buckling problems of stiffened panels with partial anisotropy considering discrete stiffeners. Regarding local buckling problems of stiffened panels, a new methodology is developed for estimating the critical local buckling load with the use of energy methods. The approach considers the stiffened panel segment located between two stiffeners, while the remaining panel is replaced by equivalent transverse and rotational springs of varying stiffness, which act as elastic edge supports. The buckling analysis of the segment provides an accurate and conservative prediction of the panel local buckling behavior. Consequently, the developed methodology is extended in the prediction of post-buckling response of stiffened panels where skin has undergone local buckling. The developed methodologies for calculating the critical buckling load are applied for sizing the wing members of an unconventional wing (multispar configuration) from composite materials. An efficient methodology based on fast Finite Element (FE) stress analysis combined to analytically formulated design criteria is presented for the initial sizing of a large scale composite component. A detailed comparison between optimized designs of conventional (2-spar) and three alternative wing configurations which comprise 4-, 5-, and 6-spars for the wing construction is performed. In order to understand the effect of different material properties, as well as the variation of maximum strain level allowed in the total wing mass, parametric analyses are performed for all wing configurations considered. It arises that under certain conditions the multispar configuration demonstrates significant advantages over the conventional design. This would lead to a mass reduction of 12%.

碩士
建國科技大學
機械工程系暨製造科技研究所
100
A compact aerodynamics model combined computer-aided engineering analysis software was used to investigate the stresses of small scale vertical wind turbine blades subjected to the aerodynamics loads in this paper. The airfoil types used in the paper are NACA0012, NACA0015 and NACA0018 with the rotational speed 150 rpm and wind velocity 4 m/s. The model used to decide the aerodynamics loads on the blades was formulated with the blade element momentum method (BEM) and multiple stream-tubes method. The relations between lift-coefficient/drag-coefficient and angle of attack from X-FOIL are compared to experimental data to validate the data. Although the maximal errors are up to 50% in some ranges of Reynolds’ number and angle of attack, the comparison data shows that the two data are agree with well when Reynolds’ number is greater than 1.6 x 105 and the angle of attack less than 10°. The angle of attack, lift coefficient and drag coefficient are varied with the angles of rotation of blades, which change the stresses distribution of the blades. The results show that the stresses of blades varied with the angles of rotation of blades. The maximal von Mises stresses occurs when the blade chord is almost parallel to the direction of wind.

碩士
國立清華大學
動力機械工程學系
103
Wind blades are composite structure, which is mainly composed of two tough outer facesheets and with reinforcing spar structure in it. In this study, wind blades without spar, with one or three spars were analyzed using the finite element code ANSYS. The material of two facesheets and the reinforcing spar of 3kW wind blade are assumed to be aluminum alloy. The ratio of far and near airfoil was varied. The Whiffle-tree loading was adopted to simulate the wind load on the blade. The stress distribution and maximum deflection of the 3kW wind blade were evaluated. The results show that to reduce the maximum von Mises stress on the 3kW wind blade, increasing the number of reinforcing spar is better than increasing the ratio of far and near airfoil. In addition, to reduce the maximum deflection in the vertical direction, increasing the ratio of far and near airfoil is better than increasing the number of reinforcing spar. For the large-scale wind blade, the One-way FSI was used to calculate the pressure distribution on the NREL 5MW wind blade. The carbon fiber sheet/epoxy experiment was conducted for the material constants. The stress distribution and maximum deflection of the 5MW wind blade were evaluated under different structural support and material constants conditions. The results show that the NREL wind blade has the maxium stress and deflection at the attack angle of -90°, while the NREL wind blade has the minium stress and deflection at the attack angle of 0°, and the rigidity of NREL wind blade is enhanced by spar-cap. The results in both stress and deflection analysis are similar in Carbon fiber sheet/epoxy and glass fiber/vinyl epoxy, while the Carbon fiber sheet/epoxy has the lower density, which can significantly reduce the weight of the NREL blade.

碩士
國立成功大學
土木工程學系
104
Offshore wind turbine structure suffer from different types of loadings, such as wind and wave loads. The objective of this research is to analyze the stress fields of the grouted connection in the structure and predict the possible damage due to external forces. We used the finite element software ABAQUS to simulate the problem, and cohesive elements to describe the behavior of grout material. Stresses in the grouted connection due to different loading were shown. The results indicate the possibility of damage is dominated by shear. There also exists normal stress (opening) near the edges of the grouted connection when a horizontal load is acting. The damage zone may develop in the long term but further investigations and studies will be needed.

碩士
國立中央大學
機械工程學系
106
Energy issue is one of the most important topics in Taiwan, and so is the wind power. The size of wind turbine is getting larger to satisfy the needs of electric energy, but it results in an increment of cost and external loads. The wind-induced loads and stresses of Vestas V80-2MW wind turbine were analyzed. The author inspected the setting of inflow conditions first, and then investigated the wind-induced loads of the wind turbine under both normal power production and several extreme wind speeds. The influence of wind shear, turbulence intensity, and yaw misalignment were discussed. The results showed that the pitch angle of blade was at 0 position when the projected area of a blade was maximum. The wind-induced loads increased with turbulence intensity, and such increment of loads increased with mean wind speed. The yaw misalignment had an apparent effect on the lateral force acting on the tower. When the wind turbine was under extreme wind speed and parked, e.g., typhoon, the maximum lateral force occurred at 67.5 of yaw misalignment. With respect to the stress analysis of tower, the tower of wind turbine will not fail under normal power production or the specific wind speeds that were required in IEC 61400-1 standard. In modal analysis, the resonance of tower will not occur under normal power production with rated rotor speed. In stress analysis of bolt, the stress at bolt root was very high due to the effect of stress concentration, which resulted in local plastic deformation. The stress at the center of cross-section of bolt was identical to the value that was calculated by empirical formula.

The procedure of joining flexible or nonrigid parts using applied loads is called compliant assembly, and it is widely used in automotive, aerospace, electronics, and appliance manufacturing. Uncontrolled assembly processes may produce geometric errors that can exceed design tolerances and induce an increment of elastic energy in the structure due to the accumulation of internal stresses. This condition might create unexpected deformations and residual stress distributions across the structure that compromise product functionality. This thesis presents a method based on nonlinear Finite Element Analysis (FEA), metamodelling, and optimization techniques to provide accurate and on-time shimming strategies to support the definition of optimum assembly strategies. An example of the method on a typical aerospace wing box structure is demonstrated in the present study. The delivered outputs intend to support the production line by anticipating the response of the structure under a specific assembly condition and presenting alternative assembly strategies that can be applied to address eventual predicted issues on product requirements.
Graduate

碩士
國立清華大學
動力機械工程學系
104
Wind turbine blade is one of the most important parts of whole wind turbine. In order to investigate the operating condition of wind turbine, the modal analysis was carried out using finite element software ANSYS in this study, and the free vibration of NREL 5MW wind turbine blade was analyzed. Four different cases, including two materials of glass fiber/vinylester and carbon fiber fabric/epoxy, two supporting structures of reinforcing spar and box girder, were analyzed to explore their effects on the nature frequencies and mode shapes of NREL 5MW wind turbine blade. The results show that the blade made of carbon fiber fabric/epoxy composites has higher natural frequency for the torsional mode than that of blade made of glass fiber/vinylester composites. The natural frequencies of blade made of carbon fiber fabric/epoxy are higher than those of blade made of glass fiber/vinylester; these results can be used for wind blade to alter its resonance frequencies. In order to apply the wind load on the wind blade, the flow field analysis was carried out to determine the pressure distribution on the wind blade with four different pitch angles. The resulting pressure distribution was then inputted to the static structure analysis of the blade model as distributed forces on the wind blade. Wind turbine blades were also analyzed during operation, shutdown status and at the different positions to examine their stress distribution and deflection, and the Tsai-Hill failure criterion was used for failure prediction. The composite fiber orientations and blade skin thickness were investigated in the analysis. The results show that the wind blade has a maximum pressure and lift force at the pitch angle of 0°. When the pitch angle increases, the pressure on the wind blade decreases, thus the lift force reduces and the stress in the wind blade becomes smaller. The blade at the 120 degrees of right side of shaft from the vertical axis has a maximum stress and displacement. Four different fiber orientations of glass fiber/vinylester composites, 〖"[" 0_n "/±" 〖15〗_n "/" 〖90〗_n "]" 〗_s, 〖"[" 0_n "/±" 〖30〗_n "/" 〖90〗_n "]" 〗_s, 〖"[" 0_n "/±" 〖45〗_n "/" 〖90〗_n "]" 〗_s, and 〖"[" 0_n "/±" 〖60〗_n "/" 〖90〗_n "]" 〗_s, were studied and the blade with 60° fiber orientation had a minimum stress and displacement. When varying the skin thickness of wind blade, both of linear and optimal skin thickness can enhance the rigidity of blade, though increasing the weight of blade.

A simplified model is proposed to investigate the stress fields and the strain energy release rate (SERR) associated with cracks in bonded joints in wind turbine blades. The proposed two-dimensional model consists of nonparallel upper and lower shells with adhesive between the shells at the tapered end. Nonlinear finite element analysis (FEA) is performed in a systematic parametric study of material and geo- metric properties. Two failure modes and their locations are predicted at different combinations of parameters: yielding at the outside end of the adhesive and interface cracking at the inside end of the bondline. Effect of the shell curvature on the stress fields is also considered. Based on the classic beam theory and the beam-on-elastic-foundation (BOEF) assumption, stress and displacement fields of the adhesively-bonded joint were determined by a new theoretical model to support the results from the numerical computation. The failure analysis is continued by studying the effects of manufacturing defects in the adhesive bond. Single and multiple voids are embedded to simulate air bubble trapped in the interface. The numerical and analytical studies are conducted to investigate SERR associated with the voids and results are provided to illustrate the effects of void position and void size.

碩士
國立高雄海洋科技大學
輪機工程研究所
98
This research analyzes the loads and stresses of the turbine shaft of a mini-size vertical axis wind power generator. The main loads considered includes (1) the shaft lateral force due to the lift and drag forces generated by the rotation of the turbine blades, (2) the shaft lateral force due to the centrifugal force caused by unbalanced weights in turbine blades and frames, and (3) the axial resistant torque of the generator set. As the wind turbine is expected to be installed on an off-shore platform and worked with a magnet lift device, the turbine shaft must be anti-corrosive to the sea water, and should not interfere with the magnet field. Three material options, (1) pure SUS316 and (2) pure glass fiber reinforced composites, and (3) SUS316 enclosed by anti-corrosive glass fiber composites, are considered in the turbine shaft design analysis. Both fundamental solid mechanics and finite element method are used in the shaft stress analysis. The results show that the lower bearing of the turbine shaft is the component subject to the highest load. The normal stress of the turbine shaft is 3 times larger than the shear stress in the shaft, due to the high lateral forces produced by the turbine and relatively resistant torque in the generator set.

碩士
長庚大學
機械工程學系
100
According to the forecast of IEA 2010 World Energy Outlook (WEO 2010), the growth rate of the generating capacity of renewable energy (3.2 %) is more than that of the world's total generating capacity (2.4 %). In addition, having the advantages of mature technology, low maintenance costs, and promotions through legislation in many nations, wind power may become the secondary major renewable energy next to hydroelectric energy. The cost of the blades accounts for 20% in the overall cost of a wind turbine. In this study, the material used for the blades is RX06422, which is a mixture of nylon 66 with 50% of fiberglass. However, the blades are too strong when an auto-adjustment mechanism for blade-pitch is used. Therefore, we hope to adjust the fiberglass content and consequently reduce the mass, the cost, and the difficulty in manufacture of the blades through stress and strain analysis to the blade with blade-pitch adjustment mechanism. The stress and strain analysis was carried out using the commercial code ANSYS®, whose numerical model was verified with experimental data. The loadings from the flow field on the blade sections are predicted with the BEM (blade element momentum) method. The maximum stress under strong wind (25 m/s) would not exceed the strength of the blade if a small amount of blade-pitch variation is applied. Under the same wind speed, the blade-pitch variation must be large to keep the blade tip from hitting the mast due to blade bending. Therefore, the blade-tip displacement (13 cm) becomes the constrain of the analysis for choosing the blade material. For different materials, the Rule of Mixture is used to determine material properties. The Young's modului of the materials lie between 12 and 4 GPa, corresponding to 50 % and 0% of fiberglass content, respectively. The results show that for a maximum blade-pitch variation of 7, a fiberglass content of 17.23% is enough to prevent the blade from breaking and the blade-tip to hit the mast. The weight reduction of the blade with the new fiberglass content from the original one is about 26%.

Turbulent flow inside a street canyon was investigated in an open circuit wind tunnel and in a blow-down wind channel. Two geometries were used for comparison purposes: buildings with pitched roofs and with flat roofs. Both generate the flow of a different category, so the induced ventilation regimes are fundamentally different. Quadrant, Fourier and Wavelet analysis, Proper Orthogonal Decomposition (POD) and vortex detection methods are used to identify coherent structures in the flow and establish their impact on the ventilation of pollution. Two types of the organised motions are detected: the compact areas of sweep and ejection with the scale comparable to the size of building and the small vortices generated in the shear layer behind the building roof. POD identifies the most dominant modes with high coherency in the flow and evaluates the relative contributions of each mode to the overall kinetic energy of turbulence. Rigorous analysis of the correctness of the physical interpretation for such a decomposition is carried out. Wavelet analysis is applied to the time-series of the POD expansion coefficients in order to reveal control mechanism of the dynamics of the modes. Vorticity, calculated from the original velocity data, is decomposed by POD as well. Finally, the correlation between the vorticity...